USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers for Young Minds”

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,406 records · Page 26Linked to original sources
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗
USGS2026

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗
USGS2026

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗
USGS2024

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods↗
USGS2026

Generating geochemical and mineralogy distributions of soil in the conterminous United States using Bayesian hierarchical spatial models

Characterizing geochemical and mineralogical soil distributions across large spatial extents is essential for understanding mineral resources, ecosystem processes, and environmental risks. Rasters of soil geochemical distributions for the conterminous United States, however, are limited. We present a Bayesian modeling workflow and tool for generating predictive geochemical and mineralogy distribution maps for the conterminous United States using integrated nested Laplace approximation (INLA) with the stochastic partial differential equation approach. By modeling soil geostatistical data with environmental covariates (soil properties, topography, climate, and land cover), we generate predictive distributions of soil geochemistry that can be mapped or extracted for further analyses. As an example, we model the spatial distribution of trace elements in soil relevant to vertebrate health (cobalt, copper, iron, manganese, selenium, and zinc) and provide a workflow that can be used to generate and visualize predictive distributions of 39 other major and trace elements and 21 minerals of the soil survey, supporting a variety of ecological, environmental, and agricultural applications.

MethodsX↗
USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2025

Four-band image mosaic of the Colorado River Corridor downstream of Glen Canyon Dam in Arizona, derived from the May 2021 airborne image acquisition

In May 2021, the U.S. Geological Survey’s Grand Canyon Monitoring and Research Center acquired airborne multispectral high-resolution data for the Colorado River in the Grand Canyon, Arizona. The image data, which consist of four spectral bands (red, band 1; green, band 2; blue, band 3; and near infrared, band 4) with a ground resolution of 20 centimeters, are available as 16-bit unsigned-integer GeoTIFF files in Sankey and others (2024) (available online at https://doi.org/10.5066/P9BBGN6G ). The image files are projected in the State Plane Coordinate System, using the central Arizona zone (202) with the North American Datum of 1983 National Adjustment of 2011. The assessed spatial accuracy for these data is based on 47 ground-control points that were independent from the ground-control points used by the contractor for aerotriangulation and is reported at the 95-percent confidence level as 0.514 meter (m) and a root mean square error of 0.297 m. The intended uses of this dataset are primarily in support of scientific research and monitoring applications. Examples of these applications include high-resolution spatial and temporal change detection of the river channel, geomorphic landforms, riparian vegetation, and backwater and nearshore habitat, as well as other ecosystem-wide mapping. These imagery data also serve as reference material for field science mission planning, as base data for field data collection including community science activities, and as a highly detailed guide for technical boat operation during science activities such as reconnaissance for nighttime missions and navigating rapids during low flows.

Arizona, Nevada, Utah↗
USGS2025

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗
USGS2024

Peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. This work is part of that reexamination. The stationarity assumption is a concern because flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done to assess potential nonstationarity in peak flows in the north-central United States. This report summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Scientific Investigations Report↗
USGS2024

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed↗
USGS2026

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

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

Engineering Geology↗
USGS2025

Gaps in water quality modeling of hydrologic systems

This review assesses gaps in water quality modeling, emphasizing opportunities to improve next-generation models that are essential for managing water quality and are integral to meeting goals of scientific and management agencies. In particular, this paper identifies gaps in water quality modeling capabilities that, if addressed, could support assessments, projections, and evaluations of management alternatives to support ecosystem health and human beneficial use of water resources. It covers surface water and groundwater quality modeling, dealing with a broad suite of physical, biogeochemical, and anthropogenic drivers. Modeling capabilities for six constituents (or constituent categories) are explored: water temperature, salinity, nutrients, sediment, geogenic constituents, and contaminants of emerging concern. Each constituent was followed through the coupled atmospheric-hydrologic-human system, with prominent modeling gaps described for a diverse array of relevant inputs, processes, and human activities. Commonly identified modeling gaps primarily fall under three types: (1) model gaps, (2) data gaps, and (3) process understanding gaps. In addition to potential solutions for addressing specific individual modeling limitations, some broad approaches (e.g., enhanced data collection and compilation, machine learning, reduced-complexity modeling) are discussed as ways forward for tackling multiple gaps. This gap analysis establishes a framework of diverse approaches that may support improved process representation, scale, and accuracy of models for a wide range of water quality issues.

Water↗
USGS1990

Effects of storm-water runoff on local ground-water quality, Clarksville, Tennessee

Storm-related water-quality data were collected at a drainage-well site and at a spring site in Clarksville, Tennessee, to define the effects of storm-water runoff on the quality of ground water in the area. A dye-trace test verified the direct hydraulic connection between the drainage well and Mobley Spring. Samples of storm run off and spring flow were collected at these sites for nine storms during the period February to October 1988. Water samples were collected also from Mobley Spring and two other springs and two observation wells in the area during dry-weather conditions to assess the general quality of ground water in an urban karst terrain. Evaluation of the effect of storm-water runoff on the quality of local ground water is complicated by the presence of other sources of contaminants in the area Concentrations and load for most major constituents were much smaller in storm-water runoff at the drainage well than in the discharge of Mobley Spring, indicating that much of the chemical constituent load discharged from the spring comes from sources other than the drainage well. However, for some of the minor constituents associated with roadway runoff (arsenic, copper, lead, organic carbon, and oil and grease), the drainage well contributed relatively large amounts of these constituents to local ground water during storms. The close correlation between concentrations of total organic carbon and concentrations of most trace metals at the drainage-well and Mobley Spring sites indicates that these constituents are transported together. Many trace metals were flushed early during each runoff event. Mean storm loads for copper, lead, zinc, and four nutrient species (total nitrogen, ammonia nitrogen, total phosphorus, and orthophosphorus) in storm-water runoff at the drainage-well site were lower than mean storm load predicted from an existing regression model. The overprediction by the model may be a result of the small size of the drainage area relative to the range of drainage areas used in the development of the models, or to the below-normal amounts of rainfall during the period of sampling for this investigation. Loads& in storm-water runoff for 22 constituents were extrapolated from sampled storms to total loads for the period February to October 1988. Calculated loads for trace metals for the period ranged from 0.030pound.s for cadmium to 12pound.s for strontium. Loads of the primary nutrients ranged from 0.97pounds for nitrite as nitrogen to 34pounds of organic nitrogen. Storm-water quality at the drainage-well and Mobley Spring sites was compared to background water quality of the local aquifer; as characterized by dry-weather samples from three springs and two observation wells in the Clarksville area. Concentrations of total-recoverable cadmium, chromium, copper, lead, and nickel were higher in many stormwater samples from both the drainage-well and Mobley Spring sites than in samples from any other site. In addition, concentrations of total organic carbon, methylene blue active substances, and total-recoverable oil and grease were generally higher in storm-water samples from the drainage-well site than in any ground-water sample. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total recoverable iron, manganese, and methylene blue active substances in storm samples from the drainage-well site exceeded the maximum contaminant levels listed in Tennessee’s drinking-water standards (1988) by as much as 2,500 and 5,500 colonies per 100 milliliters, and 2.7, 0.29, and 0.05 milligrams per liter, respectively. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total-recoverable iron, manganese, and lead in storm samples from Mobley Spring exceeded the maximum contaminant levels by as much as 500 and 4,500 colonies per 100 milliliters, and 18.7,0.65, and 0.02 milligrams per liter, respectively. For iron, manganese, and bacteria, these undesirable levels are not necessarily attributable to storm-water recharge, because concentrations of these constituents also exceeded drinking-water standards in one or more of the dry-weather samples from selected springs and observation wells in the area.

Tennessee↗
USGS1995

Streamflow, ground-water recharge and discharge, and characteristics of surficial deposits in Buzzards Bay Basin, southeastern Massachusetts

Streamflow measurements at 14 low-flow partial-record stations and two discontinued streamflow-gaging stations were related to concurrent streamflows at six long-term streamflow-gaging stations to estimate streamflows at selected flow durations from 50 to 99 percent for water years 1967-91. At low flows, stream discharge per square mile generally increased with increasing percentage of subbasin underlain by stratified-drift deposits. At the 70-, 90-, and 99-percent flow durations, subbasins underlain primarily by stratified-drift deposits had stream discharges per square mile at least 2, 4, and 8 times greater, respectively, than subbasins underlain primarily by till and bedrock deposits. Streamflow measured three times at the Paskamanset River showed that streamflows downstream of municipal pumped wells increased at a lower rate per unit drainage area than streamflows upstream of the pumped wells. When pumpage from wells near the river was added to measured streamflow at the most downstream station, the sum was similar to the potential streamflow at the station if its streamflow per unit area was the same as that upstream from the pumped wells. Mean ground-water recharge and discharge rates were computed from continuous records of daily mean discharge during water years 1967-91 for six streamflow-gaging stations in southeastern Massachusetts and Rhode Island. Estimates of mean ground-water recharge were 19.7 to 22.6 inches per year for stations with drainage areas primarily underlain by till and bedrock deposits, and 23.8 to 25.2 inches per year for stations with drainage areas primarily underlain by stratified-drift deposits. During drought years, annual ground-water recharge can be less than 50 percent of the mean ground-water recharge rate for water years 1967-91. Estimates of ground-water discharge rates generally were 1 to 4 inches per year less than estimates of ground-water recharge rates. Ground-water discharge-duration curves were calculated from hydrograph separation of streamflow records at the six continuous streamflow-gaging stations for water years 1967-91. When streamflows were less than the median streamflow (50-percent duration), ground-water discharges averaged 85.7 percent of the total streamflow at the six streamflow-gaging stations. Ground-water discharges were estimated at the 14 low-flow partial-record stations and two discontinued streamflow-gaging stations by multiplying streamflows at the selected flow durations from 50 to 99 percent for water years 1967-91 by 85.7 percent. Data on surficial deposits in the basin were obtained by drilling 10 wells, by conducting 12 seismic-refraction surveys, and by collecting well log data for 313 wells from consulting firms and the Massachusetts Department of Environmental Management. These data indicated that a few minor corrections could be made to existing saturated thickness maps for stratified-drift deposits in Buzzards Bay Basin. Few well logs were detailed enough to estimate transmissivities; thus existing transmissivity maps were not updated for stratified-drift deposits in the basin.

Massachusetts↗
USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗
USGS2026

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS1996

Analysis of nutrients in the surface waters of the Georgia–Florida Coastal Plain study unit, 1970–91

During the early phase of the Georgia-Florida National Water Quality Assessment study, existing information on nutrients was compiled and analyzed in order to evaluate the nutrient concentrations within the 61,545 square mile study unit. Evaluation of the nutrient concentrations collected at surface-water sites between October 1, 1970, and September 30, 1991, utilized the environmental characteristics of land resource provinces, land use, and nonpoint and point-source discharges within the study unit. Long-term trends were investigated to determine the temporal distribution of nutrient concentrations. In order to determine a level of concern for nutrient concentrations, the U.S. Environmental Protection Agency (USEPA) guidelines were used: (1) for nitrate concentrations, the maximum contaminant level in public-drinking water supplies (10 mg/L); (2) for ammonia concentrations, the chronic exposure of aquatic organisms to un-ionized ammonia (2.1 mg/L); (3) for total-phosphorus concentrations, the recommended concentration in flowing water to discourage excessive growth of aquatic plants (0.1 mg/L); and (4) for kjeldahl concentrations, however, no guidelines were available. For sites within the 10 major river basins, median nutrient concentrations were enerally below USEPA guidelines, except for total-phosphorus concentrations where 45 percent of the medians exceeded the guideline. The only median ammonia concentration that exceeded the guideline occurred at the Swift Creek site (3.4 mg/L), in the Suwannee River basin, perhaps due to wastewater discharges. For all sites within the Withlacoochee, Aucilla, and St. Marys River basins, median concentrations of nitrate, ammonia, and total phosphorus were below the USEPA guidelines. Nutrient data at each monitoring site within each major basin were aggregated for comparisons of median nutrient concentrations among major basins. The Ochlockonee and Hillsborough River basins had the highest median nutrient concentrations, the Aucilla River basin had the lowest. Median concentrations of nitrate and ammonia among all major basins were below USEPA guidelines. The median total-phosphorus concentrations for the following river basins exceeded the USEPA guideline: Hillsborough, St. Johns, Suwannee, Ochlockonee, Satilla, Altamaha, and Ogeechee. Although nutrient concentrations within the study unit were low, long-term increasing trends were found in all four nutrients. All 18 study-unit wide nitrate trends had increasing slopes ranging from less than 0.01 to 0.07 (mg/L)/yr. The range in slope for the 13 ammonia trends was -0.03 to 0.01 (mg/L)/yr with 6 increasing trends in the northern part of the study unit. Of the 17 total-phosphorus trends found in the study unit, 10 were found at sites where the median concentration exceeded the USEPA guideline. At these 10 sites, 4 sites had increasing trends with slopes ranging from less than 0.01 to 0.07 (mg/L)/yr, 5 sites had decreasing trends with slopes ranging from -0.01 to -0.24 (mg/L)/yr, and one site showed a seasonal concentration trend. Median nutrient concentrations were significantly different among the four land resource provinces: Southern Piedmont, Southern Coastal Plain, Coastal Flatwoods, and Central Florida Ridge. As a result, nutrient concentrations among basins with similar nutrient inputs but located within different land resource provinces are not expected to be the same due to differences in the combination of factors such as soil permeability, runoff rates, and stream channel slopes. This concept is an important consideration in designing a surface-water quality network within the study area. For the most part, the Coastal Flatwoods showed the lowest median nutrient concentrations and the Southern Coastal Plain had the highest median nutrient concentrations. Lower median nitrate concentrations in surface-water basins were associated with the forest/wetland land-use category and higher median concentrations of nitrate and ammonia with the urban category when land-use percentages were classified into four land-use categories (agriculture, forest/wetland, mixed, and urban). These results were reasonable based on expected high nutrient inputs from urban areas and low inputs from forested and wetland areas. However, the lack of association between high nutrient concentrations and the agricultural land-use category was not expected since high nutrient inputs are generally needed for agriculture production.

Florida, Georgia↗