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 1,405 records · Page 78Linked to original sources

Estimation of groundwater recharge parameters by time series analysis

A model is proposed that relates water level fluctuations in a Dupuit aquifer to effective precipitaton at the top of the unsaturated zone. Effective precipitation, defined herein as that portion of precipitation which becomes recharge, is related to precipitation measured in a nearby gage by a two-parameter function. A second-order stationary assumption is used to connect the spectra of effective precipitation and water level fluctuations. Measured precipitation is assumed to be Gaussian, in order to develop a transfer function that relates the spectra of measured and effective precipitation. A nonlinear least squares technique is proposed for estimating parameters of the effective-precipitation function. Although sensitivity analyses indicate difficulties that may be encountered in the estimation procedure, the methods developed did yield convergent estimates for two case studies.

Water Resources Research↗

Aquifer reclamation design: The use of contaminant transport simulation combined with nonlinear programing

A simulation-management methodology is demonstrated for the rehabilitation of aquifers that have been subjected to chemical contamination. Finite element groundwater flow and contaminant transport simulation are combined with nonlinear optimization. The model is capable of determining well locations plus pumping and injection rates for groundwater quality control. Examples demonstrate linear or nonlinear objective functions subject to linear and nonlinear simulation and water management constraints. Restrictions can be placed on hydraulic heads, stresses, and gradients, in addition to contaminant concentrations and fluxes. These restrictions can be distributed over space and time. Three design strategies are demonstrated for an aquifer that is polluted by a constant contaminant source: they are pumping for contaminant removal, water injection for in-ground dilution, and a pumping, treatment, and injection cycle. A transient model designs either contaminant plume interception or in-ground dilution so that water quality standards are met. The method is not limited to these cases. It is generally applicable to the optimization of many types of distributed parameter systems.

Water Resources Research↗

A nonparametric trend test for seasonal data with serial dependence.

Statistical tests for monotonic trend in seasonal (e.g., monthly) hydrologic time series are commonly confounded by some of the following problems: nonnormal data, missing values, seasonality, censoring (detection limits), and serial dependence. An extension of the Mann-Kendall test for trend (designed for such data) is presented here. Because the test is based entirely on ranks, it is robust against nonnormality and censoring. Seasonality and missing values present no theoretical or computational obstacles to its application. Monte Carlo experiments show that, in terms of type I error, it is robust against serial correlation except when the data have strong long-term persistence (e.g., ARMA (1, 1) monthly processes with ϕ > 0.6) or short records (∼ 5 years). When there is no serial correlation, it is less powerful than a related simpler test which is not robust against serial correlation.

Water Resources Research↗

Double-porosity models for a fissured groundwater reservoir with fracture skin

Theories of flow to a well in a double-porosity groundwater reservoir are modified to incorporate effects of a thin layer of low-permeability material or fracture skin that may be present at fracture-block interfaces as a result of mineral deposition or alteration. The commonly used theory for flow in double- porosity formations that is based upon the assumption of pseudo–steady state block-to-fissure flow is shown to be a special case of the theory presented in this paper. The latter is based on the assumption of transient block-to-fissure flow with fracture skin. Under conditions where fracture skin has a hydraulic conductivity that is less than that of the matrix rock, it may be assumed to impede the interchange of fluid between the fissures and blocks. Resistance to flow at fracture-block interfaces tends to reduce spatial variation of hydraulic head gradients within the blocks. This provides theoretical justification for neglecting the divergence of flow in the blocks as required by the pseudo–steady state flow model. Coupled boundary value problems for flow to a well discharging at a constant rate were solved in the Laplace domain. Both slab-shaped and sphere-shaped blocks were considered, as were effects of well bore storage and well bore skin. Results obtained by numerical inversion were used to construct dimensionless-type curves that were applied to well test data, for a pumped well and for an observation well, from the fractured volcanic rock terrane of the Nevada Test Site.

Water Resources Research↗

Eulerian-Lagrangian solution of the convection-dispersion equation in natural coordinates

The vast majority of numerical investigations of transport phenomena use an Eulerian formulation for the convenience that the computational grids are fixed in space. An Eulerian-Lagrangian method (ELM) of solution for the convection-dispersion equation is discussed and analyzed. The ELM uses the Lagrangian concept in an Eulerian computational grid system. The values of the dependent variable off the grid are calculated by interpolation. When a linear interpolation is used, the method is a slight improvement over the upwind difference method. At this level of approximation both the ELM and the upwind difference method suffer from large numerical dispersion. However, if second-order Lagrangian polynomials are used in the interpolation, the ELM is proven to be free of artificial numerical dispersion for the convection-dispersion equation. The concept of the ELM is extended for treatment of anisotropic dispersion in natural coordinates. In this approach the anisotropic properties of dispersion can be conveniently related to the properties of the flow field. Several numerical examples are given to further substantiate the results of the present analysis.

Water Resources Research↗

On the treatment of evapotranspiration, soil moisture accounting, and aquifer recharge in monthly water balance models

Several two- to six-parameter regional water balance models are examined by using 50-year records of monthly streamflow at 10 sites in New Jersey. These models include variants of the Thornthwaite-Mather model, the Palmer model, and the more recent Thomas abcd model. Prediction errors are relatively similar among the models. However, simulated values of state variables such as soil moisture storage differ substantially among the models, and fitted parameter values for different models sometimes indicated an entirely different type of basin response to precipitation. Some problems in parameter identification are noted, including difficulties in identifying an appropriate time lag factor for the Thornthwaite-Mather-type model for basins with little groundwater storage, very high correlations between upper and lower storages in the Palmer-type model, and large sensitivity of parameter a of the abcd model to bias in estimates of precipitation and potential evapotranspiration. Modifications to the threshold concept of the Thornthwaite-Mather model were statistically valid for the six stations in northern New Jersey. The abcd model resulted in a simulated seasonal cycle of groundwater levels similar to fluctuations observed in nearby wells but with greater persistence. These results suggest that extreme caution should be used in attaching physical significance to model parameters and in using the state variables of the models in indices of drought and basin productivity.

Water Resources Research↗

Origin and distribution of carbon dioxide in the unsaturated zone of the southern High Plains of Texas

Partial pressures of CO 2 , O 2 , N 2 , and Ar were monitored at two locations in the Ogallala aquifer system on the Southern High Plains of Texas. Samples were collected monthly during parts of 1980–1981 from nine depths ranging from 0.6 to 36 meters below land surface. P CO 2 was observed to be greater at depth than in the active soil zone and thus appears to contradict the normal process in which CO 2 is generated in the soil zone and diffuses upward to the atmosphere and downward to the water table. The δ 13 C of the CO 2 gas was quite uniform and averaged −17.9 per mil. P O 2 declined with depth, suggesting in situ generation of CO 2 by the oxidation of carbon. Several hypotheses were considered to explain the origin of the CO 2 at depth. It was concluded that the most probable hypothesis was that dissolved and particulate organic carbon introduced by recharging water was oxidized to CO 2 by the aerobic microbial community that utilized oxygen diffusing in from the atmosphere. This hypothesis is consistent with the CO 2 concentration profile, calculated production profile of CO 2 , δ 13 C values of CO 2 gas, caliche, soil humic acid fraction, and dissolved carbonate in groundwater. The abundance of CO 2 , its concentration profile, and its probable origin provide information for evaluating the observed complex sequence of caliche dissolution and precipitation known to occur in the aquifer.

Texas↗

Note on the applicability of the James-Stein Estimator in regional hydrologic studies

The applicability of the James-Stein estimator in regional hydrologic studies which entail the estimation of an N -dimensional location parameter is discussed. Regional studies are frequently characterized by relatively short, generally correlated, samples drawn from nonsymmetric and bounded, i.e., nonnormal, populations. By means of computer simulation studies the James-Stein estimator, subject to the Lindley modification and adoption of the positive part rule suggested by Efron and Morris and conditioned on the assumption of independence, was shown to be robust in the case of the hydrologically plausible distribution considered here, namely, Weibull distributions with coefficient of skewness ranging from 0 to 10. However, in contrast to traditional methods of regionalization the effect of cross correlation is a diminishment of the relative risk advantage of the James-Stein estimator, even in the best case of normal variables : this is discussed and illustrated.

Water Resources Research↗

Interactions of solutes and streambed sediment: 1. An experimental analysis of cation and anion transport in a mountain stream

An experimental injection was performed to study the transport of stream water solutes under conditions of significant interaction with streambed sediments in a mountain pool-and-riffle stream. Experiments were conducted in Little Lost Man Creek, Humboldt County, California, in a period of low flow duringwhich only a part of the bank-full channel held active surface flow. The injection of chloride and several trace cations lasted 20 days. In this report we discuss the results of the first 24 hours of the injection and survey the results of the first 10 days. Solute-streambed interactions of two types were observed. First, the physical transport of the conservative tracer, chloride, was affected by intergravel flow and stagnant watt, zones created by the bed relief. Second, the transport of the cations (strontium, potassium, and lithium) was appreciably modified by sorption onto streambed sediment. In the stream the readily observable consequence of the solute-streambed interactions was an attenuation of the dissolved concentration of each of the tracers. The attenuation in the stream channel occurred concurrently with the storage of tracers in the streambed via both physical and chemical processes. All tracers were subsequently present in shallow wells dug several meters from the wetted part of the channel. Sediment samples collected approximately 3 weeks after the start of the injection contained increased concentrations of the injected cations.

California↗

Interactions of solutes and streambed sediment: 2. A dynamic analysis of coupled hydrologic and chemical processes that determine solute transport

Solute transport in streams is determined by the interaction of physical and chemical processes. Data from an injection experiment for chloride and several cations indicate significant influence of solutestreambed processes on transport in a mountain stream. These data are interpreted in terms of transient storage processes for all tracers and sorption processes for the cations. Process parameter values are estimated with simulations based on coupled quasi-two-dimensional transport and first-order mass transfer sorption. Comparative simulations demonstrate the relative roles of the physical and chemical processes in determining solute transport. During the first 24 hours of the experiment, chloride concentrations were attenuated relative to expected plateau levels. Additional attenuation occurred for the sorbing cation strontium. The simulations account for these storage processes. Parameter values determined by calibration compare favorably with estimates from other studies in mountain streams. Without further calibration, the transport of potassium and lithium is adequately simulated using parameters determined in the chloride-strontium simulation and with measured cation distribution coefficients.

Water Resources Research↗

The effect of glaciers on streamflow variations

The effect of temperate glaciers on runoff variations is examined for the North Cascade Mountains of Washington State. The principal influences of glaciers on streamflow are often unexpected contributions to streamflow volume, a delay of the maximum seasonal flow, and a decrease in annual and monthly variation of runoff. The delay of maximum flow is caused by temporary englacial storage of spring meltwater and by peak meltwater production occurring in midsummer. The englacial storage, for one case, is 54% of the potential May runoff. An algorithm is presented that calculates the coefficient of variation of runoff for any arbitrary glacier cover. The results suggest that a minimum in year-to-year variation occurs for basins about 36% glacierized. On a month-to-month basis, maximum variation occurs in July and August for basins with less than 10% glacier cover but is a minimum for basins with glacier covers greater than 30%.

Water Resources Research↗

Water balance models in one-month-ahead streamflow forecasting

T echniques are tested that incorporate information from water balance models in making 1-month-ahead streamflow forecasts in New Jersey. The results are compared to those based on simple autoregressive time series models. The relative performance of the models is dependent on the month of the year in question. The water balance models are most useful for forecasts of April and May flows. For the stations in northern New Jersey, the April and May forecasts were made in order of decreasing reliability using the water-balance-based approaches, using the historical monthly means, and using simple autoregressive models. The water balance models were useful to a lesser extent for forecasts during the fall months. For the rest of the year the improvements in forecasts over those obtained using the simpler autoregressive models were either very small or the simpler models provided better forecasts. When using the water balance models, monthly corrections for bias are found to improve minimum mean-square-error forecasts as well as to improve estimates of the forecast conditional distributions.

Water Resources Research↗

Interannual streamflow variability in the United States based on principal components

Interannual modes of streamflow variation at 106 locations across the United States during the period 1931–1978 are defined by using principal components. Five statistically significant components are found to account for more than 56% of the total streamflow variance. The first principal component represents a nationwide tendency for either above- or below-mean streamflow. The second component represents a north-south opposition in departures from mean flow, and the third, an east-west opposition. Higher-order components (fourth and fifth) geographically depict regional patterns of opposition in the sign of streamflow departures between coastal-continental areas and between the northern and southern plains, respectively. Analyses using spatially and temporally modified data sets indicate that the first three components (which explain 45% of the variance) are quite stable spatially, while only the first component is stable temporally. Time series analysis of principal component scores indicates that all but the fourth component are first-order autoregressive processes, as is mean annual nationwide streamflow. The fourth component is an autoregressive (AR)(2) process. In general, the principal components of streamflow are found to exhibit more persistence over annual time scales than the mean annual flow data themselves.

Water Resources Research↗

Nitrogen and phosphorus speciation and flux in a large Florida river wetland System

Hydrologic measurements and analyses of various nitrogen and phosphorus species were made on the Apalachicola River system in northern Florida in 1979 and 1980. Annual outflows of total nitrogen (TN) and total phosphorus (TP) were not substantially different from annual inflows. However, there was significant net import of ammonia and soluble reactive phosphorus and net export of some particulate and organic species. The TN: TP ratio ranged from 12 to 15, but the specific ratio of dissolved inorganic nitrogen: soluble reactive phosphorus was much higher (up to 40) and increased in a downstream direction; this contributed to a phosphorus-limiting situation in Apalachicola estuary. Processes within the flood plain ecosystem accounted for much of the release of organic and particulate species and retention of inorganic species. This flood plain function is probably critical for maintaining a nutrient pool in the estuary which supports secondary productivity and a detrital-based food web.

Florida↗

Unit hydrograph approximations assuming linear flow through topologically random channel networks

The instantaneous unit Hydrograph (IUH) of a drainage basin is derived in terms of fundamental basin characteristics ( Z , α, β), where α parameterizes the link (channel segment) length distribution, and β is a vector of hydraulic parameters, Z is one of three basin topological properties, N , ( N , D ), or ( N , M ), where N is magnitude (number of first-order streams), D is diameter (mainstream length), and M is order. The IUH is derived based on assumptions that the links are independent and identically distributed random variables and that the network is a member of a topologically random population. Linear routing schemes, including translation, diffusion, and general linear routing are used, and constant drainage density is assumed. By using ( N , α, β) as the fundamental basin characteristics, asymptotic (for large N ) considerations lead to a Weibull probability density function for the IUH, with time to peak given by t p = (2 N ) ½ α * /β * where α * is mean link length, and β * is a scalar hydraulic parameter (usually average celerity). This asymptotic IUH is identical for all linear routing schemes.

Water Resources Research↗

A quantitative analysis of the Lassen hydrothermal system, north central California

Our conceptual model of the Lassen system is termed a liquid-dominated hydrothermal system with a parasitic vapor-dominated zone. The essential feature of this model is that steam and steam-heated discharge at relatively high altitudes in Lassen Volcanic National Park (LVNP) and liquid discharge with high chloride concentrations at relatively low altitudes outside LVNP are both fed by an upflow of high-enthalpy two-phase fluid within the Park. Liquid flows laterally away from the upflow area toward the areas of high-chloride discharge, and steam rises through a vapor-dominated zone to feed the steam and steam-heated features. Numerical simulations show that several conditions are necessary for the development of this type of system, including (1) large-scale topographic relief; (2) an initial period of convective heating within an upflow zone followed by (3) a change in hydrologic or geologic conditions that initiates drainage of liquid from portions of the upflow zone; and (4) low-permeability barriers that inhibit the movement of cold water into the vapor zone. Simulations of thermal fluid withdrawal south of LVNP, carried out in order to determine the effects of such withdrawal on portions of the hydrothermal system within the Park, generally showed decreases in pressure and liquid saturation beneath the vapor zone which resulted in temporary increases and subsequent decreases in the rate of upflow of steam. A generalized production-injection scenario that could mitigate the effects of development on both the high-chloride and steam-fed features was identified.

California↗