Search

Mathematical models for stream depletion with stream stage decline or drawdown are developed to overcome the deficiency in existing models that typically use the constant-head (Dirichlet) or general (Robin) boundary condition and source terms at the stream-aquifer interface. Existing approaches assume a fixed stream stage during pumping, implies that the stream is an infinite water source, with depletion defined as a decrease in stream discharge. We refer to this depletion without drawdown as the ``stream depletion paradox.'' It is a glaring model limitation, ignoring the most observable adverse effect of long-term groundwater abstraction near a stream, namely stage declines that eventually lead to dry streambeds. Field data are presented to demonstrate that stream stage responds to pumping near the stream, motivating the development of an alternative theory predicts transient stream drawdown based on the concepts of finite stream storage and mass continuity at the stream-aquifer interface. Based on this alternative theory, models are developed for the cases of a non- and a fully-penetrating stream. The proposed model reduces to the fixed-stage model in the limit as stream storage becomes infinitely large and to the limiting case of confined aquifer flow with a no-flow boundary at the streambed when the stream storage vanishes. The model is applied to field observations of both aquifer and stream drawdown from tests conducted in a confined aquifer over which a shallow stream flows. Model fits and parameter estimates are obtained both aquifer and stream drawdown data. Model predicted and observed transient drawdown behavior indicate that fixed-stage models (a) underestimate late-time aquifer drawdown and (b) overestimate the available recharge from streams to pumping wells. This has significant implications for the sustainable management of water resources in hydraulically connected stream-aquifer systems with heavy groundwater abstraction.

The results reported in this work are based in part on measurements of sap flow in a few select trees on a representative riparian forest plot coupled with a forest-wide randomized sampling of tree sapwood area in a watershed located along the Pacific coast in Santa Cruz County, California. These measurements were upscaled to estimate evapotranspiration (ET) across the forest and to quantify groundwater usage by dominant phreatophyte vegetation. Canopy cover in the study area is dominated by red alder (Alnus rubra) and arroyo willow (Salix lasiolepis), deciduous phreatophyte trees from which a small sample was selected for instrumentation with sap flow sensors on a single forest plot. These localized sap flow measurements were then upscaled to the entire riparian forest to estimate forest ET using data from a survey of sapwood area on six plots scattered randomly across the entire forest. The estimated canopy-scale ET was compared to reference ET and NDVI based estimates. The results show positive correlation between sap flow based estimates and those of the other two methods, though over the winter months, sap flow-based ET values were found to significantly underestimate ET as predicted by the other two methods. The results illustrate the importance of ground-based measurements of sap flow for calibrating satellite based methods and for providing site-specific estimates and to better characterize the ET forcing in groundwater flow models.

We present an approach to uncoupling the pair of transient governing equations used in electrokinetics (i.e., streaming potential and electroosmosis). This approach allows for the solution of two u...

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 multiporosity extension of classical double and triple porosity fractured rock flow models for slightly compressible fluids is presented. The multiporosity model is an adaptation of the multirate solute transport model of Haggerty and Gorelick (1995) to viscous flow in fractured rock reservoirs. It is a generalization of both pseudo-steady-state and transient interporosity flow double porosity models. The model includes a fracture continuum and an overlapping distribution of multiple rock matrix continua, whose fracture-matrix exchange coefficients are specified through a discrete probability mass function. Semi-analytical cylindrically symmetric solutions to the multiporosity mathematical model are developed using the Laplace transform to illustrate its behavior. The multiporosity model presented here is conceptually simple, yet flexible enough to simulate common conceptualizations of double and triple porosity flow. This combination of generality and simplicity makes the multiporosity model a good choice for flow in low-permeability fractured rocks., 5 figures, 3 tables

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.

Model applicability to core-scale solute transport is evaluated using breakthrough data from column experiments conducted with conservative tracers tritium (H-3) and sodium-22, and the retarding solute uranium-232. The three models considered are single-porosity, double-porosity with single-rate mobile-immobile mass-exchange, and the multirate model, which is a deterministic model that admits the statistics of a random mobile-immobile mass-exchange rate coefficient. The experiments were conducted on intact Culebra Dolomite core samples. Previously, data were analyzed using single- and double-porosity models although the Culebra Dolomite is known to possess multiple types and scales of porosity, and to exhibit multirate mobile-immobile-domain mass transfer characteristics at field scale. The data are reanalyzed here and null-space Monte Carlo analysis is used to facilitate objective model selection. Prediction (or residual) bias is adopted as a measure of the model structural error. The analysis clearly shows single- and double-porosity models are structurally deficient, yielding late-time residual bias that grows with time. On the other hand, the multirate model yields unbiased predictions consistent with the late-time -5/2 slope diagnostic of multirate mass transfer. The analysis indicates the multirate model is better suited to describing core-scale solute breakthrough in the Culebra Dolomite than the other two models.

[1] Riparian zones in semiarid regions often exhibit high rates of evapotranspiration (ET) in spite of low-soil moisture content due to the presence of phreatophytic vegetation that is able to withdraw water from shallow aquifers. This work seeks to better define the relationship between ET, the saturated zone and the river boundary by comparing observed water table drawdown records to analytically modeled drawdown in fully penetrating wells of an unconfined aquifer in response to daily ET flux. ET at the Boise Hydrogeophysical Research Site (BHRS), a riparian zone in a temperate, semiarid environment, is calculated using a radiation-based method to provide ET values at four different wells with different vegetation densities. Analytically modeled drawdown response to ET forcing shows that drawdown magnitude increases with increasing distance from the river edge even as the surficial ET forcing remains constant. This behavior is also observed in well hydrographs and shows the buffering effect that flow from the river has on drawdown in fully penetrating riparian wells in highly permeable, unconfined aquifers. Relative contributions of river water to aquifer storage are calculated for ET-induced diurnal fluctuations of the water table at increasing distances from the river boundary. Failure to account for these spatial differences in drawdown related to the river source may explain some errors associated with estimating ET from well hydrographs alone.

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

A new semi-empirical model that can be used to predict fracture deformation behavior under normal compressive loading is presented. The development of a simple exponential model is presented first after which a modified and more general exponential model, with an additional degree of freedom in the model parameters, is obtained. The simple and the modified exponential models are then compared to available fracture closure models, namely the empirical Barton–Bandis hyperbolic model, and a power-law model based on Hertzian contact theory, to determine how good they fit the results of fracture closure experiments under monotonically increasing normal compressive loading. A new parameter called the half-closure stress, σ 1/2 , is introduced and is used, in addition to the maximum fracture closure, Δ v m , in the model fitting procedures for the Barton–Bandis and the simple exponential model. The half-closure stress is shown to be related to the initial normal stiffness, K ni , used in the original Barton–Bandis model. An additional parameter, n, is used in fitting the modified exponential model to the experimental data. Of the models presented herein, the modified exponential model was found to provide the best fit to the experimental data, for the same values of σ 1/2 and Δ v m , over the entire range of compressive stresses. The power-law model based on Hertzian contact theory was found to be unsuitable for predicting normal fracture deformation behavior.

Laboratory experiments and numerical simulations, using Particle Flow Code (PFC3D ), were performed to study the behavior of jointed blocks of model material under uniaxial loading. The effect of joint geometry parameters on the uniaxial compressive strength of jointed blocks was investigated and this paper presents the results of the experiments and numerical simulations. The fracture tensor component in a given direction is used to quantify the combined directional effect of joint geometry parameters including joint density, orientation and size distributions, and the number of joint sets. The variation of the uniaxial compressive strength of the jointed blocks of the model material with the fracture tensor component was investigated. Both the laboratory experiments and the numerical simulations showed that the uniaxial block strength decreases in a nonlinear manner with increasing values of the fracture tensor component. It was observed that joint geometry configuration controls the mode of failure of the jointed blocks and three modes of failure were identified, namely (a) tensile splitting through the intact material, (b) failure by sliding along the joint plane and/or by displacement normal to the joint plane and, (c) mixed mode failure involving both the failure mechanisms in (a) and (b). It has also been shown that with careful parameter calibration procedures, PFC3D could be used to model the strength behavior of jointed blocks of rock.

This report is a summary of the physical processes, primary governing equations, solution approaches, and historic testing related to brine migration in geologic salt. Although most information presented in this report is not new, we synthesize a large amount of material scattered across dozens of laboratory reports, journal papers, conference proceedings, and textbooks. We present a mathematical description of the governing brine flow mechanisms in geologic salt. We outline the general coupled thermal, multi-phase hydrologic, and mechanical processes. We derive these processes' governing equations, which can be used to predict brine flow. These equations are valid under a wide variety of conditions applicable to radioactive waste disposal in rooms and boreholes excavated into geologic salt.

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; 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.

International audience; An algorithm is developed to interpret self-potential (SP) data in terms of distribution of Darcy velocity of the ground water. The model is based on the proportionality existing between the streaming current density and the Darcy velocity. Because the inverse problem of current density determination from SP data is underdetermined, we use Tikhonov regularization with a smoothness constraint based on the differential Laplacian operator and a prior model. The regularization parameter is determined by the L-shape method. The distribution of the Darcy velocity depends on the localization and number of non-polarizing electrodes and information relative to the distribution of the electrical resistivity of the ground. A priori hydraulic information can be introduced in the inverse problem. This approach is tested on two synthetic cases and on real SP data resulting from infiltration of water from a ditch.