USGS Science⌕ Search

SEARCH · USGS Science

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.

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

At least 127 records · Page 7Linked to original sources

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas↗

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗

Geomagnetic storms, the Dst ring-current myth and lognormal distributions

The definition of geomagnetic storms dates back to the turn of the century when researchers recognized the unique shape of the H-component field change upon averaging storms recorded at low latitude observatories. A generally accepted modeling of the storm field sources as a magnetospheric ring current was settled about 30 years ago at the start of space exploration and the discovery of the Van Allen belt of particles encircling the Earth. The Dst global 'ring-current' index of geomagnetic disturbances, formulated in that period, is still taken to be the definitive representation for geomagnetic storms. Dst indices, or data from many world observatories processed in a fashion paralleling the index, are used widely by researchers relying on the assumption of such a magnetospheric current-ring depiction. Recent in situ measurements by satellites passing through the ring-current region and computations with disturbed magnetosphere models show that the Dst storm is not solely a main-phase to decay-phase, growth to disintegration, of a massive current encircling the Earth. Although a ring current certainly exists during a storm, there are many other field contributions at the middle-and low-latitude observatories that are summed to show the 'storm' characteristic behavior in Dst at these observatories. One characteristic of the storm field form at middle and low latitudes is that Dst exhibits a lognormal distribution shape when plotted as the hourly value amplitude in each time range. Such distributions, common in nature, arise when there are many contributors to a measurement or when the measurement is a result of a connected series of statistical processes. The amplitude-time displays of Dst are thought to occur because the many time-series processes that are added to form Dst all have their own characteristic distribution in time. By transforming the Dst time display into the equivalent normal distribution, it is shown that a storm recovery can be predicted with remarkable accuracy from measurements made during the Dst growth phase. In the lognormal formulation, the mean, standard deviation and field count within standard deviation limits become definitive Dst storm parameters.

Journal of Atmospheric and Terrestrial Physics↗

Atlas of microscopic images of biochar using reflected light microscopy in biochar characterization

Derived through the thermochemical conversion of biomass, biochar is a carbon-rich substance recognized for its significance in environmental applications and sustainable agriculture. As interest in its utilization continues to rise, it becomes crucial to comprehend how the source material and pyrolysis parameters influence the properties of biochar and, consequently, to research the suitability of various analytical methods for characterizing it. Despite the current utilization of numerous physical and chemical methods, the untapped potential of reflected light microscopy warrants further exploration. While a few recent studies suggest a correlation between certain microscopic characteristics and selected physical and chemical properties of biochar, the data are limited and difficult to compare. This is primarily due to variations in the types of original biomass used and lack of information about pyrolysis conditions. Moreover, because only a limited number of photographs taken under a reflected light microscope are publicly available to-date, it is difficult to evaluate morphological differences between various biochars and other organic materials such as inertinites from coal, charcoal, etc. To address limited availability of publicly available data, this “Atlas of Microscopic Images of Biochar” presents a collection of more than 300 images contributed by researchers from Poland, the United States, Canada, Australia, Brazil, and Denmark. These photomicrographs capture optical characteristics of a diverse array of biochar, demonstrating its unique morphological and structural features. This visual documentation can serve as a valuable resource for researchers, industry professionals, educators, and enthusiasts interested in investigating the complexities of biochar forms.

Indiana Journal of Earth Sciences↗

Debris flow runup on vertical barriers and adverse slopes

Runup of debris flows against obstacles in their paths is a complex process that involves profound flow deceleration and redirection. We investigate the dynamics and predictability of runup by comparing results from large-scale laboratory experiments, four simple analytical models, and a depth-integrated numerical model (D-Claw). The experiments and numerical simulations reveal the important influence of unsteady, multidimensional flow on runup, and the analytical models highlight key aspects of the underlying physics. Runup against a vertical barrier normal to the flow path is dominated by rapid development of a shock, or jump in flow height, associated with abrupt deceleration of the flow front. By contrast, runup on sloping obstacles is initially dominated by a smooth flux of mass and momentum from the flow body to the flow front, which precedes shock development and commonly increases the runup height. D-Claw simulations that account for the emergence of shocks show that predicted runup heights vary systematically with the adverse slope angle and also with the Froude number and degree of liquefaction (or effective basal friction) of incoming flows. They additionally clarify the strengths and limitations of simplified analytical models. Numerical simulations based on a priori knowledge of the evolving dynamics of incoming flows yield quite accurate runup predictions. Less predictive accuracy is attained in ab initio simulations that compute runup based solely on knowledge of static debris properties in a distant debris flow source area. Nevertheless, the paucity of inputs required in ab initio simulations enhances their prospective value in runup forecasting.

Journal of Geophysical Research F: Earth Surface↗

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↗

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↗

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↗

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↗

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&ntilde;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↗

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↗

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↗

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↗

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗

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↗

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↗

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↗

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↗