DNAPL Characterization Methods and Approaches, Part 2: Cost Comparisons

Contamination from the use of chlorinated solvents, often classified as dense nonaqueous phase liquids (DNAPLs) when in an undissolved state, pose environmental threats to ground water resources worldwide. DNAPL site characterization method performance comparisons are presented in a companion paper (Kram et al. 2001). This study compares the costs for implementing various characterization approaches using synthetic unit model scenarios (UMSs), each with particular physical characteristics. Unit costs and assumptions related to labor, equipment, and consumables are applied to determine costs associated with each approach for various UMSs. In general, the direct-push sensor systems provide cost-effective characterization information in soils that are penetrable with relatively shallow (less than 10 to 15 m) water tables. For sites with impenetrable lithology using direct-push techniques, the Ribbon NAPL Sampler Flexible Liner Underground Technologies Everting (FLUTe) membrane appears to be the most cost-effective approach. For all scenarios studied, partitioning interwell tracer tests (PITTs) are the most expensive approach due to the extensive pre-and post-PITT requirements. However, the PITT is capable of providing useful additional information, such as approximate DNAPL saturation, which is not generally available from any of the other approaches included in this comparison.     

Integrating Hydrogeological and Geophysical Methods for the Characterization of a Deltaic Aquifer System

Groundwater management needs detailed aquifer characterization, especially in semiarid costal aquifer systems that are under hydrological pressure. Our study area is in the Tordera delta, northeastern coast of Spain, where a detrital fluvio-deltaic aquifer system has been developed above granitic basement. The main purpose of this study is to characterize the complex lithological structure and the seawater intrusion state by combining hydrological information, audiomagnetotelluric (AMT) and seismic reflection and refraction models. This allowed us to provide spatially continuous information about aquifer properties and processes. Thus, we have determined the thickness and continuity of the aquifer units, as well as the morphology and depth to the basement. The models revealed that the main seawater intrusion main path is found in the western deltaic area that coincides with an existing buried paleochannel. This new result explains the anomalously high chlorine concentrations observed in the deep semiconfined aquifer more than 1,500 m inland.     

Use of Improved Hydrologic Testing and Borehole Geophysical Logging Methods for Aquifer Characterization

Depth-discrete aquifer in formal ion was obtained using recently developed adaptations and improvements to conventional characterization techniques. These improvements included running neutron porosity and hulk density geophysical logging tools through a cased hole, performing an enhanced point-dilution tracer test for monitoring tracer concentration as a function of Lime and depth, and using pressure derivatives for diagnostic and quantitative analysis of constant rate discharge lest data. Data results from the use of these techniques were used to develop a conceptual model of a heterogeneous aquifer. Depth-discrete aquifer information was required to effectively design field-scale deployment and monitoring of an in situ bioremediation technology.
Geophysical logging and point-dilution tracer test results provided the relative distribution of porosity and horizontal hydraulic conductivity, respectively, with depth and correlated well. Hydraulic pumping tests were conducted to estimate mean values for transmissivity and effective hydraulic conductivity, Tracer lest and geophysical logging results indicated that ground water flow was predominant in the upper approximate 10 feet of the aquifer investigated. These results were used to delineate a more representative interval thickness for estimating effective hydraulic conductivity. Hydraulic conductivity, calculated using this representative interval, was estimated lo be 73 ft/d, approximately three limes higher than that calculated using the full length of the screened test interval.     

DNAPL Characterization Methods and Approaches, Part 1: Performance Comparisons

Contamination from the use of chlorinated solvents, often classified as dense nonaqueous phase liquids (DNAPLs) when in an undissolved state, represents an environmental challenge with global implications. Mass-transfer limitations due to rate-limited dissolution can lead to long-term aquifer persistence for even small volumetric fractions. The identification of DNAPL source zones located beneath the water table is critical to ultimately achieve site remediation and aquifer restoration. This paper provides a comparison of the advantages and disadvantages of many of the methods being used for detecting and delineating DNAPL contaminant source zones. The objective is to determine which options are best to pursue based on site characteristics, method performance, and method costs. DNAPL characterization methods are grouped into approaches, which include site preparation, characterization, and data-processing activities necessary to design an effective remediation system. We compare the different approaches based on the level of chemical and hydrogeologic resolution, and the need for additional data requirements. Our findings can be used to assist with selection of appropriate site remediation management options.     

Characterizing Heterogeneity in Infiltration Rates During Managed Aquifer Recharge

Infiltration rate is the key parameter that describes how water moves from the surface into a groundwater aquifer during managed aquifer recharge (MAR). Characterization of infiltration rate heterogeneity in space and time is valuable information for MAR system operation. In this study, we utilized fiber optic distributed temperature sensing (FO‐DTS) observations and the phase shift of the diurnal temperature signal between two vertically co‐located fiber optic cables to characterize infiltration rate spatially and temporally in a MAR basin. The FO‐DTS measurements revealed spatial heterogeneity of infiltration rate: approximately 78% of the recharge water infiltrated through 50% of the pond bottom on average. We also introduced a metric for quantifying how the infiltration rate in a recharge pond changes over time, which enables FO‐DTS to be used as a method for monitoring MAR and informing maintenance decisions. By monitoring this metric, we found high‐spatial variability in how rapidly infiltration rate changed during the test period. We attributed this variability to biological pore clogging and found a relationship between high initial infiltration rate and the most rapid pore clogging. We found a strong relationship (R² = 0.8) between observed maximum infiltration rates and electrical resistivity measurements from electrical resistivity tomography data taken in the same basin when dry. This result shows that the combined acquisition of DTS and ERT data can improve the design and operation of a MAR pond significantly by providing the critical information needed about spatial variability in parameters controlling infiltration rates.     

Accounting for Aquifer Heterogeneity from Geological Data to Management Tools

A nested workflow of multiple‐point geostatistics (MPG) and sequential Gaussian simulation (SGS) was tested on a study area of 6 km² located about 20 km northwest of Quebec City, Canada. In order to assess its geological and hydrogeological parameter heterogeneity and to provide tools to evaluate uncertainties in aquifer management, direct and indirect field measurements are used as inputs in the geostatistical simulations to reproduce large and small‐scale heterogeneities. To do so, the lithological information is first associated to equivalent hydrogeological facies (hydrofacies) according to hydraulic properties measured at several wells. Then, heterogeneous hydrofacies (HF) realizations are generated using a prior geological model as training image (TI) with the MPG algorithm. The hydraulic conductivity (K) heterogeneity modeling within each HF is finally computed using SGS algorithm. Different K models are integrated in a finite‐element hydrogeological model to calculate multiple transport simulations. Different scenarios exhibit variations in mass transport path and dispersion associated with the large‐ and small‐scale heterogeneity respectively. Three‐dimensional maps showing the probability of overpassing different thresholds are presented as examples of management tools.     

Reducing Uncertainty in Calibrating Aquifer Flow Model with Multiple Scales of Heterogeneity

Modeling and calibration of natural aquifers with multiple scales of heterogeneity is a challenging task due to limited subsurface access. While computer modeling plays an essential role in aquifer studies, large uncertainty exists in developing a conceptual model of an aquifer and in calibrating the model for decision making. Due to uncertainties such as a lack of understanding of subsurface processes and a lack of techniques to parameterize the subsurface environment (including hydraulic conductivity, source/sink rate, and aquifer boundary conditions), existing aquifer models often suffer nonuniqueness in calibration, leading to poor predictive capability. A robust calibration methodology is needed that can address the simultaneous estimations of aquifer parameters, source/sink, and boundary conditions. In this paper, we propose a multistage and multiscale approach that addresses subsurface heterogeneity at multiple scales, while reducing uncertainty in estimating the model parameters and model boundary conditions. The key to this approach lies in the appropriate development, verification, and synthesis of existing and new techniques of static and dynamic data integration. In particular, based on a given set of observation data, new inversion techniques can be first used to estimate aquifer large‐scale effective parameters and smoothed boundary conditions, based on which parameter and boundary condition estimation can be refined at increasing detail using standard or highly parameterized estimation techniques.     

Single Well Tracer Tests in Aquifer Characterization

When the purpose of aquifer testing is to yield data for modeling aqueous mass transport, pumping tests and gradient measurement can only partially satisfy characterization requirements. Effective porosity, ground water flow velocity, and the vertical distribution of hydraulic conductivity within the aquifer are left as unknowns. Single well tracer methods, when added to the testing program, can be used to estimate these parameters. A drift, and pumpback test yields porosity and velocity, and point-dilution testing yields depth-discrete hydraulic information, A single emplacement of tracer into a test well is sufficient to conduct both tests. The tracer tests are facilitated by a simple method for injecting and evenly distributing the tracer solution into a wellbore, and by new ion-selective electrode instrumentation, specifically designed for submersible service, for monitoring the concentration of tracers such as bromide.     

Electrical Resistivity Characterization of Anisotropy in the Biscayne Aquifer

Electrical anisotropy occurs when electric current flow varies with azimuth. In porous media, this may correspond to anisotropy in the hydraulic conductivity resulting from sedimentary fabric, fractures, or dissolution. In this study, a 28‐electrode resistivity imaging system was used to investigate electrical anisotropy at 13 sites in the Biscayne Aquifer of SE Florida using the rotated square array method. The measured coefficient of electrical anisotropy generally ranged from 1.01 to 1.12 with values as high as 1.36 found at one site. The observed electrical anisotropy was used to estimate hydraulic anisotropy (ratio of maximum to minimum hydraulic conductivity) which ranged from 1.18 to 2.83. The largest values generally were located on the Atlantic Coastal Ridge while the lowest values were in low elevation areas on the margin of the Everglades to the west. The higher values of anisotropy found on the ridge may be due to increased dissolution rates of the oolitic facies of the Miami formation limestone compared with the bryozoan facies to the west. The predominate trend of minimum resistivity and maximum hydraulic conductivity was E‐W/SE‐NW beneath the ridge and E‐W/SW‐NE farther west. The anisotropy directions are similar to the predevelopment groundwater flow direction as indicated in published studies. This suggests that the observed anisotropy is related to the paleo‐groundwater flow in the Biscayne Aquifer.     

Hydrogeophysical Methods for Analyzing Aquifer Storage and Recovery Systems

Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.     

Aquifer Characterization Using Fiber Bragg Grating Multi-Level Monitoring System

Groundwater responses measured from multiple wells at different depths are essential for delineating the aquifer heterogeneity using hydraulic tomography (HT). In general, conducting HT requires many wells because traditional well monitoring is usually partially open at a specific depth interval or is fully penetrating. Accordingly, conducting an HT survey is typically costly and time-consuming. To tackle these issues, a new multi-level monitoring system (MLMS) for the HT survey was developed using the fiber Bragg grating (FBG) technique. This FBG MLMS could collect the depth-discrete groundwater observations from a fully penetrated 2-inch well. Three field campaigns were conducted to validate the capability of the FBG MLMS for HT surveys. The results show that the accuracy and stability of this MLMS are reliable and that FBG MLMS is beneficial for conducting an HT survey. Specifically, compared to the traditional monitoring well in an injection event, this FBG MLMS can concurrently cause an increase in the number of cross-hole tests several times and collect many more head observations than the standard methods, resulting in the observed flow fields efficiently reaching ergodic conditions and effectively improving the accuracy of the estimated hydraulic heterogeneity. Therefore, the FBG MLMS could be an alternative MLMS for efficiently and economically conducting an HT survey.     

Evaluation of Combined Direct-Push Methods Used for Aquifer Model Generation

Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.     

Characterization of Flow Dynamics and Vulnerability in a Coastal Aquifer System

Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ¹³C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ¹⁸O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.