Model evaluation of denitrification under rapid infiltration basin systems

Rapid Infiltration Basin Systems (RIBS) are used for disposing reclaimed wastewater into soil to achieve additional treatment before it recharges groundwater. Effluent from most new sequenced batch reactor wastewater treatment plants is completely nitrified, and denitrification (DNF) is the main reaction for N removal. To characterize effects of complex surface and subsurface flow patterns caused by non-uniform flooding on DNF, a coupled overland flow-vadose zone model is implemented in the multiphase flow and reactive transport simulator TOUGHREACT. DNF is simulated in two representative soils varying the application cycle, hydraulic loading rate, wastewater quality, water table depth, and subsurface heterogeneity. Simulations using the conventional specified flux boundary condition under-predict DNF by as much as 450% in sand and 230% in loamy sand compared to predictions from the coupled overland flow-vadose zone model, indicating that simulating coupled flow is critical for predicting DNF in cases where hydraulic loading rates are not sufficient to spread the wastewater over the whole basin. Smaller ratios of wetting to drying time and larger hydraulic loading rates result in greater water saturations, more anoxic conditions, and faster water transport in the vadose zone, leading to greater DNF. These results in combination with those from different water table depths explain why reported DNF varied with soil type and water table depth in previous field investigations. Across all simulations, cumulative percent DNF varies between 2 and 49%, indicating that NO3 removal in RIBS may vary widely depending on operational procedures and subsurface conditions. These modeling results improve understanding of DNF in RIBS and suggest operational procedures that may improve NO3 removal.