This study investigates the impact of aeration strategy on the performance of total nitrogen (TN) removal in a compact hybrid aerated treatment wetland (TW), called Rhizosph'air®. The system combines a single-stage French vertical flow wetland with an aerated horizontal-flow wetland, offering a unique and flexible approach for optimizing TN removal. In total, seven experimental conditions were tested, with different aeration modes, hydraulic loading rates and ammonium addition. The wetland system demonstrated high performance in terms of chemical oxygen demand removal (>85%) and solids removal (>90%), regardless of the experimental condition. However, TN removal was found to be directly impacted by operational changes. Increasing the hydraulic loading rate from 0.15 to 0.25 m/day led to an improvement in TN removal, achieving over 60%. Furthermore, when ammonium was added to the inlet and when the aeration timing was synced with the timing of the influent batch load, the environmental conditions facilitated the denitrification process, resulting in TN removal of approximately 70% and the lowest effluent NO3-N concentrations (8.70 ± 4.40 mg/L). In summary, the timing of the aeration strategy according to influent batch loading improved TN removal, suggesting its potential for optimization in future studies.

  • This study explored different experimental conditions to understand the impact of aeration on the dynamics of nitrogen removal in a hybrid aerated treatment wetland.

  • The total nitrogen removal is limited when the aeration is controlled according to time.

  • The transition from time-based to batch-hour aeration strategy increased the denitrification process, resulting in 70% removal of total nitrogen.

Treatment wetlands (TWs) have gained considerable attention for wastewater treatment due to their robustness, low operation and maintenance requirements and ability to deliver high-quality effluent (Morvannou et al. 2015). Total nitrogen (TN) in wastewater is usually removed through the use of a traditional nitrification–denitrification pathway. The process requires two steps: aerobic autotrophic nitrification to transform ammonia into nitrates and anoxic heterotrophic denitrification, to transform the nitrates into nitrogen gas. Both phases demand specific environmental conditions, foremost the oxygen concentration and organic carbon source (Saeed & Sun 2012). Passive TWs usually require two or three treatment stages to address TN removal by implementing separate aerobic and anoxic filters. Therefore, depending on the design chosen, the land requirement can vary between 4 and 10 m²/PE, which is a disadvantage for those systems (Dotro et al. 2017) compared to conventional systems (0.12–0.3 m²/PE for activated sludge).

In addition, passive TWs have some limitations regarding pollutant removal; namely, horizontal-flow (HF) wetlands do not nitrify to a great extent, and vertical flow wetlands nitrify but are limited in their capacity to remove TN (Masih 2017; Zhuang et al. 2019; Cross et al. 2021). While vertical flow (VF) TWs exhibit almost complete ammonium removal (upwards of 90%), passive TWs (either vertical or horizontal) achieve only 20–50% TN removal (Cross et al. 2021). Consequently, considering the restrictions imposed on nutrient discharge and the problem of available land area in urban areas, intensified TW designs have been developed. The term ‘intensification’ implies the use of energetic, chemical or operational modifications to increase pollutant removal (Fonder & Headley 2010). By adopting an energy-oriented approach, it becomes possible to enhance the availability of oxygen within a TW, thereby promoting the removal of oxygen-requiring compounds such as carbon and ammonium.

Aeration in TW refers to the injection of air into the system, enabling aerobic conditions in a saturated zone. Aerated treatment wetland (TW) technology was created over 25 years ago (Wallace 2001), and its implementation has witnessed a notable surge over the past decade. Today, approximately 500 aerated TW systems are operational on a global scale (Nivala et al. 2020). The main advantage of aerated TWs is their ability to completely nitrify (Masih 2017; Zhuang et al. 2019; Cross et al. 2021), and their compactness compared to passive TWs (1 m2/population equivalent (PE) versus 4–5 m2/PE) (DWA 2017); consequently, rendering them more viable for implementation where land availability is limited, such as in urban areas. Despite the disadvantage of increased energy costs, those systems operate at low pressure and consume less energy than conventional intensive wastewater treatment plants (Rochard 2017).

Moreover, intermittent or controlled aeration strategies (compared to continuous 24/7 aeration) have demonstrated the potential to enhance denitrification processes in aerated wetlands (Dotro et al. 2017). Continuous aeration has been shown to achieve TN concentration removal ranging from 15 to 60% (Nivala et al. 2019; Cross et al. 2021) and intermittent aeration has been found to optimize TN removal, yielding removal efficiencies of up to 90% (Masih 2017).

Recent configurations and adaptations have been mixing the aeration and hybrid TW to go even further in the TN removal. One of these systems, known as Rhizosph'air® (Troesch et al. 2020), allows the nitrification and denitrification processes in a single-stage TW. The aeration in this system provides increased oxygen transfer, resulting in the reduction of the system footprint from 2 m2 PE−1 (classical two-stage French vertical flow wetland) to 0.8–1.2 m2 PE−1 depending on the conditions (Dou et al. 2016; Cross et al. 2021). This hybrid system has been installed around France (18 plants) and one of them was monitored for 2 years, treating domestic wastewater for 1400 PE (Prost-Boucle & Molle 2021). Regardless of the tested aeration mode, the treatment performances for COD, BOD5 and TSS consistently remained high and stable. However, under the aeration mode of 12 h per day divided into four cycles, nitrification was successfully achieved, while denitrification was limited due to insufficient carbon availability. To enhance TN removal, a reduction in the daily aeration duration proved beneficial. By adopting a four-cycle aeration mode with a total duration of 3 h per day, the system was capable of removing 97% COD concentration and 82% TN concentration, remaining 19 mg TN/L at the outlet.

However, they tested limited aeration strategies with high variability in inflow and outflow. This current research intends to go further by testing the difference between each experimental condition in a controlled system, allowing a comparison to be made. Thus, to guarantee a stable and optimized nitrogen removal, it is crucial to adapt aeration strategies to meet specific treatment goals, ranging from complete nitrification to nearly complete removal of TN. Aeration strategies should be based on carbon and nitrification requirements, considering the absence of oxygen for denitrification. Therefore, this research aims to investigate the impact of aeration strategy in a hybrid aerated vertical/horizontal flows TW on the performance of TN removal.

TW pilot

The industrial pilot-scale TW studied is part of the REFLET research platform (https://eng-reversaal.lyon-grenoble.hub.inrae.fr/equipment-platforms/reflet), located in Craponne, France, a commune situated in the Auvergne Rhône Alpes region. The pilot system had 1 month of startup before starting the sampling campaign. Its commission started in January 2019. In total, the pilot covers an area of 20 m2 (8 × 2.5 m2). It replicates the Rhizosph'air system, which comprises a single-stage French vertical flow TW, receiving raw wastewater directly on the surface of the filter, above an aerated horizontal-flow TW. This arrangement enables the TW system to provide combined primary and secondary wastewater treatment.

The TW receives wastewater from around 4 m3/day (26 PE), collected by a combined sewer from the upstream sewage basin produced by the community. The domestic wastewater has a ratio of COD/BOD5 of approximately 2.0, showing that is a typical municipal untreated wastewater (Kadlec & Wallace 2009).

The TW is planted with Phragmites australis and the upper layer is divided into two filters (10 m2 each) by a PVC plate embedded up to the level of the filtration (0.2 m depth) and transition layers (0.15 m depth). While plants can impact nutrient transformation, sequestration and uptake (Val del Río et al. 2017) the main role of plants in this system is related to their mechanical role to avoid surface clogging. Indeed, nutrient uptake is minor compared to the load applied. This division allows for alternating feeding/resting phases of the filters (3.5 days each) to ensure mineralization of the organic deposit layer. The TW is fed by batches, and the filter that actively receives wastewater is referred to as the ‘primary filter (PF)’, while the filter at rest is referred to as the ‘secondary filter (SF)’. The drainage system of the pilot system is designed such that during the feeding phase of the PF, the outflow is directed to the outlet of the SF, which is achieved through the opening and closing of associated electro-valves, as depicted in Figure 1.
Figure 1

Hydraulic operating scheme of the hybrid aerated vertical/horizontal-flow treatment wetland pilot.

Figure 1

Hydraulic operating scheme of the hybrid aerated vertical/horizontal-flow treatment wetland pilot.

Close modal

The saturation layer (1 m depth) is shared between the two parallel filters. Aeration of the system is accomplished by employing 120 driplines positioned at the bottom of the saturation layer (60 diffuser/m2). The aeration process is facilitated by two side channel blowers (Becker SV 5.90/1). The estimated supply airflow rate for each blower is approximately 15 m3/h (equivalent to 1.5 m3/h/m2). It is worth noting that no gas flow meter or temperature sensor was installed and therefore, the airflow rate estimation is based on these specifications.

Pilot operation and experimental conditions

In order to evaluate the influence of aeration timing on nitrogen removal dynamics in the hybrid aerated TW, seven experimental conditions were tested. Table 1 and Figure 2 present the experimental condition parameters. The nomenclature for each experimental condition follows the order: T or F (aeration based on time or feeding), sum total aeration time per day (the same aeration for both filters or differentiating aeration for each filter, in h/day), hydraulic loading rate (HLR) on the PF (in m/day), the addition or not of ammonium (NH4).
Table 1

Experimental conditions parameters

Experimental condition
C1C2C3C4C5C6C7
Ammonium addition no no no no Yes yes yes 
Aeration based on timer (T) or feeding (F) 
Aeration on the primary filter Number of cycles/day 10 10 
Operation (min per cycle) 90 45 45 45 45 20 15 
Stop (min per cycle) 150 315 315 315 315 130 135 
Aeration on the secondary filter Number of cycles/day 10 
Operation (min per cycle) 90 45 45 45 45 10 
Stop (min per cycle) 150 315 315 315 315 140 
Hydraulic loading rate m/day 0.15 0.15 0.25 0.35 0.35 0.35 0.35 
Sampling period 1/19–5/19 5/19–7/19 07/19–10/19 10/19–01/20 05/20–09/20 06/21–08/21 08/21–10/21 
Duration (days) 129 73 84 113 119 41 62 
Number of samples 10 
Experimental condition
C1C2C3C4C5C6C7
Ammonium addition no no no no Yes yes yes 
Aeration based on timer (T) or feeding (F) 
Aeration on the primary filter Number of cycles/day 10 10 
Operation (min per cycle) 90 45 45 45 45 20 15 
Stop (min per cycle) 150 315 315 315 315 130 135 
Aeration on the secondary filter Number of cycles/day 10 
Operation (min per cycle) 90 45 45 45 45 10 
Stop (min per cycle) 150 315 315 315 315 140 
Hydraulic loading rate m/day 0.15 0.15 0.25 0.35 0.35 0.35 0.35 
Sampling period 1/19–5/19 5/19–7/19 07/19–10/19 10/19–01/20 05/20–09/20 06/21–08/21 08/21–10/21 
Duration (days) 129 73 84 113 119 41 62 
Number of samples 10 
Figure 2

Parameters for each experimental condition.

Figure 2

Parameters for each experimental condition.

Close modal

Aeration was based on time (‘T’ condition) for the first five experiments, e.g., hours of aeration per day spread across a number of aeration cycles. Initially, a low HLR (0.15 m/day) with the longest aeration time (9 h) was tested as a C1 condition to favor plant establishment as well as good aerobic performances: T–9 h–0.15. The aim of this condition was to not overload the system from the beginning and fasten nitrifying bacteria growth in winter. Then, the aeration time and the number of aeration cycles were decreased (C2: T–3 h–0.15). To better identify the treatment capacity of the system, the HLR was gradually increased (C3: T–3 h–0.25 and C4: T–3 h–0.35), and then the pollutant load was increased by adding ammonia to the influent (C5: T–3 h–0.35–NH4). Since the influent comes from a combined sewer and, consequently, the inlet ammonium concentration was low, ammonium addition (the goal was to reach around 80 mg NH4-N/L) allowed us to test the capacity of the system for more stringent wastewater concentrations (C5, C6 and C7 conditions).

To increase denitrification, the timing of the aeration was coordinated with pulse feeding on the PF (‘F’ condition), optimizing the use of the available carbon in raw wastewater. Additionally, the aeration was decreased on the SF (C6: F–3.3 h–1.7 h–0.35–NH4). For the last condition (C7: F–2.5 h–0 h–0.35–NH4) aeration was decreased on the PF (2.5 h) and completely stopped on an SF (0 h) to favor denitrification. In these cases, aeration was set to start 40 min before the batch and stop after 20 min for C6 and 15 min for C7 (time needed to degrade the oxygen inside the filter and remain in an anoxic environment).

Pilot monitoring – chemical parameters

The pilot operation was monitored using a combination of continuous online sensors and laboratory analysis. Eight sensors were strategically installed along the filter, with measurements recorded at a 5-min time step (Figure 3). Three of these sensors were dedicated to monitoring NH4-N and NO3-N concentrations (WTW, FDO 700IQ), positioned at the inlet, center and outlet of the pilot. The remaining sensors were oxygen sensors (WTW, VARION 700IQ), distributed within horizontal and vertical piezometers located inside the pilot.
Figure 3

Position and type of sensor at the inlet, inside and outlet of the pilot, and location of the 24-h composite sampling. The sensors were strategically placed to have a different vertical and horizontal profile and to have a view of the processes taking place in each filter and in half of the pathway.

Figure 3

Position and type of sensor at the inlet, inside and outlet of the pilot, and location of the 24-h composite sampling. The sensors were strategically placed to have a different vertical and horizontal profile and to have a view of the processes taking place in each filter and in half of the pathway.

Close modal

The laboratory analysis occurred once per week, sampling the inlet and outlet through 24-h composite samples. Each week, the filter was sampled alternated (different inlet and outlet). The following parameters were analyzed: total chemical oxygen demand (COD), filtered COD (fCOD), total suspended solids (TSS), NH4-N, NO3-N, NO2-N and Total Kjeldahl Nitrogen (TKN). The chemistry laboratory followed the European standard methods (NF EN 872, EN ISO 10304-1, NF EN ISO 14911, NF EN 25663, NF EN 6878, NF EN ISO 9963-1). TN was calculated as the sum of NO3-N, NO2-N and TKN concentrations.

A statistical analysis was conducted to compare the treatment performance of each filter since they are fed alternately and heterogeneity could change the removal processes. This comparison was performed using a linear mixed model, which is suitable for analyzing grouped observations (Zuur et al. 2009). In our case, we considered each experimental condition as a random effect and the use of a filter on the left or right side as a fixed effect. The following equations describe such a relationship:
(1)
(2)
(3)
where: i is the observation, j is the experimental condition, a is the intercept term, is the random intercept term, b is the fixed effect, is a dummy variable indicating the filter used (0 = left filter; 1 = right filter), is the residual term, N is the normal distribution and and are the residual and random effect variance, respectively. If the results showed a low effect between the filter results, it means that they were homogeneous and could be processed together. The statistical analysis for all experimental conditions (refer to the Supplementary Material for detailed results) revealed a statistically insignificant difference between filters, indicating their homogeneity. Consequently, the results were treated collectively as representative of a single treatment system.

Treatment performances and limitations: timing aeration strategy

Table 2 presents the treatment performance for the experimental conditions employing the timing aeration strategies (‘T’). Notably, the concentrations measured at the inlet exhibited variation across different conditions, evidenced by lower COD concentrations for C1 (T–9 h–0.15) (384 mg COD/L) and C4 (T–3 h–0.35) (280 mg COD/L), emphasizing the influence of seasonal factors. These periods corresponded to the winter season, characterized by wastewater dilution due to rainwater inflow and intrusion of a greater quantity of clear water in the combined sewer. Overall, the system demonstrated high performance in removing carbon and solids pollution regardless of the experimental conditions, with an average outlet concentration of less than 29 mg COD/L and 4 mg TSS/L. These results highlight the filtration capacity provided by the TW. Moreover, these findings align with previous studies on the treatment efficiency of similar full-scale systems (Petitjean et al. 2021), affirming the effectiveness and stability of Rhizosph'air® in removing carbon and solids from raw wastewater.

Table 2

Inlet, outlet concentrations and concentration removal for each condition, considering temperature inside the filter, COD, TSS, TN, NH4-N and NO3-N, where Avg denotes the average, Std dev denotes the standard deviation and n denotes the number of observations

In
Out
Concentration percent removal (%)
MedianAvg + Std devnMedianAvg + Std devnMedianAvg + Std devn
C1: T–9 h–0.15 (Temperature inside the filter: 8.4 ± 2.5 °C) 
COD (mg/L) 325 384 ± 197 10 21 24 ± 6 10 93% 92 ± 4% 10 
TSS (mg/L) 244 289 ± 146 10 <2 <2 ± 0 10 99% 99 ± 1% 10 
TN (mg/L) 36.1 33.9 ± 7.7 10 35.7 34.4 ± 5.7 10 5.5% 5.0 ± 12.7% 
NH4-N (mg/L) 34.7 32.3 ± 9.3 10 0.10 0.23 ± 0.48 10 99.7% 99.2 ± 1.7% 10 
NO3-N (mg/L) 0.83 1.33 ± 1.79 10 35.6 34.0 ± 5.9 10 – – – 
C2: T–3 h–0.15 (Temperature inside the filter: 17.2 ± 2.7 °C) 
Temperature (°C)          
COD (mg/L) 642 723 ± 202.0 25 26 ± 2 96% 96 ± 1% 
TSS (mg/L) 570 585 ± 146 <2 <2 ± 1 99.6% 99.6 ± 0.2% 
TN (mg/L) 51.7 55.8 ± 19.4 24.1 26.4 ± 5.7 45.1% 47.2 ± 19.6% 
NH4-N (mg/L) 45.5 43.0 ± 8.5 0.02 0.23 ± 0.31 99.9% 99.5 ± 0.7% 
NO3-N (mg/L) 0.51 0.71 ± 0.43 22.6 25.3 ± 6.2 – – – 
C3: T–3 h–0.25 (Temperature inside the filter: 19.9 ± 1.4 °C) 
Temperature (°C)          
COD (mg/L) 554 545 ± 212 25 24 ± 4 96% 95 ± 3% 
TSS (mg/L) 380 434 ± 275 <2 <2 ± 0 99.5% 99.4 ± 0.3% 
TN (mg/L) 35.3 41.1 ± 21.9 16.8 17.1 ± 1.8 66.4% 62.6 ± 17.7% 
NH4-N (mg/L) 38.3 36.2 ± 16.2 0.47 0.65 ± 0.62 98.4% 97.6 ± 3.0% 
NO3-N (mg/L) 0.49 0.75 ± 0.71 16.2 16.0 ± 2.4 – – – 
C4: T–3 h–0.35 (Temperature inside the filter: 12.0 ± 3.2 °C) 
Temperature (°C)          
COD (mg/L) 274 280 ± 112 <20 <20 ± 1 93% 92 ± 4% 
TSS (mg/L) 176 179 ± 91 <2 <2 ± 0 99% 99 ± 1% 
TN (mg/L) 35.9 29.8 ± 12.2 15.4 15.0 ± 1.1 56.0% 42.8 ± 25.7% 
NH4-N (mg/L) 23.3 23.9 ± 11.1 0.24 0.65 ± 0.91 98.2% 97.8 ± 2.4% 
NO3-N (mg/L) 0.74 1.17 ± 0.9 13.6 13.6 ± 0.2 – – – 
C5: T–3 h–0.35–NH4 (Temperature inside the filter: 19.8 ± 1.7 °C) 
Temperature (°C)          
COD (mg/L) 772 789 ± 242 28 29 ± 6 97% 96 ± 1% 
TSS (mg/L) 446 494 ± 145 3.7 3.9 ± 1.1 99% 99 ± 0.3% 
TN (mg/L) 106.4 108.2 ± 18.1 28.1 29.9 ± 11.1 75.2% 71.2 ± 12.8% 
NH4-N (mg/L) 83.8 86.9 ± 19.0 6.09 8.17 ± 8.92 93.7% 89.2 ± 12.5% 
NO3-N (mg/L) 0.45 0.48 ± 0.33 18.0 17.6 ± 4.9 – – – 
In
Out
Concentration percent removal (%)
MedianAvg + Std devnMedianAvg + Std devnMedianAvg + Std devn
C1: T–9 h–0.15 (Temperature inside the filter: 8.4 ± 2.5 °C) 
COD (mg/L) 325 384 ± 197 10 21 24 ± 6 10 93% 92 ± 4% 10 
TSS (mg/L) 244 289 ± 146 10 <2 <2 ± 0 10 99% 99 ± 1% 10 
TN (mg/L) 36.1 33.9 ± 7.7 10 35.7 34.4 ± 5.7 10 5.5% 5.0 ± 12.7% 
NH4-N (mg/L) 34.7 32.3 ± 9.3 10 0.10 0.23 ± 0.48 10 99.7% 99.2 ± 1.7% 10 
NO3-N (mg/L) 0.83 1.33 ± 1.79 10 35.6 34.0 ± 5.9 10 – – – 
C2: T–3 h–0.15 (Temperature inside the filter: 17.2 ± 2.7 °C) 
Temperature (°C)          
COD (mg/L) 642 723 ± 202.0 25 26 ± 2 96% 96 ± 1% 
TSS (mg/L) 570 585 ± 146 <2 <2 ± 1 99.6% 99.6 ± 0.2% 
TN (mg/L) 51.7 55.8 ± 19.4 24.1 26.4 ± 5.7 45.1% 47.2 ± 19.6% 
NH4-N (mg/L) 45.5 43.0 ± 8.5 0.02 0.23 ± 0.31 99.9% 99.5 ± 0.7% 
NO3-N (mg/L) 0.51 0.71 ± 0.43 22.6 25.3 ± 6.2 – – – 
C3: T–3 h–0.25 (Temperature inside the filter: 19.9 ± 1.4 °C) 
Temperature (°C)          
COD (mg/L) 554 545 ± 212 25 24 ± 4 96% 95 ± 3% 
TSS (mg/L) 380 434 ± 275 <2 <2 ± 0 99.5% 99.4 ± 0.3% 
TN (mg/L) 35.3 41.1 ± 21.9 16.8 17.1 ± 1.8 66.4% 62.6 ± 17.7% 
NH4-N (mg/L) 38.3 36.2 ± 16.2 0.47 0.65 ± 0.62 98.4% 97.6 ± 3.0% 
NO3-N (mg/L) 0.49 0.75 ± 0.71 16.2 16.0 ± 2.4 – – – 
C4: T–3 h–0.35 (Temperature inside the filter: 12.0 ± 3.2 °C) 
Temperature (°C)          
COD (mg/L) 274 280 ± 112 <20 <20 ± 1 93% 92 ± 4% 
TSS (mg/L) 176 179 ± 91 <2 <2 ± 0 99% 99 ± 1% 
TN (mg/L) 35.9 29.8 ± 12.2 15.4 15.0 ± 1.1 56.0% 42.8 ± 25.7% 
NH4-N (mg/L) 23.3 23.9 ± 11.1 0.24 0.65 ± 0.91 98.2% 97.8 ± 2.4% 
NO3-N (mg/L) 0.74 1.17 ± 0.9 13.6 13.6 ± 0.2 – – – 
C5: T–3 h–0.35–NH4 (Temperature inside the filter: 19.8 ± 1.7 °C) 
Temperature (°C)          
COD (mg/L) 772 789 ± 242 28 29 ± 6 97% 96 ± 1% 
TSS (mg/L) 446 494 ± 145 3.7 3.9 ± 1.1 99% 99 ± 0.3% 
TN (mg/L) 106.4 108.2 ± 18.1 28.1 29.9 ± 11.1 75.2% 71.2 ± 12.8% 
NH4-N (mg/L) 83.8 86.9 ± 19.0 6.09 8.17 ± 8.92 93.7% 89.2 ± 12.5% 
NO3-N (mg/L) 0.45 0.48 ± 0.33 18.0 17.6 ± 4.9 – – – 

Regarding ammonium removal, the nitrification process was nearly complete, regardless of the total aeration duration per day (C1: 9 h/day or C2, C3, C4: 3 h/day). The increase of HLR (C2: T–3 h–0.15 and C3: T–3 h–0.25) did not have a significant impact on the nitrification process, since the treated loads were approximately equal to the applied loads (Figure 4(a). The low output NH4-N concentrations (<1 mg NH4-N/L) and associated treated loads (<18 g NH4-N/m2/day) for conditions C1 (T–9 h–0.15), C2 (T–3 h–0.15), C3 (T–3 h–0.25) and C4 (T–3 h–0.35) demonstrated stable ammonium removal efficiency. These values fall within the lower ranges of reference values for the first stage of French VF TWs (Morvannou et al. 2015), with approximately 99% ammonium removal and production of nitrate.
Figure 4

Applied and treated loads for NH4-N (A) and TN (B), for experimental conditions C1 (T–9 h–0.15), C2 (T–3 h–0.15), C3 (T–3 h–0.25), C4 (T–3 h–0.35) and C5 (T–3 h–0.35–NH4).

Figure 4

Applied and treated loads for NH4-N (A) and TN (B), for experimental conditions C1 (T–9 h–0.15), C2 (T–3 h–0.15), C3 (T–3 h–0.25), C4 (T–3 h–0.35) and C5 (T–3 h–0.35–NH4).

Close modal

However, when ammonium was added to the influent raw wastewater (C5: T–3 h–0.35–NH4), removal efficiency decreased to 89.2 ± 12.5%, indicating that nitrification had reached its maximum. To further enhance ammonium removal, the possible strategies are to increase the number of aeration cycles per day, the duration of aeration for each cycle or to automate the aeration using online sensors. Nevertheless, the result reinforces previous observations about the ability of a hybrid aerated TW to achieve ammonium removal over 90% (Prost-Boucle & Molle 2021).

In classical wastewater treatment process engineering, the nitrification process is the first step in TN removal. The challenge of TN removal is the second part, the denitrification process (Saeed & Sun 2012). In the case of condition C1 (T–9 h–0.075), the issue for initiating denitrification was excessive oxygenation within the filter, where the aerobic periods lasted 24 h per day, as shown in Figure 5. The high dissolved oxygen concentration (in average 10 mg/L) led to the complete oxidation of ammonium (99.2 ± 1.7% removal), resulting in the reduction of ammonium nitrogen into nitrate and minimal removal of TN (concentration removal of 5.0 ± 12.7%, Table 2). Uggetti et al. (2016) tested a continuously aerated TW and found nearly complete nitrification, with 7–8 mg DO/L at the outlet.
Figure 5

Duration in hours per day of aerobic condition (>0.2 mg O2/L) and anoxic condition (<0.2 mg O2/L) for the experimental conditions employing the time-based aeration strategy. A single oxygen sensor is exhibited for a period of 1 month. The red dash lines represent the filter alternation, PF stands for primary filter and SF for secondary filter.

Figure 5

Duration in hours per day of aerobic condition (>0.2 mg O2/L) and anoxic condition (<0.2 mg O2/L) for the experimental conditions employing the time-based aeration strategy. A single oxygen sensor is exhibited for a period of 1 month. The red dash lines represent the filter alternation, PF stands for primary filter and SF for secondary filter.

Close modal

Comparing the C1 (T–9 h–0.15) and C2 (T–3 h–0.15) conditions (6 cycles versus 4 cycles of aeration/day), the C2 condition showed significantly higher denitrification, as indicated by the lower average outlet nitrate concentration (C1: 34 mg NO3-N/L; C2: 25 mg NO3-N/L), resulting in higher TN percent concentration removal for C2 of 47%. This can be attributed to the longer duration of the anoxic condition, with an average of 13.1 h per day in the PF and 3.2 h/day in the SF, which favored denitrifying bacteria activity. In Figure 5, for condition C2, whenever the filter is being fed, there is more carbon to degrade and the oxygen is rapidly consumed, while in the SF, the oxygen concentration increases, exhibiting longer aerobic conditions (3.2 h/day for the PF and 20.5 h/day for the SF).

When a similar aeration strategy (3 h/day) was applied but the HLR was increased from 0.15 m/day (C2) to 0.25 m/day (C3) and 0.35 m/day (C4), there was a decrease in the average effluent NO3-N concentration (C2: 25.3 ± 6.2 mg/L; C3: 16.0 ± 2.4 mg/L; C4: 13.6 ± 0.2 mg/L). Consequently, TN removal increased from C2 (47.2 ± 19.6%) to C3 (62.6 ± 17.7%). The lower percent concentration TN removal in C4 (T–3 h–0.35) (42.8 ± 25.7%) was attributed to the lower inlet TN concentration and lower hydraulic retention time as the HLR was increased. During conditions C2 and C3, denitrification was more efficient as the HLR were higher and higher inlet COD concentrations (about twice as high) provided higher availability of organic carbon. In addition, for these three experimental conditions, the loads applied in NH4-N (2.3–18 g/m2/day) and TN (2.4–21.7 g/m2/day) remained within the same range of values, and the treated loads (2.1–17.9 g NH4-N/m2/day and 0.7–17.7 g TN/m2/day) did not differ among the experimental conditions (Figure 4). These treated loads are higher than aerated VF (5.2 g NH4-N/m2/day) or aerated HF (7.3 g NH4-N/m2/day) (Dotro et al. 2017).

For condition C5, the condition with an ammonium addition (C5: T–3 h–0.35–NH4), the inlet average concentration was 108 ± 18 mg TN/L (between 9 and 19 g TN/m2/day of applied load), approximately 2–4 times higher than the inlet TN concentration in other conditions. This operational condition achieved a high TN removal, reaching 71%. However, effluent NO3-N concentrations were still measurable, indicating limitations in the denitrification process. The oxygen sensors reported an average of 18.1 h/day of anoxic conditions in the PF and 9.4 h/day in the SF. As Figure 5 shows, it is the first time that the duration of an anoxic condition was constantly longer than aerobic condition. Therefore, with high nitrification and a long anoxic time during this experimental condition, it became clear that the limitation of denitrification is not due to a lack of nitrates or excessive oxygenation in the saturated zone, but probably to the lack of organic carbon in the saturated zone. This issue is commonly reported in the literature (Bezbaruah & Zhang 2003; Wang et al. 2022; Yao et al. 2023).

Based on the assumptions of Morvannou et al. (2015) (80% concentration removal for COD and 50% for TKN in the unsaturated zone of French VF TWs), the COD/NO3-N ratio arriving in the saturated zone can be estimated. According to Morvannou et al. (2017), this ratio should be at least 2.0 at the entrance to the saturated zone (and a minimum residence time of 0.75 days) for efficient denitrification. Estimations in this study would indicate that only condition C5 (T–3 h–3.5–NH4) lacked sufficient biodegradable carbon for denitrification. The deficit would be approximately between 40 and 120 mg COD/L (results are presented in Supplementary material).

It was not possible to match aeration requirements to batch feeding by controlling aeration strategies on a temporal basis. When a batch arrived during a period of high dissolved oxygen concentration in the saturated zone, oxygen was used for carbon degradation resulting in a lack of carbon for denitrification during anoxic periods. Under the conditions tested from C1 to C5, the balance between nitrification and denitrification periods was not ideal for effectively transforming nitrogen forms while maintaining a sufficient carbon source for the denitrification process. Thus, to remedy this deficiency and attempt to optimize denitrification, a new aeration strategy was tested that coordinated the aeration period and batch feeding.

Treatment performance and limitations: aeration strategy according to batch feeding periods

Table 3 presents the treatment performance results for the experimental conditions employing the aeration strategy according to feeding time (‘F’). The first condition C6 (F–3.3 h–1.7 h–0.175–NH4) had similar influent concentrations to condition C4 (T–3 h–0.175), characterized by diluted wastewater and subsequently, lower COD and NH4-N influent concentrations. Despite this dilution, the system demonstrated high efficiency in removing those pollutants, even in wet weather conditions (refer to the Supplementary material for the local rainfall data), demonstrating some flexibility in the performance of the system.

Table 3

Inlet, outlet concentrations and concentration removal for each condition, considering COD, TSS, TN, NH4-N and NO3-N, where Avg denotes the average, Std dev denotes the standard deviation and n denote the number of observations

In
Out
Concentration removal (%)
MedianAvg + Std devnMedianAvg + Std devnMedianAvg + Std devn
C6: F-3.3 h–1.7 h–0.35–NH4 (Temperature inside the filter: 19.5 ± 0.5 °C) 
COD (mg/L) 286 252 ± 136 21 23 ± 4 91% 87 ± 10% 
TSS (mg/L) 201 189 ± 113 <2.0 <2.0 ± 0.0 99% 98 ± 2% 
TN (mg/L) 53.3 59.6 ± 42.8 31.5 35.3 ± 12.0 38.2% 24.7 ± 9.2% 
NH4-N (mg/L) 52.7 47.1 ± 20.9 0.08 1.35 ± 3.01 99.6% 98.0 ± 4.4% 
NO3-N (mg/L) 0.45 0.68 ± 0.55 30.5 30.3 ± 3.1 – – – 
C7: F-2.5 h-0 h-0.35-NH4 (Temperature inside the filter: 20.2 ± 0.5 °C) 
COD (mg/L) 595 548 ± 137 40 34 ± 13 94% 94 ± 2% 
TSS (mg/L) 367 384 ± 85 7 ± 1 98% 98 ± 0% 
TN (mg/L) 98.0 91.8 ± 41.6 29.6 27.0 ± 10.2 70.5% 70.7 ± 5.3% 
NH4-N (mg/L) 97.3 91.3 ± 41.4 19.5 18.3 ± 7.1 80.5% 80.0 ± 5.3% 
NO3-N (mg/L) 0.45 0.38 ± 0.16 9.76 8.70 ± 4.40 – – – 
In
Out
Concentration removal (%)
MedianAvg + Std devnMedianAvg + Std devnMedianAvg + Std devn
C6: F-3.3 h–1.7 h–0.35–NH4 (Temperature inside the filter: 19.5 ± 0.5 °C) 
COD (mg/L) 286 252 ± 136 21 23 ± 4 91% 87 ± 10% 
TSS (mg/L) 201 189 ± 113 <2.0 <2.0 ± 0.0 99% 98 ± 2% 
TN (mg/L) 53.3 59.6 ± 42.8 31.5 35.3 ± 12.0 38.2% 24.7 ± 9.2% 
NH4-N (mg/L) 52.7 47.1 ± 20.9 0.08 1.35 ± 3.01 99.6% 98.0 ± 4.4% 
NO3-N (mg/L) 0.45 0.68 ± 0.55 30.5 30.3 ± 3.1 – – – 
C7: F-2.5 h-0 h-0.35-NH4 (Temperature inside the filter: 20.2 ± 0.5 °C) 
COD (mg/L) 595 548 ± 137 40 34 ± 13 94% 94 ± 2% 
TSS (mg/L) 367 384 ± 85 7 ± 1 98% 98 ± 0% 
TN (mg/L) 98.0 91.8 ± 41.6 29.6 27.0 ± 10.2 70.5% 70.7 ± 5.3% 
NH4-N (mg/L) 97.3 91.3 ± 41.4 19.5 18.3 ± 7.1 80.5% 80.0 ± 5.3% 
NO3-N (mg/L) 0.45 0.38 ± 0.16 9.76 8.70 ± 4.40 – – – 

It is noteworthy that NH4-N concentrations at the outlet were 20 times higher for condition C7 (F–2.5 h–0 h–0.175–NH4) than for condition C6 (F–3.3 h–1.7 h–0.175–NH4). However, condition C7 also exhibited the lowest NO3-N outlet concentration among all experimental conditions (8.70 ± 4.40 mg NO3-N/L), consequently resulting in the highest TN removal. This condition demonstrated its capability to effectively treat TN loads that were nearly twice as high as those of the other experimental conditions (as depicted in Figure 6). These applied and treated loads were much higher than the previous studies with aerated TW (Zhang et al. 2010; Li et al. 2014; Uggetti et al. 2016).
Figure 6

Applied and treated loads for NH4-N and TN, for experimental conditions C1 (T–9 h–0.15), C2 (T–3 h–0.15), C3 (T–3 h–0.25), C4 (T–3 h–0.35), C5 (T–3 h–0.35–NH4), C6 (F–3.3 h–1.7 h–0.35–NH4) and C7 (F–2.5 h–0 h–0.35–NH4).

Figure 6

Applied and treated loads for NH4-N and TN, for experimental conditions C1 (T–9 h–0.15), C2 (T–3 h–0.15), C3 (T–3 h–0.25), C4 (T–3 h–0.35), C5 (T–3 h–0.35–NH4), C6 (F–3.3 h–1.7 h–0.35–NH4) and C7 (F–2.5 h–0 h–0.35–NH4).

Close modal
Figure 7 provides insight into the duration of aerobic and anoxic environments for experimental conditions with an aeration strategy based on batch time. Condition C6 (F–3.3 h–1.7 h–0.175–NH4) appeared to have a balanced distribution between anoxic and aerobic environments, with an average duration of 14.7 h/day of anoxic in the PF and 5 h/day in the SF and 7.8 h/day of aerobic in the PF and 18 h/day in the SF. This balance seemed to favor only the nitrification process since 30.3 ± 3.1 mg NO3-N/L were measured at the outlet and the TN removal efficiency was only 24.7 ± 9.2%. Nonetheless, the lack of biodegradable carbon (approximately 50 mg COD/L) was another reason for the lack of denitrification.
Figure 7

Duration in hours per day of aerobic condition (>0.2 mg O2/L) and anoxic condition (<0.2 mg O2/L) for the experimental conditions employing aeration strategy based on batch time. A single oxygen sensor is exhibited for a period of 1 month. PF, primary filter and SF, secondary filter.

Figure 7

Duration in hours per day of aerobic condition (>0.2 mg O2/L) and anoxic condition (<0.2 mg O2/L) for the experimental conditions employing aeration strategy based on batch time. A single oxygen sensor is exhibited for a period of 1 month. PF, primary filter and SF, secondary filter.

Close modal

In contrast, condition C7 (F–2.5 h–0 h–0.175–NH4) exhibited an opposite trend compared to the other conditions: the majority of the day, the saturated zone remained anoxic (17.1 h/day in the PF and 22.9 h/day in the SF), leading to an almost complete denitrification process at the saturated zone, with 30.3 ± 3.1 mg NO3-N/L effluent concentration. However, as the aerobic environment decreased, 5.7 h/day in PF and 0.9 h/day in SF, it also reduced the concentration removal of NH4-N (80.0 ± 5.3%). The dissolved oxygen concentration observed for the C7 condition (supplementary material) was similar to those observed by Uggetti et al. (2016), who also found an average removal of 66% for TN removal.

Despite the extended anoxic periods, denitrification was again hindered by organic carbon limitation. The C/N ratios for C6 and C7 averaged 1.9 and 1.3, respectively. Thus, the carbon source from the batch was not sufficient to meet the denitrification requirement, causing nitrate accumulation (Wang et al. 2020). To overcome this, it would be necessary to increase the carbon source by using an external chemical product that is readily biodegradable, such as methanol or acetic acid (Fu et al. 2022).

Figure 8 illustrates the fluctuation in dissolved oxygen concentration over the course of the day for the C7 condition (F–2.5 h–0 h–0.175–NH4). The aeration was set to start 40 min before the batch, it took around 10 min for the oxygen concentration to reach its maximum value and then when the batch occurred, there were low oxygen concentrations at the saturated layer (on average 0.3 mg/L and maximum 0.6 mg/L). Thus, at the exact moment when the influent dose arrives, the environmental conditions were not completely anoxic. Nonetheless, Figure 7 pointed out that the environment was predominantly anoxic during the day, which created favorable conditions for denitrification. Consequently, the aeration strategy employed in C7 could be optimized if the aeration was timed early (to provide carbon when denitrification is occurring). To further improve understanding of nitrogen dynamics in aerated TWs, online sensors such as pH and redox could be useful for monitoring and controlling the system through online monitoring instruments (Nivala et al. 2020). Such sensors would provide information about the characteristic patterns of biodegradation processes, including the detection of nitrification and denitrification phases. The advantage in this case is to know the timing and duration of nitrification periods or incomplete denitrification stages (Olsson 2012). Ammonium and nitrate sensors are mainly used for nitrogen control; however, the high investment and maintenance costs of these sensors can make their use impractical. To overcome these problems, other physical–chemical parameters can be an option to measure indirectly, like dissolver oxygen (DO), pH, oxidation–reduction potential (ORP), oxygen uptake rate (OUR) and OUR (secondary variables) (Zanetti et al. 2012).
Figure 8

Dissolved oxygen concentration during 24 h for multiple days for the condition C7 (F–2.5 h–0 h–0.35–NH4). A single oxygen sensor is exhibited, located at vertical 1 (84 cm depth). Dashed black lines represent the batches sent to the primary filter.

Figure 8

Dissolved oxygen concentration during 24 h for multiple days for the condition C7 (F–2.5 h–0 h–0.35–NH4). A single oxygen sensor is exhibited, located at vertical 1 (84 cm depth). Dashed black lines represent the batches sent to the primary filter.

Close modal

The Rhizosph'air hybrid aerated TW demonstrated high removal of carbon and suspended solids removal for all seven experimental conditions tested. The conditions varied with influent concentration, duration of aeration, timing of influent loading and aeration, different hydraulic loading rates and addition of ammonium in the influent. TN removal was directly impacted by each operational change. In the experimental conditions with time-based aeration (C1–C5), percent ammonium concentration removal exceeded 90%. Aerating the system for 9 h total per day (6 cycles of 150 min aeration per cycle) resulted in the lowest TN removal (5%). Increasing the HLR (from 0.15 to 0.25 m/day) with 3 h aeration per day, the TN removal improved from 47 to 62%. This is the result of a better balance between oxygen demand and supply. When ammonium was added to the influent (C5: T–3 h–0.35–NH4), nitrification was preserved while denitrification increased, reaching 71% TN concentration removal.

It is important to note the seasonal effect, particularly during winter when the influent on the pilot in this study was diluted by clear water intrusion (C4: T–3 h–0.35 and C6: F–3.3 h–1.7 h–0.35–NH4), resulting in lower influent COD and NH4-N concentrations. Consequently, this had an impact on the efficiency of the system, by decreasing the TN removal, while preserving low outlet TN concentration (around 15 mgN/L). However, despite these effects, complete nitrification still occurs due to efficient oxygen transfer.

When aeration was set according to the batch feeding (C6: F–3.3 h–1.7 h–0.35–NH4 and C7: F–2.5 h–0 h–0.35–NH4), the subsurface conditions favored denitrification with longer periods of anoxic conditions and dissolved oxygen concentrations of approximately 0.2 and 0.9 mg O2/L. The condition C7 (F–2.5h –0 h–0.35–NH4) removed 71% of TN, similar to the performance of condition C5 (T–3 h–0.35–NH4) (71%). The change in aeration strategy from time-based to batch time provided an instantaneous source of carbon for denitrifying bacteria and consequently resulted in lower NO3-N concentrations at the outlet. Nevertheless, this condition reached its maximum denitrification rates but showed an increase in outlet ammonium concentrations (18.3 ± 7.1 mg NH4-N/L effluent concentration). This underlines the importance of both oxygen supply and sequencing.

Therefore, to optimize the hybrid aerated TW system in this study and decrease effluent nitrate and TN concentrations without adding carbon, the following modifications are suggested:

  • Optimize the use of carbon for denitrification by stopping aeration earlier before the influent batch feeding to ensure anoxic environmental conditions when the system is fed.

  • Use online sensors (pH, redox or O2) to identify bending points or patterns related to nitrification and denitrification dynamics. This information could be used to control the operation of the aeration system based on the influent pollutant loads (e.g., treatment on demand).

The authors would like to thank the Office Français de la Biodiversité, Rhône–Méditerranée–Corse and Adour-Garrone water agencies as well as EUR H2O'Lyon for its financial support to the TONIC research project. In addition, we would like to thank Olivier Garcia for his contribution to the field.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Bezbaruah
A. N.
&
Zhang
T. C.
2003
Performance of a constructed wetland with a sulfur/limestone denitrification section for wastewater nitrogen removal
.
Environmental Science & Technology
37
,
1690
1697
.
https://doi.org/10.1021/es020912w
.
Cross
K.
,
Tondera
K.
,
Rizzo
A.
,
Andrews
L.
,
Pucher
B.
,
Istenic
D.
,
Karres
N.
&
McDonald
R.
2021
Nature-Based Solutions for Wastewater Treatment: A Series of Factsheets and Case Studies
.
IWA Publishing. ed
.
https://doi.org/10.2166/9781789062267.
Dotro
G.
,
Langergraber
G.
,
Molle
P.
,
Nivala
J.
,
Puigagut
J.
,
Stein
O.
&
von Sperling
M.
2017
Biological Wastewater Treatment Series – Treatment Wetlands (Volume 7)
.
IWA Publishing
.
https://doi.org/10.2166/9781780408774.
Dou
T.
,
Troesch
S.
,
Petitjean
A.
,
Gabor
P. T.
&
Esser
D.
2016
Wastewater and rainwater management in urban areas: A role for constructed wetlands
. In
Science Direct – International Conference Green Urbanism
.
DWA
2017
Standard DWA-A 262E: Principles for Dimensioning, Construction, and Operation of Wastewater Treatment Plants with Planted and Unplanted Filters for Treatment of Domestic and Municipal Wastewater
.
Fonder
N.
&
Headley
T.
2010
Systematic nomenclature and reporting for treatment wetlands
.
In: Water and Nutrient Management in Natural and Constructed Wetlands (Vymazal, J., ed.). Springer, Dordrecht, pp. 191–220
.
Fu, X., Hou, R., Yang, P., Qian, S., Feng, Z., Chen, Z., Wang, F., Yuan, R., Chen, H. & Zhou, B.
2022
Application of external carbon source in heterotrophic denitrification of domestic sewage: A review. Science of The Total Environment 817, 153061. https://doi.org/10.1016/j.scitotenv.2022.153061.
Kadlec
R. H.
&
Wallace
S. D.
2009
Treatment Wetland
.
CRC Press, New York, NY
.
Li
F.
,
Lu
L.
,
Zheng
X.
&
Zhang
X.
2014
Three-stage horizontal subsurface flow constructed wetlands for organics and nitrogen removal: Effect of aeration
.
Ecological Engineering
68
,
90
96
.
https://doi.org/10.1016/j.ecoleng.2014.03.025
.
Masih
I.
2017
The performance of the intensified constructed wetlands for organic matter and nitrogen removal: A review
.
Journal of Environmental Management
198
,
372
383
.
https://doi.org/10.1016/j.jenvman.2017.04.098
.
Morvannou
A.
,
Forquet
N.
,
Michel
S.
,
Stéphane
T.
&
Molle
P.
2015
Treatment performances of French constructed wetlands: Results from a database collected over the last 30 years
.
Water Science & Technology
71
,
1333
.
https://doi.org/10.2166/wst.2015.089
.
Morvannou
A.
,
Troesch
S.
,
Esser
D.
,
Forquet
N.
,
Petitjean
A.
&
Molle
P.
2017
Using one filter stage of unsaturated/saturated vertical flow filters for nitrogen removal and footprint reduction of constructed wetlands
.
Water Science and Technology
76
,
124
133
.
https://doi.org/10.2166/wst.2017.115
.
Nivala
J.
,
Boog
J.
,
Headley
T.
,
Aubron
T.
,
Wallace
S.
,
Brix
H.
,
Mothes
S.
,
van Afferden
M.
&
Müller
R. A.
2019
Side-by-side comparison of 15 pilot-scale conventional and intensified subsurface flow wetlands for treatment of domestic wastewater
.
Science of The Total Environment
658
,
1500
1513
.
https://doi.org/10.1016/j.scitotenv.2018.12.165
.
Nivala
J.
,
Murphy
C.
&
Freeman
A.
2020
Recent advances in the application, design, and operations & maintenance of aerated treatment wetlands
.
Water
12
,
1188
.
https://doi.org/10.3390/w12041188
.
Olsson
G.
2012
ICA and me – a subjective review
.
Water Research
46
,
1585
1624
.
https://doi.org/10.1016/j.watres.2011.12.054
.
Petitjean
A.
,
Daniau
P.
&
Troesch
S
.
2021
Constructed wetlands for domestic wastewater treatment: From classical ‘French’ systems to optimized intensified treatment wetlands
. In:
International Congress Smallwat21v – Abstract Book
.
Promedio, Badajoz, Spain
, p.
548
.
Prost-Boucle
S.
&
Molle
P.
2021
Intensified treatment wetlands: Forced aeration in Tarcenay, France. Case study
. In
Nature-Based Solutions for Wastewater Treatment. Science for Nature and People Partnership (SNAPP)
(Cross, K., Tondera, K., Rizzo, A., Andrews, L., Pucher, B., Isteni, D., Karres, N. & McDonald, R., eds.).
IWA publishing
, London
p. 344 (180–185)
.
Rochard
J.
2017
Processes of Treatment of Wineries Effluents Adapted to the Organic Wine Sector: Current Situation and Prospects
.
Troesch
S.
,
Dirk
E.
,
Wallace
S.
&
Oirschot
D. V.
2020
Wastewater Purification Device and Utilizations. U.S. Patent 10,526,222 B2 Rhizosph'air®
.
Uggetti
E.
,
Hughes-Riley
T.
,
Morris
R. H.
,
Newton
M. I.
,
Trabi
C. L.
,
Hawes
P.
,
Puigagut
J.
&
García
J.
2016
Intermittent aeration to improve wastewater treatment efficiency in pilot-scale constructed wetland
.
Science of the Total Environment
559
,
212
217
.
https://doi.org/10.1016/j.scitotenv.2016.03.195
.
Val del Río, Á., Campos Gómez, J. L. & Mosquera Corral, A.
2017
Technologies for the Treatment and Recovery of Nutrients from Industrial Wastewater:, Advances in Environmental Engineering and Green Technologies. IGI Global, Hershey, PA. https://doi.org/10.4018/978-1-5225-1037-6.
Wallace
S. D.
2001
System for Removing Pollutants From Water. U.S. Patent 6,200,469 B1
.
Wang
J.
,
Hou
J.
,
Xia
L.
,
Jia
Z.
,
He
X.
,
Li
D.
&
Zhou
Y.
2020
The combined effect of dissolved oxygen and COD/N on nitrogen removal and the corresponding mechanisms in intermittent aeration constructed wetlands
.
Biochemical Engineering Journal
153
,
107400
.
https://doi.org/10.1016/j.bej.2019.107400
.
Wang
H.
,
Liu
Y.
,
Yang
Y.
,
Fang
Y.
,
Luo
S.
,
Cheng
H.
&
Wang
A.
2022
Element sulfur-based autotrophic denitrification constructed wetland as an efficient approach for nitrogen removal from low C/N wastewater
.
Water Research
226
,
119258
.
https://doi.org/10.1016/j.watres.2022.119258
.
Yao
D.
,
Dai
N.
,
Hu
X.
,
Cheng
C.
,
Xie
H.
,
Hu
Z.
,
Liang
S.
&
Zhang
J.
2023
New insights into the effects of wetland plants on nitrogen removal pathways in constructed wetlands with low C/N ratio wastewater: Contribution of partial denitrification-anammox
.
Water Research
243
,
120277
.
https://doi.org/10.1016/j.watres.2023.120277
.
Zanetti
L.
,
Frison
N.
,
Nota
E.
,
Tomizioli
M.
,
Bolzonella
D.
&
Fatone
F.
2012
Progress in real-time control applied to biological nitrogen removal from wastewater. A short-review
.
Desalination
286
,
1
7
.
https://doi.org/10.1016/j.desal.2011.11.056
.
Zhang
L.-y.
,
Zhang
L.
,
Liu
Y.
,
Shen
Y.
,
Liu
H.
&
Xiong
Y.
2010
Effect of limited artificial aeration on constructed wetland treatment of domestic wastewater
.
Desalination
250
,
915
920
.
https://doi.org/10.1016/j.desal.2008.04.062
.
Zhuang
L.-L.
,
Yang
T.
,
Zhang
J.
&
Li
X.
2019
The configuration, purification effect and mechanism of intensified constructed wetland for wastewater treatment from the aspect of nitrogen removal: A review
.
Bioresource Technology
293
,
122086
.
https://doi.org/10.1016/j.biortech.2019.122086
.
Zuur
A. F.
,
Ieno
E. N.
,
Walker
N.
,
Saveliev
A. A.
&
Smith
G. M.
2009
Mixed Effects Models and Extensions in Ecology with R, Statistics for Biology and Health
.
Springer
,
New York, NY
.
https://doi.org/10.1007/978-0-387-87458-6
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-NC-ND 4.0), which permits copying and redistribution for non-commercial purposes with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Supplementary data