USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Environment”

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 76Linked to original sources
USGS2025

U.S. Geological Survey Groundwater Climate Response Network, 2024

As of October 2024, the U.S. Geological Survey (USGS) operated 588 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells selected to monitor the effects of climate variability, such as droughts, on groundwater levels nationwide. The CRN includes nearly 500 locations with real-time data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.

General Information Product↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Earthquake cycle mechanics during caldera collapse: Simulating the 2018 Kīlauea eruption

In multiple observed caldera-forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero-dimensional slider-block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non-vertical ring fault geometries often inferred at volcanoes. We present two-dimensional axisymmetric mafic piston-like collapse earthquake cycle models that include rate-and-state friction, fully-dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth-varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well-monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate-and-state friction parameterizations.

Hawaii↗
USGS1976

Guidelines for collection and field analysis of ground-water samples for selected unstable constituents

The unstable nature of many chemical and physical constituents in ground water requires special collection procedures and field analysis immediately after collection. This report describes the techniques and equipment commonly used in the collection and field analysis of samples for pH, temperature, carbonate, bicarbonate, specific conductance, Eh, and dissolved oxygen.

Open-File Report↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2026

King Rails in southeastern Arkansas are majority resident

King Rails ( Rallus elegans ), a declining species of secretive marsh bird, have a mixed migratory strategy, with some individuals migrating and some being year-round residents. Determining the distribution of migratory versus resident King Rails is critical to guiding habitat management and assessing risks to populations. We examined the migratory status of adult King Rails at the Choctaw West Wildlife Management Area in southeastern Arkansas, the only known site in Arkansas with multiple breeding pairs. Of the 23 adult King Rails we tagged in 2022 and 2023, two birds migrated away from the local area, one in each year. One migrated 145 km, and the other 450 km. Our findings suggest that the line between the migratory and resident range of Arkansas King Rails could be north of our site.

Arkansas↗
USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2025

Seamless Geologic Map Database for the Intermountain West, United States: A foundational dataset for mineral systems analysis

The Intermountain West has a complex geologic history, resulting in the formation of a diverse array of mineral deposits. Effective mineral exploration requires understanding the spatial and temporal relationships among geologic processes and events, a key focus of the mineral systems approach to exploration. This paper presents the Intermountain West Seamless Geologic Map Database, a unified dataset designed to support mineral exploration. Integrating geologic provinces, structural settings, and hydrothermal alteration, the database leverages the Seamless Integrated Geologic Mapping (SIGMa) extension to the USGS Geologic Map Schema (GeMS) to standardize geologic data from varied sources. SIGMa's hierarchical stratigraphic organization and feature-level metadata enhance data interoperability and reusability, enabling seamless query, analysis, and visualization of lithology, structural features, mineral deposits, geochronology, hydrothermal alteration. volcanic activity, and By providing a regionally consistent and dynamically evolving geologic map, this database provides a foundational framework for mineral exploration and geologic research. It also allows for an efficient workflow that expedites the publication of integrated geologic map databases.

Conference Paper↗
USGS2024

Automated, near real-time ground-motion processing at the U.S. Geological Survey

We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.

Seismological Research Letters↗
USGS2024

Effective site coefficients for the 2024 International Building Code (IBC)

The U.S. National Seismic Hazard Models (NSHMs), developed by the U.S. Geological Survey (USGS), have long been the scientific foundation of seismic design guidelines and have been used to compute design ground motions for construction of new buildings and retrofit of existing buildings in the United States and its territories. The 2018 NSHM is adopted by the 2024 International Building Code (IBC). Prior to the 2018 NSHM update, hazard calculations were required only at one reference site condition defined by V S30 =760 m/s (where V S30 is the time-averaged shear wave velocity from the surface to a depth of 30 m) and three periods (peak ground acceleration, PGA, and pseudo spectral accelerations at periods of 0.2 and 1 s, S s and S 1 ). Site coefficients, F PGA , F a , and F v , were then defined by the Building Seismic Safety Council (BSSC) Provisions Update Committee (PUC) in the site-specific procedures of National Earthquake Hazard Reduction Program (NEHRP) Recommended Seismic Provisions to calculate ground motions for other site classes with different V S30 values at the given periods. The design ground motions at other periods were then estimated using a generic spectral shape that was also defined by the BSSC PUC in NEHRP provisions. In recent years, the engineering community has realized there were deficiencies with the F a and F v site coefficients and the generic spectral shape. To avoid potentially dangerous underestimations of design ground motions for long period structures on soft site conditions in high seismicity regions, the BSSC PUC recommended the use of multi-period response spectra (MPRS) in 2017. As a result, the USGS produced multi-period and multi-V S30 response spectral values in the 2018 NSHM for calculations of design ground motions and the site coefficients F a and F v were eliminated from the 2020 NEHRP Provisions. As these site coefficients were widely used inside and outside of the United States, in this study we back-calculate the “effective” site coefficients F a,eff , and F v,eff by comparing MPRS for various site classes with the MPRS for the reference site condition, and discuss the changes that are observed in the 2024 IBC compared to its previous version in 2021. The effective site coefficients are presented for test site locations and their dependence on various factors including period, ground motion intensity, and regional models are discussed. Ratio maps between the new effective site coefficients and the old ones are then presented for soft site classes and for short and long periods. For soft site classes at short periods, the new effective site coefficients are lower than the old site coefficients for high seismicity regions and higher for low seismicity regions. As it was expected, for soft site classes at long periods and high seismicity regions, the new effective site coefficients are much larger than the old site coefficients without imposing the 50% increase as a penalty that was suggested in the 2021 IBC, whereas they could be much smaller if the 50% increase would have been imposed particularly around New Madrid and Charleston high seismicity regions. For low seismicity regions, the long period effective site coefficients can be smaller or larger by 20% compared to the 2021 IBC coefficients.

Conference Paper↗
USGS2026

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States↗
USGS2026

Egg thiamine concentrations reflect dietary differences between four hatchery strains, including two morphotypes, of Lake Trout in Lake Ontario

Objective Test whether dietary differences among hatchery strains of Lake Trout Salvelinus namaycush were associated with variation in egg thiamine concentrations and assess the potential for using strain-specific traits to mitigate thiamine deficiency complex (TDC) in Lake Ontario restoration efforts. Methods Stomach content analysis and stable isotope mixing models were used to compare diets among four hatchery strains of Lake Trout, representing lean and humper morphotypes. Egg thiamine concentrations were measured to test whether dietary differences correlated with susceptibility to TDC. Results Significant differences in diet and egg thiamine concentrations were observed among Lake Trout strains. A humper-like morphotype strain exhibited elevated thiamine levels, which were associated with reduced consumption of prey known to degrade thiamine levels when consumed, Alewife Alosa pseudoharengus and Rainbow Smelt Osmerus mordax , and increased reliance on Round Goby Neogobius melanostomus or sculpin species (Slimy Sculpin Cottus cognatus and Deepwater Sculpin Myoxocephalus thompsonii ). These findings support the link between the consumption of Alewife and Rainbow Smelt and thiamine deficiency. Conclusions Our findings support a link between strain-specific dietary preferences and susceptibility to TDC in Lake Trout. Preferential stocking of strains less prone to TDC, such as the strain derived from a humper morph from Lake Superior, could enhance restoration success in Lake Ontario. This strategy may be especially effective given the continued dominance of Alewife in the lake's prey base and the current management priorities supporting fisheries for Pacific salmon Oncorhynchus spp.

North American Journal of Fisheries Management↗
USGS2026

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗
USGS1991

Wetland sedimentation and vegetation patterns near selected highway crossings in West Tennessee

Wetland sedimentation and vegetation patterns at 11 highway crossings in West Tennessee were studied from 1987 to 1989. The purpose of the study was to investigate potential adverse effects of highway crossings on wetlands. Sedimentation rates, determined from root-burial depths, were highly variable. Average rates offine-grained deposition ranged from 0.005 to 0.033 foot per year for stations in locally ponded areas and from -0.002 to 0.039 foot per year for stations in drained areas. Sedimentation rates upstream from highway crossings were not significantly different from downstream rates at 8 of the 11 study sites. Three study sites had significantly greater sedimentation rates downstream. Sand splays were observed downstream from bridges at most study sites. Vegetation patterns and tree growth appear most strongly related to hydropetiod, defined as the average length of time an area is covered by water each year. The influence of sedimentation on tree growth is difficult to separate from the influence of hydroperiod because areas with high sedimentation rates typically have long hydroperiods. Estimated hydroperiod increased no more than 1 percent because of backwater from the highway crossings at the 11 study sites, while the estimated average depth of flood-plain inundation increased by an average of 6 percent.

Tennessee↗
USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗
USGS2025

Assessment of density pattern retention of generalized data for 1:100,000-scale United States topographic maps

Cartographic generalization reduces the complexity of geographic data to produce legible, smaller-scale displays that retain essential information and logical geographic patterns. Generalization is a vital process in topographic map production. An important challenge in this process is managing and evaluating consistency across scale in the density and spatial distribution of map features such as buildings, roads, streams, water bodies, and elevation contours. Density patterns in these features reflect underlying physiographic conditions, which include factors such as bedrock geology, tectonics, climate, and landforms. Assessments of an acceptable level of change in feature density patterns are critical to ensuring the readability, usability, and accuracy of generalized maps and data. Preserving realistic density patterns across mapping scales also supports sustainable development goals in cartography, by helping to prioritize and communicate the relative reliability of geospatial data at specific scales.

Conference Paper↗
USGS2024

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗