Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Near Surface Geophysics”

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.

735 records · Page 37Linked to original sources
USGS2026

Assessing future hydrologic extremes using an integrated hydrology and river operations model in the Russian River watershed

Study region The Russian River watershed, situated in coastal, northern California, experiences hydrologic extremes, including periodic droughts and flooding. Water managers are working to maintain sustainable water supplies and environmental flows, while mitigating flood risks. Study focus This paper introduces an integrated hydrology and river operations model for the Russian River watershed. This model is distinct from models in previous studies because it represents surface-groundwater interactions and uses climate forcings to estimate dynamic water use demands that are superimposed onto both reservoir operations and water supply constraints. The model was used to examine three historical (1990–2015) and eight future (2016–2099) water use and climate change scenarios. New hydrological insights for the region The direct connection between streams and aquifers facilitated both annual aquifer replenishment by high winter streamflows and streamflow depletion by groundwater wells (19 % of pumped groundwater in alluvial aquifers from stream leakage) during critical low flow periods. Simulated streamflow changes included 59 % longer and 54 % more severe streamflow droughts, 26 % lower seasonal low streamflows, and up to 125 % higher peak streamflows, averaged over future climate and water use scenarios, suggesting increased future flood and water availability risks. Results showed the importance of reservoir operations for mitigating the impacts of increased hydroclimatic volatility, despite a decrease in reservoir reliability at Lake Mendocino, suggesting that reservoir management may be used to decrease future risks.

California↗
USGS2025

Pathways for potential exposure to onshore oil and gas wastewater: What we need to know to protect human health

Produced water is a chemically complex waste stream generated during oil and gas development. Roughly four trillion liters were generated onshore in the United States in 2021 (ALL Consulting, 2022, https://www.gwpc.org/wp-content/uploads/2021/09/2021_Produced_Water_Volumes.pdf ). Efforts are underway to expand historic uses of produced water to offset freshwater needs in water-stressed regions, avoid induced seismic activity associated with its disposal, and extract commodities. Understanding the potential exposures from current and proposed produced water uses and management practices can help to inform health-protective practices. This review summarizes what is known about potential human exposure to produced water from onshore oil and gas development in the United States. We synthesize 236 publications to create a conceptual model of potential human exposure that illustrates the current state of scientific inquiry and knowledge. Exposure to produced water can occur following its release to the environment through spills or leaks during its handling and management. Exposure can also arise from authorized releases, including permitted discharges to surface water, crop irrigation, and road treatment. Knowledge gaps include understanding the variable composition and toxicity of produced water released to the environment, the performance of treatment methods, migration pathways through the environment that can result in human exposure, and the significance of the exposures for human and ecosystem health. Reducing these uncertainties may help in realizing the benefits of produced water use while simultaneously protecting human health.

Environmental Health Perspectives↗
USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗
USGS2026

Relating glassy rind thicknesses to ambient air temperatures at the Lost Jim flow field in the Imuruk Lake volcanic field, Alaska

The Lost Jim flow field, in the Imuruk Lake volcanic field, Alaska, extends west ~34 km from a single vent, crossing subarctic tundra and currently touches several lakes and streams. The weighted mean of five 36Cl cosmogenic exposure ages from the Lost Jim pāhoehoe flow is 7.73 ± 0.37 ka, indicating this eruption occurred substantially after the eruption of the underlying Camille flow, which was emplaced at 39.7 ± 1.3 ka. Paleoclimate records indicate the period when the Lost Jim flow field was emplaced was after deglaciation, and the climate was similar to today. We propose that the emplacement of lava in these cold subarctic conditions can lead to faster cooling of the lava surface compared to lava emplaced in warmer locations such as mid- latitude cold deserts. Glass abundance in the outermost rinds at the Lost Jim flow field was on average 74 % with 6.4 mm thick rims, compared to 60 % with 2.9 mm rims for cold mid-latitude desert samples. We interpret increased glass content as a proxy for rapid cooling likely occurring partly during winter. Glassiness values varied less across vent, margin, and mid-flow locations when compared to the mid-latitude flows suggesting the Lost Jim flow field was broadly impacted by the subarctic climate as opposed to responding to local microclimates. Our results indicate that lava glassiness may be a useful environmental indicator of cooler (in this case subarctic) conditions.

Alaska↗
USGS2023

Preliminary geologic map of the Southern Santa Rosa Mountains and Borrego Badlands, San Diego County, Southern California

This investigation delineates the geologic framework of an area of 75 square kilometers (km 2 ) located west of the Salton Sea in southern California (fig. 1, on sheet 1). The study area encompasses the south flank of the Santa Rosa Mountains and the eastern part of the Borrego Badlands (sheet 1). In this study area, regionally important stratigraphic and structural elements collectively inform the late Cenozoic geologic evolution of the Anza-Borrego sector of the Salton Trough province. Critical stratigraphic and structural elements in the map area include the following: The well exposed sequence of late Cenozoic, nonmarine sedimentary rocks that filled the Anza-Borrego subbasin (fig. 1) of the Salton Trough; A tectonic boundary that—in the southern Santa Rosa Mountains—separates the sedimentary strata from underlying crystalline rocks of Peninsular Ranges type. This tectonic boundary, named the West Salton Detachment Fault System by Axen and Fletcher (1998), is projected to underlie all late Cenozoic sedimentary strata in the Anza-Borrego subbasin of the Salton Trough; A variety of transpressional, transtensional, and strike-slip structures that have deformed the late Cenozoic sedimentary strata and collectively guided syntectonic and posttectonic depositional events within the Anza-Borrego subbasin of the Salton Trough; and The southeasternmost surface expression of the Clark Fault, a major strand of the dextral San Jacinto Fault Zone. Geologic mapping and analysis for this investigation focused on clarifying geologic relations among these four stratigraphic and structural aspects in the map area.

California↗
USGS2025

Beyond the mangroves: A global synthesis of tidal forested wetland types, drivers and future information opportunities

There is increasing awareness of the global diversity of tidal forested wetlands (TFWs) and their significance in the provision of ecosystem services. These ecosystems, including mangrove forests, tidal freshwater forested wetlands, supratidal forests and transitional forests together span multiple climatic zones, geomorphic settings, and inundation and salinity regimes. We utilise case studies across five continents to demonstrate the state of knowledge among TFWs. Intertidal mangroves are the best-defined of the TFWs thanks to decades of research on their geomorphology, hydrology and ecology across their broad distribution. Non-mangrove forest settings, however, demonstrate more diverse hydrological, biochemical and vegetation conditions. In many cases, non-mangrove forests are situated at upper intertidal or supratidal elevations, where surface waters and groundwater are subject to interactions between tides freshwater inputs. Salinity datasets show variations ranging from tidal freshwater forested wetlands and ‘low-salinity mangroves’ to mesohaline or marine salinities, often with high temporal variability. While the floristic composition of non-mangrove forests vary among biogeographic regions, locally dominant TFW species are commonly distributed beyond the tidal niche into non-tidal wetland and upland forests. This presents challenges for traditional remote sensing approaches to ecosystem mapping, which are mostly lacking for non-mangrove forests. Geomorphic approaches and developments in machine learning offer opportunities to address this.

Preprint↗
USGS2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi↗
USGS2024

Using surficial geologic maps, vegetation, and monitoring to address erosion impacts from grazing in Channel Islands National Park, California

Employing a map-unit classification scheme based on geomorphic process and age, the U.S. Geological Survey (USGS) mapped Quaternary surficial deposits of the five islands comprising Channel Islands National Park (CHIS), California, as no such maps previously existed. Mapping was motivated through an agreement with the National Park Service (NPS) to aid natural resource assessments, including post-grazing disturbance recovery and identification of mass wasting and tectonic hazards. The resulting detailed (1:12,000 scale) maps portray areas of upland erosional transport processes and alluvial, fluvial, eolian, beach, marine terrace, mass wasting, and mixed depositional processes at a scale commensurate with the process signal on the landscape. Detailed, GPS-located observations of sedimentology, pedogenic development, landscape position, and geomorphology constrain map unit identification. In addition, map-unit boundary delineation was determined in a GIS using a variety of high-resolution (sub-meter) aerial imagery and lidar-based DEMs complemented by derivative raster products. The GIS geodatabase schema used is GeMS compliant. Resolvable map-unit widths are as small as 10-15 meter. Locational accuracy of unit boundaries is estimated to be as low as 10 m or better. Detailed mapping is at a scale allowing determination of geomorphic transport processes operating on the landscape and conveys geologic information to land managers at a resolution commensurate with park resources. Map unit ages are constrained by relative soil chronostratigraphy, vertebrate and invertebrate fossils, and materials collected for dating control using radiogenic, luminescence, and cosmogenic radionuclide techniques. Associated topical studies examined connections between geology, surface water, and vegetation distribution and restoration.

California↗
USGS2026

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences↗
USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2026

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California↗
USGS2025

A 10-year continuous daily simulation of chloride flux from a suburban watershed in Fairfax County, Virginia, USA

Increasing levels of chloride in surface water are associated with detrimental effects on water quality, aquatic ecosystems, infrastructure, and human health. Numerous mass-balance studies have inferred watershed transport processes by interpreting chloride inputs and outputs, but few represent internal dynamics explicitly. We constructed a coupled water/chloride mass balance model to gain insights into storage, residence time, and transport processes in a 10-km 2 urban watershed. The model, which operates over a 10-year period at a daily time scale, represents storage in a dynamic soil-moisture reservoir, quick-flow runoff from storm events, and slow-flow runoff that sustains streamflow in dry weather. The calibrated model accurately represented (a)the observed transition from a streamflow enrichment regime in cold months to a dilution regime in warmer months, (b) the observed tendency for late-summer concentrations to be higher after winters with heavy snowfall, and (c) a period-of-record downward trend in chloride concentration likely associated with a downward trend in annual snowfall. Estimated chloride inputs averaged 195 metric tons per year, while the average output was 270 metric tons per year. In contrast, estimated storage was only 107 metric tons. The estimated mean residence time in groundwater was 1.27 years. This short residence time indicates that efforts to reduce inputs will manifest as decreased concentrations in streamflow on a management-relevant time scale of several years. The coupled mass balance model yielded insights into internal watershed dynamics that would not be possible from simple input/output analysis; such models can be useful tools for gaining insight into small watershed hydrology and pollutant transport.

Virginia↗
USGS2025

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona↗
USGS1992

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee↗
USGS2024

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗
USGS2025

Invasion of perennial sagebrush steppe by shallow-rooted exotic cheatgrass reduces stable forms of soil carbon in a warmer but not cooler ecoregion

Soil organic carbon ('SOC') in drylands comprises nearly a third of the global SOC pool and has relatively rapid turnover and thus is a key driver of variability in the global carbon cycle. SOC is also a sensitive indicator of longer-term directional change and disturbance-responses of ecosystem C storage. Biome-scale disruption of the dryland carbon cycle by exotic annual grass invasions (mainly Bromus tectorum, 'Cheatgrass') threatens carbon storage and corresponding benefits to soil hydrology and nutrient retention. Past studies on cheatgrass impacts mainly focused on total C, and of the few that evaluated SOC, none compared the very different fractions of SOC, such as relatively unstable particulate organic carbon (POC) or relatively stable, mineral-associated organic carbon (MAOC). We measured SOC and its POC and MAOC constituents in the surface soils of sites that had sagebrush canopies but differed in whether their understories had been invaded by cheatgrass or not, in both warm and relatively colder ecoregions of the western USA. MAOC stocks were 36.1% less in the 0–10 cm depth and 46.1% less in the 10–20 cm depth in the cheatgrass-invaded stands compared to the uninvaded stands of the warmer Colorado Plateau, but not in the cooler and more carbon-rich Wyoming Basin ecoregion. In plots where cheatgrass increased SOC, it was via unstable POC. These findings indicate that cheatgrass effects on the distribution of soil carbon among POC and MAOC fractions may vary among ecoregions, and that cheatgrass can reduce forms of carbon that are otherwise considered stable and 'secure', i.e. sequestered.

Environmental Research Communications↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗