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 289 records · Page 16Linked to original sources
USGS1999

A new multistage groundwater transport inverse method: presentation, evaluation, and implications

More computationally efficient methods of using concentration data are needed to estimate groundwater flow and transport parameters. This work introduces and evaluates a three‐stage nonlinear‐regression‐based iterative procedure in which trial advective‐front locations link decoupled flow and transport models. Method accuracy and efficiency are evaluated by comparing results to those obtained when flow‐ and transport‐model parameters are estimated simultaneously. The new method is evaluated as conclusively as possible by using a simple test case that includes distinct flow and transport parameters, but does not include any approximations that are problem dependent. The test case is analytical; the only flow parameter is a constant velocity, and the transport parameters are longitudinal and transverse dispersivity. Any difficulties detected using the new method in this ideal situation are likely to be exacerbated in practical problems. Monte‐Carlo analysis of observation error ensures that no specific error realization obscures the results. Results indicate that, while this, and probably other, multistage methods do not always produce optimal parameter estimates, the computational advantage may make them useful in some circumstances, perhaps as a precursor to using a simultaneous method.

Water Resources Research↗
USGS2010

The relationship between soil heterotrophic activity, soil dissolved organic carbon (DOC) leachate, and catchment-scale DOC export in headwater catchments

Dissolved organic carbon (DOC) from terrestrial sources forms the major component of the annual carbon budget in many headwater streams. In high-elevation catchments in the Rocky Mountains, DOC originates in the upper soil horizons and is flushed to the stream primarily during spring snowmelt. To identify controls on the size of the mobile soil DOC pool available to be transported during the annual melt event, we measured soil DOC production across a range of vegetation communities and soil types together with catchment DOC export in paired watersheds in Summit County, Colorado. Both surface water DOC concentrations and watershed DOC export were lower in areas where pyrite weathering resulted in lower soil p H. Similarly, the amount of DOC leached from organic soils was significantly smaller ( p < 0.01) at sites having low soil p H. Scaling point source measurements of DOC production and leaching to the two basins and assuming only vegetated areas contribute to DOC production, we calculated that the amount of mobile DOC available to be leached to surface water during melt was 20.3 g C m −2 in the circumneutral basin and 17.8 g C m −2 in the catchment characterized by pyrite weathering. The significant ( r 2 =0.91 and p < 0.05), linear relationship between over-winter CO 2 flux and the amount of DOC leached from upper soil horizons during snowmelt suggests that the mechanism for the difference in production of mobile DOC was heterotrophic processing of soil carbon in snow-covered soil. Furthermore, this strong relationship between over-winter heterotrophic activity and the size of the mobile DOC pool present in a range of soil and vegetation types provides a likely mechanism for explaining the interannual variability of DOC export observed in high-elevation catchments.

Water Resources Research↗
USGS1999

Denitrification in marine shales in northeastern Colorado

Parts of the South Platte River alluvial aquifer in northeastern Colorado are underlain by the Pierre Shale, a marine deposit of Late Cretaceous age that is <1000 m thick. Ground water in the aquifer is contaminated with NO 3 ‐ , and the shale contains abundant potential electron donors for denitrification in the forms of organic carbon and sulfide minerals. Nested piezometers were sampled, pore water was squeezed from cores of shale, and an injection test was conducted to determine if denitrification in the shale was a sink for alluvial NO 3 − and to measure denitrification rates in the shale. Measured values of NO 3 − , N 2 , NH 4 + , δ 15 N[NO 3 − ], δ 15 N[N 2 ], and δ 15 N[NH 4 + ] in the alluvial and shale pore water indicated that denitrification in the shale was a sink for alluvial NO 3 − . Chemical gradients, reaction rate constants, and hydraulic head data indicated that denitrification in the shale was limited by the slow rate of NO 3 − transport (possibly by diffusion) into the shale. The apparent in situ first‐order rate constant for denitrification in the shale based on diffusion calculations was of the order of 0.04–0.4 yr −1 , whereas the potential rate constant in the shale based on injection tests was of the order of 60 yr −1 . Chemical data and mass balance calculations indicate that organic carbon was the primary electron donor for denitrification in the shale during the injection test, and ferrous iron was a minor electron donor in the process. Flux calculations for the conditions encountered at the site indicate that denitrification in the shale could remove only a small fraction of the annual agricultural NO 3 ‐ input to the alluvial aquifer. However, the relatively large potential first‐order rate constant for denitrification in the shale indicated that the percentage of NO 3 − uptake by the shale could be considerably larger in areas where NO 3 − advection.

Colorado↗
USGS2010

A procedure for classifying textural facies in gravel‐bed rivers

Textural patches (i.e., grain‐size facies) are commonly observed in gravel‐bed channels and are of significance for both physical and biological processes at subreach scales. We present a general framework for classifying textural patches that allows modification for particular study goals, while maintaining a basic degree of standardization. Textures are classified using a two‐tier system of ternary diagrams that identifies the relative abundance of major size classes and subcategories of the dominant size. An iterative procedure of visual identification and quantitative grain‐size measurement is used. A field test of our classification indicates that it affords reasonable statistical discrimination of median grain size and variance of bed‐surface textures. We also explore the compromise between classification simplicity and accuracy. We find that statistically meaningful textural discrimination requires use of both tiers of our classification. Furthermore, we find that simplified variants of the two‐tier scheme are less accurate but may be more practical for field studies which do not require a high level of textural discrimination or detailed description of grain‐size distributions. Facies maps provide a natural template for stratifying other physical and biological measurements and produce a retrievable and versatile database that can be used as a component of channel monitoring efforts.

Water Resources Research↗
USGS1999

Stochastic analysis of virus transport in aquifers

A large-scale model of virus transport in aquifers is derived using spectral perturbation analysis. The effects of spatial variability in aquifer hydraulic conductivity and virus transport (attachment, detachment, and inactivation) parameters on large-scale virus transport are evaluated. A stochastic mean model of virus transport is developed by linking a simple system of local-scale free-virus transport and attached-virus conservation equations from the current literature with a random-field representation of aquifer and virus transport properties. The resultant mean equations for free and attached viruses are found to differ considerably from the local-scale equations on which they are based and include effects such as a free-virus effective velocity that is a function of aquifer heterogeneity as well as virus transport parameters. Stochastic mean free-virus breakthrough curves are compared with local model output in order to observe the effects of spatial variability on mean one-dimensional virus transport in three-dimensionally heterogeneous porous media. Significant findings from this theoretical analysis include the following: (1) Stochastic model breakthrough occurs earlier than local model breakthrough, and this effect is most pronounced for the least conductive aquifers studied. (2) A high degree of aquifer heterogeneity can lead to virus breakthrough actually preceding that of a conservative tracer. (3) As the mean hydraulic conductivity is increased, the mean model shows less sensitivity to the variance of the natural-logarithm hydraulic conductivity and mean virus diameter. (4) Incorporation of a heterogeneous colloid filtration term results in higher predicted concentrations than a simple first-order adsorption term for a given mean attachment rate. (5) Incorporation of aquifer heterogeneity leads to a greater range of virus diameters for which significant breakthrough occurs. (6) The mean model is more sensitive to the inactivation rate of viruses associated with solid surfaces than to the inactivation rate of viruses in solution.

Water Resources Research↗
USGS1999

Natural attenuation of chlorinated volatile organic compounds in a freshwater tidal wetland: Field evidence of anaerobic biodegradation

Field evidence collected along two groundwater flow paths shows that anaerobic biodegradation naturally attenuates a plume of chlorinated volatile organic compounds as it discharges from an aerobic sand aquifer through wetland sediments. A decrease in concentrations of two parent contaminants, trichloroethylene (TCE) and 1,1,2,2‐tetrachloroethane (PCA), and a concomitant increase in concentrations of anaerobic daughter products occurs along upward flow paths through the wetland sediments. The daughter products 1,2‐dichloroethylene, vinyl chloride, 1,1,2‐trichloroethane, and 1,2‐dichloroethane are produced from hydrogenolysis of TCE and from PCA degradation through hydrogenolysis and dichloroelimination (reductive dechlorination) pathways. Total concentrations of TCE, PCA, and their degradation products, however, decrease to below detection levels within 0.15–0.30 m of land surface. The enhanced reductive dechlorination of TCE and PCA in the wetland sediments is associated with the naturally higher concentrations of dissolved organic carbon and the lower redox state of the groundwater compared to the aquifer. This field study indicates that wetlands and similar organic‐rich environments at groundwater/surface‐water interfaces may be important in intercepting groundwater contaminated with chlorinated organics and in naturally reducing concentrations and toxicity before sensitive surface‐water receptors are reached.

Water Resources Research↗
USGS1999

Soil water retention and maximum capillary drive from saturation to oven dryness

This paper provides an alternative method to describe the water retention curve over a range of water contents from saturation to oven dryness. It makes two modifications to the standard Brooks and Corey [1964] (B-C) description, one at each end of the suction range. One expression proposed by Rossi and Nimmo [1994] is used in the high-suction range to a zero residual water content. (This Rossi-Nimmo modification to the Brooks-Corey model provides a more realistic description of the retention curve at low water contents.) Near zero suction the second modification eliminates the region where there is a change in suction with no change in water content. Tests on seven soil data sets, using three distinct analytical expressions for the high-, medium-, and low-suction ranges, show that the experimental water retention curves are well fitted by this composite procedure. The high-suction range of saturation contributes little to the maximum capillary drive, defined with a good approximation for a soil water and air system as H cM = ∫ 0 ∞ K rw dh c , where k rw is relative permeability (or conductivity) to water and h c is capillary suction, a positive quantity in unsaturated soils. As a result, the modification suggested to describe the high-suction range does not significantly affect the equivalence between Brooks-Corey (B-C) and van Genuchten [1980] parameters presented earlier. However, the shape of the retention curve near “natural saturation” has a significant impact on the value of the capillary drive. The estimate using the Brooks-Corey power law, extended to zero suction, will exceed that obtained with the new procedure by 25 to 30%. It is not possible to tell which procedure is appropriate. Tests on another data set, for which relative conductivity data are available, support the view of the authors that measurements of a retention curve coupled with a speculative curve of relative permeability as from a capillary model are not sufficient to accurately determine the (maximum) capillary drive. The capillary drive is a dynamic scalar, whereas the retention curve is of a static character. Only measurements of infiltration rates with time can determine the capillary drive with precision for a given soil.

Water Resources Research↗
USGS1999

Effects of hydraulic roughness on surface textures of gravel‐bed rivers

Field studies of forest gravel‐bed rivers in northwestern Washington and southeastern Alaska demonstrate that bed‐surface grain size is responsive to hydraulic roughness caused by bank irregularities, bars, and wood debris. We evaluate textural response by comparing reach‐average median grain size ( D 50 ) to that predicted from the total bank‐full boundary shear stress (т 0 bf ), representing a hypothetical reference condition of low hydraulic roughness. For a given т 0 bf , channels with progressively greater hydraulic roughness have systematically finer bed surfaces, presumably due to reduced bed shear stress, resulting in lower channel competence and diminished bed load transport capacity, both of which promote textural fining. In channels with significant hydraulic roughness, observed values D 50 can be up to 90% smaller than those predicted from т 0 bf . We find that wood debris plays an important role at our study sites, not only providing hydraulic roughness but also influencing pool spacing, frequency of textural patches, and the amplitude and wavelength of bank and bar topography and their consequent roughness. Our observations also have biological implications. We find that textural fining due to hydraulic roughness can create usable salmonid spawning gravels in channels that otherwise would be too coarse.

Water Resources Research↗
USGS1999

Evaluation of prediction intervals for expressing uncertainties in groundwater flow model predictions

We tested the accuracy of 95% individual prediction intervals for hydraulic heads, streamflow gains, and effective transmissivities computed by groundwater models of two Danish aquifers. To compute the intervals, we assumed that each predicted value can be written as the sum of a computed dependent variable and a random error. Testing was accomplished by using a cross‐validation method and by using new field measurements of hydraulic heads and transmissivities that were not used to develop or calibrate the models. The tested null hypotheses are that the coverage probability of the prediction intervals is not significantly smaller than the assumed probability (95%) and that each tail probability is not significantly different from the assumed probability (2.5%). In all cases tested, these hypotheses were accepted at the 5% level of significance. We therefore conclude that for the groundwater models of two real aquifers the individual prediction intervals appear to be accurate.

Water Resources Research↗
USGS1999

Pumping tests in networks of multilevel sampling wells: Motivation and methodology

The identification of spatial variations in hydraulic conductivity ( K ) on a scale of relevance for transport investigations has proven to be a considerable challenge. Recently, a new field method for the estimation of interwell variations in K has been proposed. This method, hydraulic tomography, essentially consists of a series of short‐term pumping tests performed in a tomographic‐like arrangement. In order to fully realize the potential of this approach, information about lateral and vertical variations in pumping‐induced head changes (drawdown) is required with detail that has previously been unobtainable in the field. Pumping tests performed in networks of multilevel sampling (MLS) wells can provide data of the needed density if drawdown can accurately and rapidly be measured in the small‐diameter tubing used in such wells. Field and laboratory experiments show that accurate transient drawdown data can be obtained in the small‐diameter MLS tubing either directly with miniature fiber‐optic pressure sensors or indirectly using air‐pressure transducers. As with data from many types of hydraulic tests, the quality of drawdown measurements from MLS tubing is quite dependent on the effectiveness of well development activities. Since MLS ports of the standard design are prone to clogging and are difficult to develop, alternate designs are necessary to ensure accurate drawdown measurements. Initial field experiments indicate that drawdown measurements obtained from pumping tests performed in MLS networks have considerable potential for providing valuable information about spatial variations in hydraulic conductivity.

Water Resources Research↗
USGS1999

Effects of sediment supply on surface textures of gravel‐bed rivers

Using previously published data from flume studies, we test a new approach for quantifying the effects of sediment supply (i.e., bed material supply) on surface grain size of equilibrium gravel channels. Textural response to sediment supply is evaluated relative to a theoretical prediction of competent median grain size ( D 50 ′ ). We find that surface median grain size ( D 50 ) varies inversely with sediment supply rate and systematically approaches the competent value ( D 50 ′ ) at low equilibrium transport rates. Furthermore, equilibrium transport rate is a power function of the difference between applied and critical shear stresses and is therefore a power function of the difference between competent and observed median grain sizes ( D 50 ′ and D 50 ). Consequently, we propose that the difference between predicted and observed median grain sizes can be used to determine sediment supply rate in equilibrium channels. Our analysis framework collapses data from different studies toward a single relationship between sediment supply rate and surface grain size. While the approach appears promising, we caution that it has been tested only on a limited set of laboratory data and a narrow range of channel conditions.

Water Resources Research↗
USGS1999

Coupling groundwater and riparian vegetation models to assess effects of reservoir releases

Although riparian areas in the arid southwestern United States are critical for maintaining species diversity, their extent and health have been declining since Euro‐American settlement. The purpose of this study was to develop a methodology to evaluate the potential for riparian vegetation restoration and groundwater recharge. A numerical groundwater flow model was coupled with a conceptual riparian vegetation model to predict hydrologic conditions favorable to maintaining riparian vegetation downstream of a reservoir. A Geographic Information System (GIS) was used for this one‐way coupling. Constant and seasonally varying releases from the dam were simulated using volumes anticipated to be permitted by a regional water supplier. Simulations indicated that seasonally variable releases would produce surface flow 5.4–8.5 km below the dam in a previously dry reach. Using depth to groundwater simulations from the numerical flow model with conceptual models of depths to water necessary for maintenance of riparian vegetation, the GIS analysis predicted a 5‐ to 6.5‐fold increase in the area capable of sustaining riparian vegetation.

Water Resources Research↗
USGS2000

Stream chemistry modeling of two watersheds in the Front Range, Colorado

We investigated the hydrologic, geochemical, and biogeochemical controls on stream chemical composition on the Green Lakes Valley and Andrews Creek watersheds using the alpine hydrochemical model (AHM). Both sites had comparable data sets from 1994 and 1996, including high‐resolution spatial data and high‐frequency time series of hydrology, geochemistry, and meteorology. The model of each watershed consisted of three terrestrial subunits (soil, talus, and rock), with the routing between the subunits determined by spatial land cover data. Using 1994 data for model calibration and 1996 data for evaluation, AHM captured the dominant processes and successfully simulated daily stream chemical composition on both watersheds. These results confirm our procedure of using spatial and site‐specific field and laboratory data to generate an initial catchment model and then calibrating the model to calculate effective parameters for unmeasured processes. A net source of nitrogen was identified in the Andrews Creek watershed during the spring snowmelt period, whereas nitrogen was immobilized in the Green Lakes Valley. This difference was most likely due to the larger and more dominant area of talus in the Andrews Creek watershed. Our results also indicate that routing of snowmelt through either soil or talus material is sufficient for retention of H + and release of base cations but that N retention is more important on areas mapped as soil. Owing to the larger ionic pulse and larger fraction of surface runoff the Green Lakes Valley was more sensitive to a doubling of wet deposition chemistry than the Andrews Creek watershed.

Colorado↗
USGS2000

Combining binary decision tree and geostatistical methods to estimate snow distribution in a mountain watershed

We model the spatial distribution of snow across a mountain basin using an approach that combines binary decision tree and geostatistical techniques. In April 1997 and 1998, intensive snow surveys were conducted in the 6.9‐km 2 Loch Vale watershed (LVWS), Rocky Mountain National Park, Colorado. Binary decision trees were used to model the large‐scale variations in snow depth, while the small‐scale variations were modeled through kriging interpolation methods. Binary decision trees related depth to the physically based independent variables of net solar radiation, elevation, slope, and vegetation cover type. These decision tree models explained 54–65% of the observed variance in the depth measurements. The tree‐based modeled depths were then subtracted from the measured depths, and the resulting residuals were spatially distributed across LVWS through kriging techniques. The kriged estimates of the residuals were added to the tree‐based modeled depths to produce a combined depth model. The combined depth estimates explained 60–85% of the variance in the measured depths. Snow densities were mapped across LVWS using regression analysis. Snow‐covered area was determined from high‐resolution aerial photographs. Combining the modeled depths and densities with a snow cover map produced estimates of the spatial distribution of snow water equivalence (SWE). This modeling approach offers improvement over previous methods of estimating SWE distribution in mountain basins.

Colorado↗
USGS2000

Bed material transport in the Virgin River, Utah

Detailed information concerning the rate and particle size distribution of bed material transport by streamflows can be very difficult and expensive to obtain, especially where peak streamflows are brief and bed material is poorly sorted, including some very large boulders. Such streams, however, are common in steep, arid watersheds. Any computational approach must consider that (1) only the smaller particle sizes present on the streambed move even during large floods and (2) the largest bed particles exert a significant form drag on the flow. Conventional methods that rely on a single particle size to estimate the skin friction shear stress acting on the mobile fraction of the bed material perform poorly. Instead, for this study, the skin friction shear stress was calculated for the observed range of streamflows by calculating the form drag exerted on the reach‐averaged flow field by all particle sizes. Suspended and bed load transported rates computed from reach‐averaged skin friction shear stress are in excellent agreement with measured transport rates. The computed mean annual bed material load, including both bed load and suspended load, of the East Fork Virgin River for the water years 1992‐1996 was approximately 1.3×10 5 t. A large portion of the bed material load consists of sand‐sized particles, 0.062–1.0 mm in diameter, that are transported in suspension. Such particles, however, constituted only 10% of the surface bed material and less than 25% of the subsurface bed material. The mean annual quantity of bed load transported was 1060 t/yr with a median size of 15 mm.

Utah↗
USGS1999

Reactive solute transport in streams: A surface complexation approach for trace metal sorption

A model for trace metals that considers in-stream transport, metal oxide precipitation-dissolution, and pH-dependent sorption is presented. Linkage between a surface complexation submodel and the stream transport equations provides a framework for modeling sorption onto static and/or dynamic surfaces. A static surface (e.g., an iron- oxide-coated streambed) is defined as a surface with a temporally constant solid concentration. Limited contact between solutes in the water column and the static surface is considered using a pseudokinetic approach. A dynamic surface (e.g., freshly precipitated metal oxides) has a temporally variable solid concentration and is in equilibrium with the water column. Transport and deposition of solute mass sorbed to the dynamic surface is represented in the stream transport equations that include precipitate settling. The model is applied to a pH-modification experiment in an acid mine drainage stream. Dissolved copper concentrations were depressed for a 3 hour period in response to the experimentally elevated pH. After passage of the pH front, copper was desorbed, and dissolved concentrations returned to ambient levels. Copper sorption is modeled by considering sorption to aged hydrous ferric oxide (HFO) on the streambed (static surface) and freshly precipitated HFO in the water column (dynamic surface). Comparison of parameter estimates with reported values suggests that naturally formed iron oxides may be more effective in removing trace metals than synthetic oxides used in laboratory studies. The model's ability to simulate pH, metal oxide precipitation-dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between trace metal chemistry and hydrologic transport at the field scale.

Water Resources Research↗
USGS2000

Sensitivity of a high-elevation Rocky Mountain watershed to altered climate and CO2

We explored the hydrologic and ecological responses of a headwater mountain catchment, Loch Vale watershed, to climate change and doubling of atmospheric CO 2 scenarios using the Regional Hydro-Ecological Simulation System (RHESSys). A slight (2°C) cooling, comparable to conditions observed over the past 40 years, led to greater snowpack and slightly less runoff, evaporation, transpiration, and plant productivity. An increase of 2°C yielded the opposite response, but model output for an increase of 4°C showed dramatic changes in timing of hydrologic responses. The snowpack was reduced by 50%, and runoff and soil water increased and occurred 4–5 weeks earlier with 4°C warming. Alpine tundra photosynthetic rates responded more to warmer and wetter conditions than subalpine forest, but subalpine forest showed a greater response to doubling of atmospheric CO 2 than tundra. Even though water use efficiency increased with the double CO 2 scenario, this had little effect on basin-wide runoff because the catchment is largely unvegetated. Changes in winter and spring climate conditions were more important to hydrologic and vegetation dynamics than changes that occurred during summer.

Water Resources Research↗
USGS2000

Estimating the theoretical semivariogram from finite numbers of measurements

We investigate from a theoretical basis the impacts of the number, location, and correlation among measurement points on the quality of an estimate of the semivariogram. The unbiased nature of the semivariogram estimator ŷ ( r ) is first established for a general random process Z( x ). The variance of ŷ Z ( r ) is then derived as a function of the sampling parameters (the number of measurements and their locations). In applying this function to the case of estimating the semivariograms of the transmissivity and the hydraulic head field, it is shown that the estimation error depends on the number of the data pairs, the correlation among the data pairs (which, in turn, are determined by the form of the underlying semivariogram γ( r )), the relative locations of the data pairs, and the separation distance at which the semivariogram is to be estimated. Thus design of an optimal sampling program for semivariogram estimation should include consideration of each of these factors. Further, the function derived for the variance of ŷ Z ( r ) is useful in determining the reliability of a semivariogram developed from a previously established sampling design.

Water Resources Research↗