Geology topics

Browse USGS reports

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

At least 1,315 records · Page 73Linked to original sources

USGS2026

Environmental tradeoffs of urban stream restoration in Fairfax County, Virginia

Regulatory mandates to improve water quality and stream health have driven substantial investment in stream restoration. Most projects aim to improve channel-floodplain connectivity, reduce sediment erosion, and enhance habitat for aquatic organisms, yet few include adequate pre- and post-restoration monitoring to assess outcomes. Since 2007, Fairfax County, Virginia, and the U.S. Geological Survey have partnered to monitor and evaluate water-quality conditions in Flatlick Branch. In 2018, a 1.72-km reach of stream above the monitoring station was restored using a Natural Channel Design approach. This study applied the Stream Functions Pyramid (SFP) framework to evaluate restoration tradeoffs. Post-restoration, watershed hydrology remained largely unchanged, but channel modifications increased flow capacity, reduced velocity, and further disconnected the channel from the floodplain. Nutrient and sediment reductions exceeded expected amounts, but the removal of over 20 % of riparian tree canopy increased physicochemical variability and the frequency and magnitude of water temperature heatwaves. Post-restoration, state standards for low dissolved oxygen and elevated pH were exceeded 2.5 and 7.5 times more often, respectively. Gross primary production and ecosystem respiration increased and organic matter sources supporting metabolism shifted from allochthonous to autochthonous. Trends in several benthic macroinvertebrate metrics, which were improving prior to construction, have since plateaued or declined, and the fish assemblage shifted from a native minnow dominated community to non-native, warmwater tolerant taxa. This study highlights the need for comprehensive assessments of stream restoration and benefits of using the SFP to understand the consequences and possible tradeoffs of different ecosystem management decisions.

Virginia↗

USGS2026

A tool to monitor hydrologic conditions on tree islands in the Everglades

Tree islands are patchy upland forested habitats in Florida's Everglades that face degradation and disappearance due to altered hydrologic patterns. The U.S. Geological Survey coordinated with the Miccosukee Tribe of Indians of Florida and the Seminole Tribe of Florida to co-develop a decision-support tool based on tree-island hydrologic conditions. Everglades managers can use this tool to help with restoration planning and water operations decisions that affect tree-island conditions. After a series of organized workshops and meetings, a list of hydrologic metrics was selected as indicators of tree-island health, including hydroperiod, number of days since last dry, and maximum water depth at the head of the island. As a result, a web application tool, called ETree, has been developed and is publicly available online. This web application provides data on daily metrics for the current Everglades water year and annual summaries for past years, beginning in 2000.

Florida↗

USGS2026

Vulnerability of different Colorado Plateau land types to drivers of change

Public lands are challenged by a range of pressures—changing climate, increasing visitation, resource extraction—and their effects can span spatial scales, often crossing land management jurisdictional boundaries. Research approaches which explicitly span jurisdictions can support strategies to contend with regional pressures. We assess management-relevant drivers of change—aridification, livestock grazing, invasive species, surface disturbance, and fire—across a patchwork of land management units and agencies on the Colorado Plateau, focusing on southeastern Utah, USA. We use vulnerability analysis, first evaluating exposure to drivers across the landscape, then quantifying sensitivity to each driver across different land types, defined by mapped Ecological Site Groups, a system for classifying landscapes according to physical factors including climate, soils, and topographic setting. We address the questions: 1) how are drivers spatially distributed across the study region; and, 2) based on exposure and sensitivity, are certain land types more vulnerable to these drivers? We find that the study region has high exposure and sensitivity—and thus high vulnerability—to aridification and grazing, but low exposure and vulnerability to other drivers. Although more sensitive land types were not generally more exposed, identifying which areas are most sensitive can guide adaptive measures, like where new uses or disturbances would be least harmful and which areas could be prioritized for restoration. The method we demonstrate is a flexible tool for assessing landscape-scale impacts, is built on nationally available datasets, and can be tailored to different datasets and sensitivity metrics.

Utah↗

USGS2026

Satellite time series analysis to quantify changing climax ciénegas using a state and transition model approach

Ciénegas are rare wetlands in arid landscapes of the North American Southwest, historically providing critical ecological and hydrological functions but increasingly threatened by changing climate and land use pressures. This study quantifies changes in ciénega condition and floodplain dynamics using a state-and-transition model (STM) informed by expert knowledge and remote sensing. Key factors include woody plant encroachment, water availability, and soil aggradation. We mapped 31 ciénegas with high-resolution imagery and analyzed Landsat data (1985–2023) to assess vegetation health and moisture using the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Infrared Index (NDII). Results show substantial interannual variability in phenology, water stress, and soil moisture, with regional drying and elevation strongly influencing ciénega resilience. We classified ciénegas into three functional states—healthy, desiccated, and dormant—and mapped their 2023 condition. Trend analyses indicate most ciénegas exhibit greening despite drought, though localized variability underscores the need for site-specific management. None are in a stable climax (reference) state; rather, they transition among states in response to external drivers. Increasing woody plant cover and surface drying, likely linked to declining regional water tables, favor deep-rooted species over wetland grasses—a pattern mirrored in adjacent control plots. Spatially explicit analysis revealed intra-ciénega variability often masked by aggregated data, highlighting the importance of high-resolution monitoring. Seasonal and long-term trends provide context for understanding ciénega dynamics, including degradation and restoration pathways. This study emphasizes the importance of groundwater conservation and demonstrates how remote sensing supports long-term monitoring. The STM framework offers a practical tool for adaptive management to sustain freshwater resources in arid environments.

Arizona, New Mexico↗

USGS2026

Severity and certainty metric for kelp stressors: A tool for research, management, and conservation

Foundation species face multiple threats, and many are in decline worldwide. Kelps are no exception, with declines documented throughout their range. Stressors such as warming, eutrophication, and reduced water quality—often linked to human activities—are well recognized, yet their relative severity and certainty have not been systematically evaluated. Such a framework may aid scientific, regulatory, and conservation organizations in prioritizing actions and decision making. Here we build on a conceptual model of stressors to kelp in the Salish Sea by rating stressor severity and scientific certainty. We convened a panel of regional kelp experts to perform the rating exercise for two canopy taxa, Nereocystis luetkeana and Macrocystis tenuifolia , as well as a grouping for understory species, including both the sporophyte (adult) and gametophyte (juvenile) life stages for these algae. Our process produced a rich dataset which can be separated by stressor, taxa, and life stage allowing for broad applications. The panel identified temperature, low nutrients, and substrate change/benthic sedimentation as the most severe stressors to kelp, with relatively high certainty for most taxa. Certainty scores were relatively high for Nereocystis and Macrocystis sporophytes but were considerably lower for understory species and the gametophyte life stage, for most stressors, highlighting research gaps. Our data set could be applied for management prioritization, setting ecological baselines, and spatial modeling exercises. We encourage scientific, regulatory, and conservation organizations to consider both our data and our process when evaluating their needs. As foundation species continue to face threats, exercises such as this can provide scientists and managers with actionable indices that support strategies to stem degradation.

Ecological Indicators↗

USGS2026

The use of an analog sclerometer to estimate strength of hard substrates supporting coastal communities

In the marine environment, substrate differences such as ‘soft’ and ‘hard’ bottom result in different ecological communities. It is unclear, however, whether strength of hard bottom affects organism presence. To date, a limitation of many studies on organism-substrate interactions is the inability to measure rock strength in situ for marine substrates. Previous efforts involved measurements of rock compressive strength, which relates to the hardness of the interior rock and does not necessarily reflect how organisms interact with the rock surface. Here we present a method for estimating the surface strength, or friability, of marine rocks in the field. We used an analog sclerometer to scratch the surface of marine rocks at predetermined force levels across 38 intertidal sites. Friability ranged over two orders of magnitude (0.13–12.10 MPa) with significant differences among rock types: mudstone (mean = 0.56 MPa) had the lowest resistance, while granitic rocks (mean = 2.14 MPa) had the highest, consistent with theoretical expectations. Friability measurements indicated expected compressive rock strengths across rock types and captured differences within types (i.e., surface weathering). The sclerometer effectively measured surface rock strength on soft rocks, such as mudstone, and hard rocks, like granite. The ability to measure the strength of surfaces represents a new tool for benthic marine research, enabling investigations of substrate-organism interactions in future studies.

Estuarine, Coastal and Shelf Science↗

USGS2026

Cumulative effects analysis to inform public land management in the United States: Key characteristics and legal challenges

Considering potential cumulative effects of proposed actions is fundamental to environmental impact analysis. However, cumulative effects analyses historically are not robust, especially for site-specific decisions. We sought to identify opportunities to strengthen cumulative effects analysis in a large United States public land management agency, the Bureau of Land Management (BLM). We asked 1) how cumulative effects analyses were legally challenged, 2) how site-specific cumulative effects analyses aligned with policy and compared to the broader-scale analyses to which they tiered, and 3) whether characteristics of cumulative effects analyses varied with category of proposed action, type of resource, or agency office. We used thematic analysis to assess litigation and appeals case documents finalized from 2010 to 2020 and a set of document analysis questions to assess National Environmental Policy Act (NEPA) analyses for BLM decisions completed prior to 2020 in Alaska and Colorado. We found that legal challenges related to cumulative effects focused on absence of cumulative effects analysis. In NEPA analyses, cumulative effects were frequently considered, but elements recommended in policy, such as citations, methods, and scope, were rarely included. These elements were present more often in the broader analyses to which site-specific analyses tiered. Many elements of cumulative effects analyses varied by proposed action and BLM office, and analyses of potential cumulative effects on air quality were consistently more detailed than for other resources. Our results suggest that many problems that historically plagued cumulative effects analysis persist. Advances in methods, training, and guidance could strengthen the defensibility of NEPA analyses.

Alaska, Colorado↗

USGS2026

Snow simulations predict future changes in rain-on-snow events across the upper Gallatin River watershed, a Greater Yellowstone Ecosystem headwater system

Study region: The upper Gallatin River watershed, an alpine headwater system in the Greater Yellowstone Ecosystem, in Wyoming and Montana. Study focus: As global and regional air temperatures rise, mountain headwaters across the Greater Yellowstone Ecosystem (GYE) are projected to see more precipitation falling as rain. While the hydrologic effects of this snow-to-rain transition depends on a variety of factors, it can lead to an increased occurrence of rain-on-snow (RoS) events. To investigate these changes, we used high-resolution (30 m) SnowModel simulations of the upper Gallatin River watershed. Simulations were run for 2001-2013 using two scenarios: (1) historical meteorology as control and (2) pseudo global warming (PGW) where control air temperature and precipitation conditions were perturbed to represent mean end-of-century conditions under a high-emissions scenario. New hydrological insights for the region: SnowModel outputs show that changes in PGW precipitation and snow accumulation varied with elevation. Warmer air temperatures at low elevations (< 2,500 m) led to less snow accumulation and less precipitation falling as snow. Colder baseline air temperatures for elevations above 2,500 meters (m) resulted in minor reductions in winter snowfall fraction. For PGW simulations, spring (April-June) months were rainier, and elevations above 2,500 m experienced more RoS events. Snowpacks between 2,500-3,100 m generated more snowmelt during RoS events, which was reflected in the watershed average. More high-intensity melt events can affect aquatic habitat, water quality, and the accuracy of streamflow forecasts across the region.

Montana, Wyoming↗

USGS2026

VegET evapotranspiration for Africa: Continental-scale simulation, multi-product evaluation, and drought assessment

Study region Continental Africa, encompassing diverse climatic zones—tropical, arid, and temperate—and spanning major transboundary river basins such as the Nile, Niger, Congo, Volta, and Zambezi River Basins. The region exhibits pronounced hydroclimatic gradients and heterogeneous land use systems ranging from rainfed croplands and rangelands to dense tropical forests and irrigated schemes. Study focus Actual evapotranspiration (ETa) is a central component of the terrestrial water balance, governing the redistribution of water and energy between the land surface and the atmosphere. Accurate estimation of ETa at continental scale is critical for hydrological monitoring, water resource management, and climate adaptation, as well as for quantifying water, energy, and carbon fluxes that underpin sustainable development. In this study, we applied the agro-hydrologic VegET v2 model to simulate a new, high-resolution, continental-scale ETa dataset for Africa (2000–2021). The model results were benchmarked against four widely used remote sensing-based products—MODIS16 v6.1, SSEBop v6.1, WaPOR v3, and GLEAM v4.1a—across major climate zones, land use types, and River Basins, providing a comprehensive multi-product evaluation of evapotranspiration dynamics across the continent. New hydrological insights for the region Validation against eddy covariance flux tower observations at eight representative sites confirmed that VegET v2 accurately reproduces the seasonal dynamics of observed ETa, achieving a correlation (r) of 0.8 and an RMSE of 25 mm month⁻¹ —accuracy that is comparable to or higher than accuracies of satellite-based products MODIS16, SSEBop, and GLEAM. This study represents one of the first Africa-wide hydrological simulations of ETa, extending the VegET model beyond basin-scale applications. Intercomparisons reveal that VegET aligns closely with MODIS16, SSEBop, and GLEAM in humid and tropical regions (r = 0.80–0.90; RMSE < 20 mm month⁻¹), while greater discrepancies appear in arid and semi-arid zones, where WaPOR tends to overestimate ETa (RMSE ≥ 28 mm month⁻¹). Despite these differences, VegET effectively captures spatial and temporal ETa variability across rainfed croplands, forests, and savannas, supporting its utility in regional water balance assessments, water accounting, and drought monitoring. A key application of VegET v2 is the Evapotranspiration Deficit Index (ETDI), derived by integrating VegET-based ETa with potential evapotranspiration (PET) to quantify water stress. ETDI successfully captured major drought episodes across Africa, including persistent Sahelian and southern African dry spells, the 2020–2021 winter drought in the Maghreb, and the 2018–2019 austral summer drought in southern Africa, while identifying positive anomalies over central Africa indicative of recurrent wetness. These results underscore VegET’s capability as a hydrologically consistent, operational tool for continental ETa monitoring and drought assessment, offering support for basin-scale water balance studies, food security planning, and climate resilience across Africa’s diverse hydrological environments.

Journal of Hydrology: Regional Studies↗

USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗

USGS2026

From landslide susceptibility to risk assessment in the conterminous U.S.

Understanding the spatial distribution of landslide prone-areas and what consequences they may have is important for risk management and land-use planning. In the United States, although landslides occur in every state, a comprehensive landslide risk assessment is still missing. Existing efforts, such as the Federal Emergency Management Agency (FEMA)’s National Risk Index, rely on aggregated products and coarse cartographic units, limiting their geomorphological and practical accuracy. In this study, we present a methodological advance for landslide risk assessment across large areas with incomplete and sparse data. We apply our procedures to the conterminous United States by integrating geomorphologically meaningful partitions and spatial and temporal probability data-driven models. Landslide susceptibility is estimated using a Generalized Additive Mixed Model incorporating a bias capture/correction scheme to account for inventory inaccuracies (reference Area Under the Curve = 0.75). The exceedance probabilities of landslide occurrence are defined for three temporal scenarios (2, 5, and 10 year). Then, we explore the associated potential economic consequences for human settlements and agricultural areas. The findings indicate that the spatial variability of risk is primarily controlled by exposure rather than by susceptibility/hazard alone. The mean risk increases by ∼170% from the 2-year to the 10-year scenario. Beyond its quantitative outcomes, this study offers a blueprint for continental or sub-continental scale landslide risk assessments, demonstrating both the opportunities and current limitations.

Engineering Geology↗

USGS2026

Per- and polyfluoroalkyl substances and pesticides in black bears (Ursus americanus) and grizzly bears (Ursus arctos horribilis) from Montana's Cabinet-Yaak ecosystem: A baseline assessment of emerging and legacy contaminants

Per- and polyfluoroalkyl substances (PFAS) and pesticides are globally distributed contaminants that persist in terrestrial food webs, yet baseline data for large omnivores in the continental U.S. remain limited. This study quantified PFAS and pesticides in whole blood from free-ranging black bears ( Ursus americanus , n = 15) and grizzly bears ( Ursus arctos horribilis , n = 9) in the Cabinet-Yaak Ecosystem of northwestern Montana. Using high-resolution mass spectrometry-based methods (GC-MS/MS and LC-MS/MS), 52 PFAS and 186 pesticides were analyzed. PFAS were detected in all bears, with total PFAS concentrations ranging from 79.5 to 317 pg/mL. PFAS profiles were dominated by long-chain perfluoroalkyl carboxylic acids, particularly PFOA, PFDA, and PFUnDA, with minimal contribution from short-chain compounds. Species was the strongest ecological predictor of PFAS concentrations, although variability across age, sex, and capture location were observed. Only two pesticides, piperonyl butoxide and permethrin, were detected at quantifiable concentrations (2070–12,600 pg/mL), and no correlations were observed between pesticide and PFAS concentrations, indicating independent exposure pathways. The predominance of long-chain PFAS suggests diffuse environmental sources and bioaccumulation within terrestrial food webs. Although measured concentrations were generally low relative to other wildlife toxicological studies, the persistence and protein-binding properties of long-chain PFAS warrant continued monitoring, particularly for long-lived omnivores. These results establish baseline contaminant concentrations for two apex omnivores in a remote ecosystem and highlight the utility of large mammals as integrative sentinels of environmental contaminant exposure for regional monitoring and wildlife conservation. Continued monitoring in the Cabinet-Yaak and across western North America will be critical for assessing temporal trends, identifying new contaminant sources, and evaluating ecological health in bear populations.

Montana↗

USGS2026

Weathering the ride: Associations of heat, smoke, precipitation, and ozone with bus ridership in Colorado

Weather and climate variability increasingly shape urban travel behavior, yet the short-term temporal dynamics and contextual modifiers of weather-transit relationships remain poorly understood. We analyzed system-wide hourly bus ridership in the Denver metropolitan area from June 2022 through September 2023 using fixed-effects negative binomial distributed lag nonlinear models with lags up to 24 h. Nonlinear exposure-response and lag-response functions were specified for hourly Universal Thermal Climate Index (UTCI) and precipitation, with additional indicators for daily ozone exceedance and wildfire smoke. Models adjust for hour-of-day, day-of-week, and month-year fixed effects, with standard errors clustered by day. Stratified analyses assess heterogeneity by time period, season, fare policy, and shelter availability, with interaction evaluated using Wald tests. Cold thermal stress was associated with the largest, most persistent reductions in ridership, with cumulative declines of −15.7% (95% CI: −26.0%, −4.0%) over 0-24 h. Precipitation was associated with sharp but transient reductions concentrated within 3-6 h (−8.4% to −9.8%), with little evidence of longer-term displacement. Heat associations were weaker and context-dependent, with modest short-run increases, but net same-day declines. Associations varied by time of day and season, were attenuated during the free-fare period, and differed by shelter availability. In contrast, ozone exceedance and wildfire smoke exhibited limited and inconsistent associations, with measurable reductions primarily during morning commute during heavier smoke conditions. Overall, transit ridership is more strongly associated with short-term weather exposure, than with ambient air quality. Fare policy and stop-level infrastructure modify these associations, highlighting actionable strategies to enhance transit resilience under increasing climate variability.

Colorado↗

USGS2026

A catalogue of Do's and Don'ts in the modeling of environmental systems

Modeling plays a vital role in understanding and managing complex environmental systems, but its credibility and quality depend heavily on a comprehensive set of defensible model activities and practices, especially when the system of interest is plagued with uncertainties and conflicting stakeholder perspectives. This paper proposes a catalogue of Do's and Don'ts to guide modelers in addressing the many pertinent considerations through the whole modeling cycle. This practical tool provides advice on approaching modeling effectively through adhering to good modeling practice. It emphasizes model choices that align with the model purpose and context, and the justification and documentation of modeling decisions and assumptions. Managing uncertainty is a core consideration. The identification, assessment and reporting of these uncertainties is important across the entire modeling process, which spans problem framing, technical design, implementation and application phases. Such good practices are critical for transparency and reliability of the modeling.

Environmental Modelling and Software↗

USGS2026

Multiple machine-learning estimation of groundwater levels and trends for the regional Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer provides irrigation, public, and domestic water supplies across the south-central United States. Declining groundwater levels require improved characterization of changing conditions. Traditional potentiometric-surface mapping does not use all available water-level data or quantify uncertainty. To address these limitations, we developed a data-driven multiple machine-learning (MML) framework delivered through two open-source R packages. The covMRVAgen1 software assembles covariates to 155,960 monthly groundwater levels from 57,695 wells; the mmlMRVAgen1 software trains Cubist and Random Forest models, blends them, and makes 1-kilometer gridded predictions of monthly potentiometric surfaces for the period January 1980–December 2022. The MML approach provides a methodological foundation for region-scale spatiotemporal groundwater prediction and uncertainty quantification, generating 90-percent prediction limits with appropriate empirical coverage. Model performance is acceptable, with a root-mean-square error of about 4.2 feet, standard deviation of 24.82 feet, and a normalized Nash–Sutcliffe efficiency of 0.973.

Arkansas, Illinois, Louisiana, Mississippi, Missou↗

USGS2026

Patterns of rift basin development and the fidelity of the subsidence record: Insights through Bayesian modeling of rapid tectonic subsidence in a Rio Grande rift basin, Socorro, NM, U.S.A

Characterizing tectonic subsidence rates within depositional sequences provides direct insight into the driving mechanism(s) of accommodation in a basin. However, the temporal resolution of this record is often stymied by a lack of high-precision and high-resolution ages, which enable a more complete description of basin subsidence drivers. We explore the effect of high-precision and high-resolution ages in modeling accommodation for the Miocene La Jencia Basin of the central Rio Grande rift (RGR) and interpret driving mechanism behavior from these models (e.g., lithospheric thinning). We present a new geochronologic dataset of both laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS) and chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) data and use these ages in Bayesian accommodation modeling. Models constrained by high-precision and high-accuracy TIMS ages yield peak tectonic subsidence rates exceeding 220 m/Myr, and an average Miocene subsidence of ∼120 m/Myr. While timing and magnitude vary, all models suggest two pulses of rapid Miocene tectonic subsidence, which we interpret to reflect basin-bounding fault movement. Prior to peak subsidence, there was an initial period of fault linkage and organization that occurred over <1–3 Myr that produced the basin-bounding La Jencia-Cerro Colorado fault zone. A comparison of published tectonic subsidence rates to those modeled here shows that while tectonic subsidence during continental rifting is highly variable, the La Jencia Basin rates appear relatively high. However, the significant difference between peak and average La Jencia Basin rates modeled here highlights the potential for underestimation of many records of tectonic subsidence due to a lack of high-precision and high-resolution age constraints. Furthermore, age data and modeling results presented here document fault movement and consequent rates of tectonic subsidence that lower-resolution data would not, providing a high-fidelity case study of continental rift basin development.

New Mexico↗

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↗

USGS2026

Comprehensive characterization of Mauna Loa basalt rheology based on the 2022 eruption

Rheological data are crucial for the understanding of storage, transport, eruption, and emplacement of magma and lava. While viscosity estimates for lavas from Mauna Loa, the largest active volcano on Earth, exist, no direct rheological measurements have been done to date. This is especially surprising given the role Mauna Loa lavas and resulting landforms play in the interpretation of volcanic processes across our solar system and beyond. We present the first complete rheological characterization of Mauna Loa lava, based on the 2022 eruption. We constrain the melt viscosity and map the lava's rheology during crystallization at isothermal subliquidus conditions and at cooling at rates between 0.25 and 3.00 °C/minute. These experiments were performed at oxygen fugacities relevant to storage and eruption conditions at Mauna Loa (log fO 2 = -8.7). We integrate textural and thermal characteristics obtained from the experiments with those from water quenched natural samples (i.e. sampled from actively flowing lava and quenched within seconds to preserve the textural state of the lava while flowing) to reconstruct the evolution of lava rheology during emplacement. This integration suggests that the Mauna Loa lavas were emplaced at viscosities between 10 1.5 to 10 4.5 Pa s over the temperature range of 1150 – 1090 °C before entering a region of rapid viscosity increase and effective solidification caused by groundmass crystallization. We provide a detailed comparison of our data to rheological reconstructions of the 1984 eruption. The resulting rheological map can help guide physical property-based magma ascent and lava emplacement models, as well as the interpretation of flow morphologies and the conditions of formation of basaltic landforms on other planets.

Hawaii↗