USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Water Process Engineering”

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,711 records · Page 20Linked to original sources
USGS2025

Disparate groundwater responses to wildfire

Post-wildfire investigations of groundwater response reveal a range of outcomes, varying from substantial increases to notable decreases in recharge and baseflow, with some studies indicating negligible or short-lived effects. This review assesses these varied responses within five critical categories: climate, vegetation, hydrogeology, fire characteristics, and the cryosphere, examining both short-term (within 2 years) and intermediate (2–10 years post-fire) effects. Despite considerable variability, some consistent patterns emerge. For instance, in hydroclimatic settings where water input and evaporative demand cycles are out of sync, post-wildfire groundwater responses tend to be positive (i.e., increased flux or storage), whereas under low fire severity conditions or in vegetation types that quickly recover, groundwater responses tend to be negative (i.e., decreased flux or storage). We synthesize relevant findings into a compendium of testable hypotheses aimed at explaining the spatiotemporal variability in observed post-wildfire groundwater responses. A recurring theme is the critical influence of the pre-wildfire groundwater regime on expected response and recovery. We identify opportunities for specific improvements in post-wildfire monitoring and modeling that would further advance capabilities to predict groundwater response. A key area for further research is understanding how wildfire effects on snow dynamics and other cryospheric processes translate to changes in groundwater.

WIREs Water↗
USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗
USGS2026

Open-source gravity reduction workflows for geothermal resource assessment

Potential-field geophysical data such as gravity can enhance understanding of geothermal resources at all stages of the resource life cycle, including assessment, exploration, development, and monitoring, and at multiple scales, from the reservoir scale to regional scale. However, to make gravity data useful for geothermal resource characterization, several processing steps are required to isolate the effects of density variations in the Earth’s crust to enable the identification of structural features associated with geothermal resources. Although this process is well-established, standard computational implementations for processing gravity data that are FAIR (Findable, Accessible, Interoperable, and Reproduceable) are still lacking. This paper details ongoing efforts at the U.S. Geological Survey (USGS) to develop a standard set of open-source Python tools for gravity data reduction that align with the FAIR principles. This workflow makes use of existing open-source tools for geophysical data processing with the goal of maximizing opportunities for rapid improvements, interoperability, and adaptability to other types of geophysical data.

Conference Paper↗
USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗
USGS2026

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 3 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km 2 ) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 3, the subject of this scientific investigations map, has a mapped area of 185 km 2 and has water depths that range from about 30 meters (m) on the Stellwagen Bank crest to about 135 m in a basin east of South Ninety Bank, which lies off the eastern margin of Stellwagen Bank. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates on the southeastern part of Stellwagen Bank, on adjacent banks and basins in deeper water to the east, in the eastern part of Race Point Channel to the south of the bank, and on the northern slope of Cape Cod. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 309 stations were analyzed, including 279 sediment samples. The geologic substrate maps of quadrangle 3 show the distribution of 21 geologic substrates that represent a wide range of textures, such as rippled sand, immobile sand, immobile muddy sand, sand that partially veneers gravel, and boulder ridges. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 3. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that use the various substrates, and support seabed management in the region.

Quadrangle 3 of the Stellwagen Bank National Marin↗
USGS2025

Relations of groundwater quality to long-term surface disposal of produced water near the Midway-Sunset and Buena Vista Oil Fields, California, USA

Contamination of groundwater by oil-field fluids in proximity to oil and gas development has been an issue of concern to water users and regulators given long histories of development and legacy disposal practices. A robust set of geochemical tracers including petroleum hydrocarbon compounds, thermogenic gases, inorganic ion concentrations, stable isotopes, radioactive isotopes, and noble gases were used to assess if oil-field fluids mixed with groundwater near the Midway-Sunset and Buena Vista Oil Fields in California, USA. Results show evidence of mixing of oil-field fluids with groundwater within the study area from either anthropogenic or natural processes. Produced water plumes associated with modern surface disposal facilities, used since the late 1950s, extend up to 1.5 km and currently remain within the boundaries of the oil fields. Plumes associated with earlier routing of produced water down natural drainages and in large retention structures (Midway Basin and Sunset Basin) near the Buena Vista Lake Bed are present in groundwater east of the oil fields. Based on geochemical tracer evidence, aerial imagery, and aerial electromagnetic surveys, these legacy plumes reach the western portion of the Central Valley aquifer system, an important groundwater resource for agricultural and domestic supply. The legacy plume associated with Sunset Basin may further be detected downgradient in deeper groundwater beneath the southern extent of the Buena Vista Lake Bed based on the presence of thermogenic gases and petroleum hydrocarbon compounds.

California↗
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↗
USGS2024

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 5 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km 2 ) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 5, the subject of this scientific investigations map, has an area of 211 km 2 and has water depths that range from 23 meters (m) on the Stellwagen Bank crest to 105 m in the Stellwagen Basin. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates on the crest and western flank of the south-central part of Stellwagen Bank and in Stellwagen Basin to the west. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 729 stations were analyzed, including 620 sediment samples. The geologic substrate maps of quadrangle 5 show the distribution of 20 substrates that represent a wide range of textures, such as mobile and rippled sand, immobile sand, sand that partially veneers gravel, boulder ridges, and mud. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 5. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that utilize the various substrates, and to support seabed management in the region.

Massachusetts↗
USGS1999

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota↗
USGS2024

Short-term sediment dispersal on a large retreating coastal river delta via 234Th and 7Be sediment geochronology: The Mississippi River Delta Front

Many Mississippi River Delta studies have shown recent declines in fluvial sediment load from the river and associated land loss. In contrast, recent sedimentary processes on the subaqueous delta are less documented. To help address this knowledge gap, multicores were collected offshore from the three main river outlets at water depths of 25–280 m in June 2017 just after the peak river discharge period, with locations selected based on 2017 U.S. Geological Survey seabed mapping. The coring locations included the undisturbed upper foreset, mudflow lobes, gullies, and the undisturbed prodelta. Nine multicores were analyzed for Beryllium-7 activity, and four cores were analyzed for excess Thorium-234 activity via gamma spectrometry, granulometry and X-radiography. Our results indicate a general trend of declining 7 Be and 234 Th activities and inventories with increasing distance from sources and in deeper water. The core X-radiographs are graded from the predominantly physically stratified nearshore to the more bioturbated offshore, consistent with the sedimentation patterns. Sediment focusing assessed via the 7 Be and 234 Th sediment inventories shows preferential sedimentation in gully and lobe environments, whereas the upper foreset and prodelta focusing factors are relatively depleted. Overall, short-term sediment deposition from the main fluvial source remains active offshore from all three major river outlets, despite the overall declining river load.

Louisiana↗
USGS2024

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California↗
USGS2025

Status of water-quality conditions in the United States, 2010–20

Degradation of water quality can make water harmful or unusable for humans and ecosystems. Although many studies have assessed the effect of individual constituents or narrow suites of constituents on freshwater systems, no consistent, comprehensive assessment exists over the wide range of water-quality effects on water availability. Using published studies, data, and models completed at regional or national scales in the United States during 2010–20, this chapter moves towards a comprehensive assessment by summarizing how selected anthropogenic and geogenic water-quality constituents affect national-scale water availability for human and ecosystem needs. Several types of human health, agricultural, ecological, and beneficial-use standards or thresholds were used to provide context for categorizing surface-water and groundwater quality. Water availability for human and ecological use is limited by elevated concentrations of geogenic and anthropogenic constituents in surface and groundwater. Elevated concentrations of five geogenic constituents (arsenic, manganese, strontium, radium, and adjusted gross alpha) are common in groundwater and collectively affect the drinking water supply to over 30 million people. Surface water sourced drinking water supplies are impaired in about a third of assessed stream miles, most commonly because of non-mercury metals and salinity. Health-based violations at community water systems may disproportionately affect socially vulnerable communities. Ecological water uses are predominantly limited by nutrients, sediment, temperature, pathogens, salinity, and pesticides. Water availability for human and ecological use is adversely affected by human activities including human contaminant sources (for example, wastewater, agriculture), processes (for example, dredging, groundwater pumping), or permanent landscape modifications (for example, dams, urbanization). Primary contaminant sources vary spatially and include fertilizer and manure, atmospheric deposition, wastewater treatment plants, urban land, and a range of natural sources. Contaminants of emerging concern, contaminants without regulatory thresholds, and mixtures of geogenic and anthropogenic water contaminants also contribute to ecological degradation and human exposure.

conterminous United States↗
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, ↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗
USGS2024

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water↗
USGS2024

Framework for implementing damping scaling factors in U.S. Geological Survey National Seismic Hazard Models

Traditionally, probabilistic seismic hazard analysis (PSHA) has focused on calculating ground motion hazard curves for elastic, 5%-damped pseudo spectral accelerations, Sa(T,5%), which are used as the basis for engineering design parameters and targets for ground motion selection and modification. However, structures and geotechnical systems can exhibit a wide range of damping ratios both above and below the 5% level, depending on the construction material, structural system, nonstructural elements, or subsurface soil properties. When spectral parameters at such damping levels are required for certain applications, 5%-damped accelerations have traditionally been extracted from PSHA-based hazard curves and adjusted outside of the hazard integral using damping scaling factors (DSF) such as those from Newmark & Hall (1982). Recent advances in the development of more rigorous and comprehensive damping scaling models (e.g., Rezaeian et al., 2014; Rezaeian et al., 2021) have allowed for the modeling of means and standard deviations of DSFs as functions of earthquake source and path properties for crustal, intraslab, and subduction interface tectonic environments. These DSF models can be applied to ground motion model (GMM) estimates of Sa(T,5%) for a given earthquake rupture scenario to produce a corresponding mean and standard deviation Sa at a specified damping ratio β, Sa(T,β). In this study, the DSF models of Rezaeian et al. (2014) and Rezaeian et al. (2021) are implemented within the U.S. Geological Survey National Seismic Hazard Model (NSHM) PSHA framework to calculate probabilistic hazard curves for spectral accelerations at damping ratios from 0.5% to 30%. The DSF models are applied directly to the mean and standard deviation of Sa(T,5%) predictions from each GMM in the NSHM logic tree. Resulting hazard curves and uniform hazard and risk spectra for Sa(T,β) are presented for several geographic locations and compared with corresponding spectra estimated using current design practices by applying the same DSFs outside of the PSHA calculation. Key differences between the two methods for estimating Sa(T,β) are discussed, and potential strategies are presented for the implementation and usage of the hazard-consistent Sa(T,β) in building codes. Comparing the results to those from DSFs used in current design practices that are mainly based on Newmark & Hall (1982) is not explored in this study.

Conference Paper↗
USGS2024

Joint pilot fish habitat framework

This story map will take you through the process of exploring and testing methods necessary for a higher resolution, seamless fish habitat assessment across both inland and estuarine waters through the lens of our joint pilot assessment Fish habitat assessments attempt to relate past, current, or future landscape conditions to the state of fish species occurrence, distribution, abundance, or community and habitat condition in streams, rivers, or estuaries. Previous fish habitat assessments, such as the National Fish Habitat Assessment, conducted separate and disconnected assessments for inland waters and estuaries. In this project, National Oceanic and Atmospheric Administration ( NOAA ) and U.S. Geological Survey ( USGS ) researchers created a seamless spatial framework to allow assessments that integrate influences on fish habitat from headwaters to the estuary. This effort began when the Chesapeake Bay Program Fish Habitat Action Team expressed interest in a Baywide fish habitat assessment spanning tidal salt, tidal fresh, warm non-tidal, and cold non-tidal waters. However, the complexity of the myriad of implementation details to consider when developing such an assessment necessitated the need for a tributary-specific pilot assessment. To conduct this pilot assessment, a NOAA/USGS joint partnership was formed with cooperation and support from the Chesapeake Bay Agreement and Chesapeake Bay Fish Habitat Action Team (FHAT).

Maryland↗
USGS2024

Water-quality characteristics of the Red River of the North and tributaries in the Fargo-Moorhead metropolitan area, North Dakota, 2019–22

The Flood Risk Management Project was initiated in 2008 in the Fargo-Moorhead metropolitan area to reduce flood risk, flood damages, and flood protection costs in the Fargo-Moorhead metropolitan area. In cooperation with the U.S. Army Corps of Engineers, the U.S. Geological Survey initiated a water-quality monitoring study to describe the water-quality characteristics of the Red River of the North and its tributaries in the Fargo-Moorhead metropolitan area during the preconstruction period of the Flood Risk Management Project from October 1, 2019, to October 1, 2022. The monitoring study included the collection of discrete and continuous water-quality data and streamflow monitoring at selected sites that integrated and enhanced existing monitoring programs within the study area. Discrete samples collected at 10 sites in the Fargo-Moorhead metropolitan area were analyzed for major ions, trace elements, nutrients, suspended sediment, pesticides, and fecal indicator bacteria. In general, major ion concentrations were higher at sites on the tributaries (Wild Rice, Sheyenne, and Maple Rivers) compared to sites on the Red River of the North. In general, bicarbonate, calcium, magnesium, and sulfate represented most of the dissolved ions measured in samples collected at the 10 sites. Calcium, chloride, fluoride, potassium, silica, and sodium were also measured in samples, but they represented a smaller portion of the total dissolved ions. Sulfate was the most dominant dissolved ion that had the highest concentrations among the major ions measured in samples. A total of 18 trace elements were analyzed in discrete samples. Several of the trace elements had concentrations below the laboratory reporting level in all of the samples, including antimony, beryllium, cadmium, chromium, silver, and thallium. Sites on the Wild Rice River generally had the highest concentrations of arsenic, barium, boron, manganese, and nickel compared to the other sites. Nutrients analyzed in discrete samples included filtered and unfiltered concentrations of ammonia, nitrate plus nitrite, phosphorus, and organic carbon. The median filtered ammonia concentration at most sites was less than the laboratory reporting level of 0.03 milligram per liter as nitrogen except for the Sheyenne River at Harwood, North Dakota (U.S. Geological Survey [USGS] station 05060400), and Red River of the North near Georgetown, Minnesota (USGS station 05062130). The lowest median unfiltered nitrate plus nitrite concentration was measured at sites on the Red River of the North upstream from the Fargo-Moorhead metropolitan area and the highest median was at sites on the Red River of the North downstream from the Fargo-Moorhead metropolitan area compared to all other sites. The increase in nitrate plus nitrite concentrations could reflect the effect of the wastewater-treatment plant discharge that enters the Red River of the North upstream from the site located downstream from the Fargo-Moorhead metropolitan area and from urban runoff. Phosphorus (unfiltered) concentrations were generally higher at sites on the Maple and Sheyenne Rivers compared to the other sites and were higher at sites on the Red River of the North downstream from the Fargo-Moorhead metropolitan area compared to sites upstream on the Red River of the North. Suspended-sediment concentrations were generally highest at sites in the Sheyenne River and lowest in the upstream Red River of the North sites. Suspended-sediment concentration was highly variable in samples collected at the 10 sites, mostly influenced by the occurrence of snowmelt and rainfall-runoff events. The Sheyenne River near Kindred, N. Dak. (USGS station 05059000) had the largest range in sediment concentrations in samples collected at the 10 sites. For all sites other than the Sheyenne River near Kindred, N. Dak., 95 percent or more of the suspended sediment had particle diameter sizes less than 0.0625 millimeter in 50 percent of the samples (median). Of the 102 pesticides and pesticide degradates analyzed, 45 constituents had no detectable concentrations in any of the 17 samples collected at five sites. The remaining 57 pesticides had at least one detection in the samples collected at the five sites. The sites on the Wild Rice River (near Abercrombie, N. Dak., USGS station 05053000, and near St. Benedict, N. Dak., USGS station 05053500) and Sheyenne River near Kindred, N. Dak., had fewer pesticide detections compared to the Maple River below Mapleton, N.Dak. (USGS station 05060100) and the Red River of the North at Fargo, N. Dak (USGS station 05054000) and near Georgetown, Minn. Patterns in annual loads generally followed the same pattern as streamflow at the 10 sites for water years 2020–22. A water year is the 12-month period from October 1 to September 30 and is designated by the calendar year in which it ends. The greatest loads for all constituents were delivered at the two downstream sites on the Red River of the North; sites that also had the highest annual streamflows among the sites and the greatest loads were delivered in water year 2020 when the highest streamflows occurred at the sites. Likewise, the least loads for most constituents were at the Maple River and were least in 2021 compared to the other years because of low-streamflow conditions. Water-quality measurements continuously recorded at the Red River of the North at Hickson, N. Dak. (USGS station 05051522); Red River of the North at Fargo, N. Dak.; and Red River of the North near Georgetown, Minn. included water temperature, specific conductance, dissolved oxygen, pH, and turbidity. Specific conductance values were similar for the Red River of the North near Hickson, N. Dak., and Red River of the North at Fargo, N. Dak., when compared to the Red River of the North near Georgetown, Minn. that had higher values than the other two sites. Dissolved oxygen concentrations and pH were similar among the three sites on the Red River. The patterns in turbidity were mostly related to streamflow conditions and were similar among the three sites on the Red River of the North.

Minnesota, North Dakota↗