Usefulness of flux measurements during hydraulic tomographic survey for mapping hydraulic conductivity distribution in a fractured medium

Using the first-order analysis, we investigate the spatial cross-correlation between hydraulic conductivity variation and specific discharge (flux) as well as its components measured in a borehole under steady-state flow conditions during cross-hole pumping tests in heterogeneous aquifers. These spatial correlation patterns are found to be quite different from that between the hydraulic conductivity variation and the hydraulic head measurement in the same borehole. This finding suggests that a specific discharge measurement carries non-redundant information about the spatial distribution of heterogeneity, even this measurement is collected from the same location where the head measurement is taken. As such, specific discharge observations should be included in the analysis of hydraulic tomography to increase the resolution of estimated aquifer heterogeneity. Using numerical experiments, we demonstrate the effectiveness of the joint interpretation of both hydraulic heads and fluxes for mapping fracture distributions in a hypothetic geologic medium.     

Three-dimensional geostatistical inversion of synthetic tomographic pumping and heat-tracer tests in a nested-cell setup

A main purpose of groundwater inverse modeling lies in estimating the hydraulic conductivity field of an aquifer. Traditionally, hydraulic head measurements, possibly obtained in tomographic setups, are used as data. Because the groundwater flow equation is diffusive, many pumping and observation wells would be necessary to obtain a high resolution of hydraulic conductivity, which is typically not possible. We suggest performing heat tracer tests using the same already installed pumping wells and thermometers in observation planes to amend the hydraulic head data set by the arrival times of the heat signals. For each tomographic combinations of wells, we recommend installing an outer pair of pumping wells, generating artificial ambient flow, and an inner well pair in which the tests are performed. We jointly invert heads and thermal arrival times in 3-D by the quasi-linear geostatistical approach using an efficiently parallelized code running on a mid-range cluster. In the present study, we evaluate the value of heat tracer versus head data in a synthetic test case, where the estimated fields can be compared to the synthetic truth. Because the sensitivity patterns of the thermal arrival times differ from those of head measurements, the resolved variance in the estimated field is 6 to 10 times higher in the joint inversion in comparison to inverting head data only. Also, in contrast to head measurements, reversing the flow field and repeating the heat-tracer test improves the estimate in terms of reducing the estimation variance of the estimate. Based on the synthetic test case, we recommend performing the tests in four principal directions, requiring in total eight pumping wells and four intersecting observation planes for heads and temperature in each direction.     

Two-dimensional physical-based inversion of confined and unconfined aquifers under unknown boundary conditions

An inverse method is developed to simultaneously estimate multiple hydraulic conductivities, source/sink strengths, and boundary conditions, for two-dimensional confined and unconfined aquifers under non-pumping or pumping conditions. The method incorporates noisy observed data (hydraulic heads, groundwater fluxes, or well rates) at measurement locations. With a set of hybrid formulations, given sufficient measurement data, the method yields well-posed systems of equations that can be solved efficiently via nonlinear optimization. The solution is stable when measurement errors are increased. The method is successfully tested on problems with regular and irregular geometries, different heterogeneity patterns and variances (maximum K_max/K_min tested is 10,000), and error magnitudes. Under non-pumping conditions, when error-free observed data are used, the estimated conductivities and recharge rates are accurate within 8% of the true values. When data contain increasing errors, the estimated parameters become less accurate, as expected. For problems where the underlying parameter variation is unknown, equivalent conductivities and average recharge rates can be estimated. Under pumping (and/or injection) conditions, a hybrid formulation is developed to address these local source/sink effects, while different types of boundary conditions can also exert significant influences on drawdowns. Local grid refinement near wells is not needed to obtain accurate results, thus inversion is successful with coarse inverse grids, leading to high computation efficiency. Furthermore, flux measurements are not needed for the inversion to succeed; data requirement of the method is thus not much different from that of interpreting classic well tests. Finally, inversion accuracy is not sensitive to the degree of nonlinearity of the flow equations. Performance of the inverse method for confined and unconfined aquifer problems is similar in terms of the accuracy of the estimated parameters, the recovered head fields, and the solver speed.     

Three-dimensional deformation and strain induced by municipal pumping, Part 2: Numerical analysis

Field measurements consisting of water levels from a municipal well and three-dimensional surface deformations and strains from high-precision GPS measurements at various radial distances from the well were collected as part of a 62-day controlled aquifer test at Mesquite, NV. These measurements were used as observations in several numerical models and a parameter estimation code to characterize and constrain hydraulic and mechanical properties of a 400 m thick basin-fill aquifer. A parsimonious approach was used in conceptualizing the aquifer system. Nonetheless, results from the calibrated deformation and flow models accurately reproduced the observed head and deformations during the first 20 days of pumping, the time at which a new equilibrium was achieved. Surface deformations were shown to reflect hydraulic anisotropy and direction of principal conductivity. In addition, the radius of influence and cone of depression from pumping was approximated in spite of the fact that no monitoring well data existed at the site. Sensitivity analysis indicates that cyclical head values are most sensitive to changes in horizontal hydraulic conductivity, while time-dependent vertical deformations are most sensitive to changes in skeletal specific storage. This investigation shows that GPS monitoring can be used in place of costly monitoring wells to characterize aquifers for water-management purposes where skeletal deformation tends to be elastic.     

Interpreting Variations in Groundwater Flows from Repeated Distributed Thermal Perturbation Tests

To better understand the groundwater resources of southern Nye County, Nevada, a multipart distributed thermal perturbation sensing (DTPS) test was performed on a complex of three wells. These wells penetrate an alluvial aquifer that drains the Nevada National Security Site, and characterizing the hydraulic properties and flow paths of the regional groundwater flow system has proven very difficult. The well complex comprised one pumping well and two observation wells, both located 18 m from the pumping well. Using fiber‐optic cables and line heaters, DTPS tests were performed under both stressed and unstressed conditions. Each test injects heat into the water column over a period of one to two days, and observes the rising temperature during heat injection and falling temperatures after heating ceases. Aquifer thermal properties are inferred from temperature patterns in the cased section of the wells, and fluxes through the 30‐m screened section are estimated based on a model that incorporates conductive and advective heat fluxes. Vertical variations in flux are examined on a scale of tens of cm. The actively flowing zones of the aquifer change between the stressed and unstressed test, and anisotropy in the aquifer permeability is apparent from the changing fluxes between tests. The fluxes inferred from the DTPS tests are compared to solute tracer tests previously performed on the same site. The DTPS‐based fluxes are consistent with the fastest solute transport observed in the tracer test, but appear to overestimate the mean flux through the system.     

Necessary conditions for inverse modeling of flow through variably saturated porous media

Non-unique solutions of inverse problems arise from a lack of information that satisfies necessary conditions for the problem to be well defined. This paper investigates these conditions for inverse modeling of water flow through multi-dimensional variably saturated porous media. It shows that in order to obtain a unique estimate of hydraulic parameters, along each streamline of the flow field (1) spatial and temporal head observations must be given; (2) the number of spatial and temporal head observations required should be greater or equal to the number of unknown parameters; (3) the flux boundary condition or the pumping rate of a well must be specified for the homogeneous case and both boundary flux and pumping rate are a must for the heterogeneous case; (4) head observations must encompass both saturated and unsaturated conditions, and the functional relationships for unsaturated hydraulic conductivity/pressure head and for the moisture retention should be given, and (5) the residual water content value also need to be specified a priori or water content measurements are needed for the estimation of the saturated water content.For field problems, these necessary conditions can be collected or estimated but likely involve uncertainty. While the problems become well defined and have unique solutions, the solutions likely will be uncertain. Because of this uncertainty, stochastic approaches are deemed to be appropriate for inverse problems as they are for forward problems to address uncertainty. Nevertheless, knowledge of these necessary conditions is critical to reduce uncertainty in both characterization of the vadose zone and the aquifer, and prediction of water flow and solute migration in the subsurface.     

Combining 3D Hydraulic Tomography with Tracer Tests for Improved Transport Characterization

Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.     

Cross-borehole flowmeter tests for transient heads in heterogeneous aquifers

Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.     

A Potential-Based Inversion of Unconfined Steady-State Hydraulic Tomography

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.     

Effects of experimental uncertainty on the calculation of hillslope flow paths

Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Assessing Lakebed Hydraulic Conductivity and Seepage Flux by Potentiomanometer

Potentiomanometers (PMs) are commonly used to determine flux directions across interfaces between surface waters and aquifers. We describe a complementary function: estimating small‐scale hydraulic conductivity (K) in a lakebed, using the constant‐head injection test (CHIT) by Cardenas and Zlotnik (2003) with the PM designed by Winter et al. (1988). A piezometer with a small screen is inserted into the lakebed. Local head potential is obtained by measuring the head difference between the test point and the aquifer interface. The piezometer is then used for water injection. This technique is illustrated by measurements taken from Alkali Lake in the Sand Hills, Nebraska, United States. Lakebed K and seepage fluxes ranged from 0.037 to 0.090 m/d and Darcy velocities ranged from 0.004 to 0.027 m/d. Results were consistent with the supplementary data gathered using a modified CHIT and a cone penetrometer. The compact size of the device and the small volumes used for injection enable this method to estimate lakebed K values as low as 0.01 to 0.1 m/d, a range seldom explored in lake‐aquifer interface systems.     

Steady flow to a partially penetrating blind‐wall well in a confined aquifer

Wells in aquifers of loose collapsible sediment are cased so that they have a blind wall and gain water only from the bottom. The hydraulic gradient established at the bottom of these wells during pumping brings the aquifer materials in a quicksand state, which may cause abrasion of pipes and pumps and even the destruction of well structure. To examine the quicksand occurrence, an analytical solution for the steady flow to a partially penetrating blind‐wall well in a confined aquifer is developed. The validity of the proposed solution is evaluated numerically. The sensitivity of maximum vertical gradient along the well bottom in response to aquifer and well parameters is examined. The solution is presented in the form of dimensionless‐type curves and equations that can be easily used to design the safe pumping rate and optimum well geometry to protect the well against sand production. The solution incorporates the anisotropy of aquifer materials and can also be used to determine the hydraulic conductivity of the aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

Time-related capture zones for radial flow in two dimensional randomly heterogeneous media

We consider the effect of randomly heterogeneous hydraulic conductivity on the spatial location of time-related capture zones (isochrones) for a non-reactive tracer in the steady-state radial flow field due to a pumping well in a confined aquifer. A Monte Carlo (MC) procedure is used in conjunction with FFT-based spectral methods. The log hydraulic conductivity field is assumed to be Gaussian and stationary, with isotropic exponential correlation. Various degrees of domain heterogeneity are considered and stability and accuracy of the MC procedure is examined. The location of an isochrone becomes uncertain due to heterogeneity, and it is strongly influenced by hydraulic conductivity variance. The probability that a particle released at a point in the aquifer is pumped by the well within a given time is identified. We propose a new expression for the probabilistic spatial distribution of isochrones, which is formally similar to the analytical solution for a uniform medium and takes into account the effects of heterogeneity.     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Modeling Steady Sea Water Intrusion with Single‐Density Groundwater Codes

Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m³). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.     

Small‐Scale ASR Between Flow Barriers in a Saline Aquifer

Regular aquifer storage recovery, ASR, is often not feasible for small‐scale storage in brackish or saline aquifers because fresh water floats to the top of the aquifer where it is unrecoverable. Flow barriers that partially penetrate a brackish or saline aquifer prevent a stored volume of fresh water from expanding sideways, thus increasing the recovery efficiency. In this paper, the groundwater flow and mixing is studied during injection, storage, and recovery of fresh water in a brackish or saline aquifer in a flow‐tank experiment and by numerical modeling to investigate the effect of density difference, hydraulic conductivity, pumping rate, cyclic operation, and flow barrier settings. Two injection and recovery methods are investigated: constant flux and constant head. Fresh water recovery rates on the order of 65% in the first cycle climbing to as much as 90% in the following cycles were achievable for the studied configurations with constant flux whereas the recovery efficiency was somewhat lower for constant head. The spatial variation in flow velocity over the width of the storage zone influences the recovery efficiency, because it induces leakage of fresh water underneath the barriers during injection and upconing of salt water during recovery.     

Simple Hydraulic Conductivity Estimation by the Kalman Filtered Double Constraint Method

This paper presents the Kalman Filtered Double Constraint Method (DCM‐KF) as a technique to estimate the hydraulic conductivities in the grid blocks of a groundwater flow model. The DCM is based on two forward runs with the same initial grid block conductivities, but with alternating flux‐head conditions specified on parts of the boundary and the wells. These two runs are defined as: (1) the flux run, with specified fluxes (recharge and well abstractions), and (2) the head run, with specified heads (measured in piezometers). Conductivities are then estimated as the initial conductivities multiplied by the fluxes obtained from the flux run and divided by the fluxes obtained from the head run. The DCM is easy to implement in combination with existing models (e.g., MODFLOW). Sufficiently accurate conductivities are obtained after a few iterations. Because of errors in the specified head‐flux couples, repeated estimation under varying hydrological conditions results in different conductivities. A time‐independent estimate of the conductivities and their inaccuracy can be obtained by a simple linear KF with modest computational requirements. For the Kleine Nete catchment, Belgium, the DCM‐KF yields sufficiently accurate calibrated conductivities. The method also results in distinguishing regions where the head‐flux observations influence the calibration from areas where it is not able to influence the hydraulic conductivity.