Evaluation of relative importance of different microbial reactions on organic matter removal in horizontal subsurface-flow constructed wetlands using a 2D simulation model

In this study, we used a two-dimensional (2D) mechanistic mathematical model in order to evaluate the relative contribution of different microbial reactions to organic matter removal (in terms of COD) in horizontal subsurface-flow constructed wetlands that treated urban wastewater. We also used the model to analyse the effect of increasing or decreasing the organic loading rate (changing the hydraulic loading rate (HLR) at a constant influent organic matter concentration, or changing the organic matter concentration at a constant HLR) on both the removal efficiency and the relative importance of the microbial reactions. The model is based on the code RetrasoCodeBright, which we modified to include the main microbial processes related to organic matter and nitrogen transformations in the wetlands: hydrolysis, aerobic respiration, nitrification, denitrification, sulphate reduction and methanogenesis. The model was calibrated and validated with data from two wetlands (each with a surface area of 55 m²) located in a pilot plant near Barcelona (Spain). According to the simulations, anaerobic processes (methanogenesis and sulphate reduction) are more widespread in the wetlands and contribute to a higher COD removal rate (60–70%) than anoxic (denitrification) and aerobic reactions do. These model results are confirmed by experimental observations. In all the cases tested, the reaction that most contributed to COD removal was methanogenesis (33–52%). According to our simulations, decreasing the HLR (for example, from 40 to 25 mm/d) while maintaining a constant COD influent concentration has a clear positive impact on COD removal efficiency (which increases from 65% to 89%). Changing influent COD concentration (for example, from 290 to 190 mg/L) while maintaining a constant HLR has a smaller impact, causing efficiency to increase from 79% to 84%. Changes in influent COD concentration (at a constant HLR) affect the relative contribution of the microbial reactions to organic matter removal. However, this trend is not seen when the HLR changes and the COD influent concentration remains constant.