Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x     

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 (δ³⁷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.     

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.     

Lessons learned from a suite of CFB Borden experiments

This article summarizes several of many field-based studies of subsurface contaminant transport conducted over the last 30 years at the Canadian Forces Base (CFB) Borden site. The field research initially consisted of extensive monitoring of a leachate plume from an abandoned landfill and its analytical and numerical modeling. Lessons learned from these initial studies led to the execution and interpretation of a variety of tracer tests involving conservative and reactive/organic solutes tests performed at various scales. The lessons learned from these tracer tests revealed a number of deficiencies in classical theories of contaminant dispersion and reaction processes as they occur in groundwater, and thus spawned a new era of process-oriented research within the hydrogeological community. The extensively monitored tracer tests were followed by controlled spills of organic contaminants to observe their subsurface movement and distribution as well as the emplacement of a variety of contaminant sources in the saturated and unsaturated zones to study the ambient transport of contaminants. The controlled spills and emplaced sources of various contaminants were then utilized for testing various active and passive remediation technologies. These studies have led to fundamental insights and lessons learned that have significantly contributed to research on contaminant transport in both the saturated and unsaturated zones. Over the years, data generated by the University of Waterloo (UW) researchers and their collaborators continues to be examined by various research groups and has led to additional new insights on subsurface transport of various chemicals.     

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.     

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 (S_s) 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 S_s heterogeneity for fractured rocks. While hydraulic tomography has recently been shown to be a robust technique, the resolution of the K and S_s 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.     

Estimating hydraulic parameters of a heterogeneous aquitard using long-term multi-extensometer and groundwater level data

The classical aquitard-drainage model COMPAC has been modified to simulate the compaction process of a heterogeneous aquitard consisting of multiple sub-units (Multi-COMPAC). By coupling Multi-COMPAC with the parameter estimation code PEST++, the vertical hydraulic conductivity (K _v) and elastic (S _ske) and inelastic (S _skp) skeletal specific-storage values of each sub-unit can be estimated using observed long-term multi-extensometer and groundwater level data. The approach was first tested through a synthetic case with known parameters. Results of the synthetic case revealed that it was possible to accurately estimate the three parameters for each sub-unit. Next, the methodology was applied to a field site located in Changzhou city, China. Based on the detailed stratigraphic information and extensometer data, the aquitard of interest was subdivided into three sub-units. Parameters K _v, S _ske and S _skp of each sub-unit were estimated simultaneously and then were compared with laboratory results and with bulk values and geologic data from previous studies, demonstrating the reliability of parameter estimates. Estimated S _skp values ranged within the magnitude of 10^?4 m^?1, while K _v ranged over 10^?10–10^?8 m/s, suggesting moderately high heterogeneity of the aquitard. However, the elastic deformation of the third sub-unit, consisting of soft plastic silty clay, is masked by delayed drainage, and the inverse procedure leads to large uncertainty in the S _ske estimate for this sub-unit.     

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.     

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.     

Hydraulic tomography for detecting fracture zone connectivity

Fracture zones and their connectivity in geologic media are of great importance to ground water resources management as well as ground water contamination prevention and remediation. In this paper, we applied a recently developed hydraulic tomography (HT) technique and an analysis algorithm (sequential successive linear estimator) to synthetic fractured media. The application aims to explore the potential utility of the technique and the algorithm for characterizing fracture zone distribution and their connectivity. Results of this investigation showed that using HT with a limited number of wells, the fracture zone distribution and its connectivity (general pattern) can be mapped satisfactorily although estimated hydraulic property fields are smooth. As the number of wells and monitoring ports increases, the fracture zone distribution and connectivity become vivid and the estimated hydraulic properties approach true values. We hope that the success of this application may promote the development and application of the new generations of technology (i.e., hydraulic, tracer, pneumatic tomographic surveys) for mapping fractures and other features in geologic media.