Next Article in Journal
Hydraulic Responses and Flow Regulation in Multi-Demand Water Transfer Systems
Next Article in Special Issue
Toxicity Abatement of Wastewaters from Tourism Units by Constructed Wetlands
Previous Article in Journal
Analysis of Wet Soil Granular Flow down Inclined Chutes Using Discrete Element Method
Previous Article in Special Issue
Diverse Arbuscular Mycorrhizal Fungi (AMF) Communities Colonize Plants Inhabiting a Constructed Wetland for Wastewater Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The State of the Art of Clogging in Vertical Flow Wetlands

Institute of Sanitary Engineering and Water Pollution Control, University of Natural Resources and Life Sciences Vienna (BOKU), Muthgasse 18, A-1190 Vienna, Austria
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2019 / Revised: 1 November 2019 / Accepted: 8 November 2019 / Published: 15 November 2019
(This article belongs to the Special Issue Recent Progress of Constructed Wetland for Wastewater Treatment)

Abstract

:
Clogging in vertical flow (VF) wetlands is an important process influencing water purification processes. The main contributing factors are the growth of microorganisms within the filter media, the accumulation of suspended solids on top of the wetland, as well as within the filter media. Both processes lead to a decrease of the available pore space, hence changing the soil’s hydraulic properties. This will alter the water flow and cause malfunctioning of the system. This paper summarizes the state of the art of the prevailing physical, biological and chemical processes influencing clogging in VF wetlands. Different design and operational parameters are discussed to give a better understanding on their influence to prevent malfunctioning. Based on a literature review, a detailed overview on experimental as well as modelling studies carried out is presented. The main conclusions are that on the one hand, important insights on clogging processes in VF wetlands have been gained but, on the other hand, design parameters such as intermittent loading operation and the grain size of the filter media are not well represented in those studies. Clogging models use different conceptual approaches ranging from black box models to process based models.

1. Introduction

Vertical flow (VF) wetlands are nowadays referred to as a mainstream technology for wastewater treatment. Due to their proven robustness as well as low costs for operation and maintenance, this technology is implement around the world to treat various polluted waters, e.g., domestic wastewater, combined sewer overflow, stormwater and all sorts of industrial wastewaters. [1,2,3]. VF wetlands evolved from the work of Seidel [4] and have been developed due to the legal requirements of ammonia reduction in the effluent. When nitrification is required, only a system providing a positive oxygen balance is appropriate. This led to the design of VF wetlands with intermittent feeding, where oxygen can enter the filter during the resting time [5,6,7]. Using the standard design, VF wetlands are suited for treating mechanically pre-treated wastewater. For the treatment of a raw wastewater, the so called French VF wetlands were developed. This system is capable of nitrification, but also stabilizes the accumulated sludge. For detailed descriptions, the reader is referred to the listed literature [1,8,9,10].
The main process leading to the malfunction of VF wetlands is clogging. Despite specific recommendations for design and operation by national guidelines and scientific literature [9,10], clogging is still listed as the number one problem. The main reasons for clogging of well-designed systems are failure or mismanagement of the pre-treatment and overloading of organic matter (OM) as well as suspended solids (SS). This has to be considered especially during the start-up time when plants and microorganisms are not yet fully developed. Two main processes leading to clogging can be identified, namely clogging due to SS and bioclogging. SS, which can be retained within the filter matrix due to sedimentation, filtration and attachment. Bioclogging occurs due to the growth of bacteria. What has to be kept in mind is that bacteria growth within the media is required to efficiently degrade wastewater pollutants, hence there needs to be a distinction between normal operation and dysfunctional clogging [11].
In order to get a perspective on the experimental and modelling research on clogging in VF wetlands, a literature search was conducted in the Scopus and Web of Science databases using the following search string: “vertical AND wetland AND clogging”. The results (22 February 2019) provided 123 hits for Scopus and 108 for Web of sciences (18 results of the latter were not included in Scopus). Based on the abstracts, publications not having their main emphasis on clogging were eliminated. In the next step, conference papers and publications not available in English language were sorted. Based on the remaining 54 publications, the following sections are compiled, including the cited primary literature. Eleven publications were chosen to present the development of experimental studies on clogging within the last two decades. In the second part, eight publications were used to describe the different modelling approaches on clogging.

2. Clogging Processes in VF Wetlands

Clogging is defined as blockage and/or filling of the pore space which alters and hinders the movement of water and pollutants. Such behaviour is described by a decrease of the main soil’s hydraulic properties of the filter media, namely the infiltration rate (ki) at the surface, the hydraulic conductivity (kf) within the filter body, as well as the local porosity. Due to the heterogeneity of the clogging development within the media, the effects can range from minor changes of the local flow velocities to the development of preferential flow paths and the exclusion of main parts of the filter. The built up of a clogging mat on the surface of VF wetlands leads to an uneven distribution of the wastewater, followed by exclusion of filter parts, surface ponding and overflow. Consequences thereof can be the risk of human contact with wastewater and the development of bad odour [12].
Clogging can be categorized by its spatial occurrence into inner and outer blockage, and is mainly influenced by the pre-treatment, the applied organic loading rate (OLR, g COD m−2 day−1) and SS concentration (mg L−1). Other influential parameters include the grain size distribution of the filter media, the feeding intervals and resting periods. Outer blockage occurs due to the accumulation of SS, namely particulate organic and inorganic matter, at the surface, while inner blockage happens within the filter media [13,14,15]. Two main types of clogging are identified, namely bioclogging and clogging due to SS. Bioclogging describes the blockage of the pore space due to surplus sludge production based on biofilm (BF) growth and microbial uptake. A BF is composed of extracellular polymeric substances (EPS) and cells which undergo replication producing biopolymers when sufficient nutrients are present [16]. Clogging due to SS describes the accumulation of SS on and its deposition within the filter. Caselles-Osorio et al. [17] differentiate two types of accumulated solids, namely interstitial solids entrapped in the empty space between the media and adhered solids which are tightly adsorbed by the media. Interstitial solids are easily released by washing and filtering while adhered solids need ultrasonic treatment to be released by the media.
The prevailing mechanisms for SS removal from the wastewater are transport, filtering and attachment [18]. Based on the particle size, different physio-chemical processes occur. Four main filtration mechanisms can be identified, namely (i) surface filtration, (ii) straining, (iii) bridging and (iv) physical-chemical filtration [19]. Surface filtration describes the accumulation of particles on the wetland surface which are larger than the pore size of the media. Straining defines a blockage within the media when the average particle size is smaller than the pore size but, depending on the local pore-size distribution, those particles are trapped within smaller pores. Bridging occurs when several particles smaller than the pore size arrive at the same time and together block the pore. This is mainly driven by high flow velocities [20].
Particles smaller than 10 µm are subjected to physical-chemical processes, namely electrostatic and London–Van der Waals forces. Those forces determine if particle–media or particle–particle associations are favourable. This depends on the charge of the media and the particle. If particles are already attached, the capacity of the media to attract more is limited, hence only small changes of the porosity are observed. Changes in water chemistry may destabilize particles which can form aggregates and attach via London–Van der Waals forces. These larger aggregates can have an impact on the local porosity. The actual attachment is controlled by hydrodynamic forces [19].
The aggregates of the filter media, as well as the BF, provide surface area for the attachment of decomposing biological matter. Microbial BF growth decreases the local porosity and the hydraulic conductivity, because cells and their EPS plug the pores between aggregates. Biological clogging will be enhanced if nutrient loadings are high [21].

3. Design and Operational Suggestions to Avoid Clogging

3.1. Design

In recent years, several national design guidelines for VF wetlands have been updated. The main suggestion is that by fulfilling these recommendations, the risk of clogging should be eliminated [1,9]. An overview on the main design parameters is given in Table 1. Detailed information including expected treatment efficiencies can be found elsewhere [1].
Generally, three types of VF wetlands are implemented, namely the single stage system, the two-stage system and the French system. The single stage and the two-stage system are designed for the treatment of mechanically pre-treated wastewater. The single stage is designed with a surface area of 4 m2 per person equivalent (PE) and is operated in unsaturated free flow conditions.
The two-stage VF wetland system [10,22] is designed by using 2 m2 per PE. The first stage has a coarser filter media with a grain size of 1–4 mm and an impounded drainage layer (15 cm). In this faster flowing media, OM and nitrogen (after nitrification) is available in the impounded layer leading to denitrification. In the second stage (0–4 mm), the remaining COD and NH4-N is treated.
The French VF wetland system [8] treats raw wastewater. Due to accumulation of sludge on the surface, this system is fed intermittently for 3.5 days followed by a resting period of 7 days, where the sludge is mineralized and the surface is opened up by wind induced plant stem movement. The recommended surface area is 2–2.5 m2 per PE. One of the biggest VF wetlands implemented is a French system in Moldova treating the wastewater of up to 20,000 PE [23].

3.1.1. Filter Media

The grain size distribution of the media is crucial for the general operational strategy and pollutant removal performance of VF wetlands and can be related to the recommended OLR and hydraulic loading rate (HLR). Based on several studies, three different media sizes are identified, namely single stage systems using washed sand with a grain size of 0–2 mm and 0–4 mm respectively, two stage systems using a 1–4 mm sand in the first stage and 0–4 mm sand in the second stage and the French system, fed with raw wastewater having a 2–8 mm gravel in the first stage and a 0–4 mm sand in the second stage [9] (Table 1).
To prevent clogging, it is crucial to use washed sand. Parameters listed in the design guidelines to ensure the quality of the sand are the effective grain size d10 and the uniformity coefficient U = d60/d10. d10 and d60 describe the grain size under which 10% and 60% of the grains pass during sieve analysis (by weight), respectively. The German design guideline DWA-A262 [24] requests a grain size of 0–2 mm d10 = 0.2 to 0.4 mm and U < 5, for a size of 0–4 mm d10 = 0.25 to 0.4 mm and U < 5.

3.1.2. Loading

The OLR is one of the most crucial parameters. For single stage systems, the OLR can range from 20 to 27 g COD m−2 day−1 [9,13,25,26], 80 g COD m−2 day−1 for a two stage system [9,27] and 100–350 g COD m−2 day−1 for French VF wetlands [8,9]. It is important to note that the OLR is also a temperature dependent parameter, as the prevailing biological processes decrease with decreasing temperatures [28]. In order to avoid clogging, it is recommended to not only choose the appropriate OLR but also the related dosing strategy, which will be discussed in the next section.

3.1.3. Dosing Strategy

Intermittent feeding has proven to be a sufficient strategy to foster the preferred aerobic conditions in VF wetlands [25] and is compulsory based on several national design guidelines [1,9]. This operational strategy is also beneficial to hinder clogging as it supports an equilibrium in BF production [14].
During the resting period, three main processes prevail [29]. As no nutrients are available outside of the cells, microorganisms enter an endogenous respiration state which leads to a decrease in BF and EPS respectively. Oxygen can enter the filter by convection and diffusion. Convection is driven by the gradient of air pressure resulting in a vacuum produced by the effluent water. Diffusion can be related to the actual resting time. Platzer [6] gives a diffusion rate of 1 g O2 m−1 h−1. The availability of O2 supports degradation processes during the resting period as accumulated SS are hydrolysed with its products adding another source of COD which has to be further degraded by bacteria [1,30].

3.2. Pre-Treatment

Pre- or primary treatment is an important step as it prepares the influent for biological treatment. Generally, it consists of a screen, grit removal and primary sedimentation. Grit includes well settleable inorganic and organic solids like sand, gravel, coffee grounds and seeds. The aim of the primary treatment is the reduction of the solid load on the wetland and should not exceed a SS concentration of 100 mg L−1 [13,31,32]. Regular maintenance of the pre-treatment facility is of the upmost importance to avoid solids carryover which will lead to clogging [1,33].

3.3. Effects of Macrophytes

In VF wetlands, macrophytes have mainly physical functions. Developed plant roots loosen the filter media, and plant stems break up the surface area due to wind-induced movement [34]. This effect is especially important for opening up the sludge deposition layer in the French system [35]. Decomposition of plant litter on the surface can lead to an increase in the COD and SS load [36], but a relation to clogging effects is not mentioned in the literature. Maintenance recommendations [1] when cutting is carried out in fall state that the plant material shall be left on the filter to provide insulation during winter.

4. Remediation Strategies

Available methods to remediate a fully clogged system can be characterized as destructive active or passive treatment methods. As shown in Table 2, destructive methods, either excavation and replacement or washing and reuse, are associated with high costs but have proven to be highly effective [25,37,38,39].
While passive treatment is carried out by applying a resting period for the clogged filter, active treatment offers various options. Several authors investigated the use of hydrogen peroxide (H2O2), an aggressive oxidation agent [40,41,42]. Based on their findings, good results are achieved within a short time period, but the method of application is crucial. The injection of 1600 L of H2O2 directly into the media of a 670 m2 horizontal (HF) wetland was effective immediately. The use on the surface of a VF wetland (75 m2, 100 L of 35% H2O2) was reported as insufficient, as the peroxide reacted with the clog matter on top of the filter but not within the subsurface media [42].
Based on a laboratory study, three solubilisation agents, namely sodium hydroxide (NaOH), hydrochloric acid (HCl) and sodium hypochlorite (NaClO), were under investigation for the remediation of planted columns filled with gravel [43]. The clogging of the columns was accelerated by feeding with activated sludge. After a clogging state was reached (effective porosity was decreased from 27% to 5%), each agent was applied daily for one week. Results showed an increase of the effective porosity up to 15% (HCl), 18% (NaOH) and 25% (NaClO) respectively. Damage of plant roots was observed when using NaClO.
In a recent study [44], the use of hydrolytic enzymes, namely α-glucoamylase and β-glucan, to catalyse and hydrolyse large molecular polymers of the EPS was investigated to treat bioclogging. Results of column experiments show a positive effect by monitoring the hydraulic conductivity. For the remediation of a full scale treatment wetland (TW), large quantities are needed, resulting in high costs.
The application of earthworms accelerates the biodegradation of accumulated SS and can force translocation of clog matter from the pore space to the surface and thereby reduces the dry weight of clog matter on average by 56% [45]. Studies on VF wetlands [46,47] show an efficient remediation performance when using 0.5 kg m−2 of earthworms after 10 days, as well as their positive effect during ongoing operation.
Passive treatment is equal to applying an overall resting period on the clogged bed. During this time, bacteria enter an endogenous respiration state as no nutrient supply from outside the cell is occurring. These starvation conditions limit the microbial growth, hence decay of BF is the prevailing process [48,49]. Langergraber et al. [14] observed a recovery time of 14 days after clogging of an intermittent feeding system operating at an HLR of 250 mm day−1 and a TSS loading rate of 48 g TSS m−2 day−1. Based on a VF column experiment using artificial wastewater without SS, the recovery time of the filter was achieved after one to three weeks [50,51], while others suggest 10 days in summer and 20 days in winter conditions [52].

5. Experimental Studies on Clogging

5.1. Methods to Evaluate Clogging Behavior

The quantification of accumulated SS within the pore matrix is mainly carried out by extracting a media sample, washing of the clog matter and determine its dry weight after drying at 105 °C. To determine the OM content, loss-on-ignition (heating to 550 °C) is used resulting in the measured mass of volatile SS (VSS) [53]. Based on this method, the conclusion was drawn that preferential accumulation of solids takes place in the upper layer of the VF wetland [54]. Two different types of SS accumulation are distinguished, namely interstitial solids which are trapped in the pore space and adhered solids which are tightly adsorbed on the media. While interstitial solids are the solids washed off, adhered solids are determined by ultrasonic treatment [17].
A direct correlation of accumulated solids and the decrease of hydraulic conductivity as well as the porosity is difficult to assess. The form, density and water retention characteristics of the solids effects the free volume available for water movement and therefore the local water content, hence the hydraulic conductivity [25,55]. To determine the hydraulic conductivity two in situ methods can be used, namely the falling head and the constant head method [56].
The change in porosity of different layers can be determined by the drainable porosity method. This method is mainly applicable for column experiments where the height of the water table can be more easily determined. Thereby the volume of water needed to fill the system from the bottom to a certain height is measured. A similar approach can be used to determine the hydraulic conductivity using Darcy’s law by measuring the drainage volume over time [15].
In situ measurements of water content and matrix potential are done by using time domain reflectometry (TDR) probes and tensiometers. This provides detailed data describing the water flow for each feeding at different depths. Change in water content is a good indicator for the occurrence and extant of clogging [14].
The actual observation of the pore space by producing images is a promising method to directly investigate clogging processes. Kim and Forquet [57] used the “Thin Section Method” to produce 2D images of undisturbed samples. They analysed the images and estimated phase volume fractions and void size distribution. By using X-ray Computed Tomography in combination with image processing algorithms, actual 3D images of the pore space can be produced [58]. With this novel method, changes in pore morphology and retained solids’ agglomeration by drying can be illustrated.

5.2. Experimental Studies

Based on the literature research, two main types of experiments are identified, namely experiments using continuous feeding and others using intermittent feeding operations. While the intermittent feeding operation is carried out in unsaturated conditions, continuous feeding experiments are carried out in unsaturated as well as saturated conditions.
Different types of wastewaters are used based on the objective of the study, namely real or synthetic wastewater. Table 3 presents the mixtures of synthetic wastewaters used in clogging investigations. Information on the different VF wetland experiments including filter media, loading rates and feeding schemas is given in Table 4 and Table 5.
Langergraber et al. [14] investigated the effect of different HLRs on media clogging of pilot scale VF wetlands designed based on the Austrian guideline [65] and fed with real wastewater. Two different grain sizes of sand were used, namely 0–4 mm and 1–4 mm. Clogging occurred within the fine media for a HLR of 150 L m−2 day−1 after two months and in the coarse media at a HLR of 250 L m−2 day−1 after one month. Measured NH4-N effluent concentrations were low also for high HLRs with removal rates of 99% for the fine and 90% for the coarse media. For early detection of clogging monitoring the water content within the substrate showed good results. For the 0–4 mm substrate, a linear increase in the residual water content before a dose was observed two weeks before clogging occurred. Based on soil physical and structural investigations, the main cause for clogging was found to be the SS load and not BM growth.
Hua et al. carried out several studies [60,66] using continuously fed column experiments to investigate clogging behaviour. In the first study [60], the accumulation of particulate inorganic matter, namely river sand, on top and within the filter was determined. The columns were loaded with a HLR of 500 L m−2 day−1 and an SS loading rate of 300 g SS m−2 day−1 and filled with different grain sizes of gravel, namely 3, 10 and 20 mm. Based on monitoring the infiltration rate and the effective porosity at different depths, the study concludes that firstly a blanket-like deposition layer is formed on top of the filter, followed by particles accumulating within the pore space in the first 0–4 cm. In a follow-up study, Hua et al. [66] investigate the time until clogging occurs, which was observed as 75, 120 and 240 days.
In the study by Sani et al. [15], column experiments using pea gravel of 10 and 20 mm grain size were fed with real wastewater. The feeding pattern used is described as continuously with a resting period, namely 72 or 36 h of feeding and 48h of rest. Two different OLRs were used, namely 11 and 22 g COD m−2 day−1. The effluent quality, as well as the kf and the SS distribution within the filter, were monitored. A change in porosity/kf was measured but no clogging leading to malfunction of the columns occurred during the measurement campaign of one year.
Zhao et al. [59] investigated the effect of both bioclogging and clogging due to SS in column experiments using two synthetic wastewater receipts, one glucose based and one starch based. Each process was investigated by using five columns fed with different OLR. The effective porosity as well as the infiltration rate decreased with an increase in the OLR and was higher overall in the starch-fed system. The reduction of the effective porosity over the depth of the filter was as follows: for the glucose-fed system, an effect was measured over the total filter depth with a decreasing extant while the starch-fed systems were only affected within the first 25 cm. The influence of bioclogging on the infiltration rate was negligible. Therefore, the authors conclude that bioclogging is mainly accelerating already occurring clogging behaviour.
Hua et al. [50] investigated the effects of a resting period after clogging due to BF growth occurred. Therefore, glucose based synthetic wastewater was used to only trigger BF growth. The saturated columns were fed continuously for 90 days at a high HLR of 500 L m−2 day−1 and a BOD concentration of 600 mg L−1. The influence a resting time was determined after 3, 7 and 10 days in different filter depths. Main BF accumulation was found in the top 20 cm of the filter. The recovery ratio for kf was high after 7 and 10 days but started decreasing, while for the porosity it stayed constant. One explanation, therefore, was the remaining inert BF within the pore space. To evaluate the reduction of BF, polysaccharides and proteins were measured as a good correlation is stated. The anthrone–sulfuric acid method is used to measure polysaccharides and Coomassie Brilliant Blue G-250 colorimetry for the measurement of proteins. The quantity of EPS is calculated as the sum of both, and a significant decrease over the resting period is observed.
In a follow up study, Hua et al. [29] investigated the contribution of evapotranspiration (ET) as well as oxygen diffusion into the media to alleviate clogging. Oxygen transport was measured using an optical probe (Unisense, Aarhus, Denmark). Results show that oxygen enters the resting filter over time, supporting OM degradation. ET was calculated using an equation by Fujimaki and Inoue [67], which depends on the relative humidity (RH) within the media. RH was determined by measuring the water content. The results point out a contribution of ET but did not separate the evaporation and the transpiration by the plants.
Another experimental study investigated the effects of BF growth on the local porosity in a saturated VF column with a full water recirculation batch operation [62]. A constant oxygen influent concentration of 4–5 mg L−1 was supplied within the recirculation process. Every week, the column was drained and filled with new synthetic wastewater. Thereby the porosity was evaluated. BF development caused preferential flow and local porosity reduction in dead zones within the column after one year.
Song et al. [63] investigated the effect of different packing strategies on the nitrogen treatment performance and clogging behaviour in VF column experiments. Four layers with different grain sizes of sand where either in increasing or decreasing size stacked. A third column was packed universally with 5–6 mm sand. The operation followed a fill and drain pattern with an HRT of 10 h and a drainage time of 2 h. The OLR was alternated between high and low loads. BF accumulation was determined by the decrease in media water storage volume. COD and NH4-N removal rates of 39%–46% and 28%–39% respectively for the high load regime. The coarse to fine layering (from top to bottom) showed the least reduction in the water storage volume.
The effect of aeration in a saturated VF wetland on clogging is that aeration influences the distribution over the depth of accumulated SS and promotes its mineralization as well as BF growth [68]. Compared to a non-aerated system, bioclogging becomes more of a factor.
Petitjean et al. [11] conducted a column experiment to investigate the interplay between clogging and oxygen concentration to identify parameters capable of determining clogging when it occurs. Their results show that the absence of oxygen within the filter was leading to clogging one day after. This also correlated with a decrease of measured NO3-N while measured NH4-N increased only nine days after clogging occurred. Measuring the oxygen concentration within the media using a portable gas analyser (Dräger X am 7000) is easy to do, but the as the oxygen content varies quickly, the time of observation is critical and therefore not practical. Observation of the NO3-N effluent concentration combined with an observation of surface ponding can be a valid method to identify if clogging is occurring.
An investigation in the spatial and temporal distribution of accumulated insoluble inorganic particles was carried out by Yang et al. [31]. In total, 80%–90% of accumulated solids were distributed in within the first 6 cm [60] of the columns with a media height of 10 cm. In total, 15 columns were established. Every 15 days, one column was dismantled and separated into five 2 cm layers. Extracted media was washed and filtered using a microporous membrane of 0.45 µm. The dry weight represents the amount of interstitial solids. The media is than subjected to ultrasound treatment (Jiupin JP10-200) and filtered (0.45 µm) to quantify adhered solids. The particle size of accumulated solids is determined using a laser particle size analyser (Malvern MS-2000). The prevailing mechanism on particle accumulation was adsorption. Therefore, the pore space was reduced but not filled completely. At the end of the experiment, the number of large particles in the first 0–2 cm increased, leading to a significant reduction of the pore space. Based on these results, the importance of pre-treatment can be pointed out to reduce the SS influent load.
The interaction of EPS and SS is presented in a study by Zhou et al. [64]. Based on four saturated VF columns the impact of different synthetic wastewater mixtures on the clogging behaviour was under investigation. The four mixtures were glucose-based (Bioclogging), starch-based (Ip-clogging), starch with bacteriostatic (Op-clogging) and glucose with starch (C-clogging). Feeding was carried out until clogging, which occurred after 48, 24, 24 and 28 days respectively. In the following resting period, a reduced water level to 20 cm was established and porosity and kf were monitored. The main conclusions were as follows: bioclogging is a main driver while SS accelerates clogging. Undegraded starch particles are compacted by EPS and reduce the porosity. Bioclogging, Op-clogging and C-clogging can be relieved by a resting period, especially when the water level is reduced.

6. Model Concepts

In the field of treatment wetlands (TW), two main types of models can be identified, namely (i) models to design TW which is mainly done using so called black box models and (ii) the application of process-based models to describe the processes occurring within the wetland. The latter type of models are generally built up using several sub-models, usually including a water flow model, a transport model, a biokinetic model, a plant model and a clogging model [69]. Within the next section, several approaches on modelling clogging are described (including other sub-models when provided by the authors). Different approaches can be identified ranging from black box models to white box models. Black box models are data dependent and based on functions relating input to output without describing the occurring processes. White box models on the other hand describe the actual processes. Data in white box models are used for calibration and validation of the models only and are not needed for developing the model itself. In between, grey box models combine both approaches based on the existing knowledge [70].
The general approach for modelling clogging is described by determine the pore space reduction based on BF development and SS retention. This change alters the initial soil hydraulic parameters (SHP) of the initial or clean filter media, hence the water flow behaviour is changed.
The next subsections are divided based on the type of clogging modelled, namely clogging due to SS, bioclogging and both combined. The units of the described parameters are generalized using L as a length unit, M as a unit of mass and T as a unit of time. An overview on the discussed models is given in Table 6.

6.1. Clogging Due to SS

In the work of Hua et al. [66], a model to estimate the clogging time based on SS is described. The conceptual model considers the effluent concentration of SS as a function of the operation time which is influenced by the gradual reduction of the infiltration rate and the particle size distribution of SS. Clogging is described as pore volume reduction caused by the mass of SS accumulated as volume over time. A mass balance (Equation (1)) is used to determine the accumulated SS (Mss, M) within the pore space based on the flow rate Qin (L3 T−1) and the SS influent and effluent concentration (Cin, Cout, M L−3). The SS mass in the outflow ( M o u t , M) is assumed to relate to a certain particle size (d0), hence a fraction of the total mass (Equation (2)). The certain sized SS accumulated in the pores would be washed off when the flow rate is higher than that of the maximum tolerated impact speed. Therefore, M d 0 % is a function of the infiltration rate K (L T−1) (Equation (3)). The calculation of K (Equation (4)) is based on the Hagen–Poiseuille equation in combination with Darcy’s law and the porosity ε (-),
dM SS dt   =   Q in   ×   ( C in     C out )
M out   =   M ss   ×   M d 0 %
M d 0 %   =   f ( K )
K   =   ε ( ρ w g µ w ) ( d 0 2 32 )   ×   ϕ ( t )
where g is acceleration of gravity (L T−2), µ w is the viscosity of the fluid (M T−1 L−1) and ρ w is the fluid density (M L−3). The authors further assume that K is a function of the operational time ϕ ( t ) .
After n days of operation, the SS mass accumulated within the pore space (Mss,n) is calculated as follows (Equation (5)): clogging is described as pore volume reduction, hence when a certain mass is reached the pore space is used up and a fully clogged state is reached. This threshold condition is described in Equation (6),
M s s , n   =   M s s , n 1   +   Q i n , n   ×   C i n   ×   { 1     f [ ε ( ρ w g µ w ) ( d 0 2 32 )   ×   ϕ ( t ) ] }
M s s , n ρ   ×   ( 1     ω )     ε   × h c   ×   A
where ε   ×   h c describes the total pore space per unit area (L), A (L2) is the filter area, h c is the height of the clogging layer (L) and ρ and ω are the density (M L−3) and moisture content (%) of SS, respectively. Values for the last named parameters are not provided by the authors.
Sani et al. [15] present a model to describe the impact of SS on clogging. Three mechanisms are incorporated which affect particle transport in a filter, namely diffusion, sedimentation and particle adsorption. Settling and aggregation mechanisms are described in Equation (7),
φ i t = D 2 φ i z 2 ( u v i ) φ i z ± ψ i + q ( z ) A φ i , i n
where φ i is the SS concentration with particle sizes in the range of i (M L−1), D is the dispersion coefficient (L2 T−1), z is the depth (L), u is the vertically flowing water velocity (L T−1), vi is the fall velocity or settling velocity of SS of particle size i (L T−1), ψ i is the source or sink term of the SS of particle size i (M L−3) and is used to take account of the effect of the aggregation or break-up of particles, q(z) is the lateral inflow to the wetland (L3 T−1), A is the wetland area (L−2) and φ i , i n is the concentration of the SS of size i in the lateral flow (M L−3). For the described continuous flow system, Darcy’s law (Equation (8)) is used to model the water flow where u is calculated as,
u   = K H z
where K is the hydraulic conductivity (L T−1), H is the water head (L) and z is the depth (L). Within Equation (7), lateral flow is neglected and effects of aggregation and breakup of SS are reflected by the dispersion coefficient and the settling velocity. A modified mass conservation model is given in Equation (9), where R is the source or sink term of SS,
φ i t = D 2 φ i z 2 ( u v i ) φ i z + R
Four sub-models (Equations (10)–(13)) provide the needed parameters for Equation (9). Two mechanisms can be described using dispersion, namely molecular diffusion and mechanical dispersion [71]. In this model, only mechanical dispersion is respected (Equation (10)),
D m d = α × u
where D m d is the mechanical dispersion (L2 T−1), α is the dispersivity (L) and u is the velocity based on Equation (8). The settling velocity for SS is represented using Equations (11) and (12) [72].
v i = w 0 i × f i
f i = ( 1 φ ) n
where w 0 i represents the terminal settling velocity of an isolated particle of size i, f i is the hindered settling velocity, φ is the total particle concentration (M L−3) and n is an empirical parameter given the value of 5.1 [72] as this represents the physical properties of the presented work [15]. The parameter w 0 i is an adjustable parameter representing the average effect of the sedimentation velocity of varied particle sizes. The last equation to solve Equation (9) is the description of the source or sink term R (Equation (13)). As attached BF supports SS accumulation, a form of the Monod equation [73] is used to describe particle adsorption.
R = M b s s × q m × φ φ S + φ
where M b s s is the adsorbed SS concentration on the substrate grains, q m is the maximum growth rate (T−1) and φ is the total SS concentration and φ S represents the half-saturation-coefficient (M L−3) of the Monod equation. The reduction of the hydraulic conductivity based on the decreasing pore space is taken into account using the Kozeny–Carman Equation (14) [74],
K = K 0 [ ( 1 + p D v t o t ε 0 ) x ( 1 D v t o t ε 0 ) y ]
where D v t o t is the total volumetric specific deposit (L3 L−3), K0 is the hydraulic conductivity (M T−1) and ε 0 the porosity (-) of the clean filter respectively and p, x and y describe empirical parameters (-).

6.2. Bioclogging

Mostafa and van Geel [75] present a microscopic conceptual model for bioclogging in unsaturated conditions. The voids in the unsaturated soil are simulated as capillary tubes with different diameters, where each diameter represents a volume of void space which can be filled at a certain capillary pressure. The capillary pressure–saturation relationship (Equation (15)) [76] and the permeability equation [77,78] are used to determine the relationship of the flow factors with and without BF. The relationship between capillary pressure and saturation is described by,
P c = ( S e 1 m 1 ) 1 n α
where P c is the capillary pressure (L) and α (L−1), m and n are van Genuchten model parameters [76]. The radii of the capillary tube can be determined by (Equation (16)),
r i = 2 σ c o s β P c , i
where r is the radii of the capillary tube (L), σ is the water surface tension (M T−2) and β is the contact angle, which is generally set to 0. The bundle of capillary tubes representing the porous media is divided into a user specific number N. The thickness of the BF t h (L) is calculated using Equation (17) where the microbial volume V m (L3) is an input parameter. Equation (18) defines the relationship between the clean filter media and the impact of the BF thickness on the media. In a last step, the hydraulic conductivity (L T−1) is updated based on the BF thickness in Equation (19). Growth and decay of BF is described using the Monod [73] equation similar to Equation (20).
V m = i = 1 N [ N i π r i 2 N i π ( r i t h ) 2 ]
U = i = 1 N N i π ( r i t h ) 4 i = 1 N N i π r i 4
K ( S e ) = K s a t × S e [ 1 ( 1 S e 1 m ) ] 2 × U
Brovelli et al. [79] present a macro scale model describing the effect of BF growth on the hydraulic properties of saturated porous media. In their work, only BF in the immobile phase is considered. BF growth (Equation (20)) and decay (Equation (21)) are described using a Monod equation,
µ = µ m a x × C e a K e a + C e a × C S K S + C S
d r = k d × X
where µ is the growth rate (T−1), X is the BF concentration (M L−3), µ m a x is the maximum growth rate (T−1), C is the concentration of the substrate (M L−3) and K is the half saturation constant (M L−3). The subscripts ea and s stand for electron acceptor and substrate respectively. d r represents the lysis rate with k d as the first-order decay constant (T−1). As within a porous media, the space for growth of BF is limited. Brovelli et al. [79] introduced a self-limiting function (Equation (22)) [80,81],
I b i o = X m a x ( X S + X a ) X m a x
where X m a x is the maximum BF content of the porous media per mass of solids and X S and X a represent the current amount of mobile and immobile BF, respectively. I b i o is defined by an upper bound with a maximum near the porosity but generally should be smaller as with the increasing BF growth, nutrient transport becomes a diffusion-controlled process. To account for BF attachment and detachment, the authors [79] adopted the classical deep-bed filtration. Thereby, the attachment coefficient (T−1) (Equation (23)) is described as [82,83],
k a t t = 3 ( 1 n ) v P η 2 d g
where v P is the pore velocity (L T−1), d g is the characteristic grain diameter (L) and η is the collector efficiency representing the frequency of collisions between mobile BF and grain surface (-). The collector efficiency is used by the authors as fitting parameter during model calibration as its general determination is difficult to assess [83]. Detachment due to shear forces is well discussed within the literature providing different approaches from dismissal [84] to simplified laws independent from the flow velocity. Here, a semi-empirical equation proposed by Rittmann [85] (Equation (24)) is used,
k d e t = c d ( γ v p ( 1 n ) 3 d g 2 n 3 M ) 0.58
where k d e t represents the detachment rate (T−1), γ is the viscosity of water (M L−1 T−1), M is the specific surface area (L2) and c d is an empirical parameter which is dependent on the experimental setup and calibrated for the work of Brovelli et al. [79]. Rittmann [85] proposes a value of 2.29 × 10−6. When combining Equations (20) to (24) one gets the variation of BF over time,
X s t = µ s I b i o X s k d X s k d e t X s k a t t X a
X a t = µ a I b i o X a k d X a k d e t X s k a t t X a
where the subscripts s and a refer to the immobile and mobile BF, respectively. As BF consists of up to 95% water, their density is assumed equal to that of water. Furthermore, the fraction of soil grains and water phase remain constant over time, hence solid BF is expressed as pore-fluid concentration instead of concentration per unit of soil [79]. The porosity of the filter media and its reduction due to BF as calculated in Equations (27) and (28) are described as follows,
n = n 0 n b i o
where n is the current porosity, n 0 is the porosity of the clean material and n b i o the porosity when occupied by BF. The latter represents the sum of different components of the BF, namely multiple bacteria strains, EPS and macromolecules. Therefore, n b i o is calculated as,
n b i o = i = 1 n X S i ρ b ρ S i
where X S i represents the dry weight of the ith component of the immobile BF (M), ρ b is the bulk density of the clean porous media (M L−3) and ρ S i is the density of the ith component of the immobile BF (M L−3). The change of the hydraulic conductivity based on the reduction of the pore space has been investigated based on experiments providing an exponential relationship (Equation (29)) [86,87],
K r e l ( n r e l ) = n r e l p ,   p > 0 K r e l = K K 0 , n r e l = n n 0
where K r e l and n r e l are the relative hydraulic conductivity (L T−1) and porosity (-), respectively, the subscript 0 represents the value of the clean media and, based on fitted data, p = 19 / 6 . The model based on Equation (29) assumes that BF clogs bigger pores first [86]. Another model describing the reduction of the hydraulic conductivity is based on pore network simulations where the assumption is in contrast with the description above stating that BF growth occurs first in the smaller pores [88]. This so called colonies model (Equation (30)) is described as
K r e l ( n r e l ) = a × ( n r e l n r e l 0 1 n r e l 0 ) + ( 1 a ) × ( n r e l n r e l 0 1 n r e l 0 ) 2
where a and n r e l 0 are adjustable parameters and 1 n r e l 0 being interpreted as the relative volume of BF needed to get a maximum reduction of the hydraulic conductivity. This maximum reduction occurs when n r e l approaches n r e l 0 . Good agreement on the model fitting is found in the literature [84] with values between 0.4–0.9 for n r e l 0 and −1 to −1.9 for a . A third model presented assumes a single, connected layer of BF covering the wall of each pore leading to a reduction of the pore radius [88]. Equation (31) presents the relationship of the hydraulic conductivity and the porosity as follows,
K r e l ( n r e l ) = [ ( n r e l n r e l 0 1 n r e l 0 ) b + K m i n ]   ×   1 1 + K m i n
where K m i n is the lower limit of the hydraulic conductivity. A good fit of the model was reached using parameters n r e l 0 in the range of 0.2–0.4 and 6 × 10−3–6 × 10−2 for K m i n . Within the model, the update of the hydraulic conductivity is carried out using the presented Equations (29)–(31) but no detail on which is actually used, nor a comparison of modelling results, is given.
Rajabzadeh et al. [62] implemented a computational fluid dynamics (CFD) model in COMSOL, accounting for spatial and temporal dynamics in VF wetlands. The model combines five sub-models, namely a fluid transport model, a solute transport model, a biokinetic model, a BF detachment model and a clogging model. The local porosity is estimated based on a BF model considering two mechanisms, namely growth of BF due to organic pollutants and the effect of fluid shear stress on local BF detachment. Darcy’s equation provides a linear relationship between hydraulic gradient and the flow velocity in porous media, ignoring the viscosity as well as inertia effects of the fluid flow. For the coarse media used (pea gravel) (Table 5), the relationship between velocity and hydraulic gradient can be non-linear [89], hence the Brinkman model, which extends Darcy’s law, is used (Equations (32) and (33), combining the continuity equation and momentum balance equation,
t ( ε p ρ ) + × ( ρ u ) = Q b r
ρ ε p ( u t + ( u × ) u ε p ) = p + [ 1 ε p µ ( u + ( u ) T 2 3 µ ( u ) I ] ( µ K + Q b r ) u + F
where µ is the dynamic viscosity (M L−1 T−1), u is the velocity vector (L T−1), ρ is the density of the fluid (M L−3), p is the pressure (Pa), ε p is the porosity (-), K is the permeability (L2) and Q b r is defined a mass source/sink (M L−3 T−1). The influence of volume forces such as gravity is respected by the force term F (M L−2 T−2). As the column experiment was carried out in saturated conditions, a single phase flow model was used. The solute transport is described by a simple advection-diffusion Equation (34), incorporating the biokinetic model within the source/sink term R ,
c t = ( D c ) ( u × c ) + R
where c is the solute concentration (M L−3), D is the molecular diffusivity (L2 T−1) and R is the reaction rate for the solute (M L−3 T−1) derived from a biokinetic model. The biokinetic model used in cooperates the processes of mineralization, hydrolysis, growth and lysis of heterotrophic bacteria (XH), respectively and is based on the work of Langergraber and Simunek [90]. For OM, three fractions are respected, namely, readily, slowly and inter OM. Hydrolysis of slowly to readily available OM is described in Equation (35), while the growth and decay of XH on readily OM is described by Equation (18),
d S S t   =   k h   ×   C S X H k x + C S X H   ×   X H
where k h and k x are rate constants for hydrolysis (T−1) and saturation/inhibition of hydrolysis (M M−1), S S is the readily available OM concentration (M L−3) and C S is slowly available OM (M L−3).
d X H t   =   µ m a x , H   ×   S S K S + S S   ×   X H     b i n a , H   ×   X H
where µ m a x , H ist the maximum aerobic growth rate (T−1), b i n a , H is the lysis rate (T−1) and K S is the half saturation coefficient (M L−3). The removal of S S is further described by Equation (37),
d S S t   =   1 Y H µ m a x , H   ×   S S K S + S S   ×   X H
where Y H is the yield coefficient representing the stoichiometric link between BF growth and OM consumption and is defined as BF production by OM consumption (M M−1). The next sub-model describes the BF detachment rate [85] implement in this model (Equation (38)). BF detachment is caused by fluid shear stress and is based on the flow rate, and therefore increases the local porosity.
b d e t , H   =   2.29   × 10 6 ( µ u ( 1     ε p ) 3 d p 2 ε p 3 M ) 0.58
where b d e t , H (T−1) is the BF detachment rate, d p is the average spherical gravel diameter (L) and M is the specific surface area of the filter media (L−1). The parameter M is determined by using a uniform gravel size with a spherical shape.
The fifth model, namely the clogging model, describes the hydrodynamic change caused by BF development within the pore space. Therefore, the local porosity at each time step is calculated based on the estimated BF concentration (Equation (39)),
ε ( x , y , t )   =   ε 0 ( 1 X H ( x , y , t ) ρ B i o f i l m )
where ε ( x , y , t ) is the calculated local porosity (-), ε 0 is the initial porosity (-) and ρ B i o f i l m is the BF density (M L−3). With this information, the relative change in local porosity and its effect on the permeability is calculated using the Kozeney–Carman Equation (40) [91] and is used to update the velocity and pressure profile (Equations (32) and (33),
K ( x , y , t )   =   K t = 0 ( ε ( x , y , t ) ε 0 ) 3 ( 1     ε ( x , y , t ) 1     ε 0 ) 2
A parameter efficient bioclogging model is presented Hua et al. [92]. Water flow is described by using the Darcy Equation (42) and the Kozeny–Carman Equation (49) [93] for the relationship between hydraulic conductivity and porosity. The change of BF accumulated within the pore space is described by an advection-reaction (41). Only BOD as substrate for BF growth is respected.
C t = u C x     r
where C is the substrate concentration (M L−3) and u is the pore-water velocity (L T−1) described as
u i , j   =   K i , j   ×   i i , j n i , j
where i is the hydraulic gradient (-), K is the hydraulic conductivity (L T−1) and n is the porosity (-). The subscripts i and j represent the time index and the space grid index, respectively. The reaction term r is represented by the following equation,
C t   =   k   ×   C   ×   B M
where k is the degradation rate constant (T−1). The BF growth is described by
B M t   =   C t     k d   ×   B F
where k d is the lysis rate of BF (T−1). Here the authors claim using the Monod equation, but BF growth is only related to consumption of substrate as well as lysis. The substrate concentration within the micro region for the time step Δ t is described as
C i , j =   C i 1 , j e x p [ λ ( x i x i 1 ) 2   ×   ( 1 u i 1 , j )   +   1 u i , j ]
where λ represents the deposition coefficient (T−1) of BF and is determined empirically by Kretzschmar et al. [94] as,
λ   =   v L l n ( C 0 C )
where v is the flow velocity (L.T−1), L is the column length (L) and C 0 and C represent the influent and effluent BOD concentration (M.L−3), respectively. The total amount of BF deposition σ i , j (M L−3) is calculated as
σ i , j = j = 1 j λ C i , j Δ t
The change in porosity based on deposited BF and follows
n i , j + 1 = n 0 n i , j × σ i , j ρ s
where ρ s is the BF density (M L−3) and n 0 is the initial porosity. The update of the hydraulic conductivity K i , j + 1 is finally carried out using the Kozeny–Carman equation
K i , j + 1 = K i , j × ( n i , j + 1 n i , j ) 3 × ( 1 n i , j + 1 1 n i , j ) 2

6.3. Bioclogging and SS Clogging

With FITOVERT [95], a model is presented to model unsaturated water flow, transport of dissolved as well as particulate components and a biodegradation model. To describe clogging, a model is formulated for the reduction of the porosity due to BM growth and accumulation of particulate matter. The one dimensional model formulation was developed using the Matlab® (Natick, MA, US) environment. Within the model domain, different layers can be implemented which are described by van Genuchten–Mualem SHP [76,78]. The variable saturated water flow is described by the Richards Equation (50) [96],
θ t = z [ K   ( h z 1 ) ]
where θ is the volumetric water content (L3 L−3), t is the time (T), z is the spatial coordinate (L), K is the unsaturated hydraulic conductivity (L T−1) and h is the matrix potential (L). The relationship between the pressure head, hydraulic conductivity and water content are described by the van Genuchten–Mualem model (Equations (51)),
θ ( h ) = { θ r + θ s θ r [ 1 + | α × h | n ] m h < 0 θ s h 0
K ( h ) = K S × S e l × [ 1 ( 1 S e 1 m ) m ] 2 , m = 1 1 n
S e = θ θ r θ s θ r
where θ r is the residual and θ s the saturated water content (L3 L−3) as well as K S are the saturated hydraulic conductivity (L T−1). The empirical parameters α (L−1) m and n influence the shape of the functions θ (h) and K(h), and l is defined as the pore-connectivity parameter. Se equals the effective water content as shown in Equation (51a). Transport of dissolved components is described using Bresler’s equation (Equations (54)) [97],
θ × C t = z ( d × C z ) q × C z + R × θ
where C is the concentration of a single soluble component in the liquid phase (M L−3), d is the dispersion coefficient (L2 T−1), q is the specific flow rate (L T−1) and R is the reaction term. The dispersion coefficient accounts for diffusion and mechanical dispersion, but the authors assume that in the liquid phase the effect of diffusion compared to dispersion can be neglected. Under saturated conditions and constant flow, the mechanical dispersion is considered to be proportional to the average flow velocity and dependent on the dispersivity λ (L),
d = λ × q θ
The authors assume further that the same relationship is true for unsaturated conditions. The reaction term R is determined by the biokinetic model implement based on the Activated Sludge Model 1 (ASM1) [98] describing the degradation of 13 components, seven of which are dissolved and six are particulate. Detailed information is missing within the original publication. The transport and filtration of particulate components is described using the sand filtration process [99],
q T m z   =   +   θ T m t + q × f × T m
where T m is the concentration of each single particulate component in water (M L−3) and f is the filter coefficient (L−1). The reduction of the porosity is related to the total volumetric specific deposit D v t o t (L3 L−3), which is continuously updated. How this parameter is determined is not further described. The effect of the change in porosity on the hydraulic conductivity is described using the Kozeny–Carman Equation (57) [100],
K = K 0 [ ( 1 + p D v t o t ε 0 ) x ( 1 D v t o t ε 0 ) y ]
where K0 is the hydraulic conductivity (L T−1), ε 0 is the porosity (-) of the clean filter respectively and p, x and y describe empirical parameters (-).
Hua et al. [101] developed a clogging model describing the pore volume reduction based on BF growth, SS retention and plant detritus. The model calculates the mass of each substance contributing to clogging and converts it to the volume reduction within the pore space. The overall application is the calculation of the operational time until the porosity is near zero, hence full clogging is reached. Influent parameters are defined as inert SS and BOD, as well as plant roots. SS and BOD are considered within one sub-model and plant root detritus in another sub-model. The model is based on the mass balance of the named contributors,
d M d t = F i n F o u t + S
where M is the total solid mass (M), F i n / F o u t represent the influent and effluent mass (M T−1) and S the source/sink term (M T−1) representing reaction and conversation rates. This is used as analogues for the influent parameters SS ( M I S ) and BF solids ( M B S ), contributing to SS clogging and bioclogging, respectively. The accumulation of total solid mass is represented by Equation (59). Two source terms represent the contributing loads from either BF development based on BOD ( S B S , (Equations (60)) and from inert matter production ( S I S , (Equations (61)),
d M T S d t = d M I S d t + d M B S d t
S B S = Q i n × C B O D × Y H 1 + K d × θ
S I S = f p × K d × θ × S B S
where Q i n is the influent flow rate (L3 T−1), t is the time (T), C B O D the BOD influent concentration (M L−3), Y H is the observed yield for heterotrophic BF (-), K d is the heterotrophic microbial endogenous decay coefficient (T−1), θ is the mean residence time of biosolids within the system (T) and f p is the fraction of microbial BM converted to inert matter (-).
While F i n , B S is represented by the BOD influent concentration, F i n , I S is calculated using Equation (62). The effluent mass of each parameter, represented by F o u t is assumed to be proportional to their respective influent mass (Equations (63)),
F i n , I S = Q i n × C S S × ( 1 f v + f v × f n v )
F o u t = F i n × λ
where C S S is the SS influent concentration (M L−3), f v is the proportion of organic matter in SS (-), f n v is the proportion of the inert matter in the organic SS (-) and λ is a proportional factor with a value between 0 and 1. The authors are aware of the simplification representing BF detachment and SS retention. In addition, Hua et al. [101] introduce a model predicting the effect of plant roots on clogging depending on two seasons, namely the growing season (Equations (64)) and the non-growing season (Equations (65)). During the growing season plant BM is increasing while during the non-growing season plant decay is prevailing,
G r o w t h   s e a s o n :   d M p l a n t , L d t = M p l a n t , L × ( f N H + f N O b P )
N o n   g r o w t h   s e a s o n :   d M p l a n t , L d t = M p l a n t , L × b P
where M p l a n t , L is the living plant BM (M), b P (T−1) is the plant decay coefficient for living plant material, f N H is the plant growth rate depending on ammonia and f N O is the plant growth rate depending on nitrate (T−1). Those growth rates are calculated as follows,
f N H = k p l × ( C N H K p N H + C N H )
f N O = k p l × ( C N O K p N O + C N O ) × ( K p N H K p N H + C N H )
where k p l is the relative plant growth rate (T−1), K p N H and K p N O (M L−3) are the half saturation coefficients for ammonia and nitrogen respectively and C N H and C N O are the influent concentrations (g m−3) for ammonia and nitrogen respectively. Due to plant decay, dead plant BM is produced. The change of this mass over time is calculated using Equation (68),
d M p l a n t , D d t = M p l a n t , L × b p M p l a n t , D × k d e g
where M p l a n t , D represents the dead plant BM (M) and k d e g is the first order decay rate (T−1) representing the loss of BM due to physical degradation processes such as physical degradation and invertebrate consumption. The total plant BM for the growth season (Equations (69)) and non-growth season (Equations (70)) is calculated by combining Equations (64) and (65) with Equation (68).
G r o w t h   s e a s o n :   d M p l a n t d t = M p l a n t , L × ( f N H + f N O ) M p l a n t , D × k d e g
N o n   g r o w t h   s e a s o n :   d M p l a n t d t = M p l a n t , D × k d e g
The underground portion of the plant mass is then computed as follows,
M r o o t = M p l a n t × η
where η is the root shoot ratio (-). In a final step, the models are coupled to predict the total volume occupied by the total solids M T S (Equations (59)) and the root mass M r o o t (Equations (71)),
V T = M T S ρ S × ( 1 ω ) + M r o o t ρ r o o t
where ρ S and ω are the density (M L−3) and the moisture content (-) of the total solids and ρ r o o t is the density of the plant roots (M L−3). The time after full clogging is reached is determined over the operational time of the model run and the following criteria,
V T ε × h c × A
where ε is the porosity (-), h c (m) is the depth of the filtration layer and A is the surface area (L2).

7. Summary and Conclusions

7.1. Experimental Studies

Throughout the literature, not many practical definitions on the extent of clogging are given. The authors mainly use the expression “clogged” to escribe a fully clogged system. Only one study evaluated different stages of clogging. Based on the observations of 21 VF wetlands, three categories of clogging were defined as the extent of permanent surface ponding between two loadings [13] (Table 7).
As presented, most studies are well documented, but usually some information is missing and full comparability is not given (Table 4 and Table 5). Table 8 presents a recommendation on data necessary to fully describe an experimental setup of VF wetlands, not only for studies on clogging.
In most experimental studies, the system is overloaded on purpose to get a fully clogged system in a short period of time for investigation. From studies investigating both types of clogging it can be concluded that SS is the main contributing factor to clogging. Only the study by Zhou et al. [64] comes to a different conclusion. In this study, bioclogging occurred faster within the glucose fed column than in the column investigating clogging due to SS. This can be explained by the rather high OLR (163 g COD m−2 day−1) used compared to other studies (Table 4). Generally, it can be concluded that
  • A general definition of clogging in various states is missing.
  • The majority of the experimental studies in the literature used continuous loading to fast track clogging of the system.
  • Besides analyses of occurring processes, the impact of resting on the recovery is observed and its importance pointed out. This strongly supports the importance of the intermittent loading operation used in the design guidelines [9].
  • For future studies on clogging processes, using real life operation is highly recommended to close the gap between experimental studies and implemented systems.

7.2. Modelling

As stated in Table 6 several models are set up as black box models, which provide good results but are only valid for the system under investigation. Fitted parameters cannot not be transferred to other systems, hence will need new fitting based on measured data. In order to get a deeper understanding of the underlying processes, grey box or white box models need to be developed further.
Existing water flow models follow two main approaches. Often, Darcy’s law is used which is a simplification when describing water flow in unsaturated conditions as the hydraulic conductivity is depended on the water content. This is respected when using the Richards equation in combination with a soil hydraulic model such as the one by van Genuchten–Mualem [76,102]. Here it should be mentioned, that by using this particular model only uniform or matrix flow is respected. Preferential flow, when occurring, cannot be described [103,104,105].
For updating the SHP of the filter material the semi empirical Kozeny–Carman equation, providing a porosity-permeability relationship is used regularly [62,95]. This equation which is based on a power function of the effective saturation, gives a relationship between the permeability, the porosity and the grain size and assumes Darcian flow which describes laminar flow and low pore water velocity. For media used in VF wetlands, this assumption might not be true [93]. Another approach, introduced by Mostafa and van Geel [75], is based on the work of Mualem [78] describing and changing unsaturated flow parameters based on the soil water retention curve. By using the capillary model, an applicable concept for the representation of the pore space and its change due to BF development is available. Based on this conceptual model Soleimani et al. [106] set up a two-dimensional unsaturated flow and transport model. For describing the BF growth the Monod equation is used. Their results show that the time until a clogged state is reached depends on the SHP of the media described by the van Genuchten–Mualem model [76]. The most influential parameters are the saturated hydraulic conductivity, the residual water content and the parameter n. The work by Samso et al. [107,108] is mainly concerned with HF wetlands but also provides a model implementation for variable saturated porous media [109] based on the concept of Mostafa and van Geel [75] to describe the change of porosity and hydraulic conductivity. Water flow is described using the Richards equation and BM growth using the Monod equation. Results showed that the model calculates the same BF thickness in all pores, hence also within bigger pores where no water content is expected. Still, this can be also related to a residual water content in big pores [75] but should be further investigated. Generally, it can be concluded that
  • A variety of modelling concepts are available to describe clogging processes in VF wetlands, ranging from simple black box models to deterministic process based models describing the actual processes.
  • Thereby different contributing factors, alone or in combination, are included, namely clogging due to SS, bioclogging and plant detritus.
  • Interactions between those contributors are, up to now, not implemented in existing models.
  • The main sub-models used [69] follow the same concepts over all modelling studies despite different details in the implemented versions.
  • Only a small fraction of the sub-models describe water flow in unsaturated conditions.
  • Detailed information on model parameters is often missing.

Author Contributions

The contributions of the authors are as follows: conceptualization of the review, B.P. and G.L.; writing—original draft preparation, B.P.; data curation, B.P.; writing—review and editing, G.L.

Funding

This open access publication was supported by the Doctoral School “Human River Systems in the 21st Century (HR21)” of the University of Natural Resources and Life Sciences, Vienna.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BFBiofilmPEPerson equivalent
BMBiomassRHRelative humidity
CFDComputational fluid dynamicsSHPSoil hydraulic parameters
CODChemical oxygen demandSSSuspended solids
EPSExtracellular polymeric substancesTDRTime domain reflectometry
ETEvapotranspirationTSSTotal suspended solids
HFHorizontal flowTWTreatment wetland
HLRHydraulic loading rateVFVertical flow
OLROrganic loading rateVSSVolatile suspended solids
OMOrganic matterXHHeterotrophic bacteria

References

  1. Dotro, G.; Langergraber, G.; Molle, P.; Nivala, J.; Puigagut, J.; Stein, O.R.; von Sperling, M. Treatment Wetlands; Biological Wastewater Treatment Series; IWA Publishing: London, UK, 2017; Volume 7, p. 172. ISBN 9781780408774. Available online: http://wio.iwaponline.com/content/16/9781780408774 (accessed on 14 November 2019).
  2. Stefanakis, A.I. Introduction to constructed wetland technology. In Constructed Wetlands for Industrial Wastewater Treatment; Alexandros, S., Ed.; John Wiley & Sons, Ltd.: Chichester, UK, 2018; pp. 1–21. ISBN 9781119268413. [Google Scholar]
  3. Langergraber, G.; Weissenbacher, N. Survey on number and size distribution of treatment wetlands in Austria. Water Sci. Technol. 2017, 75, 2309–2315. [Google Scholar] [CrossRef] [PubMed]
  4. Seidel, K. Über die Selbstreinigung natüirlicher Gewässer. Naturwissenschaften 1967, 63, 286–291. [Google Scholar] [CrossRef] [PubMed]
  5. Haberl, R.; Perfler, R.; Mayer, H. Constructed wetlands in Europe. Water Sci. Technol. 1995, 32, 305–315. [Google Scholar] [CrossRef]
  6. Platzer, C. Design recommendations for subsurface flow constructed wetlands for nitrification and denitrification. Water Sci. Technol. 1999, 40, 257–263. [Google Scholar] [CrossRef]
  7. Brix, H.; Arias, C.A. The use of vertical flow constructed wetlands for on-site treatment of domestic wastewater: New Danish guidelines. Ecol. Eng. 2005, 25, 491–500. [Google Scholar] [CrossRef]
  8. Molle, P.; Liénard, A.; Boutin, C.; Merlin, G.; Iwema, A. How to treat raw sewage with constructed wetlands: An overview of the French systems. Water Sci. Technol. 2005, 51, 11–21. [Google Scholar] [CrossRef]
  9. Nivala, J.; van Afferden, M.; Hasselbach, R.; Langergraber, G.; Molle, P.; Rustige, H.; Nowak, J. The new German standard on constructed wetland systems for treatment of domestic and municipal wastewater. Water Sci. Technol. 2018, 78, 2414–2426. [Google Scholar] [CrossRef]
  10. Langergraber, G.; Leroch, K.; Pressl, A.; Sleytr, K.; Rohrhofer, R.; Haberl, R. High-rate nitrogen removal in a two-stage subsurface vertical flow constructed wetland. Desalination 2009, 246, 55–68. [Google Scholar] [CrossRef]
  11. Petitjean, A.; Forquet, N.; Boutin, C. Oxygen profile and clogging in vertical flow sand filters for on-site wastewater treatment. J. Environ. Manag. 2016, 170, 15–20. [Google Scholar] [CrossRef]
  12. Licciardello, F.; Aiello, R.; Alagna, V.; Iovino, M.; Ventura, D.; Cirelli, G.L. Assessment of clogging in constructed wetlands by saturated hydraulic conductivity measurements. Water Sci. Technol. 2019, 79, 314–322. [Google Scholar] [CrossRef]
  13. Winter, K.-J.; Goetz, D. The impact of sewage composition on the soil clogging phenomena of vertical flow constructed wetlands. Water Sci. Technol. 2003, 48, 9–14. [Google Scholar] [CrossRef]
  14. Langergraber, G.; Haberl, R.; Laber, J.; Pressl, A. Evaluation of substrate clogging processes in vertical flow constructed wetlands. Water Sci. Technol. 2003, 48, 25–34. [Google Scholar] [CrossRef]
  15. Sani, A.; Scholz, M.; Babatunde, A.; Wang, Y. Impact of water quality parameters on the clogging of vertical-flow constructed wetlands treating urban wastewater. Water Air Soil Pollut. 2013, 224, 1488. [Google Scholar] [CrossRef]
  16. Bryers, J.D. Biofilm formation and chemostat dynamics: Pure and mixed culture considerations. Biotechnol. Bioeng. 1984, 26, 948–958. [Google Scholar] [CrossRef]
  17. Caselles-Osorio, A.; Puigagut, J.; Segú, E.; Vaello, N.; Granés, F.; García, D.; García, J. Solids accumulation in six full-scale subsurface flow constructed wetlands. Water Res. 2007, 41, 1388–1398. [Google Scholar] [CrossRef]
  18. Yao, K.-M.; Habibian, M.T.; O’Melia, C.R. Water and waste water filtration. Concepts and applications. Environ. Sci. Technol. 1971, 5, 1105–1112. [Google Scholar] [CrossRef]
  19. McDowell-Boyer, L.M.; Hunt, J.R.; Sitar, N. Particle transport through porous media. Water Resour. Res. 1986, 22, 1901–1921. [Google Scholar] [CrossRef]
  20. Sen, T.K.; Khilar, K.C. Review on subsurface colloids and colloid-associated contaminant transport in saturated porous media. Adv. Colloid Interface Sci. 2006, 119, 71–96. [Google Scholar]
  21. Al-Isawi, R.; Scholz, M.; Wang, Y.; Sani, A. Clogging of vertical-flow constructed wetlands treating urban wastewater contaminated with a diesel spill. Environ. Sci. Pollut. Res. 2015, 22, 12779–12803. [Google Scholar] [CrossRef]
  22. Langergraber, G.; Pressl, A.; Leroch, K.; Rohrhofer, R.; Haberl, R. Comparison of single-stage and a two-stage vertical flow constructed wetland systems for different load scenarios. Water Sci. Technol. 2010, 61, 1341–1348. [Google Scholar] [CrossRef]
  23. Masi, F.; Bresciani, R.; Martinuzzi, N.; Cigarini, G.; Rizzo, A. Large scale application of French reed beds: Municipal wastewater treatment for a 20,000 inhabitant’s town in Moldova. Water Sci. Technol. 2017, 76, 134–146. [Google Scholar] [CrossRef] [PubMed]
  24. DWA-A262E. Principles for Dimensioning, Construction and Operation of Wastewater Treatment Plants with Planted and Unplanted Filters for Treatment of Domestic and Municipal Wastewater, DWA Set of Rules; German Association for Water, Wastewater and Waste (DWA): Hennef, Germany, 2017. [Google Scholar]
  25. Platzer, C.; Mauch, K. Soil clogging in vertical flow reed beds—Mechanisms, parameters, consequences and....solutions? Water Sci. Technol. 1997, 35, 175–181. [Google Scholar] [CrossRef]
  26. Martí, A.C.; Pucher, B.; Hernández-Crespo, C.; Monerris, M.M.; Langergraber, G. Numerical simulation of vertical flow wetlands with special emphasis on treatment performance during winter. Water Sci. Technol. 2018, 78, 2019–2026. [Google Scholar] [CrossRef] [PubMed]
  27. Langergraber, G.; Leroch, K.; Pressl, A.; Rohrhofer, R.; Haberl, R. A two-stage subsurface vertical flow constructed wetland for high-rate nitrogen removal. Water Sci. Technol. 2008, 57, 1881–1887. [Google Scholar] [CrossRef]
  28. Langergraber, G.; Pressl, A.; Leroch, K.; Rohrhofer, R.; Haberl, R. Experiences with a top layer of gravel to enhance the performance of vertical flow constructed wetlands at cold temperatures. Water Sci. Technol. 2009, 59, 1111–1116. [Google Scholar] [CrossRef]
  29. Hua, G.; Chen, Q.; Kong, J.; Li, M. Evapotranspiration versus oxygen intrusion: Which is the main force in alleviating bioclogging of vertical-flow constructed wetlands during a resting operation? Environ. Sci. Pollut. Res. 2017, 24, 18355–18362. [Google Scholar] [CrossRef]
  30. García, J.; Rousseau, D.P.L.; Morató, J.; Lesage, E.; Matamoros, V.; Bayona, J.M. Contaminant removal processes in subsurface-flow constructed wetlands: A review. Crit. Rev. Environ. Sci. Technol. 2010, 40, 561–661. [Google Scholar] [CrossRef]
  31. Yang, M.; Lu, M.; Sheng, L.; Wu, H. Study of the spatial and temporal distribution of accumulated solids in an experimental vertical-flow constructed wetland system. Sci. Total Environ. 2018, 509–516. [Google Scholar] [CrossRef]
  32. Kadlec, R.H.; Wallace, S.D. Treatment Wetlands, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2009; ISBN 9781566705264. [Google Scholar]
  33. Langergraber, G. Are constructed treatment wetlands sustainable sanitation solutions? Water Sci. Technol. 2013, 67, 2133–2140. [Google Scholar] [CrossRef]
  34. Brix, H. Do macrophytes play a role in constructed treatment wetlands? Water Sci. Technol. 1997, 35, 11–17. [Google Scholar] [CrossRef]
  35. Molle, P.; Liénard, A.; Grasmick, A.; Iwema, A. Effect of reeds and feeding operations on hydraulic behaviour of vertical flow constructed wetlands under hydraulic overloads. Water Res. 2006, 40, 606–612. [Google Scholar] [CrossRef]
  36. Scholz, M.; Xu, J. Performance comparison of experimental constructed wetlands with different filter media and macrophytes treating industrial wastewater contaminated with lead and copper. Bioresour. Technol. 2002, 83, 71–79. [Google Scholar] [CrossRef]
  37. Cooper, P. The performance of vertical flow constructed wetland systems with special reference to the significance of oxygen transfer and hydraulic loading rates. Water Sci. Technol. 2005, 51, 81–90. [Google Scholar] [CrossRef] [PubMed]
  38. Wallace, S.; Knight, R. Small-Scale Constructed Wetland Treatment Systems: Feasibility, Design Criteria, and O&M Requirements; Water Environment Research Foundation (WERF): Alexandria, VA, USA, 2006. [Google Scholar]
  39. Murphy, C.; Cooper, D. Reed bed refurbishment: A sustainable approach. In Proceedings of the 3rd International Symposium on “Wetland Pollutant Dynamics and Control—WETPOL 2009”, Barcelona, Spain, 20–24 September 2009. [Google Scholar]
  40. Behrends, L.L.; Bailey, E.; Houke, L.; Jansen, P.; Smith, S. Non-invasive methods for treating and removing sludge from subsurface flow constructed wetlands I. In Proceedings of the Annual National Onsite Wastewater Recycling Association (NOWRA) Meeting, Baltimore, MD, USA, 12–14 March 2007. [Google Scholar]
  41. Behrends, L.L.; Bailey, E.; Houke, L.; Jansen, P.; Smith, S. Non-invasive methods for treating and removing sludge from subsurface flow constructed wetlands II. In Proceedings of the 10th International Conference on Wetland Systems for Water Pollution Control, Lisbon, Portugal, 23–29 September 2006; pp. 1271–1281. [Google Scholar]
  42. Nivala, J.; Rousseau, D.P.L. Reversing clogging in subsurface-flow constructed wetlands by hydrogen peroxide treatment: Two case studies. Water Sci. Technol. 2009, 59, 2037–2046. [Google Scholar] [CrossRef]
  43. Hua, G.F.; Wei, Z.; Lianfang, Z.; Yunhui, Z. Applying solubilization treatment to reverse clogging in laboratory-scale vertical flow constructed wetlands. Water Sci. Technol. 2010, 61, 1479–1487. [Google Scholar]
  44. Tang, P.; Xiang, Z.; Zhou, Y.; Zhang, Y. Enzyme treatment improves the performance of laboratory-scale vertical flow constructed wetland. Bioresour. Technol. 2018, 268, 665–671. [Google Scholar] [CrossRef]
  45. Davison, L.; Headley, T.; Pratt, K. Aspects of design, structure, performance and operation of reed beds—Eight years’ experience in northeastern New South Wales, Australia. Water Sci. Technol. 2005, 51, 129–138. [Google Scholar] [CrossRef]
  46. Ye, J.; Xu, Z.; Chen, H.; Wang, L.; Benoit, G. Reduction of clog matter in constructed wetlands by metabolism of Eisenia foetida: Process and modeling. Environ. Pollut. 2018, 238, 803–811. [Google Scholar] [CrossRef]
  47. Li, H.Z.; Wang, S.; Ye, J.F.; Xu, Z.X.; Jin, W. A practical method for the restoration of clogged rural vertical subsurface flow constructed wetlands for domestic wastewater treatment using earthworm. Water Sci. Technol. 2011, 63, 283–290. [Google Scholar] [CrossRef]
  48. Kim, D.-S.; Fogler, H.S. Biomass evolution in porous media and its effects on permeability under starvation conditions. Biotechnol. Bioeng. 2000, 69, 47–56. [Google Scholar] [CrossRef]
  49. Thullner, M. Comparison of bioclogging effects in saturated porous media within one- and two-dimensional flow systems. Ecol. Eng. 2010, 36, 176–196. [Google Scholar] [CrossRef]
  50. Hua, G.; Zeng, Y.; Zhao, Z.; Cheng, K.; Chen, G. Applying a resting operation to alleviate bioclogging in vertical flow constructed wetlands: An experimental lab evaluation. J. Environ. Manag. 2014, 136, 47–53. [Google Scholar] [CrossRef] [PubMed]
  51. Hua, G.; Kong, J.; Ji, Y.; Li, M. Influence of clogging and resting processes on flow patterns in vertical flow constructed wetlands. Sci. Total Environ. 2018, 621, 1142–1150. [Google Scholar] [CrossRef] [PubMed]
  52. Bouwer, H.; Lance, J.C.; Riggs, M.S. High-rate land treatment II: Water quality and economic aspects of the flushing meadows project. Water Pollut. Control Fed. 1974, 46, 844–859. [Google Scholar]
  53. Tanner, C.; Sukias, J.P. Accumulation of organic solids in gravel-bed constructed wetlands. Water Sci. Technol. 1995, 32, 229–239. [Google Scholar] [CrossRef]
  54. Chazarenc, F.; Merlin, G. Influence of surface layer on hydrology and biology of gravel bed vertical flow constructed wetlands. Water Sci. Technol. 2005, 51, 91–97. [Google Scholar] [CrossRef]
  55. Tanner, C.C.; Sukias, J.P.; Upsdell, M.P. Organic matter accumulation during maturation of gravel-bed constructed wetlands treating farm dairy wastewaters. Water Res. 1998, 32, 3046–3054. [Google Scholar] [CrossRef]
  56. Pedescoll, A.; Knowles, P.R.; Davies, P.; García, J.; Puigagut, J. A Comparison of in situ constant and falling head permeameter tests to assess the distribution of clogging within horizontal subsurface flow constructed wetlands. Water Air Soil Pollut. 2012, 223, 2263–2275. [Google Scholar] [CrossRef]
  57. Kim, B.; Forquet, N. Pore-scale observation of deposit within the gravel matrix of a vertical flow constructed wetland. Environ. Technol. 2016, 37, 3146–3150. [Google Scholar] [CrossRef] [Green Version]
  58. Martinez-Carvajal, G.D.; Oxarango, L.; Adrien, J.; Molle, P.; Forquet, N. Assessment of X-ray computed tomography to characterize filtering media from vertical flow treatment wetlands at the pore scale. Sci. Total Environ. 2019, 658, 178–188. [Google Scholar] [CrossRef]
  59. Zhao, L.; Zhu, W.; Tong, W. Clogging processes caused by biofilm growth and organic particle accumulation in lab-scale vertical flow constructed wetlands. J. Environ. Sci. 2009, 21, 750–757. [Google Scholar] [CrossRef]
  60. Hua, G.F.; Zhu, W.; Zhao, L.F.; Huang, J.Y. Clogging pattern in vertical-flow constructed wetlands: Insight from a laboratory study. J. Hazard. Mater. 2010, 180, 668–674. [Google Scholar] [CrossRef] [PubMed]
  61. Hua, G.F.; Zhao, Z.W.; Kong, J.; Guo, R.; Zeng, Y.T.; Zhao, L.F.; Zhu, Q.D. Effects of plant roots on the hydraulic performance during the clogging process in mesocosm vertical flow constructed wetlands. Environ. Sci. Pollut. Res. 2014, 21, 13017–13026. [Google Scholar] [CrossRef] [PubMed]
  62. Rajabzadeh, A.R.; Legge, R.L.; Weber, K.P. Multiphysics modelling of flow dynamics, biofilm development and wastewater treatment in a subsurface vertical flow constructed wetland mesocosm. Ecol. Eng. 2015, 74, 107–116. [Google Scholar] [CrossRef]
  63. Song, X.; Ding, Y.; Wang, Y.; Wang, W.; Wang, G.; Zhou, B. Comparative study of nitrogen removal and bio-film clogging for three filter media packing strategies in vertical flow constructed wetlands. Ecol. Eng. 2015, 74, 1–7. [Google Scholar] [CrossRef]
  64. Zhou, Y.; Luo, S.; Yu, B.; Zhang, T.; Li, J.; Zhang, Y. A comparative analysis for the development and recovery processes of different types of clogging in lab-scale vertical flow constructed wetlands. Environ. Sci. Pollut. Res. 2018, 25, 24073–24083. [Google Scholar] [CrossRef]
  65. ÖNORM B2505. Bepflanzte Bodenfilter (Pflanzenkläranlagen)—Anwendung, Bemessung, Bau und Betrieb; Österreichisches Normeringsinstitut: Vienna, Austria, 2009. (In German) [Google Scholar]
  66. Hua, G.F.; Zhu, W.; Zhang, Y.H. A conceptual approach based on suspended solids to estimate clogging time in constructed wetlands. J. Environ. Sci. Heal. Part A 2010, 45, 1519–1525. [Google Scholar] [CrossRef]
  67. Fujimaki, H.; Inoue, M. A flux-controlled steady-state evaporation method for determining unsaturated hydraulic conductivity at low matric pressure head values. Soil Sci. 2003, 168, 385–395. [Google Scholar] [CrossRef]
  68. Ren, Y.; Gong, Q.; Zhang, T.C.; Kang, J.; Liu, D. Wavy subsurface flow and vertical flow constructed wetlands: Effects of aeration and wastewater composition on wetland clogging. J. Environ. Eng. 2015, 141, 04015018. [Google Scholar] [CrossRef]
  69. Langergraber, G. Applying process-based models for subsurface flow treatment wetlands: Recent developments and challenges. Water 2017, 9, 18. [Google Scholar] [CrossRef]
  70. Grandjean, T.; McGordon, A.; Jennings, P. Structural identifiability of equivalent circuit models for li-ion batteries. Energies 2017, 10, 90. [Google Scholar] [CrossRef] [Green Version]
  71. Rowe, R.K.; Babcock, D. Modelling the clogging of coarse gravel and tire shreds in column tests. Can. Geotech. J. 2007, 44, 1273–1285. [Google Scholar] [CrossRef] [Green Version]
  72. Richardson, J.F.; Zaki, W.N. The sedimentation of a suspension of uniform spheres under conditions of viscous flow. Chem. Eng. Sci. 1954, 3, 65–73. [Google Scholar] [CrossRef]
  73. Monod, J. The Growth of Bacterial Cultures. Annu. Rev. Microbiol. 1949, 3, 371–394. [Google Scholar] [CrossRef] [Green Version]
  74. Carman, P.C. Fluid flow through granular beds. Chem. Eng. Res. Des. 1997, 75, S32–S48. [Google Scholar] [CrossRef]
  75. Mostafa, M.; Van Geel, P.J. conceptual models and simulations for biological clogging in unsaturated soils. Vadose Zo. J. 2007, 6, 175. [Google Scholar] [CrossRef]
  76. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1. Soil Sci. Soc. Am. J. 1980, 44, 892. [Google Scholar] [CrossRef] [Green Version]
  77. Burdine, N.T. Relative permeability calculations from pore size distribution data. J. Pet. Technol. 1953, 5, 71–78. [Google Scholar] [CrossRef]
  78. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef] [Green Version]
  79. Brovelli, A.; Malaguerra, F.; Barry, D.A. Bioclogging in porous media: Model development and sensitivity to initial conditions. Environ. Model. Softw. 2009, 24, 611–626. [Google Scholar] [CrossRef] [Green Version]
  80. Zysset, A.; Stauffer, F.; Dracos, T. Modeling of reactive groundwater transport governed by biodegradation. Water Resour. Res. 1994, 30, 2423–2434. [Google Scholar] [CrossRef]
  81. Kildsgaard, J.; Engesgaard, P. Numerical analysis of biological clogging in two-dimensional sand box experiments. J. Contam. Hydrol. 2001, 50, 261–285. [Google Scholar] [CrossRef]
  82. Harvey, R.W.; Garabedian, S.P. Use of colloid filtration theory in modeling movement of bacteria through a contaminated sandy aquifer. Environ. Sci. Technol. 1991, 25, 178–185. [Google Scholar] [CrossRef]
  83. Scheibe, T.D.; Dong, H.; Xie, Y. Correlation between bacterial attachment rate coefficients and hydraulic conductivity and its effect on field-scale bacterial transport. Adv. Water Resour. 2007, 30, 1571–1582. [Google Scholar] [CrossRef]
  84. Thullner, M.; Schroth, M.H.; Zeyer, J.; Kinzelbach, W. Modeling of a microbial growth experiment with bioclogging in a two-dimensional saturated porous media flow field. J. Contam. Hydrol. 2004, 70, 37–62. [Google Scholar] [CrossRef] [PubMed]
  85. Rittman, B.E. The effect of shear stress on biofilm loss rate. Biotechnol. Bioeng. 1982, 24, 501–506. [Google Scholar] [CrossRef] [PubMed]
  86. Clement, T.P.; Hooker, B.S.; Skeen, R.S. Macroscopic models for predicting changes in saturated porous media properties caused by microbial growth. Ground Water 1996, 34, 934–942. [Google Scholar] [CrossRef]
  87. Ives, K.J.; Pienvichitr, V. Kinetics of the filtration of dilute suspensions. Chem. Eng. Sci. 1965, 20, 965–973. [Google Scholar] [CrossRef]
  88. Thullner, M.; Zeyer, J.; Kinzelbach, W. Influence of microbial growth on hydraulic properties of pore networks. Transp. Porous Media 2002, 49, 99–122. [Google Scholar] [CrossRef]
  89. Reddy, N.B.; Rao, P.R. Effect of convergence on nonlinear flow in porous media. J. Hydraul. Eng. 2006, 132, 420–427. [Google Scholar] [CrossRef]
  90. Langergraber, G.; Šimůnek, J. Reactive transport modeling of subsurface flow constructed wetlands using the HYDRUS wetland module. Vadose Zone J. 2012, 11. [Google Scholar] [CrossRef] [Green Version]
  91. Blazejewski, R.; Murat-Blazejewska, S. Soil clogging phenomena in constructed wetlands with subsurface flow. Water Sci. Technol. 1997, 35, 183–188. [Google Scholar] [CrossRef]
  92. Hua, G.F.; Shao, C.; Cheng, Y.; Kong, J.; Zhao, Z. Parameter-efficient bioclogging model: Calibration and comparison with laboratory data. Environ. Sci. Pollut. Res. 2019, 26, 3731–3740. [Google Scholar] [CrossRef] [PubMed]
  93. Carrier, W.D. Goodbye, hazen; hello, kozeny-carman. J. Geotech. Geoenviron. Eng. 2003, 129, 1054–1056. [Google Scholar] [CrossRef]
  94. Kretzschmar, R.; Borkovec, M.; Grolimund, D.; Elimelech, M. Mobile subsurface colloids and their role in contaminant transport. In Advanced in Agronomy; Academic Press: Cambridge, MA, USA, 1999; Volume 66, pp. 121–193. ISBN 0065-2113. [Google Scholar]
  95. Giraldi, D.; de Michieli Vitturi, M.; Iannelli, R. FITOVERT: A dynamic numerical model of subsurface vertical flow constructed wetlands. Environ. Model. Softw. 2010, 25, 633–640. [Google Scholar] [CrossRef]
  96. Richards, L.A. Capillary conduction of liquids through porous mediums. Physics 1931, 1, 318–333. [Google Scholar] [CrossRef]
  97. Bresler, E. Simultaneous transport of solutes and water under transient unsaturated flow conditions. Water Resour. Res. 1973, 9, 975–986. [Google Scholar] [CrossRef]
  98. Henze, M.; Gujer, W.; Mino, T.; van Loosdrecht, M.C.M. Activated Sludge Models ASM1, ASM2, ASM2d and ASM3; IWA Publishing: London, UK, 2000; ISBN 1900222248. [Google Scholar]
  99. Iwasaki, T. Some Notes on Sand Filtration. J. Am. Water Works Assoc. 1937, 29, 1591–1597. [Google Scholar] [CrossRef]
  100. Boller, M.A.; Kavanaugh, M.C. Particle characteristics and headloss increase in granular media filtration. Water Res. 1995, 29, 1139–1149. [Google Scholar] [CrossRef]
  101. Hua, G.F.; Li, L.; Zhao, Y.Q.; Zhu, W.; Shen, J.Q. An integrated model of substrate clogging in vertical flow constructed wetlands. J. Environ. Manag. 2013, 119, 67–75. [Google Scholar] [CrossRef] [Green Version]
  102. Peszyńska, M.; Yi, S.Y. Numerical methods for unsaturated flow with dynamic capillary pressure in heterogeneous porous media. Int. J. Numer. Anal. Model. 2008, 5, 126–149. [Google Scholar]
  103. Morvannou, A.; Forquet, N.; Vanclooster, M.; Molle, P. Which hydraulic model to use for vertical flow constructed wetlands? In Proceedings of the 4th International Conference HYDRUS Software Applications to Subsurface Flow and Contaminant Transport Problems, Prague, Czech Republic, 21–22 March 2013; pp. 245–255. [Google Scholar]
  104. Forquet, N.; Morvannou, A. Preferential flows in treatment wetlands. In Proceedings of the 8th International Symposium on Wetland Pollutant Dynamics and Control—WETPOL 2019, Aarhus, Denmark, 17–21 June 2019; p. 219, ISBN 978-87-971486-0-0. [Google Scholar]
  105. Pucher, B.; Langergraber, G. Simulating vertical flow wetlands using filter media with different grain sizes with the HYDRUS Wetland Module. J. Hydrol. Hydromech. 2018, 66, 227–231. [Google Scholar] [CrossRef] [Green Version]
  106. Soleimani, S.; Van Geel, P.J.; Isgor, O.B.; Mostafa, M.B. Modeling of biological clogging in unsaturated porous media. J. Contam. Hydrol. 2009, 106, 39–50. [Google Scholar] [CrossRef] [PubMed]
  107. Samsó, R.; Garcia, J. BIO_PORE, a mathematical model to simulate biofilm growth and water quality improvement in porous media: Application and calibration for constructed wetlands. Ecol. Eng. 2013, 54, 116–127. [Google Scholar] [CrossRef]
  108. Samsó, R.; Blázquez, J.; Agulló, N.; Grau, J.; Torres, R.; García, J. Effect of bacteria density and accumulated inert solids on the effluent pollutant concentrations predicted by the constructed wetlands model BIO_PORE. Ecol. Eng. 2015, 80, 172–180. [Google Scholar] [CrossRef] [Green Version]
  109. Samsó, R.; García, J.; Molle, P.; Forquet, N. Modelling bioclogging in variably saturated porous media and the interactions between surface/subsurface flows: Application to constructed wetlands. J. Environ. Manag. 2016, 165, 271–279. [Google Scholar] [CrossRef] [Green Version]
Table 1. Comparison of design parameters of VF wetland types to prevent clogging (adapted from [1]).
Table 1. Comparison of design parameters of VF wetland types to prevent clogging (adapted from [1]).
TypeSurface Area Per PEPre-TreatmentMediaOLRDosing
Single stage VF4 m2YesSand0–2 mm20>6 h
4 m2YesSand0–4 mm20–27>6 h
Two stage VF2 m2YesSand1–4 mm/0–4 mm80>3 h
French VF 2–2.5 m2NoFine gravel2–8 mm/0–4 mm100–3503.5 day
7 day rest
Table 2. Available remediation strategies for clogged wetlands.
Table 2. Available remediation strategies for clogged wetlands.
CategoryMethodTime Until Recovery
Destructive methodsExcavation and replacement using new mediaWithin one day
Several days
Excavation, washing and reuse of media
Active treatmentApplication of oxidising agentImmediate up to several hours
Addition of solubilisation agentOne week
Enzyme treatment-
Addition of earthwormsAround 10 days
Passive treatmentResting of the whole treatment wetland (TW)10–20 days
Table 3. Mixtures of synthetic wastewater used for investigation.
Table 3. Mixtures of synthetic wastewater used for investigation.
StudyBiocloggingSS CloggingNutrients Added
Zhao et al.[59]Glucose StarchKNO3, K2HPO4
Hua et al.[60]-Lake sediment-
Hua et al.[61]Glucose-(NH4)2SO4, CO(NH2)2, K2HPO4,
Rajabzadeh et al.[62]Molasses, urea-NH4H2PO4
Song et al.[63]not mentioned
Yang et al.[31]-Zeolite powder-
Zhou et al.[64]GlucoseStarchNaCl, (NH4)2SO4, CO(NH2)2, K2HPO4, NaH2PO4·12H2O, MgSO4·7H20, CaCl2
Table 4. Operational description of experimental studies.
Table 4. Operational description of experimental studies.
Authors Feeding Wastewater O L R g C O D   m 2   d a y 1 T S S g T S S   m 2   d a y 1 H L R L   m 2   d a y 1 COD mg   L 1 TSS mg   L 1
Langergraber et al.[14]intermittentevery 6 hreal5–100 *8–2940–150100–700 *50–300 *
5–120 *8–4840–250100–700 *50–300 *
Zhao et al.[59]continuously synthetic40–3434085047–404
Hua et al.[60]continuously Lake sediments 300500
300500
300500
Sani et al.[15]continuously72 h/36 h feedingreal222373301
48 h restreal diluted118151
Hua et al.[50]continuously synthetic300 (as BOD) 500
Rajabzadeh et al.[62]continuously synthetic 500
Song et al.[63]fill and drain 10 h HRT 2 h draining synthetic 112
236
113
Ren et al.[68]continuously real66288008235
Petijean et al.[11]intermittentevery 6 or 8 hreal35973483119
Yang et al.[31]continuously synthetic 18.537 500
Zhou et al.[64]continuously synthetic163137815200168
6884
* unpublished data.
Table 5. Description on the filter media used in the described experiments.
Table 5. Description on the filter media used in the described experiments.
Authors Media S i z e m m Porosity k f c m   s 1 d 10 m m d 60 m m Cu
Langergraber et al.[14]Sand0.063–4 0.1317.69
1–4 1.132.73
Zhao et al.[59]Coarse sand 0.364.85 × 10−214.44.4
Hua et al.[60]Gravel30.34 2.55.82.32
100.44 611.41.9
200.47 15251.67
Sani et al.[15]Pea gravel10
20
Hua et al.[50]Coarse sand 0.120.27842.32
Rajabzadeh et al.[62]Pea gravel10
Song et al.[63]Sand, 4 layers3–8
Sand, 4 layers8–3
Uniform5–6
Ren et al.[68]Sand
Petitjean et al.[11]Sand 0–4 0.160.0480.3
Yang et al.[31]Gravel3–40.370.458
Zhou et al.[64]Sand1–2
Table 6. Overview on the modelling studies investigated.
Table 6. Overview on the modelling studies investigated.
ReferenceDescriptionCloggingModel TypeSHP UpdateSub-Model
Hua et al.[66]Estimate clogging time based on SS influent and effluent dataSSBBNonBased on mass balance
Sani et al. [15]Describe the impact of SS on cloggingSSGBKozeny-CarmanWater flow Mechanical dispersion Settling Adsorption
Mostafa and van Geel[75]Conceptual model describing the impact of BF growth on the SHPBCWBChange of porosityWater flow Change of SHP based on reduction of pore space
Brovelli et al.[79]Simulate effect of biomass on SHP in saturated conditionsBCWBChange of porosityBiokinetic model BF attachment and detachment
Rajabzadeh et al.[62]Effect of BF on permeability in saturated conditions, including BF detachment due to shear stressBCGBKozeny-CarmanFluid transport Solute transport Biokinetic model BF detachment Clogging
Hua et al.[92]Parameter efficient model to describe bioclogging based on colloid transport modelsBCBBKozeny-CarmanWater flow Solute transport (advection) Linear growth of BF Clogging model
Giraldi et al.[95]Forecast behaviour and treatment properties due to cloggingSS, BCGBKozeny-CarmanWater flow Solute transport Biokinetic model Transport and filtration of SS
Hua et al. [101]Clogging time based on influent of SS, BF, development on BOD and Plant detritusSS, BC, PBBNonBF, SS based on mass balance Plant clogging model
SS: suspended solids, BC: bioclogging, P: plant detritus; BB: Black Box, GB: Grey Box, WB: White Box; SHP: soil hydraulic parameters.
Table 7. Identification and classification of clogging (adapted from [13]).
Table 7. Identification and classification of clogging (adapted from [13]).
Degree of CloggingDescription
Clogging>80% of the surface is permanently ponded between two loadings
Partly clogged30%–70% of the surface is ponded between two loadings
No clogging<30% of the surface is ponded between two loadings
Table 8. Recommendations on information for describing VF wetland experimental setups.
Table 8. Recommendations on information for describing VF wetland experimental setups.
Geometry of Plot/ColumnFilter MediaLoading RatesQuality Parameters
Diameter or length widthMaterial usedOLRCOD
AreaGeneral grain size distributionHLRSS
Total heightd10, d60SS LRVSS
Height of main layerPorosityLoading intervalNH4-N
Height of other layerHydraulic conductivityDuration of a single dosingNO3-N

Share and Cite

MDPI and ACS Style

Pucher, B.; Langergraber, G. The State of the Art of Clogging in Vertical Flow Wetlands. Water 2019, 11, 2400. https://0-doi-org.brum.beds.ac.uk/10.3390/w11112400

AMA Style

Pucher B, Langergraber G. The State of the Art of Clogging in Vertical Flow Wetlands. Water. 2019; 11(11):2400. https://0-doi-org.brum.beds.ac.uk/10.3390/w11112400

Chicago/Turabian Style

Pucher, Bernhard, and Guenter Langergraber. 2019. "The State of the Art of Clogging in Vertical Flow Wetlands" Water 11, no. 11: 2400. https://0-doi-org.brum.beds.ac.uk/10.3390/w11112400

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