The Hydraulic Profiling Tool for Hydrogeologic Investigation of Unconsolidated Formations

The hydraulic profiling tool (HPT) has become one of the basic tools for investigation of soils and unconsolidated formations over the last 10 years. The HPT is advanced into the subsurface using direct push methods. Clean water is injected into the formation from a small screened port on the side of the probe as it is steadily advanced into the subsurface. A downhole pressure sensor detects the pressure required to inject the water into the formation while an up-hole flowmeter monitors the water flow rate. An electrical conductivity (EC) array included in the lower end of the probe provides a simultaneous EC log of the bulk formation. The EC log, HPT pressure, and flow rate are logged and displayed onscreen as the probe is advanced. These logs enable the investigator to evaluate vertical changes in relative formation permeability at high resolution. Pressure dissipation tests may be performed at selected depths in coarse-grained materials to determine the piezometric pressure in saturated formations. This enables the operator to define the piezometric profile and determine the piezometric surface without a well. Post processing of the log in the viewing software provides for calculation of the corrected HPT pressure (Pc) and estimation of hydraulic conductivity (Est. K) within limits (~0.1 to 75 ft/d). In clean, coarse-grained materials the tandem EC log may be used to estimate groundwater specific conductance based on an Archie's Law model. Cross sections of HPT logs provide an efficient means to define hydrostratigraphy. When combined with contaminant logging tools such as the membrane interface probe (MIP) the HPT data may help to define contaminant migration pathways or contaminated low permeability zones that may result in back diffusion. The HPT can be a useful tool for many geoenvironmental investigations in unconsolidated formations.     

Controls on deep critical zone architecture: a historical review and four testable hypotheses

The base of Earth's critical zone (CZ) is commonly shielded from study by many meters of overlying rock and regolith. Though deep CZ processes may seem far removed from the surface, they are vital in shaping it, preparing rock for infusion into the biosphere and breaking Earth materials down for transport across landscapes. This special issue highlights outstanding challenges and recent advances of deep CZ research in a series of articles that we introduce here in the context of relevant literature dating back to the 1500s. Building on several contributions to the special issue, we highlight four exciting new hypotheses about factors that drive deep CZ weathering and thus influence the evolution of life‐sustaining CZ architecture. These hypotheses have emerged from recently developed process‐based models of subsurface phenomena including: fracturing related to subsurface stress fields; weathering related to drainage of bedrock under hydraulic head gradients; rock damage from frost cracking due to subsurface temperature gradients; and mineral reactions with reactive fluids in subsurface chemical potential gradients. The models predict distinct patterns of subsurface weathering and CZ thickness that can be compared with observations from drilling, sampling and geophysical imaging. We synthesize the four hypotheses into an overarching conceptual model of fracturing and weathering that occurs as Earth materials are exhumed to the surface across subsurface gradients in stress, hydraulic head, temperature, and chemical potential. We conclude with a call for a coordinated measurement campaign designed to comprehensively test the four hypotheses across a range of climatic, tectonic and geologic conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Integrity of Clay Till Aquitards to DNAPL Migration: Assessment Using Current and Emerging Characterization Tools

Field investigations were carried out to determine the occurrence of tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL), the source zone architecture and the aquitard integrity at a 30‐ to 50‐year old DNAPL release site. The DNAPL source zone is located in the clay till unit overlying a limestone aquifer. The DNAPL source zone architecture was investigated through a multiple‐lines‐of‐evidence approach using various characterization tools; the most favorable combination of tools for the DNAPL characterization was geophysical investigations, membrane interface probe, core subsampling with quantification of chlorinated solvents, hydrophobic dye test with Sudan IV, and Flexible Liner Underground Technologies (FLUTe) NAPL liners with activated carbon felt (FACT). While the occurrence of DNAPL was best determined by quantification of chlorinated solvents in soil samples supported by the hydrophobic dye tests (Sudan IV and NAPL FLUTe), the conceptual understanding of source zone architecture was greatly assisted by the indirect continuous characterization tools. Although mobile or high residual DNAPL (S _t > 1%) only occurred in 11% of the source zone samples (intact cores), they comprised 86% of the total PCE mass. The dataset, and associated data analysis, supported vertical migration of DNAPL through fractures in the upper part of the clay till, horizontal migration along high permeability features around the redox boundary in the clay till, and to some extent vertical migration through the fractures in the reduced part of the clay till aquitard to the underlying limestone aquifer. The aquitard integrity to DNAPL migration was found to be compromised at a thickness of reduced clay till of less than 2 m.     

The Utility of Multiple-Completion Monitoring Wells for Describing a Solvent Plume

At a study site in the midwestern United States, multiple-completion wells demonstrated that a vertical hydraulic gradient was responsible for the contamination pattern exhibited by chlorinated solvent plumes. The typical pattern consisted of little or no contamination in the upper portion of the aquifer with concentrations increasing with depth. When ground water contamination was discovered in an unexpected portion of the site, water level elevations and contaminant distribution data obtained from multiple-completion wells resulted in identification of the source location. The well eventually determined to be located in the source area displayed contaminant levels much higher in the upper zone of the aquifer — the opposite contamination pattern of other on-site wells. Such results indicated that the spill had occurred near this location and that solvent residing along the capillary fringe was continuing to contaminate the aquifer.     

Persistence of a Groundwater Contaminant Plume after Hydraulic Source Containment at a Chlorinated‐Solvent Contaminated Site

The objective of this study was to characterize the behavior of a groundwater contaminant (trichloroethene, TCE) plume after implementation of a source‐containment operation at a site in Arizona. The plume resides in a quasi‐three‐layer system comprising a sand/gravel unit bounded on the top and bottom by relatively thick silty clayey layers. The system was monitored for 60 months beginning at start‐up in 2007 to measure the change in contaminant concentrations within the plume, the change in plume area, the mass of the contaminant removed, and the integrated contaminant mass discharge (CMD). The concentrations of TCE in groundwater pumped from the plume extraction wells have declined significantly over the course of operation, as have concentrations for groundwater sampled from 40 monitoring wells located within the plume. The total CMD associated with operation of the plume extraction wells peaked at 0.23 kg/d, decreased significantly within 1 year, and thereafter began an asymptotic decline to a current value of approximately 0.03 kg/d. Despite an 87% reduction in contaminant mass and a comparable 87% reduction in CMD for the plume, the spatial area encompassed by the plume has decreased by only approximately 50%. This is much less than would be anticipated based on ideal flushing and mass‐removal behavior. Simulations produced with a simplified three‐dimensional (3D) numerical model matched reasonably well to the measured data. The results of the study suggest that permeability heterogeneity, back diffusion, hydraulic factors associated with the specific well field system, and residual discharge from the source zone are all contributing to the observed persistence of the plume, as well as the asymptotic behavior currently observed for mass removal and for the reduction in CMD.     

Field Application of the Combined Membrane‐Interface Probe and Hydraulic Profiling Tool (MiHpt)

The Membrane‐Interface Probe and Hydraulic Profiling Tool (MiHpt) is a direct push probe that includes both the membrane interface probe (MIP) and hydraulic profiling tool (HPT) sensors. These direct push logging tools were previously operated as separate logging systems for subsurface investigation in unconsolidated formations. By combining these two probes into one logging system the field operator obtains useful data about the distribution of both volatile organic contaminants (VOCs) and relative formation permeability in a single boring. MiHpt logging was conducted at a chlorinated VOC contaminated site in Skuldelev, Denmark, to evaluate performance of the system. Formation cores and discrete interval slug tests are used to assess use of the HPT and electrical conductivity (EC) logs for lithologic and hydrostratigraphic interpretation. Results of soil and groundwater sample analyses are compared to the adjacent MiHpt halogen specific detector (XSD) logs to evaluate performance of the system to define contaminant distribution and relative concentrations for the observed VOCs. Groundwater profile results at moderate to highly contaminated locations were found to correlate well with the MiHpt‐XSD detector responses. In general, soil sample results corresponded with detector responses. However, the analyses of saturated coarse‐grained soils at the site proved to be unreliable as demonstrated by high RPDs for duplicate samples. The authors believe that this is due to pore water drainage observed from these cores during sampling. Additionally, a cross section of HPT pressure and MiHpt‐XSD detector logs provides insight into local hydrostratigraphy and formation control on contaminant migration.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Tree Coring as a Complement to Soil Gas Screening to Locate PCE and TCE Source Zones and Hot Spots

Preliminary risk assessment for prioritisation of site investigations requires efficient screening to reveal type and level of contamination. The screening methods, tree coring and soil gas sampling were applied and compared at two forested sites contaminated with tetrachloroethylene (PCE) or trichloroethylene (TCE) to evaluate their ability to locate source zones and contaminant hot spots. One test site represented a relatively homogeneous sandy soil and aquifer, and the second a more heterogeneous geology with both sandy and less permeable clay till layers overlying a chalk aquifer. Tree cores from different tree species were sampled and analysed, and compared to soil gas measurements and existing soil gas data. Both methods were found useful as screening tools to locate hot spots of PCE and TCE in the shallow subsurface. Tree coring was found to be particularly beneficial as a complement to soil gas sampling at sites with low permeable soils, and where contamination was located in the capillary rise or shallow groundwater. The shorter time required for tree coring reduced the costs compared to soil gas sampling, but the sensitivity and precision of tree coring were lower. However, this did not affect the feasibility of using tree coring to locate the hot spots. Moreover, a combination of the two methods can help to focus any subsequent investigations like soil or groundwater sampling. The use of tree coring to complement soil gas sampling for pre‐screening is expected to result in higher certainty for revealing hot spots and source zones at contaminated sites.     

A rapid method for hydraulic profiling in unconsolidated formations

Information on vertical variations in hydraulic conductivity ( K ) can often shed much light on how a contaminant will move in the subsurface. The direct-push injection logger has been developed to rapidly obtain such information in shallow unconsolidated settings. This small-diameter tool consists of a short screen located just behind a drive point. The tool is advanced into the subsurface while water is injected through the screen to keep it clear. Upon reaching a depth at which information about K is desired, advancement ceases and the injection rate and pressure are measured on the land surface. The rate and pressure values are used in a ratio that serves as a proxy for K . A vertical profile of this ratio can be transformed into a K profile through regressions with K estimates determined using other techniques. The viability of the approach was assessed at an extensively studied field site in eastern Germany. The assessment demonstrated that this tool can rapidly identify zones that may serve as conduits for or barriers to contaminant movement.     

Vertically Integrated Hydraulic Conductivity: A New Parameter for Groundwater-Surface Water Analysis

Numerical modeling of groundwater-surface water interactions provides vital information necessary for determining the extent of nutrient transport, quantifying water budgets, and delineating zones of ecological support. The hydrologic data that drive these models are often collected at disparate scales and subsequently incorporated into numerical models through upscaling techniques such as piecewise constancy or geostatistical methods. However, these techniques either use basic interpolation methods, which often simplifies the system of interest, or utilize complex statistical methods that are computationally expensive, time consuming, and generate complex subsurface configurations. We propose a bulk parameter termed “vertically integrated hydraulic conductivity” (K_V), and defined as the depth-integrated resistance to fluid flow sensed at the groundwater-surface water interface, as an alternative to hydraulic conductivity when investigating vertical fluxes across the groundwater-surface water interface. This bulk parameter replaces complex subsurface configurations in situations dominated by vertical fluxes and where heterogeneity is not of primary importance. To demonstrate the utility of K_V, we extracted synthetic temperature time series data from a forward numerical model under a variety of scenarios and used those data to quantify vertical fluxes using the amplitude ratio method. These quantified vertical fluxes and the applied hydraulic head gradient were subsequently input into Darcy's Law and used to quantify K_V. This K_V was then directly compared to the equivalent hydraulic conductivity (K_T) assuming an infinitely extending layer. Vertically integrated hydraulic conductivity allows for more accurate and robust flow modeling across the groundwater-surface water interface in instances where complex heterogeneities are not of primary concern.     

The effects of vertical barrier walls on the hydraulic control of contaminated groundwater

We present explicit analytic solutions describing the hydraulic head and discharge vector for two-dimensional, steady groundwater flow past an impermeable barrier embedded in a regional flow field. We use the solution to investigate the effects of open vertical barriers on the flow field; in particular, we examine the hydraulic containment of contaminant plumes or source zones by combination of a vertical barrier wall and extraction wells. We quantify the local reduction in discharge rates due to the barrier wall and the local increase in the size of the capture zone of an extraction well near an open, up-gradient barrier. We find that the combination of an open vertical barrier with down-gradient extraction wells can be very effective in decreasing the well discharge rate necessary to control a contaminant plume or source area. Design charts are presented for quantifying the effects of the barrier wall on the hydraulic control of the groundwater flow field and for estimating the jump in head across a barrier. The charts are appropriate for use in the preliminary design and cost estimating of remedial systems, and for the design of dewatering systems.     

Development of an electronic seepage chamber for extended use in a river

Seepage chambers have been used to characterize the flux of water across the water-sediment interface in a variety of settings. In this work, an electronic seepage chamber was developed specifically for long-term use in a large river where hydraulic gradient reversals occur frequently with river-stage variations. A bidirectional electronic flowmeter coupled with a seepage chamber was used to measure temporal changes in the magnitude and direction of water flux across the water-sediment interface over an 8-week period. The specific discharge measured from the seepage chamber compared favorably with measurements of vertical hydraulic gradient and previous specific discharge calculations. This, as well as other supporting data, demonstrates the effectiveness of the electronic seepage chamber to accurately quantify water flux in two directions over a multimonth period in this setting. The ability to conduct multimonth measurements of water flux at a subhourly frequency in a river system is a critical capability for a seepage chamber in a system where hydraulic gradients change on a daily and seasonal basis.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

An alternative approach in establishing relation between vertical and horizontal gradients of 2D potential field

The relation in which the vertical and horizontal gradients of potential field data measured along a profile across a two‐dimensional source are a Hilbert transform pair is re‐established using complex domain mathematics. In addition, a relation between the measured field and its vertical gradient in terms of a closed‐form formula is also established. The formula is based on hypersingular or Hadamard's finite‐part integral. To estimate the vertical gradient directly from the field data, Linz's algorithm of computing Hadamard's finite‐part integral is implemented. Numerical experiments are conducted on synthetically generated total magnetic intensity data with a mild level of noise contamination. A model of a magnetically polarised vertical thin sheet buried at a finite depth within a non‐magnetic half‐space was considered in generating the synthetic response. The results from numerical experiments on the mildly noise‐contaminated synthetic response are compared with those from using classical Fourier and robust regularised Hilbert transform‐based techniques.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Microtopographical and hydrophysical controls on subsurface flow and solute transport: A continuous solute release experiment in a subarctic bog

Resource extraction and transportation activities in subarctic Canada can result in the unintentional release of contaminants into the surrounding peatlands. In the event of a release, a thorough understanding of solute transport within the saturated zone is necessary to predict plume fate and the potential impacts on peatland ecosystems. To better characterize contaminant transport in these systems, approximately 13,000 L/day of sodium chloride tracer (200 mg/L) was released into a bog in the James Bay Lowland. The tracer was pumped into a fully penetrating well (1.5 m) between July 5 and August 18, 2015. Horizontal and vertical plume development was measured via in situ specific conductance and water table depth from an adaptive monitoring network. Over the spill period, the bulk of the plume travelled a lateral distance of 100 m in the direction of the slight regional groundwater and topographical slope. The plume shape was irregular and followed the hollows, indicating preferential flow paths due to the site microtopography. Saturated transport of the tracer occurred primarily at ~25 cm below ground surface (bgs), and at a discontinuous high hydraulic conductivity layer ~125 cm bgs due to a complex and heterogeneous vertical hydraulic conductivity profile. Plume measurement was confounded by a large amount of precipitation (233 mm over the study period) that temporarily diluted the tracer in the highly conductive upper peat layer. Longitudinal solute advection can be approximated using local water table information (i.e., depth and gradient); microtopography; and meteorological conditions. Vertical distribution of solute within the peat profile is far more complex due to the heterogeneous subsurface; characterization would be aided by a detailed understanding of the site‐specific peat profile; the degree of decomposition; and the type of contaminant (e.g., reactive/nonreactive). The results of this research highlight the difficulty of tracking a contaminant spill in bogs and provide a benchmark for the characterization of the short‐term fate of a plume in these complex systems.     

From Grain to Floodplain: Evaluating heterogeneity of floodplain hydrostatigraphy using sedimentology,geophysics, and remote sensing

Stratigraphy is a fundamental component of floodplain heterogeneity and hydraulic conductivity and connectivity of alluvial aquifers, which affect hydrologic processes such as groundwater flow and hyporheic exchange. Watershed-scale hydrological models commonly simplify the sedimentology and stratigraphy of floodplains, neglecting natural floodplain heterogeneity and anisotropy. This study, conducted in the upper reach of the East River in the East River Basin, Colorado, USA, combines point-, meander-, and floodplain-scale data to determine key features of alluvial aquifers important for estimating hydrologic processes. We compare stratigraphy of two meanders with disparate geometries to explore floodplain heterogeneity and connectivity controls on flow and transport. Meander shape, orientation, and internal stratigraphy affected residence time estimates of laterally exchanged hyporheic water. Although the two meanders share a sediment source, vegetation, and climate, their divergent river migration histories resulted in contrasting meander hydrofacies. In turn, the extent and orientation of these elements controlled the effective hydraulic conductivity and, ultimately, estimates of groundwater transport and hyporheic residence times. Additionally, the meanders’ orientation relative to the valley gradient impacted the hydraulic gradient across the meanders—a key control of groundwater velocity. Lastly, we combine our field data with remotely sensed data and introduce a potential approach to estimate key hydrostratigraphic packages across floodplains. Prospective applications include contaminant transport studies, hyporheic models, and watershed models. © 2019 John Wiley & Sons, Ltd.     

Direct-Push Electrical Conductivity Logging for High-Resolution Hydrostratigraphic Characterization

Fine-scale hydrostratigraphic features often play a critical role in controlling ground water flow and contaminant transport. Unfortunately, many conventional drilling- and geophysics-based approaches are rarely capable of describing these features at the level of detail needed for contaminant predictions and remediation designs. Previous work has shown that direct-push electrical conductivity (EC) logging can provide information about site hydrostratigraphy at a scale of relevance for contaminant transport investigations in many unconsolidated settings. In this study, we evaluate the resolution and quality of that information at a well-studied research site that is underlain by highly stratified alluvial sediments. Geologic and hydrologic data, conventional geophysical logs, and particle-size analyses are used to demonstrate the capability of direct-push EC logging for the delineation of fine-scale hydrostratigraphic features in saturated unconsolidated formations. When variations in pore-fluid chemistry are small, the electrical conductivity of saturated media is primarily a function of clay content, and hydrostratigraphic features can be described at a level of detail (<2.5 cm in thickness) that has not previously been possible in the absence of continuous cores. Series of direct-push EC logs can be used to map the lateral continuity of layers with non-negligible clay content and to develop important new insights into flow and transport at a site. However, in sand and gravel intervals with negligible clay, EC logging provides little information about hydrostratigraphic features. As with all electrical logging methods, some site-specific information about the relative importance of fluid and sediment contributions to electrical conductivity is needed. Ongoing research is directed at developing direct-push methods that allow EC logging, water sampling, and hydraulic testing to be done concurrently.     

Case Histories of Shallow Time Domain Electromagnetics in Environmental Site Assessment

Time domain electromagnetic (TDEM) soundings were utilized in site assessment studies for the purpose of accomplishing two objectives: (1) measuring the migration of brines from oil and gas field evaporation pits, and (2) mapping the continuity of clay strata. Both case histories are representative of common objectives in site assessment. The data for both examples were analyzed by one-dimensional ridge regression inversions.
Both case histories illustrate that TDEM is an effective method for determining the lateral and vertical resistivity distribution (geoelectric section) of the subsurface in the depth range from about 5 m to 100 m, and it is known from other investigations that depths of investigations up to 2500 m can be achieved with TDEM.
Frequently, the lateral extent of contaminant plumes emanating from localized sources are of limited areal extent. The case history involving the use of TDEM to map a brine plume shows that a TDEM measurement has a relatively small zone of influence, so that meaningful information about the lateral and vertical extent of the plume can be obtained. Both case histories demonstrate the ability of TDEM to determine geoelectric sections below layers of low resistivity.     

Transition probability/Markov chain analyses of DNAPL source zones and plumes

At sites where a dense nonaqueous phase liquid (DNAPL) was spilled or released into the subsurface, estimates of the mass of DNAPL contained in the subsurface from core or monitoring well data, either in the nonaqueous or aqueous phase, can be highly uncertain because of the erratic distribution of the DNAPL due to geologic heterogeneity. In this paper, a multiphase compositional model is applied to simulate, in detail, the DNAPL saturations and aqueous-phase plume migration in a highly characterized, heterogeneous glaciofluvial aquifer, the permeability and porosity data of which were collected by researchers at the University of Tübingen, Germany. The DNAPL saturation distribution and the aqueous-phase contaminant mole fractions are then reconstructed by sampling the data from the forward simulation results using two alternate approaches, each with different degrees of sampling conditioning. To reconstruct the DNAPL source zone architecture, the aqueous-phase plume configuration, and the contaminant mass in each phase, one method employs the novel transition probability/Markov chain approach (TP/MC), while the other involves a traditional variogram analysis of the sampled data followed by ordinary kriging. The TP/MC method is typically used for facies and/or hydraulic conductivity reconstruction, but here we explore the applicability of the TP/MC method for the reconstruction of DNAPL source zones and aqueous-phase plumes. The reconstructed geometry of the DNAPL source zone, the dissolved contaminant plume, and the estimated mass in each phase are compared using the two different geostatistical modeling approaches and for various degrees of data sampling from the results of the forward simulation. It is demonstrated that the TP/MC modeling technique is robust and accurate and is a preferable alternative compared to ordinary kriging for the reconstruction of DNAPL saturation patterns and dissolved-phase contaminant plumes.