1. Introduction
Fluorine (F) has received increasing attentions because of its high toxicity and reactivity [1, 2]. The Ministry of Environment has established the regulatory level of F at 400 mg/kg for residential and agricultural areas in the Korean Soil Environment Conservation Act [3]. Approximately 8 tons of highly toxic anhydrous hydrofluoric acid (AHF) was accidentally released from a chemical plant in Gumi, Korea into the environment on September 2012 [4]. AHF was rapidly converted to hydrogen fluoride (HF) and it dispersed throughout the atmosphere, and spread to the surrounding environment, including soil and groundwater. The F released into the environment could induce human health problems unless properly managed. Damage to human health was caused by HF gas exposure immediately after the accident; however, chronic health impacts related to the F presenting in the soil have been exerted continuously [5, 6].
Many studies have clearly demonstrated that the role of bioavailable concentrations (i.e., water-soluble and leachable concentrations) may be more important in determining the toxic effects of F in soil than that of total F concentration [7, 8]. Plant toxicity appeared to be considerably affected by water-soluble F concentrations in soil [9]. In addition, Geretharan et al. [10] observed that water-soluble F concentration is a key factor determining the retention and mobility of F in soil, and F concentration correlates with toxicity to soil microbes. Cronin et al. [5] suggested that the critical levels of the more bioavailable concentrations of F in soil, particularly in solution or in the water-soluble fraction, are probably the most important for toxicity and groundwater pollution. Bioavailable F concentrations should be considered for the proper risk management of F-contaminated soil.
Stabilization technology is a possible option for reducing the toxicity and mobility of contaminants in environmental media, including soil [11, 12]. Arsenic leaching from solid waste was successfully stabilized by binder materials including cement and fly ash [13], and steel slags [14]. The stabilization of Cd was tested using common amendment such as hydroxyapatite and biochar, and it was effectively immobilized in agricultural soils [15]. The reduction of extractable Pb, Cd, Cu, and Zn from smelter contaminated soils by combined with lime was also investigated by Lahori et al. [16]. The term stabilization refers to a remediation process involving chemical reactions for reducing the leachability of the contaminants [17] by mixing a contaminated medium with a stabilizing agent [18]. An appropriate stabilizing agent can convert bioavailable F into chemically stable forms, thus preventing the spread of F contamination and reducing the risk to human health. Among various stabilizing agents, calcium-based stabilizing agents, such as calcium hydroxide (Ca(OH)2), are most widely used to stabilize heavy metals as well as F [19–22]. The mechanism of F stabilization is explained by the reaction of Ca(OH)2 with F ions from NaF and HF resulting in insoluble calcium fluoride (CaF2) formation, as shown in Eq. (1) [23].
CaF2 is one of the most stable forms F-containing compounds, and sparingly soluble (0.016 g/L for CaF2), compared with MgF2 (0.13 g/L), Na3AlF6 (0.42 g/L), HF, SiF4, and NaF (40 g/L) [24]. Most studies on F stabilization focused on water, such as F stabilization in groundwater, wastewater, and drinking water [25–27]. To date, the efficiency of F stabilization efficiency in soil and related changes in human health risk have not been fully demonstrated.
The study was aimed at evaluating the stabilization of F-contaminated soil in terms of human health risk, and providing insight for the risk-based management of F-contaminated soil. We investigated F stabilization in soil using Ca(OH)2, and to determine the consequent reduction of human health risk. Lab-scale batch experiments were performed to determine the stabilization efficiencies with 3% and 6% Ca(OH)2. Modified sequential extraction was conducted to determine the stabilized F concentration in soil, and X-ray diffraction (XRD) analysis was performed to identify the F-containing compounds, such as CaF2, formed owing to stabilization. A human health risk assessment was performed, and changes in risks of F in soil from various exposure pathways were estimated before and after F stabilization.
2. Experimental
2.1. Preparation and Characterization of Soil Sample
Soil samples were collected from the area within a radius of 1 km from the site of accidental HF spillage in Gumi, Korea. Soil was sampled at a depth of 0–30 cm, air-dried, and passed through a 2 mm sieve to examine its physicochemical properties. Soil pH was measured at a soil:water ratio of 1:5, and organic matter content was determined using the Walkley-Black method [28]. Three textural fractions, namely clay, silt, and sand, were used to determine soil texture using the soil texture triangle recommended by the US Department of Agriculture (USDA). The total contents of 10 major elements (i.e., Al, Ca, Fe, K, Mg, Mn, Na, P, Si, and Ti) in the soil samples were analyzed using a wavelength dispersive X-ray fluorescence spectrometer (WD-XRF, PW2404, Phillips, Netherlands) at the Korea Basic Science Institute after preparing a glass bead using Li2B4O7 (soil:Li2B4O7 ratio of 1:10 (w/w)).
The F concentration in the soil samples was determined using the standard method of the Ministry of Korea, i.e., the UV/Visible spectrometric method [29]. Briefly, 1.0 g of soil sample was passed through a 75-μm sieve and 1.0 g of soil and 5.0 g of CaO were placed in a nickel pot, in a furnace at 500°C for 5 h and 800°C for 2 h. Then, 50 mL of 70% HClO4 and 25 mL of deionized water were added to the residue obtained in the previous step, and F ions were extracted by distillation at 135 ± 2°C. Subsequently, 50 mL of distilled sample was mixed with 10 mL sodium-2-(parasulfophenylazo)-dihydroxy-3, 6-napthalene disulfonate with zirconyl acid (Zr-SPADNS solution) in a 100 mL volumetric flask, and aliquots were analyzed using a spectrophotometer at 570 nm. For determining F ion concentration in the solution, the standard method of the Ministry of Korea, i.e., ion selective electrode (ISE) method, was used [29]. After the preparation step described above, the filtrate was mixed with total ionic strength adjustment (TISAB) solution at a 1:1 ratio (v/v), and the mixture was analyzed using ISE (pH-270L, Istek, Korea).
A batch equilibrium experiment was conducted to determine the adsorption capacity of F onto the soil sample used in this study. In this, 1 g of soil sample and 25 mL of NaF solution (i.e., final F concentrations of 5, 10, 100, 250, 500, 1,000, 2,500, and 5,000 mg/L) were shaken for 24 h at room temperature. After centrifugation, the supernatant was collected, and F concentrations in equilibrium solution (mg/L) were determined using the ISE method. The weight of F adsorbed per unit weight of soil at equilibrium (mg/kg) was calculated. Saturation of F concentration adsorbed onto the soil sample in this study occurred when F solution (1,000 mg/L) was introduced thus, the soil sample was artificially contaminated with 1,000 mg/L-F solution to simulate accidental spillage, and maximize stabilization efficiency. Therefore, 40 g of soil and 40 mL of F solution were mixed, followed by aging for 7 d.
2.2. F Stabilization Experiment Using Ca(OH)2
A lab-scale batch experiment was conducted for F stabilization with 3% and 6% Ca(OH)2 (w/w), relative to the soil samples. Distilled water at field capacity (13.2%) was added to the mixture of soil samples and Ca(OH)2, and agitated for 24 h at room temperature. The mixture was dried at room temperature, and four-step sequential extraction was performed to determine the soluble F concentration in the soil samples stabilized using Ca(OH)2. The four-step sequential extraction proposed by Xu et al. [30], was adopted to determine the stabilized F concentration in the soil samples. In the first step of the extraction of the water-soluble F fraction (fraction 1), 1 g of soil sample was stirred with 30 mL of distilled water for 15 h at 25 ± 0.5°C. Then, the residue was agitated with 30 mL of 20% (v/v) acetic acid/water for 15 h at 25 ± 0.5°C to extract the fraction of F bound to carbonates (fraction 2). Next, to determine the fraction of F bound to sulfides (fraction 3), the residue obtained from the second step was treated with 30 mL of 15% (v/v) nitric acid/water in a covered vessel for 1 h at room temperature and for 1 h at 85°C. Finally, the residual F fraction (fraction 4) was calculated by subtracting the F concentration obtained in steps 1 to 3 from the total F concentration determined using UV/Visible spectrometry. The sample obtained after each extraction step was centrifuged at approximately 5,000 g, and filtered through a 0.45 μm GHP filter. The F concentrations in the extracted solutions were determined using the ISE method. All experiments were performed in triplicate.
2.3. Human Health Risk Assessment
2.3.1. Exposure assessment and toxicity information
A residential scenario was considered for the human health risk assessment, assuming that the contaminated site would be developed for residential and other anthropogenic activities after F stabilization. Direct exposure pathways, including soil ingestion, inhalation of fugitive dusts from soil, and water ingestion were considered. Additional exposure pathways (e.g., dermal contact) were ignored, because of the absence of a chemical-specific dermal absorption fraction (ABSd) for F. Similarly, the inhalation of volatile F in soils was neglected, because inorganic contaminants, excluding mercury, are regarded non-volatile compounds [31]. The exposure parameters for child and adult used in this study were obtained from the Korean soil contaminant risk assessment guidance [32] and Korean exposure factors handbook [33] (Table 1). Table 2 describes certain parameters used for calculating the particulate emission factor (PEF) for fugitive soil dust inhalation using Eq. (2).
Default values provided by the U.S. EPA were used for the site-specific parameters (e.g., width of contaminated area (LS), area of contamination (A), and fraction of vegetative cover (G)) and diffusion height (D)). The equivalent threshold value of wind speed at 700 m (Ut) and F(x) were determined following the method of Cowherd et al. [34] Mean annual wind speed (Um) was obtained from the Korea Meteorological Administration, and wind speed in the mixing zone (V) was calculated from half of Um.
Toxicity values for the human health risk assessment, such as oral reference dose (RfDo) and reference concentration (RfC) of F were obtained from the Integrated Risk Information System [35] and are shown in Table 1. Because F is classified as a non-carcinogenic element, carcinogenic toxicity values, including that of oral slope factor (SF), are not provided.
2.3.2. Risk characterization
The human health risk assessment for F were conducted following the U.S. EPA [36, 37]. Eq. (3), (4), and (5) presented below were used to estimate average daily F doses from soil ingestion, inhalation of fugitive dust from soil, and water ingestion, respectively.
where, ADDs is average daily dose from soil ingestion (mg kg−1d−1), ADDw is average daily dose for water ingestion (mg L−1 d−1), ADE is average daily exposure (mg/m3), Cs is the exposure concentration of F in soil (mg/kg), Cw is the exposure concentration of F in water (mg/L), IRS is soil ingestion rate (mg/d), IRW is water ingestion rate (L/d), EF is exposure frequency (d/y), ED is exposure duration (y), AT is averaging time (d), BW is body weight (kg), ETrs is the soil exposure time for soil inhalation (h/d), and PEF is the particulate emission factor.
F concentration in groundwater was not measured in this study. Therefore, Cw, used to calculate the ADDw from water ingestion, was estimated using the soil-water phase partitioning equation for inorganic contaminants as follows in Eq. (6) [31, 38].
where, Cw is the exposure concentration of F in groundwater (mg/L), Cs is the exposure concentration of F in soil (mg/kg), Kd is the soil-water partition coefficient (L/kg), θw is water-filled soil porosity, θa is air-filled soil porosity, ρb is the dry soil bulk density (kg/L), and H’ is Henry’s constant. θw was estimated using Eq. (7), and θa was calculated as the difference between total porosity (n) and θw. n was expressed in terms of density according to Eq. (8).
where, ρb is the dry soil bulk density (kg/L), Ks is the mean annual wind speed (m/y), I is infiltration rate (m/y), and 1/(2b + 3) is a constant (0.07), and ρs is soil particle density (kg/L). Chemical-specific Kd values were obtained from the generic table on regional screening level provided by U.S. EPA. Site-specific values such as ρb, Ks, and the constant [i.e., 1/(2b + 3)] were selected according to soil texture (loam) [39]; ρs was used as the Korean soil-specific value provided by the Rural Development Administration (RDA). Infiltration rate was 14.44% of the average annual precipitation, based on data obtained from the Korea Meteorological Administration, considering the groundwater recharge rate in Korea. The parameters used for calculating F concentration in groundwater are summarized in Table 3. A dilution attenuation factor was not considered for conservative risk calculation in this study.
The estimated exposure to F (i.e., ADD and ADE) was divided by the toxicity values (i.e., RfDo and RfC). Total non-carcinogenic risks of F were determined based on the sum of the risk from soil ingestion, inhalation of fugitive dust from soil, and water ingestion and expressed as the hazard quotient (HQ; Eq. (9)).
where, HQ is the hazard quotient, RfDo is the oral reference dose (mg kg−1 d−1), RfC is the reference concentration (mg/m3). If HQ exceeded one, it was generally considered that there might be a concern of non-carcinogenic effects.
3. Results and Discussion
3.1. Soil Characterization
The physicochemical characteristics of the soil samples used in this study are summarized in Table 4. pH was 6.64, and organic matter content was 13.1%. The texture was loamy, with 44.9%, 33.3%, and 21.8% sand, silt, and clay, respectively. Elemental composition analysis using XRF showed that SiO2 content was the highest in the soil samples (76.9%) among all detected compounds, followed by Al2O3 (10.7%), Fe2O3 (2.68%), and K2O (2.63%) contents. Total F concentration in the soil sample was 1,654 ± 7.36 mg/kg. The modified sequential extraction results showed that the majority of F in the soil sample existed in a water-soluble form (i.e., a F concentration of 711 ± 23.4 mg/kg on average, and 43.0% in fraction 1) and as residue (i.e., 942 ± 23.6 mg/kg on average and 57.0% in fraction 4). On average, only 0.060% and 0.00002% F was bound to carbonates (1.00 ± 0.171 mg/kg in fraction 2) and sulfides (0.00004 ± 0.00004 mg/kg in fraction 3), respectively. Many researchers have investigated the relationship between the bioavailability (i.e., leachability or mobility) of heavy metals, as well as for F, and their chemical forms in soils. It is widely recognized that heavy metals, as well as F, are most toxic and mobile in their water-soluble form (fraction 1), and it is unlikely that in the residue form (fraction 4), it is accessible to living organisms, or exhibits any toxic effects [8, 40–42]. The leachable F concentration for stabilization in the soil samples used in this study was 711 ± 23.4 mg/kg (i.e., water-soluble F in fraction 1); furthermore, this concentration exceeded the regulatory soil F level set in the Korean Soil Environment Conservation Act (i.e., 400 mg/kg) [3].
3.2. F Stabilization in Soil Using Ca(OH)2
Table 5 showed that when 3% and 6% Ca(OH)2 were added to the soil samples, the chemical forms of F changed, compared with the untreated soil sample (i.e., control soil sample). Importantly, the bioavailable F concentration related to fraction 1 decreased to 3.7% and 3.5% on average (61.2 ± 1.65 and 57.9 ± 1.27 mg/kg in the 3% and 6% Ca(OH)2-treated soil samples, respectively) from 43.0% (711 ± 23.4 mg/kg in the control soil sample). Nearly all of the total F was present in the residue form (i.e., fraction 4) in the Ca(OH)2-treated soil samples. The ratio of F as residue increased markedly to 95.8% and 95.9% on average (1,584 ± 2.33 and 1,586 ± 1.31 mg/kg in the 3% and 6% Ca(OH)2-treated soil samples, respectively), compared with 57.0% on average (942 ± 23.6 mg/kg in the control soil sample). This indicated that most of the leachable F was converted into stable forms, such as CaF2 in the Ca(OH)2-treated soil samples; thus, Ca(OH)2 successfully acted as the stabilizing agent for F in the soil sample. The leachable concentrations F in the Ca(OH)2-treated soil samples were observed to be within the Korean soil regulatory level (i.e., 400 mg/kg) [3]. Meanwhile, there was no significant difference in stabilization between the 3% and 6% Ca(OH)2-treated soil samples.
Ca and F are known to form precipitates, such as CaF2, which could contribute to F stabilization [23, 43]. CaF2 formation in the Ca(OH)2-treated soil samples was detected by XRD. XRD patterns were compared between an untreated soil sample and a 6% Ca(OH)2-treated soil sample contaminated with 1.64% of F, and CaF2 mineral served as the control. Because the low F concentration in the soil used in this study (i.e., 1,654 ± 7.36 mg/kg of F contaminated soil sample) resulted in low sensitivity, it could not be detected by XRD; the soil sample contaminated with a high F concentration (1.64%), considering the solubility of NaF, was used for the noticeable detection of CaF2 formation. Fig. 1 showed that a significantly higher extent of CaF2 formation was observed in the Ca(OH)2-treated soil samples than in the untreated soil samples. Because the soil samples used in this study had low Ca contents, most of the CaF2 produced resulted from Ca(OH)2 addition. Aldaco et al. [44] observed that while adding doses of calcium ions to fluoride present in wastewater, the ions (Ca2+ and F−) form an insoluble product (CaF2). Additionally, CaF2 formation confirmed the possibility of CaF2 precipitate formation in the reaction between calcium ions and soluble F present in soil using PHREEQC, a geochemical model [45].
3.3. Effect of F Stabilization on Human Health Risk Reduction
A human health risk assessment was performed to evaluate the effect of F stabilization in soil on risk reduction. Because similar stabilization efficiencies were observed for the 3% and 6% Ca(OH)2-treated soil samples, only 6% Ca(OH)2-treated soil samples, having a relatively higher stabilization efficiency compared with that of 3% Ca(OH)2, were used for risk calculation in this study. Water-soluble F concentration, in fraction 1 determined by sequential extraction in this study, was considered the exposure concentration (e.g., Cs) for calculating ADDs and ADE. The exposure concentration used for risk calculation is presented in Table 6. The estimated HQ of F to non-carcinogenic human health risk under the stabilization conditions used is shown in Table 7. Non-carcinogenic risk was expected in the control soil samples (i.e., untreated soil sample), however the total risk reduced to less than 1.0 when the soil samples were treated with Ca(OH)2 for both children and adults, suggesting that the potential risk reduced. Total HQ values were the highest in the control soil samples for children (4.2), followed by the control soil samples for adults (3.0), Ca(OH)2-treated soil samples for children (0.34) and Ca(OH)2-treated soil samples for adults (0.25). Majority of the non-carcinogenic risk to children and adults in the control soil samples resulted from water ingestion, corresponding to 3.9 and 3.0, respectively. When the soil sample were treated with Ca(OH)2, the risk from water ingestion decreased sharply to 0.32 and 0.24 for children and adults, respectively. Interestingly, non-carcinogenic risk owing to soil ingestion was not expected in the control soil sample, although the exposure concentrations of F in soil (711 mg/kg) exceeded the Korean soil F regulatory level et in the Korean Soil Environment Conservation Act (i.e., 400 mg/kg) [3]. The HQ values resulting from soil ingestion were 0.23 and 0.018 in the control soil samples for children and adults, respectively. The HQ values of the Ca(OH)2-treated soil samples were lower than those of the control soil samples (0.019 and 0.0015 for children and adults, respectively). This indicated that the Korean soil F regulatory level might be excessively conservative. The regional screening levels (RSLs) superfund sites is set by the U.S. EPA as 3,100 mg/kg in residential soil for children based on non-carcinogenic risk, i.e., HQ < 1 [35]. The contribution of inhalation of fugitive dust from soil to non-carcinogenic risk was negligible for children and adults, i.e., 1.9 × 10−7 and 1.6 × 10−8 in the Ca(OH)2-treated and control soil samples, respectively. Most importantly, our results revealed that the F leaching from soil into groundwater is a critical pathway contributing to the high risk of contaminated water ingestion; therefore, the investigation of groundwater at a nearby site is essential for risk management for F-contaminated soil. Additionally, the reduction of non-carcinogenic risk from all exposure pathways upon soil treatment with Ca(OH)2 demonstrated that stabilization is an effective technique for the remediation of F-contaminated soil.
4. Conclusions
In this study, F stabilization was conducted using Ca(OH)2, and the reduction in human health risk was estimated by conducting a human health risk assessment. The bioavailable F concentration decreased remarkably, and this resulted from conversion of water-soluble F into stable CaF2. The non-carcinogenic risks (HQ) of F in soil to children and adults decreased to less than 1.0 after stabilization. The ingestion of F-contaminated water, caused by F leaching from soil into groundwater, contributes more to the risk than do other exposure types, such as soil ingestion and inhalation of fugitive dust from soil. The risk from soil ingestion made a smaller contribution to the total non-carcinogenic risk, although there is a stringent soil F regulatory level in Korea. These findings indicated that stabilization could be an effective option for reducing risk from F-contaminated soil, and a more realistic and reasonable management strategy in terms of soil regulatory level and investigation of contaminated sites should be developed based on the human health risk posed by the contaminant.