Search USGS⌕ Search

SEARCH · Search USGS

Results for “Water Research”

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.

At least 991 records · Page 55Linked to original sources
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS2026

Potential flood events in the lower Missouri River basin over multiple centuries identified using tree-ring based multi-model streamflow reconstructions

The Missouri River basin (MRB), the largest river basin in the United States, presents major water management challenges due to its complex topology and extensive infrastructure designed to manage high annual flows. Severe hydroclimatic events, particularly floods in 1993, 2011, and 2019, have highlighted vulnerabilities, leading the U.S. Army Corps of Engineers to enhance flood risk assessments in 2023. This study focuses on the full MRB, using tree-ring based streamflow reconstructions to extend the Basin's historic streamflow perspective. We developed annual streamflow reconstructions for the Hermann, Missouri streamgage from water years 1360–1900 using K-nearest neighbor and principal component regression methods. These reconstructions incorporate 183 tree-ring chronologies (1360–2020) and modeled natural streamflow (1901–‍2020). Reconstructions were evaluated for consistency using statistical measures and wavelet analysis. Both methods produced similar results, confirming their reliability in reconstructing multi-century streamflows. This novel methodology, guided by wavelet analysis, was developed to identify potential flood events in the lower Missouri River basin (LMRB), resulting in the identification of 172 potential flood events from water years 1360–2020, with 126 events common to the two versions of the MRB reconstructions. The reconstructed streamflow records showed an increase in the number of events until the late-18th- or mid-19th-century. These findings provide insights into evolving flood characteristics in the LMRB. The novel methodology developed in this study can be used to identify potential flood events from tree-ring based annual streamflow reconstructions and generally be applicable to other river basins.

Colorado, Iowa, Kansas, Missouri, Nebraska, Wyomin↗
USGS2026

Ecogeomorphologic response of a tidally restricted salt marsh to tidal restoration and sea-level rise

Diked marshes with decades of restricted tidal flow exhibit fundamentally altered hydrologic and ecologic dynamics that result in limited vertical accretion, land subsidence, and underdeveloped tidal channels. These changes complicate marsh restoration efforts and constrain the marsh's ability to build elevation in response to future sea levels. In this study, we developed an integrated hydrodynamic and marsh accretion model to conduct the first assessment of the restoration potential of a diked marsh in response to reintroduction of tidal flow and sea-level rise (SLR). We applied our model to evaluate changes in biomass production in the Herring River estuary, Cape Cod, MA, USA—which has experienced significant ecological shifts due to more than a century of tidal restriction—in response to a range of SLR scenarios, tidal flow conditions, management interventions, and biomass production rates through 2100. We found that with restoration of tidal flow, the simulated marsh extent expanded in areas with efficient drainage, increasing biomass production by up to 4.4 times. However, when SLR was imposed, parts of the marsh gradually transitioned to open water, leading to a decline in marsh coverage from 19% to as low as 12%. Management interventions aimed at enhancing marsh elevation and improving the drainage capacity of tidal channels slowed marsh loss, supporting marsh area through 2100. Our findings provide insight into the key drivers of marsh resilience and loss, and our approach can be applied across tidally restricted and restored systems to develop targeted, adaptive restoration strategies that achieve desired ecological outcomes.

Massachusetts↗
USGS2026

Global performance of remote sensing-based and reanalysis-driven models to estimate open water evaporation

Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.

Water Resources Research↗
USGS2026

Seasonal estimates of nutrient loading, sources, and impacts on water availability in streams across the conterminous United States, 2000–2020

Excess riverine nutrients can limit water availability for humans and ecosystems. Nutrient monitoring on individual stream reaches provides important information about current and past conditions, yet most river reaches are unmonitored. Additionally, monitored data integrate effects of upstream sources, sinks, and pathways by which nutrients get into streams, making explicit planning for effective nutrient reduction strategies challenging. We developed and applied dynamic SPAtially Referenced Regression On Watershed attributes (SPARROW) models to quantify total nitrogen (TN) and total phosphorus (TP) delivery and source-specific loading in all stream reaches of the conterminous United States (CONUS) at a seasonal timestep from 2000 through 2020. Between one third and one half of non-point source nutrients lagged more than one season in their delivery to CONUS streams over the study period. The largest sources of nutrients that were applied and delivered to streams within a season were agricultural fertilizers, treated wastewater, atmospheric deposition, and manure. Increases in loads were widespread across the CONUS and predominantly associated with increased loading from lagged delivery of non-point source inputs and current season fertilizer application. Model results also identified extensive areas where predicted seasonal concentrations frequently exceeded the National Rivers and Streams Assessment “poor” criteria. Increased loading occurred in areas with high yields and delivery to downstream receiving waterbodies, particularly across the Midwest, raising the risk of negative impacts to downstream water availability. The modeling approach presented here provides a continental-scale framework for identifying when, where, and why nutrients are generated and how that may ultimately impact receiving waterbodies.

conterminous United States↗
USGS2026

Refinement of a framework for Moving Aircraft River Velocimetry (MARV) and application to particle tracking along Alaskan rivers

Information on river velocities enhances understanding flood hazards, evaluating habitat conditions, and predicting the transport of floating materials. In this follow-up study, we used data from two new sites, one with a more complex morphology and the other with a lower suspended sediment concentration, to provide further evidence that Moving Aircraft River Velocimetry (MARV) can yield accurate velocity estimates ( R 2 up to 0.87 when compared to field measurements) for long segments of large, turbid rivers. The MARV workflow is packaged in freely available software and is robust to implementation details; neither buffering to mitigate edge effects nor a new approach to aggregating velocity vectors improved performance. MARV was not sensitive to parameters used to establish overlapping image sequences, but combining a long window with a short jump between consecutive windows was the optimal configuration. Although accuracy varied from one cross section to the next, agreement between remotely sensed velocities and those measured in the field was independent of position within a frame range. As an initial step toward application of the approach to help address practical problems, we showed how MARV can drive particle tracking models. Our first-order simulations suggest that channel morphology and flow velocity are the primary controls on travel time and particle fate, with diffusive processes playing a lesser role. Although MARV can be used to characterize an instantaneous flow field, a more comprehensive framework that accounts for other physical processes would be required to model specific types of events like oil spills.

Alaska↗
USGS2008

Comment on “Flow and tracer transport in a single fracture: A stochastic model and its relation to some field observations” by L. Moreno et al.

Moreno et al. [1988] (hereinafter referred to as MT) used a particle-tracking scheme to investigate the physics of solute movement in a variable-aperture planar fracture. The spatially heterogeneous fluid velocity was assumed to be the only mechanism of solute movement; local or pore scale dispersion and molecular diffusion were assumed to be negligible. The particle-tracking scheme used by MT consisted of routing particles from node to node in a finite difference grid. In this scheme, the direction of an individual particle is randomly selected and the probability associated with the particle movement in a given direction is proportional to the fluid flux in that direction. The same method was used by Desbarats [1990] to investigate advective transport in aquifers composed of two porous media of different hydraulic conductivities.

Water Resources Research↗
USGS2008

Sampling design for groundwater solute transport: Tests of methods and analysis of Cape Cod tracer test data

Tests of a one-dimensional sampling design methodology on measurements of bromide concentration collected during the natural gradient tracer test conducted by the U.S. Geological Survey on Cape Cod, Massachusetts, demonstrate its efficacy for field studies of solute transport in groundwater and the utility of one-dimensional analysis. The methodology was applied to design of sparse two-dimensional networks of fully screened wells typical of those often used in engineering practice. In one-dimensional analysis, designs consist of the downstream distances to rows of wells oriented perpendicular to the groundwater flow direction and the timing of sampling to be carried out on each row. The power of a sampling design is measured by its effectiveness in simultaneously meeting objectives of model discrimination, parameter estimation, and cost minimization. One-dimensional models of solute transport, differing in processes affecting the solute and assumptions about the structure of the flow field, were considered for description of tracer cloud migration. When fitting each model using nonlinear regression, additive and multiplicative error forms were allowed for the residuals which consist of both random and model errors. The one-dimensional single-layer model of a nonreactive solute with multiplicative error was judged to be the best of those tested. Results show the efficacy of the methodology in designing sparse but powerful sampling networks. Designs that sample five rows of wells at five or fewer times in any given row performed as well for model discrimination as the full set of samples taken up to eight times in a given row from as many as 89 rows. Also, designs for parameter estimation judged to be good by the methodology were as effective in reducing the variance of parameter estimates as arbitrary designs with many more samples. Results further showed that estimates of velocity and longitudinal dispersivity in one-dimensional models based on data from only five rows of fully screened wells each sampled five or fewer times were practically equivalent to values determined from moments analysis of the complete three-dimensional set of 29,285 samples taken during 16 sampling times.

Water Resources Research↗
USGS2008

Comment on the treatment of residual water content in “A consistent set of parametric models for the two-phase flow of immiscible fluids in the subsurface” by L. Luckner et al.

Luckner et al. [1989] (hereinafter LVN) present a clear summary and generalization of popular formulations used for convenient representation of porous media fluid flow characteristics, including water content ( θ ) related to suction ( h ) and hydraulic conductivity ( K ) related to θ or h. One essential but problematic element in the LVN models is the concept of residual water content ( θ r ; in LVN, θ w , r ). Most studies using θ r determine its value as a fitted parameter and make the assumption that liquid flow processes are negligible at θ values less than θ r . While the LVN paper contributes a valuable discussion of the nature of θ r , it leaves several problems unresolved, including fundamental difficulties in associating a definite physical condition with θ r , practical inadequacies of the models at low θ values, and difficulties in designating a main wetting curve.

Water Resources Research↗
USGS2008

Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 1. Experimental design and observed tracer movement

A large-scale natural gradient tracer experiment was conducted on Cape Cod, Massachusetts, to examine the transport and dispersion of solutes in a sand and gravel aquifer. The nonreactive tracer, bromide, and the reactive tracers, lithium and molybdate, were injected as a pulse in July 1985 and monitored in three dimensions as they moved as far as 280 m down-gradient through an array of multilevel samplers. The bromide cloud moved horizontally at a rate of 0.42 m per day. It also moved downward about 4 m because of density-induced sinking early in the test and accretion of areal recharge from precipitation. After 200 m of transport, the bromide cloud had spread more than 80 m in the direction of flow, but was only 14 m wide and 4–6 m thick. The lithium and molybdate clouds followed the same path as the bromide cloud, but their rates of movement were retarded about 50% relative to bromide movement because of sorption onto the sediments.

Massachusetts↗
USGS2008

Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 2. Analysis of spatial moments for a nonreactive tracer

A large-scale natural gradient tracer test was conducted to examine the transport of reactive and nonreactive tracers in a sand and gravel aquifer on Cape Cod, Massachusetts. As part of this test the transport of bromide, a nonreactive tracer, was monitored for about 280 m and quantified using spatial moments. The calculated mass of bromide for each sampling date varied between 85% and 105% of the injected mass using an estimated porosity of 0.39, and the center of mass moved at a nearly constant horizontal velocity of 0.42 m per day. A nonlinear change in the bromide longitudinal variance was observed during the first 26 m of travel distance, but afterward the variance followed a linear trend, indicating the longitudinal dispersivity had reached a constant value of 0.96 m. The transverse dispersivities were much smaller; transverse horizontal dispersivity was 1.8 cm, and transverse vertical dispersivity was about 1.5 mm.

Massachusetts↗
USGS2008

Centrifugal techniques for measuring saturated hydraulic conductivity

Centrifugal force is an alternative to large pressure gradients for the measurement of low values of saturated hydraulic conductivity ( K sat ). With a head of water above a porous medium in a centrifuge bucket, both constant-head and falling-head measurements are practical at forces up to at least 1800 times normal gravity. Darcy's law applied to the known centrifugal potential leads to simple formulas for K sat that are analogous to those used in the standard gravity-driven constant- and falling-head methods. Both centrifugal methods were tested on several fine-textured samples of soil and ceramic with K sat between about 10 −10 and 10 −9 m/s. The results were compared to falling-head gravity measurements. The comparison shows most measurements agreeing to within 20% for a given sample, much of the variation probably resulting from run-to-run changes in sample structure. The falling-head centrifuge method proved to be especially simple in design and operation and was more accurate than the constant-head method. With modified apparatus, K sat measurements less than 10 −10 m/s should be attainable.

Water Resources Research↗
USGS2008

Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 1. Geochemical assessment

A study was undertaken to quantitatively evaluate the hydrologic processes affecting the chemical and isotopic composition of drain lateral water in a drained agricultural field in the western San Joaquin Valley, California. The results elucidate the process of mixing of deep and shallow groundwater (below and within 6 m from land surface) entering the drain laterals. The deep groundwater was subject to evapoconcentration prior to drainage system installation and has been displaced downward (to depths greater than 6 m) in the groundwater system. The proportions of deep and shallow groundwater entering the drain laterals was calculated from the end-member oxygen 18 compositions determined in groundwater samples. The percentage of total drain lateral flow which is deep groundwater flow is about 30% for the shallow drain lateral (1.8 m below land surface) (drain lateral 1)) and 60% for the deep drain lateral (2.7 m below land surface (drain lateral 2)). During irrigation, the percentages of deep groundwater flow decrease to 0 and 30% for the shallow and deep drain laterals, respectively. Selenium concentrations in drain lateral waters decrease during irrigation but selenium loads increase. Total estimated annual loads were 1.1 and 5.4 kg of selenium for drain laterals 1 and 2, respectively. Substantial percentages of the annual load occurred during 8 days of irrigation, 23 and 9% for drain laterals 1 and 2, respectively.

California↗
USGS2008

Groundwater flow and solute movement to drain laterals, western San Joaquin Valley, California: 2. Quantitative hydrologic assessment

Groundwater flow modeling was used to quantitatively assess the hydrologic processes affecting ground water and solute movement to drain laterals. Modeling results were used to calculate the depth distribution of groundwater flowing into drain laterals at 1.8 m (drain lateral 1) and 2.7 m (drain lateral 2) below land surface. The simulations indicated that under nonirrigated conditions about 89% of the flow in drain lateral 2 was from groundwater originating from depths greater than 6 m below land surface. The deep groundwater has higher selenium concentrations than shallow groundwater. Simulation of irrigated conditions indicates that as recharge (deep percolation) increases, the proportional contribution of deep groundwater to drain lateral flow decreases. Groundwater flow paths and travel times estimated from the simulation results indicate that groundwater containing high concentrations of selenium (greater than 780 μg L −1 ) probably will continue to enter drain lateral 2 for decades.

California↗
USGS2008

Preliminary analysis of the hydrologic and geochemical controls on acid‐neutralizing capacity in two acidic seepage lakes in Florida

In late 1988, parallel studies of Lake Five‐O ( p H 5.14) in the Florida panhandle and Lake Barco ( p H 4.50) in north central Florida were initiated to develop hydrologic and major ion budgets of these lakes as part of an overall effort to improve understanding of the hydrologic, depositional, and biogeochemical factors that control acid‐neutralizing capacity (ANC) in seepage lakes. Preliminary findings from these studies indicate that earlier perceptions of lake hydrology and mechanisms of ANC regulation in Florida seepage lakes may have to be revised. The traditional perspective of seepage lakes in the Florida panhandle views these systems as dominated by precipitation inputs and that ANC regulation is due largely to in‐lake processes. Our results for Lake Five‐O show modest to steep hydraulic gradients almost entirely around the lake. In addition, the horizontal hydraulic conductivity of the surficial aquifer is high (8–74 m day −1 ), indicating that large quantities of groundwater flow into Lake Five‐O. Calculations of net groundwater flow from hydrologic budgets also indicate that groundwater may comprise at least 38 to 46% of the total inflow. For Lake Barco, net flow estimates of the minimum groundwater inflow range from 5 to 14% of total inflow. Enrichment factor and ion flux calculations for Lake Five‐O and Lake Barco indicate that terrestrial as well as in‐lake processes contribute significantly to ANC regulation. The extent that terrestrial processes contribute to ANC generation is directly related to the magnitude of groundwater inflow as well as the degree of ion enrichment or depletion that occurs in the surficial aquifer. Net ANC generation in both study lakes was dominated by anion retention (NO 3 − and SO 4 2− ). Where previous studies concluded that in‐lake reduction was the primary sink for SO 4 2− , our preliminary calculations show that adsorption of SO 4 2− within the watershed is perhaps twice as important as in‐lake reduction as a source of ANC. Net base cation enrichment in both lakes was negligible.

Florida↗
USGS2008

Evaluation of unsaturated zone air permeability through pneumatic tests

Predicting the steady state distribution of air pressure in the unsaturated zone resulting from a pneumatic test provides a method for determining air-phase permeability. This technique is analogous to the inverse problem of well hydraulics; however, air flow is more complicated than ground water flow because of air compressibility, the Klinkenberg effect, variations in air density and viscosity that result from temperature fluctuations in the unsaturated zone and the possibility of inducing water movement during the pneumatic test. An analysis of these complicating factors reveals that, when induced water movement can be neglected, a linear version of the airflow equation can provide an appropriate approximation for the purpose of determining air-phase permeability. Two analytical solutions for steady state, two-dimensional, axisymmetric airflow to a single well partially screened in the unsaturated zone are developed. One solution applies where there is a stratum of relatively low air permeability, separating the stratum in which the well is completed, from the atmosphere. The other solution applies where there is no separating stratum between the domain and atmosphere. In both situations the water table forms the lower horizontal boundary. Applications of both solutions to determine air permeability from data collected during pneumatic tests are presented.

Water Resources Research↗
USGS2008

Convergent radial dispersion: A note on evaluation of the Laplace transform solution

A numerical inversion algorithm for Laplace transforms that is capable of handling rapid changes in the computed function is applied to the Laplace transform solution to the problem of convergent radial dispersion in a homogeneous aquifer. Prior attempts by the author to invert this solution were unsuccessful for highly advective systems where the Peclet number was relatively large. The algorithm used in this note allows for rapid and accurate inversion of the solution for all Peclet numbers of practical interest, and beyond. Dimensionless breakthrough curves are illustrated for tracer input in the form of a step function, a Dirac impulse, or a rectangular input.

Water Resources Research↗
USGS2010

Gravity-driven groundwater flow and slope failure potential: 1. Elastic effective-stress model

Hilly or mountainous topography influences gravity-driven groundwater flow and the consequent distribution of effective stress in shallow subsurface environments. Effective stress, in turn, influences the potential for slope failure. To evaluate these influences, we formulate a two-dimensional, steady state, poroelastic model. The governing equations incorporate groundwater effects as body forces, and they demonstrate that spatially uniform pore pressure changes do not influence effective stresses. We implement the model using two finite element codes. As an illustrative case, we calculate the groundwater flow field, total body force field, and effective stress field in a straight, homogeneous hillslope. The total body force and effective stress fields show that groundwater flow can influence shear stresses as well as effective normal stresses. In most parts of the hillslope, groundwater flow significantly increases the Coulomb failure potential Φ, which we define as the ratio of maximum shear stress to mean effective normal stress. Groundwater flow also shifts the locus of greatest failure potential toward the slope toe. However, the effects of groundwater flow on failure potential are less pronounced than might be anticipated on the basis of a simpler, one-dimensional, limit equilibrium analysis. This is a consequence of continuity, compatibility, and boundary constraints on the two-dimensional flow and stress fields, and it points to important differences between our elastic continuum model and limit equilibrium models commonly used to assess slope stability.

Water Resources Research↗