Next Article in Journal
Reverse Osmosis Membrane Zero Liquid Discharge for Agriculture Drainage Water Desalination: Technical, Economic, and Environmental Assessment
Previous Article in Journal
The Steady State Characteristics of Multicomponent Diffusion in Micro- and Mesoporous Media for Adsorbable and Nonadsorbable Species
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Review of the Anti-Pollution Performance of Triple-Layer GM/GCL/AL Composite Liners

School of Water Conservancy and Civil Engineering, Zhengzhou University, Zhengzhou 450002, China
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(10), 922; https://doi.org/10.3390/membranes12100922
Submission received: 11 August 2022 / Revised: 12 September 2022 / Accepted: 16 September 2022 / Published: 23 September 2022
(This article belongs to the Section Membrane Analysis and Characterization)

Abstract

:
Landfill leachates contain several types of pollutants and complex components, which pollute soils and groundwater. To compensate for the limitations of single-layer and double-layer liners, a triple-layer liner system composed of a geomembrane (GM), geosynthetic clay liner (GCL), and attenuation layer (AL) was invented and widely used in landfill anti-pollution systems. In this paper, the available literature on triple-layer GM/GCL/AL composite liners is summarized. First, the four main transport processes of pollutants through the composite liner, including convection, diffusion, adsorption, and degradation, were analyzed, and the anti-pollution performances were evaluated. According to this, the pollutant transport model considering the transport activity and transport state was classified, and the solution methods were summarized. Finally, the breakthrough time expressions of the composite liners were determined, which provided a base for evaluating their long-term performance and predicting the service life. The purpose of this literature review is to scientifically evaluate the anti-pollution performance of GM/GCL/AL and provide a scientific base and theoretical guidance for extending its application.

1. Introduction

Soil and groundwater pollution is a global problem, with more than five million contaminated sites around the world [1], and the municipal solid waste (MSW) landfills constitute a potential major threat. A wide range of toxic inorganic and volatile organic compounds (VOCs) are contained in MSW leachate, which pose an environmental threat to groundwater sources. Accordingly, the liner system is one of the most important elements of the modern engineered landfill.
The liner system is made of impermeable materials, such as geomembrane (GM), geosynthetic clay liner (GCL), compacted clay liner (CCL), and attenuation layer (AL). Among them, GM is one of the most widely used geosynthetic materials because of its excellent anti-seepage performance. GM can effectively inhibit the leakage of inorganic pollutants, but is not efficient in blocking the transport of organic pollutants [2]. In order to overcome the defects, some of the liners are combined with GM to form composite liners. In the past, CCL was often used as the GM underlying liner, but due to the high cost and inconvenient construction of CCL, GCL has more recently emerged as a substitute. However, some studies have shown that GM/GCL has a mediocre barrier effect to certain pollutants, but this can be improved by attaching AL, so GM/GCL/AL is proving to be a viable liner system (Figure 1), and the study of its anti-pollution performance is of great theoretical and practical significance to further applications.
This paper reviews the relevant studies and analyzes the research trends and progress in recent decades. As can be seen from Figure 2a–d, the number of studies and researchers showed an overall upward trend, and with a peak in the past 5 years. During this period, scholars have proposed many examined approaches. Moreover, the research of China, Canada, and the United States accounts for the majority of the current research, which is due to the large number of landfills in these countries, and the related policy norms are also more complete. In terms of anti-pollution liner types, most of the research has focused on GM single-layer and double-layer liners, while research on triple-layer composite liners, especially GM/GCL/AL, is relatively lacking.
Therefore, the anti-pollution performance of GM/GCL/AL composite liners was comprehensively evaluated in this paper, and the research process was presented in the form of a flow chart (Figure 3). In general, the anti-pollution performance of liners was evaluated in three aspects, which are pollutants’ transport mechanism, transport model, and breakthrough time. The research methods and important parameters are also summarized in Table 1. Through the evaluation of the above indicators, some suggestions and perspectives were put forward, which provides a basis for the application of the liners.

2. Convection of Pollutants through GM/GCL/AL

The convection of pollutants through the composite liner is mainly through leakage through GM defects, which can be assessed by the leakage rate. Although a small amount of pollutants are transported through GCL/Al, the two layers are still effective in inhibiting the convection of pollutants through the GM.

2.1. GM Defects

Because pollutants leak mainly through the GM, its integrity is the decisive factor in controlling the leakage. However, in an actual landfill, GM defects cannot be avoided because of transportation or construction inadequacies and wrinkles formed during GM installation. Rowe [2] proposed that the leakage rate of pollutants in the defective GM is a sum of: (1) the leakage rate through the holes on the flat GM and (2) the leakage rate through the holes connected to the wrinkles on the GM.
The formula to calculate the leakage rate has been derived in several studies [3,4,5,21,22], most of which were based on experience and experimental data, whereas few were derived theoretically. Giroud [3] classified the formulae of the leakage rate calculation according to the interface condition, log(Ki/Ks), where Ki and Ks represent the permeability coefficients of the interface and underlying soil cushion or GCL, respectively. Table 2 lists the applicable conditions and calculation methods for the various leakage rate formula.
The method available in this paper allows the calculation of leakage rates for conditions beyond those for which experimental verification exists. The reader is therefore urged to use the equations carefully and to apply judgement in interpreting calculation results. At the same time, researchers are encouraged to conduct experiments to quantify leakage rates through composite liners and thereby contribute to the body of available data.
However, it is worth noting that the above formulae can only be applied to a flat GM. When the GM is exposed to solar radiation during installation, it bends upward to form wrinkles [23,24,25,26,27,28]. Because of the large number of wrinkles on GM (the coverage rate of wrinkles can reach 30% or even higher) and most of them being hydraulically connected, the defects most likely combine with them [29,30,31,32,33,34]. To calculate the leakage rate of the hole combined with the wrinkle (Figure 4), the length of the wrinkle is assumed to be significantly greater than its width, and that the leachate first moves laterally in the transfer layer, beyond the boundary of the wrinkle to a wetting radius, after which it moves longitudinally in the underlying soil layer. Rowe et al. [21] provided a formula to calculate leakage rate as:
Q c = 2 L k b + ( k D θ ) 0.5 h d / D
where L (m) is the length of the wrinkle; 2b (m) is the width of the wrinkle; k (m2/s) is the permeability of the clay liner; D (m) is the thickness of the clay liner; θ (m2/s) is the interfacial hydraulic conductivity between the GM layer and clay liner; hd (m) is the water head loss in the composite liner.
Rowe et al. [35] calculated the leakage rate of the leachate in a composite liner made of GM and clay, and showed that the effect of wrinkles is significant, and the length of hydraulically connected wrinkles significantly affects the leakage rate. Zhang [36] further studied the influence of the length of hydraulically connected wrinkles on the leakage rate of composite liners, and showed that the leakage rate of the hole with the wrinkle was 10–5000 times more than that of the hole on the flat GM. Therefore, wrinkles are important because of the increased potential for contaminant migration through a hole in the GM at or near the wrinkle. There is also increased potential for the development of future holes due to stress cracking at points of high tensile stress in the wrinkle. While it is not essential to eliminate all of the wrinkles, the presence of wrinkles should be considered in calculating leakage through composite liners in the design. The construction procedure should be adopted to limit the length of connected wrinkles present at the time the GM is covered.

2.2. Interface Transmissivity

Although the leakage rate is mainly controlled by GM defects, the interface transmissivity between the GM and the underlying medium also influences the leakage rate of composite liners. Chai et al. [37] described two methods for calculating the value of transmissivity: one is by inverse calculation from the results of flow test experiments and the other is by directly measuring the GM/underlying layer interface transmissivity. For example, Barroso [38] measured the leakage rate of a GM (with holes)/GCL/soil-based composite liner through an indoor leakage test. The interface transmissivity was obtained using the formula proposed by Rowe and Jabin [39] for the leakage rate of leachate in the wetted radius of the composite liner. The obtained value of the GM/GCL was 3.2 × 10−11 to 5 × 10−11 m2/s.
Touze-Foltz et al. [40] directly measured the interface transmissivity of a GM with holes and a GCL using an indoor test. The interface transmissivity formula for the GM/GCL was deduced as follows:
θ = a l n R 2 / r 0 l n h 0 / h 1 2 π Δ t
where r0 (m) is the radius of the circular hole; R2 (m) is the outer diameter of the sample; a (m2) is the cross-sectional area of the burette; h0 (m) is the initial water-head difference; h1 (m) is the final water-head difference; and Δt (s) is the time interval. The obtained interface transmissivity was 6 × 10−12 to 2 × 10−10 m2/s.
Chai et al. [41] conducted large-scale and long-term leakage tests on the leachate using a composite liner made of defective GM and clay layers. They showed that when the defect was a circular hole of 5 mm radius, the range of interface transmissivity under different loads was 1.3 × 10−6 to 1 × 10−7 m2/s; when the defect was a 1 × 30 mm rectangular hole, the range became 1.3 × 10−10 to 1.4 × 10−8 m2/s. Considering the influence of load p′ on the leakage rate Q of the GM with holes, the empirical formula of the effect of p′ on θ is given as follows:
θ = θ 0 1 1 + p / p a 2.5
where θ0 (m2/s) is the interface transmissivity when p′ = 0 and pa is the atmospheric pressure.
Additionally, the conditions of contact between the GM and underlying layer play a decisive role in the interface transmissivity. Imperfect contact allows for a transmissive zone between the two liners, and hence, once fluid has migrated through the defect, it can then move laterally in the interface between the two liners before fully percolating through the clay liner. Giroud [3] defined the good and poor contact conditions between the GM and the underlying soil layer; where good contact indicates that the GM is almost wrinkle-free, the underlying soil layer is fully compact, and the surface is flat, while poor contact indicates that the GM has a certain number of wrinkles, and the underlying soil layer is uneven and not fully compact.
Rowe [42] used the empirical formula derived by Giroud [3] to obtain the empirical relationship between the interface transmissivity of the composite liner and permeability coefficient of the underlying GM liner:
① Good contact: log 10 θ = 0.07 + 1.036 log 10 k L + 0.0180 log 10 k L 2
② Poor contact: log 10 θ = 1.15 + 1.092 log 10 k L + 0.0207 log 10 k L 2 where KL (m2/s) is the permeability coefficient of the underlying layer. When KL is 1 × 10−9 m2/s for a typical CCL, the corresponding interface transmissivity (θ) is 1.6 × 10−8 m2/s for a good contact and 1 × 10−7 m2/s for a poor contact. It was observed that the interface transmissivity of a good contact is lower by nearly one order of magnitude than that of a poor contact.

2.3. GCL/AL Consolidation Characteristics

As mentioned above, pollutants mainly leak through the defective GM of a composite liner, and the defect type, liner contact conditions, external load, and properties of the underlying layer would affect the leakage rate. Therefore, for the triple-layer GM/GCL/AL composite liner reviewed in this paper, the characteristics of GM, GCL, and AL need to be considered when calculating the leakage rate, in which the consolidation characteristics are the most important factors.
AL is a natural soil layer. When consolidated, water flows out of the soil, and the solid particles and fluids move vertically through the upper and lower surfaces of the unit. Following the mass conservation principle for the transport of pollutants in soil, Peter and Smith [43] obtained the mass conservation equation of the liquid and solid phases through the unit. The Darcy rate for the consolidated soil can be obtained by the following derivation:
q = k v ρ f g u h z + p z
where kv (m2/s) is the permeability coefficient of the soil; ρf (kg/m3) is the density of the liquid phase; μ is the viscosity coefficient; g (m/s2) is the acceleration due to gravity; uh (kPa) is the total pressure of the liquid phase; and p (kPa) is the excess pore water pressure caused by consolidation. The first term on the right of the equation represents the Darcy velocity caused by the hydraulic gradient, and the second term represents the Darcy velocity caused by consolidation.
Based on this, Peters and Smith [43] proposed consolidation equations for small and large deformations using the Terzaghi effective stress theory. For small deformations, the stress change in the soil can be assumed as constant without considering the effect of velocity caused by the hydraulic gradient:
1 + e a v k v ρ f g p z z + σ t = p t
where e is the porosity ratio; av is compression coefficient; and σ (kPa) is the total vertical pressure.
In case of consolidation by large deformations, porosity varies with time and depth (porosity varies uniformly with depth for small deformations). Consolidation can be expressed by the equation for change in the porosity ratio:
e t = ± 1 + e 0 ρ s ρ f 1 d d e k v 1 + e e z + 1 + e 0 z k v ρ f g a v 1 + e 0 1 + e e z
where e0 is the initial porosity ratio and ρs (kg/m3) is the density of the solid phase.
Many scholars have studied the consolidation characteristics of GCL. Kang and Shackelford [44] tested GCL repeat samples in a flexible-walled instrument and measured the relevant parameters. They showed that the consolidation coefficient of the GCL ranged from 5.2 × 10−10 m2/s to 2.1 × 10−9 m2/s, and the hydraulic conductivity coefficient k was ≤5.0 × 10−9 cm/s, indicating that the consolidation behavior of the GCL is consistent with that in small deformations. The influence of consolidation characteristics of the GCL on the semi-permeable membrane behavior has been further studied [45]. Such behavior reflects the relative barrier effect of liquid solutes, which is usually quantified by membrane efficiency (ω), where a larger ω indicates a smaller liquid flow through the GCL. The experimental results showed that the membrane behavior of the GCL was significantly enhanced by consolidation. Lu et al. [46] evaluated the effect of pore water salinity on the consolidation characteristics of the GCL using one-dimensional consolidation testing equipment, and showed that the compression and expansion indexes decreased with increasing liquid salinity, while the characteristics of the elastic hysteretic volume change of the GCL were observed in the permeable consolidation test results.
Thus, it is inferred that the convective transport of pollutants in the composite liner is mainly influenced by leakage through the GM defects, and can be evaluated by estimating the leakage rate. The interface transmissivity and consolidation characteristics are also important factors affecting the leakage rate. The anti-pollution performance of the liner can be preliminarily evaluated by the leakage rate. It is suggested to consider the above factors and select an appropriate formula to calculate the leakage rate for evaluation.
Besides, the study of the convection characteristics of pollutants through the composite liner is lacking in some aspects: (1) The leakage rate of the GM with different defects was calculated based on a single defect, and it is necessary to consider the combined effect of various defects on the leakage rate. (2) Presently, research on interface transmissivity is scarce Most studies are being conducted on the GM double-layer composite liner, with little on the GM/GCL/AL liner. (3) The influence of the consolidation characteristics of GCL/AL consolidation on the leakage rate of the composite liner has not been systematically studied.

3. Diffusion of Pollutants through GM/GCL/AL

The transport of pollutants, such as volatile organic compounds (VOC), is dominated by diffusion, which can be evaluated by the diffusion coefficient. The diffusion mechanisms of pollutants differ in different media.

3.1. Diffusion of Pollutants through GM

The diffusion of pollutants in the GM occurs through the pore space between polymers [47]. This process was simulated using Fick’s law:
J G = D g C g z
where JG [kg/(m2·s)] is the mass flux of pollutants through the GM and Dg (m2/s) is the diffusion coefficient of pollutants in the GM. The transient diffusion equation can be obtained by applying the law of mass conservation:
C g z , t t = D g 2 C g z , t z 2
The diffusion of organic pollutants through the GM has been studied to reveal that the diffusion coefficients of polar molecules, such as water molecules and inorganic ions in the GM, are very small, while those of the organic pollutants are quite large [47]. Eun et al. [48] and Mcwatters [6] showed that the diffusion rate of organic compounds in the GM is quite high. Mcwatters et al. [7] further showed that the breakthrough time of organic chemical transport through a GM was only 24 h, and the pollutant flux reached a steady state after approximately one week. In contrast, Saheli et al. [8] showed that inorganic pollutants cannot diffuse through the GM. Therefore, the diffusion of organic pollutants through GM is the focus of most studies.
The diffusion coefficients of organic pollutants in the GM mainly depend on the pollutant properties. The affinity between the pollutants and GM is often expressed by the n-octanol/water partition coefficient (log Kow); log Kow refers to the ratio of the concentration of an organic pollutant in the n-octanol and water phases. The range of log Kow is generally between −3 and 7, and compounds with lower log Kow values are more hydrophilic. For organic pollutants, Dg increases with log Kow. Rowe et al. [49] obtained a simple empirical formula for the geometric diffusion coefficient based on previous experimental data:
log D g = 12.3624 + 0.9205   l o g K ow   0.3424 log K ow 2
Through the above formula, the approximate diffusion coefficient of organic pollutants can be calculated, and then the diffusion characteristics can be evaluated. However, in order to measure the diffusion coefficient more accurately, it is still recommended to measure it through indoor tests. The technique that best matches actual situations uses a two-compartment diffusion cell, where contaminant diffusion is monitored from the source to the receptor. When equilibrium is reached, the value of Sgf can be deduced from mass balance considerations. Scholars have tested the diffusion coefficients of some organic pollutants and discussed their influencing factors.
Rowe et al. [49] experimentally studied the influence of surface fluorination on GM performance. Through adsorption and diffusion tests of VOCs on traditional 1.5 mm thick high density polyethylene (HDPE) GM and surface-fluorinated GM (f-HDPE), it was found that the diffusion coefficient and permeability of f-HDPE were 1.5–4.5 times lower than those of the traditional HDPE GM. Rowe et al. [50] studied the effect of surface fluorination on the diffusion characteristics of benzene, toluene, ethylbenzene, xylene, and other organic substances (BTEX). It was observed that at 22 °C and 6 °C, the permeability of the BTEX solution in the f-HDPE GM is 2.4 times and 1.8 times lower than that in the ordinary HDPE GM, respectively.
Therefore, the characteristics of the pollutants are the main factors affecting the diffusion characteristics; additionally, the GM properties and ambient temperature are influential factors, and the influence is significant under specific conditions. It is suggested that the above factors should be taken into account as far as possible in practical application.

3.2. Diffusion of Pollutants through GCL

The diffusion of pollutants in the GCL is described by Fick’s law:
J G C L = n e D * C z = D p C z
where JGCL is the flux of pollutants through the GCL and DP = neD* is the diffusion coefficient of porous media. Lake and Rowe [51] found that DP is a linear function of the final porosity ratio eB of the GCL for Na+ and Cl:
D p 1.02 e B 0.89 × 10 10
The diffusion coefficient linearly increases with the increasing pore ratio because of the increase in the free pore space. Moreover, the concentration of the pollutant solution, vertical pressure exerted by GCL hydration, and ambient temperature can all affect the diffusion coefficient of pollutants in the GCL. For example, Lake and Rowe [51] showed that when the concentration of NaCl solution increases from 4.6 g/L to 114 g/L, the diffusion coefficient of Cl through the GCL increased from 1 × 10−10 m2/s to 3 × 10−10 m2/s. Rowe et al. [52] showed that when the temperature increased from 7 to 22 °C, the effective diffusion coefficient of benzene in GCL increased from 2.2 × 10−10 to (3.7–4.0) × 10−10 m2/s.
Previous studies suggest that the diffusion coefficient of pollutants in the GCL is similar to that in soil, so GCL and AL can be considered as a whole in the diffusion effect. In addition, the diffusion coefficient of organic pollutants ((3.2–5.6) × 10−10 m2/s) is greater than that of inorganic pollutants ((1.5–3.0) × 10−10 m2/s).

3.3. Diffusion of Pollutants through AL

The pollutants diffuse by molecular diffusion in AL [53], which is caused by their irregular thermal movement (Brownian motion), and they migrate from the high-concentration area to the low-concentration area.
Steady state diffusion in free solution is usually described by Fick’s first law:
J D = D 0 C z
where JD [kg/(m2·s)] is the diffusion flux and D0 (m2/s) is the diffusion coefficient in the free solution.
The transient molecular diffusion problem is described by Fick’s second law:
C t = z D 0 C z
Compared to diffusion in a free solution, the diffusion of pollutants through AL is more complex and relatively slow. The presence of solids reduces the cross-sectional area of the water flow as the solute diffuses only in the soil pores. However, because of the tortuosity of the surface of soil particles, the seepage path of solutes is more tortuous because of the surface tortuosity of the soil particles, which increases the diffusion distance of the pollutants. To accurately characterize the molecular diffusion mechanism of pollutants in AL, Shackelford [54] proposed a practical definition of the effective diffusion coefficient:
D * = τ D 0
where D* (m2/s) is the effective diffusion coefficient and τ is the tortuosity factor, which reflects the influence of soil particles, skeleton, and pore tortuosity on pollutant diffusion.
The apparent tortuosity factor of soil is usually 0.01–0.07 [55], and varies with soil types. For example, the tortuosity factor of montmorillonite is the smallest (0.15–0.20), followed by those of kaolin and clay that are similar (0.24–0.50 and 0.28–0.46, respectively), silty clay (0.08~0.30), sandy soil (0.28), and sand, which is the largest (0.47–0.51). The tortuosity factor of sand–bentonite has the largest variation range (0.01–0.84) [56,57,58,59,60].
The available methods to measure the effective diffusion coefficient of pollutants in soil include field test data inversion, vertical diffusion test, and horizontal diffusion test, which have been used to measure the effective diffusion coefficients of different pollutants through different types of soil. Xie [61] classified these values and demonstrated that the effective diffusion coefficient of inorganic pollutant ions through soil is in the range of (1–15) × 10−10 m2/s (average = (2–8) × 10−10 m2/s). The effective diffusion coefficient of organic compounds is generally up to 9 × 10−10 m2/s (average = (1.5–6) × 10−10 m2/s). Thus, it was inferred that the diffusion coefficient of organic pollutants in soil is smaller than that of inorganic ions.
In summary, the study of the diffusion of pollutants in a composite liner is mainly focused on organic pollutants, which can be evaluated by the diffusion coefficient. The method to obtain the diffusion coefficient is usually an indoor diffusion test, and the anti-pollution performance of the liner can be preliminarily evaluated by the diffusion coefficient. In addition, the chemical properties of the pollutants and the consolidation characteristics of the soil are also important factors, which should be considered in the actual situation.
However, these studies lack in the following aspects: (1) Studies on the diffusion characteristics of other pollutants are rare, and such studies are necessary to establish a database of the diffusion behavior of different pollutants. (2) Single-layer or double-layer composite liners are mainly studied, and diffusion coefficient of pollutants through the GM/GCL/AL triple-layer liner is rarely measured.

4. Adsorption of Pollutants through GM/GCL/AL

Adsorption is also a common mechanism in pollutant transport, which is mainly evaluated by adsorption coefficient. The adsorption of each layer is described below.

4.1. Adsorption of Pollutants through GM

Generally, contaminants are transported in non-defective GM by three continuous processes: ① distribution of pollutants between GM and leachate; ② diffusion of pollutants in GM; ③ distribution of pollutants between GM and the underlying media [62].
When the leachate and GM are in contact, the organic pollutants in the leachate are adsorbed by the GM, redistributing the organic pollutant concentration in the leachate on the surface between the leachate and the GM. This is essentially a process of adsorption of pollutants in the leachate by GM.
When the GM and leachate are in contact for a sufficient period of time, a definite relationship forms between the equilibrium concentration Cg (g/L) of the leachate in GM and the equilibrium concentration Cf (g/L) of the adjacent fluid. This relationship is usually expressed in terms of the Nernst distribution equation in a linear form (Henry’s law):
C g = S g f C f
where Sgf is the partition coefficient, which is mainly affected by the chemical composition of the fluid, molecular structure of the GM, and the ambient temperature.
Henry’s law can also be used for the linear simulation of the distribution of pollutants between GM and the pore water of the underlying medium, i.e., the desorption of pollutants in the GM:
C g = S g f C f
where Sgf represents the distribution coefficient of pollutants in the GM and the pore fluid of the underlying medium. Currently, there are relatively few reports on the value of Sgf, and Sgf and Sgf are generally considered equal [63].
The adsorption of pollutants in the GM is usually evaluated by the partition coefficient using different experimental techniques. Sangam and Rowe [64] used the immersion/adsorption method, and measured the partition coefficients of three chlorinated and four aromatic hydrocarbons through 2.0 mm thick HDPE GM using gas chromatography. Rowe et al. [65] studied the transport of six types of polychlorinated biphenyls homologues in HDPE GM by the immersion/adsorption test, and found that the partition coefficient ranged from 0.7 × 106 to 7.5 × 106. Zhang [36] also obtained the partition coefficient of bisphenol A (1.4) through 1.5 mm thick GM by the batch adsorption test. All of these obtained coefficients are valuable because they serve as a database for the evaluation of the adsorption performance of GM for different pollutants.
The similarity between the pollutants and GM mainly influences adsorption, and is expressed by log Kow, with minor influence from the molecular weight (Mw) of the pollutant. Rowe [63] summarized the empirical formula of the distribution coefficient based on the data from the previous tests (Table 3). It was shown that the partition coefficient of a VOC could be predicted using the available octanol/water partition coefficient (Kow). This provides a convenient means of estimating the partition coefficient of a VOC since the octanol/water partition coefficients of VOCs are readily available.

4.2. Adsorption of Pollutants through GCL

Pollutants are adsorbed through the GCL by the upper geotextile (capping layer), bentonite, and lower geotextile (load-bearing layer) layers [51]:
K d e q = m d b K d b + m d b G T 1 K d G T 1 + m d G T 1 K d G T 1 + m d G T 2 K d G T 2 m d b + m d b G T 1 + m d G T 1 + m d G T 2
where mdb, mdbGT1, mdGT1, and mdGT2 (g) represent the masses of bentonite, geotextile contained in the bentonite soil layer per unit area, geotextile capping layer, and geotextile load-bearing layer, respectively. The corresponding Kdb, KdGT1, and KdGT2 represent the partition coefficients of the bentonite, capping, and load-bearing layers, respectively.
Pinto et al. [9] carried out batch adsorption tests of ammonium and bisphenol A solutions on bentonite and geotextiles, and showed that the adsorption capacity of bentonite for the two pollutants was significantly higher than that of geotextiles. Lake and Rowe [66] experimentally studied the adsorption characteristics of aromatic VOCs on each layer of a GCL. They showed that the adsorption capacity of the bentonite layer was almost negligible, and the overall adsorption capacity depended on the adsorption capacity of pollutants in the upper and lower geotextile layers. Saheli et al. [67] obtained the distribution coefficient of bisphenol A in GCL (10–16 mL/g).
Additionally, Zhang et al. [68] performed a chemical analysis of landfill leachate before and after infiltration and showed that GCL has a strong adsorption capacity for cations and organic matter at the beginning of infiltration, which gradually reaches saturation and loses its adsorption capacity. Therefore, GCL is used as an anti-seepage barrier mainly for its low permeability rather than its adsorption of organic matter and harmful ions.
The above studies showed that adsorption of GCL is affected by several factors, including the properties of pollutants, bentonite, and geotextiles, and the ambient temperature. The relationship between the adsorption parameters of GCL on pollutants and the change in temperature is [69]:
K d T 1 K d T 2 = θ T 1 T 2
where KdT1 and KdT2 represent the distribution coefficients at temperatures, T1 and T2 (°C), respectively, and θ is determined by the following formula:
θ = X R T 1 T 2
where X (K) is the heat of the solution. θ is generally between 1.0 and 1.1. Rowe et al. [63] obtained θ values between 1.03 and 1.04 by fitting experimental data.

4.3. Adsorption of Pollutants through AL

When solid pollutants migrate to the soil pore space, they come in contact with soil particles and are captured by them through physical adsorption, ion exchange, and interface precipitation. The adsorption models can be divided according to various adsorption states into equilibrium adsorption, non-equilibrium adsorption, desorption, and competitive adsorption [10,11,70,71,72,73,74,75,76,77,78]. The classical adsorption models for different modes are summarized in Table 4.
The above adsorption models have been widely used in various studies. Among them, the linear adsorption model is the most convenient and often used in the derivation and solution of analytical solutions. However, the linear adsorption model assumes that there is no upper limit to adsorption, which is not consistent with the actual situation. Similarly, the Freundlich adsorption model also has this defect. The Langmuir model and the two-stage adsorption model are based on the upper limit of soil adsorption capacity, and compared to the Freundlich model, the expression of two-stage adsorption model is relatively easy, which can be used for the derivation of the analytical solution of the breakthrough time of pollutants. In addition, the competitive adsorption models are usually applied to consider the simultaneous adsorption of multiple pollutants. For different pollutants, their adsorption types and characteristics can be determined according to empirical methods, and can also be determined through indoor tests or field tests.
In addition, numerous experiments have been conducted on the adsorption properties of different soil types for different solutes using batch and soil column experiments. Xie et al. [79] used batch experiments to study the influence of ammonium ion concentration, temperature, and pH on the adsorption performance of ammonia nitrogen in loess, and obtained the maximum adsorption capacity of ammonium nitrogen in loess by the Langmuir model fitting (72.3 mg/g). Zhao et al. [80] used batch experiments to study the adsorption kinetics and isothermal adsorption of Bisphenol A on three soil surfaces, and analyzed the effects of different soil–water ratios on Bisphenol A adsorption. They showed that the second-order kinetic equation could describe the adsorption process of Bisphenol A. The adsorption rate of Bisphenol A increased till the soil–water ratio reached 1:10; a further increase in the soil–water ratio did not affect it. Therefore, when evaluating the adsorption performance, it is necessary to select a suitable adsorption model for the description so that the results will be more reasonable.
In this section, the adsorption mechanisms and influencing factors of the GM/GCL/AL liner were discussed and the adsorption model of the AL liner was systematically summarized. The anti-pollution performance of the liner can be preliminarily evaluated by the diffusion coefficient. However, the current study on the adsorption characteristics of pollutants through composite liners lacks in the following aspects: (1) The adsorption parameters are currently determined in single-and double-layer liners, and similar studies are required for the three-layer GM/GCL/AL composite liners. (2) In case of triple-layer composite liners, the studies mostly discuss each layer; therefore, systematic and reasonable adsorption models can be established by comprehensively considering the characteristics of each layer and the interface between the liners.

5. Degradation of Pollutants through GM/GCL/AL

In contrast to the above mechanism, in the GM/GCL/AL liners, organic pollutants are degraded mainly through microbial degradation in the AL. Even though the properties of the GCL and clay liner are similar, the degradation of pollutants in GCL is negligible owing to its thinness.
Biodegradation is a controlled process that determines the persistence and mobility of organic matter in the soil. Studies have shown that organic pollutants, such as aliphatic and aromatic compounds, can be effectively degraded in the soil [81]. Ren et al. [12,13] found that degradation in municipal solid-waste landfills is the main cause of attenuation of VOCs. The degradation of organic pollutants in soil can effectively reduce landfill leachate pollution into groundwater [82,83]. Therefore, the impact of degradation must also be considered when evaluating the anti-pollution performance of landfill liner systems [84].
The decay of organic matter in the landfill leachate can be divided into two stages: (1) decay in the AL and (2) attenuation in the aquifer. The decay of AL is influenced by ion exchange, precipitation, dissolution, and biodegradation. However, in an aquifer, pollutants migrate with groundwater and are attenuated by filtration and precipitation in the aquifer medium. The organic pollutants in porous media are degraded by microorganisms, as they metabolize the useful compounds, while the compounds that cannot be used in microbial metabolism remain in the leachate. The biodegradation of organic matter is mainly divided into aerobic and anaerobic degradation. To date, numerous experimental and theoretical studies on biodegradation have been published [85,86,87,88]. For the experimental study, most of these studies conducted the tank and column experiment to analyze biodegradation, which is the upper (influent) reservoir concentrations decreased as a function of time even though the influent bag concentrations were kept nearly constant. For most of the VOCs, the concentration decreased faster in the influent reservoir as the soil/water partition coefficient increased. Several factors suggested that VOC degradation was significant even though disinfectants were added to the influent.
For the theoretical study, biodegradation was found to follow a first order reaction; biodegradation is usually expressed by the following equation [88]:
C t = λ C
where C (g/L) is the concentration of pollutants and λ is the first-order degradation coefficient, which can be obtained by the following formula:
C C 0 = 1 2 = e λ t   1 / 2
where λ = ln 2 t 1 / 2 = 0.693 t 1 / 2 ; t1/2 is the half-life time of pollutant degradation.
The studies have shown that the anaerobic degradation of organic pollutants occurs in the bottom liner of the landfill, which is the most common degradation mode in soil–water systems. Biological barriers are an effective way of increasing the biodegradation of organic pollutants in barrier systems, thereby reducing the entry of organic pollutants into aquifers.
Li et al. [89] studied the biodegradation of perchloroethylene (PCE) in a biological barrier through batch experiments and analyzed the effects of temperature and pH on it. They showed that the optimum pH for the biodegradation of PCE in the biological barrier was 7, the first-order degradation coefficient of PCE was 0.2489 d−1 at 20 °C, but decreased to 0.1233 d−1 when the temperature reduced to 12 °C. Careghini et al. [90] conducted soil column and batch experiments and obtained first-order degradation coefficients of benzene, toluene, ethylbenzene, o-xylene, p-xylene, and methyl tert-butyl (0.62, 0.88, 1.67, 0.69, 1.23, and 0.34 d−1, respectively). Guillén et al. [91] studied the experimental biodegradation of herbicides’ (mainly composed of picloram) reaction walls and found that the first-order degradation coefficient of herbicides in a biological barrier was 0.924–4.23 d−1.
Moreover, Krupa et al. [81] studied the effectiveness of biological barriers for pollutants in landfill leachate through soil column experiments. They showed that the concentration of organic pollutants in the leachate can be significantly reduced through biodegradation, suggesting that the microbial absorption and degradation of organic pollutants in biological barriers is sustainable and effective. Dong et al. [92] showed that the adsorption and biodegradation in the groundwater of landfills were far greater than hydrolysis and photolysis. Wu et al. [93] studied the influence of the degradation half-life time of pollutants in GM/GCL/AL composite liners on the pollutant flux and concentration at the bottom of the liner (Figure 5), which were observed to increase with the extension and increase in the degradation half-life time.
The above studies indicate that biodegradation must be considered when evaluating the antifouling performance of liner systems. Biodegradation is the main form of pollutant degradation and the degradation coefficient can be calculated by above formula. However, there are only a few studies with respect to the degradation of pollutants through composite liners, and most of them involve biological barriers. Such studies need to be conducted to deepen the understanding of the degradation mechanism under the action of different factors.

6. Analytical Model for Pollutant Transport

Based on the above summary, the pollutant transport characteristics through the liners can be understood; thus, in order to further evaluate the anti-pollution performance of liners, it is necessary to adopt the numerical methods or analytical solutions for analysis.
Numerical methods [94,95,96] have been previously used to study the transport model of pollutants in composite liners. The finite layer method developed by Smith et al. [97] and the finite difference method proposed by Foose et al. [98] are used to analyze organic pollutant transport through composite liners. El-zein [99] adopted the equivalent boundary method to transform the GM into a certain boundary condition considering both convection and diffusion, and simplified the calculation by showing that the equivalent boundary calculation is more reasonable. Leo [100] studied pollutant transport in heterogeneous media in landfills using the boundary element method. Guyonnet et al. [101] used the Laplace transform numerical inversion method to study the hydrodynamic equivalence between multilayer mineral liners. Cooke and Rowe [102] studied the two-dimensional transport of pollutants in a barrier system at the bottom of a landfill, using the finite element method. The studies show that numerical models effectively simulate pollutant transport through composite liners.
The above numerical models are important for the simulation of pollutant transport problems through composite liners. However, the complexity of the numerical methods often far exceeds the complexity of the relevant data available for practical situations. In such cases, a simplified analytical solution can often be an economical and effective alternative to numerical methods in many ways. The analytical solutions can be used to understand the basic mechanisms of pollutant transport, conduct sensitivity analysis as dimensionless variables, analyze the experimental results, and easily simulate the transport behavior on a larger spatiotemporal scale [103,104].
Foose et al. [105] established a one-dimensional diffusion model of pollutants in a double-layer soil and obtained an analytical solution under a semi-infinite space. Huysmans and Dassargues [106] studied the diffusion of pollutants in layered media by combining the equivalent diffusion coefficient with a one-dimensional diffusion analytical solution for pollutants in homogeneous media. Cleall and Li [107] developed a one-dimensional diffusion model of organic pollutant transport in a composite liner under semi-infinite and finite-thickness bottom boundary conditions and obtained the corresponding analytical solution. Chen et al. [108] studied pollutant transport through a composite liner with defective GM under a high waterhead. They established a convection–diffusion model of heavy metals and organic pollutants through GM/CCL and GCL/AL composite liners and obtained the analytical solutions.
Although several studies propose the convection–diffusion analytical model of solutes through layered media, the analytical model of pollutant transport in liners has been confined to single-layer or double-layer models for a long time. In recent years, studies have been conducted on three-layer liner models. Wu et al. [93] proposed an analytical solution for the steady-state diffusion of degradable organic pollutants through a triple-layer composite liner; however, this analytical solution did not consider the leakage effect of pollutants through GM defects. Xie et al. [15,109] compensated for the previous shortcomings and developed an analytical model that comprehensively considered convection, diffusion, and degradation. The transport state of the obtained model in GM and GCL was steady; however, the transport of organic pollutants in all of the layers was transient. Based on this, Feng et al. [16,17] derived an analytical solution of degradable organic pollutant transport through a GM/GCL/AL liner using the separation of variables method, and accurately described the transient pollutant transport of pollutants in the entire composite liner.
Notably, although the one-dimensional transport model is relatively uncomplicated compared to the multi-dimensional one, it can also describe the transport characteristics to a large extent. Therefore, several detailed studies have been conducted on one-dimensional transport models. Such studies mainly focus on the transport state and transport mechanism. Control equations and boundary conditions were established to obtain the transport model, which was then solved using an algorithm to obtain the analytical solution of the one-dimensional transport model.

6.1. Governing Equation of One-Dimensional Transport Model

The form of the governing equation is related to two factors: the transport state (steady or transient) of pollutants in the composite liner and the transport mechanism of pollutants through the GM/GCL/AL liner. The four main transport mechanisms were discussed in the first four sections. The governing equations of the different calculation models from various studies are summarized in this section (Table 5).
As can be seen in Table 5, the governing equation is mainly determined by the transport mechanism through each layer. For instance, Feng et al. [16,17] assumes that the transport of pollutants through the GM layer is by diffusion and convection, and the transport through the GCL and Al layers is by diffusion, convection, and degradation. Then, based on the law of conservation of mass, select the relevant parameters Di, Ki, and λi to establish the governing equation of each layer. Therefore, it is recommended to consider a more consistent transport mechanism according to the actual situation and establish the governing equation.

6.2. Boundary Condition of One-Dimensional Transport Model

The boundary conditions include the top and bottom boundary conditions, and the GM/GCL and GCL/AL contact surface continuity conditions. The expression of the boundary conditions differs according to the various transport states. The presently studied transport states include semi-steady, semi-transient, and fully transient. For the GM/GCL/AL, semi-steady denotes steady state transport in the GM and GCL, semi-transient denotes transient state transport in the AL, and fully transient denotes transient state transport in all of the layers.
① The top boundary is the interface between the leachate and GM. The concentration redistribution of pollutants in leachate at the infiltration interface of GM can be described according to different states as follows:
Semi-steady–semi-transient state:
C gm 0 = C 0 K g
Fully transient state:
C gm 0 , t = C 0 K g
② The bottom boundary of the composite liner is assumed to be a zero-concentration (Dirichlet boundary) or zero-flux (Neumann boundary) boundary:
C al L cl , t = 0   ( Dirichlet   boundary   condition )
C al L cl , t z = 0   ( Neumann   boundary   condition )
where Lcl is the total thickness of the triple-layer liner.
③ Continuity interface conditions include pollutant concentration continuity and continuous flux conditions.
Semi-steady–semi-transient state:
C g m L g m K g = C g c l L g m
D g m C g m b L g m z = n g c l D g c l C g c l L g m z
C g c l L g m + L g c l = C a l L g m + L g c l , t
n g c l D g c l C g c l L g m + L g c l z = n a l D a l C a l L g m + L g c l , t z
where K′g is the distribution coefficient between GM and GCL.
Fully transient state:
C g m L g m b , t K g = C g c l L g m , t
D g m b C g m b L g m b , t z = n g c l D g c l C g c l L g m b , t z   or
v a C gm L gm , t K g D gm C gm L gm , t z = n gl v gcl C gcl L gm , t n gcl D gcl C gcl L gm , t z
C g c l L g m + L g c l , t = C a l L g m + L g c l , t
n g c l D g c l C g c l L g m + L g c l , t z = n a l D a l C a l L g m + L g c l , t z
Additionally, the initial conditions are:
C gm z , 0 C gcl z , 0 C al z , 0 = 0

6.3. Analytical Solution of One-Dimensional Transport Model

The one-dimensional transport model can be solved by combining the governing equations and boundary conditions. However, the model cannot be solved directly because the governing equations are ordinary differential equations, and the boundary conditions have several parameters. Therefore, analytical solution methods have been developed. The currently available solution methods for the one-dimensional transport model of GM/GCL/AL are mainly divided into the following two types:
① First, the governing equations and boundary conditions are generalized to make them applicable to each layer of the liner; then, the problem is decomposed into two sub-problems according to the transport state, and different conditions are substituted into the sub-problems to solve them individually [105,106]. When convection, diffusion, and degradation are comprehensively considered and the transport state is fully transient, it can be expressed as:
Generalized governing equation:
C i z , t t = D i R d , i 2 C i z , t z 2 v i R d , i C i z , t z λ i C i z , t , i = 1 , 2 , 3
Generalized continuous condition:
C i z , t z = L i = C i + 1 z , t z = L i , i = 1 , 2
n i D i C i z , t z z = L i = n i + 1 D i + 1 C i + 1 z , t z z = L i , i = 1 , 2
Decomposition subproblem:
C i z , t = w i z , t + u i z C 0 , i = 1 , 2 , 3
Sub problem 1: w i z , t t = D i R d , i 2 w i z , t z 2 v i R d , i w i z , t z λ i w i z , t , i = 1 , 2 , 3
Sub problem 2: D i R d , i d 2 u i z d z 2 v i R d , i d u i z d z λ i u i z = 0 , i = 1 , 2 , 3 where i = 1, 2, and 3 represent the GM, GCL, and AL, respectively;
② Governing equations of each layer of the liner are solved separately to obtain the general solution with parameters, which is then substituted into the different boundary conditions to obtain the parameters and their corresponding analytical solution [15, 109]. Considering the solving process of Xie et al. [15] as an example:
General solution of GM: C g m z = A 1 + A 2 exp v a z D g m b
General solution of GCL: C g c l z = A 3 + A 4 exp v a z n g c l D g c l
Notably, AL has two different bottom boundary conditions; therefore, the solution may differ depending on them. The Laplace change method is often used to solve the problem for the Neumann boundary condition, while the separation of variables method is often used for the Dirichlet boundary condition, which is the same as the ① solution, where the two sub-functions are solved individually. Considering the solving process of Xie et al. [15] as an example.
For the Neumann boundary condition:
C a l z , p ¯ = 0 = e p t C a l z , t d t
where p is the Laplace variation parameter.
The governing equation of AL under such conditions is:
p C a l z , p ¯ C i = D a l R d 2 C a l z , p ¯ z 2 v a l R d C a l z , p ¯ z λ C a l z , p ¯
For the Dirichlet boundary condition: C a l z , t = θ z , t + ϕ z where θ(z,t) and Φ(z) are the variables satisfying the homogeneous and non-homogeneous differential equation, respectively.
Sub problem 1: θ z , t t = D a l R d 2 θ z , t z 2 v a l R d θ z , t z λ θ z , t
Sub problem 2: D a l R d 2 ϕ z z 2 v a l R d ϕ z z λ ϕ z = 0
Summarizing the previous studies on the analytical solution for the one-dimensional transport model of GM/GCL/AL, it is inferred that the transport model is generally established through the control equations and boundary conditions and then solved to obtain the analytical solution.
In conclusion, the analytical solution is an effective way to evaluate the anti-pollution performance of composite liners. The key point of this method is to determine the governing equations and boundary conditions of the pollutant transport; thus, the readers can select different transport mechanisms’ combinations according to actual conditions to establish the governing equations and select boundary conditions based on the transport state. The analytical solution obtained by solving the sub-problem can calculate the pollutant concentration and flux at the bottom of the liners, and then evaluate the anti-pollution performance of the liners.
However, the above mentioned studies are lacking in the following aspects under different circumstances: (1) Owing to its complexity, the studies of transport models based on the finite element method are currently at the bottleneck stage. Thus, it is necessary to develop a simplified and reliable numerical simulation method suitable for the actual situation. (2) The transport models established by the analytical method are mostly one-dimensional, and can be expanded into two-dimensional and three-dimensional problems of pollutant transport in composite liners. (3) The analytical solution method of the transport model must be studied in detail. (4) To simulate the actual conditions of the landfill site, it might be necessary to study the transport model considering more comprehensive transport effects and states.

7. Breakthrough Time

Based on the transport model mentioned before, the expression of breakthrough time can be further derived. The breakthrough time of pollutants is one of the important indexes to evaluate the anti-pollution performance of liners, which means the transport time required for the pollutant concentration in the leachate (C0) to reach a certain level (Cb) at the bottom of the liner system, where it becomes the breakthrough concentration. Most countries refer to the contaminant threshold concentration in water quality standards as the breakthrough concentration value.
To study the different indicator pollutants based on different breakthrough concentration standards, Foose [110] used cadmium and dichlorotoluene for the steady-state calculation of transport–dispersion, and the breakthrough time was judged by the maximum contaminant level in the American drinking water standard. Xie et al. [111] used cadmium and benzene to calculate the analytical solution of one-dimensional transport–dispersion, which was evaluated by the pollutant concentration in the national primary drinking water regulations. Zhan et al. [112] used a chloride ion centrifuge to simulate the transport–dispersion of breakthrough pollutants under a high waterhead and selected 10% of the initial concentration to determine the breakthrough time. Zhu et al. [110] established a one-dimensional convection–dispersion finite element model using the organic pollutant Chemical Oxygen Demand (COD) as an indicator, and calculated the breakthrough time based on the concentration standard of COD in China’s groundwater and surface water environment.
The above studies suggest that the calculation methods of breakthrough time mainly include the finite element and analytical solution methods. These two methods help to establish the transport model of pollutants in the composite liner, calculate the pollutant concentration at different moments according to the solution of the model, and finally determine the breakthrough time according to the corresponding concentrations of the breakthrough standard. Considering these aspects, mathematical expressions have been developed for the breakthrough time of pollutants in different composite liners under different effects (Table 6).
As observed from the above table, the breakthrough time of GM/GCL/AL was mainly calculated using the pollutant concentration index of the bottom layer (AL), and the calculation was derived from the analytical solution of the transport model. Xie et al. [15] provided a dimensionless expression for parameters related to the relative concentration (CN) and relative flux (FN) under semi-infinite boundary conditions and drew a dimensionless design curve of CN and FN for TR.
The dimensionless parameters in the expression were defined as follows:
K = K g K g
P 1 = v a l z L g c l L g m b D a l
P 2 = η v a l D a l
T R = D a l t R d z L g c l L g m 2
Q = R d λ z L g c l L g m 2 D a l
where K is the ratio of the distribution coefficient; P1 represents the relative ratio of convection and diffusion in AL; P2 reflects the transport characteristics of pollutants in the GM and GCL; and Q represents the relative ratio of degradation and diffusion in AL. Additionally, the value of time factor TR was determined by two dimensionless Peclet numbers. The specific relationship is shown in Figure 6 and Figure 7.
Furthermore, the factors influencing the breakthrough time have also been discussed. Zhu et al. [18] demonstrated that the initial concentration of pollutants is one of the influencing factors. Different initial concentrations lead to different pollutant concentrations during chemical diffusion. Therefore, a higher initial concentration generally denotes a greater concentration gradient and a shorter breakthrough time. Additionally, adsorption hinders pollutant transport, resulting in a prolonged breakthrough time; therefore, the retardation factors of the liner also affect the degree of extension of the breakthrough time. Chen [114] found that the permeability coefficient and retardation factor are the key factors affecting the breakthrough time, while the hydrodynamic dispersion coefficient and waterhead height have a minor influence on it.
In summary, the breakthrough time is an important indicator for evaluating the long-term performance of a composite liner and predicting its service life. It is determined using indicator pollutants and the corresponding breakthrough concentration standards. Although there are many methods to calculate the breakthrough time, it is simpler and more effective to deduce and calculate the breakthrough time through the analytical solution. Comparing the pollutant concentration at the bottom of the liner with the concentration specified in the standard, it is proved that the GM/GCL/AL composite liners have a good anti-pollution performance.
This section summarized the calculation formulas of breakthrough times under different conditions and discussed its influencing factors. However, the available studies lack in the following aspects: (1) Breakthrough time is mainly calculated from the transport model, which lacks independent and systematic research. (2) The factors influencing the breakthrough time of composite liners are relatively single, and the effects of multiple factors and their combinations on the breakthrough time must be considered. (3) Few studies have been conducted on the breakthrough time of GM/GCL/AL.

8. Summary and Prospects

The research improvements and developments of GM/GCL/AL triple-layer composite liners in recent decades were summarized and the anti-pollution performance of the liners was comprehensively evaluated in this paper. The main conclusions are as follows.
(1)
The convection of pollutants through GM/GCL/AL is mainly leakage through GM defects, which can be evaluated by calculating the leakage rate. It is recommended to determine the type of defect according to the actual situation on site and calculate the leakage rate according to the formula. Additionally, factors such as the interface transmissivity coefficient and properties of the underlying medium also affect convection, which are also necessary to be considered. However, the studies on pollutant convection through composite liners have mainly focused on single defects rather than multiple or combination of defects, and lack systematic research on the GM/GCL/AL characteristics;
(2)
The diffusion of pollutants through GM/GCL/AL mainly consider the transport of organic pollutants. The diffusion coefficient is the parameter to evaluate the diffusion characteristics, and it is recommended to measure it by an indoor diffusion test. In addition, the chemical properties of the pollutants and the consolidation characteristics of the soil are also important factors, which should be considered in the actual situation. In future research, the diffusion behavior of other pollutants in the composite liner could be studied to establish the diffusion coefficient database of different pollutants through the composite liner;
(3)
The adsorption of pollutants through the GM/GCL/AL was mainly evaluated by the adsorption coefficient. The adsorption coefficient is recommended to be obtained through an indoor adsorption test. In addition, the adsorption model is also the basis to evaluate the adsorption characteristics. It is recommended to select a suitable adsorption model according to the actual situation; the Langmuir model and the two-stage adsorption model are the more common and reasonable models in the current study. However, current research on triple-layer liners is limited, and the interaction and interface characteristics between liners are not comprehensively considered;
(4)
The degradation of organic pollutants in soil can effectively reduce landfill leachate pollution into groundwater, so it is suggested to consider the impact of degradation when evaluating the anti-pollution performance of landfill liner systems. Biodegradation is the main form of pollutant degradation, and the degradation coefficient can be calculated by the summarized formula. However, it is necessary to strengthen the study of the degradation characteristics of pollutants in the triple-layer GM/GCL/AL composite liners and deepen the research on degradation mechanisms under different conditions;
(5)
The research methods for pollutant transport models mainly include numerical methods and analytical method. An analytical solution is more recommended because of its simplicity and effectiveness, and this study mainly discussed the analytical solution of a one-dimensional pollutant transport model. Determining the governing equations and boundary conditions are the keys to solving the problem. It is suggested to establish appropriate governing equations according to the transport action, select boundary conditions based on the transport state, and then obtain analytical solutions by the subproblem method. However, in the current studies, the available transport models lack simplified and practical numerical methods, with few studies on multi-dimensional transport models. Furthermore, the simplicity of the transport mechanisms and states could be considered;
(6)
The breakthrough time of pollutants is one of the important indexes to evaluate the anti-pollution performance of liners. There are many methods to calculate the breakthrough time, and in this paper, it is recommended to use the analytical solution of the transport model for derivation and calculation. By comparing the pollutant concentration at the bottom of the liner with the concentration specified in the standard, it is proved that the GM/GCL/AL composite liners have good anti-pollution performance. However, there is still a lack of separate and systematic studies on breakthrough times, such as the use of a single evaluation index of indicator pollutants. Thus, studies need to be conducted on the coupling effect of multiple factors.

Author Contributions

J.L.: resources, and project administration; C.H.: writing—original draft; J.Z.: funding acquisition and supervision; Z.Z.: formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (51909243), National Natural Science Foundation of China (No. 52279113) and the Key Scientific and Technological Project of Henan Province (222102320097).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wuana, R.A.; Okieimen, F.E. Heavy Metals in Contaminated Soils: A Review of Sources, Chemistry, Risks and Best Available Strategies for Remediation. ISRN Ecol. 2011, 2011, 402647. [Google Scholar] [CrossRef]
  2. Rowe, R.K. Short-and long-term leakage through composite liners. The 7th Arthur Casagrande Lecture. Can. Geotech. J. 2012, 49, 141–169. [Google Scholar] [CrossRef]
  3. Giroud, J.-P. Leakage Control using Geomembrane Liners. Soils Rocks 2016, 39, 213–235. [Google Scholar] [CrossRef]
  4. Badu-Tweneboah, K.; Giroud, J.P. Discussion of “Leakage through Holes in Geomembranes below Saturated Tailings” by R. Kerry Rowe, Prabeen Joshi, R.W.I. Brachman, and H. McLeod. J. Geotech. Geoenviron. Eng. 2018, 144, 07018001. [Google Scholar] [CrossRef]
  5. Cunningham, H. Evaluating the Rate of Leakage through Defects in a Geomembrane. Ph.D. Thesis, University of Saskatchewan, Saskatoon, SK, Canada, 2019. [Google Scholar]
  6. Mcwatters, R. Diffusive Transport of Volatile Organic Compounds through Geomembranes. Ph.D. Thesis, Queen’s University, Kingston, ON, Canada, 2015. [Google Scholar]
  7. McWatters, R.S.; Rutter, A.; Rowe, R.K. Geomembrane applications for controlling diffusive migration of petroleum hydrocarbons in cold region environments. J. Environ. Manag. 2016, 181, 80–94. [Google Scholar] [CrossRef]
  8. Saheli, P.T.; Rowe, R.K.; Petersen, E.J.; O’Carroll, D.M. Diffusion of multiwall carbon nanotubes through a high-density polyethylene geomembrane. Geosynth. Int. 2017, 24, 184–197. [Google Scholar] [CrossRef]
  9. Pinto, S.M.D.O.; Bila, D.M.; Ritter, E. Ammonium and BPA sorption for GCL: Towards a sustainable geoenvironment. In Proceedings of the 8th International Congress on Environmental Geotechnics, Hangzhou, China, 28 October–1 November 2018; pp. 462–468. [Google Scholar] [CrossRef]
  10. Wang, Z.; Li, Y.; Guo, P.; Meng, W. Analyzing the Adaption of Different Adsorption Models for Describing the Shale Gas Adsorption Law. Chem. Eng. Technol. 2016, 39, 1921–1932. [Google Scholar] [CrossRef]
  11. Seaman, J.C. Competitive Sorption and Transport of Heavy Metals in Soils and Geological Media. Soil Sci. Soc. Am. J. 2013, 77, 2216. [Google Scholar] [CrossRef]
  12. Ren, X.; Zeng, G.; Tang, L.; Wang, J.; Wan, J.; Wang, J.; Deng, Y.; Liu, Y.; Peng, B. The potential impact on the biodegradation of organic pollutants from composting technology for soil remediation. Waste Manag. 2017, 72, 138–149. [Google Scholar] [CrossRef]
  13. Ren, X.; Zeng, G.; Tang, L.; Wang, J.; Wan, J.; Liu, Y.; Yu, J.; Yi, H.; Ye, S.; Deng, R. Sorption, transport and biodegradation–An insight into bioavailability of persistent organic pollutants in soil. Sci. Total Environ. 2017, 610–611, 1154–1163. [Google Scholar] [CrossRef]
  14. Singh, P.; Borthakur, A. A review on biodegradation and photocatalytic degradation of organic pollutants: A bibliometric and comparative analysis. J. Clean. Prod. 2018, 196, 1669–1680. [Google Scholar] [CrossRef]
  15. Xie, H.; Zhang, C.; Feng, S.; Wang, Q.; Yan, H. Analytical Model for Degradable Organic Contaminant Transport through a GMB/GCL/AL System. J. Environ. Eng. 2018, 144, 04018006. [Google Scholar] [CrossRef]
  16. Feng, S.-J.; Peng, M.-Q.; Chen, Z.-L.; Chen, H.-X. Transient analytical solution for one-dimensional transport of organic contaminants through GM/GCL/SL composite liner. Sci. Total Environ. 2018, 650, 479–492. [Google Scholar] [CrossRef] [PubMed]
  17. Feng, S.-J.; Peng, M.-Q.; Chen, H.-X.; Chen, Z.-L. Fully transient analytical solution for degradable organic contaminant transport through GMB/GCL/AL composite liners. Geotext. Geomembr. 2019, 47, 282–294. [Google Scholar] [CrossRef]
  18. Zhu, W.; Shu, S.; Wang, S.; Xu, H.; Fan, X. Investigation of indicating pollutant for landfill leachate breaking through anti-seepage system. Chin. J. Geotech. Eng. 2016, 38, 619–626. [Google Scholar] [CrossRef]
  19. Chen, Y.; Xie, H.; Ke, H.; Chen, R. An analytical solution for one-dimensional contaminant diffusion through multi-layered system and its applications. Environ. Earth Sci. 2008, 58, 1083–1094. [Google Scholar] [CrossRef]
  20. Shu, S.; Zhu, W.; Shi, J. A new simplified method to calculate breakthrough time of municipal solid waste landfill liners. J. Clean. Prod. 2019, 219, 649–654. [Google Scholar] [CrossRef]
  21. Rowe, R.K.; Fan, J. A general solution for leakage through geomembrane defects overlain by saturated tailings and underlain by highly permeable subgrade. Geotext. Geomembr. 2022, 50, 694–707. [Google Scholar] [CrossRef]
  22. Rowe, R.K.; Joshi, P.; Brachman, R.W.I.; McLeod, H. Leakage through Holes in Geomembranes below Saturated Tailings. J. Geotech. Geoenviron. Eng. 2016, 143, 04016099. [Google Scholar] [CrossRef]
  23. Rowe, R.; Yu, Y.; Hosney, M. Interaction between clogging of a leachate drainage layer and leakage through a composite liner. In Proceedings of the Coupled Phenomena in Environmental Geotechnics (CPEG 2013), Torino, Italy, 1–3 July 2013; pp. 283–290. [Google Scholar] [CrossRef]
  24. Chen, C.; Zhan, L.-T.; Wang, Y.; Zhang, S.; Chen, Y.-M. Stochastic evaluation of leakages through holes in wrinkle networks of composite liners. Geotext. Geomembr. 2020, 48, 284–296. [Google Scholar] [CrossRef]
  25. Lake, C.B.; Rowe, R.K. A comparative assessment of volatile organic compound (VOC) sorption to various types of potential GCL bentonites. Geotext. Geomembr. 2005, 23, 323–347. [Google Scholar] [CrossRef]
  26. Xie, H.; Chen, Y.; Ke, H.; Tang, X.; Chen, R. Analysis of diffusion-adsorption equivalency of landfill liner systems for organic contaminants. J. Environ. Sci. 2009, 21, 552–560. [Google Scholar] [CrossRef]
  27. Kandris, K.; Pantazidou, M. Landfill Base Liners: Assessment of Material Equivalency and Impact to Groundwater. Geotech. Geol. Eng. 2011, 30, 27–44. [Google Scholar] [CrossRef]
  28. Pu, H.; Qiu, J.; Zhang, R.; Zheng, J. Analytical solutions for organic contaminant diffusion in triple-layer composite liner system considering the effect of degradation. Acta Geotech. 2020, 15, 907–921. [Google Scholar] [CrossRef]
  29. Zhan, T.L.T.; Chen, Y.-M.; Tang, X.-W. A Preliminary Understanding on Performance of Barriers for MSW Landfills in Southern China. In Advances in Environmental Geotechnics; Springer: Berlin/Heidelberg, Germany; pp. 380–394. [CrossRef]
  30. Chen, R.; Ge, Y.; Chen, Z.; Liu, J.; Zhao, Y.; Li, Z. Analytical solution for one-dimensional contaminant diffusion through unsaturated soils beneath geomembrane. J. Hydrol. 2018, 568, 260–274. [Google Scholar] [CrossRef]
  31. Peng, C.-H.; Feng, S.-J.; Zheng, Q.-T.; Ding, X.-H.; Chen, Z.-L.; Chen, H.-X. A two-dimensional analytical solution for organic contaminant diffusion through a composite geomembrane cut-off wall and an aquifer. Comput. Geotech. 2020, 119, 103361. [Google Scholar] [CrossRef]
  32. Yan, H.; Xie, H.; Ding, H.; Yang, D.; Chen, C.; Chen, Y.; Chen, Y. Analytical solution for one-dimensional steady-state contaminant transport through a geomembrane layer (GMBL)/compacted clay layer (CCL)/attenuation layer (AL) composite liner considering consolidation. Int. J. Numer. Anal. Methods Geéomeéch. 2022, 46, 1046–1063. [Google Scholar] [CrossRef]
  33. Take, W.; Brachman, R.; Rowe, R. Observations of bentonite erosion from solar-driven moisture migration in GCLs covered only by a black geomembrane. Geosynth. Int. 2015, 22, 78–92. [Google Scholar] [CrossRef]
  34. Zhan, L.-T.; Chen, C.; Bouazza, A.; Chen, Y.-M. Evaluating leakages through GMB/GCL composite liners considering random hole distributions in wrinkle networks. Geotext. Geomembr. 2018, 46, 131–145. [Google Scholar] [CrossRef]
  35. Rowe, R. Systems engineering: The design and operation of municipal solid waste landfills to minimize contamination of groundwater A systems engineering approach to minimizing leachate leakage from landfills. Geosynth. Int. 2011, 18, 391–404. [Google Scholar] [CrossRef]
  36. Zhang, C. Mechanisms for Contaminant Transport in Landfill Composite Liners under Thermal Effect and Its Optimization Design Method. Ph.D. Thesis, Zhejiang University, Hangzhou, China, 2018. [Google Scholar]
  37. Chai, J.-C.; Hayashi, S.; Khalili, N. Effect of effective overburden pressure on geomembrane/soil interface transmissivity. Geosynth. Int. 2008, 15, 31–42. [Google Scholar] [CrossRef]
  38. Barroso, M. Fluid Migration through Geomembrane Seams and through the Interface between Geomembrane and Geosynthetic clay Liner. Ph.D. Thesis, Université Joseph-Fourier—Grenoble I, Saint Martin D’hères, France, 2005. [Google Scholar]
  39. Rowe, R.K.; Jabin, F. Effect of prehydration, permeant, and desiccation on GCL/Geomembrane interface transmissivity. Geotext. Geomembr. 2021, 49, 1451–1469. [Google Scholar] [CrossRef]
  40. Touze-Foltz, N.; Bannour, H.; Barral, C.; Stoltz, G. A review of the performance of geosynthetics for environmental protection. Geotext. Geomembr. 2016, 44, 656–672. [Google Scholar] [CrossRef]
  41. Chai, J.-C.; Miura, N.; Hayashi, S. Large-scale tests for leachate flow through composite liner due to geomembrane defects. Geosynth. Int. 2005, 12, 134–144. [Google Scholar] [CrossRef]
  42. Rowe, R. Environmental geotechnics: Looking back, looking forward. Talian Geotech. J. Riv. Ital. Di Geotec. 2018, 4, 8–40. [Google Scholar]
  43. Peters, G.; Smith, D. One-Dimensional Contaminant Transport through a Consolidating Soil: Application to Contaminant Transport through a Geocomposite Liner; CRC Press: Boca Raton, FL, USA, 2020; pp. 481–486. [Google Scholar] [CrossRef]
  44. Kang, J.-B.; Shackelford, C.D. Consolidation of a Geosynthetic Clay Liner under Isotropic States of Stress. J. Geotech. Geoenviron. Eng. 2010, 136, 253–259. [Google Scholar] [CrossRef]
  45. Kang, J.-B.; Shackelford, C.D. Consolidation enhanced membrane behavior of a geosynthetic clay liner. Geotext. Geomembr. 2011, 29, 544–556. [Google Scholar] [CrossRef]
  46. Lu, Y.; Abuel-Naga, H.; Leong, E.-C. Mechanical and osmotic consolidation of geosynthetic clay liners: A laboratory study. Geosynth. Int. 2020, 27, 561–569. [Google Scholar] [CrossRef]
  47. Rowe, R.K. Long-term performance of contaminant barrier systems. Geotechnique 2005, 55, 631–678. [Google Scholar] [CrossRef]
  48. Eun, J.; Tinjum, J.M.; Benson, C.H.; Edil, T.B. Equivalent Transport Parameters for Volatile Organic Compounds in Coextruded Geomem-brane-Containing Ethylene-Vinyl Alcohol. J. Geotech. Geoenviron. Eng. 2018, 144, 04018040. [Google Scholar] [CrossRef]
  49. Rowe, R.K.; Jones, D.D.; Rutter, A. Polychlorinated biphenyl diffusion through HDPE geomembrane. Geosynth. Int. 2016, 23, 408–421. [Google Scholar] [CrossRef]
  50. Rowe, R.K.; Mukunoki, T.; Lindsay, H. Effect of Temperature on BTEX Permeation Through HDPE and Fluorinated HDPE Geomembranes. Soils Found. 2011, 51, 1103–1114. [Google Scholar] [CrossRef]
  51. Lake, C.; Rowe, R. Diffusion of sodium and chloride through geosynthetic clay liners. Geotext. Geomembr. 2000, 18, 103–131. [Google Scholar] [CrossRef]
  52. Rowe, R.K.; Mukunoki, T.; Sangam, H.P. Benzene, Toluene, Ethylbenzene, m & p-Xylene, o-Xylene Diffusion and Sorption for a Geosynthetic Clay Liner at Two Temperatures. J. Geotech. Geoenviron. Eng. 2005, 131, 1211–1221. [Google Scholar] [CrossRef]
  53. Filipović, V.; Ondrašek, G.; Filipović, L. Modelling water dynamics, transport processes and biogeochemical reactions in soil vadose zone. In Groundwater—Contaminant and Resource Management; Javaid, M.S., Ed.; Intech: Rijeka, Croatia, 2016; pp. 133–162. [Google Scholar] [CrossRef]
  54. Shackelford, C.D. The ISSMGE Kerry Rowe Lecture: The role of diffusion in environmental geotechnics. Can. Geotech. J. 2014, 51, 1219–1242. [Google Scholar] [CrossRef]
  55. Carey, G.R.; McBean, E.A.; Feenstra, S. Estimating Tortuosity Coefficients Based on Hydraulic Conductivity. Ground Water 2016, 54, 476–487. [Google Scholar] [CrossRef] [PubMed]
  56. Gillham, R.W.; Robin, M.J.L.; Dytynyshyn, D.J.; Johnston, H.M. Diffusion of nonreactive and reactive solutes through fine-grained barrier materials. Can. Geotech. J. 1984, 21, 541–550. [Google Scholar] [CrossRef]
  57. Wang, Y.; Chen, Y.; Xie, H.; Zhang, C.; Zhan, L. Lead adsorption and transport in loess-amended soil-bentonite cut-off wall. Eng. Geol. 2016, 215, 69–80. [Google Scholar] [CrossRef]
  58. Fu, X.-L.; Wu, H.-L.; Zhang, R.; Jiang, Z.-Y.; Reddy, K.R.; Du, Y.-J. Heavy metals containment by vertical cutoff walls backfilled with novel reactive magnesium-activated slag-bentonite-sand: Membrane and diffusion behavior. J. Clean. Prod. 2021, 328, 129623. [Google Scholar] [CrossRef]
  59. Sample-Lord, K.M.; Ahmed, M.; Malusis, M.A. Diffusion through soil-bentonite backfill from a constructed vertical cutoff wall. Soils Found. 2021, 61, 429–443. [Google Scholar] [CrossRef]
  60. Yang, Y.-L.; Reddy, K.R.; Zhang, T.; Fan, R.-D.; Fu, X.-L.; Du, Y.-J. Enhanced contaminant retardation by novel modified calcium bentonite backfill in slurry trench cutoff walls. Constr. Build. Mater. 2022, 320, 126285. [Google Scholar] [CrossRef]
  61. Xie, H.; Chen, Y.; Thomas, H.R.; Sedighi, M.; Masum, S.A.; Ran, Q. Contaminant transport in the sub-surface soil of an uncontrolled landfill site in China: Site investigation and two-dimensional numerical analysis. Environ. Sci. Pollut. Res. 2015, 23, 2566–2575. [Google Scholar] [CrossRef] [PubMed]
  62. Li, Y.-C.; Cleall, P.J. Analytical solutions for advective-dispersive solute transport in double-layered finite porous media. Int. J. Numer. Anal. Methods Geéomeéch. 2011, 35, 438–460. [Google Scholar] [CrossRef]
  63. Rowe, R.K. Advances and remaining challenges for geosynthetics in geoenvironmental engineering applications. Geotecnia 2007, 30, 3–30. [Google Scholar] [CrossRef]
  64. Sangam, H.P.; Rowe, R. Migration of dilute aqueous organic pollutants through HDPE geomembranes. Geotext. Geomembr. 2001, 19, 329–357. [Google Scholar] [CrossRef]
  65. Rowe, R.K.; Saheli, P.T.; Rutter, A. Partitioning and diffusion of PBDEs through an HDPE geomembrane. Waste Manag. 2016, 55, 191–203. [Google Scholar] [CrossRef]
  66. Lake, C.; Rowe, R. Migration of leachate constituents through a geosynthetic clay liner by diffusion. In Clay Geosynthetic Barriers; CRC Press: Boca Raton, FL, USA, 2021; pp. 177–186. [Google Scholar] [CrossRef]
  67. Saheli, P.T.; Rowe, R.K.; Rutter, A. Diffusion of bisphenol-A (BPA) through an HDPE geomembrane. Geosynth. Int. 2016, 23, 452–462. [Google Scholar] [CrossRef]
  68. Zhang, J.; Xue, G.H.; Liu, D. Adsorption ability of geosynthetic clay liner to cations and organics in landfill leachate. Hydrogeol. Eng. Geol. 2011, 38, 77–81. [Google Scholar] [CrossRef]
  69. Xiao, X.; Chen, B.; Chen, Z.; Zhu, L.; Schnoor, J.L. Insight into Multiple and Multilevel Structures of Biochars and Their Potential Environmental Applications: A Critical Review. Environ. Sci. Technol. 2018, 52, 5027–5047. [Google Scholar] [CrossRef]
  70. Ayawei, N.; Ebelegi, A.N.; Wankasi, D. Modelling and Interpretation of Adsorption Isotherms. J. Chem. 2017, 2017, 1–11. [Google Scholar] [CrossRef]
  71. Wang, J.; Guo, X. Adsorption isotherm models: Classification, physical meaning, application and solving method. Chemosphere 2020, 258, 127279. [Google Scholar] [CrossRef] [PubMed]
  72. Liu, G.; Li, L.; Huang, X.; Zheng, S.; Xu, X.; Liu, Z.; Zhang, Y.; Wang, J.; Lin, H.; Xu, D. Adsorption and removal of organophosphorus pesticides from environmental water and soil samples by using magnetic multi-walled carbon nanotubes @ organic framework ZIF-8. J. Mater. Sci. 2018, 53, 10772–10783. [Google Scholar] [CrossRef]
  73. Yaseen, Z.M. An insight into machine learning models era in simulating soil, water bodies and adsorption heavy metals: Review, challenges and solutions. Chemosphere 2021, 277, 130126. [Google Scholar] [CrossRef]
  74. Xie, H.; Chen, Y.; Lou, Z.; Zhan, L.; Ke, H.; Tang, X.; Jin, A. An analytical solution to contaminant diffusion in semi-infinite clayey soils with piecewise linear adsorption. Chemosphere 2011, 85, 1248–1255. [Google Scholar] [CrossRef]
  75. Borys, P.; Grzywna, Z.J. On the fractality of the Freundlich adsorption isotherm in equilibrium and non-equilibrium cases. Phys. Chem. Chem. Phys. 2016, 18, 20784–20789. [Google Scholar] [CrossRef] [PubMed]
  76. Gao, R.; Liu, D.; Huang, Y.; Li, G. Preparation of diatomite-modified wood ceramics and the adsorption kinetics of tetracycline. Ceram. Int. 2020, 46, 19799–19806. [Google Scholar] [CrossRef]
  77. Riahi, K.; Chaabane, S.; Ben Thayer, B. A kinetic modeling study of phosphate adsorption onto Phoenix dactylifera L. date palm fibers in batch mode. J. Saudi Chem. Soc. 2017, 21, S143–S152. [Google Scholar] [CrossRef]
  78. Smart, S.E.; Mazziotti, D.A. Accelerated Convergence of Contracted Quantum Eigensolvers through a Quasi-Second-Order, Locally Parameterized Optimization. J. Chem. Theory Comput. 2022, 18, 5286–5296. [Google Scholar] [CrossRef]
  79. Xie, H.; Wang, S.; Qiu, Z.; Jiang, J. Adsorption of NH 4 + -N on Chinese loess: Non-equilibrium and equilibrium investigations. J. Environ. Manag. 2017, 202, 46–54. [Google Scholar] [CrossRef]
  80. Zhao, M.; Chen, H.; Li, S.; Xue, A. Sorption behaviors of Bisphenol A on soil samples. Hubei Agric. Sci. 2016, 55, 1678–1681. [Google Scholar] [CrossRef]
  81. Krupa, J.T.; Gandhimathi, R.; Ramesh, S.T.; Nidheesh, P.V.; Singh, T.S.A. Investigation of Biobarrier for Leachate Containment through Batch and Continuous Flow Studies. J. Environ. Eng. 2016, 142, C4015006. [Google Scholar] [CrossRef]
  82. Textor, S.R.; Guillemette, F.; Zito, P.A.; Spencer, R.G.M. An Assessment of Dissolved Organic Carbon Biodegradability and Priming in Blackwater Systems. J. Geophys. Res. Biogeosci. 2018, 123, 2998–3015. [Google Scholar] [CrossRef]
  83. Saeed, T.; Miah, J.; Majed, N.; Alam, K.; Khan, T. Effect of effluent recirculation on nutrients and organics removal performance of hybrid constructed wetlands: Landfill leachate treatment. J. Clean. Prod. 2021, 282, 125427. [Google Scholar] [CrossRef]
  84. Xie, H.; Lou, Z.; Chen, Y.; Jin, A.; Zhan, T.L.; Tang, X. An analytical solution to organic contaminant diffusion through composite liners considering the effect of degradation. Geotext. Geomembr. 2013, 36, 10–18. [Google Scholar] [CrossRef]
  85. Beyer, C.; Chen, C.; Gronewold, J.; Kolditz, O.; Bauer, S. Determination of First-Order Degradation Rate Constants from Monitoring Networks. Ground Water 2007, 45, 774–785. [Google Scholar] [CrossRef]
  86. Singh, D.K. Biodegradation and bioremediation of pesticide in soil: Concept, method and recent developments. Indian J. Microbiol. 2008, 48, 35–40. [Google Scholar] [CrossRef]
  87. Varank, G.; Demir, A.; Top, S.; Sekman, E.; Akkaya, E.; Yetilmezsoy, K.; Bilgili, M.S. Migration behavior of landfill leachate contaminants through alternative composite liners. Sci. Total Environ. 2011, 409, 3183–3196. [Google Scholar] [CrossRef]
  88. Insel, G.; Dagdar, M.; Dogruel, S.; Dizge, N.; Cokgor, E.U.; Keskinler, B. Biodegradation characteristics and size fractionation of landfill leachate for integrated membrane treatment. J. Hazard. Mater. 2013, 260, 825–832. [Google Scholar] [CrossRef]
  89. Li, Y.; Li, J.; Du, Y.; Liu, F. Anaerobic Dechlorination of Perchloroethylene and Factors Influence under Iron Reducing Condition. In Proceedings of the 2009 International Conference on Environmental Science and Information Application Technology, Wuhan, China, 4–5 July 2009; pp. 312–316. [Google Scholar] [CrossRef]
  90. Careghini, A.; Saponaro, S.; Sezenna, E.; Daghio, M.; Franzetti, A.; Gandolfi, I.; Bestetti, G. Lab-scale tests and numerical simulations for in situ treatment of polluted groundwater. J. Hazard. Mater. 2015, 287, 162–170. [Google Scholar] [CrossRef]
  91. Guillén, M.Y.L.O.; Galíndez-Mayer, C.J.; Ordaz, N.R.; Juárez-Ramírez, C.; Santoyo-Tepole, F.; Ramos-Monroy, O. Evaluating the degradation of the herbicides picloram and 2,4-D in a compartmentalized reactive biobarrier with internal liquid recirculation. Environ. Sci. Pollut. Res. 2014, 21, 8765–8773. [Google Scholar] [CrossRef]
  92. Dong, W.H.; Zhang, P.; Lin, X.Y.; Zhang, Y.; Tabouré, A. Natural attenuation of 1,2,4-trichlorobenzene in shallow aquifer at the Luhuagang’s landfill site, Kaifeng, China. Sci. Total Environ. 2015, 505, 216–222. [Google Scholar] [CrossRef] [PubMed]
  93. Wu, X.; Shi, J.; He, J. Analytical solutions for diffusion of organic contaminant through GCL triple-layer composite liner considering degradation in liner. Environ. Earth Sci. 2016, 75, 1371. [Google Scholar] [CrossRef]
  94. Foose, G.J.; Benson, C.H.; Edil, T.B. Comparison of Solute Transport in Three Composite Liners. J. Geotech. Geoenviron. Eng. 2002, 128, 391–403. [Google Scholar] [CrossRef]
  95. Park, M.-G.; Edil, T.B.; Benson, C.H. Modeling Volatile Organic Compound Transport in Composite Liners. J. Geotech. Geoenviron. Eng. 2012, 138, 641–657. [Google Scholar] [CrossRef] [Green Version]
  96. El-Zein, A.; McCarroll, I.; Touze, N. Three-dimensional finite-element analyses of seepage and contaminant transport through composite geosynthetics clay liners with multiple defects. Geotext. Geomembr. 2012, 33, 34–42. [Google Scholar] [CrossRef]
  97. Smith, D.; Rowe, R.K.; Booker, J.R. The analysis of pollutant migration through soil with linear hereditary time-dependent sorption. Int. J. Numer. Anal. Methods Geéomeéch. 1993, 17, 255–274. [Google Scholar] [CrossRef]
  98. Foose, G.; Benson, C.; Edil, T. Analytical Equations for Predicting Concentration and Mass Flux from Composite Liners. Geosynth. Int. 2001, 8, 551–575. [Google Scholar] [CrossRef]
  99. El-Zein, A. A general approach to the modelling of contaminant transport through composite landfill liners with intact or leaking geomembranes. Int. J. Numer. Anal. Methods Geéomeéch. 2007, 32, 265–287. [Google Scholar] [CrossRef]
  100. Leo, C.; Booker, J. A boundary element method for analysis of contaminant transport in fractured and non-fractured porous media. Comput. Geotech. 1998, 23, 165–181. [Google Scholar] [CrossRef]
  101. Guyonnet, D.; Touze, N.; Norotte, V.; Pothier, C.; Didier, G.; Gailhanou, H.; Blanc, P.; Warmont, F. Performance-based indicators for controlling geosynthetic clay liners in landfill applications. Geotext. Geomembr. 2009, 27, 321–331. [Google Scholar] [CrossRef]
  102. Cooke, A.J.; Rowe, R.K. 2D modelling of clogging in landfill leachate collection systems. Can. Geotech. J. 2008, 45, 1393–1409. [Google Scholar] [CrossRef]
  103. Kocabas, I.; Bulbul, M. Modeling solute/contaminant transport in heterogeneous aquifers. Environ. Sci. Pollut. Res. 2014, 22, 3298–3313. [Google Scholar] [CrossRef] [PubMed]
  104. Leij, F.J.; Bradford, S.A.; Sciortino, A. Analytic solutions for colloid transport with time- and depth-dependent retention in porous media. J. Contam. Hydrol. 2016, 195, 40–51. [Google Scholar] [CrossRef] [PubMed]
  105. Foose, G.J. Transit-Time Design for Diffusion through Composite Liners. J. Geotech. Geoenviron. Eng. 2002, 128, 590–601. [Google Scholar] [CrossRef]
  106. Huysmans, M.; Dassargues, A. Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium). Appl. Hydrogeol. 2009, 17, 1901–1911. [Google Scholar] [CrossRef]
  107. Cleall, P.J.; Li, Y.-C. Analytical Solution for Diffusion of VOCs through Composite Landfill Liners. J. Geotech. Geoenviron. Eng. 2011, 137, 850–854. [Google Scholar] [CrossRef]
  108. Chen, Y.; Wang, Y.; Xie, H. Breakthrough time-based design of landfill composite liners. Geotext. Geomembr. 2015, 43, 196–206. [Google Scholar] [CrossRef]
  109. Xie, H.; Thomas, H.R.; Chen, Y.; Sedighi, M.; Zhan, T.L.; Tang, X. Diffusion of organic contaminants in triple-layer composite liners: An analytical modeling approach. Acta Geotech. 2013, 10, 255–262. [Google Scholar] [CrossRef]
  110. Foose, G.J. A steady-state approach for evaluating the impact of solute transport through composite liners on groundwater quality. Waste Manag. 2010, 30, 1577–1586. [Google Scholar] [CrossRef]
  111. Xie, H.; Lou, Z.; Chen, Y.; Jin, A.; Chen, P. An analytical solution to contaminant advection and dispersion through a GCL/AL liner system. Chin. Sci. Bull. 2011, 56, 811–818. [Google Scholar] [CrossRef]
  112. Zhan, L.; Liu, W.; Zeng, X.; Chen, Y. Parametric study on breakthrough time of vertical cutoff wall for MSW landfills and simplified design formula for wall thickness. Chin. J. Geotech. Eng. 2013, 35, 1988–1996. [Google Scholar]
  113. Lin, F.; Ren, H.; Wei, T.; Kang, B.; Li, Y.; Tao, Y. Research on the Contaminant Breakthrough Time Algorithm Based on Thermal Penetration Theory. IEEE Access 2021, 9, 161991–161998. [Google Scholar] [CrossRef]
  114. Chen, Y. A fundamental theory of environmental geotechnics and its application. Chin. J. Geotech. Eng. 2014, 36, 1–46. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of GM/GCL/AL (geomembrane/geosynthetic clay liner/attenuation layer) triple-layer liners.
Figure 1. Schematic diagram of GM/GCL/AL (geomembrane/geosynthetic clay liner/attenuation layer) triple-layer liners.
Membranes 12 00922 g001
Figure 2. Research trend and progress in recent decades: (a) The number of related papers for last 20 years; (b) The number of researchers and examined approaches; (c) The number of related papers in different countries; (d) The percentage of pollution-sealing approaches.
Figure 2. Research trend and progress in recent decades: (a) The number of related papers for last 20 years; (b) The number of researchers and examined approaches; (c) The number of related papers in different countries; (d) The percentage of pollution-sealing approaches.
Membranes 12 00922 g002aMembranes 12 00922 g002b
Figure 3. Flow chart of the research process.
Figure 3. Flow chart of the research process.
Membranes 12 00922 g003
Figure 4. Schematic diagram of pollutant leakage from the hole connected with the wrinkle. Redrawn from Ref. [35].
Figure 4. Schematic diagram of pollutant leakage from the hole connected with the wrinkle. Redrawn from Ref. [35].
Membranes 12 00922 g004
Figure 5. Effect of degradation on pollutant flux (a) and concentration (b) at the bottom of GM/GCL/AL liners. Reprinted with permission from Ref. [93], Copyright 2016 Springer.
Figure 5. Effect of degradation on pollutant flux (a) and concentration (b) at the bottom of GM/GCL/AL liners. Reprinted with permission from Ref. [93], Copyright 2016 Springer.
Membranes 12 00922 g005
Figure 6. Variation curve of CN and FN of pollutant with P1 (P2 = 0.05, Q = 0.1): (a) dimensionless concentration CN; (b) dimensionless flux FN. Reprinted with permission from Ref. [15], Copyright 2018 ASCE Library.
Figure 6. Variation curve of CN and FN of pollutant with P1 (P2 = 0.05, Q = 0.1): (a) dimensionless concentration CN; (b) dimensionless flux FN. Reprinted with permission from Ref. [15], Copyright 2018 ASCE Library.
Membranes 12 00922 g006
Figure 7. Variation curve of CN and FN of pollutant with P2 (P1 = 5, Q = 2): (a) dimensionless concentration CN; (b) dimensionless flux FN. Reprinted with permission from Ref. [15], Copyright 2018 ASCE Library.
Figure 7. Variation curve of CN and FN of pollutant with P2 (P1 = 5, Q = 2): (a) dimensionless concentration CN; (b) dimensionless flux FN. Reprinted with permission from Ref. [15], Copyright 2018 ASCE Library.
Membranes 12 00922 g007
Table 1. Key evaluation indicators and research contents.
Table 1. Key evaluation indicators and research contents.
Evaluation IndicatorMain ContentEvaluation ParametersMain Research MethodsImportant References
Pollutants’ transport mechanismConvectionLeakage rateExperimental research; Empirical formula[3,4,5]
DiffusionDiffusion coefficientExperimental research; Empirical formula[6,7,8]
AdsorptionAdsorption coefficientExperimental research; Adsorption model[9,10,11]
DegradationDegradation coefficientExperimental research; Empirical formula[12,13,14]
Pollutants’ transport modelGoverning equationPollutant concentration and fluxAnalytical solution[15,16,17]
Boundary conditions
Pollutants’ breakthrough timeBreakthrough timeIndicator pollutant concentrationTransport model derivation[18,19,20]
Table 2. Calculation method of leakage rate of different defects.
Table 2. Calculation method of leakage rate of different defects.
Defects TypeApplicable ConditionsCalculation Method
Circular hole [3]Good contact; contact interface conditions: log(Ki/Ks) = 1 Q c = F c K s h t r
F c = 4 + 3.35 r L s
Circular hole [3]Soil composite liner is relatively thick (thickness ≥ 61 cm); contact interface condition: log(Ki/Ks) > 4.5 Q c , G = β c 1 + 0.1 h w L s 0.95 a 0.1 h w 0.9 K s 0.7
Circular hole [4]Soil composite liner is relatively thick (thickness ≥ 61 cm); contact interface condition: log(Ki/Ks) < 4.5Numerical model method is recommended
Circular hole [4]GCL composite liner; contact interface condition: log(Ki/Ks) > 4 Q c , R = π K s r 2 h t L s + 2 h t L s Δ 1 + 2 h t L s Δ 2 2 h w L s Δ 2
Circular hole [4]GCL composite liner; contact interface condition: log(Ki/Ks) < 4 Q a = 2.85 Q R log K i K s 0.73
Rectangular hole [5]Good contact; contact interface conditions: log(Ki/Ks) = 1 Q l = F l K s h t
F l = 1 0.52 0.76 log w L s
Rectangular hole [5]Soil composite liner is relatively thick (thickness ≥ 61 cm); contact interface condition: log(Ki/Ks) > 4.5 Q l = F l K s h t
F l = 2 L s w 2 + L s T K s
Rectangular hole [21]Soil composite liner is relatively thick (thickness ≥ 61 cm); contact interface condition: log(Ki/Ks) < 4.5Numerical model method is recommended
Rectangular hole [21]GCL composite liner; contact interface condition: log(Ki/Ks) > 4 Q a = 2.85 Q R log K i K s 0.73
Rectangular hole [22]GCL composite liner; contact interface condition: log(Ki/Ks) < 4 Q a = 2.85 Q R log K i K s 0.73
Notes: Fc is the dimensionless parameter of the circular hole defect; Ks (m/s) is the saturated permeability coefficient of the soil liner; ht (m) is the water head difference in the composite liner; r (m) is the radius of the defect hole; Ls is the thickness of the underlying liner (CCL, GCL, or Al); hw (m) is the water head above the GM layer; a (m2) is the area of the defective hole; the value of βc depends on the condition of the contact interface, βc = 0.21 denotes good contact and βc = 1.15 denotes poor contact; Ks (m2/s) is the saturated permeability coefficient of the underlying liner; Δ1 and Δ2 are expressions related to the Bessel function, defect radius (r), and hydration radius (R); Qa is the adjusted leakage rate of the GCL composite liner in the one-dimensional case, Qa = Qc, R; Fl is the dimensionless parameter of the rectangular hole defect; and w is the width of the wrinkle.
Table 3. Empirical formula of the partition coefficient of organic pollutants through HDPE GM [63].
Table 3. Empirical formula of the partition coefficient of organic pollutants through HDPE GM [63].
Influencing FactorsEmpirical FormulaFitting Degree R2
Log KowLog Sgf =−1.1523 + 1.2355 Log Kow0.97
MWOxide: Log Sgf =−3.8883 + 10.0363 Log Mw0.81
Chloride: Log Sgf =−2.0467 + 10.0305 Log Mw0.94
Aromatic: Log Sgf =−0.0776 + 10.0322 Log Mw0.95
Aliphatic: Log Sgf =−0.1107 + 10.0442 Log Mw0.91
Table 4. Classical adsorption models of various adsorption modes.
Table 4. Classical adsorption models of various adsorption modes.
Adsorption ModesAdsorption ModelsModel Expressions
Equilibrium adsorptionHenry model [10] S = K d C
Freundlich model [10] S = K F C 1 / n
Langmuir model [11] S = S max K L C 1 + K L C
Two-stage adsorption model [74] S = K C ,   C < C max
S = S max ,   C   C max
Dubinin-Radushkevich (D-R) model [70] ln S = ln S max k ε 2
ε = R T ln 1 + 1 C E = 1 2 k
Non-equilibrium adsorptionLinear irreversible model [71] d S d t = λ 1 C
Linear reversible model [71] d S d t = λ 2 C λ 3 S
Freundlich model [72] d S d t = λ 4 C N λ 5 S
Langmuir model [72] d S d t = λ 6 C S max S λ 7 S
Quasi-first-order kinetic model [73] d S d t = k 1 S e S
Quasi-second-order kinetic model [75] d S d t = k 2 S e S 2
Intraparticle diffusion kinetic model [76] S = k int t   0.5 + S 0
DesorptionFirst-order kinetic model [77] S t = S 0 e k d e l   t
Competitive adsorptionSheindorf–Rebhun–Sheintuch (SRS) model [78] S i = K i C i j = 1 l α i , j C j N i 1
Langmuir model [78] S i = K i C i S max 1 + j K j C j
Notes: ① Equilibrium adsorption: Kd is the partition coefficient, KF and N are Freundlich model constants, S (g) represents the mass of solute adsorbed per unit mass of soil at equilibrium, Smax represents the maximum adsorption capacity, KL is the Langmuir model constant, K is the partition coefficient of static equilibrium adsorption before saturation, C (g/L) is the solution concentration, Cmax (g/L) is the critical solution concentration at the maximum adsorption capacity, k is the parameter of the D-R model, ε is the Polanyi potential, R is the ideal gas constant, and t (°C) is the thermodynamic temperature. ② Non-equilibrium adsorption: λ1 is the reaction rate constant; λ2, λ4, and λ6 are the first-order positive reaction rate constants; λ3, λ5, and λ7 are the first-order negative reaction rate constants; k1 is the quasi-first-order kinetic adsorption rate constant; k2 is the quasi-second-order kinetic adsorption rate constant; kint is the intraparticle diffusion rate constant; S0 is the intercept, which is related to the thickness of the boundary layer. ③ Desorption: St (g) is the residual amount of solute on the adsorbent at time t, S0 (g) is the mass of the adsorbed solute at the beginning of desorption, Kdel is the first-order kinetic adsorption rate constant, and t is the desorption time. ④ Competitive adsorption: subscripts i and j denote the i and j solutes, respectively; l is the total number of solutes, Ki and Kj are the adsorption model constants of a single adsorption system when the i and j solutes exist alone, respectively, Ci and Cj (g/L) are the concentrations of solutes i and j, respectively, and αi,j is the dimensionless competitive adsorption parameter of component i in the presence of component j.
Table 5. Governing equations of the different transport models.
Table 5. Governing equations of the different transport models.
Transport Mechanism and StateSchematic Diagram of Transport ModelGoverning Equations
GM/GCL: diffusion (steady state)
AL: diffusion (transient state)
Reprinted with permission from Ref. [109], Copyright 2013 Springer
Membranes 12 00922 i001 D gm d 2 C gm z d z 2 = 0
D gcl d 2 C gcl z d z 2 = 0
C al z , t t = D al R d , al 2 C al z , t z 2
GM:diffusion (steady state)
GCL/AL: diffusion + degradation (transient state)
Reprinted with permission from Ref. [15], Copyright 2016 Springer
Membranes 12 00922 i002 D gm d 2 C gm z d z 2 = 0
C gcl z , t t = D gcl R d , gcl 2 C gcl z , t z 2 λ gcl C gcl z , t
C a l z , t t = D a l R d , a l 2 C a l z , t z 2 λ C a l z , t
GM/GCL: diffusion + convection (steady state)
AL: diffusion + convection + degradation (transient state)
Reprinted with permission from Ref. [104], Copyright 2018 ASCE Library
Membranes 12 00922 i003 D g m 2 C g m z 2 v a C gm z = 0
D g c l 2 C gcl z 2 v a n g c l C gcl z = 0
C a l z , t t = D a l R d , a l 2 C a l z , t z 2 v a l R d , a l C a l z , t z λ C a l z , t
GM/GCL/AL: diffusion + convection (transient)
Reprinted with permission from Ref. [16], Copyright 2018 Elsevier
Membranes 12 00922 i004 C gm z , t t = D gm 2 C gm z , t z 2 v a K g C gm z , t z
C gcl z , t t = D gcl R d , gcl 2 C gcl z , t z 2 v gcl R d , gcl C gcl z , t z
C a l z , t t = D a l R d , a l 2 C a l z , t z 2 v a l R d , a l C a l z , t z
GM: diffusion + convection (transient)
GCL/Al: diffusion + convection + degradation (transient)
Reprinted with the permission from Ref. [17], Copyright 2019 Elsevier
Membranes 12 00922 i005 C gm z , t t = D gm 2 C gm z , t z 2 v a C gm z , t z
C gcl z , t t = D gcl R d , gcl 2 C gcl z , t z 2 ν gcl R d , gcl C gcl z , t z λ gcl C gcl z , t
C a l z , t t = D a l R d , a l 2 C a l z , t z 2 v a l R d , a l C a l z , t z λ C a l z , t
Notes: Cgm and Cgcl (g/L) are the pollutant concentrations in GM and GCL, respectively; Cgm(z,t), Cgcl(z,t), and Cal(z,t) (g/L) are the pollutant concentrations in GM, GCL, and AL at time t and coordinate z, respectively. ① Diffusion: Dgm, Dgcl, and Dal (m2/s) are the effective diffusion coefficients of the pollutants in GM, GCL, and AL, respectively; Rd,gcl, Rd,al are the blocking factors of AL and GCL, respectively, and their expressions are R d , gcl = 1 + ρ d , gcl K d , gcl n gcl and R d , al = 1 + ρ d , al K d , al n al , where ρd, n, and Kd represent the dry density, porosity, and distribution coefficient of each layer, respectively; Kg is the partition coefficient of the GM and leachate interface. ② Convection: va (m/s) is the Darcy flow velocity of convective leakage through holes, and the expressions are v a = m h Q / A , Q = 2 h w L w l k b + k l θ , k = L gcl + L al L gcl / k gcl + L al / k al , where mh is the number of holes on GM, A (m2) is the area of liner, Q is the leakage rate through the hole connected with a wrinkle, Lw (m) is the length of connecting wrinkle, 2b (m) is the width of connecting wrinkle, hw (m) is the head difference of leachate, l (m) is the sum of the thickness of GCL and AL, K (m2/s) is the harmonic average permeability coefficient of GCL and AL, θ (m2/s) is the hydraulic conductivity between GM and GCL, and vgcl and val (m/s) are the seepage velocities in GCL and AL, respectively. ③ Degradation: λ is the degradation rate constant, and the expressions are λ gcl = ln 2 / t gcl , 1 / 2 , λ al = ln 2 / t al , 1 / 2 , where tgcl,1/2 and tal,1/2 represent the half-life times of organic pollutants in GCL and AL, respectively.
Table 6. Expressions of breakthrough time of pollutants through composite liners.
Table 6. Expressions of breakthrough time of pollutants through composite liners.
Various StudiesExpressions of Breakthrough Time
① Chen et al. [19] deduced the breakthrough time expression of GM/GCL/AL considering the effects of convection and diffusion. t b = T R   R d , a l   z L g c l L g m v a l
where TR is the time factor, Rd,al is the retardation factor of AL, Lgcl (m) is the thickness of GCL, Lgm (m) is the thickness of GM, and val (m/s) is the seepage velocity of the AL layer.
② Xie et al. [111] deduced the breakthrough time expression of GM/GCL/AL considering the effects of convection, diffusion, and degradation. t b = T R R d , a l z L g c l L g m b 2 D a l
Other parameters are the same as in the ② expression; Dal (m2/s) is the effective diffusion coefficient of pollutants in AL.
③ Shu et al. [20] deduced the breakthrough time expression of GM/CCL composite liner with COD as the indicator pollutant index. t b = n R d c l g m + l c 2 h + l g m + l c k s
a C 0 C A b h + l g m + l c k s D h c n c ln C 0 C A + d
where n is the porosity of CCL, Rdc is the retardation factor of CCL, lg and lc (m) are the GM and CCL thickness, respectively; m is defined as DAc/DAg, which is the ratio of the diffusion coefficients of CCL and GM; h (m) is the head of leachate; ks (m2/s) is the hydraulic conductivity of the equivalent composite liner; Dhc (m2/s) is hydrodynamic dispersion coefficient; C0/CA is the relative concentration ratio; a, b, c, and d are the state parameters with values of 0.3781, −0.2968, 0.05135, and 0.3305, respectively.
④ Lin et al. [113] deduced the general expression of anti-pollution barrier breakthrough time according to the thermal penetration theory. t b = x v + 2 u 2 D x v + 2 u 2 D 2 v 2 x 2 v 2
where x (m) is the distance from the measuring point to the pollution source, v (m/s) is the seepage velocity of groundwater (m/d), u is a function of (x, t), and D (m2/s) is the hydrodynamic dispersion coefficient of pollutants in porous media.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, J.; Huang, C.; Zhang, J.; Zhang, Z. Review of the Anti-Pollution Performance of Triple-Layer GM/GCL/AL Composite Liners. Membranes 2022, 12, 922. https://doi.org/10.3390/membranes12100922

AMA Style

Li J, Huang C, Zhang J, Zhang Z. Review of the Anti-Pollution Performance of Triple-Layer GM/GCL/AL Composite Liners. Membranes. 2022; 12(10):922. https://doi.org/10.3390/membranes12100922

Chicago/Turabian Style

Li, Jia, Chuhao Huang, Jingwei Zhang, and Zhanguang Zhang. 2022. "Review of the Anti-Pollution Performance of Triple-Layer GM/GCL/AL Composite Liners" Membranes 12, no. 10: 922. https://doi.org/10.3390/membranes12100922

APA Style

Li, J., Huang, C., Zhang, J., & Zhang, Z. (2022). Review of the Anti-Pollution Performance of Triple-Layer GM/GCL/AL Composite Liners. Membranes, 12(10), 922. https://doi.org/10.3390/membranes12100922

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