USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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.

1,687 records · Page 68Linked to original sources
USGS2026

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗
USGS2026

Streamflow, base-flow, and precipitation trends and simulated effects of groundwater withdrawals from the North Fork Red River aquifer on base flows upgradient from Lake Altus, western Oklahoma, 1980–2022

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used five scenarios created from a previously published numerical groundwater-flow model (1980–2013) and historical streamflow records (1980–2022) to investigate the relation between groundwater withdrawals from the North Fork Red River aquifer and inflows to Lake Altus from the North Fork Red River in western Oklahoma. The five scenarios were (1) a scaled-equal-proportionate-share (EPS) groundwater-withdrawal scenario, (2) a study-area-scaled-reported groundwater-withdrawal scenario, (3) a zonal-scaled-reported groundwater-withdrawal scenario, (4) a historical drought-threshold scenario, and (5) a base-flow and evapotranspiration depletion scenario. For the scaled-EPS groundwater-withdrawal scenario, EPS groundwater withdrawals were often much higher than reported groundwater withdrawals and greatly decreased base flows for most scale factors. For the study-area-scaled-reported groundwater-withdrawal scenario, base flows were reduced more but by smaller percentages during wet periods than during dry periods when scaling simulated reported groundwater withdrawals. For the zonal-scaled-reported groundwater-withdrawal scenario, scaling simulated reported groundwater withdrawals within selected zones with more groundwater withdrawals did not always affect base flows more than scaling reported groundwater withdrawals within zones with less groundwater withdrawals. For the historical drought-threshold scenario, curtailing groundwater withdrawals at the drought thresholds increased annual base flows to Lake Altus by about 1,169 to 3,665 acre-feet. For the base-flow and evapotranspiration depletion scenario, the distance between a groundwater well and a stream was a major factor affecting base flow to the North Fork Red River when increasing groundwater withdrawals; however, spatially variable hydrologic properties and saturated-zone evapotranspiration could also affect the relation between base flows and groundwater withdrawals.

Oklahoma↗
USGS2026

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States↗
USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗
USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

New Mexico↗
USGS2026

The Sedimentary Geochemistry and Paleoenvironments Project Phase 2 data release: An open data resource for the study of Earth's environmental history

Geochemical data from sedimentary rocks are the primary source of information regarding Earth's surface evolution through time, including its air and water envelopes and interactions with life and deep Earth processes. The Sedimentary Geochemistry and Paleoenvironments Project (SGP) is a scientific consortium centered around open data and community-driven development of cyberinfrastructure tools and resources for sedimentary geochemistry and Earth history. Here we describe the SGP Phase 2 data release, which focused on incorporating Paleoproterozoic and Mesoproterozoic (2500–1000 million years ago) data and better accommodating carbonate data. This data release was built through the involvement of >200 researchers worldwide in academia, government, and industry, and provides the largest available public data resource for our user community in the academic fields of geochemistry, sedimentology, tectonics, paleontology, Earth history, and paleoclimate, as well as the petroleum and minerals industries. The dataset now encompasses 126,006 samples and 4,132,371 geochemical analyses. In addition to direct entry by SGP Team Members, we have ingested and incorporated datasets from the Geoscience Australia OZCHEM database, the Alberta Geological Survey, and the Deep-Time Marine Sedimentary Element Database (DM-SED) compilation. This paper details sampling in the Phase 2 dataset with respect to age, geography, lithology, and other geological characteristics, documents access via our search website and API, discusses possible issues and/or biases in the dataset that could impact analyses, describes plans for governance and stewardship of data from Indigenous lands, and serves as the citable reference paper for the data release.

Chemical Geology↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2026

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California↗
USGS2025

Synergy between geology and geophysics in graphite mineral resource assessment

Graphite is designated as a critical mineral by the U.S. Government due to its essential role in modern technology and its vulnerability to supply chain disruption. To evaluate domestic graphite resources, the U.S. Geological Survey (USGS) conducted two airborne electromagnetic (AEM) surveys as part of the Earth Mapping Resources Initiative (Earth MRI) over prospective areas in the southern Appalachians of Alabama and the Seward Peninsula of Alaska. These surveys aimed to refine areas of interest by revealing subsurface electrical resistivity patterns indicative of graphite-bearing rocks. The results of the AEM surveys were integrated with existing geological knowledge and supplemented by field observations to assess the likelihood of undiscovered resources in these regions. This multidisciplinary approach combined geophysical techniques with traditional geological methods, enhancing the understanding of graphite potential in areas that may have been previously overlooked. The findings from this initiative are expected to contribute significantly to the ongoing efforts to secure a reliable domestic supply of graphite, thereby reducing import reliance and supporting national interests in mineral resource management.

Conference Paper↗
USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗
USGS2024

Wave propagation in layered soil deposits

Recent advances in the general theory of viscoelastic waves and rays in layered media provide a rigorous mathematical framework for site-specific, soil-response models used for earthquake resistant design. The advances provide general closed-form anelastic solutions for the classic problems of the response of a stack of soil layers to S and P waves, ray theory for reflected and refracted waves, Rayleigh- and Love-Type surface waves, and head waves. These general solutions valid for anelastic media regardless of the amount of material damping yield new insights regarding the characteristics of seismic waves and their ray paths that are not provided by conventional models. They provide corresponding numerical ground-response models and ray-tracing computation algorithms that account for changes in velocity and attenuation of anelastic waves associated with changes in inhomogeneity of the waves induced by anelastic soil and soil-rock boundaries. These anelastic effects manifest themselves as variations in amplitude response, amplitude attenuation, raypath location, and travel time as observed at the Earth’s surface. Implications of these anelastic effects for soil-response models used for earthquake resistant design are provided herein.

Conference Paper↗
USGS2023

Stream sediment geochemistry in mineral exploration: A review of fine-fraction, clay-fraction, bulk leach gold, heavy mineral concentrate and indicator mineral chemistry

Stream sediment surveys support early-stage reconnaissance mineral exploration and regional assessment programmes, enhanced by recent improvements in analytical method detection limits, continuously improving mineral chemistry, and new approaches to the interpretation of geochemical data. Sediment surveys may be used to predict catchment basin lithology, mineralization type based on pathfinder geochemistry, and geological features based on indicator mineral chemistry. Sediment surveys that target a finer-fraction sediment sample led to the discovery of the La Colosa gold deposit, Colombia. The Batu Hijau porphyry Cu–Au deposit in Indonesia was discovered based on an anomalous clay-sized fraction sample 12 km downstream. In an arid region with poorly developed drainages and minor topographic relief, the Ag-base-metal Navidad District in Argentina was discovered with clay-fraction sediment geochemistry. Heavy mineral concentrate (HMC) sediment surveys that include mineral chemistry determinations have led to global diamond discoveries. HMC surveys contributed to discovery of the Ring of Fire Ni–Cu–PGE and chromite district, Ontario, Canada. Discoveries and geochemical mapping can assist advancement of the application of stream sediment geochemistry in those global areas for which lithologies and deposits are exposed. Stream sediment surveys continue to be one of the most cost-effective geochemical methods for covering large areas for mineral exploration.

Geochemistry: Exploration, Environment, Analysis↗
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2025

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming↗
USGS2024

Monitoring marine eruptions

Introduction Submarine volcanoes produce much of the same seismicity and eruptive activity as subaerial volcanoes and can pose hazards to society. Although they can be monitored with similar techniques and methods as described in other chapters of this volume, their submerged location brings unique challenges. This chapter addresses these challenges and provides recommendations for monitoring volcanoes fully or partly in marine environments to meet the capabilities described in other chapters of this volume. The United States and its territories host dozens of submarine volcanoes with most (around 60) in the Commonwealth of the Northern Mariana Islands. Approximately 20 of the Northern Mariana Islands submarine volcanoes are known to be hydrothermally active, and 10 have confirmed eruptions since the 1950s (for example, Baker and others, 2008; Tepp and others, 2019a). Nine of those volcanoes were considered by the National Volcanic Threat Assessment (Ewert and others, 2018) to have a combination of eruptive type and summit depth that poses a higher risk of hazardous eruptions, although only one was listed as a moderate (level 3) threat. Other notable submarine volcanoes of interest to the United States that have historically erupted are Axial Seamount off the Washington State coast, Kamaʻehuakanaloa in Hawaiʻi, and Vailuluʻu seamount in American Samoa. All of these, however, have a low risk of hazards because of their depth (greater than 600 meters below sea level) and eruptive type and so are not included in the National Volcanic Threat Assessment. In addition to submarine volcanoes, the submerged flanks of island volcanoes can also be a source of hazardous submarine eruptions—for example, the 1877 eruption of Mauna Loa, Hawai‘i, in Kealakekua Bay (Wanless and others, 2006). The most notable submarine eruption in recent times was the 2022 eruption of Hunga Tonga–Hunga Haʻapai in Tonga, which was one of the largest eruptions on Earth in the past 100 years. It created a massive volcanic plume, unprecedented shock waves, and far-reaching tsunami (Lynett and others, 2022). Other recent submarine eruptions in the Pacific Ocean Basin have produced subaerial plumes that reached aircraft heights (Carey and others, 2014) and large pumice rafts that can affect marine traffic and harbors (for example, Jutzeler and others, 2014; Kornei, 2019). These examples illustrate the potential hazards of major submarine eruptions. Yet, submarine volcanoes are largely unmonitored, and many eruptions occur that are unnoticed or only identified hours or days afterward. Within U.S. territory, submarine volcanoes in the Northern Mariana Islands have been known to produce eruptive activity that can affect society. Reports from fishermen and other marine vessels in the Northern Mariana Islands have noted underwater explosions, sea-surface discoloration, and bubbling water, all of which are known to be signs of submarine volcanic activity. South Sarigan seamount, located about 160 kilometers (km) north of Saipan, erupted in 2010 from greater than 150 meters below the sea surface, resulting in a gas and ash plume that reached more than 11.9 km into the atmosphere (for example, Searcy, 2013; Embley and others, 2014), high enough to affect international air traffic. Precursory and co-eruptive seismicity was detected on the regional Northern Mariana Islands seismic network (Searcy, 2013) and on global monitoring instruments (Green and others, 2013). Monitoring of submarine volcanoes is best accomplished with marine-based instrumentation, which is also useful for monitoring small island volcanoes that may not have the land area necessary for comprehensive subaerial monitoring. The primary marine-based instrumentation used for submarine volcanoes includes ocean-bottom pressure sensors to assess sea-floor deformation, ocean-bottom seismometers (OBSs) to detect seismicity, and both moored and ocean-bottom hydrophones to detect submarine explosions. Other sensors offer important monitoring data, such as turbidity, temperature, and chemistry of hydrothermal emissions. Marine-based instruments are typically deployed in campaign-style networks with no real-time telemetry owing to cost considerations and technical limitations. However, when necessary, marine instruments can be operated in real time using cables to transmit data to land-based facilities; other technologies for this purpose are in use or in development, such as acoustic transmission from the instrument to a moored buoy (Matsumoto and others, 2016) and a winch-based system with a satellite antenna that is part of the instrument mooring (Matsumoto and others, 2019). Emerging technologies for marine-based monitoring may be considered as part of a future monitoring plan. These technologies include ocean gliders and floats with on-board hydrophones that have been used to record earthquakes and submarine eruptions (for example, Matsumoto and others, 2013; Sukhovich and others, 2015) and fiber-optic cables that have been used as strainmeters to detect earthquakes (for example, Marra and others, 2018; Lindsey and others, 2019). Land-based instruments and satellites can also provide some capability for monitoring submarine volcanoes, but they provide more limited observations than marine-based instrumentation.

Scientific Investigations Report↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2024

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↗
USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗