Leaching Kinetics of As, Mo, and Se from Acidic Coal Fly Ash Samples

Annually, coal-fired electric power plants produce large volumes of potentially hazardous coal combustion products (CCPs) including fly ash. Since majority of the coal fly ash and other CCPs deposited in dry land fills or wet lagoons, they pose risk of contamination to local environment. In this study, we present results of leaching kinetics for As, Mo, and Se from three acidic fly ash samples collected from coal-fired power plants in the southeastern United States. This study shows that the leachate concentrations of As, Mo, and Se increase over time. Three kinetics equations, pseudo-second order, Elovich, and power-function, are able to adequately describe the experimental leaching kinetics data. Experimental leaching data and modeling results indicate that the rate limiting leaching of As, Mo, and Se is controlled by the diffusional process responsible for transferring these elements from interior to the surface of the particles as well as the dissolution of the fly ash particles. Therefore, it is important to adopt effective containment/treatment schemes to avoid potential and persistent dispersion of trace elements from ash disposal facilities to surrounding environment for a long time.


Introduction
Coal-fired electric power plants generate large volumes of coal combustion products (CCPs) such as fly ash, bottom ash, boiler slag, flue gas desulfurization (FGD) materials, and various gases.Fly ash is the most voluminous fraction, accounting almost 53.0%(in 2013) and projected to increase to 55.0% (2022) of the total of CCPs produced in the United Sates [1].While recycling of fly ash for beneficial use in the United States alone has increased from 30% in 2000 [2] to nearly 54% in 2015 [3], the majority of the ash is impounded in lagoons and landfills located throughout the country [1] [4].Fly ash contains large quantities of major element constituents as oxides, hydroxides and sulfates.It also contains environmentally significant quantities of hazardous leachable trace elements such as As, Mo, and Se [5] [6].During combustion, the organic matter in coal is destroyed and, as a result, the concentrations of trace elements in fly ash are enhanced relative to the source coal [7].
The release of trace elements to the environment is of concern because of their potential toxicity.Several studies show that coal fly ash with elevated concentration of trace elements can readily release these elements into the environment [5] [6] [7].The leaching behavior of major and trace elements in fly ash varies with the properties of the ash, such as pH, composition, and the leachant.Owing to the production and deposition of large volumes of fly ash in landfills each year, the risk associated with potential release of fly ash into local environment is real.For example, an incident that was occurred on December 22, 2008 with the rupture of a containment structure spilled over 3.7 million cubic meters of wet coal ash at the Kingston coal-fired power plant of Tennessee Valley Authority (TVA) in Tennessee [8].Additionally, seepage and leakage of toxic elements such as As, Mo, Se etc. from coal ash ponds have been recently reported, and, in some cases, contaminating nearby water sources and creating health hazard risk to numerous parts of the southern US [9].
Total leachable amount as well as overall leaching behavior of hazardous trace elements from fly ash samples is important for determining the environmental consequences of potential release of fly ash into local environment.A large number of studies on dissolution kinetics have been conducted in the past for several minerals [10]; however, there have been very few studies on leaching kinetics of trace elements from contaminated geomedia [11] [12].Even there is no uniformity among several regulatory leaching protocols for leachant-solid material interaction period [13].Similarly, most of the laboratory leaching schemes also use a leaching period in the range of 3 to 168 h [14].In this situation, leaching kinetics, particularly, leaching behavior over time and factors controlling the rate limiting leaching of trace elements are very important for understanding the leaching persistency of trace elements from coal fly ash.In this study, we present results evaluating the leaching kinetics of As, Mo, and Se from three acidic fly ash samples derived from the Eastern Bituminous coals in the United States.

Fly Ash Samples and Characterization
Fresh acidic fly ash samples (HA, HB and MA) were collected from three electric power plants located in the southeastern USA.These fly ash samples were the combustion products of the Eastern Bituminous coals.All three fly ash samples were collected dry and were homogenized in the laboratory before using them for kinetics experiments.
The pH of the acidic fly ash samples are reported in the range of slightly acidic to neutral.Details on physical and chemical properties of these fly ash were described elsewhere [6].In general, these fly samples are more acidic at the surface than their bulk composition.The specific surface area of these fly ash samples measured with single-point Brunauer-Emmet-Teller (BET) method was reported to be 3.18, 2.11, and 3.46 m 2 /g for HA, HB, and MA, respectively.Similarly, the carbon contents in these fly ash samples are reported to be 6.74%, 6.37%, and 14.7% for HA, HB, and MA, respectively.Concentrations of As, Mo, and Se along with some major elements determined with microwave assisted acid digestion [15] are presented in Table 1.The concentrations of several other trace elements and their leachabilities have been reported elsewhere [6] [16] [17] [18] [19].

Kinetics Experiments
Leaching kinetics of As, Mo, and Se from acidic fly ash samples were conducted using jar leaching and batch leaching experiments.For both schemes of experiments, Barnstead nanopure water (18.2MΩ) was used as leachant.
Jar leaching experiments for kinetic study were performed using 2 L high density polyethylene (HDPE) bottles at a 1:30 solid:liquid ratio.For these experiments, 60 g of each fly ash was mixed with 1.8 L of nanopure water and agitated on an orbital platform shaker at 200 rpm.About 20 mL of the leachate solution was withdrawn at each sampling event from each experimental bottle at 1, 4, 8, 12, 24, 36, 48, 72, 96, and 120 h.The leachate supernatants were separated by centrifugation at 8500 rpm for 10 minutes and filtration through 0.2 μm syringe filters.After the separation, leachate solutions were acidified to 2% with ultrapure OPTIMA nitric acid and stored in refrigerator until chemical analysis with a Perkin Elmer Optima 3000DV inductively coupled plasma optical emission spectrometer (ICP-OES).
In compliment to the Jar leaching experiment, a separate experimental scheme  S2.

Kinetic Modeling of Leaching
Leaching of trace elements from fly ash could be explained by assuming that this process is similar to desorption/dissolution of elements from solid surface/materials (Ash et al., 2013).Desorption kinetics of different sorbates from several sorbents including soils, metal oxides, and others have been previously described using zero order, first order, second order, and their derivative equations [20] [21] [22].An additional suite of kinetic models such as Elovich equation [23], power function equation [24] [25], and parabolic diffusion equation [20] were also used previously for modeling desorption kinetics.The selection of kinetic models could depend on several factors such as simplicity to determine fitting parameters and its ability to describe the experimental data.For example, Evans and Jurinak (1976) used multiple simultaneous first order reactions while Chien and Clayton (1980) used Elovich equation to adequately describe the release and sorption kinetics of phosphate in soils [26] [27].
In this study, leaching kinetics of As, Mo, and Se are analyzed using three kinetic equations: pseudo-second order, Elovich, and power-function equations.
The attempt of using other kinetic models failed because of their inability to adequately describe the experimental data.Kinetic equations with their linear forms and fitting parameters are presented in Table 2. Kinetic parameters for pseudo-second order equation were obtained by plotting t versus t/q t .The q e was the reciprocal of the slope while k and h were obtained from intercept.If Table 2. Kinetic equations, their linear forms, and model fitting parameters.

Kinetic equations linear form Kinetic parameters
Pseudo-second order a ( ) a, n a Ho and McKay (1999) [21]; b Low (1960) [23]; c Aharoni and Sparks (1991) [25]; q e is equilibrium concentration of leaching element, q t is concentration of leaching element at any time t.The units of various model parameters (e.g., k, h, α, β, a, and n) are given in Tables 3-5.
leaching kinetics of As were to obey the Elovich equation, the plot of q t vs. ln(t) should yield a linear relation with a slope of (1/β) and an intercept of (1/β) ln(αβ).Similarly, if log transformed kinetic data plotted as log(q t ) vs. log(t) result in a linear line, then the leaching kinetic was said to be in compliance with power function equation.Power function rate constant and its order were then obtained from slope (n) and intercept (logα) of the fit.The coefficient of determination (R 2 ) for each model was obtained by using experimental and model derived data.The goodness of model fit was also evaluated by calculating two additional parameters, normalized deviation (ND) and normalized standard deviation (NSD) using Equations ( 1) and ( 2), respectively [28].
where n is number of experimental measurements, q t(exp) is the experimental concentration of element at time t, and q t(model) is the model predicted concentration of element at time t.The smaller the values of ND and NSD, the better is the fit of experimental data for the kinetic model.The leachate pH for HA fly ash decreased slightly during the early stage of leaching and increased to about 7 after 30 hours.For all fly ash samples, leaching of As, Mo, and Se was relatively faster at the early stage; however, their leaching trend became slightly sluggish or plateaued in the later part of the experiment (Figure 1).

Results of Batch Leaching Kinetic Experiments
The leachate pH and concentrations of As, Mo, and Se mobilized during longterm leaching kinetic experiments are presented in Figure 2. All fly ash samples showed increasing leachate pH until 10 weeks and it remained near-neutral over time.Although the trend is visible in results after Jar leaching tests (Figure 1), the results of long-term leaching experiments clearly illustrate the early rapid release followed by a slow but persistent increase in leachate concentration of these elements until 18 weeks (Figure 2).At the latter part of the experiment (after 18 week, leaching trend became slightly sluggish, and or plateaued.
The early relatively rapid leaching of As, Mo, and Se from these acidic fly ash As Se Mo pH generally associates with fast dissolution of sub-micron sized particles in the sample during mineral-water interaction [29].However, slow but persistent increases in the concentrations of these elements with increasing contact time could be related to slow dissolution of or solid-medium diffusion from the particle interior.In coal fly ash, As, Mo, and Se are reported to be present as surface enrichment and remain within solid particles either as impurities in crystalline phases or as heterogeneous/homogeneous mixture in glassy phases [5] [30].
Glassy aluminosilicate particles in fly ash are less stable phases in aqueous environment than mullite and quartz [31] [32].Such amorphous particles in fly ash samples could gradually dissolved during leaching experiments releasing trace elements over time.In our previous study, we observed that fly ash leachates produced with comparable interaction period are undersaturated with respect to quartz, mullite, and amorphous silica [6].Therefore, the increasing leachate concentrations of As, Mo, and Se with time measured in this study is likely a result from the slow dissolution of fly ash particles.

Modeling Results
The pseudo-second order, Elovich, and power-function kinetic models were used to describe the leaching behavior of As, Mo, and Se from acidic fly ash samples.Linear plots for these kinetic models for both leaching schemes are shown in Figure 3 and Figure 4. Similarly, the kinetic parameters for these models are presented in Tables 3-5 along with coefficient of determination (R 2 ), ND, and NSD.Results suggest that all three models were able to describe the general behavior of As, Mo, and Se from acidic fly ash samples.The pseudo-second order kinetic model appears to fit the leaching data for both leaching schemes for all fly ash samples (Table 3).For the adsorption phenomenon, the good compliance of kinetic data with pseudo-second order equation is attributed to a sorption phenomenon in which the rate-limiting step is associated with chemisorption involving valence forces through sharing or exchange of electrons between sorbent and sorbate [21].Conversely, in the case of leaching, the excellent fitting of kinetic data could be attributed to a bond breaking mechanism such as dissolution of amorphous fly ash particles.Fly ash samples contain glassy materials such as amorphous aluminosilicates particles [33], and as mentioned earlier, these glassy particles in fly ash are least stable in aqueous environment [32].
It is important to note that the leaching of any solute (e.g.trace elements) from a solid phase (e.g., fly ash) involves displacement of former from the latter into the leachant.This process is assumed to be consisted of multiple steps in- Figure 4. Kinetic leaching modeling using pseudo-second order (a.HA, b.HB, and c.MA), Elovich (d.HA, e. HB, and f.MA), and power function (g.HA, h.HB, and i. MA) models for batch leaching of As, Mo, and Se from acidic fly ash samples.at the surface by dissolution of host material, it is assumed that they instantaneously mix with the leachant owing to the consistent agitating environment during experiment.Therefore, the rate limiting leaching kinetics during the experiment is related only to the either gradual dissolution of fly ash particles or diffusional transport from particle interior to the surface.However, in natural environments where a comparable agitating condition is absent, the transport related processes active in the exterior milieu of the particles are also important for rate limiting leaching of trace elements from fly ash.In general, diffusive transport within the solid particles and/or dissolution of particles at the surface, transfer of elements from solid-leachant interface to leachant, and dispersive as well as diffusive transport in the leachant are important rate limiting processes in natural environment.
Similarly, Elovich equation also shows a strong ability to describe the leaching behavior of these trace elements from fly ash samples during batch leaching (Figure 4).However, the fits of jar leaching data with Elovich equation for Mo and Se from HB and As and Se from MA are relatively weak (Table 4).Nevertheless, this kinetic model fits the jar leaching data for Mo and Se from HB and As and Se from MA far better than other kinetic equations (e.g., first order and second order kinetic equations).The Elovich equation, which is based on assumption that the adsorption energy increases with surface coverage [23], has been previously also used to describe the kinetics of desorption of oxyanions from soils and soil constituents [27].In their study, Chein and Clayton (1980) suggested that an increase in α with or without change in β would indicate a relatively rapid reaction [27].In this perspective, the α values are relatively larger for the same elements during long-term leaching experiments, and could indicate a relatively faster leaching than during jar leaching experiments [27].However, these two experiments were conducted at different solid: liquid ratio and different time scales, and could not be directly comparable.Furthermore, such inference drawn from nature of α and β variation could be questionable at different situations [22].
The power-function equation is also able to describe most of the leaching data for both leaching schemes with few exceptions (Figure 3 and Figure 4; Table 5).
Particularly, the model fits for Mo (from HA) and As and Se (from MA) are relatively poor for jar leaching data.Similarly, this model resulted in weak fits for As leached from HA and Se leached from HB (Table 5).
The good agreement of experimental data with the model equations indicate the presence of one or multiple rate limiting mechanisms controlling the release of trace elements from fly ash samples [21] [25].As indicated above, one of such potential mechanism could be dissolution of glassy fly ash particles where these trace elements are disseminated.Similarly, another likely mechanism for rate-limiting process is migration of these elements from particle interior to the surface.A similar diffusion controlled rate limiting kinetics has been previously reported for late-stage sorption/desorption of phosphate and arsenate to/from ferric hydroxide particles [34] [35] [36].However, because these fly ash samples are rich in glassy constituents and our previous observation of increasing SiO 2 (aq) concentration over time [6], we propose that long-term kinetics of As, Mo, and Se release from these fly ash samples occur largely via dissolution of fly ash particles.During both (jar and long-term) experiments, the leachate pH shows, in general, an increasing trend with time.The proton consumption during dissolution of aluminosilicates is a well-documented phenomena in water-rock interactions [29] [37].Therefore, the increasing pH with increasing duration of fly ash-leachant interaction is a manifestation of ongoing dissolution of aluminosilicate fly ash particles, and such dissolution is attributed to be the major factor for slow but persistent increase in leachate concentrations of As, Mo, and Se over time.

Summary and Conclusions
The  Submit or recommend next manuscript to SCIRP and we will provide best service for you: Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.A wide selection of journals (inclusive of 9 subjects, more than 200 journals) Providing 24-hour high-quality service User-friendly online submission system Fair and swift peer-review system Efficient typesetting and proofreading procedure Display of the result of downloads and visits, as well as the number of cited articles Maximum dissemination of your research work Submit your manuscript at: http://papersubmission.scirp.org/Or contact jwarp@scirp.org

Figure 1
Figure 1 shows solution pH trends and concentrations of As, Mo, and Se released during jar leaching tests of the fly ash samples.Leachate pH slightly increased until 20 hours and stabilized thereafter for HB and MA fly ash samples.

Figure 2 .
Figure 1.Leachate pH trends and concentrations of As, Mo, and Se elements mobilized during fly ash jar leaching tests: a) HA fly ash, b) HB fly ash, and c) MA fly ash.

Figure 3 .
Figure 3. Kinetic leaching modeling using pseudo-second order (a.HA, b.HB, and c.MA), Elovich (d.HA, e. HB, and f.MA), and power function (g.HA, h.HB, and i. MA) models for jar leaching of As, Mo, and Se from acidic fly ash samples.
cluding chemical interactions such as dissociation of chemical bonds and transport of slackened trace elements from solid phase to the leachant.Once the ions of the trace elements reach to the particle surface from interior or become loose leaching tests of three acidic fly ash samples collected from coal-fired power plants in the southeastern United States show the increasing leachate concentrations of As, Mo, and Se with time.The jar leaching experiments show that the leachate concentrations of these elements are relatively low during the early

Table 1 .
Environmentally available concentrations (mg/kg) of major elements, As, Mo, and Se in fly ash samples.
a ICP-OES limit of detection (LOD) values are given in mg/L.

Table 3 .
Leaching kinetic parameters for pseudo-second order model.

Table 5 .
Leaching kinetic parameters for power function model

Table S1 .
Concentrations of As, Mo, and Se during Jar leaching tests.

Table S2 .
Concentrations of As, Mo, and Se during long-term batch leaching tests.