# Numerical Evaluation and Management Suggestions for Heavy Metal Pollution Risks in a Sludge Landfill: A Case Study from Fuyong Landfill, Shenzhen, China

Authors: Muyun Sun; Kaiyuan He; Shi Shu
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

Despite the extensive attention paid to the transport of heavy metals in sludge landfills, the processes of transporting these pollutants from a landfill to the underground environment are quite complicated and subject to significant uncertainty. In this study, the transport of typical heavy metal pollutants in a sludge landfill through saturated and unsaturated soil zones during rainfall was investigated via numerical modeling. The objectives of the study are to evaluate the heavy metal pollution risk from a sludge landfill under rainfall infiltration conditions and to propose several management suggestions. The results indicate that, during rainfall, heavy metal concentrations at the top of the unsaturated sludge layer decrease rapidly, but they decrease more gradually at the bottom of the layer. The maximum concentration appears in vertical distribution and decreases gradually through the saturated zone. Nickel is the first heavy metal pollutant to break through the low-permeability natural silt barrier. The transport parameters not only influence the simulated time for heavy metal pollutants to break through the silt layer and cause underground environmental pollution but also affect the extent to which the heavy metal pollutants in pore water exceed the guidelines. On the basis of these results, for dredged sludge with heavy metal concentrations significantly exceeding the standard, the concentration of heavy metals in pore water should be reduced before the sludge is landfilled, and a covering layer should be established on the sludge surface to control rainfall infiltration.


---

## Body

## 1. Introduction

Solidification/stabilization technology is widely used to treat sludge dredged from rivers and lakes before they are transported to landfills [1]. However, pollutants in the dredged sludge, such as heavy metals, cannot be removed through solidification/stabilization [2, 3]. These heavy metal pollutants may migrate to underground environments via mechanisms such as convection, diffusion, dispersion, and adsorption, especially in simple landfills common in China [4, 5]. These landfills generally do not have either a bottom liner system or a closure cover system, and their antiseepage properties are dependent on the low permeability of natural clay layers [6, 7]. Under the effect of rainfall infiltration, the dredged sludge in landfills is associated with potential environmental risks. Therefore, the transport of heavy metal pollutants in dredged sludge has already become a research focus in the field of sludge processing and disposal.Previous studies have proposed a number of methods for modeling pollutant transport in underground environments [7–9]. In contrast with pollutant transport in groundwater, pollutant transport in landfills generally occurs in unsaturated porous media. Under rainfall, what the laws of heavy metal concentration in unsaturated porous media are and how the pollutions in landfill migrate into the surrounding area, the simulation study on this phenomenon is still a novel issue. In addition, the pollution level of landfill for the underground environment is closely related to the site address. Therefore, investigations that use representative site conditions will possess more practical significance.In this study, a sludge landfill in the south of China was chosen as the research location, and the transport of typical heavy metal pollutants in the landfill during rainfall was investigated using a saturated-unsaturated pollutant migration numerical model. The objectives of the study are to evaluate the heavy metal pollution risk from landfilled sludge under rainfall infiltration conditions and to propose several management suggestions. This study provides a reference for the prediction and evaluation of the long-term safety of landfill environments as well as landfill design and management.


## 2. Overview about Fuyong Sludge Landfill

Fuyong sludge landfill is located on the east side of Lingdingyang Bay in the Pearl River estuary, west of Shenzhen, China. The landfill has a land area of 230,000 m2 (Figure 1) and the quantity of landfilled sludge after solidification/stabilization is up to 500,000 m3. A 4 m high dam has been built around the landfill. The natural silt layer is used as the impervious barrier, and no closure cover is built. The antiseepage barrier is taken at the wing dam only, as shown in Figure 2.Figure 1
Plan view of Fuyong sludge landfill.Figure 2
Schematic of Fuyong sludge landfill.Sludge samples were collected from the landfill and analyzed (Figure3). The typical heavy metal pollutants and their mean concentrations in the solidified sludge pore water samples from the Fuyong sludge landfill are summarized in Table 1. According to the threshold concentrations for heavy metal pollutants in the Chinese national standards for groundwater quality [10], Pb and Zn concentrations in pore water samples reached the Class IV standard for underground water, and Cd, Cu, and Ni concentrations exceeded the Class V standard. Therefore, the overall pore water of solidified sludge in this landfill was Class V. If the Class III standards for underground water were applied as the pollution threshold for heavy metals, concentrations of the five heavy metals would exceed the standards by between 2.6 and 167 times (Table 1).Figure 3
Sample collection in Fuyong sludge landfill.Table 1
Measured concentrations of heavy metal pollutants in the landfill.
Heavy metal pollutantsCdCuNiPbZnMean concentration (mg/L)0.033.513.340.092.60Class III standard of underground water (mg/L) [10]0.0051.000.020.011.00Excess multiple6.003.51167.009.002.60


## 3. Materials and Methods

### 3.1. Numerical Model

Heavy metal transport was simulated over long periods and under rainfall conditions using a numerical model of solute transport in flowing saturated-unsaturated fluid using the finite element method (FEM). The FEM programDtransu 2D, designed for calculating fluid flow with mass transport via coupled solute advection-dispersion, was employed [9]. This program has been verified in numerous applications pertaining to seepage flow and contaminant migration [4, 5, 11, 12]. The equation governing solute transport and dispersion, taking saturated-unsaturated fluid flow into account, is as follows [9]:(1)Rdθρ∂c∂t=∂∂xiθρDij∂c∂xj−∂∂xjθρVic−θρλRdc−Qc,where Rd (unitless) is the retardation factor, θ (unitless) is the moisture content, ρ (kg/L) is the fluid density, c (kg/L) is the solute concentration, Dij (m2/s) is the hydrodynamic dispersion coefficient, Vi (m/s) is the seepage velocity, λ (unitless) is the attenuation coefficient, and Qc is the source term. The soil-water characteristic curves for soil were described using the Van Genuchten model [13].


### 3.2. Model Generation

In accordance with the stratigraphic distribution in the sludge landfill (Figure2), a saturated-unsaturated calculation model with a height of 17 m was developed, as shown in Figure 4. Some parameters for the model layers were obtained according to the geological investigation report, as shown in Table 2. The longitudinal dispersivities, diffusion coefficients, and retardation factors used to model the transport of the different heavy metal pollutants through the different materials were taken from previous publications and are shown in Table 2 [14–17].Figure 4
Schematic diagram of the calculation model.Table 2
Parameters of the calculation model.
ParametersSolidified sludgeNatural siltSandSilty clayHydraulic conductivity (m/s)1.0 × 10−71.0 × 10−81.850 × 10−51.1 × 10−8Porosity (%)50.060.033.053.0VG–Parametersa (m−1)1.01.05.751.0n (–)1.21.21.61.2θs (–)0.5410.5410.5410.541θr (–)0.20.20.00.2Longitudinal dispersivity (m)0.10.11.00.1Diffusion coefficient (×10−10 m2/s)Cd8.908.908.908.90Cu7.907.907.907.90Ni0.590.590.590.59Pb3.603.603.603.60Zn3.303.303.303.30Retardation factor (–)Cd2.02.01.02.0Cu3.33.31.03.3Ni2.02.01.02.0Pb30.030.01.030.0Zn3.03.01.03.0


### 3.3. Initial and Boundary Conditions

The mean heavy metal concentrations of the solidified sludge pore water samples from the sample survey were set as the initial concentrations (Table1). The initial concentration (c0) of each pollutant at the nodes of the solidified sludge was as follows:(2)c¯xi,0=c0xi,0,0≤xi≤4m.The background concentrations of the heavy metal pollutants were assumed to be negligible, so the initial pollutant concentrations at the nodes of the other layers were 0 mg/L, as shown in the following equation:(3)c¯xi,0=0,4<xi≤17m.The initial groundwater level was set at 1 m below the base of the solidified sludge. The solidified sludge layer and the top meter of the natural silt layer were modelled as unsaturated zones. The model top was set as the rainfall infiltration boundary. Because the annual average rainfall of the landfill area is about 1965 mm, 900 mm was deducted as evaporation. The infiltration rate was set at 1065 mm/y, modelled as a continuous infiltration of 3 mm/d.(4)q¯bxj,t=3mm/d,xi=0m,where qb (m2/s) is the volumetric flow rate at the top boundary. The heads at the bottom boundary were held constant at 12 m, as shown in the following equation:(5)h¯xi,t=12m,xi=17m.


## 3.1. Numerical Model

Heavy metal transport was simulated over long periods and under rainfall conditions using a numerical model of solute transport in flowing saturated-unsaturated fluid using the finite element method (FEM). The FEM programDtransu 2D, designed for calculating fluid flow with mass transport via coupled solute advection-dispersion, was employed [9]. This program has been verified in numerous applications pertaining to seepage flow and contaminant migration [4, 5, 11, 12]. The equation governing solute transport and dispersion, taking saturated-unsaturated fluid flow into account, is as follows [9]:(1)Rdθρ∂c∂t=∂∂xiθρDij∂c∂xj−∂∂xjθρVic−θρλRdc−Qc,where Rd (unitless) is the retardation factor, θ (unitless) is the moisture content, ρ (kg/L) is the fluid density, c (kg/L) is the solute concentration, Dij (m2/s) is the hydrodynamic dispersion coefficient, Vi (m/s) is the seepage velocity, λ (unitless) is the attenuation coefficient, and Qc is the source term. The soil-water characteristic curves for soil were described using the Van Genuchten model [13].


## 3.2. Model Generation

In accordance with the stratigraphic distribution in the sludge landfill (Figure2), a saturated-unsaturated calculation model with a height of 17 m was developed, as shown in Figure 4. Some parameters for the model layers were obtained according to the geological investigation report, as shown in Table 2. The longitudinal dispersivities, diffusion coefficients, and retardation factors used to model the transport of the different heavy metal pollutants through the different materials were taken from previous publications and are shown in Table 2 [14–17].Figure 4
Schematic diagram of the calculation model.Table 2
Parameters of the calculation model.
ParametersSolidified sludgeNatural siltSandSilty clayHydraulic conductivity (m/s)1.0 × 10−71.0 × 10−81.850 × 10−51.1 × 10−8Porosity (%)50.060.033.053.0VG–Parametersa (m−1)1.01.05.751.0n (–)1.21.21.61.2θs (–)0.5410.5410.5410.541θr (–)0.20.20.00.2Longitudinal dispersivity (m)0.10.11.00.1Diffusion coefficient (×10−10 m2/s)Cd8.908.908.908.90Cu7.907.907.907.90Ni0.590.590.590.59Pb3.603.603.603.60Zn3.303.303.303.30Retardation factor (–)Cd2.02.01.02.0Cu3.33.31.03.3Ni2.02.01.02.0Pb30.030.01.030.0Zn3.03.01.03.0


## 3.3. Initial and Boundary Conditions

The mean heavy metal concentrations of the solidified sludge pore water samples from the sample survey were set as the initial concentrations (Table1). The initial concentration (c0) of each pollutant at the nodes of the solidified sludge was as follows:(2)c¯xi,0=c0xi,0,0≤xi≤4m.The background concentrations of the heavy metal pollutants were assumed to be negligible, so the initial pollutant concentrations at the nodes of the other layers were 0 mg/L, as shown in the following equation:(3)c¯xi,0=0,4<xi≤17m.The initial groundwater level was set at 1 m below the base of the solidified sludge. The solidified sludge layer and the top meter of the natural silt layer were modelled as unsaturated zones. The model top was set as the rainfall infiltration boundary. Because the annual average rainfall of the landfill area is about 1965 mm, 900 mm was deducted as evaporation. The infiltration rate was set at 1065 mm/y, modelled as a continuous infiltration of 3 mm/d.(4)q¯bxj,t=3mm/d,xi=0m,where qb (m2/s) is the volumetric flow rate at the top boundary. The heads at the bottom boundary were held constant at 12 m, as shown in the following equation:(5)h¯xi,t=12m,xi=17m.


## 4. Calculation Results

### 4.1. Transport Process of a Representative Pollutant in Saturated and Unsaturated Zones

The long-term unsteady flow was calculated assuming a continuous rainfall intensity of 3 mm/d, and the transport behavior of a representative heavy metal was determined first. Figure5 shows the variation of the Cd concentration distribution in the saturated-unsaturated model over time. Owing to the effect of rainfall infiltration, the Cd concentration on the surface of the unsaturated solidified sludge decreases rapidly, while it decreases more gradually at the bottom of the solidified sludge. The peak concentration appears in vertical distribution, and the pollution scope expands gradually. The Cd concentration declines gradually in the solidified sludge layer, and the decreasing range decreases with depth and time. The Cd concentration in the saturated layers rises gradually, and the increasing range also decreases gradually with depth and time. In places closer to the top of the unsaturated zone, the peak Cd concentration appears earlier and the peak value is higher. The concentrations at various vertical points decrease gradually after reaching peak values. Therefore, the peak concentration decreases from 0.03 mg/L to 0.01 mg/L during the simulation period of 50 years. Additionally, the location of the peak concentration moves downward, from 3 m from the top of the unsaturated zone in the 5th year of the simulation to 8.5 m from the top of the unsaturated zone in the 50th year.Figure 5
Variation of Cd concentration distribution with time.As for the reasons for the above results, when rainwater recharges the top of the model domain, the moisture content at the top of unsaturated solidified sludge increases rapidly while the influence in the saturated zone is relatively weak and the impact time also lags. Then, the flow velocity in the unsaturated solidified sludge pore water increases, and the advection rate rises accordingly, which promotes faster transport of the pollutants. With the gradual infiltration of rain, the pollutant transport depth also increases. The results also demonstrate the rationality of the calculations.


### 4.2. Breakthrough Behaviors of Various Heavy Metal Pollutants

In the simulated landfill, the natural silt layer was modelled as a low-permeability barrier and the sand layer as a permeable layer. Therefore point P, at the junction of the natural silt layer and the sand layer (see Figure4), was chosen as the location to report the concentration changes of the various heavy metal pollutants. Figure 6 shows the breakthrough curves of the heavy metals Cd, Cu, Ni, Pb, and Zn in the sludge landfill. The concentrations of these five heavy metals at point P all show an upward trend. However, the rates of increase vary. The Cd concentration at point P increases slowly with time in the beginning, and an inflection point appears on the breakthrough curve after approximately 10 years of simulated infiltration. Then, the concentration increases linearly and exceeds the pollution threshold concentration after 29.58 years. This indicates that the breakthrough time of Cd is 29.58 years. The Cd concentration then rises more slowly after 45 years and reaches 0.0097 mg/L after 50 years, which is 32.43% of its initial concentration (Figure 6(a)). The Cu concentration reaches 0.68 mg/L (19.37% of the initial concentration) after 50 years but does not exceed the pollution threshold concentration (Figure 6(b)). The breakthrough curve of Ni is similar to that of Cd. However, Ni exceeds the pollution threshold concentration after 16.20 years and reaches 1.14 mg/L (33.99% of the initial concentration) after 50 years (Figure 6(c)). Owing to the large retardation factor, the Pb concentration rises slowly, and it is only 2.55 × 10–8 mg/L in the 50th year (2.84 × 10–5% of the initial concentration). The pollution threshold concentration is not reached by Pb within 100 years (Figure 6(d)). The Zn concentration reaches 0.59 mg/L (22.68% of the initial concentration) after 50 years but also does not exceed its pollution threshold concentration (Figure 6(e)). Zn does not break through the natural silt layer. Therefore, Ni is the first heavy metal to break through the natural silt layer after 16.2 years.Figure 6
Breakthrough curves of various heavy metal pollutants in the sludge landfill. (a) Cd, (b) Cu, (c) Ni, (d) Pb, (e) Zn.
(a)(b)(c)(d)(e)


## 4.1. Transport Process of a Representative Pollutant in Saturated and Unsaturated Zones

The long-term unsteady flow was calculated assuming a continuous rainfall intensity of 3 mm/d, and the transport behavior of a representative heavy metal was determined first. Figure5 shows the variation of the Cd concentration distribution in the saturated-unsaturated model over time. Owing to the effect of rainfall infiltration, the Cd concentration on the surface of the unsaturated solidified sludge decreases rapidly, while it decreases more gradually at the bottom of the solidified sludge. The peak concentration appears in vertical distribution, and the pollution scope expands gradually. The Cd concentration declines gradually in the solidified sludge layer, and the decreasing range decreases with depth and time. The Cd concentration in the saturated layers rises gradually, and the increasing range also decreases gradually with depth and time. In places closer to the top of the unsaturated zone, the peak Cd concentration appears earlier and the peak value is higher. The concentrations at various vertical points decrease gradually after reaching peak values. Therefore, the peak concentration decreases from 0.03 mg/L to 0.01 mg/L during the simulation period of 50 years. Additionally, the location of the peak concentration moves downward, from 3 m from the top of the unsaturated zone in the 5th year of the simulation to 8.5 m from the top of the unsaturated zone in the 50th year.Figure 5
Variation of Cd concentration distribution with time.As for the reasons for the above results, when rainwater recharges the top of the model domain, the moisture content at the top of unsaturated solidified sludge increases rapidly while the influence in the saturated zone is relatively weak and the impact time also lags. Then, the flow velocity in the unsaturated solidified sludge pore water increases, and the advection rate rises accordingly, which promotes faster transport of the pollutants. With the gradual infiltration of rain, the pollutant transport depth also increases. The results also demonstrate the rationality of the calculations.


## 4.2. Breakthrough Behaviors of Various Heavy Metal Pollutants

In the simulated landfill, the natural silt layer was modelled as a low-permeability barrier and the sand layer as a permeable layer. Therefore point P, at the junction of the natural silt layer and the sand layer (see Figure4), was chosen as the location to report the concentration changes of the various heavy metal pollutants. Figure 6 shows the breakthrough curves of the heavy metals Cd, Cu, Ni, Pb, and Zn in the sludge landfill. The concentrations of these five heavy metals at point P all show an upward trend. However, the rates of increase vary. The Cd concentration at point P increases slowly with time in the beginning, and an inflection point appears on the breakthrough curve after approximately 10 years of simulated infiltration. Then, the concentration increases linearly and exceeds the pollution threshold concentration after 29.58 years. This indicates that the breakthrough time of Cd is 29.58 years. The Cd concentration then rises more slowly after 45 years and reaches 0.0097 mg/L after 50 years, which is 32.43% of its initial concentration (Figure 6(a)). The Cu concentration reaches 0.68 mg/L (19.37% of the initial concentration) after 50 years but does not exceed the pollution threshold concentration (Figure 6(b)). The breakthrough curve of Ni is similar to that of Cd. However, Ni exceeds the pollution threshold concentration after 16.20 years and reaches 1.14 mg/L (33.99% of the initial concentration) after 50 years (Figure 6(c)). Owing to the large retardation factor, the Pb concentration rises slowly, and it is only 2.55 × 10–8 mg/L in the 50th year (2.84 × 10–5% of the initial concentration). The pollution threshold concentration is not reached by Pb within 100 years (Figure 6(d)). The Zn concentration reaches 0.59 mg/L (22.68% of the initial concentration) after 50 years but also does not exceed its pollution threshold concentration (Figure 6(e)). Zn does not break through the natural silt layer. Therefore, Ni is the first heavy metal to break through the natural silt layer after 16.2 years.Figure 6
Breakthrough curves of various heavy metal pollutants in the sludge landfill. (a) Cd, (b) Cu, (c) Ni, (d) Pb, (e) Zn.
(a)(b)(c)(d)(e)


## 5. Discussion

The five heavy metal pollutants modelled in this study have different initial concentrations, pollution threshold concentrations, and transport parameters. As a result, the different heavy metal pollutants caused underground environmental pollution at different times under the same site conditions and rainfall conditions. In other words, the breakthrough times for the different heavy metals varied. Among these five heavy metal pollutants, Ni transport was the fastest and its breakthrough time was the shortest. Although the initial concentration of Ni is not the highest, its retardation factor is the lowest, so its absolute (1.14 mg/L) and relative concentrations (33.99%) after 50 years are the greatest reported at point P, and it exceeds the pollutant threshold concentration most significantly (66.8 times). Cu has the highest initial concentration, but it exceeds the pollutant threshold at point P by only 3.51 times. Pb transport is the slowest because its retardation factor is more than 10 times greater than the other heavy metals. The silt layer adsorbs Pb in pore water on a large scale, greatly slowing down the transport of Pb through the layer. Therefore, the transport parameters not only influence the time for heavy metal pollutants to break through the low-permeability natural silt barrier and cause underground environmental pollution. They also affect the extent to which the simulated heavy metal pollutants in pore water exceed the guidelines.For dredged sludge with heavy metal concentrations that greatly exceed the applicable standards, the concentration of heavy metals in the pore water should be reduced before the sludge is landfilled. Where rainfall in the region of a landfill is abundant, a cover layer should be arranged on the sludge surface following the sludge disposal process [18]. Landfill covering should be conducted after the sludge is processed to control rainwater infiltration with an impervious barrier [19]. Finally, closure management should be conducted for landfills and leachate should be collected until the landfill becomes stable [20].The management of landfill involves a combination of environmental, economic, and social factors [21]. The indexing method used in this study is a commonly used method for landfill risk assessment in the world [22]. The time for heavy metal pollutants breakthrough the natural clay liner was calculated in this study, but the complexity of the landfill determines that there are still many uncertainties [23]. Therefore, in the future, factors such as regional characteristics, urban planning, and the conditions of the landfill itself can be further considered to manage the landfill throughout its life cycle [24].


## 6. Conclusions

In this study, the vertical transport of typical heavy metal pollutants through an unsaturated soil zone in a sludge landfill to the saturated soils underneath was investigated using a numerical modeling approach. The results can provide a reference for the prediction and evaluation of the long-term safety of landfill environments as well as landfill design and management. Several conclusions and suggestions could be drawn and are shown below.Under the effect of rainfall, the concentrations of heavy metal pollutants on the surface of an unsaturated sludge layer reduce rapidly while concentrations at the bottom of the unsaturated sludge layer decrease more gradually. The peak concentration appeared in vertical distribution and dropped gradually in the saturated zone. The location of the peak concentrations also declined.Under the effects of rainfall infiltration, Ni was the first pollutant to break through the impervious barrier of natural silt. Therefore, the transport parameters not only influence the time for heavy metal pollutants to break through the low-permeability natural silt barrier and cause underground environmental pollution, but they also affect the extent to which the simulated heavy metal pollutants in pore water exceed the guidelines.



---

*Source: 8854313-2020-09-15.xml*