Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers: II. Model application to a field site

A multicomponent reactive transport model as presented by MacQuarrie and Sudicky [MacQuarrie, K.T.B., Sudicky, E.A., this volume. Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers: I. Model formulation and performance, J. Contam. Hydrol.] is applied to a well-studied wastewater plume in a sandy aquifer near Cambridge, Ontario. Domestic wastewater is released into the unsaturated zone via a drain field at a depth of about 0.8 m. The physical transport parameters for the model are obtained by simulating a non-reactive solute, while kinetic input data for the nitrogen and carbon reaction network are obtained from the literature. The model shows that the wastewater-loading rate has little influence on the moisture content in the unsaturated zone, thus oxygen diffusion in the air phase is an important transport mechanism. The model results are in general agreement with the field-determined moisture and oxygen profiles near the drain field. The simulation results show that oxidation of ammonium and dissolved organic carbon (DOC) goes to completion in the 1.5-m distance between the drain field and the water table, and that calcite dissolution limits the pH reduction to about 0.2 units. The model-predicted nitrate concentrations in the core of the plume are in the range of 20–25 mg N/l and are in good agreement with the field data. Overall, the results for the major reactive species from the model simulation agree well with the geochemical data obtained below the drain field and it is concluded that the major physical and biochemical processes have been correctly captured in the current model formulation.