SEARCH · Search USGS

Results for “Journal of Physics of the Earth”

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 109 records · Page 6Linked to original sources

USGS2018

Quantifying postfire aeolian sediment transport using rare earth element tracers

Grasslands, which provide fundamental ecosystem services in many arid and semiarid regions of the world, are undergoing rapid increases in fire activity and are highly susceptible to postfire-accelerated soil erosion by wind. A quantitative assessment of physical processes that integrates fire-wind erosion feedbacks is therefore needed relative to vegetation change, soil biogeochemical cycling, air quality, and landscape evolution. We investigated the applicability of a novel tracer technique—the use of multiple rare earth elements (REE)—to quantify soil transport by wind and to identify sources and sinks of wind-blown sediments in both burned and unburned shrub-grass transition zone in the Chihuahuan Desert, NM, USA. Results indicate that the horizontal mass flux of wind-borne sediment increased approximately threefold following the fire. The REE tracer analysis of wind-borne sediments shows that the source of the horizontal mass flux in the unburned site was derived from bare microsites (88.5%), while in the burned site it was primarily sourced from shrub (42.3%) and bare (39.1%) microsites. Vegetated microsites which were predominantly sinks of aeolian sediments in the unburned areas became sediment sources following the fire. The burned areas showed a spatial homogenization of sediment tracers, highlighting a potential negative feedback on landscape heterogeneity induced by shrub encroachment into grasslands. Though fires are known to increase aeolian sediment transport, accompanying changes in the sources and sinks of wind-borne sediments may influence biogeochemical cycling and land degradation dynamics. Furthermore, our experiment demonstrated that REEs can be used as reliable tracers for field-scale aeolian studies.

Journal of Geophysical Research: Biogeosciences↗

USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

USGS2014

An enhanced model of land water and energy for global hydrologic and earth-system studies

LM3 is a new model of terrestrial water, energy, and carbon, intended for use in global hydrologic analyses and as a component of earth-system and physical-climate models. It is designed to improve upon the performance and to extend the scope of the predecessor Land Dynamics (LaD) and LM3V models by better quantifying the physical controls of climate and biogeochemistry and by relating more directly to components of the global water system that touch human concerns. LM3 includes multilayer representations of temperature, liquid water content, and ice content of both snowpack and macroporous soil–bedrock; topography-based description of saturated area and groundwater discharge; and transport of runoff to the ocean via a global river and lake network. Sensible heat transport by water mass is accounted throughout for a complete energy balance. Carbon and vegetation dynamics and biophysics are represented as in LM3V. In numerical experiments, LM3 avoids some of the limitations of the LaD model and provides qualitatively (though not always quantitatively) reasonable estimates, from a global perspective, of observed spatial and/or temporal variations of vegetation density, albedo, streamflow, water-table depth, permafrost, and lake levels. Amplitude and phase of annual cycle of total water storage are simulated well. Realism of modeled lake levels varies widely. The water table tends to be consistently too shallow in humid regions. Biophysical properties have an artificial stepwise spatial structure, and equilibrium vegetation is sensitive to initial conditions. Explicit resolution of thick (>100 m) unsaturated zones and permafrost is possible, but only at the cost of long (≫300 yr) model spinup times.

Journal of Hydrometeorology↗

USGS2013

Seismicity around Parkfield correlates with static shear stress changes following the 2003 M w 6.5 San Simeon earthquake

Earthquakes trigger other earthquakes, but the physical mechanism of the triggering is currently debated. Most studies of earthquake triggering rely on earthquakes listed in catalogs, which are known to be incomplete around the origin times of large earthquakes and therefore missing potentially triggered events. Here we apply a waveform matched-filter technique to systematically detect earthquakes along the Parkfield section of the San Andreas Fault from 46 days before to 31 days after the nearby 2003 M w 6.5 San Simeon earthquake. After removing all possible false detections, we identify ~8 times more earthquakes than in the Northern California Seismic Network catalog. The newly identified events along the creeping section of the San Andreas Fault show a statistically significant decrease following the San Simeon main shock, which correlates well with the negative static stress changes (i.e., stress shadow) cast by the main shock. In comparison, the seismicity rate around Parkfield increased moderately where the static stress changes are positive. The seismicity rate changes correlate well with the static shear stress changes induced by the San Simeon main shock, suggesting a low friction in the seismogenic zone along the Parkfield section of the San Andreas Fault.

California↗

USGS2012

GFDL's ESM2 global coupled climate-carbon Earth System Models. Part I: physical formulation and baseline simulation characteristics

We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory's previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.

Journal of Climate↗

USGS2016

Sediment supply versus local hydraulic controls on sediment transport and storage in a river with large sediment loads

The Rio Grande in the Big Bend region of Texas, USA, and Chihuahua and Coahuila, Mexico, undergoes rapid geomorphic changes as a result of its large sediment supply and variable hydrology; thus, it is a useful natural laboratory to investigate the relative importance of flow strength and sediment supply in controlling alluvial channel change. We analyzed a suite of sediment transport and geomorphic data to determine the cumulative influence of different flood types on changing channel form. In this study, physically based analyses suggest that channel change in the Rio Grande is controlled by both changes in flow strength and sediment supply over different spatial and temporal scales. Channel narrowing is primarily caused by substantial deposition of sediment supplied to the Rio Grande during tributary-sourced flash floods. Tributary floods have large suspended-sediment concentrations, occur for short durations, and attenuate rapidly downstream in the Rio Grande, depositing much of their sediment in downstream reaches. Long-duration floods on the mainstem have the capacity to enlarge the Rio Grande, and these floods, released from upstream dams, can either erode or deposit sediment in the Rio Grande depending upon the antecedent in-channel sediment supply and the magnitude and duration of the flood. Geomorphic and sediment transport analyses show that the locations and rates of sand erosion and deposition during long-duration floods are most strongly controlled by spatial changes in flow strength, largely through changes in channel slope. However, spatial differences in the in-channel sediment supply regulate sediment evacuation or accumulation over time in long reaches (greater than a kilometer).

Big Bend National Park↗

USGS2022

Toward next-generation lava flow forecasting: Development of a fast, physics-based lava propagation model

During effusive volcanic crises, the eruption and propagation of lava flows pose a significant hazard to nearby populations, homes, and infrastructure. Consequently, timely lava flow forecasts are a critical need for volcano observatory and emergency management operations. Previous lava flow modeling tools are typically either too slow to produce timely forecasts, or are fast, but lack critical aspects of lava physics or important forecasting outputs. In particular, the strong thermal stratification present in laminar, high-Prandtl number flows has generally been neglected. Bulk rheological changes have previously been computed from cell-averaged temperatures, assuming that the flow is thermally mixed. Here, we detail the development and initial testing of Lava2d, a new two-dimensional depth-averaged finite volume model of lava flow propagation over natural terrain which accounts for bulk rheological changes due to thermorheological stratification. We use a novel approach to energy conservation based on tracking cooling and solidifying at the flow base and at the moving flow surface, allowing for the estimation of more realistic vertical thermorheological profiles, while maintaining computational efficiency, producing very timely model runs. We validate our approach with three examples: comparison with theoretical propagation of crust-dominated lava flows, comparison with a large-scale molten basalt experiment from the Syracuse University Lava Project, and efficiency testing and comparison with the initial phase of the 1984 Mauna Loa lava flows. Our model is shown to produce rapid, realistic forecasts, making it a good candidate for operationalization in active volcanic regions such as in Hawai'i.

Journal of Geophysical Research - Solid Earth↗

USGS2026

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE↗

USGS2012

Application of continuum models to deformation of the Aleutian Island Arc

Continuum models were constructed to describe large‐scale deformation of the Aleutian Island Arc over the past 5 m.y. These models consider the island arc as a continuum in the horizontal plane with the velocity boundary condition at the Pacific edge stated as a fraction of Pacific plate convergence transferred to the arc. First, a simple model of uniformly distributed strain is formulated to illustrate the mechanics of continuous deformation. Lineaments along the arc massif rotated about a vertical axis are matched by small‐element rotation calculated from the model. However, this model does not predict across‐arc variations in deformation and produces an unrealistic amount of crustal thickening after 5 m.y. A physically more meaningful model of deformation is the thin viscous sheet model based on averages of stress and rheology throughout the lithosphere. The amount of motion transferred from the Pacific plate to the arc is constrained by the rotated lineaments, while the effective stress‐strain exponent ( n ) and the ability the lithosphere has to sustain crustal thickness contrasts (the Argand number) are independent variables. Primarily, bathymetry, earthquake focal mechanisms, and styles of faulting are used to evaluate the models. The preferred model is one where the amount of motion transferred from the Pacific plate is greater in an arc‐parallel direction than in an arc‐normal direction, producing stresses consistent with strike‐slip faulting at the far western end of the arc and tensional stresses consistent with transverse normal faulting elsewhere in the arc massif. This model agrees with observations of slip vectors by Ekström and Engdahl (1989), who conclude that a portion of the arc‐parallel component of relative plate motion is taken up in the overriding plate. This model implies that compressive stress transferred to the arc is small in comparison to along‐arc shear stress and that stresses conducive to strikeslip faulting are prevalent throughout the arc.

Alaska↗

USGS2021

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗

USGS1996

40Ar/39Ar geochronology and paleomagnetism of Independence volcano, Absaroka volcanic supergroup, Beartooth mountains, Montana

Independence volcano is a major volcanic complex in the lower part of the Absaroka Volcanic Supergroup (AVS) of Montana and Wyoming. Recently reported Rb–Sr mineral dates from the complex give apparent ages of 91 and 84 Ma, whereas field relationships and the physical and compositional similarity of the rocks with other dated parts of the AVS indicate an Early to Middle Eocene age for eruption and deposition. To resolve the conflict between age assignments based on stratigraphic correlations and Rb–Sr dates, we report new paleomagnetic data and 40 Ar/ 39 Ar dates for Independence volcano. Paleomagnetic data for the stock and an andesite plug that cuts the stock are well grouped, of reverse polarity, and yield a virtual geomagnetic pole that is essentially identical to Late Cretaceous and Tertiary reference poles. The reverse polarity indicates that the magnetization of these rocks is probably younger than the Cretaceous normal superchron, or less than about 83.5 Ma. Hornblende from a volcanic breccia near the base of the volcanic pile gives a 40 Ar/ 39 Ar age of 51.57 Ma, whereas biotites from a dacite sill and a granodiorite stock that forms the core of the volcano give dates that range from 49.96 to 48.50 Ma. These dates record the age of eruption and intrusion of these rocks and clearly show that the age of Independence volcano is Early to Middle Eocene, consistent with stratigraphic relations. We suggest that the Rb–Sr mineral dates from the Independence stock and related intrusions are unreliable.

Montana, Wyoming↗

USGS2012

Effects of river regulation on aeolian landscapes, Colorado River, southwestern USA

Connectivity between fluvial and aeolian sedimentary systems plays an important role in the physical and biological environment of dryland regions. This study examines the coupling between fluvial sand deposits and aeolian dune fields in bedrock canyons of the arid to semiarid Colorado River corridor, southwestern USA. By quantifying significant differences between aeolian landscapes with and without modern fluvial sediment sources, this work demonstrates for the first time that the flow- and sediment-limiting effects of dam operations affect sedimentary processes and ecosystems in aeolian landscapes above the fluvial high water line. Dune fields decoupled from fluvial sand supply have more ground cover (biologic crust and vegetation) and less aeolian sand transport than do dune fields that remain coupled to modern fluvial sand supply. The proportion of active aeolian sand area also is substantially lower in a heavily regulated river reach (Marble–Grand Canyon, Arizona) than in a much less regulated reach with otherwise similar environmental conditions (Cataract Canyon, Utah). The interconnections shown here among river flow and sediment, aeolian sand transport, and biologic communities in aeolian dunes demonstrate a newly recognized means by which anthropogenic influence alters dryland environments. Because fluvial–aeolian coupling is common globally, it is likely that similar sediment-transport connectivity and interaction with upland ecosystems are important in other dryland regions to a greater degree than has been recognized previously.

Colorado↗

USGS1996

Quantitative investigations of the Missouri gravity low: A possible expression of a large, Late Precambrian batholith intersecting the New Madrid seismic zone

Analysis of gravity and magnetic anomaly data helps characterize the geometry and physical properties of the source of the Missouri gravity low, an important cratonic feature of substantial width (about 125 km) and length (> 600 km). Filtered anomaly maps show that this prominent feature extends NW from the Reelfoot rift to the Midcontinent Rift System. Geologic reasoning and the simultaneous inversion of the gravity and magnetic data lead to an interpretation that the gravity anomaly reflects an upper crustal, 11-km-thick batholith with either near vertical or outward dipping boundaries. Considering the modeled characteristics of the batholith, structural fabric of Missouri, and relations of the batholith with plutons and regions of alteration, a tectonic model for the formation of the batholith is proposed. The model includes a mantle plume that heated the crust during Late Precambrian and melted portions of lower and middle crust, from which the low-density granitic rocks forming the batholith were partly derived. The batholith, called the Missouri batholith, may be currently related to the release of seismic energy in the New Madrid seismic zone (earthquake concentrations occur at the intersection of the Missouri batholith and the New Madrid seismic zone). Three qualitative mechanical models are suggested to explain this relationship with seismicity. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS1976

Improving estimates of streamflow characteristics by using Landsat-1 imagery

Imagery from the first Earth Resources Technology Satellite (renamed Landsat-1) was used to discriminate physical features of drainage basins in an effort to improve equations used to estimate streamflow characteristics at gaged and ungaged sites. Records of 20 gaged basins in the Delmarva Peninsula of Maryland, Delaware, and Virginia were analyzed for 40 statistical streamflow characteristics. Equations relating these characteristics to basin characteristics were obtained by a technique of multiple linear regression. A control group of equations contains basin characteristics derived from maps. An experimental group of equations contains basin characteristics derived from maps and imagery. Characteristics from imagery were forest, riparian (streambank) vegetation, water, and combined agricultural and urban land use. These basin characteristics were isolated photographically by techniques of film-density discrimination. The area of each characteristic in each basin was measured photometrically. Comparison of equations in the control group with corresponding equations in the experimental group reveals that for 12 out of 40 equations the standard error of estimate was reduced by more than 10 percent. As an example, the standard error of estimate of the equation for the 5-year recurrence-interval flood peak was reduced from 46 to 32 percent. Similarly, the standard error of the equation for the mean monthly flow for September was reduced from 32 to 24 percent, the standard error for the 7-day, 2-year recurrence low flow was reduced from 136 to 102 percent, and the standard error for the 3-day, 2-year flood volume was reduced from 30 to 12 percent. It is concluded that data from Landsat imagery can substantially improve the accuracy of estimates of some streamflow characteristics at sites in the Delmarva Peninsula.

Maryland, Delaware, Virginia↗

USGS1981

Progradational sequences in Miocene shoreline deposits, southeastern Caliente Range, California

An exceptionally well exposed marine-nonmarine transition in middle Miocene strata exists in the southeastern Caliente Range, California. About 50 individual progradational sequences form a succession that ranges in thickness from approximately 1000 m (where predominantly nonmarine) to more than 2500 m (where predominantly marine). Paleogreographic evidence in basalt flows near the top of the succession and in overlying fluvial deposists indicates that these middle Miocene strata were deposited across a north-northwest trending shoreline. A complete progradational sequence typically is several meters to a few tens of meters thick and includes strata that represent three intertonguing stratigraphic units. Individual sequences generally rest on a thin gravel deposit interpreted as a transgressive lag on an erosional surface. The gravel is overlain by structureless siltstone or fine-grained sandstone deposited at water depths where the rate of faunal mixing exceeded that of production of structures by physical processes. These rocks grade upward into bedded fine sandstone deposited closer to shore where physical processes exceeded bioturbation. Crossbedded lenses of coarse sand or fine gravel in the upper part of this facies suggest the presence of failry long-period surface waves. The bedded fine sandstone is sharply overlain by a crossbedded coarse sandstone facies that is interpreted as a combined offshore bar-rip channel-surf zone assemblage. Cross-strata dip dominantly offshore, suggesting substantial deposition from rip currents. A secondary, shore=parallel mode of cross-strata direction suggests longshore currents produced by surface waves from the northwest. The crossbedded coarse-grained sandstone grades upward into planar-bedded medium-grained sandstone that is interpreted as a beach foreshore. This facies grades upward through structureless medium-grained sandstone into nonmarine or lagoonal red and green mudstone of the Caliente Formation. The middle Miocene succession was deposited in a subsiding basin that was otherwise remarkably stable tectonically; the position of the strand line differed no more than a few kilometers through a period of 1 to 3 m.y. The average duration of the transgressive-regressive cycles, a few tens of thousands of years, together with their distribution in groups of three or four in the lower two-thirds of the succession, is consistent with the pattern of long-term climatic cycles produced by periodicity of the earth's solar orbit and may be related to eustatic sea level changes attendant to the development of the Antarctic ice cap. Changes in the pattern of progradation in the upper part of the succession and nearby basaltic eruptions may hav been precursors to the onset of movement along the San Andreas fault in this area 12-14 m.y. ago.

California↗

USGS2006

The changing shapes of active volcanoes: History, evolution, and future challenges for volcano geodesy

At the very heart of volcanology lies the search for the 'plumbing systems' that form the inner workings of Earth’s active volcanoes. By their very nature, however, the magmatic reservoirs and conduits that underlie these active volcanic systems are elusive; mostly they are observable only through circumstantial evidence, using indirect, and often ambiguous, surficial measurements. Of course, we can infer much about these systems from geologic investigation of materials brought to the surface by eruptions and of the exposed roots of ancient volcanoes. But how can we study the magmatic processes that are occurring beneath Earth’s active volcanoes? What are the geometry, scale, physical, and chemical characteristics of magma reservoirs? Can we infer the dynamics of magma transport? Can we use this information to better forecast the future behavior of volcanoes? These questions comprise some of the most fundamental, recurring themes of modern research in volcanology. The field of volcano geodesy is uniquely situated to provide critical observational constraints on these problems. For the past decade, armed with a new array of technological innovations, equipped with powerful computers, and prepared with new analytical tools, volcano geodesists have been poised to make significant advances in our fundamental understanding of the behavior of active volcanic systems. The purpose of this volume is to highlight some of these recent advances, particularly in the collection and interpretation of geodetic data from actively deforming volcanoes. The 18 papers that follow report on new geodetic data that offer valuable insights into eruptive activity and magma transport; they present new models and modeling strategies that have the potential to greatly increase understanding of magmatic, hydrothermal, and volcano-tectonic processes; and they describe innovative techniques for collecting geodetic measurements from remote, poorly accessible, or hazardous volcanoes. To provide a proper context for these studies, we offer a short review of the evolution of volcano geodesy, as well as a case study that highlights recent advances in the field by comparing the geodetic response to recent eruptive episodes at Mount St. Helens. Finally, we point out a few areas that continue to challenge the volcano geodesy community, some of which are addressed by the papers that follow and which undoubtedly will be the focus of future research for years to come.

Journal of Volcanology and Geothermal Research↗

USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Slipstream: an early Holocene slump and turbidite record from the frontal ridge of the Cascadia accretionary wedge off western Canada and paleoseismic implications

Slipstream Slump, a well-preserved 3 km wide sedimentary failure from the frontal ridge of the Cascadia accretionary wedge 85 km off Vancouver Island, Canada, was sampled during Canadian Coast Guard Ship (CCGS) John P. Tully cruise 2008007PGC along a transect of five piston cores. Shipboard sediment analysis and physical property logging revealed 12 turbidites interbedded with thick hemipelagic sediments overlying the slumped glacial diamict. Despite the different sedimentary setting, atop the abyssal plain fan, this record is similar in number and age to the sequence of turbidites sampled farther to the south from channel systems along the Cascadia Subduction Zone, with no extra turbidites present in this local record. Given the regional physiographic and tectonic setting, megathrust earthquake shaking is the most likely trigger for both the initial slumping and subsequent turbidity currents, with sediments sourced exclusively from the exposed slump face of the frontal ridge. Planktonic foraminifera picked from the resedimented diamict of the underlying main slump have a disordered cluster of 14 C ages between 12.8 and 14.5 ka BP. For the post-slump stratigraphy, an event-free depth scale is defined by removing the turbidite sediment intervals and using the hemipelagic sediments. Nine 14 C dates from the most foraminifera-rich intervals define a nearly constant hemipelagic sedimentation rate of 0.021 cm/year. The combined age model is defined using only planktonic foraminiferal dates and Bayesian analysis with a Poisson-process sedimentation model. The age model of ongoing hemipelagic sedimentation is strengthened by physical property correlations from Slipstream events to the turbidites for the Barkley Canyon site 40 km south. Additional modelling addressed the possibilities of seabed erosion or loss and basal erosion beneath turbidites. Neither of these approaches achieves a modern seabed age when applying the commonly used regional marine 14 C reservoir age of 800 years (marine reservoir correction Δ R = 400 years). Rather, the top of the core appears to be 400 years in the future. A younger marine reservoir age of 400 years (Δ R = 0 years) brings the top to the present and produces better correlations with the nearby Effingham Inlet paleo-earthquake chronology based only on terrestrial carbon requiring no reservoir correction. The high-resolution dating and facies analysis of Slipstream Slump in this isolated slope basin setting demonstrates that this is also a useful type of sedimentary target for sampling the paleoseismic record in addition to the more studied turbidites from submarine canyon and channel systems. The first 10 turbidites at Slipstream Slump were deposited between 10.8 and 6.6 ka BP, after which the system became sediment starved and only two more turbidites were deposited. The recurrence interval for the inferred frequent early Holocene megathrust earthquakes is 460 ± 140 years, compatible with other estimates of paleoseismic megathrust earthquake occurrence rates along the subduction zone.

Slipstream Slump↗