# Source and Ecological Risk Characteristics of PAHs in Sediments from Qinhuai River and Xuanwu Lake, Nanjing, China

Authors: Zhenhua Zhao; Zhirui Qin; JingJing Cao; Liling Xia
Journal: Journal of Chemistry
 (2017)

---

## Abstract

In order to investigate the residual characteristics, sources, and ecological risk of PAHs in sediment from urban rivers, the sediments of 15 typical sites from Qinhuai River and Xuanwu Lake, which are typical urban rivers and lake, were collected from October 2015 to July 2016; the sources of PAHs in sediment were also identified by several methods. Results showed that ∑PAHs concentration in sediment ranged from 796.2 ng/g to 10,470 ng/g with an average of 2,713.8 ng/g. High molecular weight PAHs with 4-5 rings were most prominent in the sediment during all four seasons. Source characterization studies based on the analysis of diagnostic ratio (triangular plot method), cluster analysis, and positive factor matrix analysis suggested that the PAHs of Qinhuai River Basin were mainly from pyrogenic origin (biomass and coal combustion and vehicular emission), and the petroleum source also cannot be ignored (specially in summer). Most individual PAHs occasionally affect the aquatic organisms. The highest benzo[a]pyrene-equivalent doses (BaPeq dose) appear at the sites of sewage discharge and heavy traffic. So, the PAHs pollution sources of urban water body have obvious seasonal-dependent and human activities-dependent characteristics.


---

## Body

## 1. Introduction

Polycyclic aromatic hydrocarbons (PAHs), a group of typical chemicals containing two or more aromatic rings, are widely distributed in water [1], soil, sediment [2–5], and air [6]. PAHs have been attracting increasing attention, because many are known to be carcinogenic, mutagenic, or teratogenic [7, 8], which are serious health threats to organisms. Hence, United States Environmental Protection Agency (US EPA) has listed sixteen PAHs as priority pollutants. In recent years, PAHs research has focused on the distribution and source apportionment of sediments. PAHs may come from natural sources, but anthropogenic activity is generally considered the major source [9]. Anthropogenic PAHs primarily originate from pyrogenic and petrogenic sources. Pyrogenic PAHs are mainly from incomplete combustion of various fossil fuels (such as coal, oil, and natural gas), residential heating, industrial activity, wood, tobacco, and other hydrocarbons [10]. Petroleum sources also generate PAHs [11], including crude oil and refined petroleum products such as oil leakage [12].PAHs, once produced, enter the aquatic environment via wastewater, surface runoff, atmospheric deposition, and oil leakage [13]. Due to their low water solubility, PAHs are easily adsorbed by particles and colloids when transferred into the water and sediment [14]. After settling into the sediment, PAHs will rarely photochemically and biologically oxidize and, therefore, will accumulate to high concentrations in sediments [15]. The concentration of PAHs in sediments is related to sediment type and size, geographical location, water migration, and other factors. When environmental conditions change, adsorbed PAHs may be rereleased into the water by chemical and biological processes, leading to secondary pollution of the surrounding environment [16, 17]. Sediments act both as a reservoir and a secondary source for PAHs contamination [18].Until recently, researches on surface sediment PAHs mainly concentrated in the estuary area of some rivers and lakes [19–25]. Sofowote et al. (2008) investigated the profiles and relative contributions of pollution sources to the Hamilton harbor using principal component analyses (PCA) and positive matrix factorization (PMF) methods, and results showed both methods identified four factors and gave excellent correlation coefficients between predicted and measured levels of 25 aromatic compounds; both methods predicted similar contributions from coal tar/coal combustion sources to the harbor (19 and 26%, resp.), but PMF afforded better source identification than PCA [26]. Feng et al. (2007) studied the distribution, source, and ecological risk assessment of PAHs in suspended solids and sediment from the Yangtze River in Wuhan [27]. Chen et al. (2007) have researched the Qiantang River in the Yangtze River Delta region, mainly to analyze the temporal and spatial distribution and sources of PAHs in its water and sediments [28]. So far, relatively little research has focused on the seasonal distribution, source analysis, and ecological risk assessment of PAHs in a typical urban river [29]. Qinhuai River, a tributary of the Yangtze River, is located in the southwest of Jiangsu Province. With a total length of 110 km and the basin area of 2630 km2, it is the largest regional river in Nanjing. Internal Qinhuai River is a typical urban river, flowing across the Nanjing city, which directly affected the sustainable economic development of Nanjing city. Xuanwu Lake is the largest shallow lake of Nanjing (3.7 km2), a typical urban shallow lake, which has been subjected to a series of human disturbances, such as the inflow of domestic sewage, surface rain runoff from the surrounding residents and mountain, water diversion from drink water treatment plant (1.8 × 105 t/d), tourist entertainment, and rubbish in lake park. The objective of this article is to explore the applicability of source appointment technologies and the effects of city activities on the distribution and ecological risk of PAHs in surface sediments, provide theoretical support for the pollution sources controlling and desilting treatment, and make most of PAHs removed from the river, to reduce the ecology risk of urban river and lake.


## 2. Material and Methods

### 2.1. Sampling Region and Sites

A map showing the location of the sampling sites is provided in Figure1. Corresponding functional regional characteristics are shown in Table 1. Sediment sampling was carried out along the inner Qinhuai River and Xuanwu Lake in four seasons. The sampling time, based on local seasons, is grouped as follows: autumn (October, 2015), winter (January, 2016), spring (April, 2016), and summer (July, 2016). The middle month of different season was chosen as the typical month of sampling. According to the statistical data estimated by local government, the sedimentation rate is about 10 cm–60 cm a year for different river. Moreover, the sampling tool is the grab sampler, so the sediment samples can be used as representative samples of different seasons. Surface sediment sampling was collected by grab sampler. Each sample consisted of central parts of 5 subsamples randomly collected at 5 points within the study area. The samples were placed into the sampling bottles, sealed and stored at −20°C in the refrigerator for further analysis.Table 1
Sampling sites and the corresponding regional characteristics.
Reach
Code
Position
Sedimentation rate a year
Regional characteristics
Northern section of IQR
1#
Zhenzhu Bridge
30–40 cm
Near the main road and downtown of city, there are floating contaminations.
Eastern section of IQR
2#
Yixian Bridge
Middle of IQR
3#
Wenjin Bridge
50–60 cm
Locates at Chaotian Temple, receiving sewage discharge, with slow water flow due to gate controlling.
4#
Wenzheng Bridge
Receives incoming water from eastern section of Internal Qinhuai River, representing water quality of Qinhuai scenic zone.
Southern section of IQR
5#
Xiafu Bridge
40–50 cm
Receives sewage discharges and commercial contaminations
Nanhu Lake Area
6#
Saihong Bridge
50–60 cm
A typical black and odorous river near urban arterial road.
Yunliang River Area
7#
Yunliang River
10–30 cm
Receiving the tail water of sewage treatment plant, it is the local Sewage River.
8#
Wetland Park
10–20 cm
Locates between the downstream of Yunliang River and External Qinhuai River, belonging to the main body water of Qiqiaoweng Wetland Park
Upstream of Qinhuai River
9#
External Qinhuai River
20–30 cm
The upstream incoming water of urban water body.
External Qinhuai River
10#
Fengtai Bridge
20–40 cm
Subject to the disturbance of traffic and commercial activities with stable flow volume and concrete revetment.
11#
Hanzhongmen Bridge
12#
Yuhua Bridge
13#
Inside of the Sancha estuary
The estuary of External Qinhuai River into the Yangtze River.
14#
Outside of the Sancha estuary
Xuanwu Lake
15#
Xuanwu Lake
20–30 cm
Receiving the water diversion from drink water treatment plant (1.8 × 105 t/d).
Note. IQR: Internal Qinhuai River.Figure 1
Sampling sites in the study area in Nanjing City, China. The codes for each sampling site are given in Table1.


### 2.2. PAH Extraction, Cleanup, and Determination

For PAHs analysis, sediments were freeze-dried, crushed into fine powders, and passed through a 100-mesh sieve. Five grams of sieved samples was mixed with 2 g anhydrous sodium sulfate and ultrasonic-extracted for 10 min with 15 mLn-hexane: acetone (2 : 1V/V) after soaking for 1 h. The supernatant was transferred to a pear-shaped flask; then the residual sediment was ultrasonic-extracted again with 15 mLn-hexane: dichloromethane (3 : 1V/V) as described above. All the supernatant was combined to a 50 mL pear-shaped flask and concentrated to about 1 mL by vacuum rotary evaporation apparatus (RE-3000, Yarong Biochemical Instrument Plant, Shanghai, China). For purification, the concentrate was transferred to solid phase extraction column (upper: 1.0 g anhydrous sodium sulfate; the second layer is 0.5 g copper, the third layer is 1 cm silica gel, and 1 cm neutral alumina is at bottom) (Borui Co., Ltd., Tianjin, P. R. China) rinsed with 5 mLn-hexane, which had been rinsed with 5 mLn-hexane with a drop-wise flow rate. The column was eluted with a mixture of 25 mLn-hexane: dichloromethane (3 : 1,V :V). The eluent was collected into a pear-shaped flask and concentrated to 1 mL for HPLC determination.Among the reagents, acetone,n-hexane, and copper powder were analytical reagent; anhydrous sodium sulfate was guaranteed reagent grade (Nanjing Chemical Reagent Limited Company, Nanjing, China); neutral alumina was chromatographically pure (Shanghai Ludu Chemical Reagent Plant, Shanghai, China).Test equipment was the high-performance liquid chromatography (HPLC) of Agilent 1100 HPLC with fluorescence (FLD) and UV-adsorption detector (Agilent Technologies, Palo Alto, USA). The column was Agilent ZORBAX EclipseXDB-C18 (4.6 × 250 mm, 5μm). Chromatographic conditions: mobile phase was acetonitrile and water; the initial conditions was acetonitrile: water (6 : 4,V :V) and then linear elution to 100% acetonitrile at the rate of 0.75 mL/min within 60 min. Detection of wavelength was carried out: time programming with FLD signals, 0–9 min: Exλ = 260 nm and Emλ = 380 nm; 9–16 min: Exλ = 260 nm and Emλ = 340 nm; 16–18 min: Exλ = 260 nm and Emλ = 380 nm; 18–21 min: Exλ = 260 nm and Emλ = 380 nm; 21–23 min: Exλ = 289 nm and Emλ = 462 nm; 23–30 min: Exλ = 320 nm and Emλ = 380 nm; 30–36 min: Exλ = 266 nm and Emλ = 403 nm; 36–52 min: Exλ = 294 nm and Emλ = 430 nm; 52 min: Exλ = 290 nm and Emλ = 500 nm. PAHs were quantified by using external standard solutions obtained from Dr. Ehrenstorfer GmbH (Augsburg, Germany).The investigated PAHs in this study included naphthalene (Nap), acenaphthylene (Ace), fluorene (Flu), phenanthrene (Phe), anthracene (Ant), fluoranthene (Fla), pyrene (Pyr), benzo[a]anthracene (BaA), chrysene (Chr), benzo[b]fluoranthene (BbF), benzo[k]fluoranthene (BkF), benzo[a]pyrene (BaP), dibenzo[a,h]anthracene (DahA), and benzo[g,h,i]perylene (BghiP).


### 2.3. Quality Control

In order to ensure the accuracy and precision of experimental data, replicate samples, certified reference materials HS-5 (sediments, provided by NRCIMB of Canada), matrix spikes samples (sediment samples of known PAH levels spiked with a mixture consisting of 1μg each of PAHs), and procedural blanks were used as quality control procedures. Analyses were run in batches of 10 samples plus four quality controls (QCs) including one reagent blank, one matrix blank, one matrix spikes sample, and one random sample in duplicate. No detectable amounts of target PAHs are detected in the blanks. The recovery efficiency of PAHs ranged from 89% to 111% for HS-5 and from 78% to 105% for the spiked samples with relative standard deviation (RSD) 3.2%–10.6%. The method detection limits were in the range of 0.5–2.0 ng g−1 dry weight. The accuracy and precision with PAHs are within 8% variability for both 10 ng/ml and 200 ng/ml concentrations. PAHs standard solution was diluted with acetonitrile; six concentrations of the mixed standard were configured to quantify peak area by establishing a working curve with external standard method, and the correlation of the standard curve was above 0.99. All of the measurements were performed in triplicate, the results were not corrected according to the percent recovery of HS-5, spiked samples, and blanks, and the means were used for calculations.


### 2.4. Positive Matrix Factorization (PMF)

The PMF model, developed by Paatero and Tapper [30], is one of the typical factor analysis models. PMF is a useful factorization method that can calculate source profile and contribution. Unlike the traditional factor analysis models, PMF model utilizes nonnegativity constraints for obtaining physically realistic meanings [31].The principle of the PMF model is to factor the initialn ×m data matrixX (n: number of samples;m: number of species) into (n ×p) matricesG (source contributions) and (p ×m) matricesF (source profiles), as well as the “residual matrix”E (n ×m) firstly, as follows:(1)E=X-GF,whereG is the source contribution matrix,F is the source profile matrix, andp is the number of sources [31]. Certainly, the elements ofG andF are constrained to nonnegative values.The elementseij of matrixE (n×m) indicate the residual value not accounted for by the modeled data value xi.(2)xij=∑k=1pgikfkj+eij,where xij is thejth species concentration measured in the ith sample, gik is the kth source’s contribution to sample i, and fkj is thejth element’s concentration in sourcek. eij is the residual associated with thejth species concentration measured in the ith sample.The objective function(Q) related to the residual and uncertainty is minimized using weighted least-squares by PMF, which is defined as(3)QE=∑i=1m∑j=1neijσij2,wheren andm are the number of samples and species, respectively;p represents the number of factors extracted; i=1,2,…,n samples; j=1,2,…,m species; k=1,2,…,p sources. eij is the difference between the observations and the model; uij is the uncertainty for each observation. The PMF solution minimizes the object function Q based upon the given uncertainty u [32]. The uncertainties for each sample were calculated using measurement uncertainties (MU percent) and method detection limits (MDLs). A discussion of calculating uncertainties is provided by Reff et al. (2007) [33] and USEPA (2008) [21]. If the concentration ⩽ MDL, the uncertaintyu is calculated as [34](4)u=56×MDL,and when the concentration > MDL, u is calculated as(5)u=MU×concentration2+0.5×MDL2.In this study, the EPA PMF 5.0 model was used [21].


## 2.1. Sampling Region and Sites

A map showing the location of the sampling sites is provided in Figure1. Corresponding functional regional characteristics are shown in Table 1. Sediment sampling was carried out along the inner Qinhuai River and Xuanwu Lake in four seasons. The sampling time, based on local seasons, is grouped as follows: autumn (October, 2015), winter (January, 2016), spring (April, 2016), and summer (July, 2016). The middle month of different season was chosen as the typical month of sampling. According to the statistical data estimated by local government, the sedimentation rate is about 10 cm–60 cm a year for different river. Moreover, the sampling tool is the grab sampler, so the sediment samples can be used as representative samples of different seasons. Surface sediment sampling was collected by grab sampler. Each sample consisted of central parts of 5 subsamples randomly collected at 5 points within the study area. The samples were placed into the sampling bottles, sealed and stored at −20°C in the refrigerator for further analysis.Table 1
Sampling sites and the corresponding regional characteristics.
Reach
Code
Position
Sedimentation rate a year
Regional characteristics
Northern section of IQR
1#
Zhenzhu Bridge
30–40 cm
Near the main road and downtown of city, there are floating contaminations.
Eastern section of IQR
2#
Yixian Bridge
Middle of IQR
3#
Wenjin Bridge
50–60 cm
Locates at Chaotian Temple, receiving sewage discharge, with slow water flow due to gate controlling.
4#
Wenzheng Bridge
Receives incoming water from eastern section of Internal Qinhuai River, representing water quality of Qinhuai scenic zone.
Southern section of IQR
5#
Xiafu Bridge
40–50 cm
Receives sewage discharges and commercial contaminations
Nanhu Lake Area
6#
Saihong Bridge
50–60 cm
A typical black and odorous river near urban arterial road.
Yunliang River Area
7#
Yunliang River
10–30 cm
Receiving the tail water of sewage treatment plant, it is the local Sewage River.
8#
Wetland Park
10–20 cm
Locates between the downstream of Yunliang River and External Qinhuai River, belonging to the main body water of Qiqiaoweng Wetland Park
Upstream of Qinhuai River
9#
External Qinhuai River
20–30 cm
The upstream incoming water of urban water body.
External Qinhuai River
10#
Fengtai Bridge
20–40 cm
Subject to the disturbance of traffic and commercial activities with stable flow volume and concrete revetment.
11#
Hanzhongmen Bridge
12#
Yuhua Bridge
13#
Inside of the Sancha estuary
The estuary of External Qinhuai River into the Yangtze River.
14#
Outside of the Sancha estuary
Xuanwu Lake
15#
Xuanwu Lake
20–30 cm
Receiving the water diversion from drink water treatment plant (1.8 × 105 t/d).
Note. IQR: Internal Qinhuai River.Figure 1
Sampling sites in the study area in Nanjing City, China. The codes for each sampling site are given in Table1.


## 2.2. PAH Extraction, Cleanup, and Determination

For PAHs analysis, sediments were freeze-dried, crushed into fine powders, and passed through a 100-mesh sieve. Five grams of sieved samples was mixed with 2 g anhydrous sodium sulfate and ultrasonic-extracted for 10 min with 15 mLn-hexane: acetone (2 : 1V/V) after soaking for 1 h. The supernatant was transferred to a pear-shaped flask; then the residual sediment was ultrasonic-extracted again with 15 mLn-hexane: dichloromethane (3 : 1V/V) as described above. All the supernatant was combined to a 50 mL pear-shaped flask and concentrated to about 1 mL by vacuum rotary evaporation apparatus (RE-3000, Yarong Biochemical Instrument Plant, Shanghai, China). For purification, the concentrate was transferred to solid phase extraction column (upper: 1.0 g anhydrous sodium sulfate; the second layer is 0.5 g copper, the third layer is 1 cm silica gel, and 1 cm neutral alumina is at bottom) (Borui Co., Ltd., Tianjin, P. R. China) rinsed with 5 mLn-hexane, which had been rinsed with 5 mLn-hexane with a drop-wise flow rate. The column was eluted with a mixture of 25 mLn-hexane: dichloromethane (3 : 1,V :V). The eluent was collected into a pear-shaped flask and concentrated to 1 mL for HPLC determination.Among the reagents, acetone,n-hexane, and copper powder were analytical reagent; anhydrous sodium sulfate was guaranteed reagent grade (Nanjing Chemical Reagent Limited Company, Nanjing, China); neutral alumina was chromatographically pure (Shanghai Ludu Chemical Reagent Plant, Shanghai, China).Test equipment was the high-performance liquid chromatography (HPLC) of Agilent 1100 HPLC with fluorescence (FLD) and UV-adsorption detector (Agilent Technologies, Palo Alto, USA). The column was Agilent ZORBAX EclipseXDB-C18 (4.6 × 250 mm, 5μm). Chromatographic conditions: mobile phase was acetonitrile and water; the initial conditions was acetonitrile: water (6 : 4,V :V) and then linear elution to 100% acetonitrile at the rate of 0.75 mL/min within 60 min. Detection of wavelength was carried out: time programming with FLD signals, 0–9 min: Exλ = 260 nm and Emλ = 380 nm; 9–16 min: Exλ = 260 nm and Emλ = 340 nm; 16–18 min: Exλ = 260 nm and Emλ = 380 nm; 18–21 min: Exλ = 260 nm and Emλ = 380 nm; 21–23 min: Exλ = 289 nm and Emλ = 462 nm; 23–30 min: Exλ = 320 nm and Emλ = 380 nm; 30–36 min: Exλ = 266 nm and Emλ = 403 nm; 36–52 min: Exλ = 294 nm and Emλ = 430 nm; 52 min: Exλ = 290 nm and Emλ = 500 nm. PAHs were quantified by using external standard solutions obtained from Dr. Ehrenstorfer GmbH (Augsburg, Germany).The investigated PAHs in this study included naphthalene (Nap), acenaphthylene (Ace), fluorene (Flu), phenanthrene (Phe), anthracene (Ant), fluoranthene (Fla), pyrene (Pyr), benzo[a]anthracene (BaA), chrysene (Chr), benzo[b]fluoranthene (BbF), benzo[k]fluoranthene (BkF), benzo[a]pyrene (BaP), dibenzo[a,h]anthracene (DahA), and benzo[g,h,i]perylene (BghiP).


## 2.3. Quality Control

In order to ensure the accuracy and precision of experimental data, replicate samples, certified reference materials HS-5 (sediments, provided by NRCIMB of Canada), matrix spikes samples (sediment samples of known PAH levels spiked with a mixture consisting of 1μg each of PAHs), and procedural blanks were used as quality control procedures. Analyses were run in batches of 10 samples plus four quality controls (QCs) including one reagent blank, one matrix blank, one matrix spikes sample, and one random sample in duplicate. No detectable amounts of target PAHs are detected in the blanks. The recovery efficiency of PAHs ranged from 89% to 111% for HS-5 and from 78% to 105% for the spiked samples with relative standard deviation (RSD) 3.2%–10.6%. The method detection limits were in the range of 0.5–2.0 ng g−1 dry weight. The accuracy and precision with PAHs are within 8% variability for both 10 ng/ml and 200 ng/ml concentrations. PAHs standard solution was diluted with acetonitrile; six concentrations of the mixed standard were configured to quantify peak area by establishing a working curve with external standard method, and the correlation of the standard curve was above 0.99. All of the measurements were performed in triplicate, the results were not corrected according to the percent recovery of HS-5, spiked samples, and blanks, and the means were used for calculations.


## 2.4. Positive Matrix Factorization (PMF)

The PMF model, developed by Paatero and Tapper [30], is one of the typical factor analysis models. PMF is a useful factorization method that can calculate source profile and contribution. Unlike the traditional factor analysis models, PMF model utilizes nonnegativity constraints for obtaining physically realistic meanings [31].The principle of the PMF model is to factor the initialn ×m data matrixX (n: number of samples;m: number of species) into (n ×p) matricesG (source contributions) and (p ×m) matricesF (source profiles), as well as the “residual matrix”E (n ×m) firstly, as follows:(1)E=X-GF,whereG is the source contribution matrix,F is the source profile matrix, andp is the number of sources [31]. Certainly, the elements ofG andF are constrained to nonnegative values.The elementseij of matrixE (n×m) indicate the residual value not accounted for by the modeled data value xi.(2)xij=∑k=1pgikfkj+eij,where xij is thejth species concentration measured in the ith sample, gik is the kth source’s contribution to sample i, and fkj is thejth element’s concentration in sourcek. eij is the residual associated with thejth species concentration measured in the ith sample.The objective function(Q) related to the residual and uncertainty is minimized using weighted least-squares by PMF, which is defined as(3)QE=∑i=1m∑j=1neijσij2,wheren andm are the number of samples and species, respectively;p represents the number of factors extracted; i=1,2,…,n samples; j=1,2,…,m species; k=1,2,…,p sources. eij is the difference between the observations and the model; uij is the uncertainty for each observation. The PMF solution minimizes the object function Q based upon the given uncertainty u [32]. The uncertainties for each sample were calculated using measurement uncertainties (MU percent) and method detection limits (MDLs). A discussion of calculating uncertainties is provided by Reff et al. (2007) [33] and USEPA (2008) [21]. If the concentration ⩽ MDL, the uncertaintyu is calculated as [34](4)u=56×MDL,and when the concentration > MDL, u is calculated as(5)u=MU×concentration2+0.5×MDL2.In this study, the EPA PMF 5.0 model was used [21].


## 3. Results and Discussion

### 3.1. Temporal and Spatial Characteristics

#### 3.1.1. Temporal and Spatial Characteristics of the ∑PAHs in Sediments

Table2 shows the residue characteristic of 14 PAHs in sediments from Qinhuai River. ∑PAHs concentrations in sediments ranged between 796.2 ng/g and 10,470 ng/g, with a mean of 2,713.8 ng/g. Compared to several studies reported elsewhere (Table 3), mean concentration of ∑PAHs in sediments from Qinhuai River was found to be lower than that from Calabar River (9,370 ng/g), Nigeria [35], and Haihe River, Tianjin, China (27,074.1 ng/g) [36]. The maximum concentration in this study was lower than that in sediments of Haihe River (16,901 ng/g) in 2014 [29], but mean concentration was higher than other rivers. Owing to extensive economic and industrial development in Nigeria, the concentration of ∑PAHs in Calabar River was higher than that in Qinhuai River [35]. The average content (27,074.1 ng/g) of ∑PAHs in the sediment of Haihe River is about ten times that of Qinhuai River (2,713.8 ng/g); this result could be attributed to the fact that Tianjin was an old and established industrial city with the largest port in north China [37], along with rapid development of the petrochemical industry, the iron and steel industry, and other heavy industries and rapid urbanization [36]. In addition to Luan River as a water source with an average of 115.3 ng/g, ∑PAHs levels of other rivers were relatively high, which were closely related to the economic development and industrialization. Qinhuai River is a typical urban river through the densely populated city of Nanjing; sewage discharges and industrial pollution emissions lead to higher mean concentration of ∑PAHs than other rivers apart from Calabar River and Haihe River.Table 2
Mean PAHs concentrations (ng/g), range in surface sediments from Qinhuai River during various seasons.
PAHs (rings)
Autumn
Winter
Spring
Summer
Mean
Range
Mean
Range
Mean
Range
Mean
Range
Nap (2)
892.7
±
534.0
74.3–1486.8
162.4
±
70.7
69.0–272.5
540.0
±
301.8
0–903.0
163.8
±
325.1
0–785.1
Ace (3)
172.5
±
115.3
0–379.9
8.2
±
11.2
0–34.0
143.9
±
316
11.0–1145.2
44.3
±
42.7
0–126.7
Flu (3)
100.3
±
136.1
0.9–529.5
132.1
±
48.7
72.5–211.9
91.0
±
57.9
23.0–248.7
23.7
±
20.8
7.9–67.6
Phe (3)
517.2
±
454.5
71.3–1853.8
336.1
±
190.5
83.0–752.6
369.8
±
479.7
36–1787.7
81.3
±
135.7
0–488.8
Ant (3)
168.1
±
204.2
28.1–814.9
82.1
±
76.9
12.6–284.3
91.2
±
178.1
18.5–653.2
23.4
±
25.6
3.6–105.1
Fla (4)
84.1
±
71.6
0–293.7
21.9
±
24.2
0–88.5
32.5
±
66.2
0–219.7
11.7
±
12.3
0–32.6
Pyr (4)
16.2
±
10.8
0–35.7
12.6
±
8.1
0–32.3
29.4
±
54.3
0–199.9
51.7
±
62.2
0–178
BaA (4)
168.4
±
207.2
28.2–815.7
129.5
±
74.9
25.9–267.5
127.8
±
90.1
22.4–271.1
159.9
±
152.6
33.2–581.6
Chr (4)
543.4
±
426.4
91.6–1,625.6
132.3
±
76.6
27.6–267.9
131.4
±
140.7
19.7–519.4
250.7
±
217.5
0–810.9
BbF (5)
590.6
±
287.4
200.1–1,102.2
313.3
±
237.9
0–726.0
548.4
±
656
92.8–2213
217.1
±
208.1
32.1–846.1
BkF (5)
55.7
±
41.0
12.1–141.9
116.8
±
60.4
33.2–231.6
141.7
±
183.1
7.6–680.8
38.5
±
41.0
6.5–170.2
BaP (5)
132.3
±
117.6
17.3–396.6
223.5
±
126.5
67.8–492.6
243.2
±
332.8
21.1–1,225.1
52.3
±
28.3
10.1–112.2
DahA (5)
216.6
±
156.1
48.4–552.8
173.5
±
81.9
0–300.0
248.3
±
306.3
2.1–1,162.9
161.2
±
129.1
23.4–465.8
BghiP(6)
362.2
±
203.1
49.9–806.6
400.5
±
154.6
158.0–647.5
430.3
±
568.9
0–2,154.9
141.5
±
89.4
36.4–352.6
∑PAHs
4
,
020.1
±
2024.5
2,041.2–9,323.5
2
,
244.9
±
1
,
092.3
796.2–4,269.3
3169.1
±
2
,
492.9
1,416.4–10,470.0
1421.1
±
626
387–2,446.4
∑LMWPAHs
1
,
850.7
±
1
,
012.3
183.9–3,957.6
721.0
±
355.2
250.7–1,396.9
1236
±
566.5
679.8–2,752.2
336.6
±
322.1
65.8–925.5
∑HMWPAHs
2
,
169.4
±
1
,
215.7
608.8–5,365.8
1523.9
±
768.1
545.5–2,909.9
1933.1
±
1
,
991.6
394.1–7,717.8
1084.5
±
540.1
257.3–2,142.7Table 3
Sediment PAHs concentrations reported from elsewhere.
Study areas
Numberof PAHs
∑PAHs Conc.(ng/g)
Mean(ng/g)
References
Qinhuai River, Nanjing, China
14
796.2–10,470
2,713.8
This study
Qiantang River, China
15
91.3–614.4
313.3
Chen et al. 2007 [28]
Songhua River, China
15
20.5–632.0
178.0
Ma et al. 2013 [25]
Songhua River, China
16
68.25–654.15
234.15
Zhao et al. 2014 [23]
Arc River, France
16
151.0–1,257.0
585.3
Kanzari et al. 2012 [40]
Haihe River, China
16
257.0–16,901.0
—
Chen et al. 2014 [29]
Calabar River, Nigeria
17
1,670.0–20,100.0
9,370
Oyo-Ita et al. 2013 [35]
Tuhai-Majia River, China
16
311.7–3,736.3
914.4
Liu et al. 2012 [37]
Haihe River, Tianjin, China
16
774.8–255,371.9
27,074.1
Jiang et al. 2007 [36]
Nilufer Creek, Bursa, Turkey
12
15.0–9,600.0
—
Karaca and Tasdemir, 2014 [41]
Luan River, China
16
20.9–287.0
115.3
Cao et al. 2010 [42]
Yangtze River, Wuhan, China
16
72.4–3,995.2
1,334.5
Feng et al. 2007 [27]
Tonghui River, Beijing, China
16
127.1–927.7
540.4
Zhang et al. 2004 [22]
Rivers in Shanghai, China
17
107.0–1,707.0
823.0
Liu et al. 2008 [24]
Boston Harbor, USA
4
0.4–6300
—
McGroddy et al. 1995 [43]
Malacca Strait, Malaysia
16
347.05–6207.5
—
Keshavarzifard et al. 2017 [44]
Djibouti-city coastal area, Djibouti
16
2.65–3760.11
387.87
Mahdi ahmed et al. [45]
The coastal lagoons and estuaries of New York and New Jersey, USA
—
7.6–6766.1
468.4
Potapova et al. 2016 [46]
Pengkalan Chepa River, Malaysia
16
36.7–560.7
266.5
Noor et al. 2015 [47]
Sydney Harbor, Australia
—
Bd–460336
21782
Birch 2017 [48]As shown in Table2 and Figure 2(b), average concentrations of ∑PAHs in sediments were found to be as follows: autumn (4020.1 ng/g) > spring (3,169.1 ng/g) > winter (2244.9 ng/g) > summer (1421.1 ng/g), indicating that the worst sediment contamination by PAHs appeared in autumn. The uneven seasonal distribution of PAHs in sediments could be attributed to the difference of rainfall and temperature in different seasons. Strong rainfall can scour the surface, carrying the PAHs into the river [38]; high temperature could lead to their subsequent evaporation and degradation [39]. So the least concentration of PAHs appeared in summer.Figure 2
Temporal (b and c) and spatial variation (a) of ∑PAHs in sediments.
(a)
(b)
(c)From Table1 and Figure 2(a), PAHs of different sampling sites were quite spatially different. High levels of PAHs at sampling sites 1#, 3#, 4#, 5#, and 7# may be related to the inputs of surrounding pollution [49]. For example, 7# is located in Yunliang River, receiving the tail water from eastern-city’s sewage treatment plant, while 3# and 4#, located in the middle of Internal Qinhuai River, received polluted incoming water from eastern section of Internal Qinhuai River, contributing to the high PAHs concentrations detected. High concentration PAHs was also found at site 5#, which was in the southern section of Internal Qinhuai River, and the discharges of sewage and markets sewage were observed during the sampling period. The low concentration in site 8#, located in Qiqiaoweng Wetland Park between the downstream of Yunliang River and External Qinhuai River, might be related to the clean incoming water of upper reaches from External Qinhuai River and natural self-cleaning function of wetlands in park. Overall, mean concentrations of ∑PAHs in sediments from Internal Qinhuai River (from sites 1# to 5#) were higher than that from External Qinhuai River, and because Internal Qinhuai River flows through the urban center, frequent human activities lead to more pollution input and the deposition of PAHs in sediments caused by the low-flow rate due to the gate-control and stagnation. It is worth mentioning that the relatively low concentration ∑PAH found at the estuary (site 14#) may be attributed to the fact that this site is far from anthropogenic pressure and the probable tidal flushing of pollutants adsorbed on sediment particles. So, human activities are the main pollution source of PAHs for urban water body. This result is consistent with the view of Oyo-Ita et al. [35].


#### 3.1.2. Temporal and Spatial Variation of ∑PAHs Compositions in Sediments

Figure2(c) shows the seasonal average concentrations of 14 kinds of PAHs in sediments. Mean concentrations of 2-ring Nap, 3-ring Phe, and 5-ring BbF were very high during autumn and summer, while 4-ring Fla and Pyr were low during the four seasons. Overall, the PAHs with 3 and 5 rings were the dominant residual species in the spring and winter, and the concentrations of 3- and 5-ring PAHs were 59.25% and 61.73% of ∑PAHs in spring and winter, respectively; the PAHs with 4 and 5 rings were the dominant residual species in summer, and the percentage contribution was 66.35%, while in autumn they were 2–5 ring PAHs, accounting for 90.99%. The concentration of high molecular weight PAHs (∑HMWPAHs) showed lower variability than that of low molecular weight PAHs (∑LMWPAHs) along with the season change, because the higher solubility of LMWPAHs makes them easily enter the aquatic environment, more susceptible to be affected by the environmental conditions, such as long-distance transmission and evaporation. High temperatures and heavy rainfall in summer contributed to the conduct of this mechanism [39]. The HMWPAHs are relatively stable, are less affected by environmental conditions, and are not easily desorbed from sediments and degraded [50].


### 3.2. Source Apportionment of the PAHs

#### 3.2.1. Source Estimates from Diagnostic Ratios

Diagnostic ratios are often employed to be accurate and reliable diagnostic tools for source identification of PAHs [51, 52]. As isomers of PAHs have similar thermodynamic and kinetic characteristics and isomers from the same source undergo the same migration process, their ratio can be used as the eigenvalues of the source. The composition and content of PAHs vary from different types and conditions of combustion, so the ratio of different PAHs can determine the type of pollution sources. The ratio of low molecular weight and high molecular weight PAHs (L/H) and isomers are commonly used as a ratio method. These groups of compounds had molecular weight of 178 (phenanthrene and anthracene), 202 (fluoranthene and pyrene), 228 (benzo[a]anthracene and chrysene), and 276 (indeno[1,2,3-cd]pyrene and benzo[g,h,i]perylene). In the present study, the ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) were used to identify the PAHs sources in sediment. Ant/(Ant + Phe) ratios of >0.1 are the indicator of pyrogenic origin of PAHs, while <0.1 ratios are petrogenic origin of PAHs [51]. Meanwhile, a Fla/(Fla + Pyr) < 0.4 ratio indicates a petroleum source, a ratio between 0.4 and 0.5 indicates liquid fossil fuel combustion, and a ratio of >0.50 suggests the combustion of biomass and coal [53, 54]. In addition, BaA/(BaA + Chr) ratio is considered as a specific marker for coal combustion [55, 56]. The petroleum source is demonstrated if the ratio of BaA/(BaA + Chr) is lower than 0.20; a source from combustion coal, grass, and wood is shown if the ratio is higher than 0.35; and petroleum combustion (especially liquid fossil fuel and vehicle and crude oil spillage) is demonstrated if the ratio is between 0.20 and 0.35 [51]. For a more detailed analysis of sources, the diagnostic ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) for the four seasons are plotted in Figure 3. As shown in Figure 3, the source of each site was significantly different in summer, while the other three seasons were relatively convergent and consistent. These characteristic ratios demonstrated the PAHs in sediments of Qinhuai River mainly originated from petroleum and combustion of biomass and coal in spring, autumn, and winter is attributed to anthropogenic activities like straw burning before rice cultivation and after harvest, automobiles, and heating in winter and spring. Sources of PAHs in summer were relatively more diverse. Among the sampling sites, M14 was typical petroleum source, M13 and M15 were mixed sources of petroleum source, biomass, coal combustion, and petroleum combustion, and, additionally, M1 and M3 were biomass and coal combustion. Other sampling sites had a mixed source from pyrogenic sources and petroleum sources.Figure 3
Plot of seasonal isomeric ratios Ant/(Ant + Phe) versus BaA/(BaA + Chr) versus Fla/(Fla + Pry). A: autumn, W: winter, S: spring, and M: summer.


#### 3.2.2. Source Estimates from Cluster Analysis

Cluster analysis was performed to identify the homogeneous groups of individual PAHs in the Qinhuai River sediments. The result of cluster analysis is shown in the hierarchical dendrogram (Figure4), which distinguishes the 14 individual PAHs into three major groups. The first group is subdivided into two subgroups. The first subgroup contains BaP, BkF, and BghiP, which are the high molecular weight PAHs with 5-6 rings, usually detected in pyrogenic sources, for example, combustion of biomass, coal, wood, and petroleum. The second subgroup consists of Flu and Phe, which are 3-ring PAHs, from petrogenic sources mainly caused by petroleum spills, for example, fresh or used crankcase oil, crude, and fuel oil [57–61]. The results indicate the first group is mixed pyrogenic sources and petroleum sources. The second major group, which includes Nap, Ace, Ant, and Fla belonging to the low molecular weight PAHs with 2-3 rings, is usually detected in petrogenic sources. The third major group contains Pyr and Chr with 4 rings and BaA, BbF, and DahA with 5 rings, the same as the first subgroup of the first group, usually detected in pyrogenic sources, for example, combustion of coal, wood, vehicle fuel, and waste tire [56, 62, 63]. The results of the second and third groups are consistent with the study by Liu et al. [12], which is attributed to the fact that both Qinhuai River and Huangpu River are urban rivers affected by anthropogenic activities (heating, cooking, and vehicles).Figure 4
Hierarchical dendrogram for 14 PAHs in the Qinhuai River sediments using average linkage between groups and Pearson correlation as measure interval.


#### 3.2.3. Source Estimates from Positive Matrix Factorization (PMF)

For further ascertaining the validity of PAHs sources, fifty-two (52) objects (samples) and fourteen (14) variables (PAHs) together with their uncertainties were used as the input data. Factor numbers ranging from three to eight were considered in order to choose the “optimal solution.” “Robust” mode was adopted to simulate for eliminating the effect of individual extreme value [21, 30, 64–69]. Three factors were eventually chosen since it gave Q (robust) value ranging from 21041.1 to 23360.0, which converged. The slope and r2 values indicated that the predicted and actual concentrations were well fitted and three factors could adequately explain the data based on the background of PAH sources in the sediment.The source identification was based on the molecular markers in different sources of PAHs. As shown in Table4 and Figures 5 and 6, the first factor (f1) is responsible for 33.4% contribution rate, and it is heavily weighted in Ant, BaA, BbF, BkF, BaP, DahA, and BghiP with 4–6 rings. Venkataraman et al. (1994) thought that BkF was typical pollutant of diesel emissions [70]. Ant, BaA, Chr, and BaP were thought to be typical pollutants of biomass or coal combustion [71, 72]. Indeno[1,2,3-cd]pyrene (InP), BghiP, and DahA were typical markers of traffic emission [73, 74]. Hence, this factor represents the traffic-related source > biomass and coal combustion of PAHs, which accounts for the larger proportion of all the sources; Nanjing has a large population and also many private cars and buses. Figure 7 illustrates the seasonal contributions of vehicular emission are spring > winter > autumn > summer. The results are the same as those studied by Liu et al. [12], indicating the choice of transport is temperature-related and people tend to choose cars in cold seasons. Factor 2 (f2) is predominately composed of PAHs with 3–6 rings, including Flu, Phe, Pyr, Chr, BkF, BaP, and BghiP, which mainly come from biomass and coal combustion [75]. Some researches showed that Phe, Ant, Fla, and Pyr were indicative of coal combustion [76]. According to Khalili et al., gasoline combustion was shown to predominantly produce Nap, Flu, Acy, and Pyr [77]. Meanwhile, BghiP and BkF have been identified as tracers for gasoline engines, with Pyr found in both diesel- and gasoline-powered engines [71]. And BaA, Chr, and BaP were thought to be typical pollutants of biomass and coal combustion [71]. Therefore, factor 2 f2 suggests a mixed source of biomass and coal combustion equal to or >diesel and gasoline combustion sources, which is responsible for 25.3% of the total contribution rate throughout the year. Figure 7 shows that winter has the biggest contribution from mixed sources, followed by spring, autumn, and summer. Because winter and spring temperature are relatively low, people will use coal combustion, air conditioners, and other electrical facilities for heating. Nanjing municipal electric power supply is mainly derived from coal-burning power plants. Assigning this factor to coal combustion is demonstrated and shows a seasonal variation. The third factor f3 is predominately composed of Nap, Ace, Phe, Ant, Fla, and Chr with higher loadings. The higher level of Nap and Ace relative to other PAHs suggests a petroleum source [77]. Existing research showed that Ace and acenaphthylene (Acy) are the foremost product of coke burning [78]. And Phe, Ant, Fla, and Chr were indicative of coal combustion [71, 76]. Therefore, factor 3 (f3) was identified as petroleum source > coal combustion [79, 80]. Petroleum source and coal combustion are a major source in the Qinhuai River basin, contributing about 34.9% of PAHs throughout the year. From Figure 7, we can see its seasonal contributions: autumn > spring > winter > summer. Petroleum source may be related to the Yangzi Petrochemical Plant located in the main urban area of Nanjing city.Table 4
Factor profiles (F matrix) of individual PAHs in each of the three PMF factors identified in the analysis of Qinhuai River sediment data set.
PAHs
Factor profiles obtained from PMF (% of species)
Factor 1
Factor 2
Factor 3
Nap
0.000
0.178
0.822
Ace
0.013
0.000
0.987
Flu
0.313
0.580
0.107
Phe
0.106
0.567
0.327
Ant
0.457
0.157
0.386
Fla
0.055
0.076
0.869
Pyr
0.000
1.000
0.000
BaA
0.655
0.163
0.182
Chr
0.043
0.333
0.625
BbF
0.673
0.015
0.312
BkF
0.493
0.500
0.007
BaP
0.477
0.523
0.000
DahA
0.845
0.000
0.155
BghiP
0.547
0.353
0.100
Estimated source
Traffic-related, biomass > coal combustion
Biomass and coal combustion, = or > diesel and gasoline combustion sources
Petroleum source > coal combustion
Rate (%)
33.4
25.3
34.9Figure 5
Factor fingerprints of three PMF factors in each of the individual PAHs.Figure 6
Concentrations and percentage of individual PAHs in each of the three PMF factors.Figure 7
Seasonal contributions of PMF factor profiles.From above we can see that source characterization studies based on the analysis of diagnostic ratio (triangular plot method), cluster analysis, and positive factor matrix analysis [81] suggested that the PAHs of Qinhuai River Basin were mainly from pyrogenic origin (biomass and coal combustion and vehicular emission), and the petroleum source also can not be ignored (specially in summer).


### 3.3. Ecological Risk Assessments of PAHs in Surface Sediments

#### 3.3.1. Sediment Quality Guidelines (SQGs)

PAHs often accumulate in various environmental systems including coastal estuaries and marine sediments [82]. Sediments release PAHs into water as an intrinsic source, which harms aquatic organisms, and then it is harmful to human health, so PAHs of sediments are a potential threat to aquatic ecosystems. Therefore, ecological risk assessment of PAHs in surface sediments from the city-river like Qinhuai River is very important to protect human health and for ecological environment security. The need for chemical guidelines that could be used to predict adverse biological effects in contaminated sediments led to the development of sediment quality guidelines (SQGs) during the last decade [83]. SQGs developed on the basis of Biological Effects Database for Sediments (BEDS) are very useful for sediment assessments in freshwater and marine environments. BEDS contain data that can be used to establish links between the concentration of a given chemical and its biological effect. One of the most widely applied SQGs for estuarine and marine stations is effects range-low (ERL) and effects range-median (ERM) [84]. Threshold effects levels (TELs) and probable effects levels (PELs) are also considered to be widely applicable [85]. These two sets of SQGs can be used to identify three ranges of chemical concentrations, including a low range, within which adverse biological effects are unlikely to occur: below the TEL or ERL values; a middle range, within which adverse biological effects are possible: between the TEL or ERL and PEL or ERM values; and an upper range, within which adverse biological effects are likely to occur: above the PEL or ERM values [42, 84, 85].As shown in Tables5 and 6, all the sampling sites are below ERL values except for 13# and 9# in spring. However, concentrations of Nap, Flu, Phe, Ant, and Chr in autumn; Flu and Phe in winter; Nap, Ace, Flu, Phe, Ant, and BaP in spring; and Nap in summer are all above the PEL values, indicating that these sampling sites are likely to have adverse biological effects on the environment, and the sequence is as follows: autumn > spring > winter > summer. For ∑PAHs, the concentrations in all sites are lower than corresponding PEL values. Relative percentages of sites are 23%, 85%, 67%, and 86% in autumn, winter, spring, and summer, respectively, of which ∑PAHs concentrations are below TEL values, suggesting that the adverse biological effects may be unlikely to occur. In this study, Low molecular weight PAHs (LMWPAHs) concentrations were in the intermediate range (⩾ERL and<ERM), indicating that occasionally adverse biological effects were likely to occur for LMWPAHs at most sites of Qinhuai River during four seasons, which is the same as that studied by Liu et al. [83]. Concentrations levels of high molecular weight PAHs (HMWPAHs) except for Fla and Pyr were in the intermediate range (⩾TEL and<PEL) in most sites, indicating that occasionally adverse biological effects were likely to occur for HMWPAHs. Additionally, there are no corresponding ERL/ERM or TEL/PEL values for BbF, BkF, DahA, and BghiP, showing there will be adverse biological effects on aquatic organisms as long as these four PAHs exist. These four PAHs were all detected at every sampling site of Qinhuai River in four seasons, showing they are a threat to the environment.Table 5
The first set of SQGs values (ERL and ERM) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
ERL
ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
Nap
160
2100
23
77
—
54
46
—
17
83
—
79
21
—
Ace
44
640
23
77
—
100
—
—
33
59
8
57
43
—
Flu
19
540
8
92
—
—
100
—
—
100
—
64
36
—
Phe
240
1500
15
85
—
23
77
—
58
34
8
86
14
—
Ant
85
1100
46
54
—
62
38
—
92
8
—
93
7
—
LMW PAHs
552
3160
23
77
—
23
77
—
—
100
—
71
29
—
Fla
600
5100
100
—
—
100
—
—
100
—
—
100
—
—
Pyr
665
2600
100
—
—
100
—
—
100
—
—
100
—
—
BaA
261
1600
85
15
—
92
8
—
83
17
—
79
21
—
Chr
384
2800
46
54
—
100
—
—
92
8
—
79
21
—
BaP
430
1600
100
—
—
85
15
—
83
17
—
100
—
—
HMW PAHs
1700
9600
92
8
—
100
—
—
92
8
—
100
—
—
Total PAHs
4022
44792
85
15
—
100
—
—
92
8
—
100
—
—Table 6
The second set of SQGs values (TEL and PEL) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
TEL
PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
Nap
34.6
391
—
31
69
—
100
—
17
8
75
79
—
21
Ace
5.87
128
23
—
—
54
46
—
—
92
8
21
79
—
Flu
21.2
144
8
84
8
—
62
38
—
92
8
64
36
—
Phe
86.7
544
8
54
38
8
77
15
8
84
8
71
29
—
Ant
46.9
245
8
77
15
31
69
—
58
34
8
93
7
—
Fla
113
1494
85
15
—
100
—
—
92
8
—
100
—
—
Pyr
153
1398
100
—
—
100
—
—
92
8
—
86
14
—
BaA
75
693
31
69
—
23
77
—
25
75
—
29
71
—
Chr
108
846
8
77
15
31
69
—
67
33
—
21
79
—
BaP
88.8
763
46
54
—
23
77
—
33
59
8
93
7
—
Total PAHs
1684
11770
23
77
—
85
15
—
67
33
—
86
14
—


#### 3.3.2. Sediment Potential Human Toxicity and Biological Effects Based on PAHs

Many PAHs, especially their degradation products, are carcinogenic, while benzo[a]pyrene is the only PAH with sufficient toxicological data for derivation of a carcinogenic potency factor [86]. Therefore, the toxic equivalency factor (TEF) was used to quantify the carcinogenicity of other PAHs relative to benzo[a]pyrene. In the present study, the potential toxicity of sediment was evaluated using the total toxic benzo[a]pyrene equivalent (TEQ) of five PAHs and calculated using the following equation: (6)TEQ=∑Ci×TEFi,where Ci is the concentration of the PAH congener and TEFi is the toxic equivalency factor. According to US EPA [87], TEFs of BaA, BbF, BkF, BaP, and Chr are 0.1, 0.1, 0.1, 1, and 0.001, respectively [88]. Total TEQ values calculated for samples varied from 54.6 (site 13#) to 558.6 (site 7#) ng TEQ•g−1 with an average of 214.3 ng TEQ•g−1 in autumn, 86.7 (station 8#) to 600.4 (station 7#) ng TEQ•g−1, an average of 279.6 ng TEQ•g−1 in winter, 48.6 (site 13#) to 583.0 (site 4#) ng TEQ•g−1, an average of 235.1 ng TEQ•g−1 in spring, and 22.1 (site 14#) to 188.2 (site 1#) ng TEQ•g−1, with an average of 94.1 ng TEQ•g−1. Seasonal mean TEQ values ranged winter > spring > autumn > summer. The results showed that winter had largest impact on the ecological health and the maximum potential toxicity. This is probably due to people’s choice of cars for travelling because of low temperatures and incomplete combustion of fossil fuels. The maximum TEQ value was found at site 7# (Yunliang River) which received the tail water of sewage treatment plant, during autumn and winter. Therefore, direct sewage pollution was the main reason for the distribution of PAHs during autumn and winter in this study. It is reasonable to find that the maximum TEQ value in spring was at site 4#, which received incoming water from eastern section of Internal Qinhuai River, where there are much vehicle traffic and domestic sewage. The maximum TEQ value was found at site 1# in summer; this is attributed to the fact that site 1# is on the main roadside of Nanjing with constant heavy traffic, indicating that the distribution of PAHs of site 1# was affected by traffic. Therefore, we can analyze potential toxicity and ecological effects of the Qinhuai River in different seasons based on the TEQ value and then develop mitigation options based on the results.


## 3.1. Temporal and Spatial Characteristics

### 3.1.1. Temporal and Spatial Characteristics of the ∑PAHs in Sediments

Table2 shows the residue characteristic of 14 PAHs in sediments from Qinhuai River. ∑PAHs concentrations in sediments ranged between 796.2 ng/g and 10,470 ng/g, with a mean of 2,713.8 ng/g. Compared to several studies reported elsewhere (Table 3), mean concentration of ∑PAHs in sediments from Qinhuai River was found to be lower than that from Calabar River (9,370 ng/g), Nigeria [35], and Haihe River, Tianjin, China (27,074.1 ng/g) [36]. The maximum concentration in this study was lower than that in sediments of Haihe River (16,901 ng/g) in 2014 [29], but mean concentration was higher than other rivers. Owing to extensive economic and industrial development in Nigeria, the concentration of ∑PAHs in Calabar River was higher than that in Qinhuai River [35]. The average content (27,074.1 ng/g) of ∑PAHs in the sediment of Haihe River is about ten times that of Qinhuai River (2,713.8 ng/g); this result could be attributed to the fact that Tianjin was an old and established industrial city with the largest port in north China [37], along with rapid development of the petrochemical industry, the iron and steel industry, and other heavy industries and rapid urbanization [36]. In addition to Luan River as a water source with an average of 115.3 ng/g, ∑PAHs levels of other rivers were relatively high, which were closely related to the economic development and industrialization. Qinhuai River is a typical urban river through the densely populated city of Nanjing; sewage discharges and industrial pollution emissions lead to higher mean concentration of ∑PAHs than other rivers apart from Calabar River and Haihe River.Table 2
Mean PAHs concentrations (ng/g), range in surface sediments from Qinhuai River during various seasons.
PAHs (rings)
Autumn
Winter
Spring
Summer
Mean
Range
Mean
Range
Mean
Range
Mean
Range
Nap (2)
892.7
±
534.0
74.3–1486.8
162.4
±
70.7
69.0–272.5
540.0
±
301.8
0–903.0
163.8
±
325.1
0–785.1
Ace (3)
172.5
±
115.3
0–379.9
8.2
±
11.2
0–34.0
143.9
±
316
11.0–1145.2
44.3
±
42.7
0–126.7
Flu (3)
100.3
±
136.1
0.9–529.5
132.1
±
48.7
72.5–211.9
91.0
±
57.9
23.0–248.7
23.7
±
20.8
7.9–67.6
Phe (3)
517.2
±
454.5
71.3–1853.8
336.1
±
190.5
83.0–752.6
369.8
±
479.7
36–1787.7
81.3
±
135.7
0–488.8
Ant (3)
168.1
±
204.2
28.1–814.9
82.1
±
76.9
12.6–284.3
91.2
±
178.1
18.5–653.2
23.4
±
25.6
3.6–105.1
Fla (4)
84.1
±
71.6
0–293.7
21.9
±
24.2
0–88.5
32.5
±
66.2
0–219.7
11.7
±
12.3
0–32.6
Pyr (4)
16.2
±
10.8
0–35.7
12.6
±
8.1
0–32.3
29.4
±
54.3
0–199.9
51.7
±
62.2
0–178
BaA (4)
168.4
±
207.2
28.2–815.7
129.5
±
74.9
25.9–267.5
127.8
±
90.1
22.4–271.1
159.9
±
152.6
33.2–581.6
Chr (4)
543.4
±
426.4
91.6–1,625.6
132.3
±
76.6
27.6–267.9
131.4
±
140.7
19.7–519.4
250.7
±
217.5
0–810.9
BbF (5)
590.6
±
287.4
200.1–1,102.2
313.3
±
237.9
0–726.0
548.4
±
656
92.8–2213
217.1
±
208.1
32.1–846.1
BkF (5)
55.7
±
41.0
12.1–141.9
116.8
±
60.4
33.2–231.6
141.7
±
183.1
7.6–680.8
38.5
±
41.0
6.5–170.2
BaP (5)
132.3
±
117.6
17.3–396.6
223.5
±
126.5
67.8–492.6
243.2
±
332.8
21.1–1,225.1
52.3
±
28.3
10.1–112.2
DahA (5)
216.6
±
156.1
48.4–552.8
173.5
±
81.9
0–300.0
248.3
±
306.3
2.1–1,162.9
161.2
±
129.1
23.4–465.8
BghiP(6)
362.2
±
203.1
49.9–806.6
400.5
±
154.6
158.0–647.5
430.3
±
568.9
0–2,154.9
141.5
±
89.4
36.4–352.6
∑PAHs
4
,
020.1
±
2024.5
2,041.2–9,323.5
2
,
244.9
±
1
,
092.3
796.2–4,269.3
3169.1
±
2
,
492.9
1,416.4–10,470.0
1421.1
±
626
387–2,446.4
∑LMWPAHs
1
,
850.7
±
1
,
012.3
183.9–3,957.6
721.0
±
355.2
250.7–1,396.9
1236
±
566.5
679.8–2,752.2
336.6
±
322.1
65.8–925.5
∑HMWPAHs
2
,
169.4
±
1
,
215.7
608.8–5,365.8
1523.9
±
768.1
545.5–2,909.9
1933.1
±
1
,
991.6
394.1–7,717.8
1084.5
±
540.1
257.3–2,142.7Table 3
Sediment PAHs concentrations reported from elsewhere.
Study areas
Numberof PAHs
∑PAHs Conc.(ng/g)
Mean(ng/g)
References
Qinhuai River, Nanjing, China
14
796.2–10,470
2,713.8
This study
Qiantang River, China
15
91.3–614.4
313.3
Chen et al. 2007 [28]
Songhua River, China
15
20.5–632.0
178.0
Ma et al. 2013 [25]
Songhua River, China
16
68.25–654.15
234.15
Zhao et al. 2014 [23]
Arc River, France
16
151.0–1,257.0
585.3
Kanzari et al. 2012 [40]
Haihe River, China
16
257.0–16,901.0
—
Chen et al. 2014 [29]
Calabar River, Nigeria
17
1,670.0–20,100.0
9,370
Oyo-Ita et al. 2013 [35]
Tuhai-Majia River, China
16
311.7–3,736.3
914.4
Liu et al. 2012 [37]
Haihe River, Tianjin, China
16
774.8–255,371.9
27,074.1
Jiang et al. 2007 [36]
Nilufer Creek, Bursa, Turkey
12
15.0–9,600.0
—
Karaca and Tasdemir, 2014 [41]
Luan River, China
16
20.9–287.0
115.3
Cao et al. 2010 [42]
Yangtze River, Wuhan, China
16
72.4–3,995.2
1,334.5
Feng et al. 2007 [27]
Tonghui River, Beijing, China
16
127.1–927.7
540.4
Zhang et al. 2004 [22]
Rivers in Shanghai, China
17
107.0–1,707.0
823.0
Liu et al. 2008 [24]
Boston Harbor, USA
4
0.4–6300
—
McGroddy et al. 1995 [43]
Malacca Strait, Malaysia
16
347.05–6207.5
—
Keshavarzifard et al. 2017 [44]
Djibouti-city coastal area, Djibouti
16
2.65–3760.11
387.87
Mahdi ahmed et al. [45]
The coastal lagoons and estuaries of New York and New Jersey, USA
—
7.6–6766.1
468.4
Potapova et al. 2016 [46]
Pengkalan Chepa River, Malaysia
16
36.7–560.7
266.5
Noor et al. 2015 [47]
Sydney Harbor, Australia
—
Bd–460336
21782
Birch 2017 [48]As shown in Table2 and Figure 2(b), average concentrations of ∑PAHs in sediments were found to be as follows: autumn (4020.1 ng/g) > spring (3,169.1 ng/g) > winter (2244.9 ng/g) > summer (1421.1 ng/g), indicating that the worst sediment contamination by PAHs appeared in autumn. The uneven seasonal distribution of PAHs in sediments could be attributed to the difference of rainfall and temperature in different seasons. Strong rainfall can scour the surface, carrying the PAHs into the river [38]; high temperature could lead to their subsequent evaporation and degradation [39]. So the least concentration of PAHs appeared in summer.Figure 2
Temporal (b and c) and spatial variation (a) of ∑PAHs in sediments.
(a)
(b)
(c)From Table1 and Figure 2(a), PAHs of different sampling sites were quite spatially different. High levels of PAHs at sampling sites 1#, 3#, 4#, 5#, and 7# may be related to the inputs of surrounding pollution [49]. For example, 7# is located in Yunliang River, receiving the tail water from eastern-city’s sewage treatment plant, while 3# and 4#, located in the middle of Internal Qinhuai River, received polluted incoming water from eastern section of Internal Qinhuai River, contributing to the high PAHs concentrations detected. High concentration PAHs was also found at site 5#, which was in the southern section of Internal Qinhuai River, and the discharges of sewage and markets sewage were observed during the sampling period. The low concentration in site 8#, located in Qiqiaoweng Wetland Park between the downstream of Yunliang River and External Qinhuai River, might be related to the clean incoming water of upper reaches from External Qinhuai River and natural self-cleaning function of wetlands in park. Overall, mean concentrations of ∑PAHs in sediments from Internal Qinhuai River (from sites 1# to 5#) were higher than that from External Qinhuai River, and because Internal Qinhuai River flows through the urban center, frequent human activities lead to more pollution input and the deposition of PAHs in sediments caused by the low-flow rate due to the gate-control and stagnation. It is worth mentioning that the relatively low concentration ∑PAH found at the estuary (site 14#) may be attributed to the fact that this site is far from anthropogenic pressure and the probable tidal flushing of pollutants adsorbed on sediment particles. So, human activities are the main pollution source of PAHs for urban water body. This result is consistent with the view of Oyo-Ita et al. [35].


### 3.1.2. Temporal and Spatial Variation of ∑PAHs Compositions in Sediments

Figure2(c) shows the seasonal average concentrations of 14 kinds of PAHs in sediments. Mean concentrations of 2-ring Nap, 3-ring Phe, and 5-ring BbF were very high during autumn and summer, while 4-ring Fla and Pyr were low during the four seasons. Overall, the PAHs with 3 and 5 rings were the dominant residual species in the spring and winter, and the concentrations of 3- and 5-ring PAHs were 59.25% and 61.73% of ∑PAHs in spring and winter, respectively; the PAHs with 4 and 5 rings were the dominant residual species in summer, and the percentage contribution was 66.35%, while in autumn they were 2–5 ring PAHs, accounting for 90.99%. The concentration of high molecular weight PAHs (∑HMWPAHs) showed lower variability than that of low molecular weight PAHs (∑LMWPAHs) along with the season change, because the higher solubility of LMWPAHs makes them easily enter the aquatic environment, more susceptible to be affected by the environmental conditions, such as long-distance transmission and evaporation. High temperatures and heavy rainfall in summer contributed to the conduct of this mechanism [39]. The HMWPAHs are relatively stable, are less affected by environmental conditions, and are not easily desorbed from sediments and degraded [50].


## 3.1.1. Temporal and Spatial Characteristics of the ∑PAHs in Sediments

Table2 shows the residue characteristic of 14 PAHs in sediments from Qinhuai River. ∑PAHs concentrations in sediments ranged between 796.2 ng/g and 10,470 ng/g, with a mean of 2,713.8 ng/g. Compared to several studies reported elsewhere (Table 3), mean concentration of ∑PAHs in sediments from Qinhuai River was found to be lower than that from Calabar River (9,370 ng/g), Nigeria [35], and Haihe River, Tianjin, China (27,074.1 ng/g) [36]. The maximum concentration in this study was lower than that in sediments of Haihe River (16,901 ng/g) in 2014 [29], but mean concentration was higher than other rivers. Owing to extensive economic and industrial development in Nigeria, the concentration of ∑PAHs in Calabar River was higher than that in Qinhuai River [35]. The average content (27,074.1 ng/g) of ∑PAHs in the sediment of Haihe River is about ten times that of Qinhuai River (2,713.8 ng/g); this result could be attributed to the fact that Tianjin was an old and established industrial city with the largest port in north China [37], along with rapid development of the petrochemical industry, the iron and steel industry, and other heavy industries and rapid urbanization [36]. In addition to Luan River as a water source with an average of 115.3 ng/g, ∑PAHs levels of other rivers were relatively high, which were closely related to the economic development and industrialization. Qinhuai River is a typical urban river through the densely populated city of Nanjing; sewage discharges and industrial pollution emissions lead to higher mean concentration of ∑PAHs than other rivers apart from Calabar River and Haihe River.Table 2
Mean PAHs concentrations (ng/g), range in surface sediments from Qinhuai River during various seasons.
PAHs (rings)
Autumn
Winter
Spring
Summer
Mean
Range
Mean
Range
Mean
Range
Mean
Range
Nap (2)
892.7
±
534.0
74.3–1486.8
162.4
±
70.7
69.0–272.5
540.0
±
301.8
0–903.0
163.8
±
325.1
0–785.1
Ace (3)
172.5
±
115.3
0–379.9
8.2
±
11.2
0–34.0
143.9
±
316
11.0–1145.2
44.3
±
42.7
0–126.7
Flu (3)
100.3
±
136.1
0.9–529.5
132.1
±
48.7
72.5–211.9
91.0
±
57.9
23.0–248.7
23.7
±
20.8
7.9–67.6
Phe (3)
517.2
±
454.5
71.3–1853.8
336.1
±
190.5
83.0–752.6
369.8
±
479.7
36–1787.7
81.3
±
135.7
0–488.8
Ant (3)
168.1
±
204.2
28.1–814.9
82.1
±
76.9
12.6–284.3
91.2
±
178.1
18.5–653.2
23.4
±
25.6
3.6–105.1
Fla (4)
84.1
±
71.6
0–293.7
21.9
±
24.2
0–88.5
32.5
±
66.2
0–219.7
11.7
±
12.3
0–32.6
Pyr (4)
16.2
±
10.8
0–35.7
12.6
±
8.1
0–32.3
29.4
±
54.3
0–199.9
51.7
±
62.2
0–178
BaA (4)
168.4
±
207.2
28.2–815.7
129.5
±
74.9
25.9–267.5
127.8
±
90.1
22.4–271.1
159.9
±
152.6
33.2–581.6
Chr (4)
543.4
±
426.4
91.6–1,625.6
132.3
±
76.6
27.6–267.9
131.4
±
140.7
19.7–519.4
250.7
±
217.5
0–810.9
BbF (5)
590.6
±
287.4
200.1–1,102.2
313.3
±
237.9
0–726.0
548.4
±
656
92.8–2213
217.1
±
208.1
32.1–846.1
BkF (5)
55.7
±
41.0
12.1–141.9
116.8
±
60.4
33.2–231.6
141.7
±
183.1
7.6–680.8
38.5
±
41.0
6.5–170.2
BaP (5)
132.3
±
117.6
17.3–396.6
223.5
±
126.5
67.8–492.6
243.2
±
332.8
21.1–1,225.1
52.3
±
28.3
10.1–112.2
DahA (5)
216.6
±
156.1
48.4–552.8
173.5
±
81.9
0–300.0
248.3
±
306.3
2.1–1,162.9
161.2
±
129.1
23.4–465.8
BghiP(6)
362.2
±
203.1
49.9–806.6
400.5
±
154.6
158.0–647.5
430.3
±
568.9
0–2,154.9
141.5
±
89.4
36.4–352.6
∑PAHs
4
,
020.1
±
2024.5
2,041.2–9,323.5
2
,
244.9
±
1
,
092.3
796.2–4,269.3
3169.1
±
2
,
492.9
1,416.4–10,470.0
1421.1
±
626
387–2,446.4
∑LMWPAHs
1
,
850.7
±
1
,
012.3
183.9–3,957.6
721.0
±
355.2
250.7–1,396.9
1236
±
566.5
679.8–2,752.2
336.6
±
322.1
65.8–925.5
∑HMWPAHs
2
,
169.4
±
1
,
215.7
608.8–5,365.8
1523.9
±
768.1
545.5–2,909.9
1933.1
±
1
,
991.6
394.1–7,717.8
1084.5
±
540.1
257.3–2,142.7Table 3
Sediment PAHs concentrations reported from elsewhere.
Study areas
Numberof PAHs
∑PAHs Conc.(ng/g)
Mean(ng/g)
References
Qinhuai River, Nanjing, China
14
796.2–10,470
2,713.8
This study
Qiantang River, China
15
91.3–614.4
313.3
Chen et al. 2007 [28]
Songhua River, China
15
20.5–632.0
178.0
Ma et al. 2013 [25]
Songhua River, China
16
68.25–654.15
234.15
Zhao et al. 2014 [23]
Arc River, France
16
151.0–1,257.0
585.3
Kanzari et al. 2012 [40]
Haihe River, China
16
257.0–16,901.0
—
Chen et al. 2014 [29]
Calabar River, Nigeria
17
1,670.0–20,100.0
9,370
Oyo-Ita et al. 2013 [35]
Tuhai-Majia River, China
16
311.7–3,736.3
914.4
Liu et al. 2012 [37]
Haihe River, Tianjin, China
16
774.8–255,371.9
27,074.1
Jiang et al. 2007 [36]
Nilufer Creek, Bursa, Turkey
12
15.0–9,600.0
—
Karaca and Tasdemir, 2014 [41]
Luan River, China
16
20.9–287.0
115.3
Cao et al. 2010 [42]
Yangtze River, Wuhan, China
16
72.4–3,995.2
1,334.5
Feng et al. 2007 [27]
Tonghui River, Beijing, China
16
127.1–927.7
540.4
Zhang et al. 2004 [22]
Rivers in Shanghai, China
17
107.0–1,707.0
823.0
Liu et al. 2008 [24]
Boston Harbor, USA
4
0.4–6300
—
McGroddy et al. 1995 [43]
Malacca Strait, Malaysia
16
347.05–6207.5
—
Keshavarzifard et al. 2017 [44]
Djibouti-city coastal area, Djibouti
16
2.65–3760.11
387.87
Mahdi ahmed et al. [45]
The coastal lagoons and estuaries of New York and New Jersey, USA
—
7.6–6766.1
468.4
Potapova et al. 2016 [46]
Pengkalan Chepa River, Malaysia
16
36.7–560.7
266.5
Noor et al. 2015 [47]
Sydney Harbor, Australia
—
Bd–460336
21782
Birch 2017 [48]As shown in Table2 and Figure 2(b), average concentrations of ∑PAHs in sediments were found to be as follows: autumn (4020.1 ng/g) > spring (3,169.1 ng/g) > winter (2244.9 ng/g) > summer (1421.1 ng/g), indicating that the worst sediment contamination by PAHs appeared in autumn. The uneven seasonal distribution of PAHs in sediments could be attributed to the difference of rainfall and temperature in different seasons. Strong rainfall can scour the surface, carrying the PAHs into the river [38]; high temperature could lead to their subsequent evaporation and degradation [39]. So the least concentration of PAHs appeared in summer.Figure 2
Temporal (b and c) and spatial variation (a) of ∑PAHs in sediments.
(a)
(b)
(c)From Table1 and Figure 2(a), PAHs of different sampling sites were quite spatially different. High levels of PAHs at sampling sites 1#, 3#, 4#, 5#, and 7# may be related to the inputs of surrounding pollution [49]. For example, 7# is located in Yunliang River, receiving the tail water from eastern-city’s sewage treatment plant, while 3# and 4#, located in the middle of Internal Qinhuai River, received polluted incoming water from eastern section of Internal Qinhuai River, contributing to the high PAHs concentrations detected. High concentration PAHs was also found at site 5#, which was in the southern section of Internal Qinhuai River, and the discharges of sewage and markets sewage were observed during the sampling period. The low concentration in site 8#, located in Qiqiaoweng Wetland Park between the downstream of Yunliang River and External Qinhuai River, might be related to the clean incoming water of upper reaches from External Qinhuai River and natural self-cleaning function of wetlands in park. Overall, mean concentrations of ∑PAHs in sediments from Internal Qinhuai River (from sites 1# to 5#) were higher than that from External Qinhuai River, and because Internal Qinhuai River flows through the urban center, frequent human activities lead to more pollution input and the deposition of PAHs in sediments caused by the low-flow rate due to the gate-control and stagnation. It is worth mentioning that the relatively low concentration ∑PAH found at the estuary (site 14#) may be attributed to the fact that this site is far from anthropogenic pressure and the probable tidal flushing of pollutants adsorbed on sediment particles. So, human activities are the main pollution source of PAHs for urban water body. This result is consistent with the view of Oyo-Ita et al. [35].


## 3.1.2. Temporal and Spatial Variation of ∑PAHs Compositions in Sediments

Figure2(c) shows the seasonal average concentrations of 14 kinds of PAHs in sediments. Mean concentrations of 2-ring Nap, 3-ring Phe, and 5-ring BbF were very high during autumn and summer, while 4-ring Fla and Pyr were low during the four seasons. Overall, the PAHs with 3 and 5 rings were the dominant residual species in the spring and winter, and the concentrations of 3- and 5-ring PAHs were 59.25% and 61.73% of ∑PAHs in spring and winter, respectively; the PAHs with 4 and 5 rings were the dominant residual species in summer, and the percentage contribution was 66.35%, while in autumn they were 2–5 ring PAHs, accounting for 90.99%. The concentration of high molecular weight PAHs (∑HMWPAHs) showed lower variability than that of low molecular weight PAHs (∑LMWPAHs) along with the season change, because the higher solubility of LMWPAHs makes them easily enter the aquatic environment, more susceptible to be affected by the environmental conditions, such as long-distance transmission and evaporation. High temperatures and heavy rainfall in summer contributed to the conduct of this mechanism [39]. The HMWPAHs are relatively stable, are less affected by environmental conditions, and are not easily desorbed from sediments and degraded [50].


## 3.2. Source Apportionment of the PAHs

### 3.2.1. Source Estimates from Diagnostic Ratios

Diagnostic ratios are often employed to be accurate and reliable diagnostic tools for source identification of PAHs [51, 52]. As isomers of PAHs have similar thermodynamic and kinetic characteristics and isomers from the same source undergo the same migration process, their ratio can be used as the eigenvalues of the source. The composition and content of PAHs vary from different types and conditions of combustion, so the ratio of different PAHs can determine the type of pollution sources. The ratio of low molecular weight and high molecular weight PAHs (L/H) and isomers are commonly used as a ratio method. These groups of compounds had molecular weight of 178 (phenanthrene and anthracene), 202 (fluoranthene and pyrene), 228 (benzo[a]anthracene and chrysene), and 276 (indeno[1,2,3-cd]pyrene and benzo[g,h,i]perylene). In the present study, the ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) were used to identify the PAHs sources in sediment. Ant/(Ant + Phe) ratios of >0.1 are the indicator of pyrogenic origin of PAHs, while <0.1 ratios are petrogenic origin of PAHs [51]. Meanwhile, a Fla/(Fla + Pyr) < 0.4 ratio indicates a petroleum source, a ratio between 0.4 and 0.5 indicates liquid fossil fuel combustion, and a ratio of >0.50 suggests the combustion of biomass and coal [53, 54]. In addition, BaA/(BaA + Chr) ratio is considered as a specific marker for coal combustion [55, 56]. The petroleum source is demonstrated if the ratio of BaA/(BaA + Chr) is lower than 0.20; a source from combustion coal, grass, and wood is shown if the ratio is higher than 0.35; and petroleum combustion (especially liquid fossil fuel and vehicle and crude oil spillage) is demonstrated if the ratio is between 0.20 and 0.35 [51]. For a more detailed analysis of sources, the diagnostic ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) for the four seasons are plotted in Figure 3. As shown in Figure 3, the source of each site was significantly different in summer, while the other three seasons were relatively convergent and consistent. These characteristic ratios demonstrated the PAHs in sediments of Qinhuai River mainly originated from petroleum and combustion of biomass and coal in spring, autumn, and winter is attributed to anthropogenic activities like straw burning before rice cultivation and after harvest, automobiles, and heating in winter and spring. Sources of PAHs in summer were relatively more diverse. Among the sampling sites, M14 was typical petroleum source, M13 and M15 were mixed sources of petroleum source, biomass, coal combustion, and petroleum combustion, and, additionally, M1 and M3 were biomass and coal combustion. Other sampling sites had a mixed source from pyrogenic sources and petroleum sources.Figure 3
Plot of seasonal isomeric ratios Ant/(Ant + Phe) versus BaA/(BaA + Chr) versus Fla/(Fla + Pry). A: autumn, W: winter, S: spring, and M: summer.


### 3.2.2. Source Estimates from Cluster Analysis

Cluster analysis was performed to identify the homogeneous groups of individual PAHs in the Qinhuai River sediments. The result of cluster analysis is shown in the hierarchical dendrogram (Figure4), which distinguishes the 14 individual PAHs into three major groups. The first group is subdivided into two subgroups. The first subgroup contains BaP, BkF, and BghiP, which are the high molecular weight PAHs with 5-6 rings, usually detected in pyrogenic sources, for example, combustion of biomass, coal, wood, and petroleum. The second subgroup consists of Flu and Phe, which are 3-ring PAHs, from petrogenic sources mainly caused by petroleum spills, for example, fresh or used crankcase oil, crude, and fuel oil [57–61]. The results indicate the first group is mixed pyrogenic sources and petroleum sources. The second major group, which includes Nap, Ace, Ant, and Fla belonging to the low molecular weight PAHs with 2-3 rings, is usually detected in petrogenic sources. The third major group contains Pyr and Chr with 4 rings and BaA, BbF, and DahA with 5 rings, the same as the first subgroup of the first group, usually detected in pyrogenic sources, for example, combustion of coal, wood, vehicle fuel, and waste tire [56, 62, 63]. The results of the second and third groups are consistent with the study by Liu et al. [12], which is attributed to the fact that both Qinhuai River and Huangpu River are urban rivers affected by anthropogenic activities (heating, cooking, and vehicles).Figure 4
Hierarchical dendrogram for 14 PAHs in the Qinhuai River sediments using average linkage between groups and Pearson correlation as measure interval.


### 3.2.3. Source Estimates from Positive Matrix Factorization (PMF)

For further ascertaining the validity of PAHs sources, fifty-two (52) objects (samples) and fourteen (14) variables (PAHs) together with their uncertainties were used as the input data. Factor numbers ranging from three to eight were considered in order to choose the “optimal solution.” “Robust” mode was adopted to simulate for eliminating the effect of individual extreme value [21, 30, 64–69]. Three factors were eventually chosen since it gave Q (robust) value ranging from 21041.1 to 23360.0, which converged. The slope and r2 values indicated that the predicted and actual concentrations were well fitted and three factors could adequately explain the data based on the background of PAH sources in the sediment.The source identification was based on the molecular markers in different sources of PAHs. As shown in Table4 and Figures 5 and 6, the first factor (f1) is responsible for 33.4% contribution rate, and it is heavily weighted in Ant, BaA, BbF, BkF, BaP, DahA, and BghiP with 4–6 rings. Venkataraman et al. (1994) thought that BkF was typical pollutant of diesel emissions [70]. Ant, BaA, Chr, and BaP were thought to be typical pollutants of biomass or coal combustion [71, 72]. Indeno[1,2,3-cd]pyrene (InP), BghiP, and DahA were typical markers of traffic emission [73, 74]. Hence, this factor represents the traffic-related source > biomass and coal combustion of PAHs, which accounts for the larger proportion of all the sources; Nanjing has a large population and also many private cars and buses. Figure 7 illustrates the seasonal contributions of vehicular emission are spring > winter > autumn > summer. The results are the same as those studied by Liu et al. [12], indicating the choice of transport is temperature-related and people tend to choose cars in cold seasons. Factor 2 (f2) is predominately composed of PAHs with 3–6 rings, including Flu, Phe, Pyr, Chr, BkF, BaP, and BghiP, which mainly come from biomass and coal combustion [75]. Some researches showed that Phe, Ant, Fla, and Pyr were indicative of coal combustion [76]. According to Khalili et al., gasoline combustion was shown to predominantly produce Nap, Flu, Acy, and Pyr [77]. Meanwhile, BghiP and BkF have been identified as tracers for gasoline engines, with Pyr found in both diesel- and gasoline-powered engines [71]. And BaA, Chr, and BaP were thought to be typical pollutants of biomass and coal combustion [71]. Therefore, factor 2 f2 suggests a mixed source of biomass and coal combustion equal to or >diesel and gasoline combustion sources, which is responsible for 25.3% of the total contribution rate throughout the year. Figure 7 shows that winter has the biggest contribution from mixed sources, followed by spring, autumn, and summer. Because winter and spring temperature are relatively low, people will use coal combustion, air conditioners, and other electrical facilities for heating. Nanjing municipal electric power supply is mainly derived from coal-burning power plants. Assigning this factor to coal combustion is demonstrated and shows a seasonal variation. The third factor f3 is predominately composed of Nap, Ace, Phe, Ant, Fla, and Chr with higher loadings. The higher level of Nap and Ace relative to other PAHs suggests a petroleum source [77]. Existing research showed that Ace and acenaphthylene (Acy) are the foremost product of coke burning [78]. And Phe, Ant, Fla, and Chr were indicative of coal combustion [71, 76]. Therefore, factor 3 (f3) was identified as petroleum source > coal combustion [79, 80]. Petroleum source and coal combustion are a major source in the Qinhuai River basin, contributing about 34.9% of PAHs throughout the year. From Figure 7, we can see its seasonal contributions: autumn > spring > winter > summer. Petroleum source may be related to the Yangzi Petrochemical Plant located in the main urban area of Nanjing city.Table 4
Factor profiles (F matrix) of individual PAHs in each of the three PMF factors identified in the analysis of Qinhuai River sediment data set.
PAHs
Factor profiles obtained from PMF (% of species)
Factor 1
Factor 2
Factor 3
Nap
0.000
0.178
0.822
Ace
0.013
0.000
0.987
Flu
0.313
0.580
0.107
Phe
0.106
0.567
0.327
Ant
0.457
0.157
0.386
Fla
0.055
0.076
0.869
Pyr
0.000
1.000
0.000
BaA
0.655
0.163
0.182
Chr
0.043
0.333
0.625
BbF
0.673
0.015
0.312
BkF
0.493
0.500
0.007
BaP
0.477
0.523
0.000
DahA
0.845
0.000
0.155
BghiP
0.547
0.353
0.100
Estimated source
Traffic-related, biomass > coal combustion
Biomass and coal combustion, = or > diesel and gasoline combustion sources
Petroleum source > coal combustion
Rate (%)
33.4
25.3
34.9Figure 5
Factor fingerprints of three PMF factors in each of the individual PAHs.Figure 6
Concentrations and percentage of individual PAHs in each of the three PMF factors.Figure 7
Seasonal contributions of PMF factor profiles.From above we can see that source characterization studies based on the analysis of diagnostic ratio (triangular plot method), cluster analysis, and positive factor matrix analysis [81] suggested that the PAHs of Qinhuai River Basin were mainly from pyrogenic origin (biomass and coal combustion and vehicular emission), and the petroleum source also can not be ignored (specially in summer).


## 3.2.1. Source Estimates from Diagnostic Ratios

Diagnostic ratios are often employed to be accurate and reliable diagnostic tools for source identification of PAHs [51, 52]. As isomers of PAHs have similar thermodynamic and kinetic characteristics and isomers from the same source undergo the same migration process, their ratio can be used as the eigenvalues of the source. The composition and content of PAHs vary from different types and conditions of combustion, so the ratio of different PAHs can determine the type of pollution sources. The ratio of low molecular weight and high molecular weight PAHs (L/H) and isomers are commonly used as a ratio method. These groups of compounds had molecular weight of 178 (phenanthrene and anthracene), 202 (fluoranthene and pyrene), 228 (benzo[a]anthracene and chrysene), and 276 (indeno[1,2,3-cd]pyrene and benzo[g,h,i]perylene). In the present study, the ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) were used to identify the PAHs sources in sediment. Ant/(Ant + Phe) ratios of >0.1 are the indicator of pyrogenic origin of PAHs, while <0.1 ratios are petrogenic origin of PAHs [51]. Meanwhile, a Fla/(Fla + Pyr) < 0.4 ratio indicates a petroleum source, a ratio between 0.4 and 0.5 indicates liquid fossil fuel combustion, and a ratio of >0.50 suggests the combustion of biomass and coal [53, 54]. In addition, BaA/(BaA + Chr) ratio is considered as a specific marker for coal combustion [55, 56]. The petroleum source is demonstrated if the ratio of BaA/(BaA + Chr) is lower than 0.20; a source from combustion coal, grass, and wood is shown if the ratio is higher than 0.35; and petroleum combustion (especially liquid fossil fuel and vehicle and crude oil spillage) is demonstrated if the ratio is between 0.20 and 0.35 [51]. For a more detailed analysis of sources, the diagnostic ratios of Ant/(Ant + Phe), Fla/(Fla + Pyr), and BaA/(BaA + Chr) for the four seasons are plotted in Figure 3. As shown in Figure 3, the source of each site was significantly different in summer, while the other three seasons were relatively convergent and consistent. These characteristic ratios demonstrated the PAHs in sediments of Qinhuai River mainly originated from petroleum and combustion of biomass and coal in spring, autumn, and winter is attributed to anthropogenic activities like straw burning before rice cultivation and after harvest, automobiles, and heating in winter and spring. Sources of PAHs in summer were relatively more diverse. Among the sampling sites, M14 was typical petroleum source, M13 and M15 were mixed sources of petroleum source, biomass, coal combustion, and petroleum combustion, and, additionally, M1 and M3 were biomass and coal combustion. Other sampling sites had a mixed source from pyrogenic sources and petroleum sources.Figure 3
Plot of seasonal isomeric ratios Ant/(Ant + Phe) versus BaA/(BaA + Chr) versus Fla/(Fla + Pry). A: autumn, W: winter, S: spring, and M: summer.


## 3.2.2. Source Estimates from Cluster Analysis

Cluster analysis was performed to identify the homogeneous groups of individual PAHs in the Qinhuai River sediments. The result of cluster analysis is shown in the hierarchical dendrogram (Figure4), which distinguishes the 14 individual PAHs into three major groups. The first group is subdivided into two subgroups. The first subgroup contains BaP, BkF, and BghiP, which are the high molecular weight PAHs with 5-6 rings, usually detected in pyrogenic sources, for example, combustion of biomass, coal, wood, and petroleum. The second subgroup consists of Flu and Phe, which are 3-ring PAHs, from petrogenic sources mainly caused by petroleum spills, for example, fresh or used crankcase oil, crude, and fuel oil [57–61]. The results indicate the first group is mixed pyrogenic sources and petroleum sources. The second major group, which includes Nap, Ace, Ant, and Fla belonging to the low molecular weight PAHs with 2-3 rings, is usually detected in petrogenic sources. The third major group contains Pyr and Chr with 4 rings and BaA, BbF, and DahA with 5 rings, the same as the first subgroup of the first group, usually detected in pyrogenic sources, for example, combustion of coal, wood, vehicle fuel, and waste tire [56, 62, 63]. The results of the second and third groups are consistent with the study by Liu et al. [12], which is attributed to the fact that both Qinhuai River and Huangpu River are urban rivers affected by anthropogenic activities (heating, cooking, and vehicles).Figure 4
Hierarchical dendrogram for 14 PAHs in the Qinhuai River sediments using average linkage between groups and Pearson correlation as measure interval.


## 3.2.3. Source Estimates from Positive Matrix Factorization (PMF)

For further ascertaining the validity of PAHs sources, fifty-two (52) objects (samples) and fourteen (14) variables (PAHs) together with their uncertainties were used as the input data. Factor numbers ranging from three to eight were considered in order to choose the “optimal solution.” “Robust” mode was adopted to simulate for eliminating the effect of individual extreme value [21, 30, 64–69]. Three factors were eventually chosen since it gave Q (robust) value ranging from 21041.1 to 23360.0, which converged. The slope and r2 values indicated that the predicted and actual concentrations were well fitted and three factors could adequately explain the data based on the background of PAH sources in the sediment.The source identification was based on the molecular markers in different sources of PAHs. As shown in Table4 and Figures 5 and 6, the first factor (f1) is responsible for 33.4% contribution rate, and it is heavily weighted in Ant, BaA, BbF, BkF, BaP, DahA, and BghiP with 4–6 rings. Venkataraman et al. (1994) thought that BkF was typical pollutant of diesel emissions [70]. Ant, BaA, Chr, and BaP were thought to be typical pollutants of biomass or coal combustion [71, 72]. Indeno[1,2,3-cd]pyrene (InP), BghiP, and DahA were typical markers of traffic emission [73, 74]. Hence, this factor represents the traffic-related source > biomass and coal combustion of PAHs, which accounts for the larger proportion of all the sources; Nanjing has a large population and also many private cars and buses. Figure 7 illustrates the seasonal contributions of vehicular emission are spring > winter > autumn > summer. The results are the same as those studied by Liu et al. [12], indicating the choice of transport is temperature-related and people tend to choose cars in cold seasons. Factor 2 (f2) is predominately composed of PAHs with 3–6 rings, including Flu, Phe, Pyr, Chr, BkF, BaP, and BghiP, which mainly come from biomass and coal combustion [75]. Some researches showed that Phe, Ant, Fla, and Pyr were indicative of coal combustion [76]. According to Khalili et al., gasoline combustion was shown to predominantly produce Nap, Flu, Acy, and Pyr [77]. Meanwhile, BghiP and BkF have been identified as tracers for gasoline engines, with Pyr found in both diesel- and gasoline-powered engines [71]. And BaA, Chr, and BaP were thought to be typical pollutants of biomass and coal combustion [71]. Therefore, factor 2 f2 suggests a mixed source of biomass and coal combustion equal to or >diesel and gasoline combustion sources, which is responsible for 25.3% of the total contribution rate throughout the year. Figure 7 shows that winter has the biggest contribution from mixed sources, followed by spring, autumn, and summer. Because winter and spring temperature are relatively low, people will use coal combustion, air conditioners, and other electrical facilities for heating. Nanjing municipal electric power supply is mainly derived from coal-burning power plants. Assigning this factor to coal combustion is demonstrated and shows a seasonal variation. The third factor f3 is predominately composed of Nap, Ace, Phe, Ant, Fla, and Chr with higher loadings. The higher level of Nap and Ace relative to other PAHs suggests a petroleum source [77]. Existing research showed that Ace and acenaphthylene (Acy) are the foremost product of coke burning [78]. And Phe, Ant, Fla, and Chr were indicative of coal combustion [71, 76]. Therefore, factor 3 (f3) was identified as petroleum source > coal combustion [79, 80]. Petroleum source and coal combustion are a major source in the Qinhuai River basin, contributing about 34.9% of PAHs throughout the year. From Figure 7, we can see its seasonal contributions: autumn > spring > winter > summer. Petroleum source may be related to the Yangzi Petrochemical Plant located in the main urban area of Nanjing city.Table 4
Factor profiles (F matrix) of individual PAHs in each of the three PMF factors identified in the analysis of Qinhuai River sediment data set.
PAHs
Factor profiles obtained from PMF (% of species)
Factor 1
Factor 2
Factor 3
Nap
0.000
0.178
0.822
Ace
0.013
0.000
0.987
Flu
0.313
0.580
0.107
Phe
0.106
0.567
0.327
Ant
0.457
0.157
0.386
Fla
0.055
0.076
0.869
Pyr
0.000
1.000
0.000
BaA
0.655
0.163
0.182
Chr
0.043
0.333
0.625
BbF
0.673
0.015
0.312
BkF
0.493
0.500
0.007
BaP
0.477
0.523
0.000
DahA
0.845
0.000
0.155
BghiP
0.547
0.353
0.100
Estimated source
Traffic-related, biomass > coal combustion
Biomass and coal combustion, = or > diesel and gasoline combustion sources
Petroleum source > coal combustion
Rate (%)
33.4
25.3
34.9Figure 5
Factor fingerprints of three PMF factors in each of the individual PAHs.Figure 6
Concentrations and percentage of individual PAHs in each of the three PMF factors.Figure 7
Seasonal contributions of PMF factor profiles.From above we can see that source characterization studies based on the analysis of diagnostic ratio (triangular plot method), cluster analysis, and positive factor matrix analysis [81] suggested that the PAHs of Qinhuai River Basin were mainly from pyrogenic origin (biomass and coal combustion and vehicular emission), and the petroleum source also can not be ignored (specially in summer).


## 3.3. Ecological Risk Assessments of PAHs in Surface Sediments

### 3.3.1. Sediment Quality Guidelines (SQGs)

PAHs often accumulate in various environmental systems including coastal estuaries and marine sediments [82]. Sediments release PAHs into water as an intrinsic source, which harms aquatic organisms, and then it is harmful to human health, so PAHs of sediments are a potential threat to aquatic ecosystems. Therefore, ecological risk assessment of PAHs in surface sediments from the city-river like Qinhuai River is very important to protect human health and for ecological environment security. The need for chemical guidelines that could be used to predict adverse biological effects in contaminated sediments led to the development of sediment quality guidelines (SQGs) during the last decade [83]. SQGs developed on the basis of Biological Effects Database for Sediments (BEDS) are very useful for sediment assessments in freshwater and marine environments. BEDS contain data that can be used to establish links between the concentration of a given chemical and its biological effect. One of the most widely applied SQGs for estuarine and marine stations is effects range-low (ERL) and effects range-median (ERM) [84]. Threshold effects levels (TELs) and probable effects levels (PELs) are also considered to be widely applicable [85]. These two sets of SQGs can be used to identify three ranges of chemical concentrations, including a low range, within which adverse biological effects are unlikely to occur: below the TEL or ERL values; a middle range, within which adverse biological effects are possible: between the TEL or ERL and PEL or ERM values; and an upper range, within which adverse biological effects are likely to occur: above the PEL or ERM values [42, 84, 85].As shown in Tables5 and 6, all the sampling sites are below ERL values except for 13# and 9# in spring. However, concentrations of Nap, Flu, Phe, Ant, and Chr in autumn; Flu and Phe in winter; Nap, Ace, Flu, Phe, Ant, and BaP in spring; and Nap in summer are all above the PEL values, indicating that these sampling sites are likely to have adverse biological effects on the environment, and the sequence is as follows: autumn > spring > winter > summer. For ∑PAHs, the concentrations in all sites are lower than corresponding PEL values. Relative percentages of sites are 23%, 85%, 67%, and 86% in autumn, winter, spring, and summer, respectively, of which ∑PAHs concentrations are below TEL values, suggesting that the adverse biological effects may be unlikely to occur. In this study, Low molecular weight PAHs (LMWPAHs) concentrations were in the intermediate range (⩾ERL and<ERM), indicating that occasionally adverse biological effects were likely to occur for LMWPAHs at most sites of Qinhuai River during four seasons, which is the same as that studied by Liu et al. [83]. Concentrations levels of high molecular weight PAHs (HMWPAHs) except for Fla and Pyr were in the intermediate range (⩾TEL and<PEL) in most sites, indicating that occasionally adverse biological effects were likely to occur for HMWPAHs. Additionally, there are no corresponding ERL/ERM or TEL/PEL values for BbF, BkF, DahA, and BghiP, showing there will be adverse biological effects on aquatic organisms as long as these four PAHs exist. These four PAHs were all detected at every sampling site of Qinhuai River in four seasons, showing they are a threat to the environment.Table 5
The first set of SQGs values (ERL and ERM) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
ERL
ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
Nap
160
2100
23
77
—
54
46
—
17
83
—
79
21
—
Ace
44
640
23
77
—
100
—
—
33
59
8
57
43
—
Flu
19
540
8
92
—
—
100
—
—
100
—
64
36
—
Phe
240
1500
15
85
—
23
77
—
58
34
8
86
14
—
Ant
85
1100
46
54
—
62
38
—
92
8
—
93
7
—
LMW PAHs
552
3160
23
77
—
23
77
—
—
100
—
71
29
—
Fla
600
5100
100
—
—
100
—
—
100
—
—
100
—
—
Pyr
665
2600
100
—
—
100
—
—
100
—
—
100
—
—
BaA
261
1600
85
15
—
92
8
—
83
17
—
79
21
—
Chr
384
2800
46
54
—
100
—
—
92
8
—
79
21
—
BaP
430
1600
100
—
—
85
15
—
83
17
—
100
—
—
HMW PAHs
1700
9600
92
8
—
100
—
—
92
8
—
100
—
—
Total PAHs
4022
44792
85
15
—
100
—
—
92
8
—
100
—
—Table 6
The second set of SQGs values (TEL and PEL) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
TEL
PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
Nap
34.6
391
—
31
69
—
100
—
17
8
75
79
—
21
Ace
5.87
128
23
—
—
54
46
—
—
92
8
21
79
—
Flu
21.2
144
8
84
8
—
62
38
—
92
8
64
36
—
Phe
86.7
544
8
54
38
8
77
15
8
84
8
71
29
—
Ant
46.9
245
8
77
15
31
69
—
58
34
8
93
7
—
Fla
113
1494
85
15
—
100
—
—
92
8
—
100
—
—
Pyr
153
1398
100
—
—
100
—
—
92
8
—
86
14
—
BaA
75
693
31
69
—
23
77
—
25
75
—
29
71
—
Chr
108
846
8
77
15
31
69
—
67
33
—
21
79
—
BaP
88.8
763
46
54
—
23
77
—
33
59
8
93
7
—
Total PAHs
1684
11770
23
77
—
85
15
—
67
33
—
86
14
—


### 3.3.2. Sediment Potential Human Toxicity and Biological Effects Based on PAHs

Many PAHs, especially their degradation products, are carcinogenic, while benzo[a]pyrene is the only PAH with sufficient toxicological data for derivation of a carcinogenic potency factor [86]. Therefore, the toxic equivalency factor (TEF) was used to quantify the carcinogenicity of other PAHs relative to benzo[a]pyrene. In the present study, the potential toxicity of sediment was evaluated using the total toxic benzo[a]pyrene equivalent (TEQ) of five PAHs and calculated using the following equation: (6)TEQ=∑Ci×TEFi,where Ci is the concentration of the PAH congener and TEFi is the toxic equivalency factor. According to US EPA [87], TEFs of BaA, BbF, BkF, BaP, and Chr are 0.1, 0.1, 0.1, 1, and 0.001, respectively [88]. Total TEQ values calculated for samples varied from 54.6 (site 13#) to 558.6 (site 7#) ng TEQ•g−1 with an average of 214.3 ng TEQ•g−1 in autumn, 86.7 (station 8#) to 600.4 (station 7#) ng TEQ•g−1, an average of 279.6 ng TEQ•g−1 in winter, 48.6 (site 13#) to 583.0 (site 4#) ng TEQ•g−1, an average of 235.1 ng TEQ•g−1 in spring, and 22.1 (site 14#) to 188.2 (site 1#) ng TEQ•g−1, with an average of 94.1 ng TEQ•g−1. Seasonal mean TEQ values ranged winter > spring > autumn > summer. The results showed that winter had largest impact on the ecological health and the maximum potential toxicity. This is probably due to people’s choice of cars for travelling because of low temperatures and incomplete combustion of fossil fuels. The maximum TEQ value was found at site 7# (Yunliang River) which received the tail water of sewage treatment plant, during autumn and winter. Therefore, direct sewage pollution was the main reason for the distribution of PAHs during autumn and winter in this study. It is reasonable to find that the maximum TEQ value in spring was at site 4#, which received incoming water from eastern section of Internal Qinhuai River, where there are much vehicle traffic and domestic sewage. The maximum TEQ value was found at site 1# in summer; this is attributed to the fact that site 1# is on the main roadside of Nanjing with constant heavy traffic, indicating that the distribution of PAHs of site 1# was affected by traffic. Therefore, we can analyze potential toxicity and ecological effects of the Qinhuai River in different seasons based on the TEQ value and then develop mitigation options based on the results.


## 3.3.1. Sediment Quality Guidelines (SQGs)

PAHs often accumulate in various environmental systems including coastal estuaries and marine sediments [82]. Sediments release PAHs into water as an intrinsic source, which harms aquatic organisms, and then it is harmful to human health, so PAHs of sediments are a potential threat to aquatic ecosystems. Therefore, ecological risk assessment of PAHs in surface sediments from the city-river like Qinhuai River is very important to protect human health and for ecological environment security. The need for chemical guidelines that could be used to predict adverse biological effects in contaminated sediments led to the development of sediment quality guidelines (SQGs) during the last decade [83]. SQGs developed on the basis of Biological Effects Database for Sediments (BEDS) are very useful for sediment assessments in freshwater and marine environments. BEDS contain data that can be used to establish links between the concentration of a given chemical and its biological effect. One of the most widely applied SQGs for estuarine and marine stations is effects range-low (ERL) and effects range-median (ERM) [84]. Threshold effects levels (TELs) and probable effects levels (PELs) are also considered to be widely applicable [85]. These two sets of SQGs can be used to identify three ranges of chemical concentrations, including a low range, within which adverse biological effects are unlikely to occur: below the TEL or ERL values; a middle range, within which adverse biological effects are possible: between the TEL or ERL and PEL or ERM values; and an upper range, within which adverse biological effects are likely to occur: above the PEL or ERM values [42, 84, 85].As shown in Tables5 and 6, all the sampling sites are below ERL values except for 13# and 9# in spring. However, concentrations of Nap, Flu, Phe, Ant, and Chr in autumn; Flu and Phe in winter; Nap, Ace, Flu, Phe, Ant, and BaP in spring; and Nap in summer are all above the PEL values, indicating that these sampling sites are likely to have adverse biological effects on the environment, and the sequence is as follows: autumn > spring > winter > summer. For ∑PAHs, the concentrations in all sites are lower than corresponding PEL values. Relative percentages of sites are 23%, 85%, 67%, and 86% in autumn, winter, spring, and summer, respectively, of which ∑PAHs concentrations are below TEL values, suggesting that the adverse biological effects may be unlikely to occur. In this study, Low molecular weight PAHs (LMWPAHs) concentrations were in the intermediate range (⩾ERL and<ERM), indicating that occasionally adverse biological effects were likely to occur for LMWPAHs at most sites of Qinhuai River during four seasons, which is the same as that studied by Liu et al. [83]. Concentrations levels of high molecular weight PAHs (HMWPAHs) except for Fla and Pyr were in the intermediate range (⩾TEL and<PEL) in most sites, indicating that occasionally adverse biological effects were likely to occur for HMWPAHs. Additionally, there are no corresponding ERL/ERM or TEL/PEL values for BbF, BkF, DahA, and BghiP, showing there will be adverse biological effects on aquatic organisms as long as these four PAHs exist. These four PAHs were all detected at every sampling site of Qinhuai River in four seasons, showing they are a threat to the environment.Table 5
The first set of SQGs values (ERL and ERM) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
ERL
ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
<ERL
ERL–ERM
>ERM
Nap
160
2100
23
77
—
54
46
—
17
83
—
79
21
—
Ace
44
640
23
77
—
100
—
—
33
59
8
57
43
—
Flu
19
540
8
92
—
—
100
—
—
100
—
64
36
—
Phe
240
1500
15
85
—
23
77
—
58
34
8
86
14
—
Ant
85
1100
46
54
—
62
38
—
92
8
—
93
7
—
LMW PAHs
552
3160
23
77
—
23
77
—
—
100
—
71
29
—
Fla
600
5100
100
—
—
100
—
—
100
—
—
100
—
—
Pyr
665
2600
100
—
—
100
—
—
100
—
—
100
—
—
BaA
261
1600
85
15
—
92
8
—
83
17
—
79
21
—
Chr
384
2800
46
54
—
100
—
—
92
8
—
79
21
—
BaP
430
1600
100
—
—
85
15
—
83
17
—
100
—
—
HMW PAHs
1700
9600
92
8
—
100
—
—
92
8
—
100
—
—
Total PAHs
4022
44792
85
15
—
100
—
—
92
8
—
100
—
—Table 6
The second set of SQGs values (TEL and PEL) for PAHs (ng/g d.w.) and relative percentage of sites including the range of Sediment Quality Guidelines.
Compound
SQG (ng/g)
Autumn (% of sites)
Winter (% of sites)
Spring (% of sites)
Summer (% of sites)
TEL
PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
<TEL
TEL–PEL
>PEL
Nap
34.6
391
—
31
69
—
100
—
17
8
75
79
—
21
Ace
5.87
128
23
—
—
54
46
—
—
92
8
21
79
—
Flu
21.2
144
8
84
8
—
62
38
—
92
8
64
36
—
Phe
86.7
544
8
54
38
8
77
15
8
84
8
71
29
—
Ant
46.9
245
8
77
15
31
69
—
58
34
8
93
7
—
Fla
113
1494
85
15
—
100
—
—
92
8
—
100
—
—
Pyr
153
1398
100
—
—
100
—
—
92
8
—
86
14
—
BaA
75
693
31
69
—
23
77
—
25
75
—
29
71
—
Chr
108
846
8
77
15
31
69
—
67
33
—
21
79
—
BaP
88.8
763
46
54
—
23
77
—
33
59
8
93
7
—
Total PAHs
1684
11770
23
77
—
85
15
—
67
33
—
86
14
—


## 3.3.2. Sediment Potential Human Toxicity and Biological Effects Based on PAHs

Many PAHs, especially their degradation products, are carcinogenic, while benzo[a]pyrene is the only PAH with sufficient toxicological data for derivation of a carcinogenic potency factor [86]. Therefore, the toxic equivalency factor (TEF) was used to quantify the carcinogenicity of other PAHs relative to benzo[a]pyrene. In the present study, the potential toxicity of sediment was evaluated using the total toxic benzo[a]pyrene equivalent (TEQ) of five PAHs and calculated using the following equation: (6)TEQ=∑Ci×TEFi,where Ci is the concentration of the PAH congener and TEFi is the toxic equivalency factor. According to US EPA [87], TEFs of BaA, BbF, BkF, BaP, and Chr are 0.1, 0.1, 0.1, 1, and 0.001, respectively [88]. Total TEQ values calculated for samples varied from 54.6 (site 13#) to 558.6 (site 7#) ng TEQ•g−1 with an average of 214.3 ng TEQ•g−1 in autumn, 86.7 (station 8#) to 600.4 (station 7#) ng TEQ•g−1, an average of 279.6 ng TEQ•g−1 in winter, 48.6 (site 13#) to 583.0 (site 4#) ng TEQ•g−1, an average of 235.1 ng TEQ•g−1 in spring, and 22.1 (site 14#) to 188.2 (site 1#) ng TEQ•g−1, with an average of 94.1 ng TEQ•g−1. Seasonal mean TEQ values ranged winter > spring > autumn > summer. The results showed that winter had largest impact on the ecological health and the maximum potential toxicity. This is probably due to people’s choice of cars for travelling because of low temperatures and incomplete combustion of fossil fuels. The maximum TEQ value was found at site 7# (Yunliang River) which received the tail water of sewage treatment plant, during autumn and winter. Therefore, direct sewage pollution was the main reason for the distribution of PAHs during autumn and winter in this study. It is reasonable to find that the maximum TEQ value in spring was at site 4#, which received incoming water from eastern section of Internal Qinhuai River, where there are much vehicle traffic and domestic sewage. The maximum TEQ value was found at site 1# in summer; this is attributed to the fact that site 1# is on the main roadside of Nanjing with constant heavy traffic, indicating that the distribution of PAHs of site 1# was affected by traffic. Therefore, we can analyze potential toxicity and ecological effects of the Qinhuai River in different seasons based on the TEQ value and then develop mitigation options based on the results.


## 4. Conclusions

In this study, we chose 15 typical sites in the Qinhuai River and Xuanwu Lake during four seasons. The distribution, sources, and risk assessment of PAHs in sediment are studied using different methods. The comparisons of different methods (diagnostic ratio, cluster analysis, and positive factor matrix analysis) for source analysis confirmed the PAHs of Qinhuai River are primarily from pyrogenic origin (biomass and coal combustion and vehicular emission), and the petroleum source also can not be ignored (specially in summer), and the PAHs pollution sources of urban water body have obvious seasonal-dependent and human activities-dependent characteristics. The result of ecological risk assessment shows the highest benzo[a]pyrene-equivalent dose (BaPeq dose) appears at the sites of sewage discharge and heavy traffic, which should be the focus of continued study.



---

*Source: 3510796-2017-11-15.xml*