Assessing Heavy Metals in the Sele River Estuary: An Overview of Pollution Indices in Southern Italy

Rapid industrialization, coupled with a historical lack of understanding in toxicology, has led in an increase in estuary pollution, frequently resulting in unexpected environmental situations. Therefore, the occurrence of heavy metals (HMs) constitutes a major environmental issue, posing a serious risk both to aquatic ecosystems and public health. This study aimed to evaluate the levels of eight HMs (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn) in water, suspended particles, and sediment near the Sele River estuary (Italy) in order to assess their environmental impacts on the sea and health risks for humans. The results revealed an increasing order of HM concentration according to the scheme suspended particulate matter (SPM) > sediment (SED) > dissolved phase (DP) and a moderate contamination status in sediment. The health risk assessment indicated that the non-carcinogenic risk was negligible. Carcinogenic risk, expressed as the incremental lifetime cancer risk (ILCR), was negligible for Cd and Ni and within tolerable limits for As, Pb, and Cr. The findings suggested that, even if there are currently no specific limits for chemical parameters in the transitional waters of Italy, monitoring systems should be implemented to determine pollution levels and implement effective steps to improve river water quality and reduce human health risks.


Introduction
Human communities have always been clustered around rivers and estuaries for easy access to freshwater and transportation [1].Moreover, rapid urbanization and industrialization have caused an increase in pollution in coastal areas, especially along the lotic systems that run from rivers and estuaries to offshore coastal regions [2].Estuaries are positioned at the confluence of land and sea, and they are impacted by both, making their environmental conditions unpredictable [3,4].The transitional nature of estuaries presents challenges in providing a precise definition, yet they are commonly acknowledged as coastal water bodies strongly affected by tides.Predominantly enclosed by land, estuaries are characterized by a noticeable dilution of seawater caused by freshwater inputs from rivers and runoff [5].Consequently, the wide array of physical-chemical, climatic, and morphological factors contribute to a significant level of variability.Coupled with pronounced spatial heterogeneity and complexity, each estuarine area emerges as a unique environment with characteristics that defy straightforward generalization and classification [6].Thus, industrialization has resulted in an increase in the amount of waste dumped in estuaries and river systems, posing issues for both ecosystem and environmental management [7].
Rivers provide a variety of ecological and human services, including water transportation, ecotourism, aquaculture, ecological habitats, and ecological defense [8][9][10].However, they have been continuously threatened by various contaminants, among which some of the world's most hazardous contaminants found in riverine aquatic environments are heavy

Study Area
The Sele River is located in the Campania Region, in southern Italy (Figure 1).With a length of 64 km, it is the second most important river in the region in terms of average water volume, with a drainage basin of 3235 km 2 and an annual mean flow rate of 69 m 3 /s [31,32].The river rises in the province of Avellino near the Caposele, and receives different tributaries, including the Tanagro, the Calore and the Lucano.Then, the river enters the floodplain known as the Sele River plain [34].

Sampling
Water samples were collected in four different seasons to assess the spatio-temporal trends in HM concentrations near the Sele River estuary.The sampling was carried out during 2020-2021 in summer (July), autumn (November), winter (February), and spring (April) [33,35,36].For sediment, samples were collected only in one season (April).The sampling locations are indicated in Figure 1 and the characteristics of each sampling site are detailed in Table 1.Samples were collected at 10 stations, from the river mouth to different distances from it and, in detail, north, west, and south from the estuary.Three water samples aliquots for each sampling location (500 mL) were collected using polypropylene bottles, transported to the laboratory, refrigerated, and stored at +4 °C until analysis, which was carried out within 6 h of sampling.From each sampling site, three sediment samples, each weighing approximately 5 g, were collected at a depth of 0-5 cm using a Van Veen grab sampler.These samples were placed in polyethylene bags and promptly transported to the laboratory within 1 h of collection.Upon arrival, they were cooled at −20 °C until analysis, which was conducted within 1 day of sampling.

Sampling
Water samples were collected in four different seasons to assess the spatio-temporal trends in HM concentrations near the Sele River estuary.The sampling was carried out during 2020-2021 in summer (July), autumn (November), winter (February), and spring (April) [33,35,36].For sediment, samples were collected only in one season (April).The sampling locations are indicated in Figure 1 and the characteristics of each sampling site are detailed in Table 1.Samples were collected at 10 stations, from the river mouth to different distances from it and, in detail, north, west, and south from the estuary.Three water samples aliquots for each sampling location (500 mL) were collected using polypropylene bottles, transported to the laboratory, refrigerated, and stored at +4 • C until analysis, which was carried out within 6 h of sampling.From each sampling site, three sediment samples, each weighing approximately 5 g, were collected at a depth of 0-5 cm using a Van Veen grab sampler.These samples were placed in polyethylene bags and promptly transported to the laboratory within 1 h of collection.Upon arrival, they were cooled at −20 • C until analysis, which was conducted within 1 day of sampling.

Metal Extraction and Analysis
The methodology used for preparation and extraction was previously described [35,37].Briefly, water samples were filtered using GF/F glass fiber filters (47 mm × 0.7 µm; Whatman, Maidstone, UK) in order to assess the levels of HMs in the dissolved phase (DP), while filters were used for suspended particulate matter (SPM).Then, samples were acidified with 1% nitric acid (HNO 3 )/chloridric acid (HCl) mixture (v/v).For sediment, the method consisted of digesting aliquots of 5 g, prior to which samples were air-dried and homogenized [38][39][40][41].
Then, the digestion of both sediment and SPM samples was carried out using HNO 3 / HCl (3:1 v/v) mixture solution with a microwave digestion system (Milestone 1200).Subsequently, the samples were diluted to 25 mL with deionized water and filtered.The analysis was performed by inductively coupled plasma mass spectrometry (ICP-MS) (Thermo Scientific TM ICAP TM RQ) according to method 6020 B proposed by the US EPA (2014) [42], and the amounts of eight HMs (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn) were evaluated.The tune solution (Ba, Bi, Ce, Co, In, Li, and U, each 1.00 µg/L; (Thermo Scientific, Bremen, Germany)) was used for optimization of the instrumental operating conditions and the kinetic energy discrimination (KED) mode with helium collision gas was applied.Deionized water (18.2MΩ), HNO 3 (65% m/m), and HCl (37% m/m) were used.The concentrations of the studied compounds were estimated by the construction of calibration curves (CertiPUR ® , Merck, Darmstadt, Germany).To check the method, all samples were analyzed in duplicate and each sample was measured in triplicate by Q-ICP-MS detection.

Quality Control and Quality Assurance
In order to evaluate quality assurance and quality control, duplicates, method blanks, and common reference materials were used.To reduce the background contamination, all glassware, plastics, and quartz were cleaned in 10% HNO 3 before being washed in deionized water.Triplicate samples and spiked samples were both used to assess the precision of the instrumental procedures.The limits of detection and quantification (LOD and LOQ, respectively) were evaluated in 3 and 10 times the signal/noise ratio for each analyte, respectively.In detail, for water, the LODs were 0.20 µg/L for As, 0.07 µg/L for Hg, 0.15 µg/L for Cd, 0.19 µg/L for Cr, 0.47 µg/L for Cu, 0.52 µg/L for Ni, 0.47 µg/L for Pb, and 0.65 µg/L for Zn.For sediment, the LODs were 0.32 mg/kg dw for As, 0.21 mg/kg dw for Hg, 0.08 mg/kg dw for Cd, 0.29 mg/kg dw for Cr, 0.74 mg/kg dw for Cu, 0.80 mg/kg dw for Ni, 0.56 mg/kg dw for Pb, and 1.04 mg/kg dw for Zn.The LOQs were 0.67 µg/L for As, 0.23 µg/L for Hg, 0.50 µg/L for Cd, 0.63 µg/L for Cr, 1.57 µg/L for Cu, 1.73 µg/L for Ni, 1.57 µg/L for Pb, and 2.17 µg/L for Zn for water.The LOQs for sediment were 1.07 mg/Kg dw for As, 0.70 mg/kg dw for Hg, 0.25 mg/kg dw for Cd, 0.97 mg/kg dw for Cr, 2.47 mg/kg dw for Cu, 2.67 mg/kg dw for Ni, 1.87 mg/kg dw for Pb, and 3.47 mg/kg dw for Zn.For quantification, five-point calibration curves were prepared for each analyte (R 2 > 0.990).The quality control guidelines of the US EPA were closely followed during each analysis.Routine instrument calibrations were conducted for every analyzed batch or more frequently as needed.Periodic checks were conducted on the method detection limits.For each batch of analysis, reagent blanks, duplicate samples, and sample spike recoveries were all examined for every 10 samples.Certified criteria were established for spiked recoveries (80-120%) and the relative percent difference of duplicate samples (10%).Heavy metal recoveries in spiked samples were in the range of 88-113%, which complied with quality control standards.
For sediment samples, certified reference material (SD-MEDPOL-1/TM of the International Laboratory of Marine Radioactivity, IAEA, Monaco, Germany) was run after each batch of samples to test the accuracy and precision.Recoveries obtained for sedimentcertified reference material samples were in the range of 90-105%.

Statistical Analysis and Heavy Metal Input
Data reported by the study were analyzed using STATA MP v14.0 statistical software (College Station, TX, USA).Recognizing pollution origins and allocating their contribution play an essential role in heavy metal contamination analysis.Techniques like principal component analysis (PCA) are employed to pinpoint sources of heavy metal pollution [35].PCA, a linear method aimed at reducing complexity, transforms various environmental factors into distinct sources utilizing only a handful of principal components.Its application yields significant insights that encapsulate the entire dataset without compromising the original information.Typically, PCA is employed to uncover relationships among heavy metal variables, aiding in identifying the sources within water surfaces and sediment.Principal component analysis (PCA) was executed.The requirements were previously described [35].Briefly, the first components with eigenvalues > 1 were considered for the amount of explained variance, and the components had to explain the 70-80% of the variance.A p-value < 0.05 was considered to be the significance level.For lower LOQ levels, an LOQ half-value was used.The flow-averaged mean concentration (C aw ) and annual contaminant loads (F annual ) were obtained according to UNEP guidelines [43], based on Equations ( 1) and ( 2), respectively: where Q T represents the total river input for the study period, assessed using the sum of the monthly averaged water flow; C i is the instantaneous amount; and Q i is the load, expressed as the daily water flow average in kg/year [33,35].In detail, the cumulative loads of heavy metals (HMs) introduced into the Mediterranean Sea from the Sele River were computed by taking into account the concentration values of individual contaminants at the river mouth during the four months of sampling.Subsequently, the average of the total concentrations was multiplied by the annual mean flow rate (m 3 /year) of the Sele River.

Heavy Metal Pollution Indices in Sediment
The contamination level of HMs in sediment from the Sele River estuary was assessed using the geoaccumulation index (I geo ), the contamination factor (CF), the pollution index (PLI), and the potential ecological risk index (PERI).
I geo was used to estimate the contamination status of the environment by comparison against geochemical background concentrations, which was evaluated according to Equation (3) [35,44]: where C n and B n are the amount detected in the sediment sample and the geochemical background concentration for an individual HM (n), respectively.Specifically, the background values considered for each heavy metal under investigation were: 5.00 mg/kg dw (As), 0.07 mg/kg dw (Hg), 0.05 mg/kg dw (Cd), 90.00 mg/kg dw (Cr), 45.00 mg/kg dw (Cu), 68.00 mg/kg dw (Ni), 20.00 mg/kg dw (Pb), and 95.00 mg/kg dw (Zn) [45,46].A factor of 1.5 indicated the background matrix correlation value, which represented the possible variation due to the lithogenic effect and weathering [47].As a qualitative scale of pollution intensity, I geo includes seven classes on the basis of the degree of contamination.Therefore, based on the I geo value, a site can be classified from unpolluted to extremely polluted, as showed in Table 2.Moreover, the status of HM contamination in sediment was further assessed using the contamination factor (CF), contamination degree (CD), and pollution load index (PLI), which resulted in the cumulative indication of the overall level of heavy metal pollution in sediment from the research area.Therefore, the CF, CD and PLI indices were evaluated using Equations ( 4)-( 6) below [44,[48][49][50]: where n represents the number of HMs evaluated.The potential ecological risk index (PERI), applied to estimate the potential ecological risk of the studied HMs in sediment samples collected in the study area, allows an ecological risk assessment to be carried out by combining ecological and environmental effects with toxicological data, and it was calculated using Equations ( 7) and (8) [35,51]: where T r i indicates the toxic response factor of individual HMs, as shown in Table 3.

Ecological Risk Index for Sediment
The Sediment Quality Guidelines (SQGs) were used to estimate the potential risk for aquatic organisms and water systems associated with HMs in sediment [53].In detail, threshold effect level (TEL), probable effect level (PEL), effect range low (ERL), and effect range median (ERM) values were considered [54][55][56].These values are indicated in Table 4. Therefore, when amounts were lower than the ERL and TEL values, no adverse effects would be observed.In addition, HM concentrations were considered rarely to be toxic for values below the ERL and above the TEL.In contrast, when levels were above the ERM and PEL values, adverse effects on the ecosystem could be predicted [50,57].

National Recommended Water Quality Criteria
The toxicity of HMs in water samples from the Sele River estuary was assessed according to the National Recommended Water Quality Criteria by the US EPA [58].Therefore, when the HM amount in the water phase exceeds its corresponding criterion continuous concentration (CCC) value, chronic toxicity may occur in the aquatic ecosystem.In contrast, acute toxicity may occur when the HM concentration detected in the water phase exceeds its corresponding criterion maximum concentration (CMC) value [59].

Health Risk Assessment: Carcinogenic and Non-Carcinogenic Risks
A health risk assessment estimates the probability of any given amount of adverse health effects occurring during a certain time period [60,61].A contaminant's health risk is generally based on an estimation of the risk level and is categorized as carcinogenic or non-carcinogenic health risk [62,63].
Toxics 2024, 12, 38 According to Mohammadi et al., Table 5 shows the parameters used for risk assessment [61].Also, the hazard quotient (HQ) and hazard index (HI) were used to assess the potential non-carcinogenic health risk induced by heavy metal levels in water via ingestion and dermal absorption [61].The USEPA proposed guidelines for assessing the non-carcinogenic and carcinogenic risks via ingestion and dermal absorption pathways [64,67].Therefore, the evaluation was conducted according to Equations (10) and (11).To estimate the HQ, the reference oral doses (RfDs) through oral ingestion and dermal absorption routes were used, as indicated in Table 6.When HQ > 1, adverse non-carcinogenic effects on health may occur, while for HQ values below 1, no risk is predicted [60,64].The HI, expressed as the sum of individual HQs, is used to evaluate the total noncarcinogenic risk due to exposure to a mixture of potentially toxic chemicals in water systems [70,71].When HI > 1, non-carcinogenic dangers are predicted; however, HI < 1 denotes that exposed people are unlikely to have obvious negative health effects [60,61].
HQ =HQ As + HQ Hg + HQ Cd + HQ Cr + HQ Cu + HQ Ni + HQ Pb + HQ Zn (12) Toxics 2024, 12, 38 9 of 26 Furthermore, the incremental lifetime cancer risk (ILCR) was used to assess the probable cancer risks caused by total exposure (ingestion and dermal contact) to HM concentrations in water.It was estimated according to Equation (12).
where CSF represents the cancer slope factor (CSF), expressed as kg/day/mg, which correlates with the average exposure concentration leading to an increase in the probability of developing cancer.The CSFs applied in this study are displayed in Table 7.

Levels of Heavy Metals in the Water Dissolved Phase (DP)
The concentrations of HMs found in the dissolved phase (DP) ranged from nd to 14.28 µg/L (with a mean value of 1.74 ± 2.88 µg/L) for As, from nd to 2.35 µg/L (with an average of 0.27 ± 0.53 µg/L) for Hg, from nd to 1.69 µg/L (average of 0.28 ± 0.43 µg/L) for Cd, from nd to 14.97 µg/L (average of 2.63 ± 3.77 µg/L) for Cr, from nd to 13.23 µg/L (mean value of 1.34 ± 2.56 µg/L) for Cu, from nd to 11.15 µg/L (average value of 1.56 ± 2.53 µg/L) for Ni, from nd to 14.41 µg/L (mean value of 0.88 ± 2.50 µg/L) for Pb, and from nd to 5.46 µg/L (average of 1.02 ± 1.45 µg/L) for Zn.The results obtained for water dissolved phase (DP) samples collected at the sampling stations near the river during the research period (2020-2021) are shown in Table S1.

Levels of Heavy Metals in the Suspended Particulate Matter (SPM)
The concentrations of HMs detected in the suspended particulate matter (SPM) samples collected from the sampling locations near the Sele River estuary ranged from nd to 116.52 µg/L (with a mean value of 17.77 ± 24.82 µg/L) for As, from nd to 10.03 µg/L (with an average of 1.77 ± 2.00 µg/L) for Hg, from nd to 4.57 µg/L (average of 0.50 ± 0.86 µg/L) for Cd, from nd to 33.47 µg/L (average of 4.94 ± 6.28 µg/L) for Cr, from nd to 139.85 S2.

Levels of Heavy Metals in the Sediment
The levels of individual and total HMs detected in sediment samples collected at 10 sampling locations near the Sele River estuary during four seasons (2020-2021) are illustrated in Table S3.Regarding individual HMs, the levels found were in the range of nd to 12.56 mg/kg dw(with a mean value of 3.68 ± 3.60 mg/kg dw) for As and nd-1.10 mg/kg dw(with a mean value of 0.27 ± 0.44 mg/kg dw) for Hg.Moreover, the concentrations ranged from nd to 0.98 mg/kg dw (mean value of 0.19 ± 0.40 mg/kg dw), from 1.22 to 29.10 mg/kg dw (mean value of 9.22 ± 8.84 mg/kg dw), and from nd to 7.32 mg/kg dw (mean value of 2.22 ± 2.34 mg/kg dw) for Cd, Cr, and Cu, respectively.Additionally, the mean value detected for Ni was 4.26 ± 7.80 mg/kg dw (range of nd-24.89mg/kg dw) and 6.72 ± 10.31 mg/kg dw (range of nd -34.57mg/kg dw) for Pb.Finally, the levels found for Zn were in the range of 5.92-42.00mg/kg dw (with a mean value of 17.62 ± 12.29 mg/kg dw).

Heavy Metal Inputs and Spatio-Temporal Distribution
The total annual HM load from the Sele River estuary into the Tyrrhenian Sea (central Mediterranean Sea) was 483.97 kg/year.Specifically, regarding the individual contribution of each HM, the load was 81.69 kg/year for As, 8.98 kg/year for Hg, 3.37 kg/year for Cd, 32.95 kg/year for Cr, 111.02 kg/year for Cu, 60.25 kg/year for Ni, 112.81 kg/year for Pb, and 72.90 kg/year for Zn.
The spatial distributions of HMs for DP, SPM, and SED samples assessed at different sampling locations during four seasons are shown in Figures 2-4, respectively.
Figure 2 displays the spatio-temporal distribution of HMs detected in DP samples collected during the research period.Higher levels were detected at the mouth and 500 m to the south (500S) during the warmer and drier seasons (July > April).In detail, the highest concentrations were recorded at the mouth for Hg (2.35 µg/L), Cr (14.97 µg/L), Ni (11.15 µg/L), and Pb (14.41 µg/L) in the dry season (July).For As and Cu, the highest amounts were found in the same period, with maximum values of 14.28 µg/L and 13.23 µg/L, respectively.In contrast, lower levels of HMs were detected during the winter season (February < November).Moreover, the results showed that the HM pollution degree (expressed as the total concentration) was greater at the mouth (with an average amount over four seasons of 26.28 ± 23.99 µg/L) and 500S (with an average amount over four seasons of 27.02 ± 23.17 µg/L), which decreased considerably in the northern direction, with a minimum value recorded at the 1500N location (0.34 ± 0.68 µg/L).

Heavy Metal Inputs and Spatio-Temporal Distribution
The total annual HM load from the Sele River estuary into the Tyrrhenian Sea (central Mediterranean Sea) was 483.97 kg/year.Specifically, regarding the individual contribution of each HM, the load was 81.69 kg/year for As, 8.98 kg/year for Hg, 3.          Figure 3 illustrates the spatio-temporal distribution of HMs in SPM samples from the Sele River estuary.A quick glance reveals that a higher amount was found at 500S (500 m to the south) during the hot season (July), with a maximum value of 582.64 µg/L (53,948.15mg/kg dw), expressed as the total amount recorded.Regarding individual analytes, the highest concentrations detected during this sampling period at the 500S site were 116.52 µg/L for As, 10.03 µg/L for Hg, 4.57 µg/L Cd, 33.47 µg/L for Cr, 139.85 µg/L for Cu, 88.74 µg/L for Ni, and 107.54 µg/L for Zn.Only for Pb was the greatest amount recorded at the mouth (86.66 µg/L).Therefore, a similar pattern to that observed for DP was also observed in this instance, with a higher degree of HM pollution found in the summer season (July) and near the mouth (500S), and with a decreasing trend moving in the northern direction.Particularly, regarding spatial distribution, the highest total level of HMs was found at 500S (271.49± 231.94 µg/L, expressed as the average value over four seasons).
Figure 4 shows the spatio-temporal distribution of HMs in sediment samples near the Sele River estuary.Similar to the DP and SPM phases, the highest concentration in the sediment phase, amounting to 147.65 mg/kg dw, expressed as the total amount recorded, was also detected at 500S (500 m to the south) during the sampling season (April), followed by the mouth (88.19 mg/kg dw).Considering the amounts detected for each analyte, the highest concentrations were found at the 500S location for As (12.56 mg/kg dw), Hg (1.10 mg/kg dw), Cr (29.10 mg/kg dw), Ni (24.89 mg/kg dw), Pb (34.57mg/kg dw), and Zn (42.00 mg/kg dw).However, for Cd and Cu, the highest levels were obtained at the mouth, corresponding to values of 0.28 and 7.32 mg/kg dw, respectively.Therefore, an increasing trend moving in the southern direction was also observed for sediment.
Figure 5 illustrates the comparison between HM amounts detected in the DP and SPM phases.The figure suggests that higher levels were found in SPM samples.In addition, the HM concentration values found in the SPM phase were significantly higher than those found in the sediment phase, expressed as the total amount recorded in the same period of collection (April) (Figure 6).In fact, in contrast to the much lower concentration found in sediment (147.65 mg/kg dw), the highest total amount was found in SPM (33,102.1765 mg/kg dw), recorded at the same site (500S) and during the same sampling period (April).Accordingly, the partition coefficients (K p ), assessed as the ratio between HM amounts detected in the SED and SPM samples, were evaluated.Also, the ratios of HM concentrations obtained in the SPM samples and those in the DP samples were calculated.Particularly, the results of K p indicated an increasing trend in HM amount from DP to SED and more to SPM.Therefore, the findings revealed an increasing order of HM concentration according the scheme SPM > SED > DP.
Data were also reported as the sum of DP + SPM.Specifically, the total HM concentrations found in the water matrix (as the sum of both phases) were in the range of nd to 641.29 µg/L, with a mean value of 105.78 ± 138.47 µg/L (evaluated for all sampling campaigns and sites).Moreover, the highest amount was found at location 8 (mean value of 298.51 ± 255.05 µg/L over four seasons).Moreover, the total HM levels found in the water matrix (sum of DP and SPM) for each analyte in detail were in the range of nd to 130. 80  Compared to previous studies, the HM levels found in DP samples collected around the Sele River estuary were lower than those evaluated near the Sarno River estuary and Tiber River estuary by Montuori et al., who reported that total HM concentrations were in range of 0.32-1680.39µg/L and 11.80-1359.67µg/L, respectively [26,72].Instead, a recent study carried out in an area close to that of this study, near the estuary of the Volturno River, in southern Italy, recorded total HM concentrations similar to those obtained in this study (range of 0.20-65.67µg/L) [35].For the suspended phase, again, the amounts reported by Montuori et al. for the Sarno and Tiber River estuaries were higher (ranges of 103.6-7734.6 µg/L and 18.74-2841.56µg/L, respectively) [26,72] than those found in this study, although they were comparable to those detected by De Rosa et al. near the Volturno River estuary (1.37-601.77µg/L) [35].Compared to previous research conducted in southern Italy on HM amounts in sediment samples, the results obtained in this study showed lower values than those obtained for sediments collected near the estuaries of the Sarno and Tiber Rivers (in the range of 90.7-2470.3mg/Kg dw and 42.8-1686.8mg/kg dw, respectively), but similar to those obtained around the Volturno River estuary (14.62-157.33 mg/kg dw) [26,35,72].Data were also reported as the sum of DP + SPM.Specifically, the total HM concentrations found in the water matrix (as the sum of both phases) were in the range of nd to 641.29 µg/L, with a mean value of 105.78 ± 138.47 µg/L (evaluated for all sampling campaigns and sites).Moreover, the highest amount was found at location 8 (mean value of 298.51 ± 255.05 µg/L over four seasons).Moreover, the total HM levels found in the water   Data were also reported as the sum of DP + SPM.Specifically, the total HM concentrations found in the water matrix (as the sum of both phases) were in the range of nd to 641.29 µg/L, with a mean value of 105.78 ± 138.47 µg/L (evaluated for all sampling campaigns and sites).Moreover, the highest amount was found at location 8 (mean value of 298.51 ± 255.05 µg/L over four seasons).Moreover, the total HM levels found in the water

Statistical Analysis
For statistical analysis, PCA was employed for the three phases (dissolved phase, suspended particulate matter, and sediment).Particularly, the findings of the study indicated that the HM contamination level was higher for locations sampled south of the mouth.Furthermore, a similar situation was observed for the sediment (SED) and dissolved phase (DP) data, indicating that only the greatest eigenvalue was greater than 1, and the first component contributed to 83.78% and 80.15% of the total variance, respectively (Figures S1  and S2).Furthermore, all variables contributed to the determination of the first principal component with comparable impact, indicating that the first component was an average of each metal (variable), since the presence of all studied HMs was observed at each site.For suspended particulate matter (SPM), in addition to the eight studied HMs, the weight of filter variable was considered.The results showed that this variable was orthogonal to the others, which were high correlated with one another, so there were both the components that explained the total variance of the variables (HMs) and the component that explained the variance of the weight of filter variable.The two highest eigenvalues were larger than 1 and the two first components explained 96.15% of the total variance (Figure S3).

Heavy Metal Contamination Status
The contamination status caused by HMs in the study area was estimated using I geo .Table 8 and Figures 7 and 8 show the I geo values evaluated in SED samples from the Sele River estuary, in southern Italy.Particularly, higher values of I geo max were found for As and Hg at site 8 (500S) and for Cd at site 1 (the mouth).On the other hand, lower values were recorded for Cr and Ni at locations 4 (1500N) and 2 (500N), respectively.Moreover, in order to evaluate the status of HM contamination in sediment collected in the study area, the CF, CD, PLI, and PERI indices were estimated.The results are indicated in Table 9.Moreover, in order to evaluate the status of HM contamination in sediment collected in the study area, the CF, CD, PLI, and PERI indices were estimated.The results are indicated in Table 9.

Ecological Risk Assessment
The potential ecological risk for aquatic organisms and the water system associated with the presence of HMs in sediment was estimated by comparing the concentrations detected in samples collected from the Sele River estuary with the SQGs.Therefore, the comparisons with TEL, PEL, ERL, and ERM values are displayed in Figure 10 for As, Hg, Cd, and Cr, and in Figure 11 for Cu, Ni, Pb, and Zn.

Ecological Risk Assessment
The potential ecological risk for aquatic organisms and the water system associated with the presence of HMs in sediment was estimated by comparing the concentrations detected in samples collected from the Sele River estuary with the SQGs.Therefore, the comparisons with TEL, PEL, ERL, and ERM values are displayed in Figure 10 for As, Hg, Cd, and Cr, and in Figure 11 for Cu, Ni, Pb, and Zn.

National Recommended Water Quality Criteria
Table 10 shows the comparison of HM amounts detected in DP samples collected in the study area during four seasons with the National Recommended Water Quality Criteria proposed by the US EPA [58].

National Recommended Water Quality Criteria
Table 10 shows the comparison of HM amounts detected in DP samples collected in the study area during four seasons with the National Recommended Water Quality Criteria proposed by the US EPA [58].

Health Risk Assessment
In order to estimate the potential non-cancer health risk due to ingestion and dermal absorption of HMs in water sampled near the Sele River estuary, the CDIs, HQs, and HIs were evaluated.The results obtained are indicated in Tables 11 and 12 and Figure 12.Particularly, in Table 11, minimum, maximum, and mean levels of CDIs for adults via ingestion and dermal contact routes of HMs in water samples from the study area are shown.Moreover, the minimum, maximum, and mean values of HIs evaluated in this study for HM exposure through different pathways (ingestion and dermal contact) are shown in Table 12.As given in Table 11, the mean levels of CDI dermal were higher than those of CDI ingestion , with the highest mean values recorded for Zn (9.49 × 10 −6 mg/kg/day) > Cu (1.77 × 10 −6 mg/kg/day) > Pb (1.38 × 10 −6 mg/kg/day) > As (1.29 × 10 −6 mg/kg/day).These values were used for HQs.The highest risks were estimated to be from As (1.98 × 10 −6 ) through ingestion and from Cr (5.69 × 10 −5 ) by dermal contact.
Figure 12 illustrates the graphical representation of the HQs evaluated in this study for eight HMs (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn).
Moreover, the minimum, maximum, and mean values of HIs evaluated in this study for HM exposure through different pathways (ingestion and dermal contact) are shown in Table 12.
The carcinogenic risk assessment is summarized in Table 13, which indicates the incremental lifetime cancer risk (ILCR) values obtained for total exposure caused by ingestion and dermal contact in the study area.

Discussion
The study evaluated the impact of heavy metals in the Mediterranean Sea from the Sele River, in southern Italy.Therefore, HM amounts were estimated in water (as DP and SPM) and sediment samples collected near the estuary of the river.Moreover, since only a few studies have evaluated HMs in both phases, it is difficult to compare the results obtained in this study with those of other studies.Table S4 shows a comparison of heavy metal concentrations in previous studies carried out from other river catchments and transitional waters in areas close to the study area.Therefore, compared to a few previous studies, the HM concentrations found in the water samples collected near the Sele River estuary (as the sum of DP and SPM) were lower than those detected in other Italian rivers [26,72], but comparable with those found in samples collected near the Volturno River estuary [35].
Higher amounts of HMs were observed at the river's mouth and at 500 m moving southward, possibly reflecting that HM pollution downstream is more serious than that upstream, which is typically less influenced by human activity [24].Regarding individual HMs, Cr showed the highest levels in DP, probably because Cr is not easily absorbed by suspended particles, resulting in restricted Cr input in sediment.As a result, Cr may be found mostly in bodies of water [73].In contrast, Cu was detected in the highest amount in SPM and Zn showed the highest level in the sediment phase.Accordingly, Cu and Zn, consistently found in environmentally significant amounts in sediments worldwide, mainly originate from vehicle emissions and atmospheric deposition [74].However, they may be easily accumulated in river sediment, which is the ultimate sink for HM pollution [75].Thus, these findings could be due to the fact that they are mostly collected in sediment through adsorption, sedimentation, or flocculation.Although these activities may significantly decrease the HM levels in the overlying water body, the same HMs chelated by sediments can be released into the water through the same mechanisms under specific environmental circumstances, resulting in secondary contamination [74,76].To date, there are various strategies to mitigate the presence of HMs in river water and maintain low concentrations.Accordingly, the implementation of best management practices (BMPs) in industrial and agricultural activities can significantly reduce the runoff of pollutants into rivers.Additionally, public awareness campaigns and regulatory measures play pivotal roles in encouraging responsible waste disposal practices and fostering a collective com-mitment to preserving water quality.Employing a combination of these approaches can contribute to the sustainable management of heavy metal levels in river water.The results of this study suggest that monitoring programs are required for assessing pollution levels and establishing effective measures to enhance river water quality and decrease human health threats.
The spatio-temporal distributions of HMs in DP and SPM samples collected during four seasons in 2020-2021 at 10 locations near the estuary were evaluated.For both phases, the highest levels were detected at the mouth and 500 m to the south (500S) during the warmer and drier seasons (July and April).For sediment, samples were collected only in one season (April) and the spatial distribution indicated higher levels detected 500 m to the south, followed by the mouth.The data obtained suggested that fresh inputs of these contaminants to the sea from the Sele River estuary can occur, and the SPM phase could act as carrier from the sediment to the water phase.In fact, due to their chemical characteristics, HMs can bind to the sediment phase but perturbation phenomena such as rain or earth movements can facilitate their transition to the water phase, which occurs through the suspended phase.These findings are in agreement with those of Huang et al. [75], who stated that the bioavailability of Cu and Zn in river sediment may be affected by seasonal variations and spatial distribution [75].Particularly, they assessed that seasonal variations may influence HM contamination in sediment, mainly in the dry season [75].Accordingly, recent studies carried out in Ghana and India revealed that concentrations of Cu and Zn in river sediment were higher in the dry season than in the wet season [77,78].The results of this investigation were also consistent with those by Zhao et al., who observed that seasonal changes in these pollutant levels may be caused by anthropogenic pressures and different climatic conditions.In fact, they investigated the spatio-temporal distribution patterns of HMs in water samples collected in the Yitong River (China) and found higher amounts of HMs during the summer season and in urban areas than in the winter and in suburban areas [79].On the contrary, Liu et al. stated that levels of HMs found in the surface water of Qingjiang River (China) were all higher during the wet season than in the dry season [80].
The total annual HM load from the Sele River estuary into the Mediterranean Sea was evaluated.Compared with earlier investigations conducted nearby, the total HM load from the Sele River estuary to the sea was lower than that reported by De Rosa et al., who evaluated the annual HM input at about 620.39 kg/year from the Volturno River estuary to the Mediterranean Sea [35], and drastically lower than that reported by Montuori et al. from the Sarno River estuary (13,977.53kg/year) [72].However, these findings suggest that the Sele River could be a significant point source of HM discharge into the Mediterranean Sea, even though the concentrations of the individual HMs observed, and therefore their inputs, were low.
The HM contamination status of sediment from the Sele River estuary was assessed using I geo , CF, PLI, and PERI.The average I geo values indicated that, based on the classification proposed by Hossain et al. [44], the sampling locations could be considered as unpolluted with regard to As, Cr, Cu, Ni, Pb, and Zn and unpolluted/moderately polluted for Hg and Cd.In contrast, the maximum values of I geo at sites 8 (500S) for Hg and 1 (mouth) for Cd indicated the sites were moderately and heavily polluted, respectively.Moreover, in order to evaluate the status of HM contamination in sediment collected in the study area, the CF, CD, PLI and PERI indices were estimated.According to Khadanga et al., data obtained regarding CF for Hg and Cd suggested considerable and moderate contamination, respectively [50].For CD, the results indicated a moderate contamination level for the river [48].In addition, the PLI value indicated baseline levels of HMs in sediment collected from the study area [49].According to the E r and PERI values evaluated, Hg could entail a high ecological risk for the Sele River, followed by Cd with a considerable ecological risk for the river [51].Thus, since the mean value of PERI was in the range of ≥150 to < 300, the river could be a source of HMs, mainly Hg and Cd, with considerable/moderate levels of pollution by HMs [51].Effectively, Hg, released both from natural and anthropogenic sources, mostly enters aquatic environments by riverine input.Moreover, the literature states that industrial sewage discharge and watershed and atmospheric deposition are major sources of Hg in the sediments of rivers [81,82].In addition, the development of industry and agriculture has led to a significant increase in Cd levels in sediment because of its extensive application in fertilizer and pesticide production [83].
Concerning the ecological risk for aquatic organisms and water systems associated with the presence of HMs in sediment, the amounts recorded in sediment samples collected from the Sele River estuary were compared with the SQGs.Particularly, the amount of As detected at 500S (12.56 mg/kg dw) was higher than the TEL and ERL values (Figure 11).Moreover, the amounts of Hg detected at the mouth (0.74 mg/kg dw), at 500S (1.10 mg/kg dw), and 1000S (0.81 mg/kg dw) were higher than the TEL, PEL, ERL, and ERM values.For Ni, the amount found at 500S (24.89 mg/kg dw) was above the TEL and ERL values, while the amount of Pb detected at the same site (34.57mg/kg dw) was higher than the TEL value.On the other hand, for Cd and Cr (Figure 10) and for Cu and Zn (Figure 11), no concentration value exceeded the SQGs.Therefore, according to the study results, heavy metal concentrations found near the Sele River may suggest an issue of concern.Consequently, continuous monitoring of these types of environmental contaminants is necessary to preserve the aquatic ecology of the river.
Regarding the comparison of HM amounts detected in DP samples collected in the study area during four seasons with the National Recommended Water Quality Criteria proposed by the US EPA [58], the results showed that all of the values found for As and Zn were below the CMCs and CCCs indicated for river and marine waters.On the contrary, Hg had higher concentrations than the CMC and CCC.Specifically, the levels of Hg exceeded the CMC and CCC for freshwater in 60 and 18% of the samples, respectively, and for saltwater in 6 and 19% of the samples, respectively.For Cd, the observed levels were above the CCC for freshwater in 27% of the samples, while the amount of Ni exceeded the CCC for saltwater in 13% of the samples, unlike Cr, for which 58% of the samples were over the CMC for saltwater.Moreover, the Cu amount exceeded the CMC and CCC for saltwater in 10 and 31% of the samples, respectively, while for Pb, the amount surpassed the CCC values for freshwater and saltwater in 14% and 6% of the samples, respectively.Therefore, the findings indicate that there may be a possible risk connected with the heavy metals under consideration.
Furthermore, this research computed the health risks of HMs, some of which exhibit proven toxicity to humans.Dashtizadeh et al. [84] stated that the chronic intake of HMs, such as Cd, Cr, As, and Pb, has considerable biological toxicity and is dangerous to human health.According to this study, Cd mainly accumulates in the human hepatic system and kidneys, disturbing estrogen secretion, and is also classified as carcinogenic.In addition, Cr, Cu, and Zn can cause non-carcinogenic hazards such as neurologic involvement, headache, and liver disease, while acute and chronic arsenic exposure can also cause dermal, respiratory, reproductive and carcinogenic effects [84].Moreover, Kan et al. reported that Pb exposure can cause gastrointestinal damage, child amentia, and Alzheimer's disease, while Hg exposure may result in kidney damage, central nervous system defects, arrhythmia, and respiratory problems [85].For the non-carcinogenic risk assessment, the CDI dermal values showed higher mean levels than those of CDI ingestion , with the highest mean values recorded for Zn and Cu.The values of HQ indicated that the highest risks were associated with As due to accidental ingestion and Cr caused by dermal contact.However, the HQs for all examined HMs were lower than the level of concern (HQ < 1).Furthermore, the HIs were evaluated to estimate the total potential non-carcinogenic impact, and all values were below 1 (HI < 1), suggesting that no non-carcinogenic concerns were to be expected [85,86].Accordingly, the findings of the mean HQs suggested an acceptable level of non-carcinogenic health risk in all samples taken near the Sele River estuary.Therefore, the non-carcinogenic risk from heavy metals via accidental ingestion and dermal contact with water near the Sele River estuary was revealed to be within a safe range for the exposed population.The carcinogenic risk of the studied HMs in the Sele River was estimated as the incremental lifetime cancer risk (ILCR).However, for Hg and Zn, CSF values are not available, and consequently, the carcinogenic risk assessment was not carried out [61,69].Similarly, Yang et al. used Cd, Pb, and As for both carcinogenic and non-carcinogenic risk assessments, while they included Hg only in the non-carcinogenic assessment [87].The results indicated that considering the total exposure caused both by ingestion and dermal contact, the carcinogenic risk was negligible for Cd (ILCR mean = 9.28 × 10 −7 ), and Ni (ILCR mean = 6.54 × 10 −7 ).However, it was within the acceptable range of 10 −6 to 10 −4 for As (ILCR mean = 2.94 × 10 −6 ), Pb (ILCR mean = 1.78 × 10 −5 ), and Cr (ILCR mean = 4.91 × 10 −5 ) [87].The acceptable risk limit is expected to be 10 −6 in the case of single element carcinogenic risk and 10 −4 in the event of multi-element carcinogenic risk [88].Additionally, according to Mohammadi et al. [61], for an individual HM, when the ILCR value is less than 1 × 10 −6 , the cancer risk is considered negligible.On the other hand, an ILCR value higher than 1 × 10 −4 suggests that the cancer risk is considerable.For the total of all heavy metals through all exposure routes, the acceptable level is 1 × 10 −5 [81].Among all of the studied HMs, Cr showed the highest ILCR value (ILCR max = 2.90 × 10 −4 ), implying that exposure over a long lifetime via ingestion and dermal contact could increase the probability of cancer to a significant level for exposed people.In contrast, Ni presented the lowest chance of cancer risk.

Conclusions
This study provides data on the occurrence of eight heavy metals (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn) in water and sediment, as well as their spatial and seasonal patterns, in the Sele River estuary (Italy).Furthermore, the study offers the first assessment of the health hazards to humans from HM exposure in the study region using non-carcinogenic and carcinogenic analyses.The study's findings reveal that greater levels of HMs were detected at the river's mouth and 500 m moving southward, implying that HM pollution downstream is more severe than that upstream, which is normally less impacted by human activities.Furthermore, the spatio-temporal distributions revealed that the greatest amounts were discovered during the warmer seasons (July and April), meaning that fresh inputs of these pollutants into the sea are possible, and the SPM phase might serve as a carrier from the sediment to the water phase.Furthermore, the HM contamination status of sediment from the Sele River estuary was moderate, suggesting that the river could be a source of considerable/moderate levels of pollution by HMs, mainly Hg and Cd.In terms of health risk assessment, the study's outcomes showed that there was no detectable noncarcinogenic risk for the analyzed trace elements since the HQs for all studied HMs were lower than the thresholds of concern.The carcinogenic risk, expressed as the incremental lifetime cancer risk (ILCR), was negligible for Cd and Ni and within the acceptable range for As, Pb, and Cr.This study significantly advances environmental and scientific research by comprehensively assessing heavy metal levels in Sele River water and sediment.The spatio-temporal analysis provides a detailed perspective on heavy metal concentrations, aiding the understanding of variations over time and space.Notably, it estimates heavy metal inputs from the Sele River into the Mediterranean Sea, emphasizing the interconnectedness of fluvial and marine ecosystems.The study's characterization of ecological risk in the Mediterranean Sea and assessment of human health risk contribute to a broader understanding of environmental impacts.Its multi-faceted analysis, spanning local and Mediterranean contexts, methodology, and extensive scope, positions it as a crucial reference for future studies on heavy metal contamination and associated risks to ecosystems and human health.To address heavy metal contamination in estuarine environments, key actions involve implementing continuous monitoring systems for timely detection of variations, regulating industrial and urban discharges, and engaging in ecosystem regeneration projects like wetland creation and reforestation.Community involvement is vital, requiring awareness initiatives for sustainable practices.Comprehensive environmental impact studies are essential before new activities are undertaken near estuaries to prevent negative effects on water and sediment quality.Controlling terrestrial pollution sources, fostering

Figure 1 .
Figure 1.Map of the study area and sampling locations in the Sele River estuary, southern Italy.

Figure 1 .
Figure 1.Map of the study area and sampling locations in the Sele River estuary, southern Italy.
37 kg/year for Cd, 32.95 kg/year for Cr, 111.02 kg/year for Cu, 60.25 kg/year for Ni, 112.81 kg/year for Pb, and 72.90 kg/year for Zn.The spatial distributions of HMs for DP, SPM, and SED samples assessed at different sampling locations during four seasons are shown in Figures 2, 3 and 4, respectively.

Figure 2 .
Figure 2. Spatio-temporal distribution of heavy metals (HMs) detected in the water dissolved phase (DP) at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South) from the Sele River estuary, in southern Italy, during four sampling campaigns (2020-2021).

Figure 2 .
Figure 2. Spatio-temporal distribution of heavy metals (HMs) detected in the water dissolved phase (DP) at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South) from the Sele River estuary, in southern Italy, during four sampling campaigns (2020-2021).

Figure 4 .
Figure 4. Spatial distribution of heavy metals (HMs) detected in sediment (SED) samples collected at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South) from the Sele River estuary, in southern Italy.

Figure 2
Figure 2 displays the spatio-temporal distribution of HMs detected in DP samples collected during the research period.Higher levels were detected at the mouth and 500 m to the south (500S) during the warmer and drier seasons (July > April).In detail, the highest concentrations were recorded at the mouth for Hg (2.35 µg/L), Cr (14.97 µg/L), Ni (11.15 µg/L), and Pb (14.41 µg/L) in the dry season (July).For As and Cu, the highest amounts

Figure 4 .
Figure 4. Spatial distribution of heavy metals (HMs) detected in sediment (SED) samples collected at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South) from the Sele River estuary, in southern Italy.

Figure 2
Figure 2 displays the spatio-temporal distribution of HMs detected in DP samples collected during the research period.Higher levels were detected at the mouth and 500 m to the south (500S) during the warmer and drier seasons (July > April).In detail, the highest concentrations were recorded at the mouth for Hg (2.35 µg/L), Cr (14.97 µg/L), Ni (11.15 µg/L), and Pb (14.41 µg/L) in the dry season (July).For As and Cu, the highest amounts

Figure 4 .
Figure 4. Spatial distribution of heavy metals (HMs) detected in sediment (SED) samples collected at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South) from the Sele River estuary, in southern Italy.

Toxics 2024 , 28 Figure 5 .
Figure 5. Spatial distributions of heavy metal (HM) concentrations (µg/L) compared for dissolved phase (DP) and suspended particulate matter (SPM) samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Figure 6 .
Figure 6.Spatial distributions of heavy metal (HM) concentrations (mg/kg dw) compared for SPM and SED samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Figure 5 .
Figure 5. Spatial distributions of heavy metal (HM) concentrations (µg/L) compared for dissolved phase (DP) and suspended particulate matter (SPM) samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Toxics 2024 , 28 Figure 5 .
Figure 5. Spatial distributions of heavy metal (HM) concentrations (µg/L) compared for dissolved phase (DP) and suspended particulate matter (SPM) samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Figure 6 .
Figure 6.Spatial distributions of heavy metal (HM) concentrations (mg/kg dw) compared for SPM and SED samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Figure 6 .
Figure 6.Spatial distributions of heavy metal (HM) concentrations (mg/kg dw) compared for SPM and SED samples collected in the study area at 10 sampling locations (Site 1: Mouth; Site 2: River Mouth 500 m North; Site 3: River Mouth 1000 m North; Site 4: River Mouth 1500 m North; Site 5: River Mouth 500 m West; Site 6: River Mouth 1000 m West; Site 7: River Mouth 1500 m West; Site 8: River Mouth 500 m South; Site 9: River Mouth 1000 m South; Site 10: River Mouth 1500 m South).

Figure 8 .
Figure 8. Geoaccumulation index (Igeo) values for individual heavy metals (HMs) estimated at 10 sampling locations near the Sele River estuary, in southern Italy.

Figure 8 .
Figure 8. Geoaccumulation index (I geo ) values for individual heavy metals (HMs) estimated at 10 sampling locations near the Sele River estuary, in southern Italy.
pollutionRegarding the CF, the results obtained are displayed in Figure9.The highest CF values were found for Hg and Cd, with mean values of 3.79 and 1.06, respectively.

Figure 9 .
Figure 9. Contamination factors (CFs) evaluated for sediment (SED) samples collected during the sampling campaign near the Sele River estuary.

Figure 9 .
Figure 9. Contamination factors (CFs) evaluated for sediment (SED) samples collected during the sampling campaign near the Sele River estuary.

Toxics 2024 , 28 Figure 10 .
Figure 10.Comparison of the concentrations obtained in the sediment (SED) samples for As, Hg, Cd, and Cr with threshold effect level (TEL, orange line), probable effect level (PEL, grey line), effect range low (ERL, yellow line), and effect range median (ERM, blue line).

Figure 10 .
Figure 10.Comparison of the concentrations obtained in the sediment (SED) samples for As, Hg, Cd, and Cr with threshold effect level (TEL, orange line), probable effect level (PEL, grey line), effect range low (ERL, yellow line), and effect range median (ERM, blue line).

Figure 10 .
Figure 10.Comparison of the concentrations obtained in the sediment (SED) samples for As, Hg, Cd, and Cr with threshold effect level (TEL, orange line), probable effect level (PEL, grey line), effect range low (ERL, yellow line), and effect range median (ERM, blue line).

Figure 11 .
Figure 11.Comparison of the concentrations obtained in the sediment (SED) samples for Cu, Ni, Pb, and Zn with threshold effect level (TEL, orange line), probable effect level (PEL, grey line), effect range low (ERL, yellow line), and effect range median (ERM, blue line).

Figure 11 .
Figure 11.Comparison of the concentrations obtained in the sediment (SED) samples for Cu, Ni, Pb, and Zn with threshold effect level (TEL, orange line), probable effect level (PEL, grey line), effect range low (ERL, yellow line), and effect range median (ERM, blue line).

Figure 12 .
Figure 12.Graphical representation of the hazard quotients (HQs) obtained in this study.

Figure 12 .
Figure 12.Graphical representation of the hazard quotients (HQs) obtained in this study.

Table 1 .
Sampling locations with identification number (ID), characteristics of each site, coordinates, and abbreviation.

Table 1 .
Sampling locations with identification number (ID), characteristics of each site, coordinates, and abbreviation.

Table 4 .
Threshold effect level (TEL), probable effect level (PEL), effect range low (ERL) and effect range median (ERM) values for each heavy metal.

Table 5 .
[61]meters used for risk assessment due to ingestion and dermal exposure caused by heavy metals in water[61].

Table 6 .
Reference oral doses (RfDs) (mg/kg/day) through oral ingestion and dermal contact by heavy metals in the water system.n.a.: not available.

Table 7 .
Cancer slope factors (CSFs) (mg/kg/day) used for the carcinogenic risk assessment.
µg/L (mean value of 25.36 ± 30.99 µg/L) for Cu, from nd to 88.74 µg/L (average value of 11.81 ± 17.53 µg/L) for Ni, from nd to 86.66 µg/L (with a mean value of 19.44 ± 21.72 µg/L) for Pb, and from nd to 107.54 µg/L (average value of 14.45 ± 24.20 µg/L) for Zn.The results found for SPM at 10 sampling locations near the Sele River estuary during four seasons (2020-2021) are displayed in Table Toxics 2024, 12, x FOR PEER REVIEW 10 of 28 value of 2.22 ± 2.34 mg/kg dw) for Cd, Cr, and Cu, respectively.Additionally, the mean value detected for Ni was 4.26 ± 7.80 mg/kg dw (range of nd-24.89mg/kg dw) and 6.72 ± 10.31 mg/kg dw (range of nd -34.57mg/kg dw) for Pb.Finally, the levels found for Zn were in the range of 5.92-42.00mg/kg dw (with a mean value of 17.62 ± 12.29 mg/kg dw).

Table 8 .
I geo values for HMs detected in SED samples from the Sele River estuary, in southern Italy.

Table 9 .
Pollution indices assessed for sediment (SED) samples collected during the sampling campaign near the Sele River estuary.

Table 9 .
Pollution indices assessed for sediment (SED) samples collected during the sampling campaign near the Sele River estuary.Regarding the CF, the results obtained are displayed in Figure9.The highest CF values were found for Hg and Cd, with mean values of 3.79 and 1.06, respectively.

Table 10 .
[58]y metal (HM) concentrations (µg/L) detected in dissolved phase (DP) samples collected from the study area with National Recommended Water Quality Criteria by the US EPA[58].

Table 10 .
[58]y metal (HM) concentrations (µg/L) detected in dissolved phase (DP) samples collected from the study area with National Recommended Water Quality Criteria by the US EPA[58].

Table 11 .
Minimum, maximum, and mean values of CDIs (mg/kg/day) for HMs through different pathways.

Table 12 .
Minimum, maximum, and mean values of HIs evaluated for eight HMs (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn) through different pathways.

Table 12 .
Minimum , maximum, and mean values of HIs evaluated for eight HMs (As, Hg, Cd, Cr, Cu, Ni, Pb, and Zn) through different pathways.

Table 13 .
The incremental lifetime cancer risk (ILCR) values obtained for total exposure due to both ingestion and dermal contact in the study area.