Geology Reports⌕ Search

SEARCH · Geology Reports

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 127 records · Page 7Linked to original sources
USGS2010

Confirmation of rate-dependent behavior in water retention during drainage in nonswelling porous materials

In the water retention process in porous material, changes in water content are assumed to be independent of the rate at which the matric potential varies. Periodically, researchers have presented results that do not support this assumption, indicating that water retention may be rate-dependent under certain conditions. In the present study, long-term experiments were performed on five porous materials in which core samples were drained at different rates down to matric potentials of −50 kPa. Comparisons of these results with previous results confirm that slight rate-dependent behavior occurs during drainage of porous materials in this matric potential range. Specifically, new and previous results indicate that on average 4–5% more water was retained for a slow, multistep drainage versus a fast, one-step drainage treatment. For the present study, average results for individual materials were variable, ranging from a 1 to a 10% difference in water content due to treatment. Several possible mechanisms for the observed rate-dependent behavior are discussed. The only plausible mechanism is related to variations in pore water salt concentrations that induce differences in the pore water surface tension for fast versus slow drainage. Regardless of the mechanism, results suggest that rate-dependent behavior during water retention may contribute to an enhanced temperature dependence of water retention, due to higher rates of water redistribution at higher temperatures.

Water Resources Research↗
USGS2017

Water quality data for national-scale aquatic research: The Water Quality Portal

Aquatic systems are critical to food, security, and society. But, water data are collected by hundreds of research groups and organizations, many of which use nonstandard or inconsistent data descriptions and dissemination, and disparities across different types of water observation systems represent a major challenge for freshwater research. To address this issue, the Water Quality Portal (WQP) was developed by the U.S. Environmental Protection Agency, the U.S. Geological Survey, and the National Water Quality Monitoring Council to be a single point of access for water quality data dating back more than a century. The WQP is the largest standardized water quality data set available at the time of this writing, with more than 290 million records from more than 2.7 million sites in groundwater, inland, and coastal waters. The number of data contributors, data consumers, and third-party application developers making use of the WQP is growing rapidly. Here we introduce the WQP, including an overview of data, the standardized data model, and data access and services; and we describe challenges and opportunities associated with using WQP data. We also demonstrate through an example the value of the WQP data by characterizing seasonal variation in lake water clarity for regions of the continental U.S. The code used to access, download, analyze, and display these WQP data as shown in the figures is included as supporting information.

Water Resources Research↗
USGS2011

3-D flow and scour near a submerged wing dike: ADCP measurements on the Missouri River

Detailed mapping of bathymetry and three-dimensional water velocities using a boat-mounted single-beam sonar and acoustic Doppler current profiler (ADCP) was carried out in the vicinity of two submerged wing dikes located in the Lower Missouri River near Columbia, Missouri. During high spring flows the wing dikes become submerged, creating a unique combination of vertical flow separation and overtopping (plunging) flow conditions, causing large-scale three-dimensional turbulent flow structures to form. On three different days and for a range of discharges, sampling transects at 5 and 20 m spacing were completed, covering the area adjacent to and upstream and downstream from two different wing dikes. The objectives of this research are to evaluate whether an ADCP can identify and measure large-scale flow features such as recirculating flow and vortex shedding that develop in the vicinity of a submerged wing dike; and whether or not moving-boat (single-transect) data are sufficient for resolving complex three-dimensional flow fields. Results indicate that spatial averaging from multiple nearby single transects may be more representative of an inherently complex (temporally and spatially variable) three-dimensional flow field than repeated single transects. Results also indicate a correspondence between the location of calculated vortex cores (resolved from the interpolated three-dimensional flow field) and the nearby scour holes, providing new insight into the connections between vertically oriented coherent structures and local scour, with the unique perspective of flow and morphology in a large river.

Missouri↗
USGS2001

Confidence intervals for expected moments algorithm flood quantile estimates

Historical and paleoflood information can substantially improve flood frequency estimates if appropriate statistical procedures are properly applied. However, the Federal guidelines for flood frequency analysis, set forth in Bulletin 17B, rely on an inefficient “weighting” procedure that fails to take advantage of historical and paleoflood information. This has led researchers to propose several more efficient alternatives including the Expected Moments Algorithm (EMA), which is attractive because it retains Bulletin 17B's statistical structure (method of moments with the Log Pearson Type 3 distribution) and thus can be easily integrated into flood analyses employing the rest of the Bulletin 17B approach. The practical utility of EMA, however, has been limited because no closed‐form method has been available for quantifying the uncertainty of EMA‐based flood quantile estimates. This paper addresses that concern by providing analytical expressions for the asymptotic variance of EMA flood‐quantile estimators and confidence intervals for flood quantile estimates. Monte Carlo simulations demonstrate the properties of such confidence intervals for sites where a 25‐ to 100‐year streamgage record is augmented by 50 to 150 years of historical information. The experiments show that the confidence intervals, though not exact, should be acceptable for most purposes.

Water Resources Research↗
USGS2000

A compartmentalized solute transport model for redox zones in contaminated aquifers: 2. Field‐scale simulations

This paper, the second of two parts [see Abrams and Loague , this issue], reports the field‐scale application of COMPTRAN (compartmentalized solute transport model) for simulating the development of redox zones. COMPTRAN is fully developed and described in the companion paper. Redox zones, which are often delineated by the relative concentrations of dissolved oxygen, have been observed around the globe. The distribution of other redox‐sensitive species is affected by redox zonation. At the U.S. Geological Survey's Cape Cod research site, an anoxic zone containing high concentrations of dissolved iron has been observed. Field data were abstracted from the Cape Cod site for the one‐dimensional and two‐dimensional COMPTRAN simulations reported in this paper. The purpose of the concept‐development simulations was to demonstrate that the compartmentalized approach reported by Abrams et al . [1998] can be linked with a solute transport model to simulate field‐scale phenomena. The results presented in this paper show that COMPTRAN successfully simulated the development of redox zones at the field scale, including trends in p H and alkalinity. Thermodynamic constraints were used to prevent lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium among all redox species. Empirical methods of reaction inhibition were not needed for the simulations conducted for this study. COMPTRAN can be extended easily to include additional compartments and reactions and is capable of handling complex velocity fields in more than one dimension.

Water Resources Research↗
USGS2013

Use of NMR logging to obtain estimates of hydraulic conductivity in the High Plains aquifer, Nebraska, USA

Hydraulic conductivity (K) is one of the most important parameters of interest in groundwater applications because it quantifies the ease with which water can flow through an aquifer material. Hydraulic conductivity is typically measured by conducting aquifer tests or wellbore flow (WBF) logging. Of interest in our research is the use of proton nuclear magnetic resonance (NMR) logging to obtain information about water-filled porosity and pore space geometry, the combination of which can be used to estimate K. In this study, we acquired a suite of advanced geophysical logs, aquifer tests, WBF logs, and sidewall cores at the field site in Lexington, Nebraska, which is underlain by the High Plains aquifer. We first used two empirical equations developed for petroleum applications to predict K from NMR logging data: the Schlumberger Doll Research equation (K SDR ) and the Timur-Coates equation (K T-C ), with the standard empirical constants determined for consolidated materials. We upscaled our NMR-derived K estimates to the scale of the WBF-logging K(K WBF-logging ) estimates for comparison. All the upscaled K T-C estimates were within an order of magnitude of K WBF-logging and all of the upscaled K SDR estimates were within 2 orders of magnitude of K WBF-logging . We optimized the fit between the upscaled NMR-derived K and KWBF-logging estimates to determine a set of site-specific empirical constants for the unconsolidated materials at our field site. We conclude that reliable estimates of K can be obtained from NMR logging data, thus providing an alternate method for obtaining estimates of K at high levels of vertical resolution.

Nebraska↗
USGS2010

Estimation of aquifer scale proportion using equal area grids: assessment of regional scale groundwater quality

The proportion of an aquifer with constituent concentrations above a specified threshold (high concentrations) is taken as a nondimensional measure of regional scale water quality. If computed on the basis of area, it can be referred to as the aquifer scale proportion. A spatially unbiased estimate of aquifer scale proportion and a confidence interval for that estimate are obtained through the use of equal area grids and the binomial distribution. Traditionally, the confidence interval for a binomial proportion is computed using either the standard interval or the exact interval. Research from the statistics literature has shown that the standard interval should not be used and that the exact interval is overly conservative. On the basis of coverage probability and interval width, the Jeffreys interval is preferred. If more than one sample per cell is available, cell declustering is used to estimate the aquifer scale proportion, and Kish's design effect may be useful for estimating an effective number of samples. The binomial distribution is also used to quantify the adequacy of a grid with a given number of cells for identifying a small target, defined as a constituent that is present at high concentrations in a small proportion of the aquifer. Case studies illustrate a consistency between approaches that use one well per grid cell and many wells per cell. The methods presented in this paper provide a quantitative basis for designing a sampling program and for utilizing existing data.

Water Resources Research↗
USGS2017

Steady state fractionation of heavy noble gas isotopes in a deep unsaturated zone

To explore steady state fractionation processes in the unsaturated zone (UZ), we measured argon, krypton, and xenon isotope ratios throughout a ∼110 m deep UZ at the United States Geological Survey (USGS) Amargosa Desert Research Site (ADRS) in Nevada, USA. Prior work has suggested that gravitational settling should create a nearly linear increase in heavy-to-light isotope ratios toward the bottom of stagnant air columns in porous media. Our high-precision measurements revealed a binary mixture between (1) expected steady state isotopic compositions and (2) unfractionated atmospheric air. We hypothesize that the presence of an unsealed pipe connecting the surface to the water table allowed for direct inflow of surface air in response to extensive UZ gas sampling prior to our first (2015) measurements. Observed isotopic resettling in deep UZ samples collected a year later, after sealing the pipe, supports this interpretation. Data and modeling each suggest that the strong influence of gravitational settling and weaker influences of thermal diffusion and fluxes of CO 2 and water vapor accurately describe steady state isotopic fractionation of argon, krypton, and xenon within the UZ. The data confirm that heavy noble gas isotopes are sensitive indicators of UZ depth. Based on this finding, we outline a potential inverse approach to quantify past water table depths from noble gas isotope measurements in paleogroundwater, after accounting for fractionation during dissolution of UZ air and bubbles.

Nevada↗
USGS2002

The role of bedrock topography on subsurface storm flow

We conducted a detailed study of subsurface flow and water table response coupled with digital terrain analysis (DTA) of surface and subsurface features at the hillslope scale in Panola Mountain Research Watershed (PMRW), Georgia. Subsurface storm flow contributions of macropore and matrix flow in different sections along an artificial trench face were highly variable in terms of timing, peak flow, recession characteristics, and total flow volume. The trench flow characteristics showed linkages with the spatial tensiometer response defining water table development upslope. DTA of the ground surface did not capture the observed spatial patterns of trench flow or tensiometric response. However, bedrock surface topographic indices significantly improved the estimation of spatial variation of flow at the trench. Point‐scale tensiometric data were also more highly correlated with the bedrock surface‐based indices. These relationships were further assessed for temporal changes throughout a rainstorm. Linkages between the bedrock indices and the trench flow and spatial water table responses improved during the wetter periods of the rainstorm, when the hillslope became more hydrologically connected. Our results clearly demonstrate that in developing a conceptual framework for understanding the mechanisms of runoff generation, local bedrock topography may be highly significant at the hillslope scale in some catchments where the bedrock surface acts as a relatively impermeable boundary.

Water Resources Research↗
USGS2003

Time‐lapse imaging of saline‐tracer transport in fractured rock using difference‐attenuation radar tomography

Accurate characterization of fractured‐rock aquifer heterogeneity remains one of the most challenging and important problems in groundwater hydrology. We demonstrate a promising strategy to identify preferential flow paths in fractured rock using a combination of geophysical monitoring and conventional hydrogeologic tests. Cross‐well difference‐attenuation ground‐penetrating radar was used to monitor saline‐tracer migration in an experiment at the U.S. Geological Survey Fractured Rock Hydrology Research Site in Grafton County, New Hampshire. Radar data sets were collected every 10 min in three adjoining planes for 5 hours during each of 12 tracer tests. An innovative inversion method accounts for data acquisition times and temporal changes in attenuation during data collection. The inverse algorithm minimizes a combination of two functions. The first is the sum of weighted squared data residuals. Second is a measure of solution complexity based on an a priori space‐time covariance function, subject to constraints that limit radar‐attenuation changes to regions of the tomograms traversed by high difference‐attenuation ray paths. The time series of tomograms indicate relative tracer concentrations and tracer arrival times in the image planes; from these we infer the presence and location of a preferential flow path within a previously identified zone of transmissive fractures. These results provide new insights into solute channeling and the nature of aquifer heterogeneity at the site.

New Hampshire↗
USGS2015

A Bayesian kriging approach for blending satellite and ground precipitation observations

Drought and flood management practices require accurate estimates of precipitation. Gauge observations, however, are often sparse in regions with complicated terrain, clustered in valleys, and of poor quality. Consequently, the spatial extent of wet events is poorly represented. Satellite-derived precipitation data are an attractive alternative, though they tend to underestimate the magnitude of wet events due to their dependency on retrieval algorithms and the indirect relationship between satellite infrared observations and precipitation intensities. Here we offer a Bayesian kriging approach for blending precipitation gauge data and the Climate Hazards Group Infrared Precipitation satellite-derived precipitation estimates for Central America, Colombia, and Venezuela. First, the gauge observations are modeled as a linear function of satellite-derived estimates and any number of other variables—for this research we include elevation. Prior distributions are defined for all model parameters and the posterior distributions are obtained simultaneously via Markov chain Monte Carlo sampling. The posterior distributions of these parameters are required for spatial estimation, and thus are obtained prior to implementing the spatial kriging model. This functional framework is applied to model parameters obtained by sampling from the posterior distributions, and the residuals of the linear model are subject to a spatial kriging model. Consequently, the posterior distributions and uncertainties of the blended precipitation estimates are obtained. We demonstrate this method by applying it to pentadal and monthly total precipitation fields during 2009. The model's performance and its inherent ability to capture wet events are investigated. We show that this blending method significantly improves upon the satellite-derived estimates and is also competitive in its ability to represent wet events. This procedure also provides a means to estimate a full conditional distribution of the “true” observed precipitation value at each grid cell.

Water Resources Research↗
USGS2021

Seasonality of solute flux and water source chemistry in a coastal glacierized watershed undergoing rapid change: Wolverine Glacier watershed, Alaska

As glaciers around the world rapidly lose mass, the tight coupling between glaciers and downstream ecosystems is resulting in widespread impacts on global hydrologic and biogeochemical cycling. However, a range of challenges make it difficult to conduct research in glacierized systems and our knowledge of seasonally changing hydrologic processes and solute sources and signatures is limited. This in turn hampers our ability to make predictions on solute composition and flux. We conducted a broad water sampling campaign in order to understand the present-day partitioning of water sources and associated solutes in Alaska’s Wolverine Glacier watershed. We established a relationship between electrical conductivity (EC) and streamflow at the watershed outlet to divide the melt season into four hydroclimatic periods. Across hydroclimatic periods, we observed a shift in off-glacier source waters from snowmelt-dominated overland and shallow subsurface flow paths to deeper groundwater flow paths. We also observed the shift from a low- to high-efficiency subglacial drainage network and the associated flushing of water stored sub-glacially with higher solute loads. We used calcium, the dominant dissolved ion, from watershed outlet samples to estimate solute fluxes for each hydroclimatic period across two melt seasons. We found between 40 and 55 percent of Ca2+ export occurred during the late season rainy period. This partitioning of the melt season coupled with a characterization of the chemical makeup and magnitude of solute export provides new insight into a rapidly changing watershed and creates a framework to quantify and predict changes to solute fluxes across a melt season.

Alaska↗
USGS2018

Quantifying climate-related interactions in shallow and deep storage and evapotranspiration in a forested, seasonally water-limited watershed in the Southeastern United States

The Southeastern United States experiences recurring hydrological droughts, which can reduce water availability and can result in water-limiting conditions. Long-term monitoring at Panola Mountain Research Watershed, a small, forested, seasonally water-limited watershed near Atlanta, Georgia, was used to quantify the interactions of climatic variability with shallow and deep storage and evapotranspiration. Watershed storage (WS) and actual evapotranspiration (AET) were estimated monthly from 1985 through 2015 using a water-budget approach combined with a WS-baseflow relationship. Shallow storage (SS) was assessed from a soil moisture profile. Soil moisture transitioned from recharge to surplus as SS increased from its field capacity to a nearly saturated state during the dormant season, and transitioned from utilization to climatic water deficits as SS declined from its field capacity to its wilting point during the growing season. Deeper storage was unavailable to AET during dry conditions. The majority of deeper storage recharge occurred during the dormant season and required SS to be wet. WS was an effective drought indicator. Growing season droughts typically occurred when WS was below normal at the end of the dormant season and growing season precipitation (P) was below or near normal. A hydrologic persistence analysis found that monthly-standardized WS was significantly correlated (p-value <0.05) with past monthly-standardized WS for the previous 19 months and with past monthly P for the previous 11 months, indicating the importance of past hydrologic conditions on WS. Expected climatic changes affected recharge during the dormant season and deficits during the growing season.

Water Resources Research↗
USGS2012

Exploring changes in the spatial distribution of stream baseflow generation during a seasonal recession

Relating watershed structure to streamflow generation is a primary focus of hydrology. However, comparisons of longitudinal variability in stream discharge with adjacent valley structure have been rare, resulting in poor understanding of the distribution of the hydrologic mechanisms that cause variability in streamflow generation along valleys. This study explores detailed surveys of stream base flow across a gauged, 23 km 2 mountain watershed. Research objectives were (1) to relate spatial variability in base flow to fundamental elements of watershed structure, primarily topographic contributing area, and (2) to assess temporal changes in the spatial patterns of those relationships during a seasonal base flow recession. We analyzed spatiotemporal variability in base flow using (1) summer hydrographs at the study watershed outlet and 5 subwatershed outlets and (2) longitudinal series of discharge measurements every ~100 m along the streams of the 3 largest subwatersheds (1200 to 2600 m in valley length), repeated 2 to 3 times during base flow recession. Reaches within valley segments of 300 to 1200 m in length tended to demonstrate similar streamflow generation characteristics. Locations of transitions between these segments were consistent throughout the recession, and tended to be collocated with abrupt longitudinal transitions in valley slope or hillslope-riparian characteristics. Both within and among subwatersheds, correlation between the spatial distributions of streamflow and topographic contributing area decreased during the recession, suggesting a general decrease in the influence of topography on stream base flow contributions. As topographic controls on base flow evidently decreased, multiple aspects of subsurface structure were likely to have gained influence.

Water Resources Research↗
USGS2022

Integrated modeling of dynamic marsh feedbacks and evolution under sea-level rise in a mesotidal estuary (Plum Island, MA, USA)

Around the world, wetland vulnerability to sea-level rise (SLR) depends on different factors including tidal regimes, topography, creeks and estuary geometry, sediment availability, vegetation type, etc. The Plum Island estuary (PIE) is a mesotidal wetland system on the east coast of the United States. This research applied a newly updated Hydro-MEM (integrated hydrodynamic-marsh) model to assess the impacts of intermediate-low (50 cm), intermediate (1 m), and intermediate-high (1.5 m) SLR on marsh evolution by the year 2100. Model advancements include capturing vegetation change, inorganic and below and aboveground organic matter portion of marsh platform accretion, and mudflat creation. Although the results indicate a low vulnerability marsh at the PIE, the vegetation changes from high to low marsh under all SLR scenarios (2%–22%), with the higher bounds belonging to higher rise scenarios. Lower SLR produces more productive marsh (13% gain in high productivity regions), whereas the highest SLR scenario causes increased tidal inundation, which leads to loss in productivity (12% change from high to low productivity regions), generation of mudflats (17% of the domain land), and marsh migration to higher lands. Sensitive nonlinear tidal flow changes, which may be increased or decreased with SLR as a result of mudflat creation, marsh migration, and bottom friction change, emphasize the importance of integrated modeling approaches that include dynamic marsh feedbacks in hydrodynamic modeling and varying hydrodynamic effects on the marsh system.

Massachusetts↗
USGS2013

Incorporating probabilistic seasonal climate forecasts into river management using a risk-based framework

[1] Despite the influence of hydroclimate on river ecosystems, most efforts to date have focused on using climate information to predict streamflow for water supply. However, as water demands intensify and river systems are increasingly stressed, research is needed to explicitly integrate climate into streamflow forecasts that are relevant to river ecosystem management. To this end, we present a five step risk-based framework: (1) define risk tolerance, (2) develop a streamflow forecast model, (3) generate climate forecast ensembles, (4) estimate streamflow ensembles and associated risk, and (5) manage for climate risk. The framework is successfully demonstrated for an unregulated watershed in southwest Montana, where the combination of recent drought and water withdrawals has made it challenging to maintain flows needed for healthy fisheries. We put forth a generalized linear modeling (GLM) approach to develop a suite of tools that skillfully model decision-relevant low flow characteristics in terms of climate predictors. Probabilistic precipitation forecasts are used in conjunction with the GLMs, resulting in season-ahead prediction ensembles that provide the full risk profile. These tools are embedded in an end-to-end risk management framework that directly supports proactive fish conservation efforts. Results show that the use of forecasts can be beneficial to planning, especially in wet years, but historical precipitation forecasts are quite conservative (i.e., not very “sharp”). Synthetic forecasts show that a modest “sharpening” can strongly impact risk and improve skill. We emphasize that use in management depends on defining relevant environmental flows and risk tolerance, requiring local stakeholder involvement.

Montana↗
USGS2021

Improving remotely sensed river bathymetry by image-averaging

Basic data on river bathymetry is critical for numerous applications in river research and management and is increasingly obtained via remote sensing, but the noisy, pixelated appearance of image‐derived depth maps can compromise subsequent analyses. We hypothesized that this noise originates from reflectance from an irregular water surface and introduced a framework for mitigating these effects by Inferring Bathymetry from Averaged River Images (IBARI). This workflow produces time‐averaged images from video frames stabilized to account for platform motion and/or computes a spatial average from an ensemble simulated by randomly shifting images relative to themselves. We used field observations of water depth and helicopter‐based videos from a clear‐flowing river to assess the potential of this approach to improve depth retrieval. Our results indicated that depths inferred from averaged images were more accurate and precise than those inferred from single frames; observed versus predicted regression R 2 increased from 0.80 to 0.88. In addition, IBARI significantly enhanced the texture of image‐derived depth maps, leading to smoother, more coherent representations of channel morphology. Depth retrieval improved with image sequence duration, but the number of images was more important than the length of time encompassed; shorter acquisitions at higher frame rates would economize data collection. We also demonstrated the potential to scale up the IBARI workflow by producing a mosaic of bathymetric maps derived from averaged images acquired at several hovering waypoints distributed along a 2.36 km reach. This approach is well‐suited to data collected from helicopters and small unmanned aircraft systems.

Water Resources Research↗
USGS1998

Interactions between pool geometry and hydraulics

An experimental and computational research approach was used to determine interactions between pool geometry and hydraulics. A 20-m-long, 1.8-m-wide flume was used to investigate the effect of four different geometric aspects of pool shape on flow velocity. Plywood sections were used to systematically alter constriction width, pool depth, pool length, and pool exit-slope gradient, each at two separate levels. Using the resulting 16 unique geometries with measured pool velocities in four-way factorial analyses produced an empirical assessment of the role of the four geometric aspects on the pool flow patterns and hence the stability of the pool. To complement the conclusions of these analyses, a two-dimensional computational flow model was used to investigate the relationships between pool geometry and flow patterns over a wider range of conditions. Both experimental and computational results show that constriction and depth effects dominate in the jet section of the pool and that pool length exhibits an increasing effect within the recirculating-eddy system. The pool exit slope appears to force flow reattachment. Pool length controls recirculating-eddy length and vena contracta strength. In turn, the vena contracta and recirculating eddy control velocities throughout the pool.

Water Resources Research↗