Comparison of approaches for predicting solute transport: sandbox experiments

The main purpose of this paper was to compare three approaches for predicting solute transport. The approaches include: (1) an effective parameter/macrodispersion approach (Gelhar and Axness 1983); (2) a heterogeneous approach using ordinary kriging based on core samples; and (3) a heterogeneous approach based on hydraulic tomography. We conducted our comparison in a heterogeneous sandbox aquifer. The aquifer was first characterized by taking 48 core samples to obtain local-scale hydraulic conductivity (K). The spatial statistics of these K values were then used to calculate the effective parameters. These K values and their statistics were also used for kriging to obtain a heterogeneous K field. In parallel, we performed a hydraulic tomography survey using hydraulic tests conducted in a dipole fashion with the drawdown data analyzed using the sequential successive linear estimator code (Yeh and Liu 2000) to obtain a K distribution (or K tomogram). The effective parameters and the heterogeneous K fields from kriging and hydraulic tomography were used in forward simulations of a dipole conservative tracer test. The simulated and observed breakthrough curves and their temporal moments were compared. Results show an improvement in predictions of drawdown behavior and tracer transport when the K tomogram from hydraulic tomography was used. This suggests that the high-resolution prediction of solute transport is possible without collecting a large number of small-scale samples to estimate flow and transport properties that are costly to obtain at the field scale.     

Cost Comparisons of Aquifer Heterogeneity Characterization Methods

The characterization of heterogeneity in hydraulic conductivity (K) is a major challenge for subsurface remediation projects. There are a number of field studies that compare the K estimates obtained using various techniques, but to our knowledge, no field‐based studies exists that compare the performance of estimated K heterogeneity fields or the associated characterization costs. In this paper, we compare the costs of characterizing the three‐dimensional K heterogeneity and its uncertainty estimates of a glaciofluvial aquifer‐aquitard sequence at a 15 m × 15 m × 18 m field site situated on the University of Waterloo campus. We compare geostatistical analysis of high resolution permeameter K data obtained from repacked core samples in five boreholes and hydraulic tomography analysis of four pumping tests consisting of up to 41 monitoring points per test. Aside from the comparison of costs, we also assess the performance of each method by predicting several pumping tests. Our analysis reveals that hydraulic tomography is somewhat more costly than the geostatistical analysis of high resolution permeameter K data due to the higher capital costs associated with the method. However, the equipment may be reused at other sites; hence these costs may be recovered over the life of the equipment. More significantly, hydraulic tomography is able to capture the most important features of the aquifer‐aquitard sequence leading to more accurate predictions of independent pumping tests. This suggests that more robust remediation systems may be designed if site characterization is performed with hydraulic tomography.     

Transient flowmeter test: semi-analytic crossflow model

The transient flowmeter test (TFMT) provides more information about the well–aquifer system than the traditional quasi-steady-state flowmeter test (QFMT). The TFMT duration may be much shorter than that of a QFMT, which is desirable at highly contaminated sites where the extracted water has to be treated as hazardous waste. Here we present the TFMT model that accounts for inter-layer crossflow, a thick skin surrounding the well, and wellbore storage. The model is derived under the simplifying assumptions of the pseudo-steady-state inter-layer crossflow and the uniform wellface flux within each layer. The semi-analytic solution is inverted numerically from the Laplace domain to the time domain. Layer and skin parameters are estimated from the TFMT data via the modified Levenberg–Marquardt algorithm. The estimation is robust when the initial parameter guesses are close to their true values. Otherwise, a computationally expensive search among the local minima of the objective function is necessary to find the parameter estimates. The modeling errors and the associated parameter estimation errors are evaluated in a number of synthetic TFMTs and compared to the corresponding results obtained with a general numerical model that relaxes the two simplifying assumptions. The TFMT provides reasonably accurate estimates of hydraulic conductivities for the aquifer layers and the damaged skins and order-of-magnitude estimates of layer specific storativities and hydraulic conductivities for the normal skin. The skin specific storativities should not be estimated from a TFMT. Multi-rate TFMTs with a step-variable pumping rate yield significantly more accurate parameters than constant-pumping-rate TFMTs. The calculated modeling errors may be useful in estimating the magnitude of parameter estimation errors from the TFMT. Our field tests in a coastal aquifer at the Lizzie Site in North Carolina (USA) demonstrate the feasibility of a TFMT for aquifer characterization. The downhole hydraulic conductivity profiles from our field and synthetic TFMTs are consistent with the corresponding profiles from QFMTs.     

Analysis of an unconfined aquifer subject to asynchronous dual-tide propagation

Most published solutions for aquifer responses to ocean tides focus on the one-sided attenuation of the signal as it propagates inland. However, island aquifers experience periodic forcing from the entire coast, which can lead to integrated effects of different tidal signals, especially on narrow high-permeability islands. In general, studies disregard a potential time lag as the tidal wave sweeps around the island. We present a one-dimensional analytical solution to the ground water flow equation subject to asynchronous and asymmetric oscillating head conditions on opposite boundaries and test it on data from an unconfined volcanic aquifer in Maui. The solution considers sediment-damping effects at the coastline. The response of Maui Aquifers indicate that water table elevations near the center of the aquifer are influenced by a combination of tides from opposite coasts. A better match between the observed ground water head and the theoretical response can be obtained with the proposed dual-tide solution than with single-sided solutions. Hydraulic diffusivity was estimated to be 2.3 × 10⁷ m²/d. This translates into a hydraulic conductivity of 500 m/d, assuming a specific yield of 0.04 and an aquifer thickness of 1.8 km. A numerical experiment confirmed the hydraulic diffusivity value and showed that the y -intercepts of the modal attenuation and phase differences estimated by regression can approximate damping factors caused by low-permeability units at the boundary.     

深井水位对地震波、固体潮和气压的水力响应——以范县井为例

以河南范县井为例,利用不同的水力响应模型分析了井水位对地震波、固体潮和气压的响应特征,并基于相关水力响应模型反演估算了含水层的水力参数。结果显示:在高频加载作用过程中,井-含水层系统中的水流模式以水平向为主,而在低频加载作用过程中,则为水平向和垂直向共存的混合模式;利用周期为10—102 s的高频段的地震波响应模型估算的含水层导水系数值较大,为7.20×10-3 m2/s,利用周期为3.75×104 s的低频段的固体潮响应模型估算的含水层导水系数值较小,为2.02×10-6 m2/s,而利用周期为102—104 s的中等频率段的气压响应模型得到的估算值介于二者之间,为3.44×10-5 m2/s。由此分析认为,在周期性加载作用过程中,井-含水层系统内的水流模式与加载频率有关,基于不同水力响应模型反演估算的含水层水力参数存在尺度效应。本研究取得的认识,既可为井水位动态响应的机理解释提供理论基础,也可为目标含水层水力参数的原位测量提供技术支撑。     

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.     

Redundant and Nonredundant Information for Model Calibration or Hydraulic Tomography

Drawdown data from independent pumping tests have widely been used to validate the estimated hydraulic parameters from inverse modeling or hydraulic tomography (HT). Yet, the independent pumping test has not been clearly defined. Therefore, the goal of this paper is to define this independent pumping test concept, based on the redundant or nonredundant information about aquifer heterogeneity embedded in the observed heads during cross-hole pumping tests. The definition of complete, moderate redundancy and high nonredundancy of information are stipulated using cross-correlation analysis of the relationship between the head and heterogeneity. Afterward, data from numerical experiments and field sequential pumping test campaigns reinforce the concept and the definition.     

Aquifer Characterization and Uncertainty in Multi-Frequency Oscillatory Flow Tests: Approach and Insights

Characterizing aquifer properties and their associated uncertainty remains a fundamental challenge in hydrogeology. Recent studies demonstrate the use of oscillatory flow interference testing to characterize effective aquifer flow properties. These characterization efforts relate the relative amplitude and phase of an observation signal with a single frequency component to aquifer diffusivity and transmissivity. Here, we present a generalized workflow that relates extracted Fourier coefficients for observation signals with single and multiple frequency components to aquifer flow properties and their associated uncertainty. Through synthetic analytical modeling we show that multi-frequency oscillatory flow interference testing adds information that improves inversion performance and decreases parameter uncertainty. We show increased observation signal length, sampling frequency, and pressure sensor accuracy all produce decreased parameter uncertainty. This work represents the first attempt we are aware of to quantify effective aquifer parameters and their associated uncertainty using multi-frequency oscillatory flow interference testing.     

Practical issues in imaging hydraulic conductivity through hydraulic tomography

Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (S(s)). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and S(s) using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define "K tomogram" as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.     

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.     

Application of ERT,Saline Tracer and Numerical Studies to Delineate Preferential Paths in Fractured Granites

Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.     

Image-guided inversion in steady-state hydraulic tomography

In steady-state hydraulic tomography, the head data recorded during a series of pumping or/and injection tests can be inverted to determine the transmissivity distributions of an aquifer. This inverse problem is usually under-determined and ill-posed. We propose to use structural information inferred from a guiding image to constrain the inversion process. The guiding image can be drawn from soft data sets such as seismic and ground penetrating radar sections or from geological cross-sections inferred from the wells and some geological expertise. The structural information is extracted from the guiding image through some digital image analysis techniques. Then, it is introduced into the inversion process of the head data as a weighted four direction smoothing matrix used in the regularizer. Such smoothing matrix allows applying the smoothing along the structural features. This helps preserving eventual drops in the hydraulic properties. In addition, we apply a procedure called image-guided interpolation. This technique starts with the tomogram obtained from the image-guided inversion and focus this tomogram. These new approaches are applied on four synthetic toy problems. The hydraulic distributions estimated from the image-guided inversion are closer to the true transmissivity model and have higher resolution than those computed from a classical Gauss–Newton method with uniform isotropic smoothing.     

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.     

Basin-scale transmissivity and storativity estimation using hydraulic tomography

While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion.     

Estimation of spatial variability of aquifer parameters from geophysical methods: a case study of Sindhudurg district,Maharashtra, India

Accurate estimation of aquifer parameters, especially from crystalline hard rock area, assumes a special significance for management of groundwater resources. The aquifer parameters are usually estimated through pumping tests carried out on water wells. While it may be costly and time consuming for carrying out pumping tests at a number of sites, the application of geophysical methods in combination with hydro-geochemical information proves to be potential and cost effective to estimate aquifer parameters. Here a method to estimate aquifer parameters such as hydraulic conductivity, formation factor, porosity and transmissivity is presented by utilizing electrical conductivity values analysed via hydro-geochemical analysis of existing wells and the respective vertical electrical sounding (VES) points of Sindhudurg district, western Maharashtra, India. Further, prior to interpolating the distribution of aquifer parameters of the study area, variogram modelling was carried out using data driven techniques of kriging, automatic relevance determination based Bayesian neural networks (ARD-BNN) and adaptive neuro-fuzzy neural networks (ANFIS). In total, four variogram model fitting techniques such as spherical, exponential, ARD-BNN and ANFIS were compared. According to the obtained results, the spherical variogram model in interpolating transmissivity, ARD-BNN variogram model in interpolating porosity, exponential variogram model in interpolating aquifer thickness and ANFIS variogram model in interpolating hydraulic conductivity outperformed rest of the variogram models. Accordingly, the accurate aquifer parameters maps of the study area were produced by using the best variogram model. The present results suggest that there are relatively high value of hydraulic conductivity, porosity and transmissivity at Parule, Mogarne, Kudal, and Zarap, which would be useful to characterize the aquifer system over western Maharashtra.     

Matching solute breakthrough with deterministic and stochastic aquifer models

Two different deterministic and two alternative stochastic (i.e., geostatistical) approaches to modeling the distribution of hydraulic conductivity (K) in a nonuniform (sigma2ln(K)) = 0.29) glacial sand aquifer were used to explore the influence of conceptual model selection on simulations of three-dimensional tracer movement. The deterministic K models employed included a homogeneous effective K and a perfectly stratified 14 layer model. Stochastic K models were constructed using sequential Gaussian simulation and sequential i ndicator simulation conditioned to available K values estimated from measured grain size distributions. Standard simulation software packages MODFLOW, MT3DMS, and MODPATH were used to model three-dimensional ground water flow and transport in a field tracer test, where a pulse of bromide was injected through an array of three fully screened wells and extracted through a single fully screened well approximately 8 m away. Agreement between observed and simulated transport behavior was assessed through direct comparison of breakthrough curves (BTCs) and selected breakthrough metrics at the extraction well and at 26 individual multilevel sample ports distributed irregularly between the injection and extraction wells. Results indicate that conceptual models incorporating formation variability are better able to capture observed breakthrough behavior. Root mean square (RMS) error of the deterministic models bracketed the ensemble mean RMS error of stochastic models for simulated concentration vs. time series, but not for individual BTC characteristic metrics. The spatial variability models evaluated here may be better suited to simulating breakthrough behavior measured in wells screened over large intervals than at arbitrarily distributed observation points within a nonuniform aquifer domain.     

Joint Estimation of Hydraulic and Poroelastic Parameters from a Pumping Test

The coupling of hydraulic and poroelastic processes is critical in predicting processes involving the deformation of the geologic medium in response to fluid extraction or injection. Numerical models that consider the coupling of hydraulic and poroelastic processes require the knowledge of relevant parameters for both aquifer and aquitard units. In this study, we jointly estimated hydraulic and poroelastic parameters from pumping test data exhibiting “reverse water level fluctuations,” known as the Noordbergum effect, in aquitards adjacent to a pumped aquifer. The joint estimation was performed by coupling BIOT2, a finite element, two‐dimensional, axisymmetric, groundwater model that considers poroelastic effects with the parameter estimation code PEST. We first tested our approach using a synthetic data set with known parameters. Results of the synthetic case showed that for a simple layered system, it was possible to reproduce accurately both the hydraulic and poroelastic properties for each layer. We next applied the approach to pumping test data collected at the North Campus Research Site (NCRS) on the University of Waterloo (UW) campus. Based on the detailed knowledge of stratigraphy, a five‐layer system was modeled. Parameter estimation was performed by: (1) matching drawdown data individually from each observation port and (2) matching drawdown data from all ports at a single well simultaneously. The estimated hydraulic parameters were compared to those obtained by other means at the site yielding good agreement. However, the estimated shear modulus was higher than the static shear modulus, but was within the range of dynamic shear modulus reported in the literature, potentially suggesting a loading rate effect.     

Identifying connections in a fractured rock aquifer using ADFTs

Fractured rock aquifers are difficult to characterize because of their extremely heterogeneous nature. Developing an understanding of fracture network hydraulic properties in these aquifers is difficult and time consuming, and field testing techniques for determining the location and connectivity of fractures in these aquifers are limited. In the Clare Valley, South Australia, well interference is an important issue for a major viticultural area that uses a fractured aquifer. Five fracture sets exist in the aquifer, all dipping > 25 degrees . In this setting, we evaluate the ability of steady-state asymmetric dipole-flow tests (ADFTs) to determine the connections between a test well and a set of piezometers. The procedure involves dividing a test well into two chambers using a single packer and pumping fluid from the upper chamber to the lower chamber. By conducting a series of tests at different packer elevations, an "input" signal is generated in fracture zones connected to the test well. By monitoring the "output" response of the hydraulic dipole field at piezometers, the connectivity of the fractures between the test well and piezometers can be determined. Results indicate the test well used in this study is connected in a complex three-dimensional geometry, with drawdown occurring above and below areas of potentiometric buildup. The ADFT method demonstrates that the aquifer evaluated in this study cannot be modeled effectively on the well scale using continuum flow models.     

Inherent Limitations of Hydraulic Tomography

We offer a cautionary note in response to an increasing level of enthusiasm regarding high-resolution aquifer characterization with hydraulic tomography. We use synthetic examples based on two recent field experiments to demonstrate that a high degree of nonuniqueness remains in estimates of hydraulic parameter fields even when those estimates are based on simultaneous analysis of a number of carefully controlled hydraulic tests. We must, therefore, be careful not to oversell the technique to the community of practicing hydrogeologists, promising a degree of accuracy and resolution that, in many settings, will remain unattainable, regardless of the amount of effort invested in the field investigation. No practically feasible amount of hydraulic tomography data will ever remove the need to regularize or bias the inverse problem in some fashion in order to obtain a unique solution. Thus, along with improving the resolution of hydraulic tomography techniques, we must also strive to couple those techniques with procedures for experimental design and uncertainty assessment and with other more cost-effective field methods, such as geophysical surveying and, in unconsolidated formations, direct-push profiling, in order to develop methods for subsurface characterization with the resolution and accuracy needed for practical field applications.