Abstract

The Niger Delta, Nigeria, is noted for crude oil exploration. Whereas there seems to be a handful of data on soil polycyclic aromatic hydrocarbon (PAH) levels in this area, there is a paucity of studies that have evaluated soil and vegetation PAHs simultaneously. The present study has addressed this information gap. Fresh Panicum maximum (Jacq) (guinea grass), Pennisetum purpureum Schumach (elephant grass), Zea mays (L.) (maize), and soil samples were collected in triplicate from Choba, Khana, Trans-Amadi, Eleme, Uyo, and Yenagoa. PAHs determination was carried out using GC-MS. The percentage composition of the molecular weight distribution of PAHs, the molecular ratio of selected PAHs for identification of possible sources, and the isomeric ratio and total index of soil were evaluated. Pennisetum purpureum Schumach (elephant grass) from Uyo has the highest (10.0 mg·kg−1) PAH while Panicum maximum (Jacq) (guinea grass) has the highest PAH (32.5 mg·kg−1 from Khana. Zea mays (L.) (maize) from Uyo (46.04%), Pennisetum purpureum Schumach (elephant grass) from Trans-Amadi (47.7%), guinea grass from Eleme (49.2%), and elephant grass from Choba (39.9%) contained the highest percentage of high molecular weight (HMW) PAHs. Soil samples from Yenagoa (53.5%) and Khana (55.3%) showed the highest percentage of HMW PAHs. The total index ranged 0.27–12.4 in Uyo, 0.29–8.69 in Choba, 0.02–10.1 in Khana, 0.01–5.53 in Yenagoa, 0.21–9.52 in Eleme, and 0.13–8.96 in Trans-Amadi. The presence of HMW PAHs and molecular diagnostic ratios suggest PAH pollution from pyrogenic and petrogenic sources. Some soils in the Niger Delta show RQ(NCs) values higher than 800 and require remediation to forestall ecohealth consequences.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous and persistent organic pollutants consisting of two or more fused rings [13]. Their structural stability, semivolatility, and hydrophobicity make them ubiquitous in the environment [47]. PAHs have been implicated in diverse toxicities including mutagenicity, teratogenicity, neurotoxicity, genotoxicity, and carcinogenicity in humans [812]. Given these public health concerns, 16 PAHs have been recommended for the priority control list by the United States Environmental Protection Agency (USEPA) due to their carcinogenicity, neurotoxicity, and genotoxicity [9]. Furthermore, PAHs are known to migrate and transform among different environmental matrices, including those interfacing food chains, thus affecting both animal and human health [13].

The environmental persistence and public health importance of PAHs in recent years have attracted global attention [1416]. The principal sources of PAHs in different environmental matrices (soil, water, atmosphere, and food) are natural sources and anthropogenic processes including diverse industrial activities, especially those involving incomplete combustion of coal. Similarly, crude oil and petroleum manufacturing process also produce significant amounts of PAHs [1719].

Although the uptake of PAH by leaves is mainly by gaseous deposition [20], some reports have demonstrated that leaves can accumulate PAHs from contaminated soils through their roots [2123]. Nevertheless, Su and Zhu [24] suggested that the PAHs transported from roots to shoots may be negligible, with other researchers affirming that atmospheric PAHs are dominant contributors to the total PAHs in leaves [25, 26].

Long periods of crude oil exploitation in many countries have led to complex contamination by petroleum hydrocarbons and PAHs in many cities. For instance, some studies from China have reported an average content of 16 PAHs (∑16 PAHs) of 1840 μg/kg in Daqing street dust [27]. All in all, it is known that the pollution effects of PAHs on ecological land in oil- producing cities and surrounding communities have great significance on urban ecological security and environmental health [28, 29]. The Niger Delta area of Nigeria is one of the major crude oil-exploring regions in the world with inundated cases of oil spills, aerial deposition of organic by-products originating from flared gases and massive environmental degradation by both inorganic and organic pollutants like PAHs [10]. Some researchers have reported concentrations and compositional patterns of PAHs that can be employed in understanding the effects, sources, fate, and transport of PAHs in soils, as well as environmental quality management in the Niger Delta, Nigeria [30, 31]. PAHs strongly accumulate in the food chain and are subsequently transferred to humans, thereby posing a threat to human health [32, 33]. Some skeletal surveys of PAHs in agricultural soils have been carried out in Niger Delta [3436]. Some of these surveys suggest that some agricultural soils in the Niger Delta, Nigeria, suffer from PAH pollution due mainly to point source pollution. Surface and underground water are polluted with PAHs in Nigeria [37, 38]. PAHs can transfer from soil to borehole water [39]. Although some studies have reported soil and vegetation PAH contamination, information remains sparse on PAH contamination of both soils and vegetation from the same location [40]. Understanding of the spatial distribution of PAHs in agricultural topsoil is critical for environmental management and the safety of agricultural produce.

As recently discussed [41], the One Health strategy, including environmental health and food safety, can help risk assessors and risk managers in prioritising actions for the prevention and mitigation of PAH pollution and its spread and accumulation. In the present study, we evaluated the whole ecological risk of PAHs in soil and vegetation samples from farmlands in six major cities of Niger Delta, Nigeria.

2. Materials and Methods

2.1. Study Area

The Niger Delta, Nigeria, is the third largest mangrove forest and the second largest delta in the world. It falls within the central coastlands of southern Nigeria [10]. The Niger Delta area is well known for crude oil exploration and environmental pollution. According to Okoye et al. [10], the black to greyish brown and dark grey soils are acidic with slight to moderate electrical conductivity and high organic carbon content.

2.2. Sampling and Sample Treatment
2.2.1. Detailed Sampling Locations (GPS)

Soil and fresh plant (guinea grass, elephant grass, and maize) samples were sampled in Choba-4°53′16.6″N 6°54′31.3″E, Khana-4°39′36.0″N 7°22′39.4″E, Trans-Amadi-4°48′53″N 7°2′14″E, Eleme-4°47′24.2″N 7°07′47.6″E, Uyo-5°01′59.0″N 7°56′15.8″E, and Yenegoa-4°55′29″N 6°15′51″E in the Niger Delta. Figure 1 shows the sampling sites. The sampling, done in triplicate, was carried out during the period of mass flowering of plants, specifically in January 2018. Soil samples were sampled from 0 to 5 cm depth as the most root-inhabited soil layer according to the US EPA Method 610 (U.S. EPA, 1977). The sampled soil was cleaned of plant residues and other inclusions, ground in a porcelain mortar, and passed through a sieve with a hole diameter of 1 mm. Plants were dried and ground up to a hole diameter of 1 mm for analytical analysis.

2.3. Analytical Determinations
2.3.1. Information on LOD, LOQ, and Calibration (QA/QC)

Analyses of PAHs were done using gas chromatography (6890 series and 6890 plus) equipped with a dual detector (FID-ECD), dual column, TriPlus AS autosampler with helium carrier gas, and a quadrupole mass spectrometer (Agilent 5975 MSD) based on the US EPA method 8100. This analytical procedure has been described previously [4245]. Briefly, the extraction of PAHs from the samples was done with a sonicator (ultrasonic bath, Elmsonic S40H) in accordance with US SW-846 Method 3550. Two grams of either soil or plant samples were extracted with a 50 : 50 mixture of acetone and methylene chloride (analytical grade), spiked with 1 ml of PAH internal standard, and shaken thoroughly for proper mixing before being placed in an ultrasonic bath. Thereafter, 2.00 μl of each sample extracts were injected into the GC port set at column conditions: HP-5cross-linked PH-ME siloxane, length of 30 m, I.D: 0.25 mm, thickness of 1 μm with helium carrier gas set in the spitless, constant flow mode with a 1.2 ml/min flow rate. Other GC and MS operating setups were done according to the instrument’s method of development as specified in the operating instruction manual. Identification and quantification of individual PAHs were based on an internal calibration standard containing known concentrations of the 16 PAHs [4345]. The specificity of the 16 PAHs sought in the samples was confirmed by the presence of transition ions (quantifier and qualifier) as shown by their retention times which corresponded to those of their respective standards. The measured peak area ratios of precursor to quantifier ions were in close agreement with those of the standards. The 16 PAHs analyzed were the following: napthalene (Nap) CAS No. 91-20-3, acenaphthylene (Acy) CAS No. 208-96-8, acenapthene (Ace) CAS No. 83-32-9, fluorene (Flu) CAS No. 86-73-7, phenanthrene (Phe) CAS No. 85-01-8, anthracene (Ant) CAS No. 120-12-7, fluoranthene (Flt) CAS No. 206-44-0, pyrene (Pyr) CAS No. 129-00-0, benzo[a]anthracene (BaA) CAS No. 56-55-3, chrysene (Cry) CAS No. 218-01-9, benzo[b]fluoranthene (BbF) CAS No. 205-99-2, benzo[k]fluoranthene (BkF) CAS No. 207-08-9, benzo[a]pyrene (BaP) CAS No. 50-32-8, dibenzo[ah]anthracene (DahA) CAS No. 53-70-3, benzo[ghi]perylene (BghiP) CAS No. 191-24-2, and indeno[1,2,3-cd]pyrene (Ind) CAS No. 193-39-5. The CAS No. is a unique identification number assigned by the Chemical Abstracts Service (CAS), US, to every chemical substance described in the open scientific literature.

The detection limit (LOD) is estimated as three times the background noise (IUPAC criterion). The blank samples remained always below the quantification limit (LOQ). Table S1 shows the reproducibility relative standard deviation (RSDr; n = 6), repeatability relative standard deviation (RSDr; n = 6), recoveries, linear range, LOQ, LOD, and coefficient of estimation (r2).

2.4. Data Analysis

SPSS version 17.0 software (SPSS Inc., USA) was used to perform all statistical analyses. The data were analyzed to test for the significance of observed differences in the PAH content by using analysis of variance (ANOVA), while the Tukey test was used to establish if the observed differences in the mean content of PAHs from the different cities were significant.

The sources of the soil PAHs from different cities were evaluated by using PAH isomeric ratios.

The risk quotient in ecological risk assessment is defined as the level of risk produced by a particular PAH and is estimated from the risk quotient (RQ), as shown in the following equation:where CPAHs is the concentration of certain PAHs in the soil and CQV is the corresponding quality value concentration for these PAHs in the soil. Cao et al. [46] model was adopted in order to obtain the quality value concentrations in the present study.

The negligible concentrations (NCs) and the maximum permissible concentrations (MPCs) of PAHs are two quality values employed here with corresponding risk quotients, RQNCs and RQMPCs, respectively, computed fromwhere CQV(NCs) is the quality value of the NCs and CQV(MPCs) is the quality value of the MPCs of the PAHs in the soil.

The risk caused by a combination of all 16 PAHs can be appraised by calculating RQ∑PAHs(NCs) and RQ∑PAHs(MPCs) in which the values of RQNCs and RQMPCs of the individual PAHs that are not less than one are summed, as shown as follows:where .where .

The RQ denotes as follows: RQNCs <1.0 indicates that the individual PAH compounds are probably of negligible concern, while RQMPCs ≥ 1 suggests severe contamination by the individual PAH compound that requires remediation. RQNCs ≥ 1.0 and RQMPCs <1 indicate moderate risk posed by a single PAH compound that might require some control and remediation. However, and imply moderate risk 1. and indicate that the PAHs constitute a moderate risk 2; and show a high risk of the ∑16 PAHs in the ecosystem.

3. Results and Discussion

3.1. Level of PAHs

Table 1 shows the levels of PAHs in soil and elephant grass, maize, and guinea grass samples from Uyo, Choba, Khana, Yenagoa, Eleme, and Trans-Amadi. Guinea grass from Khana had the highest level of total PAH 32.5, whereas elephant grass from Uyo had the highest level of PAH 10.0. Maize samples from Yenagoa showed the highest level of PAH (9.73). Soil samples from Uyo had the highest soil total PAH (8.80).

3.2. Source Apportionment

Figures 2 and 3 show the PAHs content in elephant grass 2(a), maize 2(b), and guinea grass 2(c) from different locations in the Niger Delta and the abundance of individual PAHs in elephant grass 3(a), maize 3(b), and guinea grass 3(c) from Choba, Eleme, Khana, Trans-Amadi, Uyo, and Yenagoa in the Niger Delta, respectively.

The levels of different PAH content in elephant grass, maize, guinea grass, and soil samples, respectively, from Choba, Eleme, Khana, Trans-Amadi, Uyo, and Yenagoa in the Niger Delta, Nigeria, are provided in Figures 25. Figure 2(a) shows that elephant grass from Uyo had the highest number and content of PAHs whereas Khana had the least number of PAHs and the highest content of BaA (4.22 mg·kg−1). The highest content of benzo[a]anthracene (BaA) in elephant grass and dibenzo[ah]anthracene (DahA) content in guinea grass seen in Khana may be an indication of the severity of PAH pollution, and its impact on ecohealth should be further investigated.

Figure 2(b) shows that the maize from Yenagoa and Eleme had a higher content of different PAHs than Khana which contained a fewer number of PAHs with BaA being the highest. The lowest levels of PAHs were seen in guinea grass 2C from Yenagoa, Eleme, Choba, Trans-Amadi, and Uyo. Guinea grass from Khana, Niger Delta, contain all the PAHs and highest levels of Acy (11.9 mg·kg−1) and BaP (11.1 mg·kg−1) in comparison to other locations in this study (Figure 2(c)).

The PAHs content in soil samples from different locations in the Niger Delta, and the abundance of individual PAHs in soil samples from different locations in the Niger Delta are shown in Figures 4 and 5, respectively. According to Figures 4 and 5, Khana had the least number of PAHs and the highest concentration of B(k)F and Flt in comparison to other locations. The soil sample from Uyo had the highest number and concentration of PAHs.

A comparison of the total soil PAH levels (mean and ranges) from different countries and present study is shown in Table 2.

The PAH concentrations of agricultural soils from Choba, Eleme, Khana, Trans-Amadi, Uyo, and Yenagoa, Niger Delta, were 7417, 6787, 6325, 7972, 8799, and 7394 μg·kg−1, respectively. These PAH levels were higher than most of the soil PAH levels from other countries [61, 69, 70].

Figure 6 shows the total PAHs for the different samples from different locations in the Niger Delta of which Uyo has the highest (10.016 mg·kg−1) in elephant grass and Khana has the highest PAH (32.508 mg·kg−1) in guinea grass and the least in all the other samples. Vegetal levels (2.698 maize-32.508 mg·kg−1 elephant grass) of PAHs were higher than soil sample levels (6.325–8.799 mg·kg−1).

Based on the number of rings or molecular structure, the priority PAHs are classified as follows: low molecular weight (LMW) PAHs, i.e., Naph, Acy, Acen, Flu, Phen, and Anth (containing two and three rings), medium molecular weight (MMW) PAHs, i.e., Flan, Pyr, Chry, and BaA (with four rings), and high molecular weight (HMW) PAHs, i.e., BbF, BkF, BaP, IP, DBahA, and BghiP (with five and six rings) [71]. The percentage composition and molecular weight distribution of PAHs in vegetation, i.e., elephant grass, maize, and guinea grass and soils from Uyo, Choba, Khana, Yenagoa, Eleme, and Trans-Amadi in the Niger Delta, Nigeria, are shown in Table 3 and Figure 7. PAHs were grouped in three classes of LMW, MMW, and HMW. Maize from Uyo (46.0%), elephant grass from Trans-Amadi (47.9%), guinea grass from Eleme (49.2%), and elephant grass from Choba (39.9%) contained the highest percentage of HMW PAHs. Soil samples from Yenagoa (53.5%) show the highest percentage of HMW PAHs. Figure 7 shows the different distribution of PAHs according to their molecular weight: for Uyo, maize has more HMW while in guinea grass, elephant grass, and soil samples, LMW was more abundant. For Choba, samples of elephant grass and maize contain more HMW while guinea grass and soil samples have more LMW. For Khana, elephant grass and maize have more MMW PAH and HMW PAH in soil samples. The same trend follows in Yenagoa, Eleme, and Trans-Amadi with HMW in soil samples, guinea grass, and elephant grass, respectively.

The existence of different homologs of PAHs (PAHs with a number of aromatic rings) and PAHs with different molecular weights in the environment suggests their likely origin or sources [7277].

Although there were higher levels of HMW PAHs in plant/vegetation samples in some cities in the present study, the fairly appreciable presence of 4-ring PAHs and 3-ring PAHs (Table 3 and Figure 7) is indicative of mixed pyrogenic sources. The LMW and MMW PAHs are known to exist both in the vapor and particulate phases [77] and usually reside within the locality of the origin or source. The observation from an additive standpoint that 3-4-rings PAHs predominated in this study suggests localized mixed sources coupled with atmospheric transport [40, 70].

PAHs with less than 4 aromatic rings (LMW PAHs) are typified by grass and industrial oil, wood combustion, and petroleum products (Liu et al. 2017), whereas PAHs with more than 4 aromatic rings (HMW PAHs) signify pyrogenic activities at high temperature including coal combustion and vehicular emissions [78]. The higher levels of HMW PAHs in maize from Uyo, elephant grass from Trans-Amadi, guinea grass from Eleme, elephant grass from Choba, and soil samples from Yenagoa and Khana may implicate vehicular traffic.

Since HMW PAHs tend to reside in closer proximity to emission sources and LMW PAHs are carried to areas far from the emission sources [40, 70], the homolog pattern of PAHs in this study with different molecular weights may be characterized by local combustion sources in addition to atmospheric transported depositions [77, 79].

In addition to using different homologs of PAHs (PAHs with the same number of aromatic rings) and PAHs with different molecular weights, molecular diagnostic ratios of selected PAHs concentrations, including Fle/(Fle + Pyr), Ant/(Ant + Phe), Flt/(Flt + Pyr) in PAH identification, BaA/(BaA + Cry), BbF/BkF, BaP/BghiP, BaP/(BaP + Cry), and Ind/(Ind + BghiP) were also employed as inferential tools in the characterization of possible sources of PAHs from plant and soils samples from Uyo, Choba, Khana, Yenagoa, Eleme, and Trans-Amadi (Table 4). The average ratio of PAHs in plants and soil suggested diverse pyrogenic sources of PAHs emission, such as Fle/(Fle + Pyr) (0.60 (0.59–0.97)) for petrol and diesel [72, 80], Ant/(Ant + Phe) (0.88 (0.61–0.70)) for petroleum and biomass combustions [81], Flt/(Flt + Pyr) (0.45 (0.77–1.67)) for biomass, coal combustion (39), BaA/(BaA + Cry) (0.27 (0.32–1.17)) for petrogenic and combustion of petroleum and biomass [80], BbF/BkF (1.16 (1.27−1.09)) for diesel engine and vehicular emissions [82], BaP/BghiP (0.76 (1.49–1.94)) for vehicular emissions and coal combustions [83], BaP/(BaP + Cry) (0.009 (0.01–0.86)) for gasoline [84], and Ind/(Ind + BghiP) (0.45 (1.30–2.91)) for petrogenic and petroleum combustions (Table 4). Molecular diagnostic ratio analysis suggests disparate combustion activities including petrol, diesel, gasoline, biomass, and coal combustions and vehicular emissions, as the principal sources. Petrogenic sources are also related to automobile workshops and accidental spillage.

The isomeric ratios, namely, Flt/(Flt + Pyr), Ant/(Ant + Phen), Phen/Ant, LMW/HMW, BaA/(BaA + Chry), and IndP/(IndP + BghiP), and the total index of soil from Uyo, Choba, Khana, Yenagoa, Eleme, and Trans-Amadi, Niger Delta, Nigeria, are shown in Table 5. Furthermore, total indexes [8587] defined as the sum of single indices (aforementioned parameters previously discussed), respectively, were calculated and normalized for the limit value (low temperature sources–high temperature sources) [80].

Usually, PAHs associated with combustion or high-temperature processes possess a total index that is greater than 4, whereas PAHs originating from petroleum products or low-temperature processes have a total index that is less than 4. In the present study, the total index ranged 0.27–12.4 in Uyo, 0.29–8.67 in Choba, 0.00–10.1 in Khana, 0.01–5.53 in Yenagoa, 0.21–9.52 in Eleme, and 0.13–8.96 in Trans-Amadi. All the sampling locations had total index values greater than 4. It can therefore be inferred from the total index values that PAHs emanated from low- and high-temperature combustion processes. These observations seem to be in agreement with previous results from Warri [46].

3.3. Ecological Risk Assessment

Table 6 shows the RQ(NCs) and RQ∑PAHs(NCs) and RQ(MPCs) and RQ∑PAHs(MPCs) for the 16 priority PAHs in soils from Uyo, Choba, Khana, Yenagoa, Eleme, and Trans-Amadi, Niger Delta, Nigeria. RQ∑PAHs(NCs) values less than 800 are indicative of low ecological risk of PAHs while RQ∑PAHs(NCs) higher than 800 connote higher ecological risk of PAHs. Although some soil samples from the various sampling locations had RQ(NCs) values greater than 800 in the present study, preponderance of the samples showed RQ(NCs) values less than 800 which is suggestive of a low ecological risk of PAHs in these soils. Some previous studies in Nigeria have reported similar RQ(NCs) values less than 800 [88]. The RQ∑PAHs(MPCs) values were greater than 1 except for soil samples from Khana. In this study, Nap, Acy, Ace, Flu, Ant, Flt, and Pyr were the major causes of the ecological risk of PAHs in soils from Uyo, Choba, Yenagoa, Eleme, and Trans-Amadi, Niger Delta, Nigeria. The ecosystem risk of PAHs in soils from Khana was majorly, whereas due to Flt, Phe contributed less to the ecosystem risk of PAHs in soils from Uyo. Similarly, Ant contributed less to the ecosystem risk of PAHs in soils from Choba, Yenagoa, Eleme, and Trans-Amadi, Niger Delta, Nigeria.

4. Conclusion

Nap, Acy, Ace, Flu, Ant, Flt, and Pyr were the major causes of the ecological risk due to PAHs in soils from Uyo, Choba, Yenagoa, Eleme, and Trans-Amadi in the Niger Delta, Nigeria. The ecological risk assessment of PAHs derived from isomeric ratios suggested that the PAHs in soils and vegetation samples from Choba, Yenagoa, Eleme, and Trans-Amadi emanated from pyrogenic processes including traffic emissions, fossil fuels, and biomass combustion as well as petrogenic sources such as occasional spills of liquid petroleum fuels and discharges from automobile workshops.

From the risk quotient standpoint, some soil samples from the Niger Delta had RQ(NCs) values greater than 800, thus indicating that soil may require remediation to forestall ecohealth consequences.

Data Availability

The data used to support the findings of this study are included within the article.

Additional Points

Highlights. (i) There is a widespread PAH pollution in the Niger Delta, Nigeria, (ii) Nap, Acy, Ace, Flu, Ant, Flt, and Pyr pose major ecological risk of PAHs in soils, (iii) PAHs from the Niger Delta emanated from pyrogenic processes and petrogenic sources, and (iv) some soil in the Niger Delta may require remediation.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

EAO contributed to sampling and bench work. ANE and ILN contributed to manuscript writing and performed supervision. CF contributed to manuscript writing and to sourcing of fund. OEO conceptualised the study, performed supervision, and wrote the manuscript.

Acknowledgments

The study was partly funded under the subproject ‘Farm animals as sentinel for environmental health and food safety in Nigeria’ of the research project ‘The role of soil in the burden of noncommunicable diseases in Africa, from farm to lifestyle’ granted by the University of Tuscia, Viterbo, Italy, to the pan-African NGO NOODLES (Nutrition & food safety and wholesomeness-Prevention, education and research, www.noodlesonlus.org) (2017–2019). Esther Amaka Okoye received tuition waiver from World Bank.