Study area
The present sample study was conducted in three districts: Buxar: 25°34'29"N-83°58'43"E, Bhojpur: 25°33'47"N-84°40'17"E and Patna: 25°36'8"N 85°7'10"E. All three districts are situated in the western part of the Bihar. This is the most fertile alluvial plain of the Ganges at an average elevation of 65 meters above sea level. These areas receive about 85% of the total rainfall (800-850 mm) from the southwest monsoon (ESSO/IMD, 2020). Both non-irrigated and irrigated lands are used for agriculture. Rice, wheat, gram and pulses are the main crops; vegetables are also grown abundantly in some parts.
Sample collection, preparation and analysis
The sampling was done in three different districts (Buxar, Bhojpur and Patna) of Bihar from February to March, 2020 (Figure 1). Composite soil was collected at a depth of about 5-15 cm from the surface of arable areas from fields situated close to the river bank. A composite sample (200-500 g) from each sampling site from 4-5 different points was taken. The aim was to ensure that the composite sample better represents the selected field soil. The coordinates of the sites were noted via GPS (Gramin Oregon 650). To eliminate the moisture from the collected soil samples, the soil was dried in the sun for 4-5 days. The soil was crushed and its gravel, organic debris and small stone pieces were removed by using a 2 mm sieve. Soil pH, EC, and salt concentration were determined in soil suspension prepared in the ratio of 1:2 (Soil-Water) with the help of Multi-Parameter (PCSTestrTM35 OAKlON UK). To find out the oxidation reduction potential (ORP) of soil, 1 g of soil was mixed with 10 ml of milli-Q water and measured with the help of Aakiro Aquasol digital (AM-ORP-01) meter. The percentage of organic carbon and available phosphorus in the soil samples were determined by the rapid titration method of Wakley and Black and Oleson method.
Wave Dispersive X-ray Fluorescence
WDXRF (Wave Dispersive X-ray Fluorescence) (S8 Tiger, Bruker, Mannheim, Germany) method was used for elemental detection. The press pellet method was used for the analysis of As (Upadhyay et al. 2021). In the press pellet method, 1 g of boric acid powder was mixed with 4 g of dried and sieved soil sample in a homogeneous manner. Round pellet form of soil samples was formed with the help of the KBr (potassium bromide) press pellet machine. Analysis of these round pellet forms of samples in XRF analyzer was done by X-ray beam with the help of existing automatic detector. To ensure qualitative accuracy, standard reference materials (SRMs) and blanks were handled alongside the samples and were analysed in triplicates as per the details given in Upadhyay et al. (2021).
Risk assessment indices
Contamination factor (CF)
The contamination factor has been widely used in many previous studies to determine the metal contamination level at the soil surface (Jiang et al. 2018; Li et al. 2018a; Li et al. 2018b). CF is the ratio of the concentration of each metal in the soil to the background value or baseline value of that metal in soil and it is calculated by the given formula:
CF = (CHM/CBV) (1)
where CHM and CBV indicate heavy metals concentration in soil and reference value/ background value of heavy metal in soil respectively. The reference value for As was taken as 6.83 mg kg-1 (Rinklebe et al. 2019). According to CF value, contamination level may be classified as low (0>CF <1), moderate (1 ≦ CF < 3), considerable (3 ≦ CF < 6) and very high (CF > 6).
Potential ecological risk assessment
The potential ecological risk (PER) index based on the heavy metal concentration, sensitivity, toxicity, type and reference value was proposed by Håkanson (1980). The risk index was used in this study to estimate the ecological hazards of HMs present in the soil. The formulas are as follows:
Cf = CF = (CHM/CBV) (2)
Er = Tr * Cf (3)
RI = ∑Er = ∑Tr * Cf (4)
Where the method of computing Cf is exactly the same as that of CF (see Eq. (1)), Er represents monomial potential ecological risk factor and HM toxic response factor is denoted by Tr (Hg=40, As=10, Cd=30, Cu=5, Pb=5, Ni=5, Cr=2, Zn=1) (Wang et al. 2018). As per PER, ecological risk can be categorized as low (Er ≤ 40), moderate (40 < Er ≤ 80), considerable (80 < Er ≤ 160), high (160 < Er ≤ 320) and very high (320 < Er).
Geo-accumulation index (Igeo)
The soil accumulation index was calculated using the formula developed by Müller (1969), which reports the extent of metal contamination in the soil (Kormoker et al. 2019; Islam et al. 2020).
Igeo = log2(Cm/1.5Bm) ………(5)
In the above equation, Cm is the concentration of the metal in the present soil sample and Bm is the background value of the same metal. The factor 1.5 has been used with reference to possible variations in the background value. For the classification of various degrees of contamination level, total 7 classes the geo-accumulation index was reported.
Values
|
Classes/Grades
|
Igeo< 0
|
Class 0 (nearly uncontaminated):
|
0 < Igeo < 1
|
Class 1 (uncontaminated to moderately contaminated)
|
1 < Igeo < 2
|
Class 2 (moderately contaminated)
|
2 < Igeo < 3
|
Class 3 (moderately to severely contaminated)
|
3 < Igeo < 4
|
Class 4 (severely contaminated)
|
4 < Igeo < 5
|
Class 5 (severely to extremely contaminated)
|
5 < Igeo
|
Class 6 (extremely contaminated)
|
Health risk assessment
Oral intake is very important for the entry of As into the human body among ingestion, oral intake, dermal contact and inhalation. In this study, the As non-carcinogenic risk assessment was done based on the average daily dose (ADD) and hazard quotient (HQ) through soil only (Masood et al. 2019).
ADD = (EPC × IR × ED × EF)/(BW×AT) ………. (6)
EPC used in the formula denotes exposure point concentration (mg kg-1) in soil, ED exposure duration (assumed 70 years), EF exposure frequency (365 day year-1; Tyagi et al. 2020), IR ingestion rate (100 mg day-1; USEPA 2011) for adult, BW body weight (70 Kg; (USEPA 2005) and AT average period of exposure (25,550 days; Masood et al. 2019) respectively. An exposed population is considered at risk when HQ is greater than one (Masood et al. 2019; Shahid et al. 2017).
The ratio of ADD and RfD represents the hazard quotient (Masood et al. 2019; Shahid et al. 2017).
HQ=ADD/RfD ……….(7)
where RfD is representing the reference dose (mg kg-1 d-1). The maximum permissible exposure to the human population is estimated through daily exposures taking into account the vulnerable groups during the lifetime. For As, Rf D is considered as 0.0003 mg kg-1 d-1 (Tyagi et al. 2020).
The carcinogenic health risk was assessed as Incremental Lifetime Cancer Risk (ILCR) by following (Tyagi et al. 2020). Using the mean As concentration in soil, the ILCR was calculated by using the following formula.
ILCR = ADD × SF ……….(8)
SF is considered as slope factor. The SF is considered as 1.5 for As (USEPA 2017). As per USEPA (2001), the ILCR between 10−6 to 10−4 is considered to be within the safe limit.
Statistical analysis
Statistical Package for the Social Sciences for Windows (SPSS) version 16.0 and Microsoft Office Excel Version 16.0 (Microsoft Corporation, Washington, USA) was used in the study for correlation analysis. To determine significance in all tests significance level of 0.01 was used. Maximum, minimum, standard error, mean and median values were calculated. One-way ANOVA and Tukey tests were performed to measure how significantly the averages As content of samples collected from different regions differed from each other.