Effects of water content on reactive transport of 85Sr in Chernobyl sand columns

It is known that under unsaturated conditions, the transport of solutes can deviate from ideal advective-dispersive behaviour even for macroscopically homogeneous porous materials. Causes may include physical non-equilibrium, sorption kinetics, non-linear sorption, and the irregular distribution of sorption sites. We have performed laboratory experiments designed to identify the processes responsible for the non-ideality of radioactive Sr transport observed under unsaturated flow conditions in an Aeolian sandy deposit from the Chernobyl exclusion zone. Miscible displacement experiments were carried out at various water contents and corresponding flow rates in a laboratory model system. Results of our experiments have shown that breakthrough curves of a conservative tracer exhibit a higher degree of asymmetry when the water content decreases than at saturated water content and same Darcy velocity. It is possible that velocity variations caused by heterogeneities at the macroscopic scale are responsible for this situation. Another explanation is that molecular diffusion drives the solute mass transfer between mobile and immobile water regions, but the surface of contact between these water regions is small. At very low concentrations, representative of a radioactive Sr contamination of the pore water, sorption and physical disequilibrium dominate the radioactive Sr transport under unsaturated flow conditions. A sorption reaction is described by a cation exchange mechanism calibrated under fully saturated conditions. The sorption capacity, as well as the exchange coefficients are not affected by desaturation. The number of accessible exchange sites was calculated on the basis that the solid remained in contact with water and that the fraction of solid phase in contact with mobile water is numerically equal to the proportion of mobile water to total water content. That means that for this type of sandy soil, the nature of mineral phases is the same in advective and non-advective domains. So sorption reaction parameters can be estimated from more easily conducted saturated experiments, but hydrodynamic behaviour must be characterized by conservative tracer experiments under unsaturated flow conditions.     

Development of a scalable model for predicting arsenic transport coupled with oxidation and adsorption reactions

Understanding the fundamentals of arsenic adsorption and oxidation reactions is critical for predicting its transport dynamics in groundwater systems. We completed batch experiments to study the interactions of arsenic with a common MnO2(s) mineral, pyrolusite. The reaction kinetics and adsorption isotherm developed from the batch experiments were integrated into a scalable reactive transport model to facilitate column-scale transport predictions. We then completed a set of column experiments to test the predictive capability of the reactive transport model. Our batch results indicated that the commonly used pseudo-first order kinetics for As(III) oxidation reaction neglects the scaling effects with respect to the MnO2(s) concentration. A second order kinetic equation that explicitly includes MnO2(s) concentration dependence is a more appropriate kinetic model to describe arsenic oxidation by MnO2(s) minerals. The arsenic adsorption reaction follows the Langmuir isotherm with the adsorption capacity of 0.053micromol of As(V)/g of MnO2(s) at the tested conditions. The knowledge gained from the batch experiments was used to develop a conceptual model for describing arsenic reactive transport at a column scale. The proposed conceptual model was integrated within a reactive transport code that accurately predicted the breakthrough profiles observed in multiple column experiments. The kinetic and adsorption process details obtained from the batch experiments were valuable data for scaling to predict the column-scale reactive transport of arsenic in MnO2(s)-containing sand columns.     

Biodegradation of petroleum hydrocarbon vapors: laboratory studies on rates and kinetics in unsaturated alluvial sand

Predictions of natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone rely critically on information about microbial biodegradation kinetics. This study aims at determining kinetic rate laws for the aerobic biodegradation of a mixture of 12 volatile petroleum hydrocarbons and methyl tert-butyl ether (MTBE) in unsaturated alluvial sand. Laboratory column and batch experiments were performed at room temperature under aerobic conditions, and a reactive transport model for VOC vapors in soil gas coupled to Monod-type degradation kinetics was used for data interpretation. In the column experiment, an acclimatization of 23 days took place before steady-state diffusive vapor transport through the horizontal column was achieved. Monod kinetic parameters Ks and vmax could be derived from the concentration profiles of toluene, m-xylene, n-octane, and n-hexane, because substrate saturation was approached with these compounds under the experimental conditions. The removal of cyclic alkanes, isooctane, and 1,2,4-trimethylbenzene followed first-order kinetics over the whole concentration range applied. MTBE, n-pentane, and chlorofluorocarbons (CFCs) were not visibly degraded. Batch experiments suggested first-order disappearance rate laws for all VOCs except n-octane, which decreased following zero-order kinetics in live batch experiments. For many compounds including MTBE, disappearance rates in abiotic batch experiments were as high as in live batches indicating sorption. It was concluded that the column approach is preferable for determining biodegradation rate parameters to be used in risk assessment models.     

Transport of Sr2+ and SrEDTA2- in partially-saturated and heterogeneous sediments

Strontium-90 has migrated deep into the unsaturated subsurface beneath leaking storage tanks in the Waste Management Areas (WMA) at the U.S. Department of Energy's (DOE) Hanford Reservation. Faster than expected transport of contaminants in the vadose zone is typically attributed to either physical hydrologic processes such as development of preferential flow pathways, or to geochemical processes such as the formation of stable, anionic complexes with organic chelates, e.g., ethylenediaminetetraacetic acid (EDTA). The goal of this paper is to determine whether hydrological processes in the Hanford sediments can influence the geochemistry of the system and hence control transport of Sr(2+) and SrEDTA(2-). The study used batch isotherms, saturated packed column experiments, and an unsaturated transport experiment in an undisturbed core. Isotherms and repacked column experiments suggested that the SrEDTA(2-) complex was unstable in the presence of Hanford sediments, resulting in dissociation and transport of Sr(2+) as a divalent cation. A decrease in sorption with increasing solid:solution ratio for Sr(2+) and SrEDTA(2-) suggested mineral dissolution resulted in competition for sorption sites and the formation of stable aqueous complexes. This was confirmed by detection of MgEDTA(2-), MnEDTA(2-), PbEDTA(2-), and unidentified Sr and Ca complexes. Displacement of Sr(2+) through a partially-saturated undisturbed core resulted in less retardation and more irreversible sorption than was observed in the saturated repacked columns, and model results suggested a significant reservoir (49%) of immobile water was present during transport through the heterogeneous layered sediments. The undisturbed core was subsequently disassembled along distinct bedding planes and subjected to sequential extractions. Strontium was unequally distributed between carbonates (49%), ion exchange sites (37%), and the oxide (14%) fraction. An inverse relationship between mass wetness and Sr suggested that sandy sediments of low water content constituted the immobile flow regime. Our results suggested that the sequestration of Sr(2+) in partially-saturated, heterogeneous sediments was most likely due to the formation of immobile water in drier regions having low hydraulic conductivities.     

Sorption of isoxaflutole diketonitrile degradate (DKN) and dicamba in unsaturated soil

When analyzing the sorption characteristics of weakly sorbing or labile pesticides, batch methods tend to yield a high margin of error attributable to errors in concentration measurement and to degradation, respectively. This study employs a recently developed unsaturated transient flow method to determine the sorption of isoxaflutole's herbicidally active diketonitrile degradate (DKN) and dicamba. A 20-cm acrylic column was packed with soils with varied texture that had been uniformly treated with 14C-labeled chemical.The antecedent solution herbicide in equilibrium with sorbed phase herbicide was displaced by herbicide-free water, which was infiltrated into the column. Sorption coefficients, Kd, were obtained from a plot of total herbicide concentration in the soil versus water content in the region where the antecedent solution accumulated. DKN Kd values were approximately 2-3 times (average Kd = 0.71 L kg-1) greater using the unsaturated transient flow method as compared to the batch equilibration method in clay loam (Kd = 0.33 L kg-1), but similar for the two methods in sand (0.12 vs 0.09 L kg-1) soils. Dicamba Kd values were 3 times greater using the unsaturated transient flow method as compared to the batch equilibration method in the clay loam soil (0.38 vs 0.13 L kg-1), however, the Kd values were the same for the two methods in the sand (approximately 0.06 L kg-1). This demonstrates that to determine sorption coefficients for labile hydrophilic pesticides, an unsaturated transient flow method may be a suitable alternative to the batch method. In fact, it may be better in cases where transport models have overpredicted herbicide leaching when batch sorption coefficients have been used.     

Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.     

The reactive transport of trichloroethene is influenced by residence time and microbial numbers

The dechlorination rate in a flow-through porous matrix can be described by the species specific dechlorination rate observed in a liquid batch unless mass transport limitations prevail. This hypothesis was examined by comparing dechlorination rates in liquid batch with that in column experiments at various flow rates (3-9-12 cm day(-1)). Columns were loaded with an inoculated sand and eluted with a medium containing 1mM trichloroethene (TCE) for 247 days. Dechlorination in the column treatments increased with decreasing flow rate, illustrating the effect of the longer residence time. Zeroth order TCE or cis-DCE degradation rates were 4-7 folds larger in columns than in corresponding batch systems which could be explained by the higher measured Geobacter and Dehalococcoides numbers per unit pore volume in the columns. The microbial numbers also explained the variability in dechlorination rate among flow rate treatments marked by a large elution of the dechlorinating species' yield as flow increased. Stop flow events did not reveal mass transport limitations for dechlorination. We conclude that flow rate effects on reactive transport of TCE in this coarse sand are explained by residence time and by microbial transport and that mass transport limitations in this porous matrix are limited.     

Evaluating equilibrium and non-equilibrium transport of bromide and isoproturon in disturbed and undisturbed soil columns

In this study, displacement experiments of isoproturon were conducted in disturbed and undisturbed columns of a silty clay loam soil under similar rainfall intensities. Solute transport occurred under saturated conditions in the undisturbed soil and under unsaturated conditions in the sieved soil because of a greater bulk density of the compacted undisturbed soil compared to the sieved soil. The objective of this work was to determine transport characteristics of isoproturon relative to bromide tracer. Triplicate column experiments were performed with sieved (structure partially destroyed to simulate conventional tillage) and undisturbed (structure preserved) soils. Bromide experimental breakthrough curves were analyzed using convective-dispersive and dual-permeability (DP) models (HYDRUS-1D). Isoproturon breakthrough curves (BTCs) were analyzed using the DP model that considered either chemical equilibrium or non-equilibrium transport. The DP model described the bromide elution curves of the sieved soil columns well, whereas it overestimated the tailing of the bromide BTCs of the undisturbed soil columns. A higher degree of physical non-equilibrium was found in the undisturbed soil, where 56% of total water was contained in the slow-flow matrix, compared to 26% in the sieved soil. Isoproturon BTCs were best described in both sieved and undisturbed soil columns using the DP model combined with the chemical non-equilibrium. Higher degradation rates were obtained in the transport experiments than in batch studies, for both soils. This was likely caused by hysteresis in sorption of isoproturon. However, it cannot be ruled out that higher degradation rates were due, at least in part, to the adopted first-order model. Results showed that for similar rainfall intensity, physical and chemical non-equilibrium were greater in the saturated undisturbed soil than in the unsaturated sieved soil. Results also suggested faster transport of isoproturon in the undisturbed soil due to higher preferential flow and lower fraction of equilibrium sorption sites.     

Sorption of Isoxaflutole Diketonitrile Degradate (DKN) and Dicamba in Unsaturated Soil

When analyzing the sorption characteristics of weakly sorbing or labile pesticides, batch methods tend to yield a high margin of error attributable to errors in concentration measurement and to degradation, respectively. This study employs a recently developed unsaturated transient flow method to determine the sorption of isoxaflutole's herbicidally active diketonitrile degradate (DKN) and dicamba. A 20-cm acrylic column was packed with soils with varied texture that had been uniformly treated with ¹⁴C-labeled chemical.

The antecedent solution herbicide in equilibrium with sorbed phase herbicide was displaced by herbicide-free water, which was infiltrated into the column. Sorption coefficients, K_d, were obtained from a plot of total herbicide concentration in the soil versus water content in the region where the antecedent solution accumulated. DKN K_d values were ～2–3 times (average K_d = 0.71 L kg^?1) greater using the unsaturated transient flow method as compared to the batch equilibration method in clay loam (K_d = 0.33 L kg^?1), but similar for the two methods in sand (0.12 vs 0.09 L kg^?1) soils. Dicamba K_d values were 3 times greater using the unsaturated transient flow method as compared to the batch equilibration method in the clay loam soil (0.38 vs 0.13 L kg^?1), however, the K_d values were the same for the two methods in the sand (～0.06 L kg^?1). This demonstrates that to determine sorption coefficients for labile hydrophilic pesticides, an unsaturated transient flow method may be a suitable alternative to the batch method. In fact, it may be better in cases where transport models have overpredicted herbicide leaching when batch sorption coefficients have been used.     

Numerical simulations of pyrite oxidation and acid mine drainage in unsaturated waste rock piles

Numerical simulations of layered, sulphide-bearing unsaturated waste rock piles are presented to illustrate the effect of coupled processes on the generation of acid mine drainage (AMD). The conceptual 2D systems were simulated using the HYDRUS model for flow and the POLYMIN model for reactive transport. The simulations generated low-pH AMD which was buffered by sequential mineral dissolution and precipitation. Sulphide oxidation rates throughout the pile varied by about two orders of magnitude (0.004-0.4 kg m-3 year-1) due to small changes in moisture content and grain size. In the fine-grained layers, the high reactive surface area induced high oxidation rates, even though capillary forces kept the local moisture content relatively high. In waste rock piles with horizontal layers, most of the acidity discharged through vertical preferential flow channels while with inclined fine grained layers, capillary diversion channeled the AMD to the outer slope boundary, keeping the pile interior relatively dry. The simulation approach will be useful for helping evaluate design strategies for controlling AMD from waste rock.     

Experimental Zn(II) retention in a sandy loam soil by very small columns

The use of column experiments, usually performed to better approximate field conditions, may provide information that is not available from batch experiments. In such experiments heavy metals are often adsorbed until saturation followed by desorption experiments. When the affinity of the metal to soil is high, the retention factor (R) could be greater than thousands and the duration of experiments can become impractically long. In order to use reasonable laboratory time, the flow rate should be increased or the column size decreased. The increase in flow rate produces undesirable kinetic and dispersion effects, so we used very small soil columns (pore volume = 0.31–0.70 ml) and relatively high flow rates (0.03–0.12 ml min⁻¹) in studies of Zn(II) adsorption and retention in soils. Conservative tracer flow column experiments under saturation conditions were carried out to determine flow parameters for different flow rates. Column pore volume (V_p), Peclet numbers (Pe) and longitudinal dispersion coefficients (D_L) were determined from breakthrough curves. The effect of type of electrolyte and ionic strength on the Zn(II) retention onto soil was determined. The influence of flow rate and bed height on the retention coefficient and on the mass transfer zone was also studied. The effect of different influent Zn(II) concentrations on the R values obtained was analyzed. Freundlich parameters from column experiments were compared with batch ones. The leaching efficiency of different electrolytes, salts of weak organic acids and EDTA was also studied.     

Experimental study and modeling of the transfer of zinc in a low reactive sand column in the presence of acetate

Nowadays, it is necessary to understand and identify the reactions governing the fate of heavy metals introduced into the environment with low complexing organic compounds, particularly when they are transferred through soils in urban areas. In this work the concomitant influence of pH and acetate on the fate of zinc on siliceous sand was studied in batch and non-saturated column experiments. Total zinc concentrations varied between 2 and 20 mg/l, and total acetate concentrations were fixed at 22, 72, 132, and 223 mM to obtain solution pHs of 4, 5, 6 and 7, respectively. Natural sand (diameter, 0.3-2 mm), mainly constituted of silica, was used. In batch adsorption experiments, zinc adsorption is insignificant at pH 4, low and linear at pH 5, and increasingly nonlinear, of the Langmuir type, at pH 6 and 7 indicating near-saturation conditions of surface sites at these high pH values. In column experiments, Zn retardation increases and the maximum outlet concentration of Zn decreases with rising pH and acetate concentrations. Previous column tracer experiments revealed the occurrence of regionalized water transport in the column. Modeling these data was based on a non-electrostatic approach. Batch and column data modeling was based on the PHREEQC code that allows concomitant resolution of chemical speciation and regionalized water transport. The speciation calculation indicates that the ZnAcetate+ species is the dominant Zn species in the solutions used. Batch experimental curves are correctly modeled assuming the formation of the three surface species triple bond SiOZn+, triple bond SiOH-Zn Acetate+ and triple bond SiO-Zn(Acetate)2-. The column data could be adequately modeled assuming a two-region water transport and the formation of the same three species with the same thermodynamic constants determined in the batch experiments. The hypothesis of the modeling leads to a slight overestimation of the quantities of zinc eluted (10%) at pH 6 and 7, mostly in the desorption phase. These results show that the methodology used facilitates the correct modeling of both batch and transport experiments and formulation of the hypothesis on the interactions between the low reactive sand and a complex solution.     

Kinetic limitations on tracer partitioning in ganglia dominated source zones

Quantification of the relationship between dense nonaqueous phase liquid (DNAPL) source strength, source longevity and spatial distribution is increasingly recognized as important for effective remedial design. Partitioning tracers are one tool that may permit interrogation of DNAPL architecture. Tracer data are commonly analyzed under the assumption of linear, equilibrium partitioning, although the appropriateness of these assumptions has not been fully explored. Here we focus on elucidating the nonlinear and nonequilibrium partitioning behavior of three selected alcohol tracers - 1-pentanol, 1-hexanol and 2-octanol in a series of batch and column experiments. Liquid-liquid equilibria for systems comprising water, TCE and the selected alcohol illustrate the nonlinear distribution of alcohol between the aqueous and organic phases. Complete quantification of these equilibria facilitates delineation of the limits of applicability of the linear partitioning assumption, and assessment of potential inaccuracies associated with measurement of partition coefficients at a single concentration. Column experiments were conducted under conditions of non-equilibrium to evaluate the kinetics of the reversible absorption of the selected tracers in a sandy medium containing a uniform entrapped saturation of TCE-DNAPL. Experimental tracer breakthrough data were used, in conjunction with mathematical models and batch measurements, to evaluate alternative hypotheses for observed deviations from linear equilibrium partitioning behavior. Analyses suggest that, although all tracers accumulate at the TCE-DNAPL/aqueous interface, surface accumulation does not influence transport at concentrations typically employed for tracer tests. Moreover, results reveal that the kinetics of the reversible absorption process are well described using existing mass transfer correlations originally developed to model aqueous boundary layer resistance for pure-component NAPL dissolution.     

Colloid transport in unsaturated porous media: the role of water content and ionic strength on particle straining

Packed column and mathematical modeling studies were conducted to explore the influence of water saturation, pore-water ionic strength, and grain size on the transport of latex microspheres (1.1 microm) in porous media. Experiments were carried out under chemically unfavorable conditions for colloid attachment to both solid-water interfaces (SWI) and air-water interfaces (AWI) using negatively charged and hydrophilic colloids and modifying the solution chemistry with a bicarbonate buffer to pH 10. Interaction energy calculations and complementary batch experiments were conducted and demonstrated that partitioning of colloids to the SWI and AWI was insignificant across the range of the ionic strengths considered. The breakthrough curve and final deposition profile were measured in each experiment indicating colloid retention was highly dependent on the suspension ionic strength, water content, and sand grain size. In contrast to conventional filtration theory, most colloids were found deposited close to the column inlet, and hyper-exponential deposition profiles were observed. A mathematical model, accounting for time- and depth-dependent straining, produced a reasonably good fit for both the breakthrough curves and final deposition profiles. Experimental and modeling results suggest that straining--the retention of colloids in low velocity regions of porous media such as grain junctions--was the primary mechanism of colloid retention under both saturated and unsaturated conditions. The extent of stagnant regions of flow within the pore structure is enhanced with decreasing water content, leading to a greater amount of retention. Ionic strength also contributes to straining, because the number of colloids that are held in the secondary energy minimum increases with ionic strength. These weakly associated colloids are prone to be translated to stagnation regions formed at grain-grain junctions, the solid-water-air triple point, and dead-end pores and then becoming trapped.     

Effects of degree of water saturation on dispersivity and immobile water in sandy soil columns

Three natural nonaggregated soil samples, with similar grain-size distributions, have been used to determine the dispersive behavior of porous media under steady, saturated and unsaturated flow conditions. Tritium was used as a tracer and was found to have no sorption on the solid matrix. Generated breakthrough curves (BTCs) for the unsaturated experiments were symmetrical with no evidence of tailing. The unsaturated experiments for two of the soils were adequately described by considering all the water in the pore volume as mobile. However, about 10% of the pore water, independent of the degree of saturation, was found to be immobile in the case of the third soil during unsaturated flow. For this soil, there was no mass transfer between the two water regions, indicating that the immobile water is essentially isolated from the flowing water fraction. For all three soils, dispersivity under unsaturated conditions was found to be higher, independent of the degree of water saturation, than the value determined for the saturated experiments. This is inconsistent with what would be expected from the simple bundle-of-capillary-tubes model and does not agree well with a more sophisticated conceptualization of the porous medium. The data, however, clearly indicate a wider range in pore-water velocities when these soils are desaturated.     

Thioarsenate formation upon dissolution of orpiment and arsenopyrite

Thioarsenates were previously determined as dominant species in geothermal and mineral waters with excess sulfide. Here, we used batch leaching experiments to determine their formation upon weathering or industrial leaching of the arsenic-sulfide minerals orpiment (As₂S₃) and arsenopyrite (FeAsS) under different pH and oxygen conditions. Under acidic conditions, as expected based on their known kinetic instability at low pH, no thioarsenates formed in either of the two mineral systems. Under neutral to alkaline conditions, orpiment dissolution yielded mono-, di- and trithioarsenate which accounted for up to 43-55% of total arsenic. Thioarsenate formation upon arsenopyrite dissolution was low at neutral (4%) but significant at alkaline pH, especially under suboxic to sulfidic conditions (20-43%, mainly as monothioarsenate). In contrast to orpiment, we postulate that recombination of arsenite and sulfide in solution is of minor importance for monothioarsenate formation during alkaline arsenopyrite dissolution. We propose instead that hydroxyl physisorption lead to formation of As-OH-S surface complexes by transposition of hydroxyl anions to arsenic or iron sites. Concurrently formed ironhydroxides could provide re-sorption sites for the freshly released monothioarsenate. However, sorption experiments with goethite showed slower sorption kinetics of monothioarsenate compared to arsenite, but comparable with arsenate. The discovery that thioarsenates are released by natural weathering and industrial leaching processes and that, once they are released, have a higher mobility than the commonly-investigated species arsenite and arsenate requires future studies to consider them when assessing arsenic release in sulfidic natural or mining-impacted environments.     

Preparation of Calcium Silicate Reagent from Fly Ash and Lime in a Flow Reactor

Abstract

The batch reaction between fly ash and hydrated lime in water to produce high surface area calcium silicates for flue gas desulfurization has been examined extensively. This paper examines the reaction in a flow reactor using two low-calcium fly ashes and introducing gypsum, calcium sulf ite hemihydrate, and calcium chloride as additives to the reaction. The flow system is compared to the batch reaction at similar operating conditions and a segregated flow model is used to approximate flow reactor behavior.

Experiments with calcium chloride and gypsum additives were modeled fairly well by the segregated flow approximation at residence times less than 12 hours. The flow reactor produced low surface area material at longer residence times when gypsum was present. Because the changing solution chemistry affected the batch reaction rate the fly ash and hydrated lime system without gypsum or calcium chloride could not be approximated using batch reaction data. In this case, the flow reactor produced higher surface area product than the batch reactor for a given residence time due to the increased calcium hydroxide availability.     

Using aliphatic alcohols as gaseous tracers in determination of water contents and air-water interfacial areas in unsaturated sands

A new type of gaseous tracer utilizing nontoxic aliphatic alcohols for the determination of water content and air-water interfacial area is tested on unsaturated sands of low water content. Alcohol vapors are generated at room temperature and passed through the experimental sand column. Breakthrough curves (BTCs) of these vapors are obtained by monitoring their effluent concentrations using GC-FID. The retardation factor with respect to each vapor transport process is obtained by optimizing BTCs data using the CXTFIT program in the reverse problem mode. The water content and the interfacial area are subsequently calculated from their retardation factors by both equilibrium and nonequilibrium transport models. Experimental results indicate that the pentanol tracer is feasible in the determination of water content at conditions when the degree of water saturation is low. In the determination of air-water interfacial area, decanol is selected due to its interfacial adsorption characteristics. By comparing to interfacial areas from theoretical predictions as well as other conventional tarcer methods, the ones determined from the decanol tracer tests are found to be close to the true interfacial areas when the water content is low.     

Thioarsenate formation upon dissolution of orpiment and arsenopyrite

Thioarsenates were previously determined as dominant species in geothermal and mineral waters with excess sulfide. Here, we used batch leaching experiments to determine their formation upon weathering or industrial leaching of the arsenic-sulfide minerals orpiment (As₂S₃) and arsenopyrite (FeAsS) under different pH and oxygen conditions. Under acidic conditions, as expected based on their known kinetic instability at low pH, no thioarsenates formed in either of the two mineral systems. Under neutral to alkaline conditions, orpiment dissolution yielded mono-, di- and trithioarsenate which accounted for up to 43–55% of total arsenic. Thioarsenate formation upon arsenopyrite dissolution was low at neutral (4%) but significant at alkaline pH, especially under suboxic to sulfidic conditions (20–43%, mainly as monothioarsenate). In contrast to orpiment, we postulate that recombination of arsenite and sulfide in solution is of minor importance for monothioarsenate formation during alkaline arsenopyrite dissolution. We propose instead that hydroxyl physisorption lead to formation of As-OH-S surface complexes by transposition of hydroxyl anions to arsenic or iron sites. Concurrently formed ironhydroxides could provide re-sorption sites for the freshly released monothioarsenate. However, sorption experiments with goethite showed slower sorption kinetics of monothioarsenate compared to arsenite, but comparable with arsenate. The discovery that thioarsenates are released by natural weathering and industrial leaching processes and that, once they are released, have a higher mobility than the commonly-investigated species arsenite and arsenate requires future studies to consider them when assessing arsenic release in sulfidic natural or mining-impacted environments.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.