Assessment of hydraulic and thermal properties of the Antarctic active layer: Insights from laboratory column experiments and inverse modeling.

To better understand the changes in the hydrologic cycle caused by global warming in Antarctica, it is crucial to improve our understanding of the groundwater flow system, which has received less attention despite its significance. Both hydraulic and thermal properties of the active layer, through which groundwater can flow during thawing seasons, are essential to quantify the groundwater flow system. However, there has been insufficient information on the Antarctic active layer. The goal of this study was to estimate the hydraulic and thermal properties of Antarctic soils through laboratory column experiments and inverse modeling. The column experiments were conducted with sediments collected from two lakes in the Barton Peninsula, Antarctica. A sand column was also operated for comparison. Inverse modeling using HydroGeoSphere (HGS) combined with Parameter ESTimation (PEST) was performed with data collected from the column experiments, including permeameter tests, saturation-drain tests, and freeze-thaw tests. Hydraulic parameters (i.e., Ks, θs, Swr, α, ß, and Ss) and thermal diffusivity (D) of the soils were derived from water retention curves and temperature curves with depth, respectively. The hydraulic properties of the Antarctic soil samples, estimated through inverse modeling, were 1.6 × 10-5-3.4 × 10-4 cm s-1 for Ks, 0.37-0.42 for θs, 6.62 × 10-3-1.05 × 10-2 for Swr, 0.53-0.58 cm-1 for α, 5.75-7.96 for ß, and 5.11 × 10-5-9.02 × 10-5 cm-1 for Ss. The thermal diffusivities for the soils were estimated to be 0.65-4.64 cm2 min-1. The soil hydraulic and thermal properties reflected the physical and ecological characteristics of their lake environments. The results of this study can provide a basis for groundwater-surface water interaction in polar regions, which is governed by variably-saturated flow and freeze-thaw processes.

Evaluation of Hydraulic Conductivity Estimates from Various Approaches with Groundwater Flow Models.

Significant efforts have been expended for improved characterization of hydraulic conductivity (K) and specific storage (Ss ) to better understand groundwater flow and contaminant transport processes. Conventional methods including grain size analyses (GSA), permeameter, slug, and pumping tests have been utilized extensively, while Direct Push-based Hydraulic Profiling Tool (HPT) surveys have been developed to obtain high-resolution K estimates. Moreover, inverse modeling approaches based on geology-based zonations, and highly parameterized Hydraulic Tomography (HT) have also been advanced to map spatial variations of K and Ss between and beyond boreholes. While different methods are available, it is unclear which one yields K estimates that are most useful for high resolution predictions of groundwater flow. Therefore, the main objective of this study is to evaluate various K estimates at a highly heterogeneous field site obtained with three categories of characterization techniques including: (1) conventional methods (GSA, permeameter, and slug tests); (2) HPT surveys; and (3) inverse modeling based on geology-based zonations and highly parameterized approaches. The performance of each approach is first qualitatively analyzed by comparing K estimates to site geology. Then, steady-state and transient groundwater flow models are employed to quantitatively assess various K estimates by simulating pumping tests not used for parameter estimation. Results reveal that inverse modeling approaches yield the best drawdown predictions under both steady and transient conditions. In contrast, conventional methods and HPT surveys yield biased predictions. Based on our research, it appears that inverse modeling and data fusion are necessary steps in predicting accurate groundwater flow behavior.

Migration behaviour of LNAPL in fractures filled with porous media: Laboratory experiments and numerical simulations.

With the increasing requirement of maintaining world energy security and strategic reserves, oil storage and transportation facilities are being built at a large scale. Taking the safe and efficient operation of petroleum storage projects as the goal, a set of experimental apparatus to investigate the migration of contaminants in fractures filled with media was developed to predict and evaluate the environmental risk of oil contaminants leakage. A multiphase numerical flow model based on COMSOL was built based on the laboratory experimental model. Specifically, the migration behaviour of Light Non-aqueous Phase Liquid (LNAPL) through a sand-filled fractured medium was studied by laboratory experiments and numerical simulations. Image and chemical analyses methods were used to monitor and study LNAPL migration behaviour for varying grain sizes of porous medium filling the fractures and varying groundwater table elevations. Laboratory experimental results showed that the LNAPL migration velocity in filled fracture network was significantly faster than that in adjoining porous media during the initial stage of infiltration. The migration velocity increased with the relative permeability of filled sand, which was closely related to the Van Genuchten (VG) model parameters α and n. LNAPL migrated downward with the falling groundwater table and became entrapped with the rising groundwater table, and the amount of entrapment depended on VG model parameters. Hydrogeological parameters were calibrated and LNAPL migration in filled fractured media was predicted using the calibrated numerical model. Simulation results revealed that fracture inclination had an important influence on LNAPL migration in filled fractured media and its migration velocity decreased with a decrease in fracture inclination. These research results can be applied to the control and remediation of oil-contaminated sites in fractured rock settings, such as at underground oil storage tanks and caverns, as well as at underground oil pipelines.

Review of Laboratory Scale Models of Karst Aquifers: Approaches, Similitude, and Requirements.

This review focuses on investigations of groundwater flow and solute transport in karst aquifers through laboratory scale models (LSMs). In particular, LSMs have been used to generate new data under different hydraulic and contaminant transport conditions, testing of new approaches for site characterization, and providing new insights into flow and transport processes through complex karst aquifers. Due to the increasing need for LSMs to investigate a wide range of issues, associated with flow and solute migration karst aquifers this review attempts to classify, and introduce a framework for constructing a karst aquifer physical model that is more representative of field conditions. The LSMs are categorized into four groups: sand box, rock block, pipe/fracture network, and pipe-matrix coupling. These groups are compared and their advantages and disadvantages highlighted. The capabilities of such models have been extensively improved by new developments in experimental methods and measurement devices. Newer technologies such as 3D printing, computed tomography scanning, X-rays, nuclear magnetic resonance, novel geophysical techniques, and use of nanomaterials allow for greater flexibilities in conducting experiments. In order for LSMs to be representative of karst aquifers, a few requirements are introduced: (1) the ability to simulate heterogeneous distributions of karst hydraulic parameters, (2) establish Darcian and non-Darcian flow regimes and exchange between the matrix and conduits, (3) placement of adequate sampling points and intervals, and (4) achieving some degree of geometric, kinematic, and dynamic similitude to represent field conditions.

A Type-Curve Method for the Analysis of Pumping Tests with Piecewise-Linear Pumping Rates.

Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K) and specific storage (Ss ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and Ss values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and Ss values between observation wells. The difference could be attributed to the spatial heterogeneity in K and Ss . A future research topic may focus on the characterization of K and Ss heterogeneity with the currently available drawdown data from variable-rate pumping tests.

Use of steady-state hydraulic tomography to inform the selection of a chaotic advection system.

The concept of chaotic advection is a novel approach that has the potential to overcome some of the challenges associated with mixing of reagents that commonly occur when injection based in situ treatment techniques are used. The rotated potential mixing (RPM) flow system is one configuration which has been theorized to achieve chaotic advection in porous media, and enhance reagent mixing by periodically re-oriented dipole pumping at a series of radial wells. Prior to field implementation of chaotic advection, the selection of an RPM flow protocol will likely require a numerical model that can adequately represent groundwater flow within the zone of interest. As expected, the hydraulic conductivity (K) field is the most critical input requirement for the selected groundwater flow model. Hydraulic tomography (HT) is an innovative characterization approach that has shown potential to provide information on a K field. In this investigation, we explored whether the same well system required to invoke chaotic advection can also be applied in a HT analysis, and evaluated the use of the generated K tomogram for the selection of RPM flow parameters that can enhance reagent mixing. A series of dipole pumping tests were conducted within an area of interest as defined by the limits of the circular network of eight injection/extraction wells used to invoke chaotic advection. Hydraulic head data collected from independent dipole pumping tests were used in an inverse model to perform steady-state hydraulic tomography (SSHT) analysis to generate a K tomogram. Both the K tomogram and an effective parameter approach (i.e., a single K value assigned across the entire spatial domain as determined by single well pumping and slug tests) produced estimates of hydraulic head that closely resembled those observed due to the relative homogeneous nature of the aquifer and the small spatial scale of the area of interest. In contrast, particle tracking results showed that incorporating a heterogeneous K field significantly enhanced the spatial distribution of particle trajectories indicative of reagent mixing. These findings support the hypothesis that the same well system used to invoke chaotic advection can be combined with SSHT analysis as a viable site characterization tool for delineating the spatial variability of K. Incorporating this K tomogram in a groundwater flow model with a particle tracking engine can be used as a design tool to aid in the selection of a site-specific RPM flow protocol to achieve enhanced reagent mixing.

A Joint Analytic Method for Estimating Aquitard Hydraulic Parameters.

The vertical hydraulic conductivity (Kv ), elastic (Sske ), and inelastic (Sskv ) skeletal specific storage of aquitards are three of the most critical parameters in land subsidence investigations. Two new analytic methods are proposed to estimate the three parameters. The first analytic method is based on a new concept of delay time ratio for estimating Kv and Sske of an aquitard subject to long-term stable, cyclic hydraulic head changes at boundaries. The second analytic method estimates the Sskv of the aquitard subject to linearly declining hydraulic heads at boundaries. Both methods are based on analytical solutions for flow within the aquitard, and they are jointly employed to obtain the three parameter estimates. This joint analytic method is applied to estimate the Kv , Sske , and Sskv of a 34.54-m thick aquitard for which the deformation progress has been recorded by an extensometer located in Shanghai, China. The estimated results are then calibrated by PEST (Doherty 2005), a parameter estimation code coupled with a one-dimensional aquitard-drainage model. The Kv and Sske estimated by the joint analytic method are quite close to those estimated via inverse modeling and performed much better in simulating elastic deformation than the estimates obtained from the stress-strain diagram method of Ye and Xue (2005). The newly proposed joint analytic method is an effective tool that provides reasonable initial values for calibrating land subsidence models.

Compound-Specific Isotope Analyses to Assess TCE Biodegradation in a Fractured Dolomitic Aquifer.

The potential for trichloroethene (TCE) biodegradation in a fractured dolomite aquifer at a former chemical disposal site in Smithville, Ontario, Canada, is assessed using chemical analysis and TCE and cis-DCE compound-specific isotope analysis of carbon and chlorine collected over a 16-month period. Groundwater redox conditions change from suboxic to much more reducing environments within and around the plume, indicating that oxidation of organic contaminants and degradation products is occurring at the study site. TCE and cis-DCE were observed in 13 of 14 wells sampled. VC, ethene, and/or ethane were also observed in ten wells, indicating that partial/full dechlorination has occurred. Chlorine isotopic values (δ37 Cl) range between 1.39 to 4.69‰ SMOC for TCE, and 3.57 to 13.86‰ SMOC for cis-DCE. Carbon isotopic values range between -28.9 and -20.7‰ VPDB for TCE, and -26.5 and -11.8‰ VPDB for cis-DCE. In most wells, isotopic values remained steady over the 15-month study. Isotopic enrichment from TCE to cis-DCE varied between 0 and 13‰ for carbon and 1 and 4‰ for chlorine. Calculated chlorine-carbon isotopic enrichment ratios (ϵCl /ϵC ) were 0.18 for TCE and 0.69 for cis-DCE. Combined, isotopic and chemical data indicate very little dechlorination is occurring near the source zone, but suggest bacterially mediated degradation is occurring closer to the edges of the plume.

An Application of Hydraulic Tomography to a Large-Scale Fractured Granite Site, Mizunami, Japan.

While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer-scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non-redundant datasets about the fracture and fault zones can be collected. Only with these non-redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large-scale fracture and fault zone distributions.

Determination of rate constants and branching ratios for TCE degradation by zero-valent iron using a chain decay multispecies model.

The applicability of a newly-developed chain-decay multispecies model (CMM) was validated by obtaining kinetic rate constants and branching ratios along the reaction pathways of trichloroethene (TCE) reduction by zero-valent iron (ZVI) from column experiments. Changes in rate constants and branching ratios for individual reactions for degradation products over time for two columns under different geochemical conditions were examined to provide ranges of those parameters expected over the long-term. As compared to the column receiving deionized water, the column receiving dissolved CaCO3 showed higher mean degradation rates for TCE and all of its degradation products. However, the column experienced faster reactivity loss toward TCE degradation due to precipitation of secondary carbonate minerals, as indicated by a higher value for the ratio of maximum to minimum TCE degradation rate observed over time. From the calculated branching ratios, it was found that TCE and cis-dichloroethene (cis-DCE) were dominantly dechlorinated to chloroacetylene and acetylene, respectively, through reductive elimination for both columns. The CMM model, validated by the column test data in this study, provides a convenient tool to determine simultaneously the critical design parameters for permeable reactive barriers and natural attenuation such as rate constants and branching ratios.

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.

Comparison of hydraulic tomography with traditional methods at a highly heterogeneous site.

Over the past several decades, different groundwater modeling approaches of various complexities and data use have been developed. A recently developed approach for mapping hydraulic conductivity (K) and specific storage (Ss ) heterogeneity is hydraulic tomography, the performance of which has not been compared to other more "traditional" methods that have been utilized over the past several decades. In this study, we compare seven methods of modeling heterogeneity which are (1) kriging, (2) effective parameter models, (3) transition probability/Markov Chain geostatistics models, (4) geological models, (5) stochastic inverse models conditioned to local K data, (6) hydraulic tomography, and (7) hydraulic tomography conditioned to local K data using data collected in five boreholes at a field site on the University of Waterloo (UW) campus, in Waterloo, Ontario, Canada. The performance of each heterogeneity model is first assessed during model calibration. In particular, the correspondence between simulated and observed drawdowns is assessed using the mean absolute error norm, (L1 ), mean square error norm (L2 ), and correlation coefficient (R) as well as through scatterplots. We also assess the various models on their ability to predict drawdown data not used in the calibration effort from nine pumping tests. Results reveal that hydraulic tomography is best able to reproduce these tests in terms of the smallest discrepancy and highest correlation between simulated and observed drawdowns. However, conditioning of hydraulic tomography results with permeameter K data caused a slight deterioration in accuracy of drawdown predictions which suggests that data integration may need to be conducted carefully.

Hydraulic tomography offers improved imaging of heterogeneity in fractured rocks.

Fractured rocks have presented formidable challenges for accurately predicting groundwater flow and contaminant transport. This is mainly due to our difficulty in mapping the fracture-rock matrix system, their hydraulic properties and connectivity at resolutions that are meaningful for groundwater modeling. Over the last several decades, considerable effort has gone into creating maps of subsurface heterogeneity in hydraulic conductivity (K) and specific storage (Ss ) of fractured rocks. Developed methods include kriging, stochastic simulation, stochastic inverse modeling, and hydraulic tomography. In this article, I review the evolution of various heterogeneity mapping approaches and contend that hydraulic tomography, a recently developed aquifer characterization technique for unconsolidated deposits, is also a promising approach in yielding robust maps (or tomograms) of K and Ss heterogeneity for fractured rocks. While hydraulic tomography has recently been shown to be a robust technique, the resolution of the K and Ss tomograms mainly depends on the density of pumping and monitoring locations and the quality of data. The resolution will be improved through the development of new devices for higher density monitoring of pressure responses at discrete intervals in boreholes and potentially through the integration of other data from single-hole tests, borehole flowmeter profiling, and tracer tests. Other data from temperature and geophysical surveys as well as geological investigations may improve the accuracy of the maps, but more research is needed. Technological advances will undoubtedly lead to more accurate maps. However, more effort should go into evaluating these maps so that one can gain more confidence in their reliability.

A semi-analytical solution for simulating contaminant transport subject to chain-decay reactions.

We present a set of new, semi-analytical solutions to simulate three-dimensional contaminant transport subject to first-order chain-decay reactions. The aquifer is assumed to be areally semi-infinite, but finite in thickness. The analytical solution can treat the transformation of contaminants into daughter products, leading to decay chains consisting of multiple contaminant species and various reaction pathways. The solution in its current form is capable of accounting for up to seven species and four decay levels. The complex pathways are represented by means of first-order decay and production terms, while branching ratios account for decay stoichiometry. Besides advection, dispersion, bio-chemical or radioactive decay and daughter product formation, the model also accounts for sorption of contaminants on the aquifer solid phase with each species having a different retardation factor. First-type contaminant boundary conditions are utilized at the source (x=0 m) and can be either constant-in-time for each species, or the concentration can be allowed to undergo first-order decay. The solutions are obtained by exponential Fourier, Fourier cosine and Laplace transforms. Limiting forms of the solutions can be obtained in closed form, but we evaluate the general solutions by numerically inverting the analytical solutions in exponential Fourier and Laplace transform spaces. Various cases are generated and the solutions are verified against the HydroGeoSphere numerical model.

Field study of subsurface heterogeneity with steady-state hydraulic tomography.

Remediation of subsurface contamination requires an understanding of the contaminant (history, source location, plume extent and concentration, etc.), and, knowledge of the spatial distribution of hydraulic conductivity (K) that governs groundwater flow and solute transport. Many methods exist for characterizing K heterogeneity, but most if not all methods require the collection of a large number of small-scale data and its interpolation. In this study, we conduct a hydraulic tomography survey at a highly heterogeneous glaciofluvial deposit at the North Campus Research Site (NCRS) located at the University of Waterloo, Waterloo, Ontario, Canada to sequentially interpret four pumping tests using the steady-state form of the Sequential Successive Linear Estimator (SSLE) (Yeh and Liu 2000). The resulting three-dimensional (3D) K distribution (or K-tomogram) is compared against: (1) K distributions obtained through the inverse modeling of individual pumping tests using SSLE, and (2) effective hydraulic conductivity (K(eff) ) estimates obtained by automatically calibrating a groundwater flow model while treating the medium to be homogeneous. Such a K(eff) is often used for designing remediation operations, and thus is used as the basis for comparison with the K-tomogram. Our results clearly show that hydraulic tomography is superior to the inversions of single pumping tests or K(eff) estimates. This is particularly significant for contaminated sites where an accurate representation of the flow field is critical for simulating contaminant transport and injection of chemical and biological agents used for active remediation of contaminant source zones and plumes.

Numerical simulation of DNAPL emissions and remediation in a fractured dolomitic aquifer.

This study presents a numerical model of a large aqueous phase plume of a mixture of chlorinated solvents that has penetrated the fractured dolomitic bedrock near Smithville, Ontario, Canada several decades ago which, since 1989 has been hydraulically controlled by a pump-and-treat remediation system. A multiphase compositional model CompFlow is first applied to simulate the migration of DNAPLs in a discretely fractured porous medium with hydrostratigraphy representing the Smithville site. Results from CompFlow are used to estimate the pure-phase DNAPL distribution in the discrete fractures and rock matrix. Next, CompFlow results are employed to define the source term for a regional-scale transport simulation using HydroGeoSphere (HGS) by treating the layered, fractured dolomitic rocks as an equivalent porous continuum. Transport simulations are conducted both prior to and after the operation of the pump-and-treat system. Results reveal that considerable agreement with the observed mass removal data and TCE plume can be achieved by modifying the composition of the DNAPL source and by reducing the hydraulic conductivity (K) in the source zone region to account for preferential flow around it. Our transport model results support the conceptual model of TCE contamination which posits a mixed source (2 to 4%) of DNAPL with limited contact with actively flowing groundwater that is undergoing equilibrium dissolution. Model results also reveal that the pump-and-treat system has neither been effective in stabilizing the plume nor removing a significant amount of contaminant mass, but that the stability of the plume is instead due to first-order degradation.

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.

Field study of hydrogeologic characterization methods in a heterogeneous aquifer.

Hydraulic conductivity (K) and specific storage (S(s)) are required parameters when designing transient groundwater flow models. The purpose of this study was to evaluate the ability of commonly used hydrogeologic characterization approaches to accurately delineate the distribution of hydraulic properties in a highly heterogeneous glaciofluvial deposit. The metric used to compare the various approaches was the prediction of drawdown responses from three separate pumping tests. The study was conducted at a field site, where a 15 m × 15 m area was instrumented with four 18-m deep Continuous Multichannel Tubing (CMT) wells. Each CMT well contained seven 17 cm × 1.9 cm monitoring ports equally spaced every 2 m down each CMT system. An 18-m deep pumping well with eight separate 1-m long screens spaced every 2 m was also placed in the center of the square pattern. In each of these boreholes, cores were collected and characterized using the Unified Soil Classification System, grain size analysis, and permeameter tests. To date, 471 K estimates have been obtained through permeameter analyses and 270 K estimates from empirical relationships. Geostatistical analysis of the small-scale K data yielded strongly heterogeneous K fields in three-dimensions. Additional K estimates were obtained through slug tests in 28 ports of the four CMT wells. Several pumping tests were conducted using the multiscreen and CMT wells to obtain larger scale estimates of both K and S(s). The various K and S(s) estimates were then quantitatively evaluated by simulating transient drawdown data from three pumping tests using a 3D forward numerical model constructed using HydroGeoSphere (Therrien et al. 2005). Results showed that, while drawdown predictions generally improved as more complexity was introduced into the model, the ability to make accurate drawdown predictions at all CMT ports was inconsistent.

Estimating hydraulic parameters when poroelastic effects are significant.

For almost 80 years, deformation-induced head changes caused by poroelastic effects have been observed during pumping tests in multilayered aquifer-aquitard systems. As water in the aquifer is released from compressive storage during pumping, the aquifer is deformed both in the horizontal and vertical directions. This deformation in the pumped aquifer causes deformation in the adjacent layers, resulting in changes in pore pressure that may produce drawdown curves that differ significantly from those predicted by traditional groundwater theory. Although these deformation-induced head changes have been analyzed in several studies by poroelasticity theory, there are at present no practical guidelines for the interpretation of pumping test data influenced by these effects. To investigate the impact that poroelastic effects during pumping tests have on the estimation of hydraulic parameters, we generate synthetic data for three different aquifer-aquitard settings using a poroelasticity model, and then analyze the synthetic data using type curves and parameter estimation techniques, both of which are based on traditional groundwater theory and do not account for poroelastic effects. Results show that even when poroelastic effects result in significant deformation-induced head changes, it is possible to obtain reasonable estimates of hydraulic parameters using methods based on traditional groundwater theory, as long as pumping is sufficiently long so that deformation-induced effects have largely dissipated.

Hydraulic/partitioning tracer tomography for DNAPL source zone characterization: small-scale sandbox experiments.

Dense nonaqueous phase liquids (DNAPL) are prevalent at a large number of sites throughout the world. The variable release history, unstable flow, and geologic heterogeneity make the spatial distribution of DNAPLs complex. This causes difficulties in site remediation contributing to long-term groundwater contamination for decades to centuries. We present laboratory experiments to demonstrate the efficacy of Sequential Successive Linear Estimator (SSLE) algorithm that images DNAPL source zones. The algorithm relies on the fusion of hydraulic and partitioning tracer tomography (HPTT) to derive the best estimate of the K heterogeneity, DNAPL saturation (S(N)) distribution, and their uncertainty. The approach is nondestructive and can be applied repeatedly. Results from our laboratory experiments show that S(N) distributions compare favorably with DNAPL distributions observed in the sandbox but not so with local saturation estimates from core samples. We also found that the delineation of K heterogeneity can have a large impact on computed S(N) distributions emphasizing the importance of accurate delineation of hydraulic heterogeneity.