USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources 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 757 records · Page 42Linked to original sources
USGS2026

Landscape controls river continuum dissolved organic carbon in northeastern Alaska

Rivers draining continuous permafrost landscapes on the North Slope of Alaska span from the mountains to the coast, integrating gradients in vegetation, soil organic matter, and hydrologic flow. We investigated how these landscape transitions influence dissolved organic carbon (DOC) concentrations, radiocarbon ages, and export along the river continuum over seasonal timeframes in northeastern Alaska. Water samples were collected from the Canning River and its tributaries from the spring freshet through late summer via multiple campaigns during 2011–2024. A simulation from the Permafrost Water Balance Model (PWBM) estimated terrain-specific contributions to freshwater and DOC export. River DOC concentrations peaked during the spring and declined in summer, with episodic increases linked to summer rainfall. Summer river DOC concentrations from the Brooks Range remained low (0.41–0.64 mg C L −1 ) but increased downstream (0.5–3.0 mg C L −1 ), becoming younger from the mountains (805 yBP) to the river delta (525 yBP) as flow from DOC-enriched streams draining the foothills (6–18 mg C L −1 ; Modern yBP) entered the Canning. The PWBM simulation revealed that while the mountains dominate the Canning's annual (90%) and spring-summer (88%–93%) freshwater export, the foothills disproportionately influence its annual (56%) and spring-summer (44%–58%) DOC export. DOC yields from the coastal plain were also high but had little influence due to its small contributing area. Our findings highlight the outsized role of the foothills in modulating river-integrated DOC export to coastal ecosystems. Considering the impact of these landscapes on carbon source and runoff will aid in anticipating future DOC fluxes from Arctic rivers.

Alaska↗
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 “Macrodispersion in sand-shale sequences” by A. J. Desbarats

Desbarats [1990] used a particle-tracking scheme to investigate the physics of three-dimensional solute transport in aquifers composed of two porous media of different hydraulic conductivities. 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 Desbarats 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 prob- ability associated with the particle movement in a given direction is proportional to the fluid flux in that direction. The same method was used by Moreno et al. [1988] to investigate advective transport in a variable-aperture planar fracture.

Water Resources Research↗
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

Response of Ned Wilson Lake watershed, Colorado, to changes in atmospheric deposition of sulfate

The Ned Wilson Lake watershed responds directly and rapidly to changes in precipitation inputs of sulfate, which has important implications for effects of acid deposition on the aquatic system. Chemistry at three precipitation collection sites and three watershed sites (a pond, a lake, and a spring) has been monitored in and near the Flattops Wilderness Area in northwestern Colorado beginning in 1981–1983. Bulk snowpack concentration of sulfate in the watershed and volume-weighted annual mean concentration of sulfate in precipitation at two nearby sites generally decreased from 1981 to 1985, were small through 1987, and increased in 1988–1989. Changes in concentration of sulfate at the watershed sites are controlled by precipitation inputs. Responsiveness of the individual sites was dependent on their position along the hydrologic flow path. The fastest response was in the pond, which has a hydrologic residence time of less than 1 year; over 90% of the variance in concentration of sulfate in the pond was explained by changes in concentration in precipitation. The lake has a hydrologic residence time of 1 to 4 years; a regression model of the concentration of sulfate in the lake, as a function of the concentration in the lake during the previous year and the concentration in precipitation, explained 87% of the variance in concentration of sulfate in the lake. The hydrologic response time of the spring is unknown; it was not responsive to changes in concentration of sulfate in precipitation. The recent increase of sulfate concentration in precipitation and in the pond and lake is evidence for a rapid rather than a delayed response, which could not be determined when only a decreasing trend in sulfate concentration was reported in 1982–1987. Watersheds of this type are sensitive to acidification (acid-neutralizing capacity less than 60 μeq L −1 ), and these results indicate conservative behavior of sulfate. This is important in predicting effects of future changes in atmospheric deposition, which could potentially be caused by anthropogenic emissions or climatic change.

Colorado↗
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↗
USGS2010

Arrival times and temporal moments of breakthrough curves for an imperfectly stratified aquifer

Using concepts developed in an earlier study, a solution in Laplace transform space is obtained for transport of resident concentration in an imperfectly but yet highly stratified porous medium. The flow field, into which an instantaneous pulse of tracer is injected, is taken to be steady and mean uniform parallel to the direction of stratification. From this transform-space solution either temporal moments can be derived by taking derivatives with respect to the Laplace parameter, or the transform-space solution can be inverted numerically to obtain breakthrough curves for the mean concentration. When compared to an equivalent solution with a Fickian dispersive flux, these temporal moments indicate the extent to which transport in heterogeneous porous media deviates from classical Fickian behavior. The numerical inversion of the Laplace transform solution gives partial breakthrough curves for the mean concentration which have the appearance of conflicting with the derived moment information. A hypothesis is put forth which resolves this apparent conflict; this hypothesis is verified by adding a component of local dispersion to the governing transport equation. On the basis of the flux-averaged concentration a form for the expected probability density function for the arrival time of a tracer particle is derived; arrival time moments and an arrival time cumulative distribution function are available as a consequence. Arrival time moments, as derived from the flux-averaged concentration, do not differ significantly from the resident moments, as derived from the resident concentration.

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↗