Contaminant degradation in physically and chemically heterogeneous aquifers

This paper examines the importance of the correlation between hydraulic conductivity (K) and degradation rate constant (k) during the transport of reactive contaminants in heterogeneous aquifers. We simulated reactive transport in an ensemble of two-dimensional heterogeneous aquifers. Two sets of transport simulations were conducted: one in which a perfect positive correlation was assumed between ln(K) and ln(k), and one in which a perfect negative correlation was assumed. We found that the sign of the correlation has important consequences for the contaminant transport. Qualitatively, a negative correlation leads to significantly more pronounced "fingering" of the contaminant plume than does a positive correlation, with potentially important consequences for downgradient receptors. Quantitatively, the expected behavior (as quantified by the contaminant mass remaining in the aquifer) is statistically different between the positive and negative cases: on average, more contaminant mass persists when K and k are negatively correlated. Also, the negative correlation leads to more variability between realizations of the ensemble, whereas a positive correlation induces relatively little variability between realizations. We discuss the implications of these findings for the management of contaminated aquifers.