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 325 records · Page 18Linked to original sources
USGS2000

Effects of nearshore recharge on groundwater interactions with a lake in mantled karst terrain

The recharge and discharge of groundwater were investigated for a lake basin in the mantled karst terrain of central Florida to determine the relative importance of transient groundwater inflow to the lake water budget. Variably saturated groundwater flow modeling simulated water table responses observed beneath two hillsides radiating outward from the groundwater flow‐through lake. Modeling results indicated that transient water table mounding and groundwater flow reversals in the nearshore region following large daily rainfall events generated most of the net groundwater inflow to the lake. Simulated daily groundwater inflow was greatest following water table mounding near the lake, not following subsequent peaks in the water level of upper basin wells. Transient mounding generated net groundwater inflow to the lake, that is, groundwater inflow in excess of the outflow occurring through the deeper lake bottom. The timing of the modeled net groundwater inflow agreed with an independent lake water budget; however, the quantity was considerably less than the budget‐derived value.

Water Resources Research↗
USGS2000

A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development

This paper, the first of two parts [see Abrams and Loague , this issue], takes the compartmentalized approach for the geochemical evolution of redox zones presented by Abrams et al . [1998] and embeds it within a solute transport framework. In this paper the compartmentalized approach is generalized to facilitate the description of its incorporation into a solute transport simulator. An equivalent formulation is developed which removes any discontinuities that may occur when switching compartments. Rate‐limited redox reactions are modeled with a modified Monod relationship that allows either the organic substrate or the electron acceptor to be the rate‐limiting reactant. Thermodynamic constraints are used to inhibit lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium on the lower‐energy reactions. The procedure used allows any redox reaction to be simulated as being kinetically limited or thermodynamically limited, depending on local geochemical conditions. Empirical reaction inhibition methods are not needed. The sequential iteration approach (SIA), a technique which allows the number of solute transport equations to be reduced, is adopted to solve the coupled geochemical/solute transport problem. When the compartmentalized approach is embedded within the SIA, with the total analytical concentration of each component as the dependent variable in the transport equation, it is possible to reduce the number of transport equations even further than with the unmodified SIA. A one‐dimensional, coupled geochemical/solute transport simulation is presented in which redox zones evolve dynamically in time and space. The compartmentalized solute transport (COMPTRAN) model described in this paper enables the development of redox zones to be simulated under both kinetic and thermodynamic constraints. The modular design of COMPTRAN facilitates the use of many different, preexisting solute transport and geochemical codes. The companion paper [ Abrams and Loague , this issue] presents examples of the application of COMPTRAN to field‐scale problems.

Water Resources Research↗
USGS2010

Shallow subsurface storm flow in a forested headwater catchment: Observations and modeling using a modified TOPMODEL

Transient, perched water tables in the shallow subsurface are observed at the South Fork Brokenback Run catchment in Shenandoah National Park, Virginia. Crest piezometers installed along a hillslope transect show that the development of saturated conditions in the upper 1.5 m of the subsurface is controlled by total precipitation and antecedent conditions, not precipitation intensity, although soil heterogeneities strongly influence local response. The macroporous subsurface storm flow zone provides a hydrological pathway for rapid runoff generation apart from the underlying groundwater zone, a conceptualization supported by the two‐storage system exhibited by hydrograph recession analysis. A modified version of TOPMODEL is used to simulate the observed catchment dynamics. In this model, generalized topographic index theory is applied to the subsurface storm flow zone to account for logarithmic storm flow recessions, indicative of linearly decreasing transmissivity with depth. Vertical drainage to the groundwater zone is required, and both subsurface reservoirs are considered to contribute to surface saturation.

Water Resources Research↗
USGS2000

Dating young groundwater with sulfur hexafluoride: Natural and anthropogenic sources of sulfur hexafluoride

Sulfur hexafluoride (SF 6 ) is primarily of anthropogenic origin but also occurs naturally. The troposphere concentration of SF 6 has increased from a steady state value of 0.054±0.009 to more than 4 parts per trillion volume during the past 40 years. An analytical procedure was developed for measuring concentrations of SF 6 to less than 0.01 fmol/L in water. Groundwater can be dated with SF 6 if it is in equilibrium with atmospheric SF 6 at the time of recharge and does not contain significant SF 6 from other sources. The dating range of SF 6 is currently 0 to 30 years. The tracer was successfully used to date shallow groundwater of the Atlantic Coastal Plain sand aquifers of the United States and springs issuing near the top of the Blue Ridge Mountains of Virginia. Significant concentrations of naturally occurring SF 6 were found in some igneous, volcanic, and sedimentary rocks and in some hydrothermal fluids.

Water Resources Research↗
USGS2010

An evaluation of the Wyoming Gauge System for snowfall measurement

The Wyoming snow fence (shield) has been widely used with precipitation gauges for snowfall measurement at more than 25 locations in Alaska since the late 1970s. This gauge's measurements have been taken as the reference for correcting wind‐induced gauge undercatch of snowfall in Alaska. Recently, this fence (shield) was tested in the World Meteorological Organization Solid Precipitation Measurement Intercomparison Project at four locations in the United States of America and Canada for six winter seasons. At the Intercomparison sites an octagonal vertical Double Fence with a Russian Tretyakov gauge or a Universal Belfort recording gauge was installed and used as the Intercomparison Reference (DFIR) to provide true snowfall amounts for this Intercomparison experiment. The Intercomparison data collected were compiled at the four sites that represent a variety of climate, terrain, and exposure. On the basis of these data sets the performance of the Wyoming gauge system for snowfall observations was carefully evaluated against the DFIR and snow cover data. The results show that (1) the mean snow catch efficiency of the Wyoming gauge compared with the DFIR is about 80–90%, (2) there exists a close linear relation between the measurements of the two gauge systems and this relation may serve as a transfer function to adjust the Wyoming gauge records to obtain an estimate of the true snowfall amount, (3) catch efficiency of the Wyoming gauge does not change with wind speed and temperature, and (4) Wyoming gauge measurements are generally compatible to the snowpack water equivalent at selected locations in northern Alaska. These results are important to our effort of determining true snowfall amounts in the high latitudes, and they are also useful for regional hydrologic and climatic analyses.

Water Resources Research↗
USGS2010

Comment on RamaRao et al. [1995] and LaVenue et al. [1995]

A method for stochastic modeling of groundwater flow systems using a combination of pilot point parameterization and conditional simulation was presented by RamaRao et. al [1995] and LaVenue et. al [1995]. (We will collectively term these two papers RLMM and term the method developed in RLMM the CS method here.) RLMM (pp. 478-479) state that the CS method is intended to provide a frequency distribution of possible alternative spatial transmissivity T distributions that (1) are statistically similar to the observed T distribution, (2) are equally likely given the calibration data, and (3) closely reproduce the measured pressures. The frequency distribution of transmissivities is then used to form frequency distribution of derived functions, such as travel times, which are summarized in the form of uncertainty measures, such as (RLMM, p. 512) "confidence (or tolerance) intervals," on the actual values, which in the case of travels times, are values that could occur sometime in the future (RLMM, p.513). Cooley [2000] analyzes the RLMM method using linearization and bootstrap theory and concludes that that the method can yield accurate uncertainty estimates but only under some limited circumstances. In this comment we use Cooley's analysis to critique the method. We also identify and discuss some statements made my by RLMM about model calibration and their methodology that appear to be misleading. It is unusual to comment on a paper so long after publication. Subsequent work has expanded the method of RLMM, and the method was used advantageously in the testing documented by Zimmerman et. al [1998], so it is clear the method has significant strengths. We go back to the 1995 papers for this comment, however, because they display most clearly the methodological difficulties with which we are concerned.

Water Resources Research↗
USGS2000

Testing a full‐range soil‐water retention function in modeling water potential and temperature

Recent work has emphasized development of full‐range water‐retention functions that are applicable under both wet and dry soil conditions, but evaluation of such functions in numerical modeling has been limited. Here we show that simulations using the Rossi‐Nimmo (RN) full‐range function compared favorably with those using the common Brooks‐Corey function and that the RN function can improve prediction of water potentials in near‐surface soil, particularly under dry conditions. Simulations using the RN function also improved prediction of temperatures throughout the soil profile. Such improvements could be important for calculations of liquid and vapor flow in near‐surface soils and in deep unsaturated zones of arid and semiarid regions.

Water Resources Research↗
USGS2000

Unsaturated‐zone wedge beneath a large, natural lake

Lake Belle Taine (480 ha) in north central Minnesota receives on average 76,000 m 3 d −1 of water from Little Sand Creek but has no outlet. Water seeps out of the lake, flows through glacial outwash, and discharges into three nearby lakes with stages 13–14 m lower than Belle Taine. Seepage‐meter data indicate water is seeping out of Belle Taine at velocities up to 263 cm d −1 . Hydraulic‐head measurements made in the lake bed indicate the sediments are unsaturated beneath a portion of the lake edge, and a wedge of unsaturated sediments extends beneath the lake bed as much as 20 m from the shoreline. At the shoreline the water table is as much as 6.7 m below the lake surface. Modeling results of a similar hypothetical setting indicate that the horizontal extent of an unsaturated zone beneath a lake depends on (1) the permeability contrast between the outwash and the lake bed, (2) anisotropy, (3) lake bed slope, and (4) thickness of the lower‐permeability lake bed sediments. Although rarely documented, unsaturated sediments beneath a lake may not be extremely uncommon. Similar, much smaller unsaturated‐zone areas also have been observed beneath two other lakes in Minnesota.

Minnesota↗
USGS2000

Variability of bed mobility in natural, gravel‐bed channels and adjustments to sediment load at local and reach scales

Local variations in boundary shear stress acting on bed‐surface particles control patterns of bed load transport and channel evolution during varying stream discharges. At the reach scale a channel adjusts to imposed water and sediment supply through mutual interactions among channel form, local grain size, and local flow dynamics that govern bed mobility. In order to explore these adjustments, we used a numerical flow model to examine relations between model‐predicted local boundary shear stress (т j ( and measured surface particle size ( D 50 ) at bank‐full discharge in six gravel‐bed, alternate‐bar channels with widely differing annual sediment yields. Values of т j and D 50 were poorly correlated such that small areas conveyed large proportions of the total bed load, especially in sediment‐poor channels with low mobility. Sediment‐rich channels had greater areas of full mobility; sediment‐poor channels had greater areas of partial mobility; and both types had significant areas that were essentially immobile. Two reach‐mean mobility parameters (Shields stress and Q *) correlated reasonably well with sediment supply. Values which can be practicably obtained from carefully measured mean hydraulic variables and particle size would provide first‐order assessments of bed mobility that would broadly distinguish the channels in this study according to their sediment yield and bed mobility.

Water Resources Research↗
USGS2000

Modeling the influence of variable pH on the transport of zinc in a contaminated aquifer using semiempirical surface complexation models

Land disposal of sewage effluent resulted in contamination of a sand and gravel aquifer (Cape Cod, Massachusetts) with zinc (Zn). The distribution of Zn was controlled by pH‐dependent adsorption; the Zn extended 15 m into the 30‐m‐thick sewage plume within approximately 100 m of the source but only 2–4 m into the plume between 100 and 400 m downgradient. A two‐dimensional vertical cross section model coupling groundwater flow with solute transport and equilibrium adsorption is used to simulate the influence of pH on Zn transport. Adsorption is described using semiempirical surface complexation models (SCM) by writing chemical reactions between dissolved Zn and mineral surface sites. SCM parameters were determined in independent laboratory experiments. A 59‐year simulation with a one‐site SCM describes the influence of pH on Zn transport well, with greater mobility at the low pH values near the upper sewage plume boundary than at the higher pH values deeper in the sewage‐contaminated zone. Simulation with a two‐site SCM describes both the sharpness and approximate location of the leading edge of the Zn‐contaminated region. Temporal variations in pH of incoming groundwater can result in large increases in Zn concentration and mobility. The influence of spatial and temporal variability in pH on adsorption and transport of Zn was accomplished much more easily with the semiempirical SCM approach than could be achieved with distribution coefficients or adsorption isotherms.

Water Resources Research↗
USGS2001

Quantifying the relative importance of flow regulation and grain size regulation of suspended sediment transport α and tracking changes in grain size of bed sediment β

To predict changes in sediment transport, it is essential to know whether transport is regulated mainly by changes in flow or by changes in grain size of sediment on the bed. In flows where changes in suspended sediment transport are regulated purely by changes in flow (grain size of bed sediment is constant), increases in flow strength cause increases in both concentration and grain size of sediment in suspension (because stronger flows are able to suspend more sediment and coarser grains). Under this constraint of constant grain size of bed sediment concentration and median diameter of suspended sediment are positively correlated. In contrast, where transport is regulated purely by changes in grain size of sediment on the bed, concentration and median diameter of suspended sediment are negatively correlated (because increasing the median diameter of the bed sediment causes the concentration to decrease while causing the median grain size in suspension to increase). Where both flow strength and grain size on the bed are free to vary, the relation between concentration and grain size in suspension can be used to quantify the importance of grain size regulation relative to flow regulation of sediment transport, a measure defined as α. To predict sediment transport in systems that are regulated dominantly by changes in grain size on the bed, it is more useful to measure sediment input events or changes in grain size on the bed than to measure changes in flow. More commonly, grain size of bed sediment may be secondary to flow in regulating transport but may, nevertheless, be important. The relative coarseness of bed sediment (β) can be measured directly or, like α, can be calculated from measurements of concentration and grain size of suspended sediment.

Water Resources Research↗
USGS2001

Analysis of streambed temperatures in ephemeral channels to determine streamflow frequency and duration

Spatial and temporal patterns in streamflow are rarely monitored for ephemeral streams. Flashy, erosive streamflows common in ephemeral channels create a series of operational and maintenance problems, which makes it impractical to deploy a series of gaging stations along ephemeral channels. Streambed temperature is a robust and inexpensive parameter to monitor remotely, leading to the possibility of analyzing temperature patterns to estimate streamflow frequency and duration along ephemeral channels. A simulation model was utilized to examine various atmospheric and hydrological upper boundary conditions compared with a series of hypothetical temperature‐monitoring depths within the streambed. Simulation results indicate that streamflow events were distinguished from changing atmospheric conditions with greater certainty using temperatures at shallow depths (e.g., 10–20 cm) as opposed to the streambed surface. Three ephemeral streams in the American Southwest were instrumented to monitor streambed temperature for determining the accuracy of using this approach to ascertain the long‐term temporal and spatial extent of streamflow along each stream channel. Streambed temperature data were collected at the surface or at shallow depth along each stream channel, using thermistors encased in waterproof, single‐channel data loggers tethered to anchors in the channel. On the basis of comparisons with site information, such as direct field observations and upstream flow records, diurnal temperature variations successfully detected the presence and duration of streamflow for all sites.

Water Resources Research↗
USGS2001

A new tracer‐density criterion for heterogeneous porous media

Tracer experiments provide information about aquifer material properties vital for accurate site characterization. Unfortunately, density‐induced sinking can distort tracer movement, leading to an inaccurate assessment of material properties. Yet existing criteria for selecting appropriate tracer concentrations are based on analysis of homogeneous media instead of media with heterogeneities typical of field sites. This work introduces a hydraulic‐gradient correction for heterogeneous media and applies it to a criterion previously used to indicate density‐induced instabilities in homogeneous media. The modified criterion was tested using a series of two‐dimensional heterogeneous intermediate‐scale tracer experiments and data from several detailed field tracer tests. The intermediate‐scale experimental facility (10.0×1.2×0.06 m) included both homogeneous and heterogeneous (σ ln k 2 = 1.22) zones. The field tracer tests were less heterogeneous (0.24 < σ ln k 2 < 0.37), but measurements were sufficient to detect density‐induced sinking. Evaluation of the modified criterion using the experiments and field tests demonstrates that the new criterion appears to account for the change in density‐induced sinking due to heterogeneity. The criterion demonstrates the importance of accounting for heterogeneity to predict density‐induced sinking and differences in the onset of density‐induced sinking in two‐ and three‐dimensional systems.

Water Resources Research↗
USGS2001

Methyl tert‐butyl ether degradation in the unsaturated zone and the relation between MTBE in the atmosphere and shallow groundwater

Atmospheric methyl tert‐butyl ether (MTBE) concentrations in southern New Jersey generally exceeded concentrations in samples taken from the unsaturated zone. A simple unsaturated zone transport model indicates that MTBE degradation can explain the attenuation with half‐lives from a few months to a couple of years. Tert‐butyl alcohol (TBA), a possible degradation product of MTBE, was detected in unsaturated‐zone samples at concentrations exceeding atmospheric levels at some sites, suggesting the possible conversion of MTBE to TBA. At sites where MTBE was detected in shallow groundwater, the concentration was typically higher than the overlying unsaturated‐zone concentration. This observation is consistent with outgassing from the aquifer and combined with the unsaturated‐zone attenuation suggests some of the MTBE detections in shallow groundwater are nonatmospheric in origin, coming from leaking tanks, road runoff, or other sources. The identification of sources of MTBE in groundwater and attenuation mechanisms through the hydrologic cycle is critical in developing an understanding of the long‐term effect of MTBE releases.

Water Resources Research↗
USGS2001

Effective matrix diffusion in kilometer‐scale transport in fractured crystalline rock

Concentrations of tritium ( 3 H) and dichlorodifluoromethane (CFC‐12) in water samples taken from glacial drift and fractured crystalline rock over 4 km 2 in central New Hampshire are interpreted to identify a conceptual model of matrix diffusion and the magnitude of the diffusion coefficient. Dispersion and mass transfer to and from fractures has affected the 3 H concentration to the extent that the peak 3 H concentration of the 1960s is no longer distinguishable. Because of heterogeneity in the bedrock the sparsely distributed chemical data do not warrant a three‐dimensional transport model. Instead, a one‐dimensional model of CFC‐12 and 3 H migration along flow lines in the glacial drift and bedrock is used to place bounds on the processes affecting kilometer‐scale transport, arid model parameters are varied to reproduce the measured relation between 3 H and CFC‐12, rather than their spatial distributions. A model of mass exchange to and from fractures that is dependent on the time‐varying concentration gradient at fracture surfaces qualitatively reproduces the measured relation between 3 H and CFC‐12 with an upper bound for the fracture dispersivity approximately equal to 250 m and a lower bound for the effective matrix diffusion coefficient equal to 1 m 2 yr −1 . The diffusion coefficient at the kilometer scale is at least 3 orders of magnitude greater than laboratory estimates of diffusion in crystalline rock. The large diffusion coefficient indicates that diffusion into an immobile fluid phase (rock matrix) is masked at the kilometer scale by advective mass exchange between fractures with large contrasts in trarismissivity. The measured transmissivity of fractures in the study area varies over more than 6 orders of magnitude. Advective mass exchange from high‐permeability fractures to low‐permeability fractures results in short migration distances of a chemical constituent in low‐permeability fractures over an extended period of time before reentering high‐permeability fractures; viewed at the kilometer scale, this process is analogous to the chemical constituent diffusing into and out of an immobile fluid phase.

Water Resources Research↗
USGS2001

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↗
USGS2001

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↗
USGS2001

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↗