USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Microbiology”

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,596 records · Page 65Linked to original sources
USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗
USGS2023

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper↗
USGS2026

Deep subsurface organic-rich shale supports abundant, diverse, and novel fungi

As Earth’s principal reservoir of organic carbon and microbial biomass, the deep subsurface hosts microorganisms capable of mobilizing this once-sequestered carbon. Contrary to standard assumptions of eukaryotic scarcity, this study documents abundant fungal communities, ranging from 4.2 × 10 3 to 6.8 × 10 3 fungal cells mL −1 , across a methane-producing organic-rich shale 247–556 meters below the surface. Although fungal:bacterial cell ratios ranged from 1:7028 to 1:713, application of biomass conversion factors developed for oceanic systems yielded a median fungal:bacterial biomass ratio of 1:4.7. 16S rRNA gene amplicons revealed bacterial and archaeal communities mirroring those found in well-characterized extremophilic, carbon-degrading environments, while sequencing of 18S rRNA gene and ITS rRNA spacer amplicons collectively identified a eukaryotic hotspot with 689 fungal OTUs across six phyla. The dominant fungal classes, Agaricomycetes and Dothideomycetes, are well-established degraders of recalcitrant carbon compounds at the surface, suggesting they may similarly contribute to organic matter degradation and ecosystem maintenance in the subsurface. Cultivation and isolation efforts yielded 205 fungal strains, including 13 candidate novel taxa, underscoring the deep subsurface as an underexplored eukaryotic habitat. Stable carbon isotopes indicate methane is predominantly generated via microbial conversion of the fossil carbon, while water isotopes suggest in situ geochemical conditions have been relatively stable since the Late Pleistocene, with subglacial recharge as a plausible mechanism for microbial introduction. Collectively, these findings suggest that fungi are underrecognized contributors to organic matter transformation and functional diversity in the deep biosphere, revealing a critical gap in our understanding of deep subsurface ecosystem processes.

Indiana, Michigan, Ohio↗
USGS2025

Stratigraphic notes—Volume 2, 2025

This is the second volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 .

Professional Paper↗
USGS2026

Stratigraphic notes—Volume 3, 2026

This is the third volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 .

Professional Paper↗
USGS2026

Production of mineral commodities and geospatial map of the mineral industries and related infrastructure of China

As part of the U.S. Geological Survey’s (USGS) mission to distribute global mineral information and analyze supply chains, this study provides a comprehensive review of the global significance of China’s mineral production and capacity in 2023. Of 77 mineral commodities in the USGS dataset, China produced 74 and was the world’s first-ranked producer for 39 of the 74. Compared to the high share of global mineral production, including up to 98 percent of global gallium production, the country’s share of global mineral reserves was relatively small, ranging from 20 percent (zinc ore) to 52 percent (tungsten ore). China’s imports of metal ores, slag, and ash accounted for 64 percent of global imports of such commodities by value. The country’s exports of base metals and articles of base metal accounted for 17 percent of the global exports. To help nongeographic information system users assess the spatial distribution of mineral mines, processing facilities, and ports for trades in China, this study created a geospatial (also called “georeferenced”) portable document format (GeoPDF) map. In addition, the GeoPDF contains mineral resource tracts (such as antimony, copper, potash, coal, and oil and gas), exploration sites, and energy infrastructure based on the preexisting USGS data.

Open-File Report↗
USGS1997

An optimized network for phosphorus load monitoring for Lake Okeechobee, Florida

Phosphorus load data were evaluated for Lake Okeechobee, Florida, for water years 1982 through 1991. Standard errors for load estimates were computed from available phosphorus concentration and daily discharge data. Components of error were associated with uncertainty in concentration and discharge data and were calculated for existing conditions and for 6 alternative load-monitoring scenarios for each of 48 distinct inflows. Benefit-cost ratios were computed for each alternative monitoring scenario at each site by dividing estimated reductions in load uncertainty by the 5-year average costs of each scenario in 1992 dollars. Absolute and marginal benefit-cost ratios were compared in an iterative optimization scheme to determine the most cost-effective combination of discharge and concentration monitoring scenarios for the lake. If the current (1992) discharge-monitoring network around the lake is maintained, the water-quality sampling at each inflow site twice each year is continued, and the nature of loading remains the same, the standard error of computed mean-annual load is estimated at about 98 metric tons per year compared to an absolute loading rate (inflows and outflows) of 530 metric tons per year. This produces a relative uncertainty of nearly 20 percent. The standard error in load can be reduced to about 20 metric tons per year (4 percent) by adopting an optimized set of monitoring alternatives at a cost of an additional $200,000 per year. The final optimized network prescribes changes to improve both concentration and discharge monitoring. These changes include the addition of intensive sampling with automatic samplers at 11 sites, the initiation of event-based sampling by observers at another 5 sites, the continuation of periodic sampling 12 times per year at 1 site, the installation of acoustic velocity meters to improve discharge gaging at 9 sites, and the improvement of a discharge rating at 1 site.

Florida↗
USGS2026

Coastal barrier resilience and resistance: Analysis and metrics for characterizing coastal state

Barrier islands are shaped by a variety of short- and long-term environmental processes such as storms and relative sea-level rise. These islands, found along the estuarine-marine interface, provide ecosystem services including storm surge and wave attenuation, erosion protection to inland marshes, habitat for fish and wildlife, and recreation. Natural resource managers require actionable information on how barrier island resilience and resistance changes over time to understand how an island’s current state relates to past conditions and to inform restoration prioritization and implementation. The U.S. Geological Survey and The Water Institute collaborated on a study to develop indicators of resilience and resistance for barrier islands in Louisiana. Here, resilience captures island persistence on yearly to decadal time scales, and resistance captures persistence on event time scales of days to weeks. The indicators fall in two categories: Tier 1 Screening Metrics, that can be readily calculated from available data, are easily interpretable as an evaluation of barrier condition, and provide a high-level snapshot of overall barrier resilience and resistance; and Tier 2 Analysis Metrics, which are detailed metrics that required specialized analysis and interpretation and are more applicable to answering specific questions managers may have about barrier state. The research team derived Tier 1 resilience indicators from subaerial land and vegetation cover calculated from publicly available maps and products based on satellite imagery. By benchmarking the total land and vegetation extent against their respective historical maxima, this metric provides a snapshot of an island’s current state in the context of its long-term trajectory. The research team developed Tier 1 resistance indicators based on subaerial island configuration and water level recurrence as a proxy for evaluating island resistance to storms, which are the primary driver of short-term change. These Tier 1 metrics can be analyzed over time to provide a high-level assessment of how an island’s resistance decreases because of elevation loss or sea-level rise or increases due to restoration or natural recovery. The research team developed Tier 2 resilience and resistance indicators and associated analyses to provide detailed information for specific time periods or applications (e.g., wildlife management). These metrics include habitat coverage from high-resolution maps, which show composition changes over time to capture the evolving resilience of specific habitat types; high tide flooding analysis, which evaluate island area relative to specified flooding thresholds to characterize resistance in the short-term or, if analyzed over time, indicate changes in resilience; and hypsometric curve analysis, which allows managers to evaluate island area changes above their own elevation benchmarks of interest and similarly characterize resistance in the short-term or indicate changes in resilience if assessed over time. The research team calculated Tier 1 metrics of the barrier islands and headlands along the coast of Louisiana for the period of 1984 through 2021 and Tier 2 metrics for select times during that period depending on data available and quality. The results were captured in a report card for each barrier, which also includes an overview of the metrics and their interpretation; a restoration and storm history; and Tier 1 and Tier 2 metric analysis, including benchmarking against coastwide and regional values as well as to an island’s pre-restoration trajectory. These report cards provide a readily digestible synthesis of barrier condition and trajectory that coastal managers can use to support restoration prioritization and other decisions.

Louisiana↗
USGS2026

Long-term intermittent connection between the western Snake River Plain and Columbia basin: A two-phased incision history of Hells Canyon

For more than a century, researchers have debated the route of the Snake River across the northern Cordillera and U.S. Pacific Northwest, including the associated incision history of Hells Canyon. Here, we use detrital zircon U-Pb provenance analysis of Miocene strata upstream and downstream of Hells Canyon to constrain the evolution in fluvial pathways over time. Downstream of Hells Canyon, we show that the ca. 10−6 Ma Clarkston Heights gravel was dominantly sourced from the nearby Salmon and Clearwater rivers but also with substantial input from drainages that now flow into the western Snake River Plain. These same western Snake River Plain drainages were contributing to time-equivalent Lake Idaho strata, upstream of Hells Canyon, indicating a Hells Canyon fluvial connection between the western Snake River Plain and Columbia basin before 6 Ma. Results from Lake Idaho strata suggest a two-phased lacustrine history in the western Snake River Plain: Before ca. 6.0 Ma, lake strata were derived from local tributaries, suggesting that the western Snake River Plain was isolated from the modern upper Snake River system. Between ca. 4.3 Ma and 2.2 Ma, the source area expanded to include eastern tributaries in conjunction with migration of high-standing topography of the Yellowstone hotspot. Our results challenge the long-held hypothesis of “capture” of the Snake River through Hells Canyon. Instead, we infer that Hells Canyon was a long-established route for outflow of lakes occupying the western Snake River Plain, possibly intermittently, during a ca. 10−2 Ma phase of slow regional incision, followed by rapid incision in Hells Canyon starting ca. 2 Ma in conjunction with erosion of the topographic barrier impounding Lake Idaho.

Idaho, Oregon↗
USGS2025

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2024

SUTRA— A code for simulation of saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of the version 4.0 enhancements—Freeze-thaw capability, saturation and relative-permeability relations, spatially varying properties, and enhanced budget and velocity outputs

Version 4.0 of the Saturated-Unsaturated Transport (SUTRA) software code provides the capability to simulate the freezing and thawing of groundwater during energy transport simulations under saturated and unsaturated conditions. In addition to the types of hydrogeologic processes that SUTRA has been able to simulate in the past, this version can be used to study the effects of the freeze-thaw process on the flow and energy dynamics of hydrogeologic systems. The freeze-thaw simulation capability accounts for the latent heat of fusion and allows thermal property values to vary with changing total-water saturation, liquid-water saturation, and ice saturation. It allows the effective permeability of the porous medium to change as a result of freezing and thawing. This version also provides several user-selectable relations for the dependence of total-water saturation on fluid pressure, the dependence of liquid-water saturation on temperature during freezing and thawing, and the dependence of relative permeability on liquid saturation, as well as three user-selectable formulae for defining the bulk thermal conductivity of a mixture of solid grains, liquid water, ice, and air. For unsaturated simulations without freezing, the selectable total-water saturation relations eliminate the need for the user to program these and their associated relative-permeability functions, as had been required in previous SUTRA versions. Optional nonlinear dependence of fluid density on temperature, which covers the range from supercooled (about −50 degrees Celsius) to superheated (about 400 degrees Celsius), is also provided. Additionally, this version makes it possible to spatially vary parameters that, in previous versions of SUTRA, were required to be spatially uniform: solid-matrix properties, adsorption parameters, and parameters for production of solute mass or energy. Spatial variation is also allowed for the newly included freeze-thaw process parameters. Additional enhancements provide (1) output of water-mass and energy budgets that include values of all component terms in the governing balance equations, and (2) output of Darcy velocities (fluid fluxes), in addition to the velocity output provided by previous SUTRA versions. These enhanced outputs allow fuller interpretation of simulation results, especially for freeze-thaw phenomena. The set of processes simulated by this version of SUTRA are useful for studying a wide range of hydrogeologic system types, conditions, and questions. For cryohydrogeologic simulations, however, this version of the code is limited in that (1) it does not simulate thermomechanical effects of freeze-thaw, (2) pressure changes due to water density change during freezing are neglected, (3) ice saturation cannot exceed the initial porosity of the simulated medium, and (4) cryosuction, the migration of liquid water toward freezing fronts, is neglected. Furthermore, this version does not account for air flow or for water vaporization and sublimation under unsaturated conditions.

Techniques and Methods↗
USGS2025

Lunar grid systems, coordinate systems, and map projections for the Artemis missions and lunar surface navigation

Foreward This document contains design specifications of a navigational standard for the Moon, including a Lunar Transverse Mercator system, a Lunar Polar Stereographic system, a Lunar Grid Reference System, and a unique coordinate structure, Artemis Condensed Coordinates, for Artemis mission navigation and lunar surface science. The National Aeronautics and Space Administration (NASA) Artemis campaign seeks to place humans on the Moon for the first time since the Apollo missions. Early Artemis missions are heavily focused on the lunar south pole, which promises to return valuable data on the Moon’s geologic record, amongst other mission objectives. Coordinate systems in use today for the lunar south pole provides crew members on the surface neither an efficient nor intuitive means to communicate their position and orientation. A novel grid coordinate system, the Lunar Grid Reference System, is proposed to address these concerns for use in real-time extravehicular activity operations on the lunar surface. The many stakeholders involved in the Artemis missions will need a common system to communicate position and orientation while astronauts are operating on the lunar surface. To that end, Artemis crew members will need that system to be efficient and intuitive to promote efficient extravehicular activity timelines and reduce confusion. In the context of this document, these characteristics are addressed on the design of lunar coordinate systems: Efficient.—The number of characters required to communicate a location within a desired precision level in both local and global contexts, and how many steps are required for a recipient or sender to interpret a location. Intuitive.—How well the system aligns with human perceptual abilities, and whether the system yields distances that have the same relationship to actual lunar surface distance in all directions from the point where a person is located. Technological systems are currently being investigated to supplement the crew members’ ability to locate and orient themselves and other assets on the lunar surface; however, it is unlikely that those systems will be fully operational for the first few landed missions. Even with future positional aids, crew members will still need an efficient and intuitive means to communicate position and orientation. In addition, if technological systems fail, the crews will require land navigation skills and have maps available, thus providing further motivation for a crew-centric coordinate system. The contents of this U.S. Geological Survey (USGS) document detail a comprehensive framework for standardizing lunar crewed surface navigation within NASA and outlines the protocols, methods, and designs necessary for achieving consistency and interoperability across relevant space mission teams and lunar surface navigators. Key components of this document include designs of map projections, projected coordinate reference systems (Lunar Transverse Mercator and Lunar Polar Stereographic systems), and a grid system (Lunar Grid Reference System and Artemis Condensed Coordinates) for the Moon. The work proposed in this document seeks to accomplish something similar to the National Geospatial-Intelligence Agency (NGA) document SIG 0012 (NGA, 2014a), but for using grid systems for the Moon. This report incorporates initial feedback and input from NASA’s Artemis Geospatial Data Team, NASA’s Flight Operations Directorate, National Geodetic Survey, USGS Astrogeology Science Center, and NGA and is intended to serve as a resource for all involved with the Artemis missions, as well as for engineers designing and operating lunar infrastructure.

Techniques and Methods↗
USGS2026

The mineral chemistry networks of tin and tungsten reflect metallogenic eras of the Mesozoic

Continental remobilization is a crucial driver of metallogenesis and the formation of ore deposits. Some of the world’s largest mineral deposits of the economically valuable elements tin (Sn), tungsten (W), gold (Au), copper (Cu), lead (Pb), and zinc (Zn) formed during the Mesozoic Era. Additionally, the chemistry and distribution of the elements Sn and W have been investigated in previous studies to understand planetary formation and differentiation processes. These two elements are largely co-located during certain South China Mesozoic metallogenic events but are not co-located during other time periods in the same regions. Here, we investigated the mineral chemistry network similarities and dissimilarities of Sn and W to understand their mineral formation and distribution during the Mesozoic Era and throughout Earth history. Mineral chemistry network community detection analysis and electronegativity associations among mineral constituent elements of Sn minerals and W minerals indicate that the elements have similar chemistry among their oxide minerals. However, Sn forms a much wider range of minerals that also contain S compared to W, which occurs in a limited number of S-containing minerals. The divergent constituent element interactions among S-containing Sn minerals and W minerals reflect the redox sensitivity and importance of oxygen (O) fugacity in Sn mineral formation. Conversely, extensive W mineral deposits are known to form at both high and low O fugacities. The similarities and differences between the mineral chemistry networks of Sn and W reflect the mineral distribution of the two elements in the Sn-W mineralization event from 160 to 139 Ma vs. the Sn–uranium (U) mineralization event from 125 to 98 million years ago (Ma). The mineral chemistry and distribution of Mesozoic Sn and W deposits illustrate the contrasting importance of redox and O fugacity on the mineral formation of different elements, and the dynamic crustal evolution that took place during this period of Earth history.

Geosciences↗
USGS2025

PEST++IES how many iterations & realizations, finding the point of diminishing returns

PEST++IES (White 2018; White et al. 2020) is widely used in the groundwater modeling community for its ability to perform computationally efficient history matching and uncertainty analysis in a highly parameterized context. One primary advantage of using an iterative ensemble smoother is that the number of model runs required per iteration depends on the number of realizations in an ensemble, not the number of parameters in each realization. However, this raises the question: what is the optimal number of realizations and iterations to use for any one model before the point of diminishing returns? Using a modified version of the Freyberg model (Freyberg 1988; Hunt et al. 2020), different parameter and observation scenarios were evaluated for four iterations and ensembles of 10, 25, 50, 100, 250, 500, 1000, and 2000 realizations. To match observations, PEST++IES altered hydraulic conductivity ( k ), both globally across the model and locally at three different pilot point densities, as well as global recharge (via a single multiplier), global river conductance, and individual well flow rates. Risk-based well capture zone results (Fienen et al. 2022a) and estimated hydraulic conductivity fields from each scenario were quantitatively and qualitatively compared against the “truth” model and its outputs. Across the cases examined, ensemble sizes of 100 to 250 realizations and two PEST++IES iterations were generally sufficient to achieve good results.

Groundwater↗
USGS2025

Risk implications of Poisson assumptions and declustering inferred from a fully time-dependent earthquake forecast

We use the Third Uniform California Earthquake Rupture Forecast Epidemic Type Aftershock Sequence model, which is fully time-dependent in terms of including spatiotemporal clustering, to evaluate the effects of the Poisson assumption and declustering algorithms on statewide loss exceedance curves. The model is simulation based, meaning it produces synthetic catalogs that exhibit realistic behavior with respect to aftershocks and multi-fault earthquakes. A Poisson version of the model was constructed by randomizing event times, and the influence of two declustering algorithms was examined as well. We demonstrate that the probability of one-or-more loss exceedances (occurrence exceedance probability) is greater for the Poisson model because it has fewer seismically quiet time windows. The discrepancy between dollar loss estimates with a given exceedance probability is up to a factor of 32% but varies depending on the loss threshold (the x-axis value) and the forecast duration (we examined a range between 24 h and 50 years, with the discrepancy for the latter being negligible). We discuss how the one-or-more loss exceedance metric is questionable because it ignores all but the maximum loss experienced in each timeframe. An alternative metric based on total aggregate loss in each time window (aggregate exceedance probability) was therefore also examined, for which the Poisson model again implies higher risk at intermediate losses but lower risk at higher losses (because large, triggered events now contribute to total aggregate losses for the fully time-dependent model). We also argue that declustering is not a scientifically justifiable way to deal with full time dependence, in agreement with a chorus from other recent studies. It is difficult to draw generally applicable conclusions from our study, in part because application specific details will likely be important, but our results highlight how full time dependence can be reckoned with once authoritative forecast models are made available.

California↗
USGS1976

Corwin Springs Known Geothermal Resources Area, Park County, Montana

The Corwin Springs Known Geothermal Resources Area (KGRA) is contiguous to Yellowstone National Park along a part of the northern boundary of the park near Gardiner, Park County, Mont. The area contains two known sites of hot-spring activity--LaDuke Spring, 2.8 km southeast of the small resort community of Corwin Springs, and Bear Creek Spring, 2.6 km east of Gardiner. LaDuke Spring issues from brecciated quartzite and has a flow rate of 380 1/min, a surface water temperature of 65°C, a silica-geothermometer temperature of 68.7°C, and a Na-K-Ca geothermometer temperature of 76.8°C. Bear Creek Spring issues from limestone and has a flow rate of 4 1/min, a surface water temperature of 32°C, a silica temperature of 46.7°C, and a Na-K-Ca temperature of 87.2°C. The springs, which are actively depositing travertine, are on or near the trace of the Gardiner fault-a high-angle reverse fault which forms the southwestern boundary of the Beartooth uplift. Features that make the Corwin Springs KGRA potentially significant include: (1) Proximity of the Corwin Springs area to significant Pliocene and Pleistocene volcanism in Yellowstone National Park and related geothermal activity, (2) existence of a large potential geothermal reservoir in the northern extension of the Sepulcher Mountain graben, (3) localization of thermal activity along the Gardiner fault in the Corwin Springs area, (4) occurrence of negative gravity anomaly centered over the northern end of the Sepulcher Mountain graben, (5) existence of negative magnetic anomaly in the Corwin Springs area, (6) location of the area within the Intermountain seismic belt, (7) evidence of recent tectonism in the area, as demonstrated by Pleistocene and Holocene faulting, and (8) observed surface temperatures and estimated geochemical temperatures of the two known hot springs (features which do not indicate high subsurface temperatures).

Montana↗
USGS2024

Trends and environmental impacts of virtual water trade

Virtual water describes water embedded in the production of goods and offers meaningful insights about the complex interplay between water, trade and sustainability. In this Review, we examine the trends, major players, traded products and key drivers of virtual water trade (VWT). Roughly 20% of water used in global food production is traded virtually rather than domestically consumed. As such, agriculture dominates VWT, with livestock products, wheat, maize, soybean, oil palm, coffee and cocoa contributing over 70% of total VWT. These products are also driving VWT growth, the volume of which has increased 2.9 times from 1986 to 2022. However, the countries leading VWT contributions (with China, the United States, the Netherlands, Germany and India accounting for 34% of the global VWT in 2022) have remained relatively stable over time, albeit with China becoming an increasingly important importer. VWT can mitigate the effects of water scarcity and food insecurity, although there are concerns about the disconnect between consumers and the environmental impacts of their choices, and unsustainable resource exploitation. Indeed, approximately 16% of unsustainable water use and 11% of global groundwater depletion are virtually traded. Future VWT analyses must consider factors such as water renewability, water quality, climate change impacts and socioeconomic implications.

Nature Reviews Earth & Environment↗
USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗