USGS Science⌕ Search

SEARCH · USGS Science

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 757 records · Page 42Linked to original sources

Historical effects of El Nino and La Nina events on the seasonal evolution of the montane snowpack in the Columbia and Colorado River Basins

Snow‐water equivalent (SWE) data measured at several hundred montane sites in the western United States are used to examine the historic effects of El Nino and La Nina events on seasonal snowpack evolution in the major subbasins in the Columbia and Colorado River systems. Results are used to predict annual runoff. In the Columbia River Basin, there is a general tendency for decreased SWE during El Nino years and increased SWE in La Nina years. However, the SWE anomalies for El Nino years are much less pronounced. This occurs in part because midlatitude circulation anomalies in El Nino years are located 35° east of those in La Nina years. This eastward shift is most evident in midwinter, at which time, SWE anomalies associated with El Nino are actually positive in coastal regions of the Columbia River Basin. In the Colorado River Basin, mean anomalies in SWE and annual runoff during El Nino years depict a transition between drier‐than‐average conditions in the north, and wetter‐than‐average conditions in the southwest. Associations during La Nina years are generally opposite those in El Nino years. SWE anomalies tend to be more pronounced in spring in the Lower Colorado River Basin. Our predictions of runoff reveal modest skill for scenarios using only historic El Nino and La Nina information. Predictions based on the water stored in the seasonal snowpack are, in almost all cases, much higher than those based on El Nino‐Southern Oscillation (ENSO) information alone. However, combining observed midwinter snow conditions with information on seasonal snowpack evolution associated with ENSO improves predictions for basins in which ENSO signals exhibit strong seasonality.

Water Resources Research↗

A minimalist probabilistic description of root zone soil water

The probabilistic response of depth‐integrated soil water to given climatic forcing can be described readily using an existing supply‐demand‐storage model. An apparently complex interaction of numerous soil, climate, and plant controls can be reduced to a relatively simple expression for the equilibrium probability density function of soil water as a function of only two dimensionless parameters. These are the index of dryness (ratio of mean potential evaporation to mean precipitation) and a dimensionless storage capacity (active root zone soil water capacity divided by mean storm depth). The first parameter is mainly controlled by climate, with surface albedo playing a subsidiary role in determining net radiation. The second is a composite of soil (through moisture retention characteristics), vegetation (through rooting characteristics), and climate (mean storm depth). This minimalist analysis captures many essential features of a more general probabilistic analysis, but with a considerable reduction in complexity and consequent elucidation of the critical controls on soil water variability. In particular, it is shown that (1) the dependence of mean soil water on the index of dryness approaches a step function in the limit of large soil water capacity; (2) soil water variance is usually maximized when the index of dryness equals 1, and the width of the peak varies inversely with dimensionless storage capacity; (3) soil water has a uniform probability density function when the index of dryness is 1 and the dimensionless storage capacity is large; and (4) the soil water probability density function is bimodal if and only if the index of dryness is <1, but this bimodality is pronounced only for artificially small values of the dimensionless storage capacity.

Water Resources Research↗

Interbasin groundwater flow in south central Nevada: A further comment on the discussion between Davisson et. al.. [1999a, 1999b] and Thomas [1999]

In their response to the comments by Thomas [1999], Davisson et al. [1999a] dismiss a large set of potentiometric measurements pertinent to an understanding of the hydrogeology of Yucca and Frenchman Flats in southcentral Nevada. This commentary is submitted to demonstrate, first, that their dismissal of this data set is unfounded and, second, that these potentiometric data call into question the central thesis of the original paper by Davisson et al. [1999b].

Nevada↗

Seasonal subsidence and rebound in Las Vegas Valley, Nevada, observed by Synthetic Aperture Radar Interferometry

Analyses of areal variations in the subsidence and rebound occurring over stressed aquifer systems, in conjunction with measurements of the hydraulic head fluctuations causing these displacements, can yield valuable information about the compressibility and storage properties of the aquifer system. Historically, stress‐strain relationships have been derived from paired extensometer/piezometer installations, which provide only point source data. Because of the general unavailability of spatially detailed deformation data, areal stress‐strain relations and their variability are not commonly considered in constraining conceptual and numerical models of aquifer systems. Interferometric synthetic aperture radar (InSAR) techniques can map ground displacements at a spatial scale of tens of meters over 100 km wide swaths. InSAR has been used previously to characterize larger magnitude, generally permanent aquifer system compaction and land subsidence at yearly and longer timescales, caused by sustained drawdown of groundwater levels that produces intergranular stresses consistently greater than the maximum historical stress. We present InSAR measurements of the typically small‐magnitude, generally recoverable deformations of the Las Vegas Valley aquifer system occurring at seasonal timescales. From these we derive estimates of the elastic storage coefficient for the aquifer system at several locations in Las Vegas Valley. These high‐resolution measurements offer great potential for future investigations into the mechanics of aquifer systems and the spatial heterogeneity of aquifer system structure and material properties as well as for monitoring ongoing aquifer system compaction and land subsidence.

Nevada↗

Development of channel organization and roughness following sediment pulses in single‐thread, gravel bed rivers

Large, episodic inputs of coarse sediment (sediment pulses) in forested, mountain streams may result in changes in the size and arrangement of bed forms and in channel roughness. A conceptual model of channel organization delineates trajectories of response to sediment pulses for many types of gravel bed channels. Channels exhibited self‐organizing behavior to various degrees based on channel gradient, presence of large in‐channel wood or other forcing elements, the size of the sediment pulse, and the number of bed‐mobilizing flows since disturbance. Typical channel changes following a sediment pulse were initial decreases in water depth, in variability of bed elevations, and in the regularity of bed form spacing. Trajectories of change subsequently showed increased average water depth, more variable and complex bed topography, and increased uniformity of bed form spacing. Bed form spacing in streams with abundant forcing elements developed at a shorter spatial scale (two to five channel widths) than in streams without such forcing mechanisms (five to 10 channel widths). Channel roughness increased as bed forms developed.

Water Resources Research↗

Predictive modeling of flow and transport in a two‐dimensional intermediate‐scale, heterogeneous porous medium

As a first step toward understanding the role of sedimentary structures in flow and transport through porous media, this work deterministically examines how small‐scale laboratory‐measured values of hydraulic conductivity relate to in situ values of simple, artificial structures in an intermediate‐scale (10 m long), two‐dimensional, heterogeneous, laboratory experiment. Results were judged based on how well simulations using measured values of hydraulic conductivities matched measured hydraulic heads, net flow, and transport through the tank. Discrepancies were investigated using sensitivity analysis and nonlinear regression estimates of the in situ hydraulic conductivity that produce the best fit to measured hydraulic heads and net flow. Permeameter and column experiments produced laboratory measurements of hydraulic conductivity for each of the sands used in the intermediate‐scale experiments. Despite explicit numerical representation of the heterogeneity the laboratory‐measured values underestimated net flow by 12–14% and were distinctly smaller than the regression‐estimated values. The significance of differences in measured hydraulic conductivity values was investigated by comparing variability of transport predictions using the different measurement methods to that produced by different realizations of the heterogeneous distribution. Results indicate that the variations in measured hydraulic conductivity were more important to transport than variations between realizations of the heterogeneous distribution of hydraulic conductivity.

Water Resources Research↗

In‐stream sorption of fulvic acid in an acidic stream: A stream‐scale transport experiment

The variation of concentration and composition of dissolved organic carbon (DOC) in stream waters cannot be explained solely on the basis of soil processes in contributing subcatchments. To investigate in‐stream processes that control DOC, we injected DOC‐enriched water into a reach of the Snake River (Summit County, Colorado) that has abundant iron oxyhydroxides coating the streambed. The injected water was obtained from the Suwannee River (Georgia), which is highly enriched in fulvic acid. The fulvic acid from this water is the standard reference for aquatic fulvic acid for the International Humic Substances Society and has been well characterized. During the experimental injection, significant removal of sorbable fulvic acid occurred within the first 141 m of stream reach. We coinjected a conservative tracer (lithium chloride) and analyzed the results with the one‐dimensional transport with inflow and storage (OTIS) stream solute transport model to quantify the physical transport mechanisms. The downstream transport of fulvic acid as indicated by absorbance was then simulated using OTIS with a first‐order kinetic sorption rate constant applied to the sorbable fulvic acid. The “sorbable” fraction of injected fulvic acid was irreversibly sorbed by streambed sediments at rates (kinetic rate constants) of the order of 10 −4 –10 −3 s −1 . In the injected Suwannee River water, sorbable and nonsorbable fulvic acid had distinct chemical characteristics identified in 13 C‐NMR spectra. The 13 C‐NMR spectra indicate that during the experiment, the sorbable “signal” of greater aromaticity and carboxyl content decreased downstream; that is, these components were preferentially removed. This study illustrates that interactions between the water and the reactive surfaces will modify significantly the concentration and composition of DOC observed in streams with abundant chemically reactive surfaces on the streambed and in the hyporheic zone.

Water Resources Research↗

Analysis of δ15N and δ18O to differentiate NO3− sources in runoff at two watersheds in the Catskill Mountains of New York

To quantify the movement of atmospheric nitrogen deposition through two forested watersheds in the Catskill Mountains of New York, dual‐isotope analysis (δ 15 N and δ 18 O) was used to differentiate NO 3 − derived from precipitation from NO 3 − derived by microbial nitrification and to quantify the contributions of these sources to NO 3 − in drainage waters. Samples of stream water, soil water, precipitation, snowmelt, and O‐horizon soil were collected during the March and April snowmelt period of 1994 and throughout an 18‐month period from August 1995 through February 1997. The mean δ 18 O‐NO 3 − value of precipitation was +50.5‰, whereas the mean values for stream water and soil water were +17.7‰ and +23.6‰, respectively. The mean δ 15 N‐NO 3 − of precipitation was −0.2‰, that of soil water was +1.4‰, and that of stream water was +2.3‰; these values showed greater overlap among the three different waters than did the δ 18 O‐NO 3 − values, indicating that δ 15 N‐NO 3 − was not as useful for source separation. Soil water δ 18 O‐NO 3 − values decreased, and δ 15 N‐NO 3 − values increased, from the O to the B and C horizons, but most of the differences among horizons were not statistically significant. Nitrate derived by nitrification in incubated soil samples had a wide range of δ 15 N‐NO 3 − values, from +1.5‰ to +16.1‰, whereas δ 18 O‐NO 3 − values ranged more narrowly, from +13.2‰ to +16.0‰. Values of δ 18 O‐NO 3 − indicated that NO 3 − in stream water is mainly derived from nitrification. Only during a high‐flow event that exceeded the annual flood was precipitation a major contributor to stream water NO 3 − . Values of δ 18 O‐NO 3 − and δ 15 N‐NO 3 − changed at differing rates as NO 3 − cycled through these watersheds because δ 18 O‐NO 3 − values change sharply through the incorporation of oxygen from ambient water and gas during nitrification, whereas δ 15 N‐NO 3 − values change only incrementally through fractionation during biocycling processes. The results of this study show that most NO 3 − is first cycled through the biota and nitrified before entering the stream.

Water Resources Research↗

Pathways for nitrate release from an alpine watershed: Determination using δ15N and δ18O

Snowpack, snowmelt, precipitation, surface water, and groundwater samples from the Loch Vale watershed in Colorado were analyzed for δ 15 N and δ 18 O of nitrate to determine the processes controlling the release of atmospherically deposited nitrogen from alpine and subalpine ecosystems. Although overlap was found between the δ 15 N (NO3) values for all water types (−4 to +6‰), the δ 18 O (NO3) values for surface water and groundwater (+10 to +30‰) were usually distinct from snowpack, snowmelt, and rainfall values (+40 to +70‰). During snowmelt, δ 18 O (NO3) indicated that about half of the nitrate in stream water was the product of microbial nitrification; at other times that amount was greater than half. Springs emerging from talus deposits had high nitrate concentrations and a seasonal pattern in δ 18 O (NO3) that was similar to the pattern in the streams, indicating that shallow groundwater in talus deposits is a likely source of stream water nitrate. Only a few samples of surface water and groundwater collected during early snowmelt and large summer rain events had isotopic compositions that indicated most of the nitrate came directly from atmospheric deposition with no biological assimilation and release. This study demonstrates the value of the nitrate double‐isotope technique for determining nitrogen‐cycling processes and sources of nitrate in small, undisturbed watersheds that are enriched with inorganic nitrogen.

Water Resources Research↗

Dynamic interpretation of slug tests in highly permeable aquifers

Considerable progress has been made in developing a theoretical framework for modeling slug test responses in formations with high hydraulic conductivity K . However, several questions of practical significance remain unresolved. Given the rapid and often oscillatory nature of test responses, the traditional hydrostatic relationship between the water level and the transducer‐measured head in the water column may not be appropriate. A general dynamic interpretation is proposed that describes the relationship between water level response and transducer‐measured head. This theory is utilized to develop a procedure for transforming model‐generated water level responses to transducer readings. The magnitude of the difference between the actual water level position and the apparent position based on the transducer measurement is a function of the acceleration and velocity of the water column, test geometry, and depth of the transducer. The dynamic approach explains the entire slug test response, including the often‐noted discrepancy between the actual initial water level displacement and that measured by a transducer in the water column. Failure to use this approach can lead to a significant underestimation of K when the transducer is a considerable distance below the static water level. Previous investigators have noted a dependence of test responses on the magnitude of the initial water level displacement and have developed various approximate methods for analyzing such data. These methods are re‐examined and their limitations clarified. Practical field guidelines are proposed on the basis of findings of this work. The soundness of the dynamic approach is demonstrated through a comparison of K profiles from a series of multilevel slug tests with those from dipole‐flow tests performed in the same wells.

Water Resources Research↗

Natural and management influences on freshwater inflows and salinity in the San Francisco Estuary at monthly to interannual scales

Understanding the processes controlling the physics, chemistry, and biology of the San Francisco Estuary and their relation to climate variability is complicated by the combined influence on freshwater inflows of natural variability and upstream management. To distinguish these influences, alterations of estuarine inflow due to major reservoirs and freshwater pumping in the watershed were inferred from available data. Effects on salinity were estimated by using reconstructed estuarine inflows corresponding to differing levels of impairment to drive a numerical salinity model. Both natural and management inflow and salinity signals show strong interannual variability. Management effects raise salinities during the wet season, with maximum influence in spring. While year‐to‐year variations in all signals are very large, natural interannual variability can greatly exceed the range of management effects on salinity in the estuary.

California↗

Interparticle collision of natural sediment grains in water

Elastohydrodynamic theory and measurements of particle impacts on an inclined glass plane in water are used to investigate the mechanics of interparticle collisions in sediment‐transporting flows. A collision Stokes number is proposed as a measure of the momentum of an interparticle collision versus the viscous pressure force in the interstitial gap between colliding particles. The viscous pressure force opposes motion of the particles on approach and rebound. A Stokes number of between 39 and 105 is estimated as the critical range below which particle impacts are completely viscously damped and above which impacts are partially elastic. The critical Stokes number is shown to roughly coincide with the Bagnold number transition between macroviscous and grain inertial debris flows and the transition between damped and partially elastic bed load transport saltation impacts. The nonspherical nature of natural particles significantly alters the motion of the center of mass after a partially elastic collision. The normal to the point of contact between the particles does not necessarily go through the center of mass. Thus normal rebound of the center of mass may not occur. A model of particle motion after rebound for particles of arbitrary shape, conserving both linear and angular momentum, is proposed.

Water Resources Research↗

Denitrification in the recharge area and discharge area of a transient agricultural nitrate plume in a glacial outwash sand aquifer, Minnesota

Recharge rates of nitrate (NO 3 − ) to groundwater beneath agricultural land commonly are greater than discharge rates of NO 3 − in nearby streams, but local controls of NO 3 − distribution in the subsurface generally are poorly known. Groundwater dating (CFC, 3 H) was combined with chemical (ions and gases) and stable isotope (N, S, and C) analyses to resolve the effects of land use changes, flow patterns, and water‐aquifer reactions on the distributions of O 2 , NO 3 − , SO 4 = , and other constituents in a two‐dimensional vertical section leading from upland cultivated fields to a riparian wetland and stream in a glacial outwash sand aquifer near Princeton, Minnesota. Within this section a “plume” of oxic NO 3 − ‐rich groundwater was present at shallow depths beneath the fields and part of the wetland but terminated before reaching the stream or the wetland surface. Groundwater dating and hydraulic measurements indicate travel times in the local flow system of 0 to >40 years, with stratified recharge beneath the fields, downward diversion of the shallow NO 3 − ‐bearing plume by semiconfining organic‐rich valley‐filling sediments under the wetland and upward discharge across the valley and stream bottom. The concentrations and δ 15 N values of NO 3 − and N 2 indicate that the NO 3 − plume section was bounded in three directions by a curvilinear zone of active denitrification that limited its progress; however, when recalculated to remove the effects of denitrification, the data also indicate changes in both the concentrations and δ 15 N values of NO 3 − that was recharged in the past. Isotope data and mass balance calculations indicate that FeS 2 and other ferrous Fe phases were the major electron donors for denitrification in at least two settings: (1) within the glacial‐fluvial aquifer sediments beneath the recharge and discharge areas and (2) along the bottom of the valley‐filling sediments in the discharge area. Combined results indicate that the shape and progress of the oxic NO 3 − plume termination were controlled by a combination of (1) historical and spatial variations in land use practices, (2) contrast in groundwater flow patterns between the agricultural recharge area and riparian wetland discharge area, and (3) distribution and abundance of electron donors in both the sand aquifer and valley‐filling sediments. The data are consistent with slow migration of redox zones through the aquifer in response to recharging oxic groundwater during Holocene time, then an order‐of‐magnitude increase in the flux of electron acceptors as a result of agricultural NO 3 − contamination in the late twentieth century, to which the redox zone configuration still may be adjusting. The importance of denitrification for NO 3 − movement through formerly glaciated terrains should depend on the source areas and depositional environments of the glacial sediments, as well as geomorphology and recent stream‐valley sediment history.

Water Resources Research↗

Extracting low‐resolution river networks from high‐resolution digital elevation models

Including a global river network in the land component of global climate models (GCMs) is necessary in order to provide a more complete representation of the hydrologic cycle. The process of creating these networks is called river network upscaling and consists of lowering the resolution of already available fine networks to make them compatible with GCMs. Fine‐resolution river networks have a level of detail appropriate for analysis at the watershed scale but are too intensive for global hydrologic studies. A river network upscaling algorithm, which processes fine‐resolution digital elevation models to determine the flow directions that best describe the flow patterns in a coarser user‐defined scale, is presented. The objectives of this study were to develop an algorithm that advances the previous work in the field by being applicable at a global scale, allowing for the upscaling to be performed in a projected environment, and generating evenly distributed flow directions.

Water Resources Research↗

Evaluation of unsaturated zone water fluxes in heterogeneous alluvium at a Mojave Basin Site

Vertical and horizontal water fluxes in the unsaturated zone near intermittent streams critically affect ecosystems, water supply, and contaminant transport in arid and semiarid regions. The subsurface near the Oro Grande Wash is typical in having great textural diversity, pronounced layer contrasts, and extremely low hydraulic conductivities associated with nearly dry media. These features prevent a straightforward application of the Darcian method for recharge estimation, which has provided high‐quality flux estimates at simpler, wetter sites. We have augmented the basic Darcian method with theoretical developments such that a small number of core sample unsaturated hydraulic property measurements, combined with additional, easily obtained data (e.g., drillers' logs) can provide useful flux estimates and knowledge of two‐dimensional water behavior beneath the wash.

Water Resources Research↗

Macroscale water fluxes 3. Effects of land processes on variability of monthly river discharge

A salient characteristic of river discharge is its temporal variability. The time series of flow at a point on a river can be viewed as the superposition of a smooth seasonal cycle and an irregular, random variation. Viewing the random component in the spectral domain facilitates both its characterization and an interpretation of its major physical controls from a global perspective. The power spectral density functions of monthly flow anomalies of many large rivers worldwide are typified by a “red noise” process: the density is higher at low frequencies (e.g., <1 y −1 ) than at high frequencies, indicating disproportionate (relative to uncorrelated “white noise”) contribution of low frequencies to variability of monthly flow. For many high‐latitude and arid‐region rivers, however, the power is relatively evenly distributed across the frequency spectrum. The power spectrum of monthly flow can be interpreted as the product of the power spectrum of monthly basin total precipitation (which is typically white or slightly red) and several filters that have physical significance. The filters are associated with (1) the conversion of total precipitation (sum of rainfall and snowfall) to effective rainfall (liquid flux to the ground surface from above), (2) the conversion of effective rainfall to soil water excess (runoff), and (3) the conversion of soil water excess to river discharge. Inferences about the roles of each filter can be made through an analysis of observations, complemented by information from a global model of the ocean‐atmosphere‐land system. The first filter causes a snowmelt‐related amplification of high‐frequency variability in those basins that receive substantial snowfall. The second filter causes a relatively constant reduction in variability across all frequencies and can be predicted well by means of a semiempirical water balance relation. The third filter, associated with groundwater and surface water storage in the river basin, causes a strong reduction in high‐frequency variability of many basins. The strength of this reduction can be quantified by an average residence time of water in storage, which is typically on the order of 20–50 days. The residence time is demonstrably influenced by freezing conditions in the basin, fractional cover of the basin by lakes, and runoff ratio (ratio of mean runoff to mean precipitation). Large lake areas enhance storage and can greatly increase total residence times (100 to several hundred days). Freezing conditions appear to cause bypassing of subsurface storage, thus reducing residence times (0–30 days). Small runoff ratios tend to be associated with arid regions, where the water table is deep, and consequently, most of the runoff is produced by processes that bypass the saturated zone, leading to relatively small residence times for such basins (0–40 days).

Water Resources Research↗

Deep arid system hydrodynamics 1. Equilibrium states and response times in thick desert vadose zones

Quantifying moisture fluxes through deep desert soils remains difficult because of the small magnitude of the fluxes and the lack of a comprehensive model to describe flow and transport through such dry material. A particular challenge for such a model is reproducing both observed matric potential and chloride profiles. We propose a conceptual model for flow in desert vadose zones that includes isothermal and nonisothermal vapor transport and the role of desert vegetation in supporting a net upward moisture flux below the root zone. Numerical simulations incorporating this conceptual model match typical matric potential and chloride profiles. The modeling approach thereby reconciles the paradox between the recognized importance of plants, upward driving forces, and vapor flow processes in desert vadose zones and the inadequacy of the downward‐only liquid flow assumption of the conventional chloride mass balance approach. Our work shows that water transport in thick desert vadose zones at steady state is usually dominated by upward vapor flow and that long response times, of the order of 10 4 –10 5 years, are required to equilibrate to existing arid surface conditions. Simulation results indicate that most thick desert vadose zones have been locked in slow drying transients that began in response to a climate shift and establishment of desert vegetation many thousands of years ago.

Water Resources Research↗

Deep arid system hydrodynamics 2. Application to paleohydrologic reconstruction using vadose zone profiles from the northern Mojave Desert

Site‐specific numerical modeling of four sites in two arid alluvial basins within the Nevada Test Site employs a conceptual model of deep arid system hydrodynamics that includes vapor transport, the role of xeric vegetation, and long‐term surface boundary transients. Surface boundary sequences, spanning 110 kyr, that best reproduce measured chloride concentration and matric potential profiles from four deep (230–460 m) boreholes concur with independent paleohydrologic and paleoecological records from the region. Simulations constrain a pluvial period associated with infiltration of 2–5 mm yr −1 at 14–13 ka and denote a shift linked to the establishment of desert vegetation at 13–9.5 ka. Retrodicted moisture flux histories inferred from modeling results differ significantly from those determined using the conventional chloride mass balance approach that assumes only downward advection. The modeling approach developed here represents a significant advance in the use of deep vadose zone profile data from arid regions to recover detailed paleohydrologic and current hydrologic information.

Mojave Desert↗