2.1 Study area
The study area is Ashaka cement plant located in Ashaka town, 9 km north of Bajoga, Gombe State, Nigeria. The map of the area is shown in Fig. 1, located between longitudes 11o 28' 30”E and 11o 29’ 30”E and latitudes 10o55'30”N and 11o56’30”N, and The climate is tropical with two distinct seasons (dry and wet). In the dry season between November-February is harmattan, and the hot period between March-April is of the dry season with temperature ranging from 31.1 to 42.1 oC. The maximum rainfall ranges from 800 mmto 900 mm (August-September) and a minimum of 250 mm to 350 mm (May-June), with relative humidity ranging from 60–80% [22].
2.2 Sample collection and preparation
A systematic sampling technique was used for sample collection with little modification [23]. The surrounding soil of a cement plant was collected at depths ranging from 0–15 cm (topsoil). Four subsamples within each sampling point were mixed to form a composite sample. A total of ninety three soil samples were collected. Approximately 50 g of soil samples were transferred into polyethylene bags, and subsequently transported to the laboratory and air-dried at room temperature for two weeks [24]. Before the determination of PTEs, substantial impurities were detached from the air-dried soil samples, pulverized using a mortar and pestle, and sieved viaa U.S. No. 10 (2 mm) mesh. Approximately 0.5 g of the soil was transferred into a Teflon cup, and a binary combination of 3.75 mL of HCl and 1.25 mL of HNO3 was added as previously reported [25]. The resulting solutions were then analyzed for PTE concentration using inductively coupled plasma optical emission spectrometry (ICP-OES, Agilent 720-ES).
2.3 Quality assurance
The precision of the analytical procedure was determined by a recovery study. The recovery study was conducted by determining the PTE concentration in triplicate in spiked and un-spiked samples. The ICP-OES (Agilent 720-ES) was first calibrated using a multi-element standard solution (QCSTD-27). The percentage of mean recovery of the PTEs ranged from 98–104%, and the limits of detection for Al, As, Ba, Cd, Co, Cr, Cu, Fe, Hg Mn, Mo, Ni, Pb, Sb, Sc, Se, Sr, Ti, V and Zn were 0.013, 0.003, 0.004, 0.0001, 0.001, 0.0004, 0.0008, 0.01, 0.005, 0.0006, 0.0001, 0.005, 0.001, 0.0001, 0.0001, 0.003, 0.003, 0.003, 0.009 and 0.0003 µg/L, respectively.
2.4 Data analysis
The data obtained were statistically evaluated (simple descriptive and inferential) using SPSS software version 25. The generated data were used to estimate the potential ecological and health risks.
2.5 Geochemical load index (GLI)
The geochemical load index (GLI) was used to evaluate soil pollution by potentially toxic elements, as reported by other relative studies [12, 26, 27].
GLI = Log2 [\(\frac{{\text{C}}_{\text{i}}}{\text{G}\text{B}\text{V} \times 1.5 }\)] (1)
Where Ci₌ measured elements (mg.kg− 1), GBV₌ geóchemical backgróund vąlue of element, 1.5 ₌ control vąlues attributed to lithogenic variątion in the soil. The evaluation parameters of the geochemical analysis are presented in Table S1 (in the supplementary information).
2.6 Ecological risk assessment
Ecological risk was used to evaluate the overall soil pollution the expression proposed by Hakanson, [28].
Cf =\(\frac{{\text{C}}_{\text{i}}}{{\text{C}}_{\text{o}}}\) (2)
$${\text{E}}_{\text{r}}^{\text{i}} = {\text{T}}_{\text{r}}\frac{{\text{C}}_{\text{i}}}{{\text{C}}_{\text{o}}}$$
3
Ri = \(\sum {\text{E}}_{\text{r}}^{\text{i}}\) (4)
Where Ci ₌ measured elements, Co ₌background value of elements, Tr ₌ toxic response factor,\({\text{E}}_{\text{r}}^{\text{i}}\)₌ potential ecological risk factor, and Ri₌ summation of \({\text{E}}_{\text{r}}^{\text{i}}\). The Tr and Co values are shown in Table S2.
2.7 Health risk assessment
The PTEs measured in soil samples were used to evaluate health risk using the United States Environmental Protection Agency (USEPA) model and other previous studies [12, 24, 29]. Exposure through the three pathways was used.
ADDinh=\(\frac{\text{C} \times \text{I}\text{n}\text{h}\text{R} \times \text{E}\text{F} \times \text{E}\text{D} }{\text{P}\text{E}\text{P} \times \text{B}\text{W}\times \text{A}\text{T}}\) (5)
ADDing =\(\frac{\text{C} \times \text{I}\text{n}\text{g}\text{R} \times \text{E}\text{F} \times \text{E}\text{D} \times \text{C}\text{F}}{\text{B}\text{W} \times \text{A}\text{T}}\) (6)
ADDderm = \(\frac{\text{C} \times \text{S}\text{L} \times \text{S}\text{A} \times \text{A}\text{B}\text{S} \times \text{E}\text{F} \times \text{E}\text{D} \times \text{C}\text{F}}{\text{B}\text{W} \times \text{A}\text{T}}\) (7)
The non-carcinogenic effect of PTEs using the hazard quotient (HQ) in soil samples using the expression below USEPA, [30].
HQ = \(\frac{\text{A}\text{D}\text{D}}{\text{R}\text{f}\text{D}}\) (8)
The hazard index (HI) represents the summation of HQ multiple routes [31]. When HQ or HI > 1, the greater is the probability of non-carcinogenic adverse health effects to occur [30].
HI =\(\sum _{\text{n}}^{\text{i}}{\text{H}\text{Q}}_{\text{i}}\) (9)
Cancer risk (CR) was evaluated using the expression below [30].
CR = ADD × CSF (10)
Where ADD ₌ average daily exposure dose of PTEs, RfD ₌ reference dose and CSF₌ cancer slope factor. The evaluation parameters for the exposure assessment are listed in Table S3.