Investigation of polycyclic aromatic hydrocarbons in soils from Caserta provincial territory, southern Italy: Spatial distribution, source apportionment, and risk assessment

[Display omitted] •Systemic data on PAH concentrations in soils from Caserta province are reported.•The high levels of PAHs observed in Caserta are linked to pollution hotspots.•Pyrogenic emissions have a dramatic influence on the PAH concentrations.•PAH contamination in soils did not damage the ecological environment.•Soils contaminated with PAHs posed health risks to local residents. The concentrations of polycyclic aromatic hydrocarbons (PAHs) in soils from Caserta provincial territory, southern Italy, were systematically investigated along with their correlations with soil properties and health risk. The concentrations of ∑16PAHs ranged from 10.0 to 4191 ng/g, with a median (1 st quartile, Q1; 3rd quartile, Q3) of 28.5 (17.5–65.0) ng/g; Four-ring PAHs were the most abundant and contributed an average of ∼50.2% of the ∑16PAHs. Significant differences in the spatial distributions of PAHs in soil were observed, with higher levels of PAH contamination found in Caserta city and the surrounding areas. According to the positive matrix factorization (PMF) model, three sources were identified: chemical production and metal smelting, vehicle emissions, and coal/biomass combustion. Soil total organic carbon was significantly correlated with the concentration of total PAHs and the concentrations of PAHs with three-, four-, and five-rings. In contrast, only the concentration of ∑4DBPs (dibenzo(a,e)pyrene, dibenzo(a,h)pyrene, dibenzo(a,i)pyrene, dibenzo(a,l)pyrene) was well correlated with population density. The soil mass inventory of ∑16PAHs was estimated to be 6.87 metric tons (geometric mean). The ecological risks posed by PAHs in the study are negligible; however, health risks of exposure to soil-borne PAHs were identified based on a probabilistic risk model.