Next Article in Journal
Simulation of Water–Energy Nexus of the Spatial Patterns of Crops and Irrigation Technologies in the Cascade Pump Station Irrigation District
Previous Article in Journal
A Hybrid Data-Driven-Agent-Based Modelling Framework for Water Distribution Systems Contamination Response during COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seasonal Enhancement of Nitrogen Removal on Domestic Wastewater Treatment Performance by Partially Saturated and Saturated Hybrid Constructed Wetland

1
Engineering and Environmental Biotechnology Group, Environmental Science Faculty & Center EULA-Chile, Universidad de Concepción, Concepción 4070386, Chile
2
Facultad de Ingeniería y Ciencias, Universidad Adolfo Ibáñez, Avda. Padre Hurtado 750, Viña del Mar 2562340, Chile
*
Author to whom correspondence should be addressed.
Water 2022, 14(7), 1089; https://doi.org/10.3390/w14071089
Submission received: 27 January 2022 / Revised: 9 March 2022 / Accepted: 20 March 2022 / Published: 30 March 2022
(This article belongs to the Topic Sustainable Environmental Technologies)

Abstract

:
The aim of this study is to evaluate seasonal enhancement of nitrogen removal on domestic wastewater treatment performance by partially saturated and saturated HBCWs. To achieve this, two HBCWs consisting of a vertical subsurface flow constructed wetland, followed by a horizontal subsurface flow constructed wetland (VSSF-HSSF) were evaluated. Two saturation levels were used: (a) partially saturated HB1:VSSF1 (0.6 m)-HSSF1 (0.15 m), (b) saturated HB2: VSSF2 (0.8 m)-HSSF2 (0.25 m). Each unit was planted with Schoenoplectus californicus and was operated for 297 days. The removal efficiencies in HB1 and HB2 were above 70%, 86%, 77% and 55% for chemical oxygen demand (COD), total suspended solids (TSS), nitrogen as ammonium (NH4+-N), and total nitrogen (TN), respectively. For VSSF, a higher level of saturation (from 0.6 to 0.8 m) meant a decrease of 17% in the TN removal efficiencies, and for HSSF, an increase from 0.15 to 0.25 m of saturation meant a decrease of 11 and 10% in the NH4+-N and TN removal efficiencies, respectively. Thus, the increase of saturation level in HBCWs reduces the transformation and/or removal of components of the wastewaters to be treated, particularly nitrogen. Through this research, the possibility of optimizing the transformation of nitrogen with partially saturated hybrids can be examined.

1. Introduction

Hybrid constructed wetlands (HBCWs) are based on the sequential utilization of several constructed wetlands (CWs), allowing the combination of each system advantage. Moreover, several configurations have been used, VSSF + VSSF, HSSF + HSSF and VSSF + HSSF, the latter being the most commonly used to treat municipal and industrial wastewater [1,2,3,4,5]. VSSFs are characterized by an oxygen transfer of 28–100 gO2/m2·d, maintaining aerobic conditions and allowing processes such as nitrification with elimination rates of 0.01–4.17 gNH4+/m2·d [6,7,8,9]. In the case of HSSF, oxygen transfer occurs at 0.3–3.2 gO2/m2·d, maintaining anaerobic conditions and allowing processes such as denitrification if nitrate is present with elimination rates of 0.47–0.83 gNO3/m2·d [6,8,9,10]. Several studies have reported that HBCWs present removal efficiencies for COD (58–97.6%), TSS (94–96%), TN (62.8–97%) and NH4+-N (81.1–99.3%) [7,11,12,13,14].
Among the design parameters, it has been stated that the wastewater saturation level impacts the different processes that occur in the HBCWs. [15,16]. Various authors have investigated the influence of partial saturation on the component’s transformation in HSSF or VSSF separately, evidencing the potentialities to further improve the elimination efficiencies of these components. That is how Sanchez-Ramos et al. [17] determined higher removal efficiencies (94% vs. 83% for COD and 74% vs. 32% for NH4+-N) in an HSSF under different saturation levels (0.25 and 0.50 m) than a system that presented a lower saturation level.
Moreover, García et al. [15] determined in a partially saturated HSSF (0.27 m) and in a saturated HSSF (0.50 m) that the best removal efficiencies for COD and NH4+-N were in the partially saturated system with values up to 23.83 and 69.35% higher than those in a saturated HSSF. Huang et al. [18] used a VSSF with different saturation levels (0.05, 0.30, 0.45 and 0.60 m) to improve nitrogen removal. The results indicated that the modification of the saturation level significantly affected the transformation of nitrogen, and the optimal level of saturation was achieved at 0.45 m, where the authors reported that transformation/removal efficiencies of TN and NH4+-N of 76.65 and 82.91% enabled the generation of other processes, such as anammox. Furthermore, Saeed and Sun [19] employed a VSSF under unsaturated and partially saturated conditions with saturation levels of 0.14, 0.29 and 0.43 m, respectively, and conducted a comparative analysis on contaminant removal. They observed that VSSF with a lower saturation level had higher removal efficiencies of BOD5 (47%), TN (52%) and NH4+-N (53%), as it presented an aerobic environment (2.2 mg/L), favoring processes such as nitrification with mean removal rates of 1.30 gNH4+/m2·d.
Specifically, partial saturation in CWs allows the development of an anoxic/anaerobic zone inside the saturated support medium in VSSF and the penetration of atmospheric oxygen toward the saturated zone through overlaid unsaturated bulk media volume without altering commonly practiced overall main mean depth of traditional fully saturated HF wetlands [20]. Despite the advances in research that have been made regarding partial saturation and single-stage CWs, investigations about transformation of nitrogen with partially saturated HBCWs is very limited. Given all that has been mentioned thus far, the main objective of this study is to evaluate the saturation level in the in the treatment of municipal wastewater through hybrid constructed wetlands. Specifically, we will evaluate the saturation of wastewater level in the transformation of nitrogen through hybrid constructed wetlands.

2. Materials and Methods

2.1. Hybrid Constructed Wetlands

The wastewater influent was collected after primary treatment from the wastewater treatment plant in Hualqui, Concepción Province, Biobío Region (Chile) (36°59′26.93″ from south latitude and 72°56′47.23″ from west latitude), which can serve 20,000 inhabitants [21,22]. The influent was transported and stored in 20 L tanks and refrigerated at 4 °C in the dark [23].
Two hybrid constructed systems (HBCWs) were used to scale at a laboratory at Universidad de Concepción. The systems were implemented in a greenhouse-type laboratory with semicontrolled temperatures with between 1–20 °C corresponding to a temperate maritime climate with Mediterranean influence. These HBCWs consisted of a vertical subsurface flow constructed wetland, followed by a horizontal subsurface flow constructed wetland (VSSF-HSSF). In addition, different saturation levels were used. For the VSSF, a partially saturated model (0.6 m) and a saturated model (0.8 m) were employed. Similarly, for the HSSF, partially saturated (0.15 m) and saturated models (0.25 m) were also employed. Therefore, the experimental units for the HBCWs consisted of (a) HB 1: VSSF 1 (0.6 m)-HSSF 1 (0.15 m) and (b) HB 2: VSSF 2 (0.8 m)-HSSF 2 (0.25 m). Each of the four units was planted with the macrophyte species (Schoenoplectus californicus). The VSSF and HSSF systems were inoculated with 4 and 8 seedlings distributed homogeneously, respectively.
Each system was fed discontinuously with 3.6–4.2 L/d for VSSF 1 and VSSF 2 and 1.1–1.7 L/d for HSSF 1 and HSSF 2, respectively. The VSSF had a hydraulic retention time (HTR) of 1 d, and the HSSF had an HRT of 7 d. The effluent was collected using pipeline systems located underneath each wetland. The main design and operational characteristics of the CW are shown in Table 1. Moreover, Figure 1 shows a schematic representation of the HBCWs.

2.2. Sampling Strategy

During the operational period (a total of 297 days), which included summer, autumn, winter, and spring, the in situ parameters, including pH temperature, electrical conductivity, dissolved oxygen (DO) and oxide reduction potential (ORP), were measured weekly using the samplers located at the center of each wetland. Physicochemical parameters, including chemical oxygen demand (COD), total organic carbon (TOC), total suspended solids (TSS), volatile suspended solids (VSS), nitrate-nitrogen (NO2-N), the nitrogen of the nitrate (NO3-N) and ammonium nitrogen (NH4+-N), were measured weekly. Total nitrogen (TN) was measured monthly [22,24]. Removal efficiencies were determined according to the Equation (1) used by Tuttolomondo et al. [25].
R(%) = ((Qi · Ci − Qe · Ce)/(Qi · Ci)) · 100
where Qi is the influent flow rate, Qe is the effluent flow rate, Ci and Ce are the influent and effluent pollutant concentration, respectively.

2.3. Analitycal Methods

To monitor the HBCWs, samples were taken of the influent and the output of each unit, which were filtered through a Whatman 0.45 mm membrane pore size. All protocols described in the standard methods were followed for the physicochemical parameters [26]. Combustion catalytic oxidation at 680 °C and detection by a nondispersive infrared sensor (NDIR) were used to determine TOC. The COD was determined using the colorimetric method (5210-B). The TSS and VSS were determined by the gravimetric method. The nutrients NO3-N were determined by UV, NO2-N and NH4+-N spectrophotometry, which were determined using the colorimetric method, and TN was determined by Spectroquant-Nova 60 (kits of Merck, Kenilworth, NJ, USA). For the in situ parameters, OAKTON portable multiparameter equipment (PC650-480485, Charleston, SC, USA) was used, and for DO, a portable oximeter (HANNA OXI 330i/set HI 9146-04, Villafranca, Italy) was used.
Statistical analysis was performed to compare the removal efficiencies of pollutants of the different CW units. Data were grouped, and the Shapiro–Wilk normality test was performed. A paired t-test was used if the data were normally distributed, and data without normal distribution were analyzed with a Wilcoxon test. The effect of the saturation level of the HBCWs was compared; if they were characterized by a normal distribution, an ANOVA test was performed, but if data were not normally distributed, the Kruskal–Wallis test was used to conduct the analysis. All statistical analyses were performed using the InfoStat statistical program with a significance level of 0.05 [27].

3. Results and Discussion

3.1. Influent Physicochemical Characterization

Table 2 shows the results of the physicochemical characterization of the influent throughout the duration of this study. The influent had an average COD value of 211.67 mg/L and a range of 127.37–401.17 mg/L, reflected in the variations in COD concentrations when considering municipal wastewater [28]. There were applied organic load rates (OLRs) in the range of 18.23–57.76 g/m2·d for HB 1 and 21.27–67.40 g/m2·d for HB 2. Similar results were obtained by García-Ávila et al. [29], who registered an influent with an average value of 222.44 mg/L for municipal wastewater. Likewise, Al-Ajalin et al. [30] reported influents with COD concentrations of 234.00 mg/L regarding municipal wastewater. Additionally, the highest concentrations of COD were recorded during autumn with an average of 296.94 mg/L. In contrast, the lowest concentrations of COD were recorded during spring with a mean of 178.51 mg/L. Thus, the highest variation between seasons was recorded during autumn and spring at 44%.
The TSS presented a range of 32.67–98.96 mg/L with a mean of 53.20 mg/L. Higher concentrations were recorded during summer with a mean of 98.96 mg/L. Moreover, the highest observed variation was recorded between summer and spring at 66%. Similar results were reported by Zurita et al. [31] with a TSS average concentration value of 57.50 mg/L for domestic wastewater.
In terms of nitrogen, the TN recorded a concentration range from 86 to 145 mg/L with higher concentrations during winter (137.67 mg/L average) and up to 23.73% higher than during summer (105.00 mg/L average). NH4+-N had a concentration range from 66.34 to 123.78 mg/L, and a variation of 21.10% was registered in summer (81.79 mg/L average) and autumn (103.66 mg/L average). The applied nitrogen load rates (NLRs) were, on average, 12.80 and 15.01 g NH4+-N/m2·d for VSSF 1 and VSSF 2, respectively. On the other hand, HSSF 1 and HSSF 2 registered applied nitrogen loads of 0.44 and 0.67 g NH4+-N/m2·d, respectively. NO2-N registered an average value of 0.40 mg/L with higher concentrations during summer (0.54 mg/L average) and lower concentrations during spring (0.25 mg/L on average). Regarding NO3-N, a variation of 83.33% during autumn and spring was registered with a concentration range of 0.28 to 2.86 mg/L. Other studies have used influent of the same source and recorded 60.9–79.9 and 76–119 mg/L for NH4+-N and TN, respectively [22,32]. This could be explained since, as it is rural wastewater, the nitrogen contribution is higher because of the different agricultural activities in the area [33].

3.2. Parameters in HBCWs

Figure 2 shows the variation of the in situ parameters for the HBCWs systems. The average temperature was 15.0 °C with a minimum of 13.9 °C in winter and a maximum of 28.2 °C in summer. These data are consistent with Leiva et al. [22] who obtained similar results. They used a laboratory-scale VSSF and reported a temperature range that varied between 8.2 and 22.9 °C. There were no significant differences between the VSSF and HSSF groups (p > 0.05) regarding temperature. In the case of the incidence of the temperature in the biodegradation processes of organic components in anaerobic conditions, studies have concluded that there are no significant effects because bacteria are capable of performing elimination processes even at temperatures as low as 5 °C. Additionally, roots and porous media allow temperatures inside the HBCWs to be 2.3 °C above the outside temperature, thus allowing microbial activity to continue to function properly [34,35,36]. The pH (Figure 2a) presented values that ranged from 5.9 to 8.9, which is considered to be an optimal range for the survival of bacteria in HBCWs (4.0–9.5) [8,37]. Furthermore, higher levels of up to 10% were found in the VSSF (6.5–8.2) than in the HSSF (5.9–7.9). Likewise, the highest values of pH were registered during the autumn-winter period (up to 8.20 and 7.93 for VSSF and HSSF), and a significant difference (p < 0.05) was found when comparing these values to those obtained during the spring-summer period (up to 7.66 and 7.15 for VSSF and HSSF). Similar results were obtained by Marzec et al. [11]. In their study, they reported a pH range of 6.68 and 8.70 for HBCWs with a VSSF-HSSF layout that treated domestic wastewater.
Regarding conductivity, ranges of up to 1814 and 1760 µS/cm were found for VSSF 1 and VSSF 2, respectively. Conversely, HSSF 1 and 2 registered up to 1423 µS/cm and 1402 µS/cm, respectively (Figure 2b). There was a reduction of up to 34% between VSSF and HSSF, which would be explained by the absorption of micro- and macroelements and ions by plants and bacteria and their elimination through adsorption by the action of the roots of plants and sedimental suspended particles [38]. Moreover, a similar decrease was observed by Kyambadde et al. [37], who used an HBCW composed of six staggered VSSFs. In their case, a decrease of 20.72% regarding conductivity concentrations was found, going from 796.3 µS/cm to 631.5 µS/cm.
The DO concentrations in the HBCWs registered a range of 0.1–0.6 mg/L. Average values of 0.22 and 0.23 mg/L were registered for VSSF 1 and VSSF 2 (Figure 2c), respectively. No significant differences (p > 0.05) were found between seasons. In relation to other investigations, the concentrations in VSSF 1 and VSSF 2 were up to 89.52% lower than those reported by Sgroi et al. [39]. They recorded DO concentrations between 0.5 and2.1 mg/L for a VSSF that treated domestic wastewater with an OLR of 40 g COD/m2·d. Regarding HSSF 1 and HSSF 2 (Figure 2c), there were average concentrations of 0.4 mg/L, values up to 78.57% lower than those obtained by Caselles–Osorio et al. [40], who reported DO concentrations of 2.1 mg/L in an HSSF that treated domestic wastewater with an OLR between 4 and 5 g COD/m2·d. Higher concentrations of DO were registered during winter with values up to 0.4 and 0.5 mg/L for VSSF 1 and VSSF 2, respectively, and 0.6 mg/L for HSSF 1 and HSSF 2. Thus, it can be determined that anaerobic conditions prevail in the HBCWs, registering concentrations < 2 mg/L of DO [41].
Regarding the measured ORP, the VSSF (Figure 2d) scored a range from 288 to −25 mV with means of −174.9 and −169.1 mV for VSSF 1 and VSSF 2, respectively. The HSSF registered values ranged from −125 to 77 mV with means of 34.7 and −35.5 mV for HSSF 1 and HSSF 2 (Figure 2d), respectively. A value higher than 100 mV is required to consider an aerobic environment, while a value lower than −100 mV indicates anaerobic conditions [42]. Consequently, the HBCWs were characterized by an anaerobic environment in both the VSSF and HSSF. The oxygen diffusion was not high enough to increase the oxygen concentration within the HBCWs and thus modify the oxidation conditions [43].

3.3. Concentrations and Removal Efficiencies of Organic Matter and Suspended Solids for the HBCWs System

In the VSSF, the average concentration of COD (Table 3) was kept in ranges between 70.25–92.57 mg/L with an average OLR of 23.84 and 27.81 g COD/m2·d for VSSF 1 and VSSF 2, respectively. Significant differences (p < 0.05) were found for VSSF 1 and VSSF 2 between autumn and the rest of the seasons, causing an influent concentration of 296.94 mg/L with concentrations up to 44% higher than those in spring. In the VSSF, the average removal efficiencies for COD (Figure 3a,b) were 56.67 and 57.94% for VSSF 1 and VSSF 2, reaching higher levels of 78.41 and 82.68%, respectively, in autumn. In a recent study, Chang et al. [44] used a laboratory-scale VSSF and reached an average removal efficiency level of 59.9 and 62.8% for an influent with a COD of 288.66 mg/L.
Previous research has reported removal efficiencies that fluctuate between 41 and 92.90% [45,46,47,48,49]. There was only a difference of 2.21% for COD between the partially saturated VSSF system and saturated VSSF with no significant differences (p > 0.05). The results are consistent with the findings of Huang et al. [18]. In their study, no significant differences were found regarding the removal efficiencies of COD when using VSSF at different saturation levels (0.05, 0.30, 0.45 and 0.60 m). Therefore, it can be observed that the removal of organic matter was not affected by the alteration of the saturation level because organic matter decomposes under aerobic/anaerobic conditions as well as by sedimentation and filtration of organic particles [50]. However, different results were reported by Bassani et al. [51]. In their study, they used an unsaturated (0.1 m) and a partially saturated VSSF (0.4 m), and significant differences (p < 0.05) were reported regarding the removal efficiencies of COD (87.9% and 94.3% for the unsaturated and saturated systems, respectively). These results differ from the present study since in the Bassani et al. [51] study, different HRTs (0.6–0.8 d for unsaturated VSSF and 2–2.5 d for partially saturated VSSF) were used. Therefore, a higher HRT allowed a longer period of microorganism interaction with organic matter, thus enabling better removal efficiencies of COD.
HSSF registered average COD concentrations in a range of 25.07–58.92 mg/L (Table 3) with an average OLR of 0.43 and 0.59 g COD/m2·d for HSSF 1 and HSSF 2, respectively. Removal efficiencies (Figure 3c,d) fluctuated between 23.2 and 74.5% with means of 49.53 and 48.56% for HSSF 1 and HSSF 2, registering removal efficiencies of up to 37% higher during winter and spring (up to 74.52 and 71.18%, respectively). These values are lower than those reported by Haddis et al. [52] who used an HSSF that treated municipal wastewater and recorded an average removal efficiency of 65%. Other studies have reported organic matter removal efficiencies in HSSF that range between 33.60 and90.22% [53,54]. A difference of 1.90% for COD was observed between the partially saturated HSSF and the saturated HSSF; therefore, no significant differences (p > 0.05) were found. These results differ from those reported by García et al. [15]. In their study, significant differences (p < 0.05) were found for COD when using a saturated HSSF (0.50 m) and a partially saturated HSSF (0.27 m). Furthermore, in their study, the partially saturated system registered removal efficiencies of up to 83% versus removal efficiencies of 65% for the saturated system. Similarly, Al-Ajalin et al. [30] reported significant differences (p < 0.05) using an HSSF with two saturation levels (0.35 and 0.45 m) with higher removal efficiencies of COD in the system with a lower saturation level (up to 96.94%). On the other hand, Benvenuti et al. [55] indicated that a DO concentration higher than 1.5 mg/L is needed to allow aerobic microbial metabolism. Therefore, it can be expected that in a VSSF and HSSF that included DO < 1.5 mg/L concentrations and a reductive environment (between −288 and 77 mV), both in the partially saturated and saturated systems, organic matter was removed mainly by anaerobic processes [56], thus explaining that there were no significant differences between the saturation levels.
The removal efficiencies of COD in the HBCWs were 79.01 for HB 1 and 79.79% for HB 2. There were no significant differences (p > 0.05) between HB 1 and HB 2 regarding removal efficiencies. In a recent study by Kraiem et al. [57], removal efficiencies of COD of 86% were reported when using an HBCW with a VSSF-HSSF layout that treated rural wastewater. In both VSSF and HSSF, higher removal efficiencies of COD were achieved during cold seasons (autumn-winter), thus concluding that temperature did not affect the removal processes of organic matter. Likewise, Rozema et al. [58] evaluated the effect of temperature in a VSSF, concluding that cold weather did not affect the VSSF treatment efficiencies; furthermore, efficiencies reached up to 99% for COD.
These were removed mainly by physical mechanisms such as sedimentation and filtration in relation to the TSS. Filtration occurred in the roots and stems of macrophyte plants and in the support medium [29]. Additionally, VSSF 1 and VSSF 2 presented removal efficiencies of 83.55 and 82.24%, respectively, and no significant differences were found (p > 0.05) between seasons and amid the partially saturated VSSF and the saturated VSSF. The values reported are higher than those reported by Abdelhakeem et al. [59]. They reported an average removal efficiency of 75% for a laboratory-scale VSSF that used municipal wastewater. HSSF 1 and HSSF 2 recorded removal efficiencies of 57.24 and 50.60%, respectively. Jácome et al. [60] registered removal efficiencies of 77% in an HSSF that treated municipal wastewater with concentrations of TSS in the influent of 41–48 mg/L (up to 22% lower than the current study). There were no significant differences (p > 0.05) between seasons, but there was a significant difference (p < 0.05) of 11.60% between the partially saturated HSSF and saturated HSSF systems, mainly because of the increase in the saturation level, which in turn generated a free flow in the upper part of the support medium. As a consequence, the removal capacity of plants decreased when removing TSS through the roots [61]. Additionally, HB 1 and HB 2 reached average removal efficiencies of 89.81 and 90.53%, respectively. Moreover, there were no significant differences (p > 0.05) between HB 1 and HB 2 regarding the removal efficiencies of TSS. The values reported are consistent with previous research that registered average removal efficiencies of 94–96% [11,12].

3.4. Concentrations and Removal Efficiencies of Nitrogen for the HBCWs System

VSSF 1 (Figure 4a) registered lower concentrations of NH4+-N during spring, at 29.97 mg/L. On the other hand, VSSF 2 (Figure 4b) NH4+-N concentrations fluctuated between 28.53 and96.18 mg/L with lower average concentrations (Table 3) during spring (34.62 mg/L). Additionally, in HSSF 1 (Figure 5c), the minimum concentration of NH4+-N was recorded during summer (6.92 mg/L). In the case of HSSF 2 (Figure 4d), the NH4+-N remained between 7.61 and 42.67 mg/L, recording the lowest concentration during spring (13.87 mg/L).
For the removal efficiencies of NH4+-N, VSSF 1 (Figure 4a) registered a range of 12.28–70.22% with an average of 34.47%. For VSSF 2, a range of 7.52–71.66% (Figure 4b) was recorded with higher efficiencies during spring (up to 71.66%) and an average of 33.40%. Furthermore, significant differences (p < 0.05) were found between VSSF 1 and VSSF in spring and the rest of the seasons. In a similar study, Abdelhakeem et al. [59] reported removal efficiencies of NH4+-N of 36% for a laboratory-scale VSSF that registered an effluent with NH4+-N concentrations of 21 mg/L (37.5% lower than the average concentrations of this study). Moreover, several studies have reported removal efficiencies of NH4+-N that have reached values up to 88.0–97.8% [62,63,64]. Additionally, a difference of only 3.01% for NH4+-N was found between the partially saturated VSSF and the saturated VSSF; thus, no significant differences (p > 0.05) were reported, and nitrification was not affected when modifying the saturation level. In contrast, Huang et al. [18] reported significant differences (p < 0.05) concerning the removal efficiencies of NH4+-N when using different saturation levels (0.05, 0.30, 0.45 and 0.60 m). Furthermore, when changing the saturation level from 0.45 to 0.60 m, a decrease of 12.71% in removal efficiencies was reported. This can be explained by the difficulty of diffusing oxygen at 0.45 m of saturation.
HSSF registered average removal efficiencies regarding NH4+-N of 73.95–65.93% for HSSF 1 and HSSF 2. For HSSF 1 (Figure 4c), removal efficiencies of NH4+-N higher than 60% were recorded, reaching 88.98% in summer. Similarly, HSSF 2 (Figure 4d) registered removal efficiencies in the range of 36.82–87.49%, reaching 87.49% in summer. Furthermore, significant differences (p < 0.05) were found for HSSF 1 and HSSF 2 during the autumn-winter period. The values reported are higher than those obtained by Zurita et al. [31]. In their study, HSSF was used, reaching NH4+-N removal efficiencies of 45.8–48.6%. Additionally, for the partially saturated and saturated HSSFs, although there was a difference of 10.85% in NH4+-N, there were no significant differences (p < 0.05). Sanchez-Ramos et al. [17] obtained similar results when using HSSF at different saturation levels (partially saturated (0.27 m) and saturated (0.50 m)). They reported significant differences (p < 0.05) in NH4+-N since the partially saturated system registered NH4+-N removal efficiencies up to 56% higher than those of the saturated system. Moreover, a higher level of saturation generated changes at the redox state level, creating more reducing conditions and changing oxygen availability as a result of a shorter diffusion of oxygen from the surface and the air-water interface, which are fundamental in the transformation processes of nitrogen components [65].
DO concentrations (<2 mg/L), particularly in the VSSF, would be caused by the oxygen demands for the organic matter oxidation of (min 1.5 mg/L DO), generating a competition of the DO, thus decreasing the ammonium-oxidizing bacteria’s performance and, as a consequence, nitrification [66,67]. According to Liu et al. [67], DO concentrations higher than 1.5 mg/L are needed for nitrification to occur. The previously mentioned factors would help explain the VSSF’s lower performance in removing NH4+-N (34.47% for VSSF vs. 73.95% for HSSF), owing to lower concentrations of available DO (0.2 mg/L for VSSF vs. 0.4 mg/L for HSSF). In addition, different studies have suggested that temperature affects the nitrification process with the optimal temperature ranging from 16 to 32 °C [54,68,69]. With regard to VSSF and HSSF, the removal efficiencies of NH4+-N were higher, reaching 70.22 and 37.00% in the warmer seasons (spring-summer), respectively. In their study, Hua et al. [68] reported a decrease in removal efficiencies of NH4+-N of up to 65% to 20% from summer to winter, caused by the reduction of activity and proliferation of nitrification microorganisms amid a decrease in the temperature [70].
The NO2-N registered values ranged from 0.26–23.53 mg/L for VSSF 1 (Figure 5a) and 0.23–35.63 mg/L for VSSF 2 (Figure 5b). In the case of HSSF 1 (Figure 5c), NO2-N concentrations lower than 2 mg/L were recorded throughout the monitoring period with values ranging between 0.035 and1.61 mg/L. On the other hand, HSSF 2 (Figure 5d) had an average value for NO2-N of 0.019 mg/L with a range of 0.01–3.35 mg/L. The NO3-N in VSSF 1 (Figure 5a) remained at values lower than 13 mg/L, in a range of 0.07–12.56 mg/L. In the case of VSSF 2 (Figure 5b), in general, NO3-N concentrations remained lower than 5 mg/L, except in summer, where concentrations up to 32 mg/L were recorded. The HSSF 1 (Figure 5c) registered concentrations of NO3-N ranged between 10.65 and 63.48 mg/L. On the other hand, HSSF 2 (Figure 5d) recorded NO3-N values that reached 51.32 mg/L in autumn. Furthermore, in the VSSF, there were lower concentrations of NO3-N than in the HSSF (up to 97% in autumn), which could have been removed through denitrification under anaerobic conditions [66]. Alternatively, the decrease in NO3-N could also be explained by simultaneous nitrification/denitrification (SND) processes [71]. This happens as the biofilm is formed in the support medium; thus, an anaerobic microenvironment is created inside, and an aerobic microenvironment is created outside, allowing the coexistence of facultative aerobic and anaerobic microorganisms [72,73].
Conversely, the highest concentrations of NO3-N in HSSF (up to 63.48 mg/L) could be caused by the suppression of denitrification, owing to the lack of an organic carbon source, since it significantly affects the efficiency of denitrification, causing an accumulation of NO3-N in the medium [8,74]. Ding et al. [74] used different C/N relationships (0, 2, 4, 6 and 9) and established that the optimal removal appeared when increasing the C/N relationship to 9, where there were concentrations of NO3-N in the effluent of only 0.03 mg/L for an influent of 19.98 mg/L.
In this study, the lowest concentrations of C/N were found in the HSSF. Furthermore, HSSF 1 registered a C/N relation of 0.8 during the monitoring period, and HSSF 2 registered C/N relations of 1.0, 0.8, 0.7 and 0.7 during summer, autumn, winter, and spring, respectively, thus explaining the accumulation of NO3-N in these systems.
In VSSF 1 (Figure 5a), the TN remained in the range of 70–132 mg/L, and in winter, lower concentrations (70 mg/L) were registered. On the other hand, in VSSF 2 (Figure 5b), TN concentrations fluctuated between 65.00–121.00 mg/L, registering a similar concentration in summer (65.00 mg/L). Additionally, in HSSF 1 (Figure 5c) and HSSF 2 (Figure 5d), the lowest concentrations of TN were recorded in summer, at 33 and 23 mg/L, respectively. Regarding removal efficiencies, VSSF recorded average efficiencies of 22.64 and 18.78% for VSSF 1 and VSSF 2, respectively. These values were up to 62.44% lower than those reported by Zurita et al. [31], who used a VSSF that treated domestic wastewater and reported removal efficiencies of 50% for an influent with concentrations of TN of 28.7 mg/L (up to 77% lower than the concentrations of TN in this study). Additionally, significant differences (p < 0.05) were found between the partially saturated VSSF and the saturated VSS with a difference of 17.05% for TN. This could be caused by changes in the saturation levels creating changes at the DO level (0.62 to 2.2 mg/L) and ORP (−5.39 to 90.93 mV), causing the different microbial communities that act simultaneously in the medium to be affected, thus altering nitrogen removal processes such as nitrification and denitrification [75,76,77]. In the case of the HSSF, average removal efficiencies of TN of 37.03 and 39.40% were recorded for HSSF 1 and HSSF 2, respectively. In addition, no significant differences (p > 0.05) were found between the systems. In a recent study, Jácome et al. [60] used HSSF and established a TN removal efficiency of 37%. On the other hand, Zhu et al. [78] reported removal efficiencies of TN that reached 89.9%; however, this efficiency was due to the TN concentration in the influent of 40 mg/L, which was up to 68% lower than the average concentration of TN in this study (126 mg/L).
In HBCWs, the removal efficiencies of NH4+-N reached 83.50 for HB 1 and 77.28% for HB 2. Saeed et al. [14], who used partially saturated HBCWs, obtained similar results. They established removal efficiencies of NH4+-N of 81.1%. On the other hand, these values are lower than those reported by Ávila et al. [79], who registered values of up to 94%. Regarding TN, average removal efficiencies of 55.12 for HB 1 and 57.40% for HB 2 were reported. Moreover, Nguyen et al. [13] reported lower values than those previously mentioned, where removal efficiencies of 79% were reported using municipal wastewater for HBCWs with a layout of VSSF-HSSF.

3.5. In Situ Spatial Variations and Physical-Chemical Parameters

Table 4 shows the in situ spatial variations and physicochemical parameters. In the case of the COD, it can be observed that during 2019 (data with influent concentrations), the greatest decrease in concentrations occurred in the upper zone with reductions of 58.64 and 53.54 mg/L for VSSF 1 and VSSF 2, namely, removal efficiencies of 34 and 31%, respectively. As the zone being measured decreased, removal efficiencies declined since only 16% was recorded in the lower zone. In a study conducted by Yi et al. [80], it was established that COD concentrations were mainly removed at the entrance with a removal efficiency of 64.23% in the upper zone. Previous research has reported that the anaerobic routes of organic matter removal are slower than aerobic routes [19], which would explain the decrease in the removal efficiency in the lower zones when dealing with systems with a low HRT.
Similarly, for NH4+-N, the highest concentration decline was recorded in the upper zone with 31.43 mg/L for VSSF 1 and 23.53 mg/L for VSSF 2, namely, removal efficiencies of 32% and 24%, respectively. Sanchez-Ramos et al. [17] found that aerobic bacteria are distributed in the first centimeters of the wetland, thus taking advantage of the atmospheric diffusion of oxygen. Therefore, it can be expected in his study that nitrification bacteria were developed in the first centimeters of the wetland, thus explaining the removal rates for NH4+-N despite the ORP and DO values obtained, which hinder the nitrification process.
If we consider the in situ parameters, the pH presented a decrease of 0.1 pH units when the measurement zone descended vertically with a range between 7.1–7.4 for VSSF 1 and VSSF 2, and no significant differences (p > 0.05) were found. Thus, the pH was not vertically affected. Regarding DO, during 2018, there were higher concentrations in the upper zone with 0.5 mg/L for VSSF 1 and 0.4 mg/L for VSSF 2. Moreover, as the measured zone descended, the DO concentrations decreased, reaching 0.3 mg/L for VSSF 1 and 0.2 mg/L for VSSF 2. However, during 2019, the DO concentration remained consistent as the vertical measurement decreased. This vertical decline could be caused by the removal process of organic matter and nitrogen transformation, which in turn would cause a rapid consumption of DO through aerobic respiration and chemical oxidation [42]. In the case of the ORP, a range of values registered from −112.4 mV in the upper zone to −162.4 mV in the lower zone of VSSF 1. On the other hand, VSSF 2 registered a range from −133.6 mV in the upper zone to −193.2 mV in the lower zone. This decline could be caused by the oxygen diffusion mainly in the upper zone of the gravel beds, decreasing as the depth increases, thus creating more reducing conditions [15].

4. Conclusions

In this study, the medium saturation level effect in the VSSF and HSSF does not show significant difference when the COD removal efficiencies are considered. However, when the saturation level of the medium increased, the TSS removal efficiency decreased in 11% in the HSSF. Moreover, this also meant a 17% decrease in the average removal efficiencies of TN in VSSF and a 10% decrease for NH4+-N at HSSF. Due to this, we conclude that when designing HBCWs, it is essential to consider the medium saturation level because it affects the transformation and/or removal of wastewater’s components to be treated. Even more so, the increase of saturation level in HBCWs reduces the transformation of the NH4+-N. This research shows the possibility of optimizing the transformation of nitrogen with partially saturated hybrids constructed wetlands.
As future prospects for this work, it is important to enhance a mass balance of the different forms of nitrogen (NH4+-N, NO3-N NO2-N), considering the microbiological transformations and plant assimilation. In this same line, in future research it would be important to evaluate the gas production and emissions such as CH4, CO2 or N2O due to the use of partially saturated wetlands and also optimize the system to control said emissions.

Author Contributions

Conceptualization, D.L. and G.V.; methodology, D.L. and J.C.; software, G.G.; validation, G.V., D.L. and J.C.; formal analysis, J.C. and D.L.; investigation, J.C. and D.L.; resources, G.V.; data curation, J.C. and G.G.; writing—original draft preparation, J.C. and G.G.; writing—review and editing, G.V., D.L. and G.G.; visualization, G.G.; supervision, G.V.; project administration, G.V. and D.L.; funding acquisition, G.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by ANID/FONDAP/15130015 and Grant No. 3170295 from ANID-FONDECYT (Chile).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ye, F.; Li, Y. Enhancement of nitrogen removal in towery hybrid constructed wetland to treat domestic wastewater for small rural communities. Ecol. Eng. 2009, 35, 1043–1050. [Google Scholar] [CrossRef]
  2. Vera, I.; García, J.; Sáez, K.; Moragas, L.; Vidal, G. Performance evaluation of eight years experience of constructed wetland systems in Catalonia as alternative treatment for small communities. Ecol. Eng. 2011, 37, 364–371. [Google Scholar] [CrossRef]
  3. Ayaz, Ç.; Aktaş, Ö.; Akça, L.; Fındık, N. Effluent quality and reuse potential of domestic wastewater treated in a pilot-scale hybrid constructed wetland system. J. Environ. Manag. 2015, 156, 115–120. [Google Scholar] [CrossRef] [PubMed]
  4. Barco, A.; Borin, M. Treatment performance and macrophytes growth in a restored hybrid constructed wetland for municipal wastewater treatment. Ecol. Eng. 2017, 107, 160–171. [Google Scholar] [CrossRef]
  5. Vymazal, J. The use of hybrid constructed wetlands for wastewater treatment with special attention to nitrogen removal: A review of a recent development. Water Res. 2013, 47, 4795–4811. [Google Scholar] [CrossRef]
  6. Butterworth, E.; Dotro, G.; Jones, M.; Richards, A.; Onunkwo, P.; Narroway, Y.; Jefferson, B. Effect of artificial aeration on tertiary nitrification in a full-scale subsurface horizontal flow constructed wetland. Ecol. Eng. 2013, 54, 236–244. [Google Scholar] [CrossRef]
  7. Gajewska, M.; Skrzypiec, K.; Jóźwiakowski, K.; Bugajski, P. Kinetics of pollutants removal in hybrid treatment wetlands-Case study comparison. Ecol. Eng. 2018, 120, 222–229. [Google Scholar] [CrossRef]
  8. Saeed, T.; Sun, G. A review on nitrogen and organics removal mechanisms in subsurface flow constructed wetlands: Dependency on environmental parameters, operating conditions and supporting media. J. Environ. Manag. 2012, 112, 429–448. [Google Scholar] [CrossRef]
  9. Vymazal, J.; Kröpfelová, L. Multistage hybrid constructed wetland for enhanced removal of nitrogen. Ecol. Eng. 2015, 84, 202–208. [Google Scholar] [CrossRef]
  10. Hickey, A.; Arnscheidt, J.; Joyce, E.; O’Toole, J.; Galvin, G.; O’Callaghan, M.; Walsh, K. An assessment of the performance of municipal constructed wetlands in Ireland. J. Environ. Manag. 2018, 210, 263–272. [Google Scholar] [CrossRef]
  11. Marzec, M.; Gizińska-Górna, M.; Jóźwiakowski, K.; Pytka-Woszczyło, A.; Kowalczyk-Juśko, A.; Gajewska, M. The efficiency and reliability of pollutant removal in a hybrid constructed wetland with giant miscanthus and Jerusalem artichoke in Poland. Ecol. Eng. 2019, 127, 23–35. [Google Scholar] [CrossRef]
  12. Melián, J.; Rodríguez, A.; Araña, J.; Díaz, O.; Henríquez, J. Hybrid constructed wetlands for wastewater treatment and reuse in the Canary Islands. Ecol. Eng. 2010, 36, 891–899. [Google Scholar] [CrossRef]
  13. Nguyen, X.; Chang, S.; Nguyen, T.; Ngo, H.; Kumar, G.; Banu, J.; Nguyen, D. A hybrid constructed wetland for organic-material and nutrient removal from sewage: Process performance and multi-kinetic models. J. Environ. Manag. 2018, 222, 378–384. [Google Scholar] [CrossRef]
  14. Pat-Espadas, A.M.; Rene Loredo Portales, R.; Amabilis-Sosa, L.E.; Gómez, G.; Vidal, G. Review on constructed wetlands for acid mine drainage treatment. Water 2018, 10, 1685. [Google Scholar] [CrossRef] [Green Version]
  15. García, J.; Aguirre, P.; Barragán, J.; Mujeriego, R.; Matamoros, V.; Bayona, J. Effect of key design parameters on the efficiency of horizontal subsurface flow constructed wetlands. Ecol. Eng. 2005, 25, 405–418. [Google Scholar] [CrossRef]
  16. Mander, Ü.; Dotro, G.; Ebie, Y.; Towprayoon, S.; Chiemchaisri, C.; Nogueira, S.; Mitsch, W. Greenhouse gas emission in constructed wetlands for wastewater treatment: A review. Ecol. Eng. 2014, 66, 19–35. [Google Scholar] [CrossRef] [Green Version]
  17. Sanchez-Ramos, D.; Agulló, N.; Samsó, R.; García, J. Effect of key design parameters on bacteria community and effluent pollutant concentrations in constructed wetlands using mathematical models. Sci. Total Environ. 2017, 584, 374–380. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, M.; Wang, Z.; Qi, R. Enhancement of the complete autotrophic nitrogen removal over nitrite process in a modified single-stage subsurface vertical flow constructed wetland: Effect of saturated zone depth. Bioresour. Technol. 2017, 233, 191–199. [Google Scholar] [CrossRef]
  19. Saeed, T.; Sun, G. Pollutant removals employing unsaturated and partially saturated vertical flow wetlands: A comparative study. Chem. Eng. J. 2017, 325, 332–341. [Google Scholar] [CrossRef]
  20. Saeed, T.; Haque, I.; Khan, T. Organic matter and nutrients removal in hybrid constructed wetlands: Influence of saturation. Chem. Eng. J. 2019, 371, 154–165. [Google Scholar] [CrossRef]
  21. Burgos, V.; Araya, F.; Reyes-Contreras, C.; Vera, I.; Vidal, G. Performance of ornamental plants in mesocosm subsurface constructed wetlands under different organic sewage loading. Ecol. Eng. 2017, 99, 246–255. [Google Scholar] [CrossRef]
  22. Leiva, A.; Núñez, R.; Gómez, G.; López, D.; Vidal, G. Performance of ornamental plants in monoculture and polyculture horizontal subsurface flow constructed wetlands for treating wastewater. Ecol. Eng. 2018, 120, 116–125. [Google Scholar] [CrossRef]
  23. López, D.; Fuenzalida, D.; Vera, I.; Rojas, K.; Vidal, G. Relationship between the removal of organic matter and the production of methane in subsurface flow constructed wetlands designed for wastewater treatment. Ecol. Eng. 2015, 83, 296–304. [Google Scholar] [CrossRef]
  24. López, D.; Sepúlveda-Mardones, M.; Ruiz-Tagle, N.; Sossa, K.; Uggetti, E.; Vidal, G. Potential methane production and molecular characterization of bacterial and archaeal communities in a horizontal subsurface flow constructed wetland under cold and warm seasons. Sci. Total Environ. 2019, 648, 1042–1051. [Google Scholar] [CrossRef] [PubMed]
  25. Tuttolomondo, T.; Leto, C.; La Bella, S.; Leone, R.; Virga, G.; Licata, M. Water balance and pollutant removal efficiency when considering evapotranspiration in a pilot-scale horizontal subsurface flow constructed wetland in Western Sicily (Italy). Ecol. Eng. 2016, 87, 295–304. [Google Scholar] [CrossRef]
  26. American Public Health Association (APHA). Standard Methods for the Examination of Wastewater, 20th ed.; APHA: Washington, DC, USA, 1998. [Google Scholar]
  27. Di Rienzo, J.; Casanoves, F.; Balzarini, M.; Gonzalez, I.; Tableda, M.; Robledo, C. Infostat Statistical Software; FCA National University of Córdoba: Córdoba, Argentina, 2011; p. 336. Available online: https://www.researchgate.net/publication/233842986_InfoStat_Version_2011 (accessed on 19 March 2022).
  28. Al-Isawi, R.; Ray, S.; Scholz, M. Comparative study of domestic wastewater treatment by mature vertical-flow constructed wetlands and artificial ponds. Ecol. Eng. 2017, 100, 8–18. [Google Scholar] [CrossRef]
  29. García-Ávila, F.; Patiño-Chávez, J.; Zhinín-Chimbo, F.; Donoso-Moscoso, S.; Del Pino, L.F.; Avilés-Añazco, A. Performance of Phragmites Australis and Cyperus Papyrus in the treatment of municipal wastewater by vertical flow subsurface constructed wetlands. Int. Soil Water Conserv. Res. 2019, 7, 286–296. [Google Scholar] [CrossRef]
  30. Al-Ajalin, F.; Idris, M.; Abdullah, S.; Kurniawan, S.; Imron, M. Effect of wastewater depth to the performance of short-term batching-experiments horizontal flow constructed wetland system in treating domestic wastewater. Environ. Technol. Innov. 2020, 20, 101106. [Google Scholar] [CrossRef]
  31. Zurita, F.; De Anda, J.; Belmont, M. Treatment of domestic wastewater and production of commercial flowers in vertical and horizontal subsurface-flow constructed wetlands. Ecol. Eng. 2009, 35, 861–869. [Google Scholar] [CrossRef]
  32. Sepúlveda, R.; Leiva, A.M.; Vidal, G. Performance of Cyperus papyrus in constructed wetland mesocosms under different levels of salinity. Ecol. Eng. 2020, 151, 105820. [Google Scholar] [CrossRef]
  33. Li, W.; Zhai, L.; Lei, Q.; Wollheim, W.M.; Liu, J.; Liu, H.; Liu, S. Influences of agricultural land use composition and distribution on nitrogen export from a subtropical watershed in China. Sci. Total Environ. 2018, 642, 21–32. [Google Scholar] [CrossRef] [PubMed]
  34. Akratos, C.; Tsihrintzis, V. Effect of temperature, HRT, vegetation and porous media on removal efficiency of pilot-scale horizontal subsurface flow constructed wetlands. Ecol. Eng. 2007, 29, 173–191. [Google Scholar] [CrossRef]
  35. Steer, D.; Fraser, L.; Boddy, J.; Seibert, B. Efficiency of small constructed wetlands for subsurface treatment of single-family domestic effluent. Ecol. Eng. 2002, 18, 429–440. [Google Scholar] [CrossRef] [Green Version]
  36. Steinmann, C.R.; Weinhart, S.; Melzer, A. A combined system of lagoon and constructed wetland for an effective wastewater treatment. Water Res. 2003, 37, 2035–2042. [Google Scholar] [CrossRef]
  37. Kyambadde, J.; Kansiime, F.; Gumaelius, L.; Dalhammar, G. A comparative study of Cyperus papyrus and Miscanthidium violaceum-based constructed wetlands for wastewater treatment in a tropical climate. Water Res. 2004, 38, 475–485. [Google Scholar] [CrossRef] [PubMed]
  38. Kyambadde, J.; Kansiime, F.; Dalhammar, G. Nitrogen and phosphorus removal in substrate-free pilot constructed wetlands with horizontal surface flow in Uganda. Water Air Soil Pollut. 2005, 165, 37–59. [Google Scholar] [CrossRef]
  39. Sgroi, M.; Pelissari, C.; Roccaro, P.; Sezerino, P.; García, J.; Vagliasindi, F.; Ávila, C. Removal of organic carbon, nitrogen, emerging contaminants and fluorescing organic matter in different constructed wetland configurations. Chem. Eng. J. 2018, 332, 619–627. [Google Scholar] [CrossRef] [Green Version]
  40. Caselles-Osorio, A.; Vega, H.; Lancheros, J.; Casierra-Martínez, H.; Mosquera, J. Horizontal subsurface-flow constructed wetland removal efficiency using Cyperus articulatus L. Ecol. Eng. 2017, 99, 479–485. [Google Scholar] [CrossRef]
  41. Vymazal, J.; Kröpfelová, L. Is concentration of dissolved oxygen a good indicator of processes in filtration beds of horizontal flow constructed wetlands? In Wastewater Treatment, Plant Dynamics and Management, 1st ed.; Vymazal, J., Ed.; Springer: Dordrecht, The Netherlands, 2008; pp. 311–317. [Google Scholar]
  42. Ding, Y.; Wang, W.; Song, X.; Wang, Y. Spatial distribution characteristics of environmental parameters and nitrogenous compounds in horizontal subsurface flow constructed wetland treating high nitrogen-content wastewater. Ecol. Eng. 2014, 70, 446–449. [Google Scholar] [CrossRef]
  43. Carballeira, T.; Ruiz, I.; Soto, M. Effect of plants and surface loading rate on the treatment efficiency of shallow subsurface constructed wetlands. Ecol. Eng. 2016, 90, 203–214. [Google Scholar] [CrossRef]
  44. Chang, J.; Wu, S.; Dai, Y.; Liang, W.; Wu, Z. Treatment performance of integrated vertical-flow constructed wetland plots for domestic wastewater. Ecol. Eng. 2012, 44, 152–159. [Google Scholar] [CrossRef]
  45. Abou-Elela, S.; Golinielli, G.; Abou-Taleb, E.; Hellal, M. Municipal wastewater treatment in horizontal and vertical flows constructed wetlands. Ecol. Eng. 2013, 61, 460–468. [Google Scholar] [CrossRef]
  46. Martín, M.; Gargallo, S.; Hernández-Crespo, C.; Oliver, N. Phosphorus and nitrogen removal from tertiary treated urban wastewaters by a vertical flow constructed wetland. Ecol. Eng. 2013, 61, 34–42. [Google Scholar] [CrossRef]
  47. Sylla, A. Domestic wastewater treatment using vertical flow constructed wetlands planted with Arundo donax, and the intermittent sand filters impact. Ecohydrol. Hydrobiol. 2020, 20, 48–58. [Google Scholar] [CrossRef]
  48. Zhao, Y.; Liu, B.; Zhang, W.; Ouyang, Y.; An, S. Performance of pilot-scale vertical-flow constructed wetlands in responding to variation in influent C/N ratios of simulated urban sewage. Bioresour. Technol. 2010, 101, 1693–1700. [Google Scholar] [CrossRef]
  49. Zheng, X.; Zhuang, L.; Zhang, J.; Li, X.; Zhao, Q.; Song, X.; Liao, J. Advanced oxygenation efficiency and purification of wastewater using a constant partially unsaturated scheme in column experiments simulating vertical subsurface flow constructed wetlands. Sci. Total Environ. 2020, 703, 135480. [Google Scholar] [CrossRef]
  50. Saeed, T.; Paul, B.; Afrin, R.; Al-Muyeed, A.; Sun, G. Floating constructed wetland for the treatment of polluted river water: A pilot scale study on seasonal variation and shock load. Chem. Eng. J. 2016, 287, 62–73. [Google Scholar] [CrossRef]
  51. Bassani, L.; Pelissari, C.; da Silva, A.; Sezerino, P. Feeding mode influence on treatment performance of unsaturated and partially saturated vertical flow constructed wetland. Sci. Total Environ. 2020, 754, 142400. [Google Scholar] [CrossRef]
  52. Haddis, A.; Van der Bruggen, B.; Smets, I. Constructed wetlands as nature based solutions in removing organic pollutants from wastewater under irregular flow conditions in a tropical climate. Ecohydrol. Hydrobiol. 2020, 20, 38–47. [Google Scholar] [CrossRef]
  53. Kaseva, M. Performance of a sub-surface flow constructed wetland in polishing pre-treated wastewater-a tropical case study. Water Res. 2004, 38, 681–687. [Google Scholar] [CrossRef]
  54. Ma, Y.; Zhai, Y.; Zheng, X.; He, S.; Zhao, M. Rural domestic wastewater treatment in constructed ditch wetlands: Effects of influent flow ratio distribution. J. Clean. Prod. 2019, 225, 350–358. [Google Scholar] [CrossRef]
  55. Benvenuti, T.; Hamerski, F.; Giacobbo, A.; Bernardes, A.; Zoppas-Ferreira, J.; Rodrigues, M. Constructed floating wetland for the treatment of domestic sewage: A real-scale study. J. Environ. Chem. Eng. 2018, 6, 5706–5711. [Google Scholar] [CrossRef]
  56. Huang, Y.; Ortiz, L.; Aguirre, P.; García, J.; Mujeriego, R.; Bayona, J. Effect of design parameters in horizontal flow constructed wetland on the behaviour of volatile fatty acids and volatile alkylsulfides. Chemosphere 2005, 59, 769–777. [Google Scholar] [CrossRef] [PubMed]
  57. Kraiem, K.; Kallali, H.; Wahab, M.; Fra-vazquez, A.; Mosquera-Corral, A.; Jedidi, N. Comparative study on pilots between annamox favored conditions in a partially saturated vertical flow constructed wetland and a hybrid system for rural wastewater treatment. Sci. Total Environ. 2019, 670, 644–653. [Google Scholar] [CrossRef]
  58. Rozema, E.; Rozema, L.; Zheng, Y. A vertical flow constructed wetland for the treatment of winery process water and domestic sewage in Ontario, Canada: Six years of performance data. Ecol. Eng. 2016, 86, 262–268. [Google Scholar] [CrossRef]
  59. Abdelhakeem, S.; Aboulroos, S.; Kamel, M. Performance of a vertical subsurface flow constructed wetland under different operational conditions. J. Adv. Res. 2016, 7, 803–814. [Google Scholar] [CrossRef] [Green Version]
  60. Jácome, J.; Molina, J.; Suárez, J.; Mosqueira, G.; Torres, D. Performance of constructed wetland applied for domestic wastewater treatment: Case study at Boimorto (Galicia, Spain). Ecol. Eng. 2016, 95, 324–329. [Google Scholar] [CrossRef]
  61. Stecher, M.; Weaver, R. Effects of umbrella palms and wastewater depth on wastewater treatment in a subsurface flow constructed wetland. Environ. Technol. 2003, 24, 471–478. [Google Scholar] [CrossRef]
  62. Dan, T.; Chiem, N.; Brix, H. Treatment of high-strength wastewater in tropical constructed wetlands planted with Sesbania sesban: Horizontal subsurface flow versus vertical downflow. Ecol. Eng. 2011, 37, 711–720. [Google Scholar] [CrossRef] [Green Version]
  63. Dąbrowski, W.; Karolinczak, B.; Gajewska, M.; Wojciechowska, E. Application of subsurface vertical flow constructed wetlands to reject water treatment in dairy wastewater treatment plant. Environ. Technol. 2017, 38, 175–182. [Google Scholar] [CrossRef]
  64. Maina, C.; Mutua, B.; Oduor, S. Evaluating performance of vertical flow constructed wetland under various hydraulic loading rates in effluent polishing. J. Water Sanit. Hyg. Dev. 2011, 1, 144–151. [Google Scholar] [CrossRef]
  65. Nivala, J.; Wallace, S.; Headley, T.; Kassa, K.; Brix, H.; van Afferden, M.; Müller, R. Oxygen transfer and consumption in subsurface flow treatment wetlands. Ecol. Eng. 2013, 61, 544–554. [Google Scholar] [CrossRef]
  66. Fan, J.; Zhang, B.; Zhang, J.; Ngo, H.; Guo, W.; Liu, F.; Wu, H. Intermittent aeration strategy to enhance organics and nitrogen removal in subsurface flow constructed wetlands. Bioresour. Technol. 2013, 141, 117–122. [Google Scholar] [CrossRef]
  67. Liu, L.; Zhao, X.; Zhao, N.; Shen, Z.; Wang, M.; Guo, Y.; Xu, Y. Effect of aeration modes and influent COD/N ratios on the nitrogen removal performance of vertical flow constructed wetland. Ecol. Eng. 2013, 57, 10–16. [Google Scholar] [CrossRef]
  68. Hua, Y.; Peng, L.; Zhang, S.; Heal, K.V.; Zhao, J.; Zhu, D. Effects of plants and temperature on nitrogen removal and microbiology in pilot-scale horizontal subsurface flow constructed wetlands treating domestic wastewater. Ecol. Eng. 2017, 108, 70–77. [Google Scholar] [CrossRef]
  69. Katayon, S.; Fiona, Z.; Megat Mohd Noor, M.; Abdul Halim, G.; Ahmad, J. Treatment of mild domestic wastewater using subsurface constructed wetlands in Malaysia. Int. J. Environ. Stud. 2008, 65, 87–102. [Google Scholar] [CrossRef]
  70. Faulwetter, J.; Gagnon, V.; Sundberg, C.; Chazarenc, F.; Burr, M.; Brisson, J.; Stein, O. Microbial processes influencing performance of treatment wetlands: A review. Ecol. Eng. 2009, 35, 987–1004. [Google Scholar] [CrossRef]
  71. Hu, Y.; Zhao, Y.; Rymszewicz, A. Robust biological nitrogen removal by creating multiple tides in a single bed tidal flow constructed wetland. Sci. Total Environ. 2014, 470, 1197–1204. [Google Scholar] [CrossRef]
  72. Jiang, X.; Tian, Y.; Ji, X.; Lu, C.; Zhang, J. Influences of plant species and radial oxygen loss on nitrous oxide fluxes in constructed wetlands. Ecol. Eng. 2020, 142, 105644. [Google Scholar] [CrossRef]
  73. Zhuang, L.; Yang, T.; Zhang, J.; Li, X. The configuration, purification effect and mechanism of intensified constructed wetland for wastewater treatment from the aspect of nitrogen removal: A review. Bioresour. Technol. 2019, 293, 122086. [Google Scholar] [CrossRef]
  74. Ding, Y.; Song, X.; Wang, Y.; Yan, D. Effects of dissolved oxygen and influent COD/N ratios on nitrogen removal in horizontal subsurface flow constructed wetland. Ecol. Eng. 2012, 46, 107–111. [Google Scholar] [CrossRef]
  75. Adrados, B.; Sánchez, O.; Arias, C.; Becares, E.; Garrido, L.; Mas, J.; Morató, J. Microbial communities from different types of natural wastewater treatment systems: Vertical and horizontal flow constructed wetlands and biofilters. Water Res. 2014, 55, 304–312. [Google Scholar] [CrossRef] [Green Version]
  76. Pelissari, C.; Guivernau, M.; Viñas, M.; García, J.; Velasco-Galilea, M.; Souza, S.; Ávila, C. Effects of partially saturated conditions on the metabolically active microbiome and on nitrogen removal in vertical subsurface flow constructed wetlands. Water Res. 2018, 141, 185–195. [Google Scholar] [CrossRef]
  77. Sun, G.; Zhu, Y.; Saeed, T.; Zhang, G.; Lu, X. Nitrogen removal and microbial community profiles in six wetland columns receiving high ammonia load. Chem. Eng. J. 2012, 203, 326–332. [Google Scholar] [CrossRef]
  78. Zhu, H.; Yan, B.; Xu, Y.; Guan, J.; Liu, S. Removal of nitrogen and COD in horizontal subsurface flow constructed wetlands under different influent C/N ratios. Ecol. Eng. 2014, 63, 58–63. [Google Scholar] [CrossRef]
  79. Ávila, C.; Garfí, M.; García, J. Three-stage hybrid constructed wetland system for wastewater treatment and reuse in warm climate regions. Ecol. Eng. 2013, 61, 43–49. [Google Scholar] [CrossRef]
  80. Yi, X.H.; Jing, D.D.; Wan, J.; Ma, Y.; Wang, Y. Temporal and spatial variations of contaminant removal, enzyme activities, and microbial community structure in a pilot horizontal subsurface flow constructed wetland purifying industrial runoff. Environ. Sci. Pollut. Res. 2016, 23, 8565–8576. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of hybrid constructed wetlands.
Figure 1. Schematic representation of hybrid constructed wetlands.
Water 14 01089 g001
Figure 2. In situ parameters variation: (a) pH; (b) electrical conductivity; (c) dissolved oxygen; and (d) oxide reduction potential. Where VSSF 1 (△); VSSF 2 (▲); HSSF 1(○); and HSSF 2 (●).
Figure 2. In situ parameters variation: (a) pH; (b) electrical conductivity; (c) dissolved oxygen; and (d) oxide reduction potential. Where VSSF 1 (△); VSSF 2 (▲); HSSF 1(○); and HSSF 2 (●).
Water 14 01089 g002
Figure 3. Average of effluent COD concentration (bars) and removal efficiencies (points): (a) VSSF 1: (b) VSSF 2: (c) HSSF 1; and (d) HSSF 2. OLR: organic loading rate.
Figure 3. Average of effluent COD concentration (bars) and removal efficiencies (points): (a) VSSF 1: (b) VSSF 2: (c) HSSF 1; and (d) HSSF 2. OLR: organic loading rate.
Water 14 01089 g003
Figure 4. Effluent average NH4+-N concentration (bars) and removal efficiencies (points): (a) VSSF 1; (b) VSSF 2; (c) HSSF 1; and (d) HSSF 2.NLR: nitrogen loading rate.
Figure 4. Effluent average NH4+-N concentration (bars) and removal efficiencies (points): (a) VSSF 1; (b) VSSF 2; (c) HSSF 1; and (d) HSSF 2.NLR: nitrogen loading rate.
Water 14 01089 g004
Figure 5. Concentrations of different forms of nitrogen (mg/L) NO2-N (○), NO3-N (●), NH4+-N (△) and TN (▼): (a) VSSF 1; (b) VSSF 2; (c) HSSF 1; and (d) HSSF 2.
Figure 5. Concentrations of different forms of nitrogen (mg/L) NO2-N (○), NO3-N (●), NH4+-N (△) and TN (▼): (a) VSSF 1; (b) VSSF 2; (c) HSSF 1; and (d) HSSF 2.
Water 14 01089 g005
Table 1. Main characteristics of a hybrid constructed wetland.
Table 1. Main characteristics of a hybrid constructed wetland.
CharacteristicsUnitVSSFHSSF
Design parameters
Surface aream20.0250.17
Average heightm0.850.30
Water table heightm0.6–0.80.15–0.25
Total volumem30.020.05
Support medium
Type-Sand/GravelGravel
Sizemm1–4/19–2519–25
Porosity-0.2/0.40.4
Operation parameters
HRTd17
HLRL/d3.6–4.21.1–1.7
HRT: Hydraulic retention time; HLR: Hydraulic loading rate.
Table 2. Physicochemical characterization of the influent.
Table 2. Physicochemical characterization of the influent.
Concentration (mg/L) ± SD
SummerAutumnWinterSpring
Average ± SDRangeAverage ± SDRangeAverage ± SDRangeAverage ± SDRange
COD200.67 ± 66.49139.93–314.73296.94 ± 66.75193.93–401.17183.53 ± 42.77132.35–251.49178.51 ± 69.48127.37–257.61
TOC45.72 ± 22.3116.35–85.1387.09 ± 21.4152.45–116.0049.64 ± 11.5735.80–68.0248.28 ± 18.7334.45–69.68
TSS98.96 ± 11.2320.94–106.9242.94 ± 12.6725.20–61.0038.25 ± 13.8227.00–59.3332.67 ± 2.5230.00–35.00
VSS24.25 ± 9.3118.75–35.0040.26 ± 11.9723.20–58.0033.25 ± 14.1815.50–57.3328.20 ± 1.0827.00–29.10
NO2-N0.54 ± 0.040.51–0.570.42 ± 0.180.17–0.720.31 ± 0.050.25–0.490.36 ± 0.120.27–0.49
NO3-N1.69 ± 0.331.26–2.151.98 ± 0.531.34–2.861.22 ± 0.260.58–1.830.33 ± 0.070.28–0.38
NH4+-N81.79 ± 12.4466.34–96.38103.66 ± 7.9192.71–123.7895.81 ± 9.9081.73–106.0397.23 ± 15.7480.05–110.98
TN105.00 ± 19.0086.00–124.00130.67 ± 5.69126.00–137.00137.67 ± 6.66132.00–145.00130.00 ± 2.83128.00–132.00
COD: chemical oxygen demand; TOC: total organic carbon; TSS: total suspended solids; VSS: volatile suspended solids; NO2-N: nitrite nitrogen; NO3-N: nitrate nitrogen; NH4+-N: ammonium nitrogen; TN: total nitrogen.
Table 3. Physicochemical characterization of the effluent for different CW.
Table 3. Physicochemical characterization of the effluent for different CW.
Concentration ± SD
ParameterPeriodVSSF 1VSSF 2HSSF 1HSSF 2
CODSummer81.23 ± 25.3170.25 ± 17.7058.92 ± 18.2357.38 ± 16.16
Autumn81.74 ± 18.3892.57 ± 24.6935.97 ± 8.8441.25 ± 12.12
Winter84.33 ± 17.3976.45 ± 8.3635.25 ± 9.3133.40 ± 6.11
Spring85.63 ± 13.7884.69 ± 21.1035.98 ± 3.7925.07 ± 3.06
TOCSummer30.91 ± 4.5427.29 ± 5.4730.25 ± 13.4413.40 ± 2.31
Autumn18.22 ± 4.1026.96 ± 7.195.20 ± 2.015.56 ± 1.83
Winter18.80 ± 3.8822.26 ± 2.444.57 ± 1.214.72 ± 0.86
Spring19.09 ± 3.0724.67 ± 6.144.67 ± 0.493.54 ± 0.43
TSSSummer7.27 ± 2.075.41 ± 1.414.71 ± 0.423.76 ± 0.71
Autumn7.92 ± 2.688.33 ± 3.233.48 ± 0.533.66 ± 0.61
Winter5.63 ± 1.005.79 ± 1.502.93 ± 0.793.25 ± 0.94
Spring5.60 ± 0.106.10 ± 0.103.50 ± 0.203.27 ± 0.25
VSSSummer6.87 ± 1.558.19 ± 3.513.24 ± 0.853.04 ± 1.12
Autumn7.37 ± 2.636.03 ± 1.473.35 ± 1.912.73 ± 0.53
Winter4.56 ± 1.095.59 ± 1.272.38 ± 0.712.15 ± 0.77
Spring4.53 ± 0.064.73 ± 0.252.67 ± 0.152.73 ± 0.21
Summer63.50 ± 8.2654.55 ± 14.3613.06 ± 5.2516.07 ± 7.62
Autumn83.25 ± 7.9283.91 ± 10.1818.26 ± 6.8424.10 ± 10.40
Winter68.87 ± 15.2778.89 ± 12.3322.30 ± 5.3234.43 ± 6.33
Spring33.35 ± 3.0634.62 ± 7.438.86 ± 1.3513.87 ± 3.09
NH4+-NSummer63.50 ± 8.2654.55 ± 14.3613.06 ± 5.2516.07 ± 7.62
Autumn83.25 ± 7.9283.91 ± 10.1818.26 ± 6.8424.10 ± 10.40
Winter68.87 ± 15.2778.89 ± 12.3322.30 ± 5.3234.43 ± 6.33
Spring33.35 ± 3.0634.62 ± 7.438.86 ± 1.3513.87 ± 3.09
NO2-NSummer0.89 ± 0.281.72 ± 0.530.06 ± 0.020.19 ± 0.09
Autumn0.34 ± 0.090.31 ± 0.080.21 ± 0.070.24 ± 0.06
Winter9.24 ± 1.974.16 ± 1.980.49 ± 0.100.32 ± 0.05
Spring22.79 ± 3.5525.21 ± 9.381.26 ± 0.480.34 ± 0.06
Summer7.44 ± 2.657.76 ± 2.1716.48 ± 1.258.82 ± 0.99
Autumn0.64 ± 0.140.69 ± 0.2525.94 ± 10.1116.53 ± 5.97
Winter2.65 ± 0.912.08 ± 0.6939.37 ± 12.9324.36 ± 11.15
Spring5.18 ± 2.755.49 ± 2.1262.14 ± 1.2735.63 ± 5.35
TNSummer82.00 ± 10.4681.20 ± 17.4430.40 ± 3.1328.60 ± 4.34
Autumn96.67 ± 11.5098.67 ± 11.6857.00 ± 5.2951.33 ± 13.32
Winter114.00 ± 16.37113.67 ± 6.4374.33 ± 7.5173.33 ± 6.66
Spring96.00 ± 5.66103.00 ± 11.3170.00 ± 5.6668.50 ± 9.19
VSSF 1: partially saturated vertical constructed wetland; VSSF 2: saturated vertical constructed wetland; HSSF 1: partially saturated horizontal constructed wetland; HSSF 2: saturated horizontal constructed wetland.
Table 4. Spatial variations of in situ and physicochemical parameters.
Table 4. Spatial variations of in situ and physicochemical parameters.
YearCWLocationConcentration (mg/L)
COD
(mg/L)
NH4+-N (mg/L)Temperature (°C)pHORP
(mV)
DO
(mg/L)
EC
(µS/cm)
2018VSSF1Upper--10.9 ± 1.37.2 ± 0.2−166.9 ± 44.60.5 ± 0.21247.3 ± 393.1
Middle--10.7 ± 1.47.1 ± 0.1−194.1 ± 47.50.3 ± 0.11232.2 ± 360.5
Bottom--10.6 ± 1.27.1 ± 0.1−162.4 ± 127.90.3 ± 0.01034.0 ± 79.2
VSSF2Upper--10.9 ± 0.57.4 ± 0.2−162.9 ± 99.30.4 ± 0.21170.3 ± 167.2
Middle--10.5 ± 0.67.2 ± 0.3−164.3 ± 127.10.3 ± 0.11139.3 ± 111.1
Bottom--10.0 ± 0.57.2 ± 0.2−193.2 ± 83.20.2 ± 0.11031.3 ± 122.7
2019VSSF1Upper111.9 ± 14.565.8 ± 12.611.1 ± 1.17.4 ± 0.2−112.4 ± 24.80.3 ± 0.11339.0 ± 95.3
Middle104.7 ± 19.468.9 ± 13.511.0 ± 1.27.4 ± 0.2−135.5 ± 36.60.3 ± 0.11341.0 ± 98.7
Bottom87.6 ± 18.958.8 ± 13.210.9 ± 1.27.3 ± 0.1−141.2 ± 33.30.3 ± 0.11332.0 ± 111.2
VSSF2Upper117.0 ± 17.073.7 ± 7.111.2 ± 1.17.4 ± 0.2−133.6 ± 14.90.3 ± 0.11361.7 ± 65.6
Middle92.2 ± 20.262.9 ± 17.711.1 ± 1.27.4 ± 0.2−154.7 ± 11.90.3 ± 0.11354.0 ± 54.6
Bottom77.2 ± 20.255.9 ± 15.310.9 ± 1.27.1 ± 0.0−153.9 ± 11.90.3 ± 0.11182.3 ± 38.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Contreras, J.; López, D.; Gómez, G.; Vidal, G. Seasonal Enhancement of Nitrogen Removal on Domestic Wastewater Treatment Performance by Partially Saturated and Saturated Hybrid Constructed Wetland. Water 2022, 14, 1089. https://doi.org/10.3390/w14071089

AMA Style

Contreras J, López D, Gómez G, Vidal G. Seasonal Enhancement of Nitrogen Removal on Domestic Wastewater Treatment Performance by Partially Saturated and Saturated Hybrid Constructed Wetland. Water. 2022; 14(7):1089. https://doi.org/10.3390/w14071089

Chicago/Turabian Style

Contreras, José, Daniela López, Gloria Gómez, and Gladys Vidal. 2022. "Seasonal Enhancement of Nitrogen Removal on Domestic Wastewater Treatment Performance by Partially Saturated and Saturated Hybrid Constructed Wetland" Water 14, no. 7: 1089. https://doi.org/10.3390/w14071089

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop