The mean value for the samples ranged from 0.004 to 0.24 mg/L. Cadmium, cobalt, nickel, chromium, manganese, zinc, and lead were not significant and below the permissible limit set by the World Health Organization. Based on the calculation of the water pollution index, the mean values of water samples from all sampling locations fall within 29.02–134.78, indicating poor water quality. The order of decreasing HQing and HQderm values was as follows: cobalt > lead > chromium > cadmium > manganese > zinc > copper > nickel > iron, for both children and adults, respectively. Thirteen distinct polycyclic aromatic hydrocarbons were detected, with mean concentrations of 0.101 ± 0.09 for sample A, 0.1010 ± 0.09 for sample B, 0.0819 ± 0.09 for sample C, 0.01735 ± 0.02 for sample D, and 0.0118 ± 0.0 for sample E. The result showed total BaPeq (mg/kg) values ranging from 0.274 for sample A, 0.280 for sample B, 0.282 for sample C, 2.64 × 10−3 for sample D, and 1.22 × 10−5 for sample E. Elevated levels of fluorene, naphthalene, and dibenzyl (a–h) anthracene were observed in all the sampled points. Irrigation quality are as follows; sodium adsorption ratio (A: 6.31, B: 4.42), magnesium adsorption ratio (B: 52.27), %Na (A: 55.61), residual sodium carbonate (B: 53.36), permeability index (A: 50.31), soluble sodium percentage (B: 52.62), and Kelly's ratio (equilibrium).

  • The study focuses on Ohaji-Egbema, Nigeria's water pollution. It assesses pollutants and health risks comprehensively.

  • Findings offer tailored recommendations for remediation.

  • Emphasizes the urgent need for monitoring and remediation efforts.

  • Crucial insights for environmental protection and public health.

  • Preserving Ohaji-Egbema's natural resources and heritage.

Groundwater pollution frequently occurs as contaminants migrate from the soil into the ground level (Cho et al. 2023). If the topsoil is polluted, it can eventually contaminate the groundwater. The geology of the soil can also impact the potential of the underlying aquifer to become polluted (Yang 2021). However, water can be polluted by oil and other natural substances that are insoluble in water, and therefore they are normally found connected to or adsorbed into sediments (Nsabimana & Li 2022; Akakuru et al. 2023). The physicochemical parameters such as pH, conductivity, total dissolved solids (TDS), and temperature point toward water pollution (Ibe et al. 2020). If the groundwater is found to be contaminated, it is important to measure the levels of certain pollutants, such as heavy metals (HMs) and polycyclic aromatic hydrocarbons (PAHs), in order to determine the extent of the pollution and assess the potential health risks of the people using the water (Ibe et al. 2018, 2020). Among nature's abundant gifts, water holds a position of extreme importance, second only to air (Rohde et al. 2017; Edokpayi et al. 2018; Ihenetu et al. 2020; Reberski et al. 2022). When it comes to irrigation waters, salinity, sodicity, and ion toxicity represent substantial quality challenges (Mohamed & Hoda 2023). To ensure the continuous monitoring of soil and water resources for variability in salt concentrations, continuous environmental assessment becomes essential (Cho et al. 2023). These changes in salt concentration are critical since they can exceed the threshold levels at which certain salts, including , , Na+, Cl, Mg2+, and various trace elements, negatively affect plant growth (Mammeri et al. 2023). The water quality index (WQI) serves as a valuable tool for assessing the quality of irrigation water.

Crude oil exploration is a significant financial venture in Nigeria, resulting in the influx of PAHs and HMs into soils and water bodies due to oil spillage (Agbozu et al. 2017; Ambade et al. 2022; Akash et al. 2023). Some studies have evaluated the concentrations of HMs and irrigation studies (Duru et al. 2017; Akakuru et al. 2022, 2023) and PAHs (Huang et al. 2023) in groundwater. Amadi et al. (2020) used the metal pollution index and factor analysis to analyze the effects of HMs in groundwater in part of the Niger Delta and reported that nickel, cadmium, lead cobalt, chromium, copper, iron, mercury, zinc, arsenic, and manganese polluted the groundwater at various concentrations (Akakuru et al. 2022). For PAHs, Amadi et al. (2020) conducted a study on the diesel contamination of groundwater in Ring Road, Jos, Plateau State, Nigeria, and found that the distribution of PAHs in the groundwater was dominated by three-fused-ring PAHs, indicating a petrogenic origin. In terms of direct ingestion, phenanthrene, acenaphthylene, and fluorene exhibited the highest exposure doses of PAHs. Also, when considering dermal absorption, benzo(a)anthracene, benzo(b)fluoranthene, and benzo(k)fluoranthene were associated with the highest exposure doses. According to Ejiogu et al. (2019), PAHs concentration in water samples can be attributed to pyrogenic sources, and some are linked to petrogenic origins. Every stage of a crude oil spill has detrimental effects on the environment, and Ohaji-Egbema is no exception. Over 40 years of oil drilling in Nigeria, more than 6,000 spills have been documented, averaging 150 spills annually (Enyoh et al. 2018; Enyoh & Isiuku 2020). Oil spills have contaminated a significant portion of farmlands in the Niger Delta region, including Ohaji-Egbema, leading to the transformation of once-productive areas into wastelands (Chika et al. 2014).

The novelty of this research lies in the application of advanced, previously unused models for the comprehensive analysis of groundwater quality, health risks, and irrigation suitability in the Niger Delta region of Imo State, Nigeria. This approach not only provides a more detailed assessment of the environmental impact of crude oil exploration but also offers new insights into the interaction between various pollutants and local groundwater resources, which have not been explored in prior studies. This study, therefore, seeks to carry out physicochemical, pollution, health risk, and irrigation quality evaluation of some groundwater in these areas. This is especially crucial as the inhabitants of this area rely on groundwater for drinking, and domestic and agricultural activities. Second, the models employed for the analysis of the results represent a departure from previous approaches, as they have not been previously used in the study of groundwater in this area. The significance of this study is as follows: (1) the study places its focus on a region heavily impacted by crude oil exploration, namely the Niger Delta area in Imo State, Nigeria. (2) Groundwater is a vital resource in this region, serving as a source of drinking water for the local population and supporting agricultural activities. (3) This study contributes significantly to the understanding of environmental sustainability in regions impacted by industrial activities and highlights the importance of proactive environmental management to mitigate pollutant impacts and safeguard public health.

Study area – climate and geology

The research was conducted in the Ohaji-Egbema area, situated in the Niger Delta region of Nigeria, as shown in Figure 1. Ohaji-Egbema harbors an expansive landmass, covering approximately 958.010 km2. Geographically, Ohaji-Egbema occupies a strategic position within Imo State, as it is bounded by Owerri to the east, Oguta Local Government Area to the north, and the Ogba/Egbema/Ndoni area in Rivers State to the southwest. This geographical juxtaposition highlights the diverse facets of the southwestern region of Imo State and emphasizes the area's role as a significant crossroads within this part of Nigeria. The geographical coordinates of Ohaji-Egbema place it at a longitude of 6° 8′ 98″ and a latitude of 5° 34′ 53″. The elevation of this area stands at approximately 126.1m above sea level. In addition to its unique geographic and administrative attributes, Ohaji-Egbema is characterized by its proximity to various prominent oil companies. This region is widely acknowledged as a pivotal hub for oil drilling and petroleum-related activities within Nigeria. The presence of these oil companies shows Ohaji-Egbema's vital role in the nation's oil industry and its significance as an area of considerable strategic importance (Chika et al. 2014; Ihenetu et al. 2017).
Figure 1

Map of Ohaji-Egbema and its environs showing sampling points.

Figure 1

Map of Ohaji-Egbema and its environs showing sampling points.

Close modal

Sample collection

Sampling was done according to our previous research (Wang et al. 2020; Ihenetu et al. 2021b). Random sampling was used for sample collection, with 45 samples collected from each of the three communities: Mmahu-Egbema, Obosima, and Umuagwo, shown as A, B, and C, respectively. In each community, the samples were composited, resulting in a total of 135 samples across all areas. Sampling was done from February 2022 to November 2022. Control samples were also collected and labeled as D and E. As a reference, control samples were collected from Owerri Municipal, a local government area in Imo State, characterized by the absence of any recorded crude oil pollution. For the sample collection, 500 mL plastic bottles were used, carefully cleaned with 10% nitric acid, thoroughly rinsed with distilled water, and oven-dried to prevent bacterial degradation. Following collection, these bottles were promptly placed in an icebox and transported to the springboard laboratory for analysis. Upon arrival, they were stored in a refrigerator at 4°C before undergoing pretreatment analysis. For HM analysis, each sample was filtered using 0.45 μm membrane filters to remove particulate matter. The filtered samples were then acidified with concentrated nitric acid to a pH of less than 2 to preserve the metal ions and prevent precipitation. For PAHs analysis, the samples were extracted using liquid–liquid extraction with dichloromethane, followed by concentration under reduced pressure. The concentrated extracts were then purified using silica gel chromatography. The concentrations of HMs were determined using a Varian AA240 atomic absorption spectrophotometer, while PAHs were analyzed using gas chromatography-mass spectrometry. Calibration standards and quality control samples were run alongside the sample analysis to ensure accuracy and precision. Analytical methods for the analysis and health risk assessment methods are found in the Supplementary Material Text S1–S2.

Quality control

Quality control measures were rigorously applied during the experiments, which utilized high-grade analytical reagents from BDH Chemical Ltd, UK, and Sigma-Aldrich Chemie GmbH, Germany. Containers were thoroughly cleaned using detergents and deionized water, followed by overnight soaking in a solution of 10% HNO3 and 1% HCl, and then extensively rinsed with deionized water. Anion concentration determinations were carried out using reagents from Hanna Instruments with the Hanna Hi 83,200 instruments. Soil samples and reagents were precisely weighed using a Shimadzu AW 320 analytical balance. The extraction processes were performed with a Yamato Scientific RE201 rotary evaporator. To ensure reliable analytical results, stringent quality assurance protocols were observed, including the proper calibration of laboratory instruments before analyses. Accuracy and precision were verified through duplicate analyses of samples against standard reference materials, and blank samples were included in the analysis process to maintain data integrity.

Physicochemical parameters of samples

The physicochemical properties of the analyzed groundwater samples were compared with World Health Organization (WHO) standards and are presented in Table 1. The water temperature falls within the range of 20–30°C. Dissolved oxygen (DO) refers to the oxygen present in dissolved form in water bodies, and in this study, DO levels ranged from 4.32 to 5.04 mg/L, respectively. The electrical conductivity (EC) measurements in this study were all below the WHO standards, ranging from 0.98 to 32.14 mg/L. pH value is a crucial factor affecting the taste of water and can also indicate potential corrosion issues arising from the dissolution of metals such as copper, zinc, and cadmium in water (Itodo et al. 2018). The slight increases observed in both DO and pH during this study may suggest that both parameters are influenced by similar human activities in the study areas. These observations are consistent with findings at the Orji Mechanic Village, where the levels of DO fell below the standards set by the WHO (WWQA 2021).

Table 1

Mean physicochemical parameters and metal concentrations in groundwater at the areas of the study location

ParametersABCDcontrolEcontrolWHO
Temp. 25.31 25.19 25.22 26.31 26.08 20–30 
DO (mg/L) 5.02 4.89 5.04 4.32 4.48 10 
EC (μS/cm) 1.10 1.08 0.98 30.71 32.14 2500 
pH 7.04 7.02 7.03 6.21 5.91 6.5–8.5 
TDS (mg/L) 0.63 0.63 0.67 20.12 20.49 500 
Col (PCU) 8.22 8.63 8.44 8.03 8.21 15 
Hardness (mg/L) 20.32 20.22 20.24 27.42 26.81 200–500 
Alkalinity (mg/L) 16.57 16.56 16.49 5.41 25.62 120 
Acidity (mg/L) 17.81 17.56 17.47 15.41 15.55 6.5–9.2 
 3.97 3.96 3.84 10.31 10.82 50 
 4.52 4.56 4.88 2.31 2.82 1.0 
 2.46 2.25 2.62 1.42 2.03 250 
Ca (mg/L) 4.05 4.02 4.06 5.02 4.87 75 
Na (mg/L) 0.18 0.19 0.17 7.11 7.32 200 
K (mg/L) 3.27 3.26 3.27 4.83 4.77 20 
Fe (mg/L) 0.005 0.004 0.005 0.24 0.21 0.3 
Cu (mg/L) 0.016 0.014 0.015 0.01 0.00 2.0 
Cd (mg/L) 0.019 0.020 0.018 0.001 0.001 0.003 
Co (mg/L) 0.051 0.049 0.052 0.002 0.002 0.5 
Ni (mg/L) 0.016 0.015 0.015 0.0 0.0 0.02 
Cr (mg/L) 0.002 0.003 0.002 0.05 0.07 0.05 
Mn (mg/L) 0.063 0.062 0.062 0.05 0.05 0.4 
Zn (mg/L) 0.556 0.562 0.554 0.0 0.0 3.0 
Pb (mg/L) 0.004 0.006 0.003 0.00 0.00 0.01 
ParametersABCDcontrolEcontrolWHO
Temp. 25.31 25.19 25.22 26.31 26.08 20–30 
DO (mg/L) 5.02 4.89 5.04 4.32 4.48 10 
EC (μS/cm) 1.10 1.08 0.98 30.71 32.14 2500 
pH 7.04 7.02 7.03 6.21 5.91 6.5–8.5 
TDS (mg/L) 0.63 0.63 0.67 20.12 20.49 500 
Col (PCU) 8.22 8.63 8.44 8.03 8.21 15 
Hardness (mg/L) 20.32 20.22 20.24 27.42 26.81 200–500 
Alkalinity (mg/L) 16.57 16.56 16.49 5.41 25.62 120 
Acidity (mg/L) 17.81 17.56 17.47 15.41 15.55 6.5–9.2 
 3.97 3.96 3.84 10.31 10.82 50 
 4.52 4.56 4.88 2.31 2.82 1.0 
 2.46 2.25 2.62 1.42 2.03 250 
Ca (mg/L) 4.05 4.02 4.06 5.02 4.87 75 
Na (mg/L) 0.18 0.19 0.17 7.11 7.32 200 
K (mg/L) 3.27 3.26 3.27 4.83 4.77 20 
Fe (mg/L) 0.005 0.004 0.005 0.24 0.21 0.3 
Cu (mg/L) 0.016 0.014 0.015 0.01 0.00 2.0 
Cd (mg/L) 0.019 0.020 0.018 0.001 0.001 0.003 
Co (mg/L) 0.051 0.049 0.052 0.002 0.002 0.5 
Ni (mg/L) 0.016 0.015 0.015 0.0 0.0 0.02 
Cr (mg/L) 0.002 0.003 0.002 0.05 0.07 0.05 
Mn (mg/L) 0.063 0.062 0.062 0.05 0.05 0.4 
Zn (mg/L) 0.556 0.562 0.554 0.0 0.0 3.0 
Pb (mg/L) 0.004 0.006 0.003 0.00 0.00 0.01 

TDS values ranged from 0.62 to 20.49 mg/L, with higher concentrations recorded in the control samples. Watercolor plays a significant role in its appearance, taste, and overall drinkability. The color of water samples at all sampling locations remained below the permissible limit, ranging from 8.03 to 8.63 Platinum-Cobalt Unit (PCU). Alkalinity measures water's acid-neutralizing capacity, acting as a buffer against abrupt shifts in pH, and its assessment depends on the specific chemical composition of the sample. Groundwater alkalinity in this study varied from 5.41 to 25.62 mg/L, contrasting with the WHO standard of 120 mg/L. Hardness, associated with high mineral content in water, occurs when water filters through calcium (Ca)- and magnesium-rich formations such as limestone, chalk, or gypsum. Water hardness ranged from 20.22 to 27.42 mg/L, indicating no pollution of groundwater with regards to hardness. Water acidity, reflecting its capacity to react with a strong base up to a specific pH value, was observed to be between 15.41 and 17.81 mg/L, consistent with previous research by Ward et al. (2020) and Birhanu et al. (2021).

, a stable form of nitrogen, was consistently detected in the groundwater samples. Nitrate levels were found to be lower than the WHO Standard for safe drinking water, ranging from 3.84 to 10.82 mg/L. Phosphorus was a crucial nutrient responsible for the fertility and productivity of fishponds, particularly in supporting plankton growth (Francis et al. 2021; Hu et al. 2014). levels in this study ranged from 2.31 to 4.56 mg/L. The higher levels observed may be attributed to the frequent use of phosphate-based fertilizers by farmers in the study area. values, ranging from 1.42 to 2.46 mg/L, were observed in this study, and are below the WHO Standards for drinking water and domestic water use. Similar findings regarding sulfate levels were recorded in a study conducted in the Niger Delta area in Nigeria (Akakuru et al. 2022). This result agrees with the outcome of the result obtained from sub-Saharan Africa (Banks et al. 2021). Fe levels observed in this study, ranging from 0.004 to 0.24 mg/L, are within the range of the WHO Standards. The presence of Fe in all the samples could be attributed to the use of Fe coagulants in the area. Cu, an essential nutrient, and potential drinking water contaminant (Joshua et al. 2016), was detected at levels below the WHO Standards for drinking water and domestic use in this study, with concentrations ranging from 0.00 to 0.015 mg/L. The low Cu levels align with findings from a study conducted in European groundwater (Bunting et al. 2021) and a similar study on Owerri municipality and environs, southeastern Nigeria (Ibe et al. 2020). Cd, Co, Ni, Ch, Mn, Zn, and Pb were all found to be below significant levels and within the permissible limits set by WHO.

A substantial positive correlation (r > 0.5) was observed between some of the metals, and anions parameters, and the relationship is shown in Figure 2. However, significant positive relations in this study are between /Fe (0.9904), /Cr (0.9853), / (0.9223), /Cu (0.7344), /Cd (0.9712), /Co (0.9832), /Ni (0.971), /Mn (0.9799), /Pb (0.8586), /Cd (0.8194), /Co (0.8522), /Ni (0.8417), /Mn (0.8414), /Zn (0.8389), /Pb (0.6505), Fe/Cr (0.9528), Cu/Cd (0.8292), Cu/Co (0.8374), Cu/Ni (0.8388), Cu/Mn (0.8393), Cu/Zn (0.8344), Cu/Pb (0.7310), Cd/Co (0.9937), Cd/Ni (0.9962), Cd/Mn (0.9956), Cd/Zn (0.9980), Cd/Pb (0.9369), Co/Ni (0.9983), Co/Mn (0.9978), Co/Zn (0.9987), Co/Pb (0.8926), Ni/Mn (0.9999), Ni/Zn (0.9986), Ni/Pb (0.9053), Mn/Zn (0.9979), Mn/Pb (0.9038), and Zn/Pb (0.9141). Once the correlation is seen as positive, the source of contamination of the positively linked metals is alike, while a negative correlation recommends dissimilar/different bases of contamination. For instance, there is a strong negative correlation between and Fe with a coefficient of −0.9904. Additionally, exhibits negative correlations with Cu, Cd, Co, Ni, Zn, and Pb, with correlation coefficients ranging from −0.8607 to −0.9985. These negative correlations imply that as the concentration of increases, the concentrations of these other parameters tend to decrease. This inverse relationship is essential to consider when assessing the complex dynamics of water quality in the study area. Understanding these correlations helps in managing groundwater resources effectively and addressing environmental concerns.
Figure 2

Correlation coefficient matrix between water quality parameters in the ground water.

Figure 2

Correlation coefficient matrix between water quality parameters in the ground water.

Close modal

Hierarchical cluster analysis (HCA)

Additionally, we performed HCA to identify groupings of physicochemical characteristics based on their square Euclidian distance. The cluster plots for physicochemical parameters in the samples are presented in Figure 3. For physicochemical parameters, three groups were identified. In group 1, the combination included all parameters except for pH and EC, and in another group then TDS, temperature, and alkalinity in the third group. Physicochemical parameters such as Ca, potassium (K), DO, phosphate (), sulfate (), pH, color, nitrate (), sodium (Na), EC, TDS, temperature, hardness, alkalinity, and acidity are analyzed. The rescaled distance cluster combined on the x-axis represents the degree of similarity or dissimilarity between parameters; the smaller the distance, the more similar the parameters are. For example, parameters that join at a lower distance are more closely associated, indicating that they may be affected by similar sources or processes in the groundwater. Similarly, in the HMs parameters, the combination includes all parameters except for Cr and Fe in group 2 and Zn in group 3. The clustering of all metals similarly indicates that their source(s) are common. For the second dendrogram, the HMs included Cd, Ni, Cu, Pb, Co, Cr, Mn, Fe, and Zn. Similarly, metals clustered together at lower distances may have a common source or similar behavior in the groundwater environment.
Figure 3

HCA for physicochemical in groundwater samples.

Figure 3

HCA for physicochemical in groundwater samples.

Close modal

Chemometric analysis for HMs

The contamination factor (CF) and pollution load index (PLI) for the groundwater samples are shown in Figure 4, offering valuable insights into the environmental quality of the sampled locations. Each set of bars in the graph likely corresponds to different HMs, and the height of the bars indicates the CF for each metal in various locations or samples. A higher bar would suggest a greater degree of contamination by that specific metal. The PLI, on the other hand, is a tool used to give a cumulative indication of the overall level of HM pollution. It is typically calculated as the nth root of the product of the CFs for n metals. A PLI value greater than 1 indicates pollution, whereas a value of 1 or less suggests no significant pollution. As the PLI serves as a comprehensive indicator of contamination, it is worth noting that all sampling points exhibit relatively low PLI values, falling within the range of 0.094–0.33. However, a more detailed examination of individual contaminants reveals intriguing patterns. Sampling points A, B, and C, which are situated in areas associated with crude oil activities, exhibit a moderate level of contamination primarily attributed to Fe levels. Conversely, the remaining sampling points, especially those in non-crude oil areas such as D and E, display higher contamination levels with regard to Fe. This discrepancy in contamination patterns highlights the complex interplay between geological and anthropogenic factors in the study area. Notably, sampling points D and E show exceptions, displaying minimal or no detectable contamination by Ni. These variations in contamination levels emphasize the role of local factors, industrial activities, and geographical differences in shaping groundwater quality.
Figure 4

Contamination factor and PLI of HMs in ground waters samples.

Figure 4

Contamination factor and PLI of HMs in ground waters samples.

Close modal

Water quality index

The WQI values reveal a spectrum of water quality ranging from 29.02 to 134.78, reflecting the diverse environmental conditions in this study area. Notably, these values collectively indicate that the overall water quality is predominantly poor, showing the pressing need for heightened environmental management efforts. Specifically, the presence of poor water quality is evident in samples A, B, and C, which are situated in regions associated with various human activities and potential sources of pollution. These activities, coupled with the possibility of oil spillage in proximity to the sampling points, contribute to the degradation of groundwater quality in these areas. The calculated WQI values for these locations fall within the higher range, emphasizing the urgency of implementing targeted mitigation measures to address the contamination. Samples D and E provide a stark contrast, boasting excellent water quality with WQI values ranging from 29.02 to 33.47. To summarize this point, the computation of the WQI highlights the varying degrees of water quality across the sampled locations. While samples A, B, and C exhibit poor water quality, attributed to a range of anthropogenic and environmental factors, samples D and E shine as beacons of excellence, showcasing the availability of uncontaminated groundwater resources. This result is shown in Table 2.

Table 2

WQI values of the sampling points

Parameterqi
S1
siwiWiABCDcontrolEcontrolABCDcontrolEcontrol
Fe 0.3 0.103 1.67 1.33 1.67 80.0 70.0 0.17 0.14 0.17 8.24 7.21 
Cu 2.0 0.068 0.80 0.70 0.75 0.50 0.00 0.05 0.05 0.05 0.03 0.03 
Cd 0.003 0.172 633.3 666.6 600.0 33.33 33.33 108.9 114.6 103.2 5.73 5.73 
Co 0.5 0.103 10.20 9.80 10.40 0.40 0.40 1.05 1.01 1.07 0.04 0.04 
Ni 0.02 0.103 80.00 75.00 75.00 0.00 0.00 8.24 7.72 7.72 0.00 0.00 
Cr 0.05 0.137 4.00 6.00 4.00 100 140 0.55 0.82 0.54 13.7 19.18 
Mn 0.4 0.103 15.75 15.50 15.50 12.50 12.50 1.62 1.59 1.59 1.28 1.28 
Zn 3.0 0.034 18.53 18.73 18.46 0.00 0.00 0.62 0.63 0.62 0.00 0.00 
Pb 0.01 0.137 40.00 60.00 30.00 0.00 0.00 5.48 8.22 4.11 0.00 0.00 
 Σwi 29           
   WQI      126.68 134.78 119.07 29.02 33.47 
Parameterqi
S1
siwiWiABCDcontrolEcontrolABCDcontrolEcontrol
Fe 0.3 0.103 1.67 1.33 1.67 80.0 70.0 0.17 0.14 0.17 8.24 7.21 
Cu 2.0 0.068 0.80 0.70 0.75 0.50 0.00 0.05 0.05 0.05 0.03 0.03 
Cd 0.003 0.172 633.3 666.6 600.0 33.33 33.33 108.9 114.6 103.2 5.73 5.73 
Co 0.5 0.103 10.20 9.80 10.40 0.40 0.40 1.05 1.01 1.07 0.04 0.04 
Ni 0.02 0.103 80.00 75.00 75.00 0.00 0.00 8.24 7.72 7.72 0.00 0.00 
Cr 0.05 0.137 4.00 6.00 4.00 100 140 0.55 0.82 0.54 13.7 19.18 
Mn 0.4 0.103 15.75 15.50 15.50 12.50 12.50 1.62 1.59 1.59 1.28 1.28 
Zn 3.0 0.034 18.53 18.73 18.46 0.00 0.00 0.62 0.63 0.62 0.00 0.00 
Pb 0.01 0.137 40.00 60.00 30.00 0.00 0.00 5.48 8.22 4.11 0.00 0.00 
 Σwi 29           
   WQI      126.68 134.78 119.07 29.02 33.47 

Assessment of health risk for HMs

Dermal and ingestion exposure, hazard quotient (HQ), hazard indices (HI)

Table 3 presents the HQ values for potential non-carcinogenic risk (CR) associated with each HM in both adults and children, as well as the cumulative HI. The HQ values are indicative of the relative health risk posed by individual HMs through dermal absorption (HQderm) and oral ingestion (HQing). For Fe, the HQ values ranged from a minimum of 7.30 × 10−9 (in adults through dermal absorption) to a maximum of 4.11 × 10−5 (in adults through oral ingestion) which indicate low hazard. Cu showed very low HQ values, often reaching zero, indicating minimal non-CR for both adults and children. Cd displayed slightly higher HQ values but remained within an acceptable range, suggesting limited non-CR. Co exhibited higher HQ values, especially in adults, implying a somewhat higher non-CR. Ni had extremely low HQ values, often at zero, indicating minimal non-CR. Cr displayed moderate HQ values, indicating a relatively higher potential non-CR. Mn exhibited HQ values within an acceptable range, signifying limited non-CR. Zn consistently showed very low HQ values, often at zero, suggesting minimal non-CR. Pb displayed HQ values within an acceptable range, indicating limited non-CR. It is important to highlight the decreasing order of HQ values for both children and adults, which ranged from Co > Pb > Cr > Cd > Mn > Zn > Co > Ni > Fe. This ranking provides insights into the relative levels of risk associated with each trace metal, with Co presenting the highest potential risk and Fe the lowest, for both age groups. The cumulative HI for all HMs in both adults and children was below one, suggesting that the combined non-CR posed by these metals was minimal and unlikely to result in significant adverse health effects in the studied populations.

Table 3

HQ for potential non-CR (HQ) and cumulative HI for each heavy metal in both adult and child populations

MetalsStatistical parameterHQderm (adult)HQderm (children)HQing (adult)HQing (children)
Fe Maximum 4.41 × 10−7 7.50 × 10−7 1.07 × 10−5 4.11 × 10−5 
 Minimum 7.30 × 10−9 1.25 × 10−8 1.78 × 10−7 6.85 × 10−7 
Cu Maximum 5.13 × 10−7 8.80 × 10−7 1.25 × 10−5 4.80 × 10−5 
 Minimum 0.00 0.00 0.00 0.00 
Cd Maximum 1.03 × 10−5 1.60 × 10−5 2.51 × 10−2 9.60 × 10−2 
 Minimum 5.14 × 10−7 8.80 × 10−7 1.25 × 10−3 4.80 × 10−3 
Co Maximum 1.07 × 10−3 1.83 × 10−3 8.15 × 10−2 3.12 × 10−1 
 Minimum 4.10 × 10−5 7.04 × 10−5 3.10 × 10−3 1.20 × 10−2 
Ni Maximum 7.61 × 10−8 1.30 × 10−7 2.57 × 10−5 9.60 × 10−5 
 Minimum 0.00 0.00 0.00 0.00 
Cr Maximum 2.40 × 10−5 4.11 × 10−5 7.33 × 10−4 2.80 × 10−3 
 Minimum 6.85 × 10−7 1.17 × 10−6 2.09 × 10−5 8.00 × 10−5 
Mn Maximum 1.68 × 10−5 2.88 × 10−5 5.87 × 10−5 3.15 × 10−4 
 Minimum 1.33 × 10−7 2.29 × 10−5 6.54 × 10−5 2.50 × 10−4 
Zn Maximum 7.23 × 10−7 1.23 × 10−6 5.86 × 10−5 2.24 × 10−4 
 Minimum 0.00 0.00 0.00 0.00 
Pb Maximum 1.46 × 10−5 2.50 × 10−5 1.34 × 10−4 5.14 × 10−4 
 Minimum 0.00 0.00 0.00 0.00 
HI Maximum 1.14 × 10−3 1.94 × 10−3 1.10 × 10−1 4.12 × 10−1 
 Minimum 4.23 × 10−5 9.53 × 10−5 4.43 × 10−3 1.71 × 10−2 
MetalsStatistical parameterHQderm (adult)HQderm (children)HQing (adult)HQing (children)
Fe Maximum 4.41 × 10−7 7.50 × 10−7 1.07 × 10−5 4.11 × 10−5 
 Minimum 7.30 × 10−9 1.25 × 10−8 1.78 × 10−7 6.85 × 10−7 
Cu Maximum 5.13 × 10−7 8.80 × 10−7 1.25 × 10−5 4.80 × 10−5 
 Minimum 0.00 0.00 0.00 0.00 
Cd Maximum 1.03 × 10−5 1.60 × 10−5 2.51 × 10−2 9.60 × 10−2 
 Minimum 5.14 × 10−7 8.80 × 10−7 1.25 × 10−3 4.80 × 10−3 
Co Maximum 1.07 × 10−3 1.83 × 10−3 8.15 × 10−2 3.12 × 10−1 
 Minimum 4.10 × 10−5 7.04 × 10−5 3.10 × 10−3 1.20 × 10−2 
Ni Maximum 7.61 × 10−8 1.30 × 10−7 2.57 × 10−5 9.60 × 10−5 
 Minimum 0.00 0.00 0.00 0.00 
Cr Maximum 2.40 × 10−5 4.11 × 10−5 7.33 × 10−4 2.80 × 10−3 
 Minimum 6.85 × 10−7 1.17 × 10−6 2.09 × 10−5 8.00 × 10−5 
Mn Maximum 1.68 × 10−5 2.88 × 10−5 5.87 × 10−5 3.15 × 10−4 
 Minimum 1.33 × 10−7 2.29 × 10−5 6.54 × 10−5 2.50 × 10−4 
Zn Maximum 7.23 × 10−7 1.23 × 10−6 5.86 × 10−5 2.24 × 10−4 
 Minimum 0.00 0.00 0.00 0.00 
Pb Maximum 1.46 × 10−5 2.50 × 10−5 1.34 × 10−4 5.14 × 10−4 
 Minimum 0.00 0.00 0.00 0.00 
HI Maximum 1.14 × 10−3 1.94 × 10−3 1.10 × 10−1 4.12 × 10−1 
 Minimum 4.23 × 10−5 9.53 × 10−5 4.43 × 10−3 1.71 × 10−2 

Chronic daily intake (CDI) and CR

The maximum CDI values for the selected metals ranged between 8.0 × 10−3 and 2.1 × 10−3 in adults, while the children's index was 4.1 × 10−2 to 3.4 × 10−3, separately. The CDI indices for HMs over the study duration for both age groups followed the sequence of Zn > Fe > Mn > Cr > Co > Cd > Ni > Cu > Pb, as seen in Figure 5. This is a sign that the groundwater has less important health effects on both adults and children via the pathways (Adimalla and Qian 2019), however, regular checks are needed to circumvent the accumulation of HMs that may pose any health problems, especially in children.
Figure 5

Chronic risk assessment (CDIing) of HMs for both adults and children.

Figure 5

Chronic risk assessment (CDIing) of HMs for both adults and children.

Close modal

The CR of Pb on the ground samples was calculated for both adults and children due to the unavailability of slope factors for other HMs in the literature. In the study, the average levels of CRing for Pb show 0.81−5 and 2.21−5 for adults and children. Under extreme regulatory programs, the CR values between 10−6 and 10−4 for an individual show potential risk; therefore, the results in this study showed that the level of Pb in the groundwater does not pose a CR to both adults and children. It is important to consider the long-term effects of chronic exposure, particularly given the cumulative nature of these contaminants. Although, the current CDI and CR values for the studied metals suggest that the immediate health risks are low, especially regarding Pb, the possibility of bioaccumulation over prolonged periods cannot be ignored. This is mostly concerning for children, who are more vulnerable to the adverse effects of HMs due to their developing systems and higher relative intake of water (Ibe et al. 2023). The absence of slope factors for other HMs in the literature should not be misconstrued as an indication of safety; rather, it highlights a critical gap in our understanding that warrants further investigation.

Characteristics of PAH concentrations and analysis

Table 4 presents the mean levels of the concentration of PAHs in the groundwater samples from different sampling points, including A, B, C, D (Control), and E (Control). Naphthalene, a common PAH, was detected in samples A, B, and C, with concentrations ranging from 0.1838 to 0.1975 mg/kg. Dibenzyl (a–h) anthracene was also present in these samples at similar concentrations. Fluorene was detected in samples A, B, and C, with concentrations around 0.25 mg/kg, while sample E had a minor presence of 0.0117 mg/kg. Acenaphthene and 1–2 benzanthracene were found in trace amounts in samples A, B, and C. Benzo (a) pyrene, a well-known carcinogenic PAH, was detected in all samples, with concentrations ranging from 0.0515 to 0.0536 mg/kg. Chrysene, another potentially harmful PAH, was present in all samples, with variable concentrations. Benzo (g-h-i) perylene and benzo (k) fluoranthene were also detected in samples A, B, and C. Acenaphthylene was detected in sample B, while pyrene and phenanthrene were found in samples C and D. Flouranthene was detected in sample E. The mean PAH levels observed in the study are in the order of Flu > Nap > DaA > BgP > Chr > Ace > BzA > Bkf > Acy > Phe > Flo > Pyr. Overall, the mean concentrations of PAHs varied among the sampling points, with samples A, B, and C generally showing higher levels of several PAH compounds compared to the control samples D and E. These findings indicate the potential presence of PAH contamination in the groundwater of the study area, with certain compounds at levels of concern, particularly benzo (a) pyrene. Monitoring and further investigation are warranted to assess the environmental and health implications of these PAHs in the groundwater. Elevated levels of fluorene, naphthalene, and dibenzyl (a–h) anthracene were detected in nearly all the sampling points. Additionally, the occurrence of pyrene (Pyr), phenanthrene (Phe), fluorene (Flu), and acenaphthylene (Acy) was observed to be relatively low.

Table 4

Mean levels of the concentration (mg/kg) of PAHs in the groundwater samples

PAH (mg/kg)ABCDcontrolEcontrol
Naphthalene 0.1838 0.1975 0.1845 Not Detected (ND) ND 
Dibenzyl (a–h) anthracene 0.2177 0.2235 0.2236 ND ND 
Fluorene 0.2641 0.2545 0.2514 ND 0.0117 
Acenaphthene 0.0343 0.0336 0.0359 ND ND 
1–2 Benzanthracene 0.0286 0.0198 0.0209 ND ND 
Benzo (a) pyrene 0.0515 0.0524 0.0536 ND ND 
Chrysene 0.0458 0.0447 0.0379 0.0425 ND 
Benzo (g-h-i) perylene 0.0728 0.0738 0.0784 ND ND 
Benzo (k) fluoranthene 0.0106 0.0110 0.0124 0.0222 ND 
Acenaphthylene ND 0.0010 ND 0.0035 ND 
Pyrene ND ND 0.0012 ND ND 
Phenanthrene ND ND 0.0021 0.0012 ND 
Flouranthene ND ND ND ND 0.001 
Mean ± Std 0.1010 ± 0.09 0.1010 ± 0.09 0.0819 ± 0.09 0.01735 ± 0.02 0.0118 ± 0.0 
PAH (mg/kg)ABCDcontrolEcontrol
Naphthalene 0.1838 0.1975 0.1845 Not Detected (ND) ND 
Dibenzyl (a–h) anthracene 0.2177 0.2235 0.2236 ND ND 
Fluorene 0.2641 0.2545 0.2514 ND 0.0117 
Acenaphthene 0.0343 0.0336 0.0359 ND ND 
1–2 Benzanthracene 0.0286 0.0198 0.0209 ND ND 
Benzo (a) pyrene 0.0515 0.0524 0.0536 ND ND 
Chrysene 0.0458 0.0447 0.0379 0.0425 ND 
Benzo (g-h-i) perylene 0.0728 0.0738 0.0784 ND ND 
Benzo (k) fluoranthene 0.0106 0.0110 0.0124 0.0222 ND 
Acenaphthylene ND 0.0010 ND 0.0035 ND 
Pyrene ND ND 0.0012 ND ND 
Phenanthrene ND ND 0.0021 0.0012 ND 
Flouranthene ND ND ND ND 0.001 
Mean ± Std 0.1010 ± 0.09 0.1010 ± 0.09 0.0819 ± 0.09 0.01735 ± 0.02 0.0118 ± 0.0 

The PAHs that were recorded were divided into three categories based on their molecular weight and distribution (Figure 6). These categories are low molecular weight (LMW), medium molecular weight (MMW), and high molecular weight (HMW). LMW PAHs tend to be formed through low-temperature processes such as burning biomass, while HMW PAHs are produced during high-temperature reactions such as burning fuel. The distribution of PAHs based on molecular weights is presented in Figure 6. The MMW PAHs had the highest total concentration in all samples (55.96%). The isomeric ratio of LMW to HMW has been used to identify the source of PAHs. If the LMW/HMW ratio is less than 1, it is considered pyrogenic, but when greater than 1, it is considered petrogenic (Akakuru et al. 2023). The LMW/HMW ratio obtained was 1.05, indicating the source of the PAHs was petrogenic. Petrogenic PAHs are those that originate from the breakdown of fossil fuels and could be linked to the study area, which is an oil-producing region in Nigeria. Most PAHs are in the HMW category, indicating larger, less volatile molecules with more than five aromatic rings, which are generally more persistent in the environment and can be more carcinogenic. This distribution pattern is critical for understanding environmental exposure risks and the development of appropriate remediation strategies, as the molecular weight of PAHs influences their environmental behavior and toxicity. Naphthalene provides an indicator of refined petroleum release or crude oil leakage (Liu et al. 2017).
Figure 6

Distribution of PAHs by weights in all samples. LMW, MMW, and HMW.

Figure 6

Distribution of PAHs by weights in all samples. LMW, MMW, and HMW.

Close modal

Carcinogenic potency of the PAHs

The results revealed total BaPeq (mg/kg) values spanning from = 2.74−1 for A, 2.80−1 for B, 2.82−1 for C, 2.64−3 for D, and 1.22−5 for sample E. It is worth noting that atmospheric sedimentation could contribute to the presence of PAH concentrations in seemingly unpolluted areas due to economic and other human activities (Bhattacharya et al. 2020; Olayinka et al. 2018). Furthermore, the elevated concentrations of PAHs in the primary study sites may be linked to recurrent local flooding events. This region experiences frequent heavy rainfall, leading to water overflow and extensive flooding, which results in significant leaching and erosion of topsoil particles containing accumulated PAHs and crude oil spillage (Ugochukwu & Ochonogor 2018). Consequently, these PAHs are distributed to floodplains and other low-lying areas. The carcinogenic potency is presented in Table 5. LMW PAHs such as naphthalene exhibit minimal carcinogenic potential (Toxic Equivalence Quality Factor (TEQF) of 0.001), while higher molecular weight PAHs such as dibenz(a,h)anthracene and benzo(a)pyrene, with five or more rings, show significantly higher potencies (TEQF of 1.00). Intermediate weights such as fluorene, acenaphthene, and chrysene have moderate TEQFs, indicating a corresponding CR. The BaPeq values, which aggregate the carcinogenic potential of all PAHs present, suggest that samples A, B, and C have a higher CR compared to the negligible levels in control samples D and E. This indicates potential environmental contamination in the study area, necessitating careful monitoring and possibly intervention to mitigate the associated health risks. Unlike previous studies of soil PAHs in the metropolitan area, the present investigation identified a high-risk level in the study area (Zheng et al. 2018; Xu et al. 2021).

Table 5

Carcinogenic potency of the PAHs in water samples in the study area

PAHNo of ringsTEQFABCD(Control)E(Control)
Naphthalene 0.001 1.83 × 10−4 1.97 × 10−4 1.84 × 10−4 0.00 0.00 
Dibenzyl (a–h) anthracene 1.00 2.17 × 10−1 2.23 × 10−1 2.24 × 10−1 0.00 0.00 
Fluorene 0.001 2.64 × 10−4 2.54 × 10−4 2.51 × 10−4 0.00 1.12 × 10−5 
Acenaphthene 0.001 3.43 × 10−5 3.36 × 10−5 3.59 × 10−5 0.00 0.00 
1–2 Benzanthracene 0.10 2.86 × 10−3 1.98 × 10−3 2.09 × 10−3 0.00 0.00 
Benzo (a) pyrene 1.00 5.15 × 10−2 5.24 × 10−2 5.36 × 10−2 0.00 0.00 
Chrysene 0.01 4.58 × 10−4 4.47 × 10−4 3.79 × 10−4 4.25 × 10−4 0.00 
Benzo (g-h-i) perylene 0.01 7.28 × 10−4 7.38 × 10−4 7.84 × 10−4 0.00 0.00 
Benzo (k) fluoranthene 0.1 1.06 × 10−3 1.10 × 10−3 1.24 × 10−3 2.22 × 10−3 0.00 
Acenaphthylene 0.001 0.00 1.00 × 10−6 0.00 3.50 × 10−6 0.00 
Pyrene 0.001 0.00 0.00 1.20 × 10−6 0.00 0.00 
Phenanthrene 0.001 0.00 0.00 2.10 × 10−6 1.20 × 10−6 0.00 
Flouranthene 0.001 0.00 0.00 0.00 0.00 1.00 × 10−6 
BaPeq (mg/kg)   2.74 × 10−1 2.80 × 10−1 2.82 × 10−1 2.64 × 10−3 1.22 × 10−5 
PAHNo of ringsTEQFABCD(Control)E(Control)
Naphthalene 0.001 1.83 × 10−4 1.97 × 10−4 1.84 × 10−4 0.00 0.00 
Dibenzyl (a–h) anthracene 1.00 2.17 × 10−1 2.23 × 10−1 2.24 × 10−1 0.00 0.00 
Fluorene 0.001 2.64 × 10−4 2.54 × 10−4 2.51 × 10−4 0.00 1.12 × 10−5 
Acenaphthene 0.001 3.43 × 10−5 3.36 × 10−5 3.59 × 10−5 0.00 0.00 
1–2 Benzanthracene 0.10 2.86 × 10−3 1.98 × 10−3 2.09 × 10−3 0.00 0.00 
Benzo (a) pyrene 1.00 5.15 × 10−2 5.24 × 10−2 5.36 × 10−2 0.00 0.00 
Chrysene 0.01 4.58 × 10−4 4.47 × 10−4 3.79 × 10−4 4.25 × 10−4 0.00 
Benzo (g-h-i) perylene 0.01 7.28 × 10−4 7.38 × 10−4 7.84 × 10−4 0.00 0.00 
Benzo (k) fluoranthene 0.1 1.06 × 10−3 1.10 × 10−3 1.24 × 10−3 2.22 × 10−3 0.00 
Acenaphthylene 0.001 0.00 1.00 × 10−6 0.00 3.50 × 10−6 0.00 
Pyrene 0.001 0.00 0.00 1.20 × 10−6 0.00 0.00 
Phenanthrene 0.001 0.00 0.00 2.10 × 10−6 1.20 × 10−6 0.00 
Flouranthene 0.001 0.00 0.00 0.00 0.00 1.00 × 10−6 
BaPeq (mg/kg)   2.74 × 10−1 2.80 × 10−1 2.82 × 10−1 2.64 × 10−3 1.22 × 10−5 

The higher BaPeq levels in samples A, B, and C, particularly those with higher concentrations of HMW PAHs such as dibenz(a,h)anthracene and benzo(a)pyrene, show the significant carcinogenic threat these areas face. The disparity between these samples and the lower BaPeq values in control samples D and E highlights the extent of limited contamination, likely aggravated by the unique environmental conditions of the region, such as recurrent flooding and oil spillage. Moreover, the presence of these contaminants in groundwater, a critical resource for both drinking and agricultural activities, raises serious public health concerns, as seen in Table 5. The results suggest that the environmental dynamics of the area, combined with anthropogenic activities, are contributing to a heightened CR that warrants immediate attention.

Irrigation quality of the investigated water sources

Sodium adsorption ratio (SAR)

Sample A exhibits a moderate SAR of 6.31. This SAR value suggests that the water contains a notable amount of Na ions relative to Ca and magnesium ions. Such conditions can potentially lead to the modification of soils, affecting soil structure and crop growth. Sample B, while still indicating a Na hazard, exhibits a relatively lower SAR of 4.42 compared to Sample A. This suggests a reduced risk of soil sodification when using this water for irrigation. In contrast, Samples C, D, and E have notably lower SAR values (2.36, 0.65, and 0.66, respectively). These lower SAR values imply a lower Na hazard in these samples, making them more suitable for irrigation without the significant risk of soil sodification.

Magnesium adsorption ratio (MAR)

Sample B presents a strikingly high MAR of 52.27. This MAR value indicates a considerable presence of magnesium in the water, which can adversely affect soil structure and permeability. Farmers using this water should be cautious about potential soil issues. In contrast, the other samples have substantially lower MAR values, signifying a diminished concern regarding excessive magnesium content in the water.

Sodium percentage (%Na)

Sample A records the highest %Na at 55.61. This figure denotes that more than half of the cations in the water are Na ions. High %Na values can have a pronounced impact on soil structure, potentially leading to soil dispersion and reduced infiltration rates.

Residual sodium carbonate (RSC)

Sample B displays an extremely negative RSC value of −53.36. This strongly negative RSC indicates a significant potential for soil alkalinity, which can adversely affect soil health and crop growth. Conversely, the other samples exhibit less negative or even positive RSC values, implying a lower risk of soil alkalinity.

Sodium permeability index (PI)

Sample A's high PI of 50.31 suggests that this water may not be suitable for sodium-sensitive soils. It implies that water with such characteristics could potentially compromise soil structure and restrict water infiltration. In contrast, the other samples possess lower PI values, indicating a better suitability for a broader range of soils.

Soluble sodium percentage (SSP)

Sample B stands out with the highest SSP of 52.62, suggesting a substantial Na hazard concerning all cations present in the water. Farmers should exercise caution when considering this water source for irrigation. The other samples present lower SSP values, indicating a comparatively lower Na hazard among all cations.

Kelly's ratio (KR)

All samples exhibit KR values close to 1. This suggests a balanced proportion of Na and Ca ions, indicating that the water does not favor one cation over the other in terms of dominance. This equilibrium in cationic composition can be advantageous for irrigation purposes, as it minimizes the risk of soil dispersion and other structural issues associated with Na dominance. All these results are seen in Figure 7.
Figure 7

Results of WDIs for irrigation.

Figure 7

Results of WDIs for irrigation.

Close modal

In this study on ecological risk assessment conducted in Imo State, Nigeria, we examined the presence of HMs and PAHs in groundwater from areas impacted by crude oil pollution. Our investigation revealed the widespread presence of these pollutants across all sampled locations, including control samples, although the control samples exhibited limited PAH presence. Our findings showed that high- fertilizers in farmlands have been used in the sampling points, which contributed to the elevation of levels that is higher than the WHO Standards. Despite being categorized as poor quality by the WDI calculations, the groundwater samples are still deemed safe for drinking purposes. The HQ and health HI, considering both ingestion and dermal exposure, indicated values within acceptable ranges for both adults and children. The notable elevation in PAH concentrations within the study locations may be closely tied to the region's susceptibility to flooding due to intense rainfall events. This recurring phenomenon leads to water overflow and substantial flooding, resulting in significant leaching and erosion of topsoil particles that carry deposited or accumulated PAHs, along with potential crude oil spillage. It is noted that Flu, Nap, and DaA levels were consistently high across all the sampling points, showing the significance of these PAHs in the studied area's environmental dynamics. The irrigation quality assessment of the investigated water sources revealed varying levels of Na and magnesium hazards, alkalinity potential, and Na dominance. These findings show the importance of considering crop and soil requirements when selecting irrigation water sources and the potential need for mitigation strategies to address water quality challenges. These findings collectively emphasize the need for continuous monitoring and proactive environmental management strategies to mitigate the impact of pollutants in the region and safeguard public health. Together, these findings emphasize the need for continuous monitoring and proactive environmental management to reduce the impact of pollutants in these areas and for public health concerns.

The authors acknowledge and appreciate Springboard Laboratories for their research assistance. We also appreciate Mr Maxwell of Alvan Ikoku College of Education for his assistance during the field analysis.

I.S.C. conceptualized the whole article, developed the methodology, software, and validated the process, rendered support in formal analysis, wrote the original draft, wrote the review and edited the article. C.A.J. rendered support in formal analysis, wrote the review and edited the article. J.N. rendered support in formal analysis, wrote the review and edited the article. C.E.E. developed the methodology, software, and validated the article, rendered support in formal analysis, wrote the review and edited the article.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Akakuru
O. C.
,
Eze
C. U.
,
Okeke
O. C.
,
Opara
A. I.
,
Usman
A. O.
,
Iheme
O. K.
,
Ibeneme
S. I.
&
Iwuoha
P. O.
(
2022
)
Hydrogeochemical evolution, water quality indices, irrigation suitability, and pollution index of groundwater
(PIG) around Eastern Niger Delta, Nigeria
,
International Journal of Energy and Water Resources
6, 389–411.
doi:10.1007/s42108-021-00162-0
.
Akakuru
O. C.
,
Adakwa
C. B.
,
Ikoro
D. O.
,
Eyankware
M. O.
,
Opara
A. I.
,
Njoku
A. O.
,
Iheme
K. O.
&
Usman
A. O.
(
2023
)
Application of artificial neural network and multi-linear regression techniques in groundwater quality and health risk assessment around Egbema, Southeastern Nigeria
,
Environmental Earth Science
. 82 (77), 1–20.
doi:10.1007/s12665-023-10753-1
.
Amadi
A. N.
,
Udoh
B. O.
,
Ozoji
T. M.
,
Akande
W. G.
,
Shaibu
I.
,
Dan-Hassan
M. A.
,
Hussaini
S. U.
&
Oguntade
A. S.
(
2020
)
Quality assessment of groundwater in parts of Niger Delta, southern Nigeria using metal pollution index and factor analysis
,
The Pacific Journal of Science and Technology
,
21
(
1
),
280
289
.
Ambade
B.
,
Sethi
S. S.
,
Giri
B.
,
Biswas
J. K.
&
Bauddh
K.
(
2022
)
Characterization, behavior, and risk assessment of polycyclic aromatic hydrocarbons (PAHs) in the estuary sediments
,
Bulletin of Environmental Contamination and Toxicology
108, 243–252.
Banks
E. W.
,
Cook
P. G.
,
Owor
M.
,
Okullo
J.
,
Kebede
S.
,
Nedaw
D.
,
Mleta
P.
,
Fallas
H.
,
Gooddy
D.
,
MacAllister
D. J.
&
Mkandawire
T.
(
2021
)
Environmental tracers to evaluate groundwater residence times and water quality risk in shallow unconfined aquifers in sub-Saharan Africa
,
Journal of Hydrology
,
598
,
12575
.
Bhattacharya
P.
,
Adhikari
S.
,
Samal
A. C.
,
Das
R.
,
Dey
D.
,
Deb
A.
,
Ahmed
S.
,
Hussein
J.
,
De
A.
,
Das
A.
&
Joardar
M.
(
2020
)
Health risk assessment of co-occurrence of toxic fluoride and arsenic in groundwater of Dharmanagar region, North Tripura (India)
,
Groundwater for Sustainable Development
,
11
,
100430
.
Birhanu
B.
,
Kebede
S.
,
Charles
K.
,
Taye
M.
,
Atlaw
A.
&
Birhane
M.
(
2021
)
Impact of natural and anthropogenic stresses on surface and groundwater supply sources of the upper awash sub-basin, central Ethiopia
,
Frontiers in Earth Science
,
9
,
321
.
Bunting
S. Y.
,
Lapworth
D. J.
,
Crane
E. J.
,
Grima-Olmedo
J.
,
Koroša
A.
,
Kuczyńska
A.
,
Mali
N.
,
Rosenqvist
L.
,
Van Vliet
M. E.
,
Togola
A.
&
Lopez
B.
(
2021
)
Emerging organic compounds in European groundwater
,
Environmental Pollution
,
269
,
115945
.
Chika
C. P.
,
Okeke
O. C.
,
Okengwo
N.
&
Onunkwo
A. A.
(
2014
)
Hydrogeochemical assessment of groundwater quality in parts of Ohaji Egbema Eastern Niger Delta, Nigeria
,
Universal Journal of Environmental Research and Technology
,
4
(
6
),
307
316
.
Duru
C. E.
,
Okoro
I. P.
&
Enyoh
C. E.
(
2017
)
Quality assessment of borehole water within orji mechanic village using pollution and contamination models
,
International Journal of Chemistry, Material and Environmental Research (IJCMER)
,
4
(
3
),
123
130
.
Edokpayi
J. N.
,
Rogawski
E. T.
,
Kahler
D. M.
,
Hill
C. L.
,
Reynolds
C.
,
Nyathi
E.
,
Smith
J. A.
,
Odiyo
J. O.
,
Samie
A.
,
Bessong
P.
&
Dillingham
R.
(
2018
)
Challenges to sustainable safe drinking water: A case study of water quality and use across seasons in rural communities in Limpopo province, South Africa
,
Water
,
10
(
2
),
159
.
Ejiogu
B. C.
,
Opara
A. I.
,
Nwosu
E. I.
,
Nwofor
O. K.
,
Onyema
J. C.
&
Chinaka
J. C.
(
2019
)
Estimates of aquifer geo-hydraulic and vulnerability characteristics of Imo state and environs, Southeastern Nigeria, using electrical conductivity data
,
Environmental Monitoring and Assessment
,
191
(
4
),
238
.
doi:10.1007/s10661-019-7335-1
.
Enyoh
C. E.
&
Isiuku
B. O.
(
2020
)
Characterization of some soils from flood basin in Amakohia, Owerri, Nigeria
,
International Journal of Environmental Analytical Chemistry
102, 3766–3785.
doi:10.1080/03067319.2020.1773455
.
Enyoh
C. E.
,
Verla
A. W.
&
Egejuru
N. J.
(
2018
)
Ph variations and chemometric assessment of borehole water in Orji, Owerri Imo State, Nigeria
.
Journal of Environmental Analytical Chemistry
,
5
(
238
),
1
9
.
doi:10.4172/2380-2391.1000238
.
Francis
C. I.
,
Chidi
E. D.
,
Benaiah
O. I.
&
Jacienta
N. A.
(
2021
)
Ecological risk assessment of the levels of polycyclic aromatic hydrocarbons in soils of the abandoned sections of Orji Mechanic Village, Owerri, Imo State, Nigeria
,
Bulletin of the National Research Centre
,
42
,
2
16
.
Hu
G.
,
Chen
S.
,
Shi
W.
,
Zhang
B.
,
Zhang
Y.
,
Huang
J.
,
Chen
J.
,
Giesy
J.
&
Yu
H.
(
2014
)
Identification of polycyclic aromatic hydrocarbons in soils in Taizhou, East China
,
Environmental Geochemistry and Health
,
1
,
13
.
Ibe
F. C.
,
Opara
A. I.
,
Ibe
B. O.
,
Adindu
B. C.
&
Ichu
B. C.
(
2018
)
Environmental and health implications of trace metal concentrations in street dust around some electronic repair workshops in Owerri, Southeastern Nigeria
,
Environmental Monitoring and Assessment
,
190
(
696
),
1
14
.
Ibe
F. C.
,
Ihenetu
S. C.
,
Njoku
V. O.
,
Ahmad
A. S.
&
Bello
Y.
(
2023
)
Physicochemical, pollution, health risk, and irrigation quality evaluation of some surface water sources in Nwangele, Southeastern, Nigeria
,
International Journal of Energy and Water Resources
23, 1–20.
doi:10.1007/s42108-023-00273-w
.
Ihenetu
S. C.
,
Ihenetu
F. C.
&
Kalu
G. I.
(
2017
)
Determination of heavy metals and physicochemical parameters of crude oil polluted soil from ohaji egbema local government in Imo state
,
Open Journal of Yangtze Gas and Oil
,
2
,
161
167
.
Ihenetu
S. C.
,
Ochule
B. I.
,
Enyoh
E. C.
,
Ibe
F. C.
,
Verla
A. W.
&
Isiuku
B. O.
(
2020
)
Pollution and health risks assessment of ground water sources around a waste disposal site in owerri west local government area of Imo state
,
Journal of Material and Environmental Science
,
11
(
9
),
1560
1573
.
Ihenetu
S. C.
,
Victor
O. N.
,
Chinweuba
A.
,
Ibe
F. C.
&
Ekeoma
N. B.
(
2021b
)
A review: Effects of air, water and land dumpsite on human health and analytical methods for determination of pollutants
,
Analytical Methods in Environmental Chemistry Journal
,
4
(
3
),
80
106
.
Itodo
A. U.
,
Sha'Ato
R.
&
Arowojolu
M. I.
(
2018
)
Polycyclic aromatic hydrocarbons in water samples from a Nigerian bitumen seepage: Gas chromatography-mass spectrometry quantification
,
Bangladesh Journal of Scientific and Industrial Research
,
53
(
4
),
319
326
.
Joshua
N. E.
,
John
O. O.
,
Oluwasein
E. P.
&
Titus
A. M. M.
(
2016
)
Determination and distribution of polycyclic aromatic hydrocarbons in rivers, sediments, and wastewater effluents in Vhembe District, South Africa
,
International Journal of Environmental Research and Public Health
,
13
,
387
.
Liu
Y.
,
Yan
C. Q.
,
Ding
X.
,
Wang
X.
,
Fu
Q. Y.
,
Zhao
Q. B.
,
Zhang
Y. H.
,
Duan
Y. S.
,
Qiu
X. H.
&
Zheng
M.
(
2017
)
Sources and spatial distribution of particulate polycyclic aromatic hydrocarbons in Shanghai, China
,
Science of the Total Environment
,
584–585
,
307
317
.
doi:10.1016/j.scitotenv.2016.12.134
.
Mammeri
A.
,
Tiri
A.
,
Belkhiri
L.
,
Salhi
H.
,
Brella
D.
,
Lakouas
E.
,
Tahraoui
H.
,
Amrane
A.
&
Mouni
L.
(
2023
)
Assessment of surface water quality using water quality index and discriminant analysis method
,
Water
,
15
,
680
.
doi:10.3390/w15040680
.
Mohamed
E. G.
&
Hoda
S.
(
2023
)
Assessing surface water uses by water quality index: Application of qalyubia governorate, southeast Nile delta, Egypt
,
Applied Water Science.
,
13
,
181
.
doi:10.1007/s13201-023-01994-3
.
Olayinka
O. O.
,
Adewu
A. A.
,
Olanrewaju
O. O.
&
Aladesida
A. A.
(
2018
)
The concentration of PAHs and estimated human health risks of water samples around Atlas Cove, Lagos, Nigeria
,
Journal of Health and Pollution
,
8
(
20
),
181210
.
doi:10.5696/2156-9614-8.20.181210
.
Reberski
J. L.
,
Terzić
J.
,
Maurice
L. D.
&
Lapworth
D. J.
(
2022
)
Emerging organic contaminants in karst groundwater: A global level assessment
,
Journal of Hydrology
,
604
,
127242
.
doi:10.1016/j.jhydrol.2021.127242
.
Ward
J. S.
,
Lapworth
D. J.
,
Read
D. S.
,
Pedley
S.
,
Banda
S. T.
,
Monjerezi
M.
,
Gwengweya
G.
&
MacDonald
A. M.
(
2020
)
Large-scale survey of seasonal drinking water quality in Malawi using in situ tryptophan-like fluorescence and conventional water quality indicators
,
Science of the Total Environment
,
744
,
140674
.
World Water Quality Alliance (WWQA)
. (
2021
)
Assessing Groundwater Quality: A Global Perspective: Importance, Methods, and Potential Data Sources. A report by the Friends of Groundwater in the World Water Quality Alliance. Information Document Annex for Display at the 5th Session of the United Nations Environment Assembly, Nairobi 2021. Available at: https://www.unigrac.org/resource/assessing-groundwater-quality-global-perspective.
Xu
Z. Y.
,
Wang
C. H.
,
Li
H. X.
,
Xu
S. D.
,
Du
J.
,
Chen
Y. J.
,
Ma
C.
&
Tang
J. H.
(
2021
)
Concentration, distribution, source apportionment, and risk assessment of surrounding soil PAHs in industrial and rural areas: A comparative study
,
Ecological Indicators
,
125
,
107513
.
doi:10.1016/j.ecolind.2021.107513
.
Zheng
H.
,
Xing
X. L.
,
Hu
T. P.
,
Zhang
Y.
,
Zhang
Y.
,
Zhang
J. Q.
,
Zhu
G. H.
,
Li
Y.
&
Qi
S. H.
(
2018
)
Biomass burning contributed most to the human cancer risk exposed to the soil-bound PAHs from Chengdu economic region, western China
,
Ecotoxicology and Environmental Safety
,
159
,
63
70
.
doi:10.1016/j.ecoenv.2018.04.065
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data