Approximately 3,500 constructed wetlands (CWs) provide raw wastewater treatment in France for small communities (<5,000 people equivalent). Built during the past 30 years, most consist of two vertical flow constructed wetlands (VFCWs) in series (stages). Many configurations exist, with systems associated with horizontal flow filters or waste stabilization ponds, vertical flow with recirculation, partially saturated systems, etc. A database analyzed 10 years earlier on the classical French system summarized the global performances data. This paper provides a similar analysis of performance data from 415 full-scale two-stage VFCWs from an improved database expanded by monitoring data available from Irstea and the French technical department. Trends presented in the first study are confirmed, exhibiting high chemical oxygen demand (COD), total suspended solids (TSS) and total Kjeldahl nitrogen (TKN) removal rates (87%, 93% and 84%, respectively). Typical concentrations at the second-stage outlet are 74 mgCOD L−1, 17 mgTSS L−1 and 11 mgTKN L−1. Pollutant removal performances are summarized in relation to the loads applied at the first treatment stage. While COD and TSS removal rates remain stable over the range of applied loads, the spreading of TKN removal rates increases as applied loads increase.

Seven treatment wetlands and a municipal wastewater treatment plant (WWTP) were weekly monitored over the course of one year for removal of conventional wastewater parameters, selected micropollutants (caffeine, ibuprofen, naproxen, benzotriazole, diclofenac, acesulfame, and carbamazepine) and biological effects. The treatment wetland designs investigated include a horizontal subsurface flow (HF) wetland and a variety of wetlands with intensification (aeration, two-stages, or reciprocating flow). Complementary to the common approach of analyzing individual chemicals, in vitro bioassays can detect the toxicity of a mixture of known and unknown components given in a water sample. A panel of five in vitro cell-based reporter gene bioassays was selected to cover environmentally relevant endpoints (AhR: indicative of activation of the aryl hydrocarbon receptor; PPARγ: binding to the peroxisome proliferator-activated receptor gamma; ERα: activation of the estrogen receptor alpha; GR: activation of the glucocorticoid receptor; oxidative stress response). While carbamazepine was persistent in the intensified treatment wetlands, mean monthly mass removal of up to 51% was achieved in the HF wetland. The two-stage wetland system showed highest removal efficacy for all biological effects (91% to >99%). The removal efficacy for biological effects ranged from 56% to 77% for the HF wetland and 60% to 99% for the WWTP. Bioanalytical equivalent concentrations (BEQs) for AhR, PPARγ, and oxidative stress response were often below the recommended effect-based trigger (EBT) values for surface water, indicating the great benefit for using nature-based solutions for water treatment. Intensified treatment wetlands remove both individual micropollutants and mixture effects more efficiently than conventional (non-aerated) HF wetlands, and in some cases, the WWTP.

Horizontal subsurface Flow Constructed Wetlands (HF CWs) are biofilters planted with aquatic macrophytes within which wastewater is treated mostly through contact with bacterial biofilms. The high concentrations of organic carbon and nutrients being transported leads to high bacterial biomass production, which decreases the flow capacity of the porous material (bioclogging). In severe bioclogging scenarios, overland flow may take place, reducing overall treatment performance. In this work we developed a mathematical model using COMSOL Multiphysics™ and MATLAB(®) to simulate bioclogging effects in HF CWs. Variably saturated subsurface flow and overland flow were described using the Richards equation. To simplify the inherent complexity of the processes involved in bioclogging development, only one bacterial group was considered, and its growth was described using a Monod equation. Bioclogging effects on the hydrodynamics were taken into account by using a conceptual model that affects the value of Mualem's unsaturated relative permeability. Simulation results with and without bioclogging were compared to showcase the impact of this process on the overall functioning of CWs. The two scenarios rendered visually different bacteria distributions, flow and transport patterns, showing the necessity of including bioclogging effects on CWs models. This work represents one of the few studies available on bioclogging in variably saturated conditions, and the presented model allows simulating the interaction between overland and subsurface flow occurring in most HF CWs. Hence, this work gets us a step closer to being able to describe CWs functioning in an integrated way using mathematical models.