USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21Linked to original sources
USGS2003

Transport and cycling of iron and hydrogen peroxide in a freshwater stream: Influence of organic acids

An in‐stream injection of two dissolved organic acids (phthalic and aspartic acids) was performed in an acidic mountain stream to assess the effects of organic acids on Fe photoreduction and H 2 O 2 cycling. Results indicate that the fate of Fe is dependent on a net balance of oxidative and reductive processes, which can vary over a distance of several meters due to changes in incident light and other factors. Solution phase photoreduction rates were high in sunlit reaches and were enhanced by the organic acid addition but were also limited by the amount of ferric iron present in the water column. Fe oxide photoreduction from the streambed and colloids within the water column resulted in an increase in the diurnal load of total filterable Fe within the experimental reach, which also responded to increases in light and organic acids. Our results also suggest that Fe(II) oxidation increased in response to the organic acids, with the result of offsetting the increase in Fe(II) from photoreductive processes. Fe(II) was rapidly oxidized to Fe(III) after sunset and during the day within a well‐shaded reach, presumably through microbial oxidation. H 2 O 2 , a product of dissolved organic matter photolysis, increased downstream to maximum concentrations of 0.25 μM midday. Kinetic calculations show that the buildup of H 2 O 2 is controlled by reaction with Fe(II), but this has only a small effect on Fe(II) because of the small formation rates of H 2 O 2 compared to those of Fe(II). The results demonstrate the importance of incorporating the effects of light and dissolved organic carbon into Fe reactive transport models to further our understanding of the fate of Fe in streams and lakes.

Water Resources Research↗
USGS2003

Comparisons of two moments‐based estimators that utilize historical and paleoflood data for the log Pearson type III distribution

The expected moments algorithm (EMA) [ Cohn et al. , 1997 ] and the Bulletin 17B [ Interagency Committee on Water Data , 1982 ] historical weighting procedure (B17H) for the log Pearson type III distribution are compared by Monte Carlo computer simulation for cases in which historical and/or paleoflood data are available. The relative performance of the estimators was explored for three cases: fixed‐threshold exceedances, a fixed number of large floods, and floods generated from a different parent distribution. EMA can effectively incorporate four types of historical and paleoflood data: floods where the discharge is explicitly known, unknown discharges below a single threshold, floods with unknown discharge that exceed some level, and floods with discharges described in a range. The B17H estimator can utilize only the first two types of historical information. Including historical/paleoflood data in the simulation experiments significantly improved the quantile estimates in terms of mean square error and bias relative to using gage data alone. EMA performed significantly better than B17H in nearly all cases considered. B17H performed as well as EMA for estimating X 100 in some limited fixed‐threshold exceedance cases. EMA performed comparatively much better in other fixed‐threshold situations, for the single large flood case, and in cases when estimating extreme floods equal to or greater than X 500 . B17H did not fully utilize historical information when the historical period exceeded 200 years. Robustness studies using GEV‐simulated data confirmed that EMA performed better than B17H. Overall, EMA is preferred to B17H when historical and paleoflood data are available for flood frequency analysis.

Water Resources Research↗
USGS2003

Radiogenic helium in shallow groundwater within a clay till, southwestern Ontario

Profiles of 4 He in pore water were measured in clay aquitards in SW Ontario. The 4 He distributions are consistent with groundwater velocities that are <6 mm yr −1 , and thus diffusion is the dominant transport mechanism for 4 He. Modeling indicates that the effective diffusion coefficient for 4 He is 6.3 ± 1.6 × 10 −6 cm 2 s −1 . Furthermore, the profiles are consistent with the internal release of 4 He from aquitard sediments at a rate of 0.03–0.13 μcc(STP) kg −1 yr −1 . These rates are also consistent with laboratory release experiments and, on average, are 600 times greater than the production of 4 He from U/Th decay. Modeling and the ratio of 21 Ne to 4 He within the sediments indicate that although the aquitards were deposited about 13 kA BP, the sediments released >70% of initial 4 He for 50 to 60 kA prior to incorporation into the till.

Water Resources Research↗
USGS2003

Discharge indices for water quality loads

Effective discharge has been used to describe the streamflow level that is responsible for transporting the most sediment over the long term. Careful inspection reveals that this concept may not have been well defined, and different interpretations have led to conflicting representations. Because total load is ultimately the quantity of interest, we define a new index, the half‐load discharge, which is that discharge above and below which half the total long‐term load is transported. The value of the half‐load discharge is derived for a reasonable model of flows and constituent concentration. The effective discharge has generally been thought to be a relatively common or frequent flood. The half‐load discharge is generally a much greater and less frequent flow than commonly used estimators of the effective discharge. Relations provided here for the frequency and magnitude of the half‐load discharge provide evidence that it is relatively rare floods that transport most of the sediment over the long term. These ideas apply to other constituents as well.

Water Resources Research↗
USGS2004

Investigating the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, using a three‐dimensional inverse flow and transport model

Flowmeter‐measured hydraulic conductivities from the heterogeneous MADE site have been used predictively in advection‐dispersion models. Resulting simulated concentrations failed to reproduce even major plume characteristics and some have concluded that other mechanisms, such as dual porosity, are important. Here an alternative possibility is investigated: that the small‐scale flowmeter measurements are too noisy and possibly too biased to use so directly in site‐scale models and that the hydraulic head and transport data are more suitable for site‐scale characterization. Using a calibrated finite element model of the site and a new framework to evaluate random and systematic model and measurement errors, the following conclusions are derived. (1) If variations in subsurface fluid velocities like those simulated in this work (0.1 and 2.0 m per day along parallel and reasonably close flow paths) exist, it is likely that classical advection‐dispersion processes can explain the measured plume characteristics. (2) The flowmeter measurements are possibly systematically lower than site‐scale values when the measurements are considered individually and using common averaging methods and display variability that obscures abrupt changes in hydraulic conductivities that are well supported by changes in hydraulic gradients and are important to the simulation of transport.

Water Resources Research↗
USGS2004

Comment on “Piezometric response in shallow bedrock at CB1: Implications for runoff generation and landsliding” by David R. Montgomery, William E. Dietrich, and John T. Heffner

Piezometric responses to rainfall on hillslopes commonly dictate the timing of landsliding. Insight to this phenomenon can be gained by evaluating the timescales for pore pressure perturbations to propagate normal and parallel to the ground surface, and these timescales can be estimated using characteristic values of hydraulic diffusivity [ Iverson , 2000 ]. However, Montgomery et al. [2002] employed an erroneous definition of hydraulic diffusivity, leading to flawed assessment of the piezometric response timescales identified by Iverson [2000] . This comment aims to rectify the errors of Montgomery et al. [2002] and clarify the methods used by Iverson [2000] to estimate piezometric response timescales at the CB1 field site studied by Montgomery et al.

Water Resources Research↗
USGS2004

Temporal scaling of hydraulic head and river base flow and its implication for groundwater recharge

Spectral analyses were conducted for hourly hydraulic head ( h ) data observed over a 4‐year period at seven monitoring wells in the Walnut Creek watershed, Iowa. The log power spectral density of the hydraulic head fluctuations versus log frequency ( f ) at all seven wells is shown to have a distinct slope or fractal dimension ( D ), indicating temporal scaling in the time series of water level fluctuations. The fractal dimension of the time series varies from well to well, and the spectrum for the average h over all seven wells has a fractal dimension of 1.46 and Hurst coefficient of 0.54. The log power spectral density of estimated base flow in the Walnut Creek and four other watersheds versus log f is shown to have two distinct slopes with a break in scaling at about 30 days. It is shown that the groundwater recharge process in a basin can be estimated from a head spectrum based on existing theoretical results. Hydraulic head in an aquifer may fluctuate as a fractal in time in response to either a white noise or fractal recharge process, depending on physical parameters (i.e., transmissivity and specific yield) of the aquifer. The recharge process at the Walnut Creek watershed is shown to have a white noise spectrum based on the observed head spectrum.

Water Resources Research↗
USGS2004

Conservative and reactive solute transport in constructed wetlands

The transport of bromide, a conservative tracer, and rhodamine WT (RWT), a photodegrading tracer, was evaluated in three wastewater‐dependent wetlands near Phoenix, Arizona, using a solute transport model with transient storage. Coupled sodium bromide and RWT tracer tests were performed to establish conservative transport and reactive parameters in constructed wetlands with water losses ranging from (1) relatively impermeable (15%), (2) moderately leaky (45%), and (3) significantly leaky (76%). RWT first‐order photolysis rates and sorption coefficients were determined from independent field and laboratory experiments. Individual wetland hydraulic profiles influenced the extent of transient storage interaction in stagnant water areas and consequently RWT removal. Solute mixing and transient storage interaction occurred in the impermeable wetland, resulting in 21% RWT mass loss from main channel and storage zone photolysis (10%) and sorption (11%) reactions. Advection and dispersion governed solute transport in the leaky wetland, limiting RWT photolysis removal (1.2%) and favoring main channel sorption (3.6%). The moderately leaky wetland contained islands parallel to flow, producing channel flow and minimizing RWT losses (1.6%).

Water Resources Research↗
USGS2004

Spatial distribution of the largest rainfall‐runoff floods from basins between 2.6 and 26,000 km2 in the United States and Puerto Rico

We assess the spatial distribution of the largest rainfall‐generated streamflows from a database of 35,663 flow records composed of the largest 10% of annual peak flows from each of 14,815 U.S. Geological Survey stream gaging stations in the United States and Puerto Rico. High unit discharges (peak discharge per unit contributing area) from basins with areas of 2.6 to 26,000 km 2 (1–10,000 mi 2 ) are widespread, but streams in Hawaii, Puerto Rico, and Texas together account for more than 50% of the highest unit discharges. The Appalachians and western flanks of Pacific coastal mountain systems are also regions of high unit discharges, as are several areas in the southern Midwest. By contrast, few exceptional discharges have been recorded in the interior West, northern Midwest, and Atlantic Coastal Plain. Most areas of high unit discharges result from the combination of (1) regional atmospheric conditions that produce large precipitation volumes and (2) steep topography, which enhances precipitation by convective and orographic processes and allows flow to be quickly concentrated into stream channels. Within the conterminous United States, the greatest concentration of exceptional unit discharges is at the Balcones Escarpment of central Texas, where maximum U.S. rainfall amounts apparently coincide with appropriate basin physiography to produce many of the largest measured U.S. floods. Flood‐related fatalities broadly correspond to the spatial distribution of high unit discharges, with Texas having nearly twice the average annual flood‐related fatalities of any other state.

Water Resources Research↗
USGS2004

Lateral mixing in the Mississippi River below the confluence with the Ohio River

Lateral dispersion coefficients for two dispersants were determined for three sections of the Mississippi River below the confluence with the Ohio River. The dispersants were the specific conductance and an industrial organic compound (trimethyltriazinetrione). Three models based on the stream tube concept were used, and lateral dispersion coefficients computed from these models were comparable. Coefficients for the two dispersants also were comparable. Lateral dispersion coefficients were consistent with expectations based on the characteristics of the river sections. Overall average values were 0.444 m 2 /s for a relatively straight section of river, 1.69 m 2 /s for a section containing two sharp bends, and 2.22 m 2 /s for a long section containing four sharp bends and several small islands. The lateral dispersion coefficients measured for the Mississippi River are consistent with literature data and a water discharge relation. Results of this study provide lateral dispersion coefficients for a water discharge not previously reported in the literature as well as new values for the Mississippi River.

Ohio River, Mississippi River↗
USGS2004

Evaluation of longitudinal dispersivity estimates from simulated forced‐ and natural‐gradient tracer tests in heterogeneous aquifers

We simulate three types of forced‐gradient tracer tests (converging radial flow, unequal strength two well, and equal strength two well) and natural‐gradient tracer tests in multiple realizations of heterogeneous two‐dimensional aquifers with a hydraulic conductivity distribution characterized by a spherical variogram. We determine longitudinal dispersivities (α L ) by analysis of forced‐gradient test breakthrough curves at the pumped well and by spatial moment analysis of tracer concentrations during the natural‐gradient tests. Results show that among the forced‐gradient tests, a converging radial‐flow test tends to yield the smallest α L , an equal strength two‐well test tends to yield the largest α L , and an unequal strength two‐well test tends to yield an intermediate value. This finding is qualitatively explained by considering the aquifer area sampled by a particular test. A converging radial‐flow test samples a small area, and thus the tracer undergoes a low degree of spreading and mixing. An equal strength two‐well test samples a much larger area, so the tracer is spread and mixed to a greater degree. Results also suggest that if the distance between the tracer source well and the pumped well is short relative to the lengths over which velocity is correlated, then the α L estimate can be highly dependent on local heterogeneities in the vicinity of the wells. Finally, results indicate that α L estimated from forced‐gradient tracer tests can significantly underestimate the α L needed to characterize solute dispersion under natural‐gradient flow. Only a two‐well tracer test with a large well separation in an aquifer with a low degree of heterogeneity can yield a value of α L that characterizes natural‐gradient tracer spreading. This suggests that a two‐well test with a large well separation is the preferred forced‐gradient test for characterizing solute dispersion under natural‐gradient flow.

Water Resources Research↗
USGS2005

Changes in the water surface profile of the Colorado River in Grand Canyon, Arizona, between 1923 and 2000

In 1923, a U.S. Geological Survey expedition surveyed the water surface profile of the Colorado River through Grand Canyon with theodolite and stadia rod. In 2000, lidar overflights collected topographic data centered on the river corridor, allowing construction of a new water surface profile and detection of change in the profile since 1923. By registering the surveys with respect to each other on the basis of 11 locations that were independently determined to have been unchanged between 1923 and 2000, 80 rapids were directly compared for change between 1923 and 2000. The average change for all measured rapids was +0.26 m, indicating net aggradation of the coarse‐grained alluvium forming the rapids throughout Grand Canyon. In addition, comparison of the two water surface profiles showed enhanced pool‐and‐rapid morphology. While 50% of the total drop of the river occurred in just 9% of the river distance in 1923, that value increased to 66% by 2000.

Arizona↗
USGS2004

Unsaturated flow and transport through a fault embedded in fractured welded tuff

To evaluate the importance of matrix diffusion as a mechanism for retarding radionuclide transport in the vicinity of a fault located in unsaturated fractured rock, we carried out an in situ field experiment in the Exploratory Studies Facility at Yucca Mountain, Nevada. This experiment involved the release of ∼82,000 L of water over a period of 17 months directly into a near‐vertical fault under both constant positive head (at ∼0.04 m) and decreasing fluxes. A mix of conservative tracers (pentafluorobenzoic acid (PFBA) and bromide (applied in the form of lithium bromide)) was released along the fault over a period of 9 days, 7 months after the start of water release along the fault. As water was released into the fault, seepage rates were monitored in a large cavity excavated below the test bed. After the release of tracers, seepage water was continuously collected from three locations and analyzed for the injected tracers. Observations of bromide concentrations in seepage water during the early stages of the experiment and bromide and PFBA concentrations in the seepage water indicate the significant effects of matrix diffusion on transport through a fault embedded in fractured, welded rock.

Water Resources Research↗
USGS2004

Inferring time‐varying recharge from inverse analysis of long‐term water levels

Water levels in aquifers typically vary in response to time‐varying rates of recharge, suggesting the possibility of inferring time‐varying recharge rates on the basis of long‐term water level records. Presumably, in the southwestern United States (Arizona, Nevada, New Mexico, southern California, and southern Utah), rates of mountain front recharge to alluvial aquifers depend on variations in precipitation rates due to known climate cycles such as the El Niño‐Southern Oscillation index and the Pacific Decadal Oscillation. This investigation examined the inverse application of a one‐dimensional analytical model for periodic flow described by Lloyd R. Townley in 1995 to estimate periodic recharge variations on the basis of variations in long‐term water level records using southwest aquifers as the case study. Time‐varying water level records at various locations along the flow line were obtained by simulation of forward models of synthetic basins with applied sinusoidal recharge of either a single period or composite of multiple periods of length similar to known climate cycles. Periodic water level components, reconstructed using singular spectrum analysis (SSA), were used to calibrate the analytical model to estimate each recharge component. The results demonstrated that periodic recharge estimates were most accurate in basins with nearly uniform transmissivity and the accuracy of the recharge estimates depends on monitoring well location. A case study of the San Pedro Basin, Arizona, is presented as an example of calibrating the analytical model to real data.

Water Resources Research↗
USGS2004

Log Pearson type 3 quantile estimators with regional skew information and low outlier adjustments

The recently developed expected moments algorithm (EMA) [ Cohn et al. , 1997 ] does as well as maximum likelihood estimations at estimating log‐Pearson type 3 (LP3) flood quantiles using systematic and historical flood information. Needed extensions include use of a regional skewness estimator and its precision to be consistent with Bulletin 17B . Another issue addressed by Bulletin 17B is the treatment of low outliers. A Monte Carlo study compares the performance of Bulletin 17B using the entire sample with and without regional skew with estimators that use regional skew and censor low outliers, including an extended EMA estimator, the conditional probability adjustment (CPA) from Bulletin 17B , and an estimator that uses probability plot regression (PPR) to compute substitute values for low outliers. Estimators that neglect regional skew information do much worse than estimators that use an informative regional skewness estimator. For LP3 data the low outlier rejection procedure generally results in no loss of overall accuracy, and the differences between the MSEs of the estimators that used an informative regional skew are generally modest in the skewness range of real interest. Samples contaminated to model actual flood data demonstrate that estimators which give special treatment to low outliers significantly outperform estimators that make no such adjustment.

Water Resources Research↗
USGS2004

Interpretation of concentration‐discharge patterns in acid‐neutralizing capacity during storm flow in three small, forested catchments in Shenandoah National Park, Virginia

Episodic concentration‐discharge ( c‐Q ) plots are a popular tool for interpreting the hydrochemical response of small, forested catchments. Application of the method involves assuming an underlying conceptual model of runoff processes and comparing observed c‐Q looping patterns with those predicted by the model. We analyzed and interpreted c‐Q plots of acid‐neutralizing capacity (ANC) for 133 storms collected over a 7‐year period from three catchments in Shenandoah National Park, Virginia. Because of their underlying lithologies the catchments represent a gradient in both hydrologic and geochemical behavior, ranging from a flashy, acidic, poorly buffered catchment to a moderate, neutral, well‐buffered catchment. The relative frequency of observed anticlockwise c‐Q loops in each catchment decreased along this gradient. Discriminant function analysis indicated that prestorm base flow ANC was an important predictor of loop rotation direction; however, the strength of the predictive relationship decreased along the same gradient. The trends were consistent with several equally plausible three‐component mixing models. Uncertainty regarding end‐member timing and relative volume and possible time variation in end‐member concentrations were key factors precluding identification of a unique model. The inconclusive results obtained on this large data set suggest that identification of underlying runoff mechanisms on the basis of a small number of c‐Q plots without additional supporting evidence is likely to be misleading.

Virginia↗
USGS2004

A resampling procedure for generating conditioned daily weather sequences

A method is introduced to generate conditioned daily precipitation and temperature time series at multiple stations. The method resamples data from the historical record “nens” times for the period of interest (nens = number of ensemble members) and reorders the ensemble members to reconstruct the observed spatial (intersite) and temporal correlation statistics. The weather generator model is applied to 2307 stations in the contiguous United States and is shown to reproduce the observed spatial correlation between neighboring stations, the observed correlation between variables (e.g., between precipitation and temperature), and the observed temporal correlation between subsequent days in the generated weather sequence. The weather generator model is extended to produce sequences of weather that are conditioned on climate indices (in this case the Niño 3.4 index). Example illustrations of conditioned weather sequences are provided for a station in Arizona (Petrified Forest, 34.8°N, 109.9°W), where El Niño and La Niña conditions have a strong effect on winter precipitation. The conditioned weather sequences generated using the methods described in this paper are appropriate for use as input to hydrologic models to produce multiseason forecasts of streamflow.

Water Resources Research↗
USGS2005

Numerical simulation of double‐diffusive finger convection

A hybrid finite element, integrated finite difference numerical model is developed for the simulation of double‐diffusive and multicomponent flow in two and three dimensions. The model is based on a multidimensional, density‐dependent, saturated‐unsaturated transport model (SUTRA), which uses one governing equation for fluid flow and another for solute transport. The solute‐transport equation is applied sequentially to each simulated species. Density coupling of the flow and solute‐transport equations is accounted for and handled using a sequential implicit Picard iterative scheme. High‐resolution data from a double‐diffusive Hele‐Shaw experiment, initially in a density‐stable configuration, is used to verify the numerical model. The temporal and spatial evolution of simulated double‐diffusive convection is in good agreement with experimental results. Numerical results are very sensitive to discretization and correspond closest to experimental results when element sizes adequately define the spatial resolution of observed fingering. Numerical results also indicate that differences in the molecular diffusivity of sodium chloride and the dye used to visualize experimental sodium chloride concentrations are significant and cause inaccurate mapping of sodium chloride concentrations by the dye, especially at late times. As a result of reduced diffusion, simulated dye fingers are better defined than simulated sodium chloride fingers and exhibit more vertical mass transfer.

Water Resources Research↗