Vadose Zone Monitoring: Experiences and Trends in the United States

The vadose zone is the portion of the geologic profile above a perennial aquifer. Inclusion of mandatory vadose zone monitoring techniques as an approach to aquifer protect ion was first proposed under the Resource Conservation and Recovery Act in the United States in 1978 and has since received increasing acceptance at federal and stale levels. The goals of a vadose zone characterization and monitoring effort are to establish background conditions, identify contaminant transport pathways, identify the extent and degree of existing contamination, establish the basis for monitoring network design, measure the parameters needed in a risk assessment, and provide detection of contaminant migration toward ground water resources. The benefits of vadose zone monitoring include early warning of contaminant migration, potential reduction of ground water monitoring efforts, reduction of contaminant spreading and volume, and reduced time and cost of remediation once a contaminant release occurs. Vadose zone characterization and monitoring techniques should be considered as critical hydrologic tools in the prevention of ground water resource degradation.     

Electrical Leak Detection System for Landfill Liners: A Case History

As landfill specifications become more stringent in the United Kingdom, the development of increasingly sophisticated monitoring methods is necessary to meet environmental protection goals. This case history describes the development of a 2-million-cubic-meter-capacity landfill located in a sandstone quarry and 1 km from a public water supply borehole, where the sensitivity of the site to ground water contamination and the proximity to a public water supply borehole are particular issues.
The landfill design incorporated a more sensitive environmental monitoring system, using a geophysical technique. The monitoring system comprises a permanent grid of electrodes installed beneath the landfill, connected by multicore cable to a computer-controlled earth resistance meter and switching unit in the site weighbridge. It was designed to detect holes in the landfill liner prior to and after covering with waste and to monitor the migration of contaminants beneath the landfill before they reach the perimeter observation boreholes, should leakage occur.
Such monitoring can enable the integrity of the landfill to be routinely reviewed; holes can be repaired if they are readily accessible and, if not, monitoring provides an early warning to enable the implementation of any additional monitoring or corrective action, based on the environmental risk posed by the site.
The system was first used as a quality assurance test once the landfill liner, which covered an area of 3 hectares, was installed. The system proved sensitive, detecting a hole consisting of two narrow knife cuts. Such sensitivity allows a high degree of confidence to be placed upon the integrity of the liner resulting in a significant contribution to public reassurance. The landfill is now operational, and monitoring using the geophysical system will be undertaken on a monthly basis for the first year, with the frequency of monitoring reviewed thereafter.     

Groundwater monitoring strategies for variable contaminant migration boundaries

Abstract

Detection efficiencies of alternative groundwater monitoring networks were evaluated in relation to distance to a buffer zone (contaminant migration) boundary. This boundary establishes a distance limit within which contaminant plumes should pass through monitoring wells, located on curvilinear segments (monitoring loci) near a waste storage facility. Alternative strategies allocated monitoring wells to loci at specified distances, measured parallel to groundwater flow, from the downgradient boundaries of a landfill. One approach constrained wells to equal spacing, measured perpendicular to groundwater flow. Compressing well locations 10% closer to the downgradient corner of the landfill rendered alternative monitoring configurations. Computations by a monitoring efficiency model indicated: (a) networks largely maintained detection efficiency for different contaminant migration boundaries; (b) one network most efficiently attained a target detection capability for all contaminant migration boundaries; and (c) compressed networks slightly outperformed equal-spaced counterparts. Compressed networks with more wells along closer monitoring loci best maintained the detection efficiency when shifting the contaminant migration boundary closer to the landfill. Procedures described in this paper may be useful for examining trade-offs between monitoring efficiency and distance limits of contaminant travel at landfills posing potential hazards to underlying groundwater.     

Relationship Between Subsurface Landfill Gas and Arsenic Mobilization into Groundwater

The maximum contaminant level for arsenic was reduced by the U.S. Environmental Protection Agency (U.S. E.P.A.) for Drinking Water Standards from 50 micrograms per liter (µg/L) to 10 µg/L, effective January 23, 2006. The subject site is a double-lined sanitary landfill facility located in the mid-Atlantic region of the United States. Arsenic was reported above the maximum contaminant level in a downgradient monitoring well (MW-18) in July 2005. Since July 2005, arsenic levels in MW-18 fluctuated above and below 10 µg/L. This research focuses on determining whether reducing conditions in groundwater, enhanced by subsurface landfill gas emissions, were causing naturally-occurring arsenic to mobilize from the native variably-saturated vadose zone soils into groundwater. The groundwater data collected from the impacted well (MW-18) were compared to an upgradient well (MW-8) to determine whether significant differences existed during the time period of April 2004 to April 2007. Linear regression analysis was also used to determine whether other parameters had a significant relationship with the arsenic concentrations detected in MW-18. The groundwater located in MW-18 was consistently more reduced than groundwater located in the upgradient/background well MW-8, and this was most likely attributed to the presence of subsurface landfill gas in the area. According to the U.S. E.P.A., oxidation-reduction potential (ORP) values in groundwater less than 50 millivolts (mv) suggest that a reducing environment may be present. The data presented in this study indicate that arsenic can mobilize into groundwater under moderately reducing conditions, with ORP measurements averaging 53 mv.     

High Purge Volume Sampling—A New Paradigm for Subslab Soil Gas Monitoring

Subslab soil gas sampling and analysis is a common line of evidence for assessing human health risks associated with subsurface vapor intrusion to indoor air for volatile organic compounds; however, conventional subslab sampling methods have generated data that show substantial spatial and temporal variability, which often makes the interpretation difficult. A new method of monitoring has been developed and tested that is based on a concept of integrating samples over a large volume of soil gas extracted from beneath the floor slab of a building to provide a spatially averaged subslab concentration. Regular field screening is also conducted to assess the trend of concentration as a function of the volume removed to provide insight into the spatial distribution of vapors at progressive distances away from the point of extraction. This approach minimizes the risk of failing to identify the areas of elevated soil vapor concentrations that may exist between discrete sample locations, and can provide information covering large buildings with fewer holes drilled through the floor. The new method also involves monitoring the extraction flow rate and transient vacuum response for mathematical analysis to help interpret the vapor concentration data and to support an optimal design for any subslab venting system that may be needed.     

Utilization of Shallow Geophysical Sensing at Two Abandoned Municipal/Industrial Waste Landfills on the Missouri River Floodplain

In the past 30 to 40 years, floodplain areas of large rivers, such as the Missouri River, have been extensively used for large industrial and municipal landfills. Many of these sites are now causing varying degrees of ground water contamination. Rapid geophysical characterization techniques have proven useful for delineation of anomalous areas indicative of potential contaminant plumes. These methods have also resulted in a cost effective approach to the location and number of monitoring wells.
An effective technique to initially characterize ground water contamination at such landfills along the Missouri River in northwestern Missouri involved a combination of electrical resistivity and electromagnetic conductivity methods. Resistivity was used to obtain soundings of the alluvium by using a modified Wenner array and to corroborate shallow electromagnetic conductivity measurements by using short Wenner array electrode spacings.
Upon confirmation of similar measurements of the upper soils for the two methods, numerous electromagnetic conductivity traverses were made at each landfill site. The data generated from these surveys were graphed and contoured to delineate anomalous areas. Based on the geophysical study, a ground water monitoring well network was then designed for each landfill.
As a result, a minimal number of wells were required to initially characterize the ground water quality at these two sites. In general, analysis of water samples from these wells displayed good correlation with the geophysical results.     

Landslide activation behaviour illuminated by electrical resistance monitoring

A common factor in landslide activation (or reactivation) is subsurface moisture and associated pore pressure variations linked to rainfall. Monitoring of these subsurface hydrogeological processes is necessary to improve our understanding of water‐induced landslide activation. Geophysical approaches, electrical methods in particular, are increasingly being applied to landslide monitoring because they provide non‐invasive spatial information in heterogeneous subsurface environments that can be difficult to characterise using surface observations or intrusive sampling alone. Electrical techniques are sensitive to changing subsurface moisture conditions, and have proven to be a useful tool for investigating the hydrogeology of natural and engineered slopes. The objectives of this investigation were to further develop electrical resistance monitoring for slope stability assessment, and to validate the approach at an intermittently‐active UK landslide system to advance the understanding of complex landslide activation mechanisms. A long‐term transfer resistance dataset was collected from a grid of electrodes to allow spatial monitoring of the landslide. These data were interpreted using a synthesis of rainfall, temperature, GPS and piezometric records. The resistance data were corrected for seasonal temperature variations and electrode movements were monitored, as these processes were shown to mask moisture related changes. Results reveal that resistance monitoring is sensitive to soil moisture accumulation, including changes in piezometric levels, and can be used to study the principal activation mechanism of slow‐moving shallow earthflows. Spatial monitoring using resistance maps was shown to be particularly valuable as it revealed the evolution of subsurface moisture distribution, in the lead up to landslide activation. Key benefits of this approach are that it provides a simple, rapid and non‐invasive means of spatially monitoring subsurface moisture dynamics linked to landslide activation at high‐temporal resolution. Crucially, it provides a means of monitoring subsurface hydraulic changes in the build‐up to slope failure, thereby contributing to early warning of landslide events. Copyright © 2018 John Wiley & Sons, Ltd.     

Dual states estimation of a subsurface flow-transport coupled model using ensemble Kalman filtering

Modeling the spread of subsurface contaminants requires coupling a groundwater flow model with a contaminant transport model. Such coupling may provide accurate estimates of future subsurface hydrologic states if essential flow and contaminant data are assimilated in the model. Assuming perfect flow, an ensemble Kalman filter (EnKF) can be used for direct data assimilation into the transport model. This is, however, a crude assumption as flow models can be subject to many sources of uncertainty. If the flow is not accurately simulated, contaminant predictions will likely be inaccurate even after successive Kalman updates of the contaminant model with the data. The problem is better handled when both flow and contaminant states are concurrently estimated using the traditional joint state augmentation approach. In this paper, we introduce a dual estimation strategy for data assimilation into a one-way coupled system by treating the flow and the contaminant models separately while intertwining a pair of distinct EnKFs, one for each model. The presented strategy only deals with the estimation of state variables but it can also be used for state and parameter estimation problems. This EnKF-based dual state-state estimation procedure presents a number of novel features: (i) it allows for simultaneous estimation of both flow and contaminant states in parallel; (ii) it provides a time consistent sequential updating scheme between the two models (first flow, then transport); (iii) it simplifies the implementation of the filtering system; and (iv) it yields more stable and accurate solutions than does the standard joint approach. We conducted synthetic numerical experiments based on various time stepping and observation strategies to evaluate the dual EnKF approach and compare its performance with the joint state augmentation approach. Experimental results show that on average, the dual strategy could reduce the estimation error of the coupled states by 15% compared with the joint approach. Furthermore, the dual estimation is proven to be very effective computationally, recovering accurate estimates at a reasonable cost.     

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.     

新型坐底式海洋可控源电磁发射系统及其海试应用

利用拖曳式海洋可控源电磁发射系统在探测埋藏较浅的天然气水合物资源时,会遇到一些新问题.如拖曳式轴向发射偶极源不能贴紧海底,发射偶极源和海底间海水的电磁衰减,使得10 Hz以上的相对高频能量难以导入至海底以下介质;发射偶极源有可能出现水平或垂向摆动;拖曳式发射时,时间窗口内叠加的数据有限;拖曳路线可能与海底构造走向平行,不利于揭示探测目标体的异常形态.新型坐底式发射系统有望解决上述问题.坐底式发射系统的硬件部分包括甲板端供电和监控单元、长距离电力和数据通信单元,以及水下发射机主体.新型坐底式发射系统的发射电极紧贴海底,没有海水层衰减,有利于相对高频的人工源电磁信号经过海底以下介质传输至接收端;通过超短基线信标更容易精确定位发射机拖体,利用姿态方位参考系统确定发射电极的供电方向,可以精确校正电偶源水平偏角的影响;增加单点供电时间,提高接收信号信噪比;提供两对电偶极源,从两个相互垂直的方向对异常体进行人工源激发.2015年海洋试验的结果表明,坐底式与拖曳式发射系统联合作业,可从多角度对海底异常体进行宽频带电磁激发,为数据采集提供丰富的场源信息.     

Configuring Detection Wells Near Landfills

A graphical method was devised for designing contaminant detection monitoring networks in aquifers. The approach eliminates bias in detection efficiency among well pairs, thereby improving the overall efficiency of a ground water monitoring network. In the equidistant configurations derived by the graphical approach, all wells are located the same distance from a landfill, but the distance is measured parallel to ground water flow, Measured perpendicular to ground water flow, there is also an equal spacing between wells in an equidistant network. A simulation model was used to compare an equidistant network to a peripheral monitoring configuration, in which wells were spaced evenly along the downgradient boundaries of a landfill. The equidistant network yielded a 12.4% higher detection efficiency and also facilitated earlier release detection. In practice, the graphical approach that yields equidistant configurations can be used to identify candidate monitoring networks to detect potential releases from landfills.     

Evaluation of Subsurface Oxygen Sensors for Remediation Monitoring

Continuous remediation monitoring using sensors is potentially a more effective and inexpensive alternative to current methods of sample collection and analysis. Gaseous components of a system are the most mobile and easiest to monitor. Continuous monitoring of soil gases such as oxygen, carbon dioxide, and contaminant vapors can provide important quantitative information regarding the progress of bioremediation efforts and the area of influence of air sparging or soil venting. Laboratory and field tests of a commercially available oxygen sensor show that the subsurface oxygen sensor provides rapid and accurate data on vapor phase oxygen concentrations. The sensor is well suited for monitoring gas flow and oxygen consumption in the vadose zone during air sparging and bioventing. The sensor performs well in permeable, unsaturated soil environments and recovers completely after being submerged during temporary saturated conditions. Calibrations of the in situ oxygen sensors were found to be stable after one year of continuous subsurface operation. However, application of the sensor in saturated soil conditions is limited. The three major advantages of this sensor for in situ monitoring arc as follows: (1) it allows data acquisition at any specified time interval; (2) it provides potentially more accurate data by minimizing disturbance of subsurface conditions; and (3) it minimizes the cost of field and laboratory procedures involved in sample retrieval and analysis.     

Target-oriented imaging using extended image volumes: a low-rank factorization approach

Imaging in geological challenging environments has led to new developments, including the idea of generating reflection responses by means of interferometric redatuming at a given target datum in the subsurface, when the target datum lies beneath a complex overburden. One way to perform this redatuming is via conventional model-based wave-equation techniques. But those techniques can be computationally expensive for large-scale seismic problems since the number of wave-equation solves is equal to two times the number of sources involved during seismic data acquisition. Also conventional shot-profile techniques require lots of memory to save full subsurface extended image volumes. Therefore, we can only form subsurface image volumes in either horizontal or vertical directions. To exploit the information hidden in full subsurface extended image volumes, we now present a randomized singular value decomposition-based approach built upon the matrix probing scheme, which takes advantage of the algebraic structure of the extended imaging system. This low-rank representation enables us to overcome both the computational cost associated with the number of wave-equation solutions and memory usage due to explicit storage of full subsurface extended image volumes employed by conventional migration methods. Experimental results on complex geological models demonstrate the efficacy of the proposed methodology and allow practical reflection-based extended imaging for large-scale five-dimensional seismic data.     

Testing the ability of surface arrays to monitor microseismic activity

Recently there has been much interest in the use of data from surface arrays in conjunction with migration‐based processing methods for passive seismic monitoring. In this study we use an example of this kind of data recorded whilst 18 perforation shots, with a variety of positions and propellant amounts, were detonated in the subsurface. As the perforation shots provide signals with known source positions and origin times, the analysis of these data is an invaluable opportunity to test the accuracy and ability of surface arrays to detect and locate seismic sources in the subsurface. In all but one case the signals from the perforation shots are not visible in the raw or preprocessed data. However, clear source images are produced for 12 of the perforation shots showing that arrays of surface sensors are capable of imaging microseismic events, even when the signals are not visible in individual traces. We find that point source locations are within typically 45 m (laterally) of the true shot location, however the depths are less well constrained (～150 m). We test the sensitivity of our imaging method to the signal‐to‐noise ratio in the data using signals embedded in realistic noise. We find that the position of the imaged shot location is quite insensitive to the level of added noise, the primary effect of increased noise being to defocus the source image. Given the migration approach, the array geometry and the nature of coherent noise during the experiment, signals embedded in noise with ratios ≥0.1 can be used to successfully image events. Furthermore, comparison of results from data and synthetic signals embedded in noise shows that, in this case, prestack corrections of traveltimes to account for near‐surface structure will not enhance event detectability. Although, the perforation shots have a largely isotropic radiation pattern the results presented here show the potential for the use of surface sensors in microseismic monitoring as a viable alternative to classical downhole methods.     

Use of LIF for Real-Time In-Situ Mixed NAPL Source Zone Detection

The site characterization and analysis cone penetrometer system (SCAPS), equipped with realtime fluorophore detection capabilities, was used to delineate subsurface contaminant releases in an area where plating shop waste was temporarily stored. Records indicated that various nonaqueous phase liquids (NAPLs) were released at the site. The investigators advanced the SCAPS laser-induced fluorescence (LIF) sensor to depths beneath the water table of the principal water-bearing zone. The water table was located approximately 6 feet (1.8 m) below ground surface (bgs) across the site. Fluorescence, attributed to fuel compounds commingled with chlorinated solvents, was observed at depths ranging from 4.0 to 11.5 feet (1.2 to 3.5 m) bgs. Fluorescence, attributed to naturally occurring organic materials (by process of elimination and spectral characteristics) commingled with chlorinated solvent constituents, was observed at depths ranging from approximately 13 to 40 feet (4.0 to 12.2 m) bgs. Fluorescence responses from compounds confirmed to be commingled with chlorinated solvents indicates that the SCAPS fluorophore detection system is capable of indirectly delineating vadose zone and subaqueous chlorinated solvents and other dense nonaqueous phase liquids (DNAPLs) at contaminant release sites. This confirmation effort represents the first documented account of the successful application of LIF to identify a mixed DNAPL/LNAPL source zone.     

Dreamlet source‐receiver survey sinking prestack depth migration

Survey sinking migration downward continues the entire surface observed multi‐shot data to the subsurface step by step recursively. Reflected energy from reflectors at current depth appear at zero time and zero offset in the extrapolated wavefield. The data (seismic records) of t > 0 at this depth are equivalent to the data acquired by a survey system deployed at this depth. This is the reason to name the process ‘survey sinking’. The records of negative time need not to be further propagated since they carry no information to image structures beneath the new survey system. In this paper, we combine survey sinking with dreamlet migration. The dreamlet migration method decomposes the seismic wavefield and one‐way wave propagator by complete time‐space localized bases. The localization on time gives flexibility on time‐varying operations during depth extrapolation. In dreamlet survey sinking migration, it only keeps the data for imaging the structures beneath the sunk survey system and gets rid of the data already used to image structures above it. The deeper the depth is, the shorter is the valid time records of the remaining data and less computation is needed for one depth step continuation. For data decomposition, in addition to time axis, dreamlet survey sinking also decomposes the data for source and receiver gathers, which is a fully localized decomposition of prestack seismic data. A three‐scatter model is first used to demonstrate the computational feature and principle of this method. Tests on the two‐dimensional SEG/EAGE salt model show that with reduced data sets the proposed method can still obtain good imaging quality on complex geology structures and a strong velocity contrast environment.     

Hot Methanol Extraction for the Analysis of Volatile Organic Chemicals in Subsurface Core Samples from Dover Air Force Base, Delaware

The evaluation of contaminant concentrations in ground water and soil is an essential aspect of most hazardous waste remedial investigations. This paper describes methods applied toward obtaining, preserving, and analyzing subsurface samples for the determination of VOC concentrations in the saturated region of an unconfined coastal plain aquifer at Dover Air Force Base (DAFB). Delaware. The described protocol involved headspace-free subsampling of cores, field preservation of subsamples in methanol, and overnight extraction of the VOCs at elevated temperature (70°C). Methanol-extracted compounds were subsequently transferred to hexane and analyzed by gas chromatography. The method was found to achieve quantitative extraction from the aquifer sands in a single step, although extraction from fine-grained and more strongly sorbing aquitard samples required multiple methanol extractions to achieve comparable recovery. An extensive set of DAFB results is presented as an indication of how these methods can be applied toward characterizing field-scale contamination with a high degree of resolution and accuracy. These data have enabled a quantitative estimate of the contaminant distribution at this site and offer valuable insight into processes of contaminant migration into the underlying aquitard. In these regards, the methods used are believed to have provided much more accurate results than could have been obtained using more commonly applied techniques.     

Geostatistical modelling of ground conditions to support the assessment of site contamination

The pattern of ground contamination across a site depends on the historical pattern of contaminant releases at the surface and the redistribution and the fate of contaminants below the surface. Using these concepts a new site assessment approach (assessment modelling) is proposed based on the development of three stochastic models: a model of the physical structure of the ground materials beneath the site; a model of the distribution of surface contaminant spills; and a model of the flow and transport of spilled material into the heterogeneous underlying ground to construct alternative, equally likely, present day contaminant distributions. Combining the models within a Monte Carlo framework can, in principle, improve the understanding of the potential for excess contamination across the site and improve decisions on remediation options and locations. A trial application has been undertaken in the UK using a particular site to assess the approach. The conditions at the site used for the trial and the first of the stochastic model developments, the geostatistical modelling of the soil heterogeneity, are presented in this paper. Non-parametric and parametric geostatistics have been employed to formulate the geostatistical models of the site soils using lithological information from 146 trial pits and boreholes. The approach to the soil modelling and the verification and validation of the results are described. The heterogeneity of the subsurface is complicated by the presence of made-ground, comprised of various inert building wastes, and the non-stationarity of the heterogeneity of the natural ground. This paper is the first of three describing the assessment modelling methodology and its trial application to the site.     

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.     

An adaptive hybrid EnKF-OI scheme for efficient state-parameter estimation of reactive contaminant transport models

Reactive contaminant transport models are used by hydrologists to simulate and study the migration and fate of industrial waste in subsurface aquifers. Accurate transport modeling of such waste requires clear understanding of the system’s parameters, such as sorption and biodegradation. In this study, we present an efficient sequential data assimilation scheme that computes accurate estimates of aquifer contamination and spatially variable sorption coefficients. This assimilation scheme is based on a hybrid formulation of the ensemble Kalman filter (EnKF) and optimal interpolation (OI) in which solute concentration measurements are assimilated via a recursive dual estimation of sorption coefficients and contaminant state variables. This hybrid EnKF-OI scheme is used to mitigate background covariance limitations due to ensemble under-sampling and neglected model errors. Numerical experiments are conducted with a two-dimensional synthetic aquifer in which cobalt-60, a radioactive contaminant, is leached in a saturated heterogeneous clayey sandstone zone. Assimilation experiments are investigated under different settings and sources of model and observational errors. Simulation results demonstrate that the proposed hybrid EnKF-OI scheme successfully recovers both the contaminant and the sorption rate and reduces their uncertainties. Sensitivity analyses also suggest that the adaptive hybrid scheme remains effective with small ensembles, allowing to reduce the ensemble size by up to 80% with respect to the standard EnKF scheme.