SEARCH · Search USGS

Results for “Earth in Space”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42Linked to original sources

USGS2023

Predicting climate change impacts on poikilotherms using physiologically guided species abundance models

Poikilothermic animals comprise most species on Earth and are especially sensitive to changes in environmental temperatures. Species conservation in a changing climate relies upon predictions of species responses to future conditions, yet predicting species responses to climate change when temperatures exceed the bounds of observed data is fraught with challenges. We present a physiologically guided abundance (PGA) model that combines observations of species abundance and environmental conditions with laboratory-derived data on the physiological response of poikilotherms to temperature to predict species geographical distributions and abundance in response to climate change. The model incorporates uncertainty in laboratory-derived thermal response curves and provides estimates of thermal habitat suitability and extinction probability based on site-specific conditions. We show that temperature-driven changes in distributions, local extinction, and abundance of cold, cool, and warm-adapted species vary substantially when physiological information is incorporated. Notably, cold-adapted species were predicted by the PGA model to be extirpated in 61% of locations that they currently inhabit, while extirpation was never predicted by a correlative niche model. Failure to account for species-specific physiological constraints could lead to unrealistic predictions under a warming climate, including underestimates of local extirpation for cold-adapted species near the edges of their climate niche space and overoptimistic predictions of warm-adapted species.

Minnesota↗

USGS1989

Seismic methods for determining earthquake source parameters and lithospheric structure

The seismologic methods most commonly used in studies of earthquakes and the structure of the continental lithosphere are reviewed in three main sections: earthquake source parameter determinations, the determination of earth structure using natural sources, and controlled-source seismology. The emphasis in each section is on a description of data, the principles behind the analysis techniques, and the assumptions and uncertainties in interpretation. Rather than focusing on future directions in seismology, the goal here is to summarize past and current practice as a companion to the review papers in this volume. Reliable earthquake hypocenters and focal mechanisms require seismograph locations with a broad distribution in azimuth and distance from the earthquakes; a recording within one focal depth of the epicenter provides excellent hypocentral depth control. For earthquakes of magnitude greater than 4.5, waveform modeling methods may be used to determine source parameters. The seismic moment tensor provides the most complete and accurate measure of earthquake source parameters, and offers a dynamic picture of the faulting process. Methods for determining the Earth’s structure from natural sources exist for local, regional, and teleseismic sources. One-dimensional models of structure are obtained from body and surface waves using both forward and inverse modeling. Forward-modeling methods include consideration of seismic amplitudes and waveforms, but lack the formal resolution estimates obtained with inverse methods. Two- and three-dimensional lithospheric models are derived using various inverse methods, but at present most of these methods consider only traveltimes of body waves. Controlled-source studies of the Earth’s structure are generally divided by method into seismic refraction/wide-angle reflection and seismic reflection studies. Seismic refraction profiles are usually interpreted in terms of two-dimensional structure by forward modeling of traveltimes and amplitudes. The refraction method gives excellent estimates of seismic velocities, but relatively low resolution of structure. Formal resolution estimates are not possible for models derived from forward modeling, but informal estimates can be obtained by perturbing the best-fitting model. Inversion methods for seismic refraction data for one-dimensional models are well established, and two- and three-dimensional methods, including tomography, have recently been developed. Seismic reflection data provide the highest resolution of crustal structure, and have provided many important geological insights in the past decade. The acquisition and processing of these data have been greatly advanced by the hydrocarbon exploration industry. However, reliable crustal velocity control is generally lacking, and the origin of deep crustal reflections remains unclear, resulting in nonunique interpretations. A new form of lithospheric seismology has recently emerged that combines the advantages of seismic refraction and seismic reflection profiles, and the distinction between the two methods is steadily diminishing. Major challenges for future work will be the collection of data that are more densely sampled in space, and the development of interpretation methods that provide quantitative estimates of the uncertainties in the calculated models.

GSA Memoirs↗

USGS2026

Near-Earth solar wind speed of fast coronal mass ejections

Using solar wind and ground magnetometer data for solar cycles 23–25, extreme-value power-law models are developed that quantify relationships between near-Earth solar wind speed and Sun-to-Earth transit times of interplanetary coronal mass ejections (ICMEs). The models of near-Earth solar wind speed, conditional on ICME transit time, are (a) approximately stable for transit times as short as 15 hr, (b) sublinear, consistent with hydrodynamic drag acting on ICMEs during transit, with linear (ballistic) models not supported by the data, (c) representative of most of the information content of the data, and (d) insensitive to the solar cycle and interplanetary preconditioning. Applying the models to the September 1859 Carrington event ( 𝑇 =17.6 hours), we estimate an ICME-average speed of 𝑉 𝑎 =102⁢1 1164 896 km/s and a 1-hr ICME-maximum speed of 𝑉 𝑚 =157⁢5 1795 1382 km/s. These estimates indicate that Carrington-class magnetic storms do not require exceptionally extreme solar wind speeds, with values lower than some previous estimates.

JGR Space Physics↗

USGS2019

Space-based imaging radar studies of U.S. volcanoes

The arrival of space-based imaging radar as a revolutionary land-surface mapping and monitoring tool little more than a quarter century ago enabled a spate of innovative volcano research worldwide. Soon after launch of European Space Agency’s ERS-1 spacecraft in 1991, the U.S. Geological Survey began SAR and InSAR studies of volcanoes in the Aleutian and Cascades arcs, in Hawai’i, and elsewhere in the western U.S. including the Yellowstone and Long Valley calderas. This paper summarizes results of that effort and presents new findings concerning: (1) prevalence of volcano deformation in the Aleutian and Cascade arcs; (2) surface-change detection and hazard assessment during eruptions at Aleutian and Hawaiian volcanoes; (3) geodetic imaging of magma storage and transport systems in Hawai’i; and (4) deformation sources and processes at the Yellowstone and Long Valley calderas. Surface deformation caused by a variety of processes is common in arc settings and could easily escape detection without systematic InSAR surveillance. Space-based SAR imaging of active lava flows and domes in remote or heavily vegetated settings, including during periods of bad weather and darkness, extends land-based monitoring capabilities and improves hazards assessments. At Kīlauea Volcano, comprehensive SAR and InSAR observations identify multiple magma storage zones beneath the summit area and along the East Rift Zone, and illuminate magma transport pathways. The same approach at Yellowstone tracks the ascent of magmatic volatiles from a mid-crustal intrusion to shallow depth and relates that process to increased hydrothermal activity at the surface. Together with recent and planned launches of highly capable imaging-radar satellites, these findings support an optimistic outlook for near-real time surveillance of volcanoes at global scale in the coming decade.

Frontiers in Earth Science↗

USGS2018

Calculation of voltages in electric power transmission lines during historic geomagnetic storms: An investigation using realistic earth impedances

Commonly, one-dimensional (1-D) Earth impedances have been used to calculate the voltages induced across electric power transmission lines during geomagnetic storms under the assumption that much of the three-dimensional structure of the Earth gets smoothed when integrating along power transmission lines. We calculate the voltage across power transmission lines in the mid-Atlantic region with both regional 1-D impedances and 64 empirical 3-D impedances obtained from a magnetotelluric survey. The use of 3-D impedances produces substantially more spatial variance in the calculated voltages, with the voltages being more than an order of magnitude different, both higher and lower, than the voltages calculated utilizing regional 1-D impedances. During the March 1989 geomagnetic storm 62 transmission lines exceed 100 V when utilizing empirical 3-D impedances, whereas 16 transmission lines exceed 100 V when utilizing regional 1-D impedances. This demonstrates the importance of using realistic impedances to understand and quantify the impact that a geomagnetic storm has on power grids.

Space Weather↗

USGS2018

Pliocene erosional pulse and glacier-landscape feedbacks in the western Alaska Range

Pliocene–Pleistocene glaciation modified the topography and erosion of most middle- and high-latitude mountain belts, because the evolution of catchment topography controls long-term glacier mass balance and erosion. Hence, characterizing how erosion rates change during repeated glaciations can help test hypothesized glacier erosion-landscape feedbacks across a range of settings. To better understand how glaciations and landscapes coevolve on geologic timescales, I quantify erosion rates in the glaciated western Alaska Range with low-temperature thermochronometric data and modeling. Zircon (U–Th)/He and apatite fission track data suggest mountain-building was underway by the early Miocene. In contrast, lower-temperature apatite (U–Th)/He age-elevation and grain age-kinetic data indicate that erosion accelerated coincident with regional Pliocene glaciation ca. 4 Ma. Furthermore, erosion rates calculated within an eroding half-space indicate slow erosion at a rate ≤0.3 km/m.y. before 4.2 Ma, an initial pulse of rapid erosion at a rate of 1.0–1.6 km/m.y. during 4.2–2.9 Ma, and more moderate erosion at a rate of 0.4–0.7 km/m.y. since 2.9 Ma. The initial erosion pulse suggests a significant transient landscape adjustment to the introduction of efficient glacial erosion. The subsequent decrease in Pleistocene erosion rates is consistent with a negative feedback between continuing glaciation and glacier size/erosivity: If glacial erosion outpaces rock uplift, glacier erosion decreases over time as topography, mass balance, valley gradients, and ice flux are reduced. These findings imply that in areas of moderate rock uplift rates, the onset of local Plio–Pleistocene glaciation may have been punctuated by an initial pulse of rapid landscape change, after which change became more gradual.

Alaska↗

USGS1988

Glaciers of Antarctica

Of all the world's continents Antarctica is the coldest, the highest, and the least known. It is one and a half times the size of the United States, and on it lies 91 percent (30,109,800 km 3 ) of the estimated volume of all the ice on Earth. Because so little is known about Antarctic glaciers compared with what is known about glaciers in populated countries, satellite imagery represents a great leap forward in the provision of basic data. From the coast of Antarctica to about 81?south latitude, there are 2,514 Landsat nominal scene centers (the fixed geographic position of the intersection of orbital paths and latitudinal rows). If there were cloud-free images for all these geographic centers, only about 520 Landsat images would be needed to provide complete coverage. Because of cloud cover, however, only about 70 percent of the Landsat imaging area, or 55 percent of the continent, is covered by good quality Landsat images. To date, only about 20 percent of Antarctica has been mapped at scales of 1:250,000 or larger, but these maps do include about half of the coastline. The area of Antarctica that could be planimetrically mapped at a scale of 1:250,000 would be tripled if the available Landsat images were used in image map production. This chapter contains brief descriptions and interpretations of features seen in 62 carefully selected Landsat images or image mosaics. Images were chosen on the basis of quality and interest; for this reason they are far from evenly spaced around the continent. Space limitations allow less than 15 percent of the Landsat imaging area of Antarctica to be shown in the illustrations reproduced in this chapter. Unfortunately, a wealth of glaciological and other features of compelling interest is present in the many hundreds of images that could not be included. To help show some important features beyond the limit of Landsat coverage, and as an aid to the interpretation of certain features seen in the images, 38 oblique aerial photographs have been included. Again, these represent only a small fraction of the large number of aerial photographs now available in various national collections. The chapter is divided into five geographic sections. The first is the Transantarctic Mountains in the Ross Sea area. Some very large outlet glaciers flow from the East Antarctic ice sheet through the Transantarctic Mountains to the Ross Ice Shelf. Byrd Glacier, one of the largest in the world, drains an area of more than 1,000,000 km 2 . Next, images from the Indian Ocean sector are discussed. These include the Lambert Glacier- Amery Ice Shelf system, so large that about 25 images must be mosaicked to cover its complex system of tributary glaciers. Shirase Glacier, a tidal outlet glacier in the sector, flows at a speed of 2.5 km a -l . About 200 km inland and 200 km west of Shirase Glacier lie the Queen Fabiola ("Yamato") Mountains, whose extensive exposures of `blue ice' lay claim to being the world's most important meteorite-collecting locality, with more than 4,700 meteorite fragments discovered since 1969. The Atlantic Ocean sector is fringed by ice shelves into which flow large ice streams like Jutulstraumen, Stancomb-Wills, Slessor, and Recovery Glaciers. Filchner and Ronne Ice Shelves together cover an area two-thirds the size of Texas. From the western margin of the Ronne Ice Shelf, the north-trending arc of the Antarctic Peninsula, with its fjord and alpine landscape and fringing ice shelves, stretches towards South America. The Pacific Ocean sector begins with the Ellsworth Mountains, which include the highest peaks (Vinson Massif at 4,897 m) in Antarctica. The area between the Ellsworth Mountains and the eastern margin of the Ross Ice Shelf is fringed with small ice shelves and some major outlet glaciers. One of these, Pine Island Glacier, was found from comparing 1973 and 1975 images to have an average ice-front velocity of 2.4 km a -1 . This part of Antarctica contains many dormant volcanoes; the summits of servers, such as Mount Takahe with its 8-km-wide summit caldera, protrude through the West Antarctic ice sheet.

Professional Paper↗

USGS1992

The Galileo Solid-State Imaging experiment

The Solid State Imaging (SSI) experiment on the Galileo Orbiter spacecraft utilizes a high-resolution (1500 mm focal length) television camera with an 800 ?? 800 pixel virtual-phase, charge-coupled detector. It is designed to return images of Jupiter and its satellites that are characterized by a combination of sensitivity levels, spatial resolution, geometric fiedelity, and spectral range unmatched by imaging data obtained previously. The spectral range extends from approximately 375 to 1100 nm and only in the near ultra-violet region (??? 350 nm) is the spectral coverage reduced from previous missions. The camera is approximately 100 times more sensitive than those used in the Voyager mission, and, because of the nature of the satellite encounters, will produce images with approximately 100 times the ground resolution (i.e., ??? 50 m lp-1) on the Galilean satellites. We describe aspects of the detector including its sensitivity to energetic particle radiation and how the requirements for a large full-well capacity and long-term stability in operating voltages led to the choice of the virtual phase chip. The F/8.5 camera system can reach point sources of V(mag) ??? 11 with S/N ??? 10 and extended sources with surface brightness as low as 20 kR in its highest gain state and longest exposure mode. We describe the performance of the system as determined by ground calibration and the improvements that have been made to the telescope (same basic catadioptric design that was used in Mariner 10 and the Voyager high-resolution cameras) to reduce the scattered light reaching the detector. The images are linearly digitized 8-bits deep and, after flat-fielding, are cosmetically clean. Information 'preserving' and 'non-preserving' on-board data compression capabilities are outlined. A special "summation" mode, designed for use deep in the Jovian radiation belts, near Io, is also described. The detector is 'preflashed' before each exposure to ensure the photometric linearity. The dynamic range is spread over 3 gain states and an exposure range from 4.17 ms to 51.2 s. A low-level of radial, third-order, geometric distortion has been measured in the raw images that is entirely due to the optical design. The distortion is of the pincushion type and amounts to about 1.2 pixels in the corners of the images. It is expected to be very stable. We discuss the measurement objectives of the SSI experiment in the Jupiter system and emphasize their relationships to those of other experiments in the Galileo project. We outline objectives for Jupiter atmospheric science, noting the relationship of SSI data to that to be returned by experiments on the atmospheric entry Probe. We also outline SSI objectives for satellite surfaces, ring structure, and 'darkside' (e.g., aurorae, lightning, etc.) experiments. Proposed cruise measurement objectives that relate to encounters at Venus, Moon, Earth, Gaspra, and, possibly, Ida are also briefly outlined. The article concludes with a description of a 'fully distributed' data analysis system (HIIPS) that SSI team members intend to use at their home institutions. We also list the nature of systematic data products that will become available to the scientific community. Finally, we append a short 'historical' note outlining the responsibilities and roles of institutions and individuals that have been involved in the 14 year development of the SSI experiment so far. ?? 1992 Kluwer Academic Publishers.

Space Science Reviews↗

USGS2025

Permafrost–wildfire interactions: active layer thickness estimates for paired burned and unburned sites in northern high latitudes

As the northern high-latitude permafrost zone experiences accelerated warming, permafrost has become vulnerable to widespread thaw. Simultaneously, wildfire activity across northern boreal forest and Arctic/subarctic tundra regions impacts permafrost stability through the combustion of insulating organic matter, vegetation, and post-fire changes in albedo. Efforts to synthesis the impacts of wildfire on permafrost are limited and are typically reliant on antecedent pre-fire conditions. To address this, we created the FireALT dataset by soliciting data contributions that included thaw depth measurements, site conditions, and fire event details with paired measurements at environmentally comparable burned and unburned sites. The solicitation resulted in 52 466 thaw depth measurements from 18 contributors across North America and Russia. Because thaw depths were taken at various times throughout the thawing season, we also estimated end-of-season active layer thickness (ALT) for each measurement using a modified version of the Stefan equation. Here, we describe our methods for collecting and quality-checking the data, estimating ALT, the data structure, strengths and limitations, and future research opportunities. The final dataset includes 48 669 ALT estimates with 32 attributes across 9446 plots and 157 burned–unburned pairs spanning Canada, Russia, and the United States. The data span fire events from 1900 to 2022 with measurements collected from 2001 to 2023. The time since fire ranges from 0 to 114 years. The FireALT dataset addresses a key challenge: the ability to assess impacts of wildfire on ALT when measurements are taken at various times throughout the thaw season depending on the time of field campaigns (typically June through August) by estimating ALT at the end-of-season maximum. This dataset can be used to address understudied research areas, particularly algorithm development, calibration, and validation for evolving process-based models as well as extrapolating across space and time, which could elucidate permafrost–wildfire interactions under accelerated warming across the high-northern-latitude permafrost zone. The FireALT dataset is available through the Arctic Data Center ( https://doi.org/10.18739/A2RN3092P , Talucci et al., 2024).

Arctic↗

USGS2015

Resolving bathymetry from airborne gravity along Greenland fjords

Recent glacier mass loss in Greenland has been attributed to encroaching warming waters, but knowledge of fjord bathymetry is required to investigate this mechanism. The bathymetry in many Greenland fjords is unmapped and difficult to measure. From 2010 to 2012, National Aeronautics and Space Administration's Operation IceBridge collected a unique set of airborne gravity, magnetic, radar, and lidar data along the major outlet glaciers and fjords in Greenland. We applied a consistent technique using the IceBridge gravity data to create 90 bathymetric profiles along 54 Greenland fjords. We also used this technique to recover subice topography where warm or crevassed ice prevents the radar system from imaging the bed. Here we discuss our methodology, basic assumptions and error analysis. We present the new bathymetry data and discuss observations in six major regions of Greenland covered by IceBridge. The gravity models provide a total of 1950 line kilometers of bathymetry, 875 line kilometers of subice topography, and 12 new grounding line depths.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Combining numerical and statistical models to predict storm-induced dune erosion

Dune erosion is an important aspect to consider when assessing coastal flood risk, as dune elevation loss makes the protected areas more susceptible to flooding. However, most advanced dune erosion numerical models are computationally expensive, which hinders their application in early-warning systems. Based on a combination of probabilistic and process-based numerical modeling, we develop an efficient statistical tool to predict dune erosion during storms. The analysis focuses on Dauphin Island, AL in the northern Gulf of Mexico, where we combine synthetic sea storms with a calibrated and validated XBeach model to develop and test a range of different surrogate models for their ability to predict barrier-island geometric parameters under storm conditions. Surrogate models are developed by combining the oceanographic forcing from 100 optimally sampled sea-storm events covering the entire multivariate parameter space (used as XBeach input) and associated changes in the dune system (XBeach output). We test four surrogate models using a k-fold approach for validation. All models perform well in predicting changes in dune elevation, barrier-island area, and width, but are less accurate in predicting alterations in the cross-shore locations of dune morphological features. Multivariate adaptive regression splines are identified as the best surrogate model based on their fast development and good performance, attaining a modified Mielke index of 0.81 for dune crest height. As demonstrated at Dauphin Island, our approach shows potential to be used in an operational framework to predict dune response (in particular crest elevation change) when water level and wave forecasts are available.

Alabama↗

USGS2009

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Identifying coseismic subsidence in tidal-wetland stratigraphic sequences at the Cascadia subduction zone of western North America

Tidal-wetland stratigraphy reveals that great plate boundary earthquakes have caused hundreds of kilometers of coast to subside at the Cascadia subduction zone. However, determining earthquake recurrence intervals and mapping the coastal extent of past great earthquake ruptures in this region are complicated by the effects of many sedimentologic, hydrographic, and oceanographic processes that occur on the coasts of tectonically passive as well as active continental margins. Tidal-wetland stratigraphy at many Cascadia estuaries differs little from that at similar sites on passive-margin coasts where stratigraphic sequences form through nonseismic processes unrelated to coseismic land level changes. Methods developed through study of similar stratigraphic sequences in Europe provide a framework for investigating the Cascadia estuarine record. Five kinds of criteria must be evaluated when inferring regional coastal subsidence due to great plate boundary earthquakes: the suddenness and amount of submergence, the lateral extent of submerged tidal-wetland soils, the coincidence of submergence with tsunami deposits, and the degree of synchroneity of submergence events at widely spaced sites. Evaluation of such criteria at the Cascadia subduction zone indicates regional coastal subsidence during at least two great earthquakes. Evidence for a coseismic origin remains equivocal, however, for the many peat-mud contacts in Cascadia stratigraphic sequences that lack (1) contrasts in lithology or fossils indicative of more than half a meter of submergence, (2) well-studied tsunami deposits, or (3) precise ages needed for regional correlation. Paleoecologic studies of fossil assemblages are particularly important in estimating the size of sudden sea level changes recorded by abrupt peat-mud contacts and in helping to distinguish erosional and gradually formed contacts from coseismic contacts. Reconstruction of a history of great earthquakes for the Cascadia subduction zone will require rigorous application of the above criteria and many detailed investigations.

Journal of Geophysical Research B: Solid Earth↗

USGS2020

Measuring, modelling and projecting coastal land subsidence

Coastal subsidence contributes to relative sea-level rise and exacerbates flooding hazards, with the at-risk population expected to triple by 2070. Natural processes of vertical land motion, such as tectonics, glacial isostatic adjustment and sediment compaction, as well as anthropogenic processes, such as fluid extraction, lead to globally variable subsidence rates. In this Review, we discuss the key physical processes driving vertical land motion in coastal areas. Use of space-borne and land-based techniques and the associated uncertainties for monitoring subsidence are examined, as are physics-based models used to explain contemporary subsidence rates and to obtain future projections. Steady and comparatively low rates of subsidence and uplift owing to tectonic processes and glacial isostatic adjustment can be assumed for the twenty-first century. By contrast, much higher and variable subsidence rates occur owing to compaction associated with sediment loading and fluid extraction, as well as large earthquakes. These rates can be up to two orders of magnitude higher than the present-day rate of global sea-level rise. Multi-objective predictive models are required to account for the underlying physical processes and socio-economic factors that drive subsidence.

Nature Reviews Earth & Environment↗

USGS2007

Forecast experiment: do temporal and spatial b value variations along the Calaveras fault portend M ≥ 4.0 earthquakes?

The power law distribution of earthquake magnitudes and frequencies is a fundamental scaling relationship used for forecasting. However, can its slope (b value) be used on individual faults as a stress indicator? Some have concluded that b values drop just before large shocks. Others suggested that temporally stable low b value zones identify future large-earthquake locations. This study assesses the frequency of b value anomalies portending M ≥ 4.0 shocks versus how often they do not. I investigated M ≥ 4.0 Calaveras fault earthquakes because there have been 25 over the 37-year duration of the instrumental catalog on the most active southern half of the fault. With that relatively large sample, I conducted retrospective time and space earthquake forecasts. I calculated temporal b value changes in 5-km-radius cylindrical volumes of crust that were significant at 90% confidence, but these changes were poor forecasters of M ≥ 4.0 earthquakes. M ≥ 4.0 events were as likely to happen at times of high b values as they were at low ones. However, I could not rule out a hypothesis that spatial b value anomalies portend M ≥ 4.0 events; of 20 M ≥ 4 shocks that could be studied, 6 to 8 (depending on calculation method) occurred where b values were significantly less than the spatial mean, 1 to 2 happened above the mean, and 10 to 13 occurred within 90% confidence intervals of the mean and were thus inconclusive. Thus spatial b value variation might be a useful forecast tool, but resolution is poor, even on seismically active faults.

California↗

USGS2016

Landsat 8: The plans, the reality, and the legacy

Landsat 8, originally known as the Landsat Data Continuity Mission (LDCM), is a National Aeronautics and Space Administration (NASA)-U.S. Geological Survey (USGS) partnership that continues the legacy of continuous moderate resolution observations started in 1972. The conception of LDCM to the reality of Landsat 8 followed an arduous path extending over nearly 13 years, but the successful launch on February 11, 2013 ensures the continuity of the unparalleled Landsat record. The USGS took over mission operations on May 30, 2013 and renamed LCDM to Landsat 8. Access to Landsat 8 data was opened to users worldwide. Three years following launch we evaluate the science and applications impact of Landsat 8. With a mission objective to enable the detection and characterization of global land changes at a scale where differentiation between natural and human-induced causes of change is possible, LDCM promised incremental technical improvements in capabilities needed for Landsat scientific and applications investigations. Results show that with Landsat 8, we are acquiring more data than ever before, the radiometric and geometric quality of data are generally technically superior to data acquired by past Landsat missions, and the new measurements, e.g., the coastal aerosol and cirrus bands, are opening new opportunities. Collectively, these improvements are sparking the growth of science and applications opportunities. Equally important, with Landsat 7 still operational, we have returned to global imaging on an 8-day cycle, a capability that ended when Landsat 5 ceased operational Earth imaging in November 2011. As a result, the Landsat program is on secure footings and planning is underway to extend the record for another 20 or more years.

Remote Sensing of Environment↗

USGS2009

Gas hydrate saturations estimated from fractured reservoir at Site NGHP-01-10, Krishna-Godavari Basin, India

During the Indian National Gas Hydrate Program Expedition 01 (NGHP-Ol), one of the richest marine gas hydrate accumulations was discovered at Site NGHP-01-10 in the Krishna-Godavari Basin. The occurrence of concentrated gas hydrate at this site is primarily controlled by the presence of fractures. Assuming the resistivity of gas hydratebearing sediments is isotropic, th?? conventional Archie analysis using the logging while drilling resistivity log yields gas hydrate saturations greater than 50% (as high as ???80%) of the pore space for the depth interval between ???25 and ???160 m below seafloor. On the other hand, gas hydrate saturations estimated from pressure cores from nearby wells were less than ???26% of the pore space. Although intrasite variability may contribute to the difference, the primary cause of the saturation difference is attributed to the anisotropic nature of the reservoir due to gas hydrate in high-angle fractures. Archie's law can be used to estimate gas hydrate saturations in anisotropic reservoir, with additional information such as elastic velocities to constrain Archie cementation parameters m and the saturation exponent n. Theory indicates that m and n depend on the direction of the measurement relative to fracture orientation, as well as depending on gas hydrate saturation. By using higher values of m and n in the resistivity analysis for fractured reservoirs, the difference between saturation estimates is significantly reduced, although a sizable difference remains. To better understand the nature of fractured reservoirs, wireline P and S wave velocities were also incorporated into the analysis.

Journal of Geophysical Research B: Solid Earth↗

USGS1978

Megalineament in southeastern Alaska marks southwest edge of Coast Range batholithic complex

The Coast Range megalineament is a prominent, nearly continuous topographic and structural feature that extends southeastward about 550 km (330 mi) from its junction with the Chatham Strait – Lynn Canal fault at Point Sherman to Tongass Passage near the mouth of Pearse Canal where it leaves southeastern Alaska. It probably extends still further southeastward into British Columbia along Work Channel and Chatham Sound – Grenville Channel.The megalineament is a zone a few hundred metres to 10 km (6 mi) wide in which closely spaced joints, foliation, compositional layering, and small faults define the megalineament trend. The zone usually coincides with topographic depressions apparently caused by selective fluvial and glacial erosion of the less resistant rocks of the zone.Studies in the Juneau, Endicott Arm, and Behm Canal areas indicate that the megalineament (1) is locally the site of lateral and (or) vertical separations of no greater than several kilometres; (2) does not mark a major structural or metamorphic discontinuity in the near-surface rocks; (3) may be located near a pre-metamorphic and pre-intrusive discontinuity; (4) is consistently associated with and parallel to steep gradients in both the gravity and aeromagnetic fields; and (5) probably is the surface expression of the western contact, at depth, of the intrusive rocks and gneisses of the Coast Range batholithic complex with the schists to the southwest.

Alaska↗