Sampling site description and sampling procedure
The North Canal belongs to the Haihe River system. It flows southward from the southern foot of Yanshan Mountain in Changping District of Beijing to the Tongzhou District. At the upstream end of Beiguan gate in Tongzhou District, it is known as Wenyu River, and it flows through Xianghe County of Hebei Province, Wuqing District of Tianjin city, and joins Haihe River at dahongqiao of Tianjin city. It is the main sewage system in Beijing, accounting for 76% of the total sewage discharge in Beijing (Yang et al. 2021).
In April 2019, ten sediment samples were collected along Beisha River (BsR), Wenyu River (WR), Beiyun River (ByR), and Haihe River (HR) (see Fig. 1). Triplicate samples were collected at each site using hand-held PVC corers (150 cm × id 80 mm) to fully represent the effects of land types on sediment quality. Individual sediment sample was packaged in a polyethylene plastic bag, which was previously rinsed by acid, and then air-dried to remove excess water. Finally, all the samples were freeze-dried, ground and sieved using a pre-cleaned 100-mesh sieve before analysis.
Chemical analysis
Samples were firstly digested with hydrofluoric and perchloric acids by microwave in Teflon vessels (Marsx press, CEM), and the obtained solution was stored at 4°C for the analysis of heavy metals. An Optima 2000DV inductively coupled plasma optical emission spectrometer (ICP-OES) (Perkin Elmer, Waltham, MA, USA), and an inductively coupled plasma mass spectrometer (ICP-MS) system (7500a Agilent Technologies, Santa Clara, CA, USA) were employed to measure the content of heavy metals in the solution, with detection limit of 0.003–0.050 mg/L and 0.025–0.200 µg/L, respectively. Certified reference material GBW07401, purchased from the Chinese Environmental Monitoring Center, was used to test the precision and accuracy of metals (i.e. As, Cr, Co, Cu, Ni, Pb, Cd, and Zn). Overall, the recoveries varied in the range of 91–102%, and the relative standard deviations (RSD) in all samples were lower than 3%.
Relevant assessment and statistical analysis
Geoaccumulation index (Igeo), shown as Eq. 1, have been extensively used to determine the metal contamination by indicating the accumulation of metals in sediment (Muller, 1969).
where Cn and Bn are the the calculated concentrations and background concentrations of heavy metals, respectively, and K is the background matrix correction factor (K = 1.5). By taking the non-polluted soil as the baseline, Igeo could be used to quantitatively measure the degree of metal enrichment or contamination and classified into seven levels (Table 1).
Table 1
Criteria of different assessment methods
Geo-accumulation index (Igeo)
|
Igeo
|
Pollution level
|
≤ 0
|
Unpolluted
|
0–1
|
Unpolluted to moderately polluted
|
1–2
|
Moderately polluted
|
2–3
|
Moderately to strongly polluted
|
3–4
|
Strongly polluted
|
4–5
|
Strongly to very strongly polluted
|
> 5
|
Very strongly polluted
|
Enrichment factor (EF)
|
EF
|
Heavy metal source
|
|
|
0.5–1.5
|
Crustal materials or natural weathering processes
|
> 1.5
|
Non-crustal or anthropogenic processes
|
Contamination factor (CF)
|
CF
|
Pollution Contamination
|
|
|
< 1
|
Low contamination
|
|
|
1–3
|
Moderate contamination
|
|
|
3–6
|
High contamination
|
|
|
> 6
|
Very high contamination
|
|
|
Pollution load index (PLI)
|
PLI
|
Pollution level
|
|
|
≤ 1
|
No pollution
|
|
|
> 1
|
Existing pollution
|
|
|
Consensus-based SQGs
|
Qm−PEC
|
Sediment quality
|
|
|
≤ 0.5
|
Not toxic
|
|
|
> 0.5
|
Toxic
|
|
|
Potential ecological risk (PER) index
|
Eir
|
PER of individual elements
|
RI
|
Comprehensive PER
|
≤ 40
|
Low
|
≤ 150
|
Low
|
40–80
|
Moderate
|
150–300
|
Moderate
|
80–160
|
Considerable
|
300–600
|
High
|
160–320
|
High
|
> 600
|
Serious
|
> 320
|
Very high
|
|
|
Enrichment factor (EF) is another extensively employed index to evaluate the potential anthropogenic inputs of observed metals. The EF values of specific heavy metals could be estimated using the Eq. 2 to assess the sources and temporal variation of metals (Bhuiyan et al., 2010):
EF= [Cn (sample) / CAl (sample)] / [Bn (baseline) / BAl (baseline)] (Eq. 2)
where CAl and BAl are the calculated concentration and background concentration of Al in the sediments, respectively. The soil background values were adopted as the baseline values in this study.
The potential ecological risk (PER) index is commonly used to determine the level of heavy metals in sediments, and shown as Eq. 3 (Håkanson 1980):
Eir= Tir×Cif =Tir×(Cis\Cin) (Eq. 3)
where Eir is the PER index of an individual metal, Cif is the single metal pollution factor. Cis and Cin are the level of the heavy metal, and the reference value of the metal, respectively. Tir is the biological toxic factor of an individual metal, where 30 was used for Cd, 2 was used for Cr, 5 was used for Cu, Co, Ni, Pb and 1 was used for Zn in this study (Guo et al., 2010).
Risk index (RI), the sum of Eir, reflects the sensitivity of the biological community to the hazardous chemicals, aiming to illustrate the comprehensive PER caused by the overall contamination.
Single contamination factor (CF) was calculated by Eq. 5:
The obtained data were used to test if there were any anthropogenic impacts. The pollution load index (PLI) was used to determine the quality of the sediments (Suresh et al., 2011):
PLI = (CF1 × CF2 × …... × CFn)1/n (Eq. 6)
CF is the contamination factor and n is the number of elements used.
To assess the effects of multiple heavy metals, the mean probable effect concentration (PEC) quotient (Qm−PEC) was proposed and be calculated using the following equation:
where Cn is the measured concentration of one heavy metal and PECn is the corresponding PEC value. The PEC benchmark values for As, Co, Cd, Cr, Cu, Ni, Pb, and Zn are 30, 4.98, 111, 149, 48.6, 128, and 459 ug/g, respectively.