The Impact of Organic Matter on Polycyclic Aromatic Hydrocarbon (PAH) Availability and Persistence in Soils

Polycyclic aromatic hydrocarbons (PAHs) exhibit persistence in soils, and most of them are potentially mutagenic/carcinogenic and teratogenic for human beings but also influence the growth and development of soil organisms. The PAHs emitted into the atmosphere are ultimately deposited (by dry or wet deposition processes) onto the soil surface where they tend to accumulate. Soil organic matter (SOM) plays an important role in the fate and transformation processes of PAHs, affecting their mobility, availability, and persistence. Therefore, the aim of this research was to investigate the influence of SOM fractional diversification (fulvic acids—FA, humic acids—HA, and humins—HN) on PAH availability and persistence in soils. Twenty soil samples (n = 20) were collected from upper horizons (0–30 cm) of agricultural soils exposed to anthropogenic emissions from industrial and domestic sources. The assessment of PAH concentrations included the determination of medium-molecular-weight compounds from the US EPA list: fluoranthene—FLA, pyrene—PYR, benz(a)anthracene—BaA, and chrysene—CHR. The assessment was conducted using the GC-MS/MS technique. Three operationally defined fractions were investigated: total extractable PAHs (TE-PAHs) fraction, available/bioavailable PAHs (PB-PAHs) fraction, and nonavailable/residual PAHs (RE-PAHs) fraction, which was calculated as the difference between total and available PAHs. TE-PAHs were analyzed by dichloromethane extraction, while PB-PAHs were analyzed with a hydrophobic β-cyclodextrin solution. SOM was characterized by total organic carbon content (Turin method) and organic carbon of humic substances including FA, HA, HN (IHSS method). Concentrations of PAHs differed between soils from 193.5 to 3169.5 µg kg−1, 4.3 to 226.4 µg kg−1, and 148.6 to 3164.7 µg kg−1 for ∑4 TE-PAHs, ∑4 PB-PAHs, and ∑4 RE-PAHs, respectively. The ∑4 PB-PAHs fraction did not exceed 30% of ∑4 TE-PAHs. FLA was the most strongly bound in soil (highest content of RE-FLA), whereas PYR was the most available (highest content of PB-PYR). The soils were characterized by diversified total organic carbon (TOC) concentration (8.0–130.0 g kg−1) and individual SOM fractions (FA = 0.4–7.5 g kg−1, HA = 0.6–13.0 g kg−1, HN = 0.9–122.9 g kg−1). FA and HA as the labile fraction of SOM with short turnover time strongly positively influenced the potential ∑4 PAH availability (r = 0.56 and r = 0.52 for FA and HA, respectively). HN, which constitutes a stable fraction of organic matter with high hydrophobicity and poor degradability, was strongly correlated with ∑4 RE-PAHs (r = 0.75), affecting their persistence in soil.

PAHs-polycyclic aromatic hydrocarbons, TE-total extractable PAHs fraction, PB-potential available PAHs fraction, RE-nonavailable/residual PAHs fraction, FLA-fluoranthene, PYR-pyrene, BaA-benz(a)anthracene, CHR-chrysene; Min-the lowest content of the test compound, above the detection limit, Max-the highest content of the analyzed compounds, Med-median, Aver-average, Lower Q-lower quartile, Upper Q-upper quartile, CoV-coefficient of variation, n-number of samples.
The highest concentrations of the available PAHs fraction were recorded for PYR (18.1 µg kg −1 ) and CHR (16.2 µg kg −1 ), while lower concentrations were recorded for FLA (14.2 µg kg −1 ) and BaA (7.8 µg kg −1 ), as shown in Table 1. The PB-PAHs (PB-FLA, PB-PYR, PB-BaA, PB-CHR) accounted for 0.2%-29.4% of TE-PAHs (TE-FLA, TE-PYR, TE-BaA, TE-CHR) with a high diversity of results (coefficient of variation, CoV, was at the level of 103% to 183% for individual compounds, Table 1, Figure 1). The potentially available fraction of FLA, PYR, BaA, and CHR did not exceed 1% in almost 25% of samples, while 15% had a total content of individual PAHs that was higher than 20%. The obtained results indicated that despite the relatively high total content of individual PAHs in soils, Molecules 2020, 25, 2470 4 of 15 the content of their available fraction is relatively small, which indicates a low risk of their transfer to living organisms and groundwater.
Molecules 2020, 25 PAHs-polycyclic aromatic hydrocarbons, TE-total extractable PAHs fraction, PB-potential available PAHs fraction, RE-nonavailable/residual PAHs fraction, FLA-fluoranthene, PYRpyrene, BaA-benz(a)anthracene, CHR-chrysene; Min-the lowest content of the test compound, above the detection limit, Max-the highest content of the analyzed compounds, Med-median, Aver-average, Lower Q-lower quartile, Upper Q-upper quartile, CoV-coefficient of variation, n-number of samples.  The contents of the residual fraction of PAHs (RE-PAHs) was very closely related to the total concentrations of these compounds (Table 1) (2018), who analyzed the total and potentially available concentrations of PAHs in soils from industrial areas subjected to long-term anthropopressure. A large array of chemical and physicochemical methods applied for the measurement of the availability of PAHs does not enable proximate comparison of the results because the final concentration of the PB-PAHs fraction results from the extraction procedure [32,37]. The methods involved either a mild organic solvent extraction or application of cyclodextrins or solid-phase hydrophobic adsorbents such as XAD resin or Tenax-TA [7,10,12,13,32,[38][39][40][41][42]. Regardless of the applied method, the quantity of available PAHs was significantly lower compared to the total concentration that indicated the low aqueous solubility of PAHs and/or their strong sorption on soil components and partitioning among various phases [14][15][16]. Cyclodextrin extraction was firstly identified by Reid et al. [38] and then confirmed by other authors [13,31,39,41,43] as a suitable technique that recognizes the readily desorbing fraction available for biota, causing a real risk in the contaminated areas.
The TE-PAH compounds were not significantly correlated with PB-PAHs, which may indicate the historical deposition of hydrocarbons and occurring aging processes. Some authors observed significant positive relationships indicating that soils with a higher concentration of TE-PAHs were characterized by a higher content of PB-PAHs [13,39]. Nevertheless, the research was conducted in soils artificially fortified by these compounds in laboratory conditions. Availability of PAHs in these soils was quite high due to a short incubation time (from one day to several months) that influenced aging processes [40,42]. In the natural soil environment, many different factors may affect PAH behavior, which is mainly associated with the repetitive depositions of contaminants and strong sorption of these pollutants by various soil components occurring over a long period of time as well as diffusion of PAH molecules within structures of organic matter [3,7,11,12,[14][15][16]32,[40][41][42]. The importance of the aging processes was confirmed by the high share of the RE-PAH fraction of individual compounds.

Assessment of Soil Organic Matter (SOM) and Its Fractional Composition (FA, HA, HN)
SOM content in the analyzed soils exhibited typical values for Polish soils of agricultural usage [44] and was in the range 8.0-130 g kg −1 (CoV = 120%) with a median value of 13.8 g kg −1 . SOM was fractionated onto individual HS components, the content of which was not dependent on the total content of organic matter ( Table 2). HS was highly diversified, and among the compounds, HN constituted the most heterogeneous fraction (CoV = 192%) compared to FA (CoV = 68%) and HA (CoV = 43%). Similarly, HN constituted a significant portion (0.9-122.9 g kg −1 ) of humic substances, while more labile fractions of humic acids occurred in the range from 0.4 to 7.5 g kg −1 and 0.6 to 13 g kg −1 for FA and HA, respectively. A general share of individual fractions (in total SOM content) followed the overall trend of HN (10.3%-94.5% of total organic carbon-TOC) > HA (4.6%-66.4% of TOC) > FA (0.3%-17.7% of TOC), as shown in Figure 2.  [44] and was in the range 8.0-130 g kg −1 (CoV = 120%) with a median value of 13.8 g kg −1 . SOM was fractionated onto individual HS components, the content of which was not dependent on the total content of organic matter ( Table 2). HS was highly diversified, and among the compounds, HN constituted the most heterogeneous fraction (CoV = 192%) compared to FA (CoV = 68%) and HA (CoV = 43%). Similarly, HN constituted a significant portion (0.9-122.9 g kg −1 ) of humic substances, while more labile fractions of humic acids occurred in the range from 0.4 to 7.5 g kg −1 and 0.6 to 13 g kg −1 for FA and HA, respectively. A general share of individual fractions (in total SOM content) followed the overall trend of HN (10.3%-94.5% of total organic carbon-TOC) > HA (4.6%-66.4% of TOC) > FA (0.3%-17.7% of TOC), as shown in Figure 2. The mutual proportions between HA and FA were above one for all soil samples (HA/FA = 1.0-18.2 with a median value of 3.7), while the ratio of (FA + HA)/HN was between 0.1 and 8.1 (with a median value of 1.9), as shown in Table 2. The proportions between HA and FA fractions are expressed as HS transformations and indicate the potential low mobility of C in the soil system and a high intensity of humification processes, while the proportion of (HA + FA)/HN shows the degree of organic matter transformation to the stable organic matter forms [24,45]. Higher values of the above two parameters suggest a lower mobility of organic matter and its weaker susceptibility to mineralization.    The mutual proportions between HA and FA were above one for all soil samples (HA/FA = 1.0-18.2 with a median value of 3.7), while the ratio of (FA + HA)/HN was between 0.1 and 8.1 (with a median value of 1.9), as shown in Table 2. The proportions between HA and FA fractions are expressed as HS transformations and indicate the potential low mobility of C in the soil system and a high intensity of humification processes, while the proportion of (HA + FA)/HN shows the degree of organic matter transformation to the stable organic matter forms [24,45]. Higher values of the above two parameters suggest a lower mobility of organic matter and its weaker susceptibility to mineralization.

Influence of SOM on Concentrations of Total Extractable PAH (TE-PAHs), Potentially Available PAHs (PB-PAHs), and Residual PAHs (RE-PAHs) in Soil
The main factor controlling the PAH content in soils is the SOM and its fractions [5,[9][10][11][12][17][18][19][20][21][22][23][24], as shown in Table 3. Although the TOC content does not affect directly the concentration of the potentially available PAHs fraction, a higher organic carbon content generally results in increased sorption of PAHs and a reduction of their accessibility. It is assumed that soils amended with exogenous organic matter increase pollutant sorption and reduce their absorption by living organisms. Many different studies have been conducted in this field [11,12,21,23]. The origin and SOM composition can profoundly influence the association and immobilization of PAHs. Different fractions of organic matter bind PAHs with different strengths and affect their half-life and resistance in soils [18][19][20][21][22]24]. Results of correlations between PB-PAHs, TE-PAHs, and RE-PAHs and FA, HA, and HN clearly indicate that availability of contaminants is promoted by labile humic substances (FA, HA), while persistence depends on the stabile carbon forms (Table 3, Figure 3). Nevertheless, the impact of individual SOM fractions on the PAH behavior may be also controlled by the properties of the contaminants and their affinity with the organic soil phase, expressed by the K o/w coefficient (Table 5).
According to correlation results, FA and HA exhibited significant positive correlations with 4 PB-PAHs (r = 0.56 and r = 0.52, respectively), while 4 TE-PAHs were strongly positively correlated with HN (r = 0.70), as shown in Table 3. HN also significantly influenced 4 RE-PAHs (r = 0.75); however, this resulted from the prevalent share of 4 RE-PAHs in the total content of these compounds. Among all PAHs, FLA indicated highest sorption affinity with FA and HA, regardless of the PAH operational fraction (r = 0.50 and r = 0.56 for PB-PAHs; r = 0.57 and r = 0.59 for TE-PAHs; r = 0.50 and r = 0.56 for RE-PAHs). The PYR and BaA occurring in potentially available forms were strongly connected with HA (r = 0.55 for PB-PYR and PB-BaA) in contrast to CHR, which was more strongly associated with FA (r = 0.73 for PB-CHR). The obtained dependencies indicate that the higher contribution of humic acids (FA, HA) to total SOM content caused the greater availability and potential mobility of hydrocarbons in the soils, while high amounts of HN affected their higher stability.  The results clearly showed (Table 3 and Figure 3) that in soils subjected to long-term anthropopressure, aging processes of the pollutants prevail. PAHs directly deposited in soils are subject to an "aging process", which often entails an initially rapid and reversible sorption process to SOM followed by a period of slow diffusion occurring over weeks, months, or even years once they are released into the soil environment [15][16][17]. This indicates stronger PAH sorption by stable SOM structures and the persistence of PAHs in soil. During aging, contaminants persevere in soils and become increasingly recalcitrant to desorption and thus much less available, possessing lower bioaccumulation and biodegradation affinity or toxicity [16,17,19,20,31,32,38,39]. Therefore, the resistant fraction content increases with time, indicating significant shares of RE-PAHs in the total PAH content.
Aging involves a number of mechanisms but is predominated by diffusion into micro-or nanopores [16,17] and partitioning into soil organic matter [18][19][20]22]. It is generally thought that For a more detailed recognition of the relationship between the effect of the SOM fractions on the PAH content, PCA analysis was used. All analyzed data are represented by three main components, explaining 89% of the total variance of the results ( Figure 3A), whereas up to 77% of variance is explained by the first two factors ( Figure 3B). The first PCA component (PCA 1), which accounted for 52% of variance, was significantly negatively correlated with individual TE-PAHs and RE-PAHs (r = −0.81 to r = −0.98 and r = −0.82 to r = −0.98, respectively), while the second PCA component (PCA 2), which represented 25% of the results' variability, was positively correlated with PB-PAHs (r = 0.64 to r = 0.85). All variables are well represented by the first two main components of the PCA coordinates, but the impact of the third PCA component (PCA 3) reached almost 12%. Although the loads were not correlated with any of the analyzed variables, its share may be related to the chemical properties of hydrocarbons or other undefined processes or properties affecting the behavior of PAHs in soils. These results more clearly indicate that individual SOM fractions determine the accumulation of PAHs and their availability in soils. PAHs were divided into two independent groups, and the distance between variables confirmed a high mutual relationship pointing to significant soil parameters affecting their contents.
The results clearly showed (Table 3 and Figure 3) that in soils subjected to long-term anthropopressure, aging processes of the pollutants prevail. PAHs directly deposited in soils are subject to an "aging process", which often entails an initially rapid and reversible sorption process to SOM followed by a period of slow diffusion occurring over weeks, months, or even years once they are released into the soil environment [15][16][17]. This indicates stronger PAH sorption by stable SOM structures and the persistence of PAHs in soil. During aging, contaminants persevere in soils and become increasingly recalcitrant to desorption and thus much less available, possessing lower bioaccumulation and biodegradation affinity or toxicity [16,17,19,20,31,32,38,39]. Therefore, the resistant fraction content increases with time, indicating significant shares of RE-PAHs in the total PAH content.
Aging involves a number of mechanisms but is predominated by diffusion into micro-or nanopores [16,17] and partitioning into soil organic matter [18][19][20]22]. It is generally thought that PAH sorption to SOM is influenced by two major domains, "soft/rubbery" and "hard/glassy" domains. Wang et al. [9] proposed the theory of sequential bonding of pollutants by SOM, which involves the following steps: Step 1, capturing the fraction of fresh pollutants through rubber FA (weak hydrogen bonds with PAH, facilitating the migration in the soil profile); Step 2, binding through HA (metastable sorbent states, inclusion of HA micellar forms); and Step 3, slow sorption of contaminants into the glassy HN structures (process aging). The noted significant positive correlations between HN and RE-PAHs as well as FA, HA, and PB-PAHs (Table 3) confirm the migration path of hydrocarbons between SOM fractions.
Luo et al. [20] have suggested that surface characteristics of SOM fractions play a critical role in interactions with organic contaminants due to better access to the surface sorption/retention domains in the first stage after they are released into the soil. This implies that the abundance of surface carbon domains of SOM fractions plays an important role in attracting PAHs and in their further interactions. According to Schaumann [29,30] and Pignatello [19], the "soft/rubbery" and "hard/glassy" phases interpenetrate in the HS structure. HA and FA are mainly of amorphic nature with mostly "soft/rubbery" phases, which very quickly interact with contaminants, but the strength of the bond is weak. HN are characterized by higher condensation of structure (benzene rings and aliphatic chains) and thus are more abundant in "hard/glassy" phases, and sorption on their structures occurs slowly, but the sorbed molecules are relatively persistent [19,29,30,46]. The "hard/glassy" phase forms regions of disordered side chains, which have a high density and low reactivity and flexibility due to the presence of numerous unsaturated bonds [29,46]; therefore, they have a high capacity for retaining PAHs by diffusing them into the tertiary structure. The structure of the condensed and rubbery phase results in more active and stable binding sites for PAHs, which allows the effective retention of these compounds in soils [17,29,30]. Although general information about chemical and physical differences between various organic matter fractions have been previously analyzed [17][18][19]22,24,29,30], the operational separated humic substances may have other specific properties that promote their interactions with PAHs.
Generally, the mechanisms of PAH partitioning during aging are thought to involve the movement of PAHs to more recalcitrant binding sites, e.g., micropores that are too small for microbial access or deeper penetration into organic matter by encapsulation (occlusion) or strong sorption by SOM components [19,46]. The positive correlation between RE-PAHs and organic matter may confirm this. Diffusion of contaminants within organic matter structure and soil micropores may also result in their migration to the more "relaxed regions", where the desorption of molecules to the soil solution is easy, or to regions where they are weakly bound by organic matter that induces the release of contaminants, resulting in an increase in the available PAH content.

Site Description and Sample Collection
The SOM and its fractions (FA, HA, HN) as well as the individual 4-ringed PAHs were isolated from soils subjected to long-term and strong industrial anthropogenic pressure (since 1910 when mining activity began on the study area). The sampling points (n = 20) were located in the Czerwionka-Leszczyny municipality (area: 115 km 2 ; population density: 366 persons km −2 ; total dust emission: 6827 kg year −1 km −2 ; dust emission from industrial sources: 1084 kg year −1 km −2 ) situated in southwestern Poland (Silesia Voivodeship, Figure 4). The research area is covered by soils of agricultural, urban, and industrial usage belonging to the Rybnik Industrial District.  The main contamination sources of soils from the Czerwionka-Leszczyny area derived from fuel combustion, road transport, and industrial activities predominantly related to the disposal of environmentally harmful mining waste and production of bituminous masses. Industrial activity affected the strong transformation of the research area, resulting in the occurrence of numerous waste The main contamination sources of soils from the Czerwionka-Leszczyny area derived from fuel combustion, road transport, and industrial activities predominantly related to the disposal of environmentally harmful mining waste and production of bituminous masses. Industrial activity affected the strong transformation of the research area, resulting in the occurrence of numerous waste dumps and storage sites for metallurgical wastes. Moreover, the single-family housing clusters were an additional source of local domestic PAH emissions. A broader description of the study area is available in the publication of Klimkowicz-Pawlas et al. [5] and Ukalska-Jaruga et al. [24].
The soil samples were collected from the surface layer (0-30 cm) of agricultural soils situated at different distances and locations from potential PAH emission sources that influenced the concentration of these compounds and their contaminant-soil contact time (aging). The location of the sampling points ( Figure 4) reflected the variability of soil physicochemical properties, which may influence differences in SOM content and its fractional composition ( Table 4). The samples were properly prepared prior to chemical analysis by air-drying, passing them through a 2-mm mesh sieve, and storing them in the dark at 12-16 • C. TOC-total organic carbon, TN-total nitrogen, TC-total carbon, CEC-cation exchange capacity; Min-the lowest value of the parameter, Max-the highest value of the parameter, Med-median, Aver-average, LQ-lower quartile, UQ-upper quartile, CoV-coefficient of variation.

Determination of PAHs
Three operationally defined PAH fractions were investigated in the study: (i) the total extractable PAHs fraction (TE-PAHs), including compounds extracted by the intensive extraction technique using organic solvents; (ii) the available PAHs fraction (PB-PAHs), which refers to the fraction of contaminants that is biologically available (bioavailable) for soil organisms, present in soil solution, and may be readily desorbed from soil under the influence of a complexing agent; and (iii) the residual/(nonavailable) PAHs fraction (RE-PAHs), which describes the bound fraction of contaminants and was calculated as the difference between the total and potentially available PAHs (RE-PAHs = TE-PAHs − PB-PAHs).
Determined PAHs included 4-ringed compounds fluoranthene-FLA, pyrene-PYR, benz(a)anthracene-BaA, and chrysene-CHR, classified as medium-molecular-weight hydrocarbons (m/z 202 and 228) and indicated as pollutants of environmental concern on the US EPA list ( Table 5). The selected PAHs exhibit various chemical properties, i.e., spatial conformation, molar mass, water solubility, vapor pressure, and the values of the octanol/water partition coefficient. Moreover, they can be degraded by soil microorganisms and sorbed by SOM. All of them can be derived from both petrogenic and pyrogenic sources and therefore may be representative compounds for the PAH group. The analytical procedure for TE-PAHs was fully described by Maliszewska-Kordybach et al. [1,2] and Ukalska-Jaruga et al. [8,24]. The analytical steps included extraction of ground soil samples (grain size ≤ 0.10 mm) with dichloromethane in an ASE200 Accelerated Solvent Extractor (Dionex Co., Sunnyvale, CA, USA). Before extraction, the samples were mixed with 2 g of diatomaceous earth dried at 550 • C and spiked with 10 µL of a recovery standard containing five deuterated PAHs (PAH-31, Dr. Ehrenstorfer GmbH, Augsburg, Germany). Clean-up of the concentrated extract was performed using hexane followed by hexane/dichloromethane (2:3 v/v) elution on activated silica gel (16 h at 135 • C).
PB-PAHs were extracted using hydroxypropyl-β-cyclodextrin (HPCD) solution according to ISO/TS 16751:2018. HPCD was applied as a complexing agent encapsulating PAH free molecules desorbing from soil to the extracting solution. Briefly, a dry soil sample of 4 ± 1 g was placed in the 100 mL centrifuge vial, and then 40 mL of 0.1 mol dm −3 HPCD solution was added. Extraction was carried out using the horizontal shaker (SM-30, Edward Bühler GmbH, Bodelshausen, Germany) equipped with a thermostat (TM-30, Edward Bühler GmbH, Bodelshausen, Germany) and a minichiller (Huber) in the dark, at 20 • C for 20 h at an amplitude of 150-180 rpm. After extraction, the soil/HPCD mixtures were centrifuged, then a certain amount of water phase (>30 mL) was collected in the separatory funnel and mixed with petroleum ether for liquid/liquid extraction. The organic phase was collected, evaporated (Vario CVC 2000, Vacuumbrand), concentrated in n-hexane and subjected to GC-MS/MS analysis (Agilent 7890B GC system; Agilent Tech., Santa Clara, CA, USA equipped with an Agilent 7000C detector and Agilent 7693 Autosampler). Soil samples were extracted in duplicates. The blank sample procedure was included in all extraction series.
The concentrations of TE-PAHs and PB-PAHs were determined by gas chromatography triple mass spectrometry on an Agilent 7890B GC system (Agilent Tech., Santa Clara, CA, USA), equipped with an Agilent 7000C detector and an Agilent 7693 Autosampler. Sample analyses were performed in multiple reaction monitoring (MRM) mode with diagnostic ions as recommended in ISO 22892:2006. Quality control measures included analysis of a certified reference material (CRM 131, ANAB Accredited Tested Laboratory), duplicate matrix samples, and a solvent blank sample. The precision expressed as a relative standard deviation (RSD) was in the range from 5% to 12%, and the recovery for individual compounds from CRM 131 was within 62% and 84%. The limit of quantification (LoQ) for individual PAH compounds ranged from 0.02 to 2.10 µg kg −1 , while the limit of detection (LoD) was within the 0.01 to 0.81 µg kg −1 range.

Determination of SOM and Its Fractions
The SOM fractionation was based on the isolation of the HS individual fractions, comprising FA, HA, and HN. The procedure was carried out using an adapted ISO 12782-4 method (ISO 12782-4, 2012, [47]), approved by the International Humic Substances Society [27]. The soil samples were acidified by 0.1 M HCl at pH 1 (1:10 w/v soil to solution ratio), shaken for 1 h in a laboratory shaker, and then centrifuged. The supernatant containing the first extract of the low molecular fulvic fraction was removed from the residue by decantation and discarded. The remains of the soil sample were neutralized with 0.1 M NaOH and then extracted with 1 M NaOH (1:10 w/v at pH > 12) overnight by continuous shaking. The decanted supernatant containing the mixture of FA and HA was acidified with 6 M HCl to gravimetrically precipitate the HA. The solutions of FA and HA were purified on DAX-8 polymer resin and analyzed for organic carbon content by a C-N analyzer (Multi N/C 2100/2100S Jena Analytics). The HN content was determined in the remaining soil solid residue after FA and HA extraction. Concentrations of organic carbon in the HN fractions were determined by a C-N analyzer (Vario Macro Cube, Elementar Analysensysteme GmbH, Langenselbold, Germany).

Statistical Analysis
Statistica software (Dell Statistica, version 13.3, TIBCO Software Inc. Greenwood Village, CO, USA) was used for statistical analysis. Basic statistical parameters such as mean, median, extreme values (min, max), lower quartile (LQ), upper quartile (UQ), and coefficient of variation (CoV) were calculated. The chi-squared test was applied to check the normal distribution of the results, and the principal component analysis (PCA) was used to provide an overview of the groups of PAHs and soil organic matter fractions. The PCA used an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of values of linearly uncorrelated variables called principal components. The number of significant PCA factors was determined based on the scree plot. All analyzed relationships between variables marked as statistically significant were measured as p < 0.05.

Conclusions
The assessment of the availability and persistence of PAHs in soils subjected to long-term anthropopressure from industrial and domestic sources indicated that the behavior of contaminants depends on the SOM content and its fractional composition. The results pointed out that availability of PAHs is promoted by labile organic matter fractions including FA and HA, while their persistence depends on the relatively stable carbon forms, such as HN. Therefore, the higher contribution of humic acids in total SOM content raises the potential mobility and bioavailability/bioaccessibility of hydrocarbons to living organisms. Therefore, the processes of SOM transformations (e.g., expressed by HA/FA proportions) significantly affect the sorption and form of PAHs in soils. This indicates that the transformation process of organic matter constitutes an important factor in the intra-soil processing responsible for the retention of PAHs in soils. The more advanced the humification processes, the greater persistence of PAHs and the lower the availability, causing lower bioaccumulation, biodegradation, and toxicity of PAHs.