Estimating spatially-variable rate constants of denitrification in irrigated agricultural groundwater systems using an Ensemble Smoother

SummaryGroundwater contaminant transport models are being used increasingly to simulate the fate and transport of reactive solutes, particularly nitrate, within aquifer systems. These models, however, often are hindered by lack of information regarding parameters, such as kinetic decay rates, that govern the subsistence of the solute within the groundwater. In an overall effort to provide accurate estimates of spatially-variable parameters in numerical reactive transport modeling we employ a data assimilation scheme, the Ensemble Smoother (ES), which yields improved estimates of spatially-variable denitrification rates within an irrigated agricultural river-aquifer system using measurements of (i) nitrate concentration in the groundwater and (ii) mass of nitrate entering the river from the aquifer via groundwater flows. Based on the Kalman Filter methodology, in which distributed uncertain model results are corrected by assimilating measurement data from a reference system, the ES incorporates uncertain parameter values, associated model results, and measurement data into an update algorithm to provide an updated, corrected model state that approaches the reference system state. As an important step in eventually employing the methodology to real-world systems, this study evaluates the parameter estimation scheme for a synthetic aquifer system approaching hydrologic (heterogeneous hydraulic conductivity, cropping patterns, irrigation recharge, canal seepage) and chemical (denitrification, leaching concentrations) complexities expected in aquifers influenced by agricultural practices. Sensitivity analyses are conducted to investigate the influence of (i) the number of measurement data assimilated and (ii) the error assigned to the measurement data. Results indicate that the spatial distribution of denitrification rates can be estimated to a satisfying degree, and when implemented in additional model runs produce (i) simulated values that coincide favorably with measurement values from the reference state, and (ii) spatial distribution of nitrate concentration comparable to that of the reference state.

► Denitrification rate constants are estimated using an Ensemble Smoother. ► Nitrate concentrations in groundwater are used to estimate the rate constants. ► Nitrate mass discharged to surface water is used to estimate the rate constants. ► The level of conditioning depends on the error assigned the measurement data.