Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirectobservations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. We demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

Dual in-aquifer and near surface processes drive arsenic mobilization in Cambodian groundwaters

Millions of people globally, and particularly in South and Southeast Asia, face chronic exposure to arsenic from reducing groundwaters in which. Arsenic release to is widely attributed largely to reductive dissolution of arsenic-bearing iron minerals, driven by metal reducing bacteria using bioavailable organic matter as an electron donor. However, the nature of the organic matter implicated in arsenic mobilization, and the location within the subsurface where these processes occur, remains debated. In a high resolution study of a largely pristine, shallow aquifer in Kandal Province, Cambodia, we have used a complementary suite of geochemical tracers (including ¹⁴C, ³H, ³He, ⁴He, Ne, δ¹⁸O, δD, CFCs and SF₆) to study the evolution in arsenic-prone shallow reducing groundwaters along dominant flow paths. The observation of widespread apparent ³H-³He ages of <55years fundamentally challenges some previous models which concluded that groundwater residence times were on the order of hundreds of years. Surface-derived organic matter is transported to depths of >30m, and the relationships between age-related tracers and arsenic suggest that this surface-derived organic matter is likely to contribute to in-aquifer arsenic mobilization. A strong relationship between ³H-³He age and depth suggests the dominance of a vertical hydrological control with an overall vertical flow velocity of ~0.4±0.1m·yr^-1 across the field area. A calculated overall groundwater arsenic accumulation rate of ~0.08±0.03μM·yr^-1 is broadly comparable to previous estimates from other researchers for similar reducing aquifers in Bangladesh. Although apparent arsenic groundwater accumulation rates varied significantly with site (e.g. between sand versus clay dominated sequences), rates are generally highest near the surface, perhaps reflecting the proximity to the redox cline and/or depth-dependent characteristics of the OM pool, and confounded by localized processes such as continued in-aquifer mobilization, sorption/desorption, and methanogenesis.

High resolution profile of inorganic aqueous geochemistry and key redox zones in an arsenic bearing aquifer in Cambodia

Arsenic contamination of groundwaters in South and Southeast Asia is a major threat to public health. In order to better understand the geochemical controls on the mobility of arsenic in a heavily arsenic-affected aquifer in northern Kandal Province, Cambodia, key changes in inorganic aqueous geochemistry have been monitored at high vertical and lateral resolution along dominant groundwater flow paths along two distinct transects. The two transects are characterized by differing geochemical, hydrological and lithological conditions. Arsenic concentrations in groundwater are highly heterogenous, and are broadly positively associated with iron and negatively associated with sulfate and dissolved oxygen. The observed correlations are generally consistent with arsenic mobilization by reductive-dissolution of iron (hydr)oxides. Key redox zones, as identified using groupings of the PHREEQC model equilibrium electron activity of major redox couples (notably ammonium/nitrite; ammonium/nitrate; nitrite/nitrate; dissolved oxygen/water) have been identified and vary with depth, site and season. Mineral saturation is also characterized. Seasonal changes in groundwater chemistry were observed in areas which were (i) sandy and of high permeability; (ii) in close proximity to rivers; and/or (iii) in close proximity to ponds. Such changes are attributed to monsoonal-driven surface-groundwater interactions and are consistent with the separate provenance of recharge sources as identified using stable isotope mixing models.

Geoelectrical monitoring of simulated subsurface leakage to support high-hazard nuclear decommissioning at the Sellafield Site, UK

A full-scale field experiment applying 4D (3D time-lapse) cross-borehole Electrical Resistivity Tomography (ERT) to the monitoring of simulated subsurface leakage was undertaken at a legacy nuclear waste silo at the Sellafield Site, UK. The experiment constituted the first application of geoelectrical monitoring in support of decommissioning work at a UK nuclear licensed site. Images of resistivity changes occurring since a baseline date prior to the simulated leaks revealed likely preferential pathways of silo liquor simulant flow in the vadose zone and upper groundwater system. Geophysical evidence was found to be compatible with historic contamination detected in permeable facies in sediment cores retrieved from the ERT boreholes. Results indicate that laterally discontinuous till units forming localized hydraulic barriers substantially affect flow patterns and contaminant transport in the shallow subsurface at Sellafield. We conclude that only geophysical imaging of the kind presented here has the potential to provide the detailed spatial and temporal information at the (sub-)meter scale needed to reduce the uncertainty in models of subsurface processes at nuclear sites.