Search

We report the results of field and laboratory investigations of stream-aquifer interactions in a watershed along the California coast to assess the impact of groundwater pumping for irrigation on stream flows. The methods used include subsurface sediment sampling using direct-push drilling, laboratory permeability and particle size analyses of sediment, piezometer installation and instrumentation, stream discharge and stage monitoring, pumping tests for aquifer characterization, resistivity surveys, and long-term passive monitoring of stream stage and groundwater levels. Spectral analysis of long-term water level data was used to assess correlation between stream and groundwater level time series data. The investigations revealed the presence of a thin low permeability silt-clay aquitard unit between the main aquifer and the stream. This suggested a three layer conceptual model of the subsurface comprising unconfined and confined aquifers separated by an aquitard layer. This was broadly confirmed by resistivity surveys and pumping tests, the latter of which indicated the occurrence of leakage across the aquitard. The aquitard was determined to be 2–3 orders of magnitude less permeable than the aquifer, which is indicative of weak stream-aquifer connectivity and was confirmed by spectral analysis of stream-aquifer water level time series. The results illustrate the importance of site-specific investigations and suggest that even in systems where the stream is not in direct hydraulic contact with the producing aquifer, long-term stream depletion can occur due to leakage across low permeability units. This has implications for management of stream flows, groundwater abstraction, and water resources management during prolonged periods of drought.

A semianalytical solution is presented for transient streaming potentials associated with flow to a pumping well in an unconfined aquifer, taking into account the effect of flow in the unsaturated zone above the water table. Flow in the unsaturated zone is modeled with a linearized form of Richards' equation using an exponential model for soil moisture retention and unsaturated hydraulic conductivity. Archie's law is invoked for unsaturated electrical conductivity. The unsaturated electrokinetic coupling coefficient is modeled with a decaying exponential, where the maximum value is at and below the water table. The coupled flow and electrokinetic problem is solved using Laplace and Hankel transforms. The results of the model predicted behavior are presented and compared to that observed in laboratory simulations of pump­ ing tests. The early time polarity reversal predicted the model is observable in the experiments. Other non- monotonic streaming potential behaviors predicted by the model are also evident in experimental measurements. The model is used to estimate hydraulic parameters from SP data and these compare well to those obtained from drawdown data. For example, a hydraulic conductivity of 3.6 3 10 24 m/s is obtained from SP data compared to 3.4 3 10 24 m/s from drawdown data.

We present transient streaming potential data collected during falling-head permeameter tests performed on samples of two sands with different physical and chemical properties. The objective of the work is to estimate hydraulic conductivity (K) and the electrokinetic coupling coefficient (Cl) of the sand samples. A semi-empirical model based on the falling-head permeameter flow model and electrokinetic coupling is used to analyze the streaming potential data and to estimate K and Cl. The values of K estimated from head data are used to validate the streaming potential method. Estimates of K from streaming potential data closely match those obtained from the associated head data, with less than 10% deviation. The electrokinetic coupling coefficient was estimated from streaming potential vs. (1) time and (2) head data for both sands. The results indicate that, within limits of experimental error, the values of Cl estimated by the two methods are essentially the same. The results of this work demonstrate that a temporal record of the streaming potential response in falling-head permeameter tests can be used to estimate both K and Cl. They further indicate the potential for using transient streaming potential data as a proxy for hydraulic head in hydrogeology applications.

A modified version of a published slug test model for unconfined aquifers is applied to cross-hole slug test data collected in field tests conducted at the Widen site in Switzerland. The model accounts for water-table effects using the linearised kinematic condition. The model also accounts for inertial effects in source and observation wells. The primary objective of this work is to demonstrate applicability of this semi-analytical model to multi-well and multi-level pneumatic slug tests. The pneumatic perturbation was applied at discrete intervals in a source well and monitored at discrete vertical intervals in observation wells. The source and observation well pairs were separated by distances of up to 4 m. The analysis yielded vertical profiles of hydraulic conductivity, specific storage, and specific yield at observation well locations. The hydraulic parameter estimates are compared to results from prior pumping and single-well slug tests conducted at the site, as well as to estimates from particle size analyses of sediment collected from boreholes during well installation. The results are in general agreement with results from prior tests and are indicative of a sand and gravel aquifer. Sensitivity analysis show that model identification of specific yield is strongest at late-time. However, the usefulness of late-time data is limited due to the low signal-to-noise ratios., 6 figures, 2 tables

Two models for slug tests conducted in unconfined aquifers are developed by (a) extending the unconfined KGS solution to oscillatory responses, yielding a model referred to herein as the unified model, and (b) replacing the constant head condition with the linearized kinematic condition at the water table. The models can be used to analyze the full range of responses from highly oscillatory to overdamped. The second model, refered to as the moving water table (MWT) model, is only applicable when effects of well bore skin are negligible. The models are validated by comparison with published solutions, and by application to a published case study of field tests conducted in wells without skin in an unconfined aquifer at the MSEA site in Nebraska. In this regard (a) the MWT model essentially yields the same results as the confined KGS model, except very close to the water table, and (b) the unified model yields slightly smaller aquifer K-values relative to the MWT model at all positions in the well. All model solutions yield comparable results when fitted to published field data obtained in an unconfined fluvial aquifer at the MSEA site in Nebraska. The unified model is fitted to field data collected in wells known to exhibit positive skin effects at the Boise Hydrogeophysical Research Site (BHRS) in Boise, Idaho. It is shown to yield hydraulic conductivity estimates of comparable magnitude to those obtained with the KGS model for overdamped responses, and the Springer–Gelhar model for oscillatory responses. Sensitivity of the MWT model to specific yield, S y , and hydraulic anisotropy, j is evaluated and the results, when plotted in log–log space and with consideration of log-scale time derivatives of the response, indicate that these two parameters should be estimable from slug test data, though challenges still remain.

Summary A recently developed unified model for partially-penetrating slug tests in unconfined aquifers ( Malama et al., in press ) provides a semi-analytical solution for aquifer response at the wellbore in the presence of inertial effects and wellbore skin, and is able to model the full range of responses from overdamped/monotonic to underdamped/oscillatory. While the model provides a unifying framework for realistically analyzing slug tests in aquifers (with the ultimate goal of determining aquifer properties such as hydraulic conductivity K and specific storage S s ), it is currently unclear whether parameters of this model can be well-identified without significant prior information and, thus, what degree of information content can be expected from such slug tests. In this paper, we examine the information content of slug tests in realistic field scenarios with respect to estimating aquifer properties, through analysis of both numerical experiments and field datasets. First, through numerical experiments using Markov Chain Monte Carlo methods for gauging parameter uncertainty and identifiability, we find that: (1) as noted by previous researchers, estimation of aquifer storage parameters using slug test data is highly unreliable and subject to significant uncertainty; (2) joint estimation of aquifer and skin parameters contributes to significant uncertainty in both unless prior knowledge is available; and (3) similarly, without prior information joint estimation of both aquifer radial and vertical conductivity may be unreliable. These results have significant implications for the types of information that must be collected prior to slug test analysis in order to obtain reliable aquifer parameter estimates. For example, plausible estimates of aquifer anisotropy ratios and bounds on wellbore skin K should be obtained, if possible, a priori. Secondly, through analysis of field data – consisting of over 2500 records from partially-penetrating slug tests in a heterogeneous, highly conductive aquifer, we present some general findings that have applicability to slug testing. In particular, we find that aquifer hydraulic conductivity estimates obtained from larger slug heights tend to be lower on average (presumably due to non-linear wellbore losses) and tend to be less variable (presumably due to averaging over larger support volumes), supporting the notion that using the smallest slug heights possible to produce measurable water level changes is an important strategy when mapping aquifer heterogeneity. Finally, we present results specific to characterization of the aquifer at the Boise Hydrogeophysical Research Site. Specifically, we note that (1) K estimates obtained using a range of different slug heights give similar results, generally within ±20%; (2) correlations between estimated K profiles with depth at closely-spaced wells suggest that K values obtained from slug tests are representative of actual aquifer heterogeneity and not overly affected by near-well media disturbance (i.e., “skin”); (3) geostatistical analysis of K values obtained indicates reasonable correlation lengths for sediments of this type; and (4) overall, K values obtained do not appear to correlate well with porosity data from previous studies.

Summary An alternative linearization of the nonlinear kinematic condition at the water table used in unconfined aquifer flow theory is proposed. It simulates the downward propagation of the water table, during water extraction via a pumping well, as a diffuse process controlled by a dimensionless linearization parameter, β D . The limiting value of β D = 0 corresponds to the solution of Neuman (1972) whereas that of β D → ∞ corresponds to the early-time Theis solution. Using the solution obtained with this linearization, data collected during a pumping test conducted at the Boise Hydrogeophysical Research Site in Idaho, US, were analyzed. The estimates of hydraulic conductivity and specific storage obtained in this study compare well with those obtained with the other models. The most significant improvement is in the estimated value of specific yield, with the proposed model yielding values that are close to those expected from known porosity values ( n ∼ 0.25) at the site. The results show that meaningful estimates of specific yield can be obtained by treating the water table as a diffuse moving material boundary, without recourse to more complex models for water table decline or inclusion of unsaturated zone flow.

Summary We investigate the response of the saturated zone (unconfined aquifer) to evapotranspiration (ET) flux at ground surface. We neglect fluid flow and storage in the unsaturated zone and treat ET as a sinusoidal forcing function at the watertable. The linearized kinematic condition is imposed at the watertable. Analytical solutions are developed for the case of flow in a domain of (a) semi-infinite extent to simulate response in a domain bounded by a river and (b) infinite lateral extent to simulate the response in a domain with no river boundaries. These solutions are fitted to observed groundwater head fluctuations recorded in observation wells at the Boise Hydrogeophysical Research Site in Idaho and the Larned Research Site in Kansas. Estimates of the amplitude of the ET flux, aquifer hydraulic conductivity, specific storage and specific yield are obtained and these compare well to published results from pumping tests conducted at the site. The field exercise is used to explore the potential for using groundwater head fluctuations to estimate ET and hydraulic parameters of unconfined aquifers.

Summary An analytical solution for flow in the vicinity of an observation well is developed. The observation is emplaced in a homogeneous aquifer of infinite radial extent that is pumped at a constant rate, and satisfies the Theis solution. We attempt to account for the effect of an observation well on the drawdown response in a neighborhood of finite radial extent centered around the well. The model of Black and Kipp [Black, J.H., Kipp Jr., K.L., 1977. Observation well response time and its effect upon aquifer tests. Journal of Hydrology 34, 297–306] only makes a correction to the drawdown response at the well location and does not modify the flow pattern in the neighborhood of the well. This may be sufficient for isolated observation wells but, for well-fields or sites with relatively closely spaced wells, the flow patterns in the neighborhoods of the wells may significantly impact drawdown response. The model we develop is applied to field data and is shown to yield an aquifer hydraulic diffusivity that is comparable to that estimated with the model of Black and Kipp [Black, J.H., Kipp Jr., K.L., 1977. Observation well response time and its effect upon aquifer tests. Journal of Hydrology 34, 297–306]. In addition to providing estimates of formation conductivity and specific storage, the solution yield estimates of the hydraulic conductivity and specific storage of the zone of influence and has the ability to predict the effect of the observation well on drawdown response in a finite region of influence in the vicinity of the well.

Summary A semi-analytical solution is presented for the problem of flow in a system consisting of unconfined and confined aquifers, separated by an aquitard. The unconfined aquifer is pumped continuously at a constant rate from a well of infinitesimal radius that partially penetrates its saturated thickness. The solution is termed semi-analytical because the exact solution obtained in double Laplace‐Hankel transform space is inverted numerically. The solution presented here is more general than similar solutions obtained for confined aquifer flow as we do not adopt the assumption of unidirectional flow in the confined aquifer (typically assumed to be horizontal) and the aquitard (typically assumed to be vertical). Model predicted results show significant departure from the solution that does not take into account the effect of leakage even for cases where aquitard hydraulic conductivities are two orders of magnitude smaller than those of the aquifers. The results show low sensitivity to changes in radial hydraulic conductivities for aquitards that are two or more orders of magnitude smaller than those of the aquifers, in conformity to findings of earlier workers that radial flow in aquitards may be neglected under such conditions. Hence, for cases were aquitard hydraulic conductivities are two or more orders of magnitude smaller than aquifer conductivities, the simpler models that restrict flow to the radial direction in aquifers and to the vertical direction in aquitards may be sufficient. However, the model developed here can be used to model flow in aquifer‐aquitard systems where radial flow is significant in aquitards. a 2008 Elsevier B.V. All rights reserved.

We develop a closed-form three-parameter model for unsaturated hydraulic conductivity associated with the Kosugi three-parameter lognormal moisture retention model. The model derivation uses a slight modification to Mualem's theory, which is nearly exact for non-clay soils. Kosugi's three-parameter lognormal moisture retention model uses physically meaningful parameters, but a corresponding closed-form relative hydraulic conductivity model has never been developed. The model is further extended to a four -parameter model by truncating the underlying pore size distribution at physically permissible minimum and maximum pore radii. The proposed closed-form models are fitted to well-known experimental data, to illustrate their utility. They have the same physical basis as Kosugi's two-parameter model, but are more general., 1 table, 5 figures

International audience; We consider the transient streaming potential response due to pumping from a confined aquifer through a fully penetrating line sink. Confined aquifer flow is assumed to occur without fluid leakage from the confining units. However, since confining units are typically clayey, and hence more electrically conductive than the aquifer, they are treated as non-insulating in our three-layer conceptual model. We develop a semi-analytical solution for the transient streaming potentials response of the aquifer and the confining units to pumping of the aquifer. The solution is fitted to field measurements of streaming potentials associated with an aquifer test performed at a site located near Montalto Uffugo, in the region of Calabria in Southern Italy. This yields an average hydraulic conductivity that compares well to the estimate obtained using only hydraulic head data. Specific storage is estimated with greater estimation uncertainty than hydraulic conductivity and is significantly smaller than that estimated from hydraulic head data. This indicates that specific storage may be a more difficult parameter to estimate from streaming potential data. The mismatch may also be due to the fact that only recovery streaming potential data were used here whereas head data for both production and recovery were used. The estimate from head data may also constitute an upper bound since head data were not corrected for pumping and observation wellbore storage. Estimated values of the electrical conductivities of the confining units compare well to those estimated using electrical resistivity tomography. Our work indicates that, where observation wells are unavailable to provide more direct estimates, streaming potential data collected at land surface may, in principle, be used to provide preliminary estimates of aquifer hydraulic conductivity and specific storage, where the latter is estimated with greater uncertainty than the former.

International audience; The importance of estimating spatially variable aquifer parameters such as transmissivity is widely recognized for studies in resource evaluation and contaminant transport. A useful approach for mapping such parameters is inverse modeling of data from series of pumping tests, that is, via hydraulic tomography. This inversion of field hydraulic tomographic data requires development of numerical forward models that can accurately represent test conditions while maintaining computational efficiency. One issue this presents is specification of boundary and initial conditions, whose location, type, and value may be poorly constrained. To circumvent this issue when modeling unconfined steady-state pumping tests, we present a strategy that analyzes field data using a potential difference method and that uses dipole pumping tests as the aquifer stimulation. By using our potential difference approach, which is similar to modeling drawdown in confined settings, we remove the need for specifying poorly known boundary condition values and natural source/sink terms within the problem domain. Dipole pumping tests are complementary to this strategy in that they can be more realistically modeled than single-well tests due to their conservative nature, quick achievement of steady state, and the insensitivity of near-field response to far-field boundary conditions. After developing the mathematical theory, our approach is first validated through a synthetic example. We then apply our method to the inversion of data from a field campaign at the Boise Hydrogeophysical Research Site. Results from inversion of nine pumping tests show expected geologic features, and uncertainty bounds indicate that hydraulic conductivity is well constrained within the central site area.

International audience; We present a semi-analytical solution for the transient streaming potential response of an unconfined aquifer to continuous constant rate pumping. We assume that flow occurs without leakage from the unit below a transverse anisotropic aquifer and neglect flow in the unsaturated zone by treating the water-table as a moving material boundary. In the development of the solution to the streaming potential problem, we impose insulating boundary conditions at land surface and the lower boundary of the lower confining unit. We solve the problem exactly in the double Laplace–Hankel transform space and obtain the inverse transforms numerically. The solution is used to analyse transient streaming potential data collected during dipole hydraulic tests conducted at the Boise Hydrogeophysical Research Site in 2007 June. This analysis yields estimates of aquifer hydraulic parameters. The estimated hydraulic parameters, namely, hydraulic conductivity, transverse hydraulic anisotropy, specific storage and specific yield, compare well to published estimates obtained by inverting drawdown data collected at the field site.

We present a nonlinear stochastic inverse algorithm that allows conditioning estimates of transient hydraulic heads, fluxes and their associated uncertainty on information about hydraulic conductivity (K) and hydraulic head (h) data collected in a randomly heterogeneous confined aquifer. Our algorithm is based on Laplace-transformed recursive finite-element approximations of exact nonlocal first and second conditional stochastic moment equations of transient flow. It makes it possible to estimate jointly spatial variations in natural log-conductivity ( Y = ln K ) , the parameters of its underlying variogram, and the variance–covariance of these estimates. Log-conductivity is parameterized geostatistically based on measured values at discrete locations and unknown values at discrete “pilot points”. Whereas prior values of Y at pilot point are obtained by generalized kriging, posterior estimates at pilot points are obtained through a maximum likelihood fit of computed and measured transient heads. These posterior estimates are then projected onto the computational grid by kriging. Optionally, the maximum likelihood function may include a regularization term reflecting prior information about Y. The relative weight assigned to this term is evaluated separately from other model parameters to avoid bias and instability. We illustrate and explore our algorithm by means of a synthetic example involving a pumping well. We find that whereas Y and h can be reproduced quite well with parameters estimated on the basis of zero-order mean flow equations, all model quality criteria identify the second-order results as being superior to zero-order results. Identifying the weight of the regularization term and variogram parameters can be done with much lesser ambiguity based on second- than on zero-order results. A second-order model is required to compute predictive error variances of hydraulic head (and flux) a posteriori. Conditioning the inversion jointly on conductivity and hydraulic head data results in lesser predictive uncertainty than conditioning on conductivity or head data alone.