Health Risk Assessment of PAHs from Estuarine Sediments in the South of Italy

Increased concerns about the toxicities of Polycyclic Aromatic Hydrocarbons (PAHs), ubiquitous and persistent compounds, as well as the associated ecotoxicology issue in estuarine sediments, have drawn attention worldwide in the last few years. The levels of PAHs in the Sele, Sarno, and Volturno Rivers sediments were evaluated. Moreover, the cancerogenic risk resulting from dermal and ingestion exposure to PAHs was estimated using the incremental lifetime cancer risk (ILCR) assessment and the toxic equivalent concentration (TEQBaP). For Sele River, the results showed that the total PAH concentration ranged from 632.42 to 844.93 ng g−1 dw, with an average value of 738.68 ng g−1 dw. ∑PAHs were in the range of 5.2–678.6 ng g−1 dw and 434.8–872.1 ng g−1 dw for the Sarno and Volturno River sediments, respectively. The cancerogenic risk from the accidental ingestion of PAHs in estuarine sediments was low at all sampling sites. However, based on the ILCRdermal values obtained, the risk of cancer associated with exposure by dermal contact with the PAHs present in the sediments was moderate, with a mean ILCRdermal value of 2.77 × 10−6. This study revealed the pollution levels of PAHs across the South of Italy and provided a scientific basis for PAH pollution control and environmental protection.


Introduction
Estuaries are the main deposits for the disposal of industrial and domestic effluents, sewage sludge, and dredged material with a significant load of contaminants, including PAHs, from pipeline discharges, vehicular emissions, atmospheric deposition, surface runoff, as well as oil spills in aquatic environments [1]. Due to their low water solubility and high lipophilicity, PAHs tend to accumulate in the sediments of aquatic systems for long periods due to their high degradation resistance and high organic carbon content [2]. In soil and sediment compartments, PAHs can undergo biodegradation processes by microorganisms. However, due to their stable physic-chemical characteristics, hydrophobicity, and a strong tendency to absorb into the soil matrix, their biodegradation rate is low. Consequently, PAHs do not degrade easily, and they can accumulate in the solid phase of the terrestrial and aquatic environment, where they persist for a long time [3,4]. Thus, sediments constitute a natural reserve of PAHs in the aquatic system [5,6]. Moreover, they can be released into the surrounding environment by means of resuspension phenomena, thus giving rise to "secondary pollution". Consequently, since high concentrations of PAHs in sediments may reveal a potential pollution risk for the environment and human health, it is essential to monitor these compounds in sediment to protect and preserve the aquatic environment and human health [7].
PAHs, as persistent organic pollutants (POPs), are widely present in the ecosystem [8][9][10]. Their spread to the environment has raised many concerns for human health, as some of them long-term exposure to low levels of Pyr and BaP has been identified as the cause of cancer in laboratory animals [32]. Animal studies have also shown adverse effects on reproduction and development due to exposure to PAHs, whereas these effects were not commonly detected in humans [33,34]. Moreover, Anyahara stated that exposure to PAHs can induce cataracts and cause kidney and liver damage and jaundice [35].
Estuaries are important aquatic systems largely affected by PAH pollution. In fact, due to their chemical properties, PAHs can persist in water where they are readily adsorbed onto particulate matter, settling in river sediments and soils. Thus, river sediments behave as the primary sink and reservoir for PAHs in the aquatic environment [36,37]. Consequently, sediments are important indicators since they can reflect the pollution status of the environment [38].
To date, no previous studies have evaluated the carcinogenic risk to human health associated with dermal and accidental ingestion exposure to PAHs from surface sediments in the South of Italy. Therefore, the main aim of this study was to assess the risk to human health from exposure to PAHs present in the sediments of surface waters in a large coastal area of the Campania Region, in the South of Italy. Specifically, the purpose of this study was to evaluate the distribution patterns of PAHs and to assess the carcinogenic risk to human health from dermal and ingestion exposure to these contaminants from the estuarine sediments of the Sarno, Volturno, and Sele Rivers, which are the main surface water streams of the Campania Region, in the South of Italy.

Study Area
This assessment of the human health risk from exposure to PAHs was carried out in a study area of approximately 3100 km 2 and included the three largest plains in the Campania Region, in the South of Italy. Particularly, the research area was close to the estuaries of the Sele, Volturno, and Sarno Rivers, which traverse the same-named plains. Figure 1 shows the three plains of interest and the respective estuaries of the rivers that cross them. skin irritants [26]. In addition, PAHs interfere with hormonal systems and, as a result, can have harmful effects on reproduction and immune function [27]. The adverse effects of exposure to PAHs have been extensively investigated, but the information currently available on human exposure to individual PAHs is scattered and incomplete, except for some accidental contact with NaP and BaP [28][29][30]. Srogi et al. stated that prolonged dermal contact with NaP may cause redness and inflammation of the skin [31]. In addition, Diggs et al. reported that long-term exposure to low levels of Pyr and BaP has been identified as the cause of cancer in laboratory animals [32]. Animal studies have also shown adverse effects on reproduction and development due to exposure to PAHs, whereas these effects were not commonly detected in humans [33,34]. Moreover, Anyahara stated that exposure to PAHs can induce cataracts and cause kidney and liver damage and jaundice [35].
Estuaries are important aquatic systems largely affected by PAH pollution. In fact, due to their chemical properties, PAHs can persist in water where they are readily adsorbed onto particulate matter, settling in river sediments and soils. Thus, river sediments behave as the primary sink and reservoir for PAHs in the aquatic environment [36,37]. Consequently, sediments are important indicators since they can reflect the pollution status of the environment [38].
To date, no previous studies have evaluated the carcinogenic risk to human health associated with dermal and accidental ingestion exposure to PAHs from surface sediments in the South of Italy. Therefore, the main aim of this study was to assess the risk to human health from exposure to PAHs present in the sediments of surface waters in a large coastal area of the Campania Region, in the South of Italy. Specifically, the purpose of this study was to evaluate the distribution patterns of PAHs and to assess the carcinogenic risk to human health from dermal and ingestion exposure to these contaminants from the estuarine sediments of the Sarno, Volturno, and Sele Rivers, which are the main surface water streams of the Campania Region, in the South of Italy.

Study Area
This assessment of the human health risk from exposure to PAHs was carried out in a study area of approximately 3100 km 2 and included the three largest plains in the Campania Region, in the South of Italy. Particularly, the research area was close to the estuaries of the Sele, Volturno, and Sarno Rivers, which traverse the same-named plains. Figure 1 shows the three plains of interest and the respective estuaries of the rivers that cross them.

Sampling
A sampling campaign was carried out during the spring season (April 2021) at 10 sampling sites near the Sele River Estuary. In detail, sediment samples were taken from the mouth of the river (Site 1) at different distances, 500 m, 1000 m, and 1500 m from the mouth, and directions, to the north, west, and south of the estuary (Figure 2). During sampling, a global positioning system (GPS) was used to locate all sampling sites.
Information on the identification number (ID), characteristics, and coordinates of each sampling location are shown in Table 1. The samples were collected at a depth of 0 to 5 cm using a scraping sampler (Van Veen Grab) and placed in aluminum containers. Then, they were transferred under refrigeration to the laboratory and stored at −20 • C until analysis.

Sampling
A sampling campaign was carried out during the spring season (April 2021) at 10 sampling sites near the Sele River Estuary. In detail, sediment samples were taken from the mouth of the river (Site 1) at different distances, 500 m, 1000 m, and 1500 m from the mouth, and directions, to the north, west, and south of the estuary (Figure 2). During sampling, a global positioning system (GPS) was used to locate all sampling sites. Information on the identification number (ID), characteristics, and coordinates of each sampling location are shown in Table 1. The samples were collected at a depth of 0 to 5 cm using a scraping sampler (Van Veen Grab) and placed in aluminum containers. Then, they were transferred under refrigeration to the laboratory and stored at −20 °C until analysis.  The PAHs levels in sediment samples from the Sarno and Volturno Rivers were evaluated previously [39,40]. Briefly, for the Sarno River, a sampling campaign was carried out during the spring of 2008 at the source of the river (site 1), just before and after the junction with Alveo Comune, at the river mouth (site 4), and in 9 sites located at different distances from the estuary ( Figure 3). More detailed information about the sediment sampling in the Sarno River is given in Table 2.  The PAHs levels in sediment samples from the Sarno and Volturno Rivers were evaluated previously [39,40]. Briefly, for the Sarno River, a sampling campaign was carried out during the spring of 2008 at the source of the river (site 1), just before and after the junction with Alveo Comune, at the river mouth (site 4), and in 9 sites located at different distances from the estuary (Figure 3). More detailed information about the sediment sampling in the Sarno River is given in Table 2.    For the Volturno River, the sampling campaign was carried out in April 2018 near the mouth of the Volturno River and in 9 sites located at different distances from it ( Figure 4). The specific details are provided in Table 3.  For the Volturno River, the sampling campaign was carried out in April 2018 near the mouth of the Volturno River and in 9 sites located at different distances from it ( Figure  4). The specific details are provided in Table 3.

Extraction Procedure and Clean-Up
The analyses were performed as described previously [41]. Briefly, for PAH extraction, the sediment samples were air-dried, crushed, sieved in 250 µm particles, and then divided into portions of 5 g. The PAH concentrations were indicated as dry weight (ng/g dw) [42,43]. The PAH extraction was performed with a Soxhlet extractor using methylene chloride as solvent. Subsequently, the extracts, first purified using a column composed of sodium sulfate/silica gel and then eluted with 70 mL of a hexane:methylene chloride (7:3, v/v) solution, were evaporated to dryness and reduced to a final volume (500 µL) with the aid of a weak current of nitrogen. Finally, the extracts were analyzed using gas chromatography coupled with mass spectrometry (GC-MS). A TOC analyzer was used to evaluate the total organic carbon (TOC) content in the sediment samples (TOC-VCPH, Shimadzu Corp., Kyoto, Japan).

Instrumental Analysis, Quality Assurance, and Quality Control
A TRACE TM 1310 gas chromatograph coupled to an ISQ TM 7000 single quadrupole mass spectrometer (GC-MS, Thermo Scientific, Waltham, MA, USA) was used, equipped with a capillary column TG-5MS (length 30 mm, inner diameter 0.25 mm, film thickness 0.25 µm) and helium as a gas carrier (constant flow of 1 mL/min), operating in the electronic ionization mode (EI) set to 70 eV. The injector operated at 280 • C, and the temperature of the detector was set to 300 • C. The acquisition was performed with the Selected Ion Monitoring (SIM) mode using two characteristic fragments for each selected analyte. A splitless injection mode was adopted with an injection volume of 1 µL. The quantification of PAHs was carried out using response factors related to the respective internal standards based on a six-point calibration curve for individual PAHs (Dr. Ehrenstorfer GmbH, Augsburg, Germany) (R 2 > 0.97). Chrysene-d 12 was used as an internal standard for sample quantification. Before the analysis, all the glassware to be used was thoroughly washed with methanol, acetone, and dichloromethane and placed in the oven at 200 • C to minimize possible sources of contamination. The column temperature was set with different gradients: from 60 • C to 200 • C with an increase of 25 • C/min (kept for 2 min), to 270 • C increasing at 10 • C min −1 (kept for 6 min), and to 310 • C with a rise of 25 • C min −1 (kept for 10 min). The single ion monitoring mode (SIM) was used for the acquisition using characteristic ions for each target analyte. The 16 priority IPA, according to the WHO and USEPA, were evaluated (Table 4) [20,44].
Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively. They were in the range of 1.5-1.9 ng g −1 and 5.1-6.3 ng g −1 , respectively. In the procedural blanks, analyzed as the samples, the PAHs showed a concentration below the LOD. Moreover, for individual PAHs, the recovery test values ranged from 80% to 97%, meeting the quality control criteria (70-130%). For the effective and reproducible detection and quantification of low concentrations of PAHs in sediments, the linear range, precision, limits of detection, and limits of quantification were performed. The precision of the method was determined using repeatability tests and was expressed as standard deviation (SD) (Tables S1-S3). The average of the results was used to estimate the precision of the method.

Human Health Risk Assessment
The human health risk assessment is useful in determining whether exposure to a chemical in a specific dose may cause an increase in the frequency of adverse effects on human health [39,42,[45][46][47]. For the population, PAH exposure represents a health and hygiene risk, which is assessed as a carcinogenic risk. Therefore, the USEPA defined carcinogenic risk as the probability that an individual may develop cancer over a lifetime from exposure to a specific substance classified as mutagenic or carcinogenic. Thus, this risk assessment consists of two phases, which are the estimation of the probability of an event occurring and the study of the likely magnitude of its adverse effect over a specific time frame [48].  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.

Crysene
Chr C 18 H 12 4 Toxics 2023, 11, x FOR PEER REVIEW 7 of 18 Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.  Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.   [20,44]. Six-point calibration curves (5-10-50-250-500-1000 ng/L), procedural blanks, and sample triplicates were carried out for every set of samples. The PAH concentrations were calculated as dry weight (ng/g dw). The limits of detection (LOD) and quantification (LOQ) were evaluated as three and ten times the noise in blank samples, respectively.

g mol −1
The incremental lifetime cancer risk (ILCR) due to exposure by direct ingestion and skin contact to PAHs present in the sediment was evaluated [49]. First, the doses of contaminants taken up by human receptors through the two different exposure pathways considered were calculated according to Equations (1) and (2) [50]: where: Dose ingestion (mg/kg-day) indicates the dose from accidental sediment ingestion; Dose dermal (mg/kg-day) is the dose from skin contact with sediment; C s (mg/kg) represents the concentration of the contaminant in the sediment; IR s (kg/day) is the rate of accidental sediment ingestion; RAF oral indicates the relative absorption factor for the gastrointestinal tract; RAF derm (dimensionless) expresses the relative absorption factor for the skin. Moreover, the dose was evaluated based on the hours per day with exposure: 0-16/16 h for accidental ingestion of sediment (D hours ); days in a week with exposure [(0-7)/7 days] (D days ); weeks in a year with exposure [(0-52)/52 weeks] (D weeks ); total years with exposure (ED years ); surface of hands (assuming only hands are exposed) (SA h (cm 2 ); SL h (kg/cm 2 -event) = Sediment load rate on exposed skin; EF (event/day) = Number of skin exposures per day; BW (kg) = receptor body weight; LE = life expectancy/average life expectancy expressed in years; CF (conversion coefficient) = 1 × 10 −6 kg/mg. Additionally, since Benzo(a)pyrene (BaP) is the most carcinogenic compound among the PAHs considered [10], all individual analyte concentrations were converted to the corresponding toxic equivalent concentrations of BaP. These concentrations are referred to as TEQ BaP or BaP eq and were obtained using the concentration product for toxic equivalence factor (TEF). The TEF factors for the 16 US EPA priority PAHs are shown in Table 5 [51]. The total concentrations of PAHs were obtained using the sum of the calculated toxic equivalents for each compound in relation to BaP. BaPeq were calculated according to Equation (3) [52]: where C s represents the average concentration of an individual PAH. Humans can encounter the PAHs present in estuarine sediments by oral ingestion and dermal contact. As a result, in addition to calculating the doses of contaminants taken up Toxics 2023, 11, 172 9 of 17 by human receptors, the incremental lifetime cancer risk by oral ingestion (ILCR ingestion ) and dermal contact (ILCR dermal ) was evaluated according to Equations (4) and (5) [49]: where: SF ingestion (Kg-day/mg) indicates the oral slope factor. SF dermal (Kg-day/mg) is the dermal slope factor. SA (cm 2 /kg) represents the area of dermal contact with the sediment. AF (mg/cm 2 ) is the skin absorption coefficient for the sediment. ABS is the skin absorption coefficient for contaminants. AT (years) is the average lifespan.
The parameters used for the calculation of the carcinogenic risk are shown in Table 6.

PAH Concentrations in Sediment from the Sele River
The concentrations of the 16 USEPA priority PAHs obtained from instrumental analyses of sediment samples taken near the mouth of the Sele River are given in Table S1. In particular, the total concentration of the PAHs ranged from 632.42 ng g −1 dw (site 10) to 844.93 ng g −1 dw (at site 1), with an average value of 738.68 ng g −1 dw. Specifically, the concentrations ranged from 2.23 to 70.64 ng g −1 dw with an average value of 36.43 ng g −1 dw for PAHs with 2 rings (NaP), from 5.45 to 51.03 ng g −1 dw for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 0.70 to 74.6 ng g −1 dw for 4-ring PAHs (Flu, Pyr, BaA, Chr), from 39.12 to 154.99 ng g −1 dw for 5-ring PAHs (BbF, BkF, BaP, DahA), and from 3.01 to 75.13 ng g −1 dw for 6-ring PAHs (BghiP, IcdP). Figure 5 shows the individual concentrations of PAHs detected in sediment samples from different sampling sites. The figure reveals that the highest concentrations of PAHs were found at the mouth of the Sele River (site 1) and 500 m from the mouth in the southerly direction (site 8). The composition profile of PAHs in the sediment is shown in Figure 6. PAHs with 5 rings were found in most test sites at 57.4% of the total PAHs in the sediment.
concentrations ranged from 2.23 to 70.64 ng g −1 dw with an average value of 36.43 ng g −1 dw for PAHs with 2 rings (NaP), from 5.45 to 51.03 ng g −1 dw for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 0.70 to 74.6 ng g −1 dw for 4-ring PAHs (Flu, Pyr, BaA, Chr), from 39.12 to 154.99 ng g −1 dw for 5-ring PAHs (BbF, BkF, BaP, DahA), and from 3.01 to 75.13 ng g −1 dw for 6-ring PAHs (BghiP, IcdP). Figure 5 shows the individual concentrations of PAHs detected in sediment samples from different sampling sites. The figure reveals that the highest concentrations of PAHs were found at the mouth of the Sele River (site 1) and 500 m from the mouth in the southerly direction (site 8). The composition profile of PAHs in the sediment is shown in Figure 6. PAHs with 5 rings were found in most test sites at 57.4% of the total PAHs in the sediment.

PAH Concentrations in Sediment from the Sarno River
Data on PAH concentrations found in the Sarno River are indicated in Table S2 [39]. The total concentration of PAHs in the sediment ranged from 5.2 ng g −1 dw at the source of the river (site 1) to 678.6 ng g −1 dw at the point 150 m to the west of the mouth (site 9), with an average value of 266.9 ng g −1 dw. The measured PAH concentrations ranged from 0.2 to 31.6 ng g −1 dw with an average of 9.7 ng g −1 dw for 2-ring PAHs (Nap), from 0.2 to 46.3 ng g −1 dw for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 0.3 to 47.2 ng g −1 dw for 4ring PAHs (Fla, Pyr, BaA, Chr), from 0.2 to 46.6 ng g −1 dw for 5-ring PAHs (BbF, BkF, BaP, DahA), and from 0.5 to 46.7 ng g −1 dw for 6-ring PAHs (BghiP, IcdP). Figure 7 shows the individual and total concentrations of PAHs found in the sediment samples taken from the different sampling sites located near the mouth of the Sarno River. The figure indicates an increase in total PAH levels at the sampling point 150 m west of the river mouth (site 8). The composition profile of PAHs in the sediment is shown in Figure 8. Three-ring

PAH Concentrations in Sediment from the Sarno River
Data on PAH concentrations found in the Sarno River are indicated in Table S2 [39]. The total concentration of PAHs in the sediment ranged from 5.2 ng g −1 dw at the source of the river (site 1) to 678.6 ng g −1 dw at the point 150 m to the west of the mouth (site 9), with an average value of 266.9 ng g −1 dw. The measured PAH concentrations ranged from 0.2 to 31.6 ng g −1 dw with an average of 9.7 ng g −1 dw for 2-ring PAHs (Nap), from 0.2 to 46.3 ng g −1 dw for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 0.3 to 47.2 ng g −1 dw for 4-ring PAHs (Fla, Pyr, BaA, Chr), from 0.2 to 46.6 ng g −1 dw for 5-ring PAHs (BbF, BkF, BaP, DahA), and from 0.5 to 46.7 ng g −1 dw for 6-ring PAHs (BghiP, IcdP). Figure 7 shows the individual and total concentrations of PAHs found in the sediment samples taken from the different sampling sites located near the mouth of the Sarno River. The figure indicates an increase in total PAH levels at the sampling point 150 m west of the river mouth (site 8). The composition profile of PAHs in the sediment is shown in Figure 8. Three-ring PAHs were found in most test sites, at a percentage of 47.3% of the total PAH amount in the sediment, followed by 5-ring PAHs at a percentage of 20.6%. of the river (site 1) to 678.6 ng g −1 dw at the point 150 m to the west of the mouth (site 9), with an average value of 266.9 ng g −1 dw. The measured PAH concentrations ranged from 0.2 to 31.6 ng g −1 dw with an average of 9.7 ng g −1 dw for 2-ring PAHs (Nap), from 0.2 to 46.3 ng g −1 dw for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 0.3 to 47.2 ng g −1 dw for 4ring PAHs (Fla, Pyr, BaA, Chr), from 0.2 to 46.6 ng g −1 dw for 5-ring PAHs (BbF, BkF, BaP, DahA), and from 0.5 to 46.7 ng g −1 dw for 6-ring PAHs (BghiP, IcdP). Figure 7 shows the individual and total concentrations of PAHs found in the sediment samples taken from the different sampling sites located near the mouth of the Sarno River. The figure indicates an increase in total PAH levels at the sampling point 150 m west of the river mouth (site 8). The composition profile of PAHs in the sediment is shown in Figure 8. Three-ring PAHs were found in most test sites, at a percentage of 47.3% of the total PAH amount in the sediment, followed by 5-ring PAHs at a percentage of 20.6%.

PAH Concentrations in Sediment from the Volturno River
Data on individual PAH concentrations found in the Volturno River are given in Table S3, while total concentrations were previously reported [40]. In detail, total concentrations were between 434.8 ng g −1 dw (site 8) and 872.1 ng g −1 dw (site 1), with an average value of 659.1 ng g −1 dw. For 2-ring PAHs (NaP), the levels ranged from 5.3 to 73.8 ng g −1 dw with an average value of 24.1 ng g −1 dw; for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 42.9 to 186.3 ng g −1 dw; for 4-ring PAHs (Fla, Pyr, BaA, Chr), from 61.7 ng g −1 dw to 199.7 ng g −1 dw; for 5-ring PAHs (BbF, BkF, BaP, DahA), from 262.7 to 507.1 ng g −1 dw; and for 6-ring PAHs (BghiP, IcdP), from 17.5 to 133.2 ng g −1 dw. Figure 9 shows the individual and total concentrations of PAHs found in the sediment samples taken at the different sampling sites near the mouth of the Volturno River. The figure shows that the highest concentrations of PAHs were found at the mouth of the river (site 1) and 500 m from the mouth in the southerly direction (site 4). The composition profile of PAHs in sediment is shown in Figure 10. Five-ring PAHs were found in most of the test sites at

PAH Concentrations in Sediment from the Volturno River
Data on individual PAH concentrations found in the Volturno River are given in Table  S3, while total concentrations were previously reported [40]. In detail, total concentrations were between 434.8 ng g −1 dw (site 8) and 872.1 ng g −1 dw (site 1), with an average value of 659.1 ng g −1 dw. For 2-ring PAHs (NaP), the levels ranged from 5.3 to 73.8 ng g −1 dw with an average value of 24.1 ng g −1 dw; for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 42.9 to 186.3 ng g −1 dw; for 4-ring PAHs (Fla, Pyr, BaA, Chr), from 61.7 ng g −1 dw to 199.7 ng g −1 dw; for 5-ring PAHs (BbF, BkF, BaP, DahA), from 262.7 to 507.1 ng g −1 dw; and for 6-ring PAHs (BghiP, IcdP), from 17.5 to 133.2 ng g −1 dw. Figure 9 shows the individual and total concentrations of PAHs found in the sediment samples taken at the different sampling sites near the mouth of the Volturno River. The figure shows that the highest concentrations of PAHs were found at the mouth of the river (site 1) and 500 m from the mouth in the southerly direction (site 4). The composition profile of PAHs in sediment is shown in Figure 10. Five-ring PAHs were found in most of the test sites at 57,4% of the total PAHs in the sediment.
value of 659.1 ng g −1 dw. For 2-ring PAHs (NaP), the levels ranged from 5.3 to 73.8 ng g −1 dw with an average value of 24.1 ng g −1 dw; for 3-ring PAHs (Acy, Ace, Flu, Phe, Ant), from 42.9 to 186.3 ng g −1 dw; for 4-ring PAHs (Fla, Pyr, BaA, Chr), from 61.7 ng g −1 dw to 199.7 ng g −1 dw; for 5-ring PAHs (BbF, BkF, BaP, DahA), from 262.7 to 507.1 ng g −1 dw; and for 6-ring PAHs (BghiP, IcdP), from 17.5 to 133.2 ng g −1 dw. Figure 9 shows the individual and total concentrations of PAHs found in the sediment samples taken at the different sampling sites near the mouth of the Volturno River. The figure shows that the highest concentrations of PAHs were found at the mouth of the river (site 1) and 500 m from the mouth in the southerly direction (site 4). The composition profile of PAHs in sediment is shown in Figure 10. Five-ring PAHs were found in most of the test sites at 57,4% of the total PAHs in the sediment.

Evaluation of the Carcinogenic Risk for Human Health from Dermal and Accidental Ingestion Exposure to the PAHs Present in the Sediments of the Surface Waters
To assess the Carcinogenic Risk from exposure to PAHs present in estuarine sediments of the Sele, Volturno, and Sarno Rivers, doses of contaminants taken up by human receptors through the different routes of exposure were evaluated, and the results obtained are reported in Table 7. Particularly, for all three rivers, the doses of each individual PAH and the total doses relating to the entire class of compounds taken up by human receptors through dermal and oral exposure were calculated. The doses from accidental ingestion of PAHs from sediments (Doseingestion) were 7.86 × 10 −4 mg Kg −1 /day for the Sele River, 8.13 × 10 −4 mg Kg −1 /day for the Volturno River, and 3.31 × 10 −4 mg Kg −1 /day for the Sarno River. On the other hand, in relation to the doses of PAHs taken up by skin contact with sediment (Dosedermal), the results were 3.23 × 10 −5 , 1.36 × 10 −5 , and 3.35 × 10 −5 mg Kg −1 /day for the Sele, Sarno, and Volturno Rivers, respectively.
Moreover, to assess the carcinogenic and mutagenic potencies of PAHs in relation to

Evaluation of the Carcinogenic Risk for Human Health from Dermal and Accidental Ingestion Exposure to the PAHs Present in the Sediments of the Surface Waters
To assess the Carcinogenic Risk from exposure to PAHs present in estuarine sediments of the Sele, Volturno, and Sarno Rivers, doses of contaminants taken up by human receptors through the different routes of exposure were evaluated, and the results obtained are reported in Table 7. Particularly, for all three rivers, the doses of each individual PAH and the total doses relating to the entire class of compounds taken up by human receptors through dermal and oral exposure were calculated. The doses from accidental ingestion of PAHs from sediments (Dose ingestion ) were 7.86 × 10 −4 mg Kg −1 /day for the Sele River, 8.13 × 10 −4 mg Kg −1 /day for the Volturno River, and 3.31 × 10 −4 mg Kg −1 /day for the Sarno River. On the other hand, in relation to the doses of PAHs taken up by skin contact with sediment (Dose dermal ), the results were 3.23 × 10 −5 , 1.36 × 10 −5 , and 3.35 × 10 −5 mg Kg −1 /day for the Sele, Sarno, and Volturno Rivers, respectively. Moreover, to assess the carcinogenic and mutagenic potencies of PAHs in relation to BaP, the most carcinogenic compound among the PAH considered [10], the average concentrations of individual analytes were converted to the corresponding toxic equivalent concentrations (BaP eq ). The equivalent toxic concentrations (BaP eq ) obtained for the Sele, Volturno, and Sarno Rivers are given in Table 7.
Furthermore, the incremental risk of developing lifelong cancer expressed as ILCR was assessed for the exposure to PAHs by ingestion (ILCR ingestion ) and dermal contact (ILCR dermal ) [65]. The ILCR ingestion and ILCR dermal values obtained for Sele, Volturno, and Sarno Rivers are shown in Table 8. The values obtained for carcinogenic risk due to exposure to PAHs by ingestion (ILCR ingestion ) and dermal contact (ILCR dermal ) were found to be comparable for the three rivers. According to the USEPA, the ILCR values were interpreted by reference to three ranges, each of which is associated with a risk of carcinogenicity: an ILCR value < 1 × 10 −6 is associated with a low or zero carcinogenic risk; ILCR values between 1 × 10 −4 and 1 × 10 −6 are indicators of a moderate carcinogenic risk; and an ILCR value higher than 1 × 10 −4 corresponds to a high carcinogenic risk associated with exposure to PAHs in sediment [66,67]. The ILCR ingestion values obtained for the three rivers were found to be much lower than the ILCR dermal , indicating that the risk of cancer associated with dermal contact exposure to PAHs present in estuarine sediments may be higher than that associated with accidental ingestion exposure. Specifically, the ILCR ingestion and ILCR dermal values obtained for all sediment samples taken near the mouth of the Sele River ranged between 6.69 × 10 −15 and 3.11 × 10 −13 and 7.27 × 10 −8 and 3.38 × 10 −6 , respectively. The ILCR ingestion and ILCR dermal values obtained for the Sarno River ranged between 2.45 × 10 −15 and 1.31 × 10 −13 and 2.66 × 10 −8 and 1.42 × 10 −6 , respectively. For the Volturno River, the ILCR ingestion values were between 8.23 × 10 −15 and 3.22 × 10 −13 , while those of ILCR dermal were between 8.95 × 10 −8 and 3.50 × 10 −6 . Thus, the risk of cancer associated with exposure to PAHs by ingestion of estuarine sediments of the Sele, Sarno, and Volturno Rivers was found to be low at all sampling sites. However, based on the ILCR dermal values obtained, the risk of cancer associated with exposure by dermal contact with the PAHs present in the sediments was found to be moderate (average ILCR dermal for the three rivers of 2.77 × 10 −6 ). In addition, the values of the ILCR ingestion and ILCR dermal indices for sediment samples taken at sites with the highest concentrations of PAHs were evaluated. For the Sele and Volturno Rivers, the assessment was carried out at the mouth (site 1), for which total PAH concentrations of 0.8449 and 0.8721 mg Kg −1 dw were found, respectively. For the Sarno River, the assessment was carried out at the sampling site 150 m to the west of the estuary (site 9), where a total concentration of PAHs of 0.6792 mg Kg −1 dw was found. The results obtained for the three rivers are shown in Table 9. The ILCR ingestion values obtained for the three rivers at the sampling sites with the highest PAH concentrations were <1 × 10 −6 (order of 10 −13 ), suggesting that the risk of cancer associated with exposure by ingestion of PAHs present in estuarine sediments was low or zero. On the other hand, the ILCR dermal values obtained for the three rivers at the sampling sites with the highest PAH concentrations were in the order of 10 −6 , suggesting that the risk of cancer associated with dermal contact exposure to PAHs present in estuarine sediments was moderate. In fact, as stated also by Cheng et al., ILCR values between 1 × 10 −4 and 1 × 10 −6 are associated with a moderate carcinogenic risk [66,67]. Thus, on the basis of the ILCR values obtained by taking into account total PAH concentrations at all sites or considering only the total concentrations recorded at the most polluted sites, the risk of cancer associated with exposure by ingestion was found to be low or zero, but the risk associated with dermal contact exposure of PAHs present in estuarine sediments of the Sele, Volturno, and Sarno Rivers is moderate. Thus, since these areas were previously considered potentially contaminated according to Italian environmental law (D. Lgs. 152/2006), and as stated by Albanese et al., who assessed an incremental lifetime cancer risk higher than 1 × 10 −5 for the city of Naples [9], a continuous monitoring of potentially hazardous substances is necessary to ensure the protection of public health.

Conclusions
This paper presents for the first time an assessment of the carcinogenic risk to human health from dermal and ingestion exposure to PAHs present in sediments of the main surface water streams of the Campania Region, southern Italy. The paper also provides information on the concentrations, spatial distribution, and composition profiles of the PAHs detected in sediments collected near the Sele, Sarno, and Volturno River estuaries. The results obtained indicate that the risk of cancer following oral exposure to PAHs in estuarine sediments, expressed as incremental lifetime cancer risk (ILCR ingestion ), is low, unlike the risk from accidental skin exposure, which was moderate with ILCR dermal values between 1 × 10 −4 and 1 × 10 −6 . This calls for ongoing assessment of the carcinogenic risk to human health posed by cutaneous and oral exposure to PAHs, as well as constant monitoring of PAH concentrations in surface water sediments in the Campania Region. In conclusion, this study represents a starting point for future studies aimed at assessing the risk of carcinogenicity to human health due to exposure to the PAHs in order to provide support for pollution prevention measures and ecological restoration strategies for rivers, as well as for the preservation of the general well-being.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/toxics11020172/s1, Table S1: PAH levels (ng g-1 dw) with SD (Standard Deviation) detected in sediment samples from Sele River; Table S2. PAH concentrations (ng g-1 dw) with SD (Standard Deviation) found in sediment samples from the Sarno River (April 2008). Table S3. PAH levels (ng g-1 dry weight) with SD (Standard Deviation) detected in sediment samples from Volturno River.