Modelling nitrogen transformations in surface flow wastewater treatment wetlands in Sweden

The purpose of this study was to describe and compare the fate of nitrogen (N) in two Swedish wastewater treatment wetlands in the cities of Oxelösund and Hässleholm. Specifically, we wanted to see if a fairly simple model, developed with regard to common data availability, could satisfactorily describe the concentration dynamics at the outlet from the wetlands. A first-order area-based model, with two alternative expressions for temperature dependence, was set up to describe three major processes: ammonification, nitrification and denitrification. The N concentration dynamics at the outlet of the Oxelösund wetland was not satisfactorily described, R2(NH4+-N) = 0.33 and R2(NO3(-)-N) = 0.10, while the modelled concentrations corresponded quite well with measured concentrations in the Hässleholm wetland, R2(NH4+-N) = 0.83 and R2(NO3(-)-N) = 0.58. The NO3(-)-N concentrations, in both wetlands, could be slightly better described when introducing a temperature coefficient as an additional free parameter. The explained variances reported above were achieved when the model was calibrated individually for the two wetlands, when the resulting (optimised) reaction rate coefficients for each of the three processes were quite different between the two wetlands. To improve model performance, the rate equations may have to be changed to include factors in addition to concentration and temperature, such as dissolved oxygen and hydraulic efficiency. It may also be important to include other processes, such as plant uptake/decay and ammonia volatilisation.