A review on numerous modeling approaches for effective, economical and ecological treatment wetlands

Review of hydraulic conditions optimization for constructed wetlands

Hydraulic conditions exert a comprehensive and vital influence on constructed wetlands (CWs). However, research on this subject is relatively limited. Hydraulic parameters can be categorized into design and operational parameters based on their properties. The design parameters are represented by the hydraulic gradient, substrate porosity, and aspect ratio, while operational parameters are represented by the hydraulic retention time, hydraulic loading rate, and water depth. These parameters directly or indirectly affect the operational lifespan and pollutant removal performance of CWs. Currently, the primary measures for optimizing the hydraulic conditions of CWs involve hydraulic structure and numerical simulation optimization methods. In this review, we aimed to elucidate the impact of hydraulic conditions on CW performance and summarize current optimization strategies. By highlighting the significance of hydraulic parameters in enhancing pollutant removal and extending operational lifespan, this review provides valuable insights for improving CW design and management. The findings will be useful for researchers and practitioners seeking to optimize CW systems and advance the application of nature-based solutions for wastewater treatment.

Mathematical modeling of subsurface flow constructed wetlands performance for arsenic removal: Review and perspectives

Constructed wetlands are nature-based solutions able to remove different pollutants from water, including arsenic. Arsenic is a pollutant of concern given its toxicity and its presence in water sources worldwide. Despite the increased interest in investigating the performance of constructed wetlands in the treatment of arsenic-contaminated water at the laboratory scale, the application of these solutions at the pilot and full scale is still limited. To understand and predict the removal of arsenic in constructed wetlands, some numerical models have been developed. Among black box models, only first-order models have been proposed, with unsuccessful results. The model that best describes arsenic retention processes in constructed wetlands is RCB-ARSENIC, a mechanistic model adapted from Retraso-CodeBright that simulates arsenic reactive transport. This model includes arsenic precipitation, arsenic sorption on supporting media, arsenic sorption on plants roots and arsenic uptake by plants; represented in the reactive term of the reactive transport equation. Thus, it includes two of the three main processes that remove arsenic in constructed wetlands: precipitation, sorption, and coprecipitation. Despite this, and what is known about arsenic geochemistry, the formulation of these reactive rates requires improvement. In addition, this model was calibrated and validated using data from a single horizontal subsurface flow constructed wetland system, which treated one type of synthetic water. Therefore, it cannot be applied to other types of arsenic-contaminated water or other constructed wetland systems. Moreover, the reactive transport of relevant species -especially iron- and their role in arsenic removal, along with relevant redox reactions associated to the presence of organic matter, oxides or bacteria-, must be included. A comprehensive mechanistic model able to simulate different design, environmental and operation conditions may be used to guide the design of constructed wetlands targeting the removal of arsenic.

Assessment of amoxicillin (AMX) removal from aqueous medium through Rhapis-based bioretention system

Antibiotics can be effectively removed from wastewater using constructed wetlands (C.W.s). However, little is known about using attractive garden plants in C.W.s to eliminate antibiotics. Thus, the current study aims to treat amoxicillin (AMX)-contaminated wastewater through a Rhapis excelsa-based bioretention system (BS). The investigation was done at 15 days hydraulic retention time (HRT) under two conditions: set-1, varied AMX 5 to 25 ppm with constant NPK (nitrogen, phosphorus, potassium) source; and set-2, varied NPK sources with constant AMX (25 ppm). During the study, it was observed that in the set-1 condition with increasing AMX concentration, the removal of AMX through BS decreased; however, in the set-2 experiment, with enhancing NPK source, the performance of the BS treating 25-ppm AMX-contaminated wastewater increased. AMX removal of 2.3%, 66.3%, 60.6%, 52.2%, 46.7%, and 44.9% was achieved for control, BS-1, BS-2, BS-3, BS-4, and BS-5, respectively, during set-1 experiment. However, in the set-2 experiment, 23.4% (control), 43.3% (BS-1), 60.3% (BS-2), 75.9% (BS-3), 88.8% (BS-4), and 99% (BS-5) AMX removal were achieved. Removing pollutants like AMX, COD, PO43--P, NO3--N, and NH4+-N followed first-order kinetics. A positive correlation of COD with AMX was observed through principal component analysis and correlation matrix. The microbial community study was also covered to prioritize the role of microbes in treating AMX through BSs. The AMX treatment through Rhapis excelsa-based BS supported plant growth and development with increasing chlorophyll content, fresh weight, and C, H, N value.

Constructed wetlands as nature-based solutions in managing per-and poly-fluoroalkyl substances (PFAS): Evidence, mechanisms, and modelling

Per- and poly-fluoroalkyl substances (PFAS) have emerged as newly regulated micropollutants, characterised by extreme recalcitrance and environmental toxicity. Constructed wetlands (CWs), as a nature-based solution, have gained widespread application in sustainable water and wastewater treatment and offer multiple environmental and societal benefits. Despite CWs potential, knowledge gaps persist in their PFAS removal capacities, associated mechanisms, and modelling of PFAS fate. This study carried out a systematic literature review, supplemented by unpublished experimental data, demonstrating the promise of CWs for PFAS removal from the influents of varying sources and characteristics. Median removal performances of 64, 46, and 0 % were observed in five free water surface (FWS), four horizontal subsurface flow (HF), and 18 vertical flow (VF) wetlands, respectively. PFAS adsorption by the substrate or plant root/rhizosphere was deemed as a key removal mechanism. Nevertheless, the available dataset resulted unsuitable for a quantitative analysis. Data-driven models, including multiple regression models and machine learning-based Artificial Neural Networks (ANN), were employed to predict PFAS removal. These models showed better predictive performance compared to various mechanistic models, which include two adsorption isotherms. The results affirmed that artificial intelligence is an efficient tool for modelling the removal of emerging contaminants with limited knowledge of chemical properties. In summary, this study consolidated evidence supporting the use of CWs for mitigating new legacy PFAS contaminants. Further research, especially long-term monitoring of full-scale CWs treating real wastewater, is crucial to obtain additional data for model development and validation.

Advances in the process-based models of constructed wetlands and a way forward for integrating emerging organic contaminants

This research examines advancements in the development of process-based models of constructed wetlands (CWs) tailored for simulating conventional water quality parameters (CWQPs). Despite the promising potential of CWs for emerging organic contaminant (EOC) removal, the available CW models do not yet integrate EOC removal processes. This study explores the need and possibility of integrating EOCs into existing CW models. Nevertheless, a few researchers have developed process-based models of other wastewater treatment systems (e.g., activated sludge systems) to simulate certain EOCs. The EOC removal processes observed in other wastewater treatment systems are analogous to those in CWs. Therefore, the corresponding equations governing these processes can be tailored and integrated into existing CW models, similarly to what was done successfully in the past for CWQPs. This study proposed the next generation of CW models, which outlines 12 areas for future work: integrating EOC removal processes; ensuring data availability for model calibration and validation; considering quantitative and sensitive parameters; quantifying microorganisms in CWs; modifying biofilm dynamics models; including pH, aeration, and redox potential; integrating clogging and plant sub-models; modifying hydraulic sub-model; advancing computer technology and programming; and maintaining a balance between simplicity and complexity. These suggestions provide valuable insights for enhancing the design and operational features of current process-based models of CWs, facilitating improved simulation of CWQPs, and integration of EOCs into the modelling framework.

Free water surface constructed wetlands: review of pollutant removal performance and modeling approaches

Free water surface constructed wetlands (FWSCWs) for the treatment of various wastewater types have evolved significantly over the last few decades. With an increasing need and interest in FWSCWs applications worldwide due to their cost-effectiveness and other benefits, this paper reviews recent literature on FWSCWs' ability to remove different types of pollutants such as nutrients (i.e., TN, TP, NH4-N), heavy metals (i.e., Fe, Zn, and Ni), antibiotics (i.e., oxytetracycline, ciprofloxacin, doxycycline, sulfamethazine, and ofloxacin), and pesticides (i.e., Atrazine, S-Metolachlor, imidacloprid, lambda-cyhalothrin, diuron 3,4-dichloroanilin, Simazine, and Atrazine) that may co-exist in wetland inflow, and discusses approaches for simulating hydraulic and pollutant removal processes. A bibliometric analysis of recent literature reveals that China has the highest number of publications, followed by the USA. The collected data show that FWSCWs can remove an average of 61.6%, 67.8%, 54.7%, and 72.85% of inflowing nutrients, heavy metals, antibiotics, and pesticides, respectively. Optimizing each pollutant removal process requires specific design parameters. Removing heavy metal requires the lowest hydraulic retention time (HRT) (average of 4.78 days), removing pesticides requires the lowest water depth (average of 0.34 m), and nutrient removal requires the largest system size. Vegetation, especially Typha spp. and Phragmites spp., play an important role in FWSCWs' system performance, making significant contributions to the removal process. Various modeling approaches (i.e., black-box and process-based) were comprehensively reviewed, revealing the need for including the internal process mechanisms related to the biological processes along with plants spp., that supported by a further research with field study validations. This work presents a state-of-the-art, systematic, and comparative discussion on the efficiency of FWSCWs in removing different pollutants, main design factors, the vegetation, and well-described models for performance prediction.

Quantitative prediction of the hydraulic performance of free water surface constructed wetlands by integrating numerical simulation and machine learning

Quantitative prediction of the design parameter-influenced hydraulic performance is significant for optimizing free water surface constructed wetlands (FWS CWs) to reduce point and non-point source pollution and improve land utilization. However, owing to limitations of the test conditions and data scale, a quantitative prediction model of the hydraulic performance under multiple design parameters has not yet been established. In this study, we integrated field test data, mechanism model, statistical regression, and machine learning (ML) to construct such quantitative prediction models. A FWS CW numerical model was established by integrating 13 groups of trace data from field tests. Subsequently, training, test and extension datasets comprising 125 (5^3), 25 (L₂₅(5⁶)) and 16 (L₁₆(4⁴)) data points, respectively, were generated via numerical simulation of multi-level value combination of three quantitative design parameters, namely, water depth, hydraulic loading rate (HLR), and aspect ratio. The short circuit index (φ₁₀), Morrill dispersion index (MDI), hydraulic efficiency (λ) and moment index (MI) were used as representative hydraulic performance indicators. Training set with large samples were analyzed to determine the variation rules of different hydraulic indicators. Based on the control variable method, φ₁₀, λ, and MI grew exponentially with increasing aspect ratio whereas MDI showed a decreasing trend; with increasing water depth, φ₁₀, λ, and MI showed polynomial decreases whereas MDI increased; with increasing HLR, φ₁₀, λ, and MI slowly increased linearly whereas MDI showed the opposite trend. Finally, we constructed models based on multivariate nonlinear regression (MNLR) and ML (random forest (RF), multilayer perceptron (MLP), and support vector regression. The coefficients of determination (R²) of the MNLR and ML models fitting the training and test sets were all greater than 0.9; however, the generalization abilities of different models in the extension set were different. The most robust MLP, MNLR without interaction term, and RF models were recommended as the preferred models to hydraulic performance prediction. The extreme importance of aspect ratio in hydraulic performance was revealed. Thus, gaps in the current understanding of multivariate quantitative prediction of the hydraulic performance of FWS CWs are addressed while providing an avenue for researching FWS CWs in different regions according to local conditions.

Experimental checking and modeling of the influence of operation conditions on the first order kinetic constants in free water surface wetlands

The most commonly used model in constructed wetlands is the first-order removal model, and first order kinetic constants (k) are the key parameters. The presumption is often made that k are constants. However, it is possible that k are functions of operating conditions, but the influence of operation conditions on k is unclear. In this study, response surface methodology was used to explore the variation patterns of k_a (area rate constants) and k_V (volume rate constants) for the removal of total nitrogen (TN) and total phosphorus (TP) in free water surface (FWS) wetlands. The experimental variables included hydraulic loading rate (HLR), water depth, and inlet concentration (C_in). The results showed that k_V was more variable than k_a, and the area-based first-order model is more suitable for simulating TN and TP in FWS wetlands. Inlet concentration (C_in) was significant for k_a; C_in and water depth were significant for k_V; HLR and the interaction between factors were insignificant. The effects of C_in on k_a and k_V can be described by an upward convex quadratic curve, while the effect of water depth on k_V demonstrates a downward convex quadratic curve. The first-order area rate constant for TN removal was given by k = -47.66 + 22.01 C_in - 1.154 C_in²; the first-order area rate constant for TP removal was given by k = -27.75 + 95.88 C_in - 30.73 C_in². Based on the variation patterns, the traditional k-C model was modified to the k^ψ-C model. The k^ψ-C model produced the best results at simulating the outlet concentration and removal efficiency (RE).

Black-Box Mathematical Model for Net Photosynthesis Estimation and Its Digital IoT Implementation Based on Non-Invasive Techniques: Capsicum annuum L. Study Case

Photosynthesis is a vital process for the planet. Its estimation involves the measurement of different variables and its processing through a mathematical model. This article presents a black-box mathematical model to estimate the net photosynthesis and its digital implementation. The model uses variables such as: leaf temperature, relative leaf humidity, and incident radiation. The model was elaborated with obtained data from Capsicum annuum L. plants and calibrated using genetic algorithms. The model was validated with Capsicum annuum L. and Capsicum chinense Jacq. plants, achieving average errors of 3% in Capsicum annuum L. and 18.4% in Capsicum chinense Jacq. The error in Capsicum chinense Jacq. was due to the different experimental conditions. According to evaluation, all correlation coefficients (Rho) are greater than 0.98, resulting from the comparison with the LI-COR Li-6800 equipment. The digital implementation consists of an FPGA for data acquisition and processing, as well as a Raspberry Pi for IoT and in situ interfaces; thus, generating a useful net photosynthesis device with non-invasive sensors. This proposal presents an innovative, portable, and low-scale way to estimate the photosynthetic process in vivo, in situ, and in vitro, using non-invasive techniques.

Evolution of Ecological Patterns of Poyang Lake Wetland Landscape over the Last One Hundred Years Based on Historical Topographic Maps and Landsat Images

Ecological pattern evolution of Poyang Lake wetland, the largest freshwater lake in China, is critical for regional ecological protection and sustainable development of migratory bird habitats; however, this information is still not fully explored. In this study, we quantitatively reconstructed the spatial distribution and landscape ecological pattern of Poyang Lake wetlands in three periods in the past 100 years based on the military topographic map in the 1930s and the Landsat satellite remote sensing image data in 1979 and 2021. Further, use the Fragstats software to analyze the ecological pattern index of wetland reconstruction results. The results show that the wetland area in the Poyang Lake region has experienced a continuous reduction process over the past 100 years, and it decreased from 3857 km2 in the 1930s to 3673 km2 in the 1970s, and then to 3624 km2 in the 2020s. The current wetland area has decreased by about 6.04% compared with the 1930s. The general trend of changes in the spatial pattern of Poyang Lake wetlands is that the surface water decreases and the open land increases. Nevertheless, the trend has certain spatial differences as a large area of wetlands disappeared in the southwest and west of Poyang Lake and the areas with enlarged wetland density values mainly appeared in the northeastern and northern parts of the study area. The NP (number of patches) in the wetlands of Poyang Lake over the past 100 years showed a downward trend during the 1930s–1970s, and an increasing trend during the 1970s–2010s. Due to the increases of constructed wetlands, the number and density of patches also increased, and PD (patch density) reached a maximum value of 0.142 in 2020s. The LPI (largest patch index) has shown a gradual downward trend in the past 100 years. Compared with the 1930s, the wetlands in 2020s dropped by about 26.64%, and the wetlands further showed a trend of fragmentation. The AI index, which indicates the concentration of wetland patches, reached the maximum value in 2020s, but the LSI (landscape shape index) showed a downward trend in general, indicating that the shape of wetland patches has been simplified over the past 100 years. The research results can provide basic data and decision-making basis for Poyang Lake wetland protection, construction of migratory bird reserve and regional sustainable development.

Machine learning exhibited excellent advantages in the performance simulation and prediction of free water surface constructed wetlands

Optimizing the design and operation parameters of free water surface constructed wetlands (FWS CWs) in runoff regulation and wastewater treatment is necessary to improve the comprehensive performance. In this study, nine machine learning (ML) algorithms were successfully developed to optimize the parameter combinations for FWS CWs. The scale effect of surface area on wetland performance was determined based on consistently smaller predictions (-6.2% to -28.9%) of the nine well-established ML algorithms. The models most suitable for FWS CW performance simulation and prediction were random forest and extra trees algorithms because of their high R² values (0.818 in both) with the training set and low mean absolute relative errors (4.7% and 3.8%, respectively) with the test set. Results from feature analysis of the six tree-based algorithms emphasized the importance of water depth and layout of inlet and outlet, and revealed the negligible effect of the aspect ratio. Feature importance and partial dependence analysis enhanced the interpretability of the tree-based algorithms. The proposed ML algorithms enabled the implementation of an extended scenario at a low cost in real time. Therefore, ML algorithms are suitable for expressing the complex and uncertain effects of the design and operation parameters on the performance of FWS CWs. Acquiring datasets consisting of more extensive, uniform, and unbiased parameter combinations is crucial for developing more robust and practical ML algorithms for the optimal design of FWS CWs.

Utilizing artificial neural network to simulate and predict the hydraulic performance of free water surface constructed wetlands

Optimizing the hydraulic performance of free water surface constructed wetlands (FWS CWs) is of great economic and ecological value. However, there is a complex nonlinear relationship between the hydraulic performance and design parameters of FWS CWs. In this study, an artificial neural network (ANN) was applied to simulate and predict the hydraulic performance corresponding to different combinations of design parameters, and orthogonal design L₉ (3⁴) was used to determine the optimal combination of the important hyperparameters of the ANN. Based on the convenient scenario prediction ability of ANN, sensitivity analysis of different design parameters was carried out by the control variate method and full factor experiment. The results showed that the combination of 3 hidden layers, 15 neural nodes in each hidden layer, 0.001 learning rate, and 8 batch sizes was optimal for the established ANN model, achieving a coefficient of determination of 0.828 in the validation set and a satisfactory prediction effect in the test set. The narrow feature distribution interval in the training set restricted the generalization ability of the ANN model to some extent. Of the four continuous design parameters, the water depth and aspect ratio had an important influence on the effective volume ratio. The layout of inlet and outlet was the most influential design parameter, as confirmed by the full factor experiment of five factors and four levels. The established ANN allowed real-time implementation in an extended scenario at a low cost. This study suggests that the ANN can simultaneously project complex and uncertain effects of several design parameters on wetland performance. In future research, acquiring further comprehensive, impartial, and unbiased experimental datasets is necessary to establish a more robust and generalizing ANN model that can guide the optimal design of FWS CWs.

Prediction of the removal efficiency of emerging organic contaminants in constructed wetlands based on their physicochemical properties

This study investigates the prediction of the removal efficiency of emerging organic contaminants (EOCs) (pharmaceuticals-PhCs, personal care products-PCPs, and steroidal hormones-SHs) in constructed wetlands based on their physicochemical properties (e.g., molecular weight-MW, octanol-water partition coefficient-Log Kow, soil organic carbon sorption coefficient-Log Koc, octanol-water distribution coefficient-Log Dow, and dissociation constant-pKa). The predictive models are formed based on statistical analysis underpinned by principle component, correlation, and regression analyses of a global data set compiled from peer-reviewed publications. The results show that the physicochemical properties of EOCs emerged as good predictors of their removal efficiency. Log Koc, Log Dow, and Log Kow are the most significant predictors, and combination with MW and/or pKa often improved the reliability of the predictions. The best performing model for PhCs was composed of MW, Log Dow, and Log Koc (coefficient of determination-R²: 0.601; probability value-p < 0.05; root mean square error-RMSE: training set: 11%; test set: 27%). Log Kow and Log Koc for PCPs (R²: 0.644; p < 0.1; RMSE: training set: 14%; test set: 14%), and a combination of MW, Log Kow, and pKa for SHs (R²: 0.941; p < 0.1; RMSE: training set: 3%; test set: 15%) formed the plausible models for predicting the removal efficiency. Similarly, reasonably good combined models could be formed in the case of PhCs and SHs or PCPs and SHs, although their individual models were comparatively better. A novel decision support tool, named as REOCW-PCP, was developed to readily estimate the removal efficiency of EOCs, and facilitate the decision-making process.

Improved solute transport and pollutant degradation model of free water surface constructed wetlands considering significant linear correlation between model parameters

To unify the structures of solute transport and pollutant degradation models and evaluate the wetland performance conveniently, a pollutant degradation model combining first-order kinetics with the hybrid solute transport model (plug flow with dispersion + continuous stirred-tank reactor, PFD + CSTR) was developed. Orthogonal tests revealed significant correlation between the model parameters, and the original models were optimized via linear substitution of parameters. The improved PFD + CSTR solute transport model exhibited a satisfactory fit with the original model, and the average relative errors of the determination coefficient (R²) and correlation coefficient were <5%. The multiple linear regressions between the hydraulic indicators and model parameters were reconstructed and exhibited consistent structures between different stages. The degradation constant k_aTN between the original and improved models exhibited high consistency (R² = 0.982). Conclusively, the improved models exhibited good consistency with the original models and allowed rapid and accurate performance evaluation.

Analysis of conservative tracer measurement results inside a planted horizontal subsurface flow constructed wetland filled with coarse gravel using Frechet distribution

We worked out a method in Maple environment to help understand the difficult transport processes in horizontal subsurface flow constructed wetlands filled with coarse gravel (HSFCW-C). With this process, the measured tracer results of the inner points of a HSFCW-C can be fitted more accurately than with the conventionally used distribution functions (Gaussian, Lognormal, Fick (Inverse Gaussian) and Gamma). This research outcome only applies for planted HSFCW-Cs. The outcome of the analysis shows that conventional solutions completely stirred series tank reactor (CSTR) model and convection-dispersion transport (CDT) model do not describe the internal transport processes with sufficient accuracy. This study may help us develop better process descriptions of very complex transport processes in HSFCW-Cs. Our results also revealed that the tracer response curves of planted HSFCW-C conservative inner points can be fitted well with Frechet distribution only if the response curve has one peak.

SURFWET: A biokinetic model for surface flow constructed wetlands

Constructed wetlands are an alternative biotechnology for wastewater treatment that have several advantages over conventional systems. In this work, a biokinetic model for surface flow constructed wetlands is presented (SURFWET). SURFWET belongs to a class of models that are not only interesting from a theoretical viewpoint, as they allow to improve the understanding of the underlying processes; but also from a practical viewpoint, because they can be useful for optimal designs of constructed wetlands, complementing current empirical methods. The proposed model is centered on the intervening physical and biochemical processes involved in pollutant removal in wastewater (organic matter, nitrogen, phosphorus, suspended solids), capturing the interplay of the main agents on contaminant removal (bacteria, macrophytes and phytoplankton). Furthermore, the hydraulic model considers water volume as a variable depending on the outlet hydraulic capacity, and dissolved oxygen has also been introduced as a key driver of reaction kinetics of wetlands. Beyond putting forward a theoretical framework, SURFWET has been applied to simulate a specific case to demonstrate its robustness, in a 12-year-interval simulation. The results show the typical seasonality of this biotechnology, highlighting the importance of dissolved oxygen, which is a key limiting factor on a large number of biochemical processes.

Biotreatment of petroleum refinery wastewater in vertical surface-flow constructed wetland vegetated with Eichhornia crassipes: lab-scale experimental and kinetic modelling

The objectives of this study were to evaluate the potential of vertical surface flow constructed wetland (VSF-CW) vegetated with Eichhornia crassipes to treat petroleum refinery secondary wastewater under tropical conditions. Also, to provide a comparative evaluation of biotreatment kinetic models (traditional first order and other alternative kinetic models) proposed to describe the removal kinetics of organics (biochemical oxygen demand (BOD) and chemical oxygen demand (COD)), nitrate-nitrogen and total petroleum hydrocarbons (TPH) in wetland systems. The refinery secondary wastewater was characterized and treated in five VSF-CWs. Eichhornia crassipes were planted in three VSF-CWs and the remaining two VSF-CWs served as the unvegetated control. The wastewater relatively had high levels of turbidity (18.30 ± 3.88 NTU), BOD (20.40 ± 2.20 mg/L), COD (86 ± 6.0 mg/L), TPH (16.6 ± 1.76 mg/L), oil and grease (18.4 ± 2.00 mg/L), heavy metals (Cadmium (0.034 ± 0.01 mg/L), Lead (0.12 ± 0.05 mg/L), Chromium (0.47 ± 0.01 mg/L), Iron (1.54 ± 0.25 mg/L) and Nickel (0.09 ± 0.01 mg/L)) and Chloride (1412 ± 9.6 mg/L). The vegetated VSF-CWs significantly performed better than the unvegetated control and resulted in the removal efficiencies of 91.5% turbidity, 94.6% BOD₅, 80.2% COD, 92.6% TPH, 90.4% oil and grease, 94% cadmium, 92.5% lead, 93% chromium, 94.8% iron, 92.2% nickel, and 57.7% chloride. The results suggested that the proposed kinetic models were applicable for understanding the removal mechanisms of organics, nitrate-nitrogen and hydrocarbons in VSF-CW systems. Thus, Eichhornia crassipes planted VSF-CW has the potential of treating refinery secondary wastewater to discharge permissible limits.

Numerical Models of Subsurface Flow Constructed Wetlands: Review and Future Development

Numerical model as a simulation tool was used to describe the pollutants transformation and degradation process in constructed wetlands (CWs). It can help provide insight into the “black box” and increase the understanding of the complex processes in CWs. In the last few decades, several process-based numerical models were developed to depict the pollutants removal processes in CWs, which include biochemical model, hydraulic model, reactive-transport model, plants model, clogging model, and coupling model combining two or more sub-models. However, there was a long way to go before fully understanding the decontamination mechanisms of CWs. On the one hand, single or a composite model coupling a small number of sub-models cannot fully reveal the decontamination processes. On the other hand, a comprehensive model including all sub-models of current cognition involves numerous parameters, most of which are interaction and cannot quantitatively determined, thus making the model complex and leading to diffuse interaction. Therefore, in order to describe the reaction processes in CWs more accurately, it is expected that all parameters should be quantified as far as possible in the future model. This study aims to provide a review of the numerical models of CWs and to reveal mechanism of decontamination. Based on the advantages and disadvantages of existing models, the study presented the improvement method and future research direction: (1) new detection/monitoring technique or computing method to quantitatively assess the parameters in CWs models, (2) correcting the simulation errors caused by the assumption of Activated Sludge Models (ASMs) and developing a complete biofilm reaction sub-model, (3) simplification of the comprehensive model, and (4) need of emerging pollutants modeling.

Application of a mathematical model to predict simultaneous reactions in anaerobic plug-flow reactors as a primary treatment for constructed wetlands

Anaerobic digestion technologies offer a set of advantages when they are implemented as a primary treatment phase prior to the use of constructed wetland systems in low cost wastewater facilities. The aim of this study is to describe a model capable of reflecting the complex functioning of anaerobic lagoons, subject to continuous flux in the feed pipe, taking into account that physicochemical properties are subject to a concentration gradient and biochemical ones to simultaneous reactions which depend on each other. Based on both Stokes and advection-diffusion-reaction equations, the proposed model includes twenty-one variables to describe hydraulic, physical, biochemical and physicochemical characteristics that take place in different points of the system and at different moments of time. Drawn up by the International Water Association, the anaerobic digestion model ADM1 is included for the purpose of incorporating the anaerobic processes in the calculation. The finite element method was used to solve the nonlinear, second order partial differential equations of the model. The calculation strategy was designed using a flowchart. Using the open-source FreeFem++ software, a simulation of the mathematical model, in bi-dimensional space, is presented to demonstrate the dynamic behaviour of the proposed model. This yields essential information about the performance of the substrate, cells, and the biochemical reaction products in each of the points within the reactor. Simulations show the potential of this methodology to carry out studies of the behaviour of each of the variables contemplated in the model, as well as comparative studies of the various possible options. In addition, this methodology can be used to help modify the behaviour of the variables based on digester geometry and the boundary values the system is subject to. From the results, it can be concluded that the proposed methodology can be a useful tool for calculating and designing the aforementioned synergistic systems of anaerobic digester plug-flow reactors and constructed wetlands.

Aeration intensity simulation in a saturated vertical up-flow constructed wetland

Simulation and performance results of a saturated vertical up-flow constructed wetland (SVU CW) operated under different operational conditions are presented. The SVU CW consists of two different systems planted with Cyperus alternifolius and Iris pseudacorus, and each system consists of three SVU beds operated in series. The SVU CW operates in continuous aeration (CA) mode using different air-water ratios from 0.5:1 to 4:1. The aerated SVU CW achieves a high (more than 85%) removal of chemical oxygen demand (COD), ammonium (NH₄⁺-N), total nitrogen (TN) and total phosphorus (TP). Furthermore, we simulate the SVU CW using the HYDRUS Wetland Module using the CWM1 biokinetic model under CA mode. According to the simulation results, aeration intensity controls the substrate distribution and growth of bacteria with depth in the SVU CW. Organic matter (OM) and nitrogen are removed in the top region (0-30 cm) of the SVU CW. The root mean square error for COD and NH₄⁺-N is >1.5, whereas R² is >0.99. A good match between observed and simulated data suggests that the CWM1 model is a suitable tool for simulating various processes and bacterial dynamics in aerated SVU CWs.

Effects of flooding regime and meteorological variability on the removal efficiency of treatment wetlands under a Mediterranean climate

The modeling of free-water surface constructed wetlands (FWS-CWs) provides an improved understanding of their processes and constitutes a useful tool for the design and management of these systems. In this work, a dynamic simulation model for FWS-CWs was developed and used to simulate the operation of a FWS-CW proposed for improving the treatment of sewage effluents entering the Tablas de Daimiel National Park in central Spain. The process-based model simulates carbon, nitrogen and phosphorus dynamics, including key hydrological processes for wetlands under a fluctuating Mediterranean semiarid climate. The model allows for the simulation of the operation of FWS-CWs with variable flooding regimes, relating the surface water level to the flooded area and the water outflow. Simulations of the proposed FWS-CW under different water management schemes and scenarios were run, and the consequences of those management strategies on the treatment efficiency were analyzed. Under the Mediterranean climate and geology of the study area, namely, high water losses through evapotranspiration and infiltration, the decrease in nutrient concentrations was higher when the flooded area was reduced in summer than when a constant flooded area was maintained. Moreover, the meteorological variability introduced in different scenarios produced different results in terms of water outflow, but differences in terms of nutrient concentrations were not significant. The ability of the model to simulate different hydrological scenarios and their consequences on water quality makes it a useful decision-support tool.

Grupos funcionales fitoplanctónicos en una laguna algal de alta tasa usada para la biorremediación de lixiviados de rellenos sanitarios

El análisis de comunidades fitoplanctónicas utilizando el enfoque de grupos funcionales (GF), es empleado para predecir la respuesta de la composición de especies fitoplanctónicas a la variación de las condiciones ambientales. El objetivo de este estudio fue evaluar por primera vez bajo condiciones de trópico americano la dinámica de los grupos funcionales de fitoplancton asociados a un sistema de laguna algal de alta tasa fotosintética (LALAT), empleada para la biorremediación de lixiviados de rellenos sanitarios, utilizando la propuesta de grupos funcionales basada en la morfología (MBFG).  Se tomaron muestras de agua para análisis biológico y fisicoquímico durante siete meses. Adicionalmente, se calculó el biovolumen del fitoplancton y se realizó un análisis de correlación canónica para explicar la relación entre las variables físicoquímicas y los grupos funcionales. La LALAT trató un lixiviado de edad intermedia caracterizado por altas concentraciones de nutrientes. Se identificaron cuatro MBFG (IV, V, VI y VII). El grupo V, conformado por fitoflagelados, fue dominante durante todo el periodo de estudio. Este grupo  puede sobrevivir y proliferar en ambientes con altas concentraciones de materia orgánica. Además, contribuye a la transformación de los contaminantes presentes en el lixiviado. El aumento del biovolumen de los grupos V y VII se relacionó con los periodos en los cuales la laguna presentó los menores valores de carbono orgánico disuelto. Finalmente, se encontró que el enfoque de los MBGF puede ser útil para comprender la dinámica de la comunidad de fitoplancton en la LALAT utilizado para el tratamiento de lixiviados en rellenos sanitarios.

Effect of key design parameters on bacteria community and effluent pollutant concentrations in constructed wetlands using mathematical models

Constructed wetlands are currently recognized as an effective environmental biotechnology for wastewater treatment, but the influence of their design parameters on internal functioning and contaminant removal efficiency is still under discussion. In this work, the effect of aspect ratio and water depth on bacteria communities as well as treatment efficiency of horizontal subsurface flow constructed wetlands (HSSF) under the Mediterranean climate was evaluated, using a mathematical model. For this purpose, experimental results from four pilot-scale wetlands of equal surface area but different aspect ratios and water depth were used. The HSSF system was fed with municipal wastewater. The experimental data were simulated using the BIO_PORE model, developed in the COMSOL Multiphysics™ platform. Simulations with the BIO_PORE model fitted well to the experimental results, showing a higher removal efficiency for the shallower HSSF for COD (93.7% removal efficiency) and ammonia nitrogen (73.8%). The aspect ratio had a weak relationship with the bacteria distribution and the removal efficiency. In contrast, the water depth was a factor. The results of the present study confirm a previous hypothesis in which depth has an important impact on the biochemical reactions causing contaminants transformation and degradation.

Modeling organic matter and nitrogen removal from domestic wastewater in a pilot-scale vertical subsurface flow constructed wetland

Constructed wetlands have become an attractive alternative for wastewater treatment. However, there is not a globally accepted mathematical model to predict their performance. In this study, the VS2DTI software was used to predict the effluent biochemical oxygen demand (BOD) and total nitrogen (TN) in a pilot-scale vertical flow constructed wetland (VFCW) treating domestic wastewater. After a 5-week adaptation period, the pilot system was monitored for another 6 weeks. Experiments were conducted at hydraulic retention times (HRTs) in the range of 2-4 days with Typha latifolia as the vegetation. The raw wastewater concentrations ranged between 144-430 and 122-283 mg L(-1) for BOD5 and TN, respectively. A first-order kinetic model coupled with the advection/dispersion and Richards' equations was proposed to predict the removal rates of BOD5 and TN from domestic wastewater. Two main physical processes were modeled in this study, porous material water flow and solute transport through the different layers of the VFCW to simulate the constructed wetland (CW) conditions. The model was calibrated based on the BOD5 and TN degradation constants. The model indicated that most of BOD and TN (88 and 92%, respectively) were removed through biological activity followed by adsorption. It was also observed that the evapotranspiration was seen to have a smaller impact. An additional data series of effluent BOD and TN was used for model validation. The residual analysis of the calibrated model showed a relatively random pattern, indicating a decent fit. Thus, the VS2DTI was found to be a useful tool for CW simulation.

Analysis of conservative tracer measurement results using the Frechet distribution at planted horizontal subsurface flow constructed wetlands filled with coarse gravel and showing the effect of clogging processes

A mathematical process, developed in Maple environment, has been successful in decreasing the error of measurement results and in the precise calculation of the moments of corrected tracer functions. It was proved that with this process, the measured tracer results of horizontal subsurface flow constructed wetlands filled with coarse gravel (HSFCW-C) can be fitted more accurately than with the conventionally used distribution functions (Gaussian, Lognormal, Fick (Inverse Gaussian) and Gamma). This statement is true only for the planted HSFCW-Cs. The analysis of unplanted HSFCW-Cs needs more research. The result of the analysis shows that the conventional solutions (completely stirred series tank reactor (CSTR) model and convection-dispersion transport (CDT) model) cannot describe these types of transport processes with sufficient accuracy. These outcomes can help in developing better process descriptions of very difficult transport processes in HSFCW-Cs. Furthermore, a new mathematical process can be developed for the calculation of real hydraulic residence time (HRT) and dispersion coefficient values. The presented method can be generalized to other kinds of hydraulic environments.

Application of divided convective-dispersive transport model to simulate conservative transport processes in planted horizontal sub-surface flow constructed wetlands

We have created a divided convective-dispersive transport (D-CDT) model that can be used to provide an accurate simulation of conservative transport processes in planted horizontal sub-surface flow constructed wetlands filled with coarse gravel (HSFCW-C). This model makes a fitted response curve from the sum of two independent CDT curves, which show the contributions of the main and side streams. The analytical solutions of both CDT curves are inverse Gaussian distribution functions. We used Fréchet distribution to provide a fast optimization mathematical procedure. As a result of our detailed analysis, we concluded that the most important role in the fast upward part of the tracer response curve is played by the main stream, with high porous velocity and dispersion. This gives the first inverse Gaussian distribution function. The side stream shows slower transport processes in the micro-porous system, and this shows the impact of back-mixing and dead zones, too. The significance of this new model is that it can simulate transport processes in this kind of systems more accurately than the conventionally used convective-dispersive transport (CDT) model. The calculated velocity and dispersion coefficients with the D-CDT model gave differences of 24-54% (of velocity) and 22-308% (of dispersion coeff.) from the conventional CDT model, and were closer to actual hydraulic behaviour.

Adaptive neuro-fuzzy inference system for real-time monitoring of integrated-constructed wetlands

Monitoring large-scale treatment wetlands is costly and time-consuming, but required by regulators. Some analytical results are available only after 5 days or even longer. Thus, adaptive neuro-fuzzy inference system (ANFIS) models were developed to predict the effluent concentrations of 5-day biochemical oxygen demand (BOD5) and NH4-N from a full-scale integrated constructed wetland (ICW) treating domestic wastewater. The ANFIS models were developed and validated with a 4-year data set from the ICW system. Cost-effective, quicker and easier to measure variables were selected as the possible predictors based on their goodness of correlation with the outputs. A self-organizing neural network was applied to extract the most relevant input variables from all the possible input variables. Fuzzy subtractive clustering was used to identify the architecture of the ANFIS models and to optimize fuzzy rules, overall, improving the network performance. According to the findings, ANFIS could predict the effluent quality variation quite strongly. Effluent BOD5 and NH4-N concentrations were predicted relatively accurately by other effluent water quality parameters, which can be measured within a few hours. The simulated effluent BOD5 and NH4-N concentrations well fitted the measured concentrations, which was also supported by relatively low mean squared error. Thus, ANFIS can be useful for real-time monitoring and control of ICW systems.

Field application of a planted fixed bed reactor (PFR) for support media and rhizosphere investigation using undisturbed samples from full-scale constructed wetlands

This study presents a novel method for investigations on undisturbed samples from full-scale horizontal subsurface-flow constructed wetlands (HSSFCW). The planted fixed bed reactor (PFR), developed at the Helmholtz Center for Environmental Research (UFZ), is a universal test unit for planted soil filters that reproduces the operational conditions of a constructed wetland (CW) system in laboratory scale. The present research proposes modifications on the PFR original configuration in order to allow its operation in field conditions. A mobile device to obtain undisturbed samples from real-scale HSSFCW was also developed. The experimental setting is presented with two possible operational configurations. The first allows the removal and replacement of undisturbed samples in the CW bed for laboratory investigations, guaranteeing sample integrity with a mobile device. The second allows the continuous operation of the PFR and undisturbed samples as a fraction of the support media, reproducing the same environmental conditions outside the real-scale system. Investigations on the hydrodynamics of the adapted PFR were carried out with saline tracer tests, validating the proposed adaptation. Six adapted PFR units were installed next to full-scale HSSFCW beds and fed with interstitial liquid pumped from two regions of planted and unplanted support media. Fourteen points were monitored along the system, covering carbon fractions, nitrogen and sulfate. The results indicate the method as a promising tool for investigations on CW support media, rhizosphere and open space for studies on CW modeling, respirometry, kinetic parameters, microbial communities, redox potential and plant influence on HSSFCW.

How to increase microbial degradation in constructed wetlands: influencing factors and improvement measures

Microorganisms play a vital role in degradation of multiple pollutants in constructed wetlands (CWs). Thus, the search for methods to improve microbial degradation in CWs is crucial. This study provides a review of critical parameters including availability of organic carbon, redox condition, temperature, pH, presence of plants, media characteristics and their influences on microbial processes. Current strategies focusing on regulation of carbon source, redox condition, and choice of substrates to enhance microbial activity in CWs are also described. A special emphasis is given to the application of bioaugmentation to enhance microbial activities in wetland in future research.

The Cartridge Theory: a description of the functioning of horizontal subsurface flow constructed wetlands for wastewater treatment, based on modelling results

Despite the fact that horizontal subsurface flow constructed wetlands have been in operation for several decades now, there is still no clear understanding of some of their most basic internal functioning patterns. To fill this knowledge gap, on this paper we present what we call "The Cartridge Theory". This theory was derived from simulation results obtained with the BIO_PORE model and explains the functioning of urban wastewater treatment wetlands based on the interaction between bacterial communities and the accumulated solids leading to clogging. In this paper we start by discussing some changes applied to the biokinetic model implemented in BIO_PORE (CWM1) so that the growth of bacterial communities is consistent with a well-known population dynamics models. This discussion, combined with simulation results for a pilot wetland system, led to the introduction of "The Cartridge Theory", which states that the granular media of horizontal subsurface flow wetlands can be assimilated to a generic cartridge which is progressively consumed (clogged) with inert solids from inlet to outlet. Simulations also revealed that bacterial communities are poorly distributed within the system and that their location is not static but changes over time, moving towards the outlet as a consequence of the progressive clogging of the granular media. According to these findings, the life-span of constructed wetlands corresponds to the time when bacterial communities are pushed as much towards the outlet that their biomass is not anymore sufficient to remove the desirable proportion of the influent pollutants.

Simulation of batch-operated experimental wetland mesocosms in AQUASIM biofilm reactor compartment

In this study, a mathematical biofilm reactor model based on the structure of the Constructed Wetland Model No.1 (CWM1) coupled to AQUASIM's biofilm reactor compartment has been used to reproduce the sequence of transformation and degradation of organic matter, nitrogen and sulphur observed in a set of constructed wetland mesocosms and to elucidate the development over time of microbial species as well as the biofilm thickness of a multispecies bacterial biofilm in a subsurface constructed wetland. Experimental data from 16 wetland mesocosms operated under greenhouse conditions, planted with three different plant species (Typha latifolia, Carex rostrata, Schoenoplectus acutus) and an unplanted control were used in the calibration of this mechanistic model. Within the mesocosms, a thin (predominantly anaerobic) biofilm was simulated with an initial thickness of 49 μm (average) and in which no concentration gradients developed. The biofilm density and area, and the distribution of the microbial species within the biofilm were evaluated to be the most sensitive biofilm properties; while the substrate diffusion limitations were not significantly sensitive to influence the bulk volume concentrations. The simulated biofilm density ranging between 105,000 and 153,000 gCOD/m(3) in the mesocosms was observed to vary with temperature, the presence as well as the species of macrophyte. The biofilm modeling was found to be a better tool than the suspended bacterial modeling approach to show the influence of the rhizosphere configuration on the performance of the constructed wetlands.

Structural and functional loss in restored wetland ecosystems

In restored wetland ecosystems with apparently natural hydrology and biological structure, biogeochemical function may remain degraded, even a century after restoration efforts.

Wetlands are among the most productive and economically valuable ecosystems in the world. However, because of human activities, over half of the wetland ecosystems existing in North America, Europe, Australia, and China in the early 20th century have been lost. Ecological restoration to recover critical ecosystem services has been widely attempted, but the degree of actual recovery of ecosystem functioning and structure from these efforts remains uncertain. Our results from a meta-analysis of 621 wetland sites from throughout the world show that even a century after restoration efforts, biological structure (driven mostly by plant assemblages), and biogeochemical functioning (driven primarily by the storage of carbon in wetland soils), remained on average 26% and 23% lower, respectively, than in reference sites. Either recovery has been very slow, or postdisturbance systems have moved towards alternative states that differ from reference conditions. We also found significant effects of environmental settings on the rate and degree of recovery. Large wetland areas (>100 ha) and wetlands restored in warm (temperate and tropical) climates recovered more rapidly than smaller wetlands and wetlands restored in cold climates. Also, wetlands experiencing more (riverine and tidal) hydrologic exchange recovered more rapidly than depressional wetlands. Restoration performance is limited: current restoration practice fails to recover original levels of wetland ecosystem functions, even after many decades. If restoration as currently practiced is used to justify further degradation, global loss of wetland ecosystem function and structure will spread.

Wetlands, which include tropical mangroves and boreal peatlands, are among the most valuable ecosystems in the world because they provide critical ecosystem goods and services, such as carbon storage, biodiversity conservation, fish production, water purification, and erosion control. As global change accelerates the loss of wetlands, attempts are increasing to restore this fragile habitat and its associated functioning. There has been no global evaluation, however, of how effective such restoration efforts have been. Here, we present a meta-analysis of the biological structure (driven mostly by plant communities) and biogeochemical functioning (driven primarily by the storage of carbon in wetland soils) of 621 wetland sites. Our analysis suggests that even a century after restoration efforts, these parameters remained on average 26% and 23% (respectively) lower in restored or created wetlands than in reference wetlands. Our results also indicate that ecosystem size and the environmental setting significantly affect the rate of recovery. Recovery may be more likely and more rapid if more than 100 contiguous hectares of habitat are restored. In warm climates, and in settings linked to riverine or tidal flows, recovery can also proceed more rapidly. In general, however, once disturbed, wetlands either recover very slowly or move towards alternative states that differ from reference conditions. Thus, current restoration practice and wetland mitigation policies will maintain and likely accelerate the global loss of wetland ecosystem functions.