Estimating spatially-variable first-order rate constants in groundwater reactive transport systems

Numerical reactive transport models are often used as tools to assess aquifers contaminated with reactive groundwater solutes as well as investigating mitigation scenarios. The ability to accurately simulate the fate and transport of solutes, however, is often impeded by a lack of information regarding the parameters that define chemical reactions. In this study, we employ a steady-state Ensemble Kalman Filter (EnKF), a data assimilation algorithm, to provide improved estimates of a spatially-variable first-order rate constant lambda through assimilation of solute concentration measurement data into reactive transport simulation results. The methodology is applied in a steady-state, synthetic aquifer system in which a contaminant is leached to the saturated zone and undergoes first-order decay. Multiple sources of uncertainty are investigated, including hydraulic conductivity of the aquifer and the statistical parameters that define the spatial structure of the parameter field. For the latter scenario, an iterative method is employed to identify the statistical mean of lambda of the reference system. Results from all simulations show that the filter scheme is successful in conditioning the lambda ensemble to the reference lambda field. Sensitivity analyses demonstrate that the estimation of the lambda values is dependent on the number of concentration measurements assimilated, the locations from which the measurement data are collected, the error assigned to the measurement values, and the correlation length of the lambda fields. (C) 2010 Elsevier B.V. All rights reserved.