A predictive model (ETLM) for arsenate adsorption and surface speciation on oxides consistent with spectroscopic and theoretical molecular evidence

The nature of adsorbed arsenate species for a wide range of minerals and environmental conditions is fundamental to prediction of the migration and long-term fate of arsenate in natural environments. Spectroscopic experiments and theoretical calculations have demonstrated the potential importance of a variety of arsenate surface species on several iron and aluminum oxides. However, integration of the results of these studies with surface complexation models and extrapolation over wide ranges of conditions and for many oxides remains a challenge. In the present study, in situ X-ray and infrared spectroscopic and theoretical molecular evidence of arsenate (and the analogous phosphate) surface speciation are integrated with an extended triple layer model (ETLM) of surface complexation, which takes into account the electrostatic work associated with the ions and the water dipoles involved in inner-sphere surface complexation by the ligand exchange mechanism.Three reactions forming inner-sphere arsenate surface species

     

A predictive model (ETLM) for As(III) adsorption and surface speciation on oxides consistent with spectroscopic data

Arsenic(III) adsorption reactions are thought to play a critical role in the mobility of arsenic in the environment. It is the nature of the As(III) surface species that must be known on a wide variety of minerals and over a range of pH, ionic strength and surface coverage in order to be able to predict adsorption behavior. EXAFS and XANES spectroscopic studies have identified bidentate, binuclear inner-sphere surface species and/or an outer-sphere species, but only a few oxides have been examined. These results need to be integrated with a predictive surface complexation model in order to ascertain the environmental conditions under which the different surface species may be important on a wide range of solids. In the present study, the surface species information from XAFS and XANES studies has been built into a recent extension of the triple-layer model (ETLM) for the formation of inner-sphere complexes of anions that takes into account the electrostatics of water dipole desorption during ligand exchange reactions. The ETLM has been applied to regress surface titration, proton coadsorption, and As(III) adsorption data over extensive ranges of pH, ionic strength, electrolyte type and surface coverage for magnetite, goethite, gibbsite, amorphous hydrous alumina, hydrous ferric oxide (HFO), ferrihydrite, and amorphous iron oxide. Two principal reactions forming inner- and outer-sphere As(III) surface species,

     

Speciation of adsorbed yttrium and rare earth elements on oxide surfaces

The distribution of yttrium and the rare earth elements (YREE) between natural waters and oxide mineral surfaces depends on adsorption reactions, which in turn depend on the specific way in which YREE are coordinated to mineral surfaces. Recent X-ray studies have established that Y³⁺ is adsorbed to the rutile (1 1 0) surface as a distinctive tetranuclear species. However, the hydrolysis state of the adsorbed cation is not known from experiment. Previous surface complexation models of YREE adsorption have suggested two to four cation hydrolysis states coexisting on oxide surfaces. In the present study, we investigate the applicability of the X-ray results to rare earth elements and to several oxides in addition to rutile using the extended triple-layer surface complexation model. The reaction producing a hydrolyzed tetranuclear surface species

     

A spectrophotometric study of Nd(III), Sm(III) and Er(III) complexation in sulfate-bearing solutions at elevated temperatures

The speciation of Nd(III), Sm(III), and Er(III) in sulfate-bearing solutions has been determined spectrophotometrically at temperatures from 25 to 250 °C and a pressure of 100 bars. The data obtained earlier on the speciation of Nd in sulfate-bearing solutions (Migdisov et al., 2006) have been re-evaluated and corrected using a more appropriate activity model and are compared with the corresponding data for Sm(III) and Er(III) and new data for Nd(III). Based on this comparison, the dominant species in the solution are interpreted to be and , with the latter complex increasing in importance at higher temperature. Equilibrium constants were calculated for the following reactions:

     

A spectrophotometric study of samarium (III) speciation in chloride solutions at elevated temperatures

The speciation of samarium (III) in chloride-bearing solutions was investigated spectrophotometrically at temperatures of 100-250 °C and a pressure of 100 bars. The simple hydrated ion, Sm³⁺, is predominant at ambient temperature, but chloride complexes are the dominant species at elevated temperatures. Cumulative formation constants for samarium chloride species were calculated for the following reactions:

     

Prediction of iodate adsorption and surface speciation on oxides by surface complexation modeling

Aqueous iodine species occur mainly as iodide (I⁻) and iodate (IO₃⁻), depending on redox conditions. The adsorption of IO₃⁻ on naturally occurring oxides under oxic conditions is of environmental concern. The adsorption behaviors of IO₃⁻ by hydrous ferric oxide (HFO), α-FeOOH, and γ-Al₂O₃ were examined in this study as functions of pH, ionic strength, and solid concentration. Adsorption data were analyzed using an extended triple-layer model (ETLM) for surface complexation modeling to infer IO₃⁻ adsorption reactions and equilibrium constants. Results of ETLM analysis suggest that adsorption of IO₃⁻ is both an outer-sphere and an inner-sphere process, as expressed by the following complexation reactions, which are consistent with the independent pressure jump kinetic results and adsorption enthalpy measurements

     

Goethite surface reactivity: III. Unifying arsenate adsorption behavior through a variable crystal face - Site density model

We developed a model that describes quantitatively the arsenate adsorption behavior for any goethite preparation as a function of pH and ionic strength, by using one basic surface arsenate stoichiometry, with two affinity constants. The model combines a face distribution-crystallographic site density model for goethite with tenets of the Triple Layer and CD-MUSIC surface complexation models, and is self-consistent with its adsorption behavior towards protons, electrolytes, and other ions investigated previously. Five different systems of published arsenate adsorption data were used to calibrate the model spanning a wide range of chemical conditions, which included adsorption isotherms at different pH values, and adsorption pH-edges at different As(V) loadings, both at different ionic strengths and background electrolytes. Four additional goethite-arsenate systems reported with limited characterization and adsorption data were accurately described by the model developed. The adsorption reaction proposed is:

     

Inhibition of calcite precipitation by orthophosphate: Speciation and thermodynamic considerations

The inhibition of heterogeneous calcite precipitation by orthophosphate was investigated under four different solution compositions using a pH-stat system. The system composition was designed to maintain a constant degree of supersaturation with respect to calcite, but with different carbonate/calcium ratios and pH values during precipitation. Inhibition in the presence of orthophosphate was found to be more effective at lower carbonate/calcium ratios and lower pH values. With the assumption that the calcite precipitation rate is proportional to the surface concentration of active crystal-growth sites, the reduction in the rate of calcite precipitation by phosphate can be explained by a Langmuir adsorption model using a conditional equilibrium constant and total phosphate concentration. Through a detailed analysis of chemical speciation in the solution phase and calcite surface speciation using chemical equilibrium computer modeling, the “conditional” equilibrium constants obtained at different solution compositions were found to converge to a single “non-conditional” value if only was considered in the adsorption reaction. This suggests that is the responsible species for inhibition of calcite precipitation because it adsorbs to the surface and blocks the active crystal-growth sites. The standard enthalpy change (ΔH⁰) and standard entropy change (TΔS⁰) of the adsorption reaction, determined by experiments performed from 15 to 45 °C, were 58.5 and 98.3 kJ/mol, respectively. The high positive values of the standard enthalpy change and the standard entropy change suggest that the adsorption reaction is an endothermic reaction, chemisorptive in nature, and driven by the entropy change, most likely resulting from the dehydration process that accompanies the adsorption of onto the calcite surface.     

A spectrophotometric study of neodymium(III) complexation in sulfate solutions at elevated temperatures

The formation constants of neodymium complexes in sulfate solutions have been determined spectrophotometrically at temperatures of 30-250 °C and a pressure of 100 bars. The dominant species in the solution are NdSO₄⁺ and Nd(SO₄)₂⁻, with the latter complex being more important at higher temperature. Equilibrium constants were calculated for the following reactions:

     

Comparative carbonate and hydroxide complexation of cations in seawater

Comparative concentrations of carbonate and hydroxide complexes in natural solutions can be expressed in terms of reactions with bicarbonate that have no explicit pH dependence (). Stability constants for this reaction with n = 1 were determined using conventional formation constant data expressed in terms of hydroxide and carbonate. Available data indicate that stability constants appropriate to seawater at 25 °C expressed in the form are on the order of 10^4.2 for a wide range of cations (M^z+) with z = +1, +2 and +3. Φ₁ is sufficiently large that species appear to substantially dominate MOH^z−1 species in seawater. Evaluations of comparative stepwise carbonate and hydroxide stability constant behavior leading to the formation of n = 2 and n = 3 complexes suggest that carbonate complexes generally dominate hydroxide complexes in seawater, even for cations whose inorganic speciation schemes in seawater are currently presumed to be strongly dominated by hydrolyzed forms (). Calculated stability constants, and , indicate that the importance of carbonate complexation is sufficiently large that carbonate and hydroxide complexes would be generally comparable even if calculated Φ₂ and Φ₃ values are overestimated by two or more orders of magnitude. Inclusion of mixed ligand species in carbonate-hydroxide speciation models allows cation complexation intensities (M_T/[M^z+]) to be expressed in the following form:

     

Oxygen isotope biogeochemistry of pore water sulfate in the deep biosphere: Dominance of isotope exchange reactions with ambient water during microbial sulfate reduction (ODP Site 1130)

Microbially mediated sulfate reduction affects the isotopic composition of dissolved and solid sulfur species in marine sediments. Experiments and field data show that the composition is also modified in the presence of sulfate-reducing microorganisms. This has been attributed either to a kinetic isotope effect during the reduction of sulfate to sulfite, cell-internal exchange reactions between enzymatically-activated sulfate (APS), and/or sulfite with cytoplasmic water. The isotopic fingerprint of these processes may be further modified by the cell-external reoxidation of sulfide to elemental sulfur, and the subsequent disproportionation to sulfide and sulfate or by the oxidation of sulfite to sulfate. Here we report values from interstitial water samples of ODP Leg 182 (Site 1130) and provide the mathematical framework to describe the oxygen isotope fractionation of sulfate during microbial sulfate reduction. We show that a purely kinetic model is unable to explain our data, and that the data are well explained by a model using oxygen isotope exchange reactions. We propose that the oxygen isotope exchange occurs between APS and cytoplasmic water, and/or between sulfite and adenosine monophosphate (AMP) during APS formation. Model calculations show that cell external reoxidation of reduced sulfur species would require up to 3000 mol/m³ of an oxidant at ODP Site 1130, which is incompatible with the sediment geochemical data. In addition, we show that the volumetric fluxes required to explain the observed data are on average 14 times higher than the volumetric sulfate reduction rates (SRR) obtained from inverse modeling of the porewater data. The ratio between the gross sulfate flux into the microbes and the net sulfate flux through the microbes is depth invariant, and independent of sulfide concentrations. This suggests that both fluxes are controlled by cell density and that cell-specific sulfate reduction rates remain constant with depth.     

A potentiometric study of the stability of aqueous yttrium-acetate complexes from 25 to 175 °C and 1-1000 bar

The stability of yttrium-acetate (Y-Ac) complexes in aqueous solution was determined potentiometrically at temperatures 25-175 °C (at P_s) and pressures 1-1000 bar (at 25 and 75 °C). Measurements were performed using glass H⁺-selective electrodes in potentiometric cells with a liquid junction. The species YAc²⁺ and were found to dominate yttrium aqueous speciation in experimental solutions at 25-100 °C (log [Ac⁻] < −1.5, pH < 5.2), whereas at 125, 150 and 175 °C introduction of into the Y-Ac speciation model was necessary. The overall stability constants β_n were determined for the reaction

     

Fluoride geochemistry of thermal waters in Yellowstone National Park: I. Aqueous fluoride speciation

Thermal water samples from Yellowstone National Park (YNP) have a wide range of pH (1-10), temperature, and high concentrations of fluoride (up to 50 mg/l). High fluoride concentrations are found in waters with field pH higher than 6 (except those in Crater Hills) and temperatures higher than 50 °C based on data from more than 750 water samples covering most thermal areas in YNP from 1975 to 2008. In this study, more than 140 water samples from YNP collected in 2006-2009 were analyzed for free-fluoride activity by ion-selective electrode (ISE) method as an independent check on the reliability of fluoride speciation calculations. The free to total fluoride concentration ratio ranged from <1% at low pH values to >99% at high pH. The wide range in fluoride activity can be explained by strong complexing with H⁺ and Al³⁺ under acidic conditions and lack of complexing under basic conditions. Differences between the free-fluoride activities calculated with the WATEQ4F code and those measured by ISE were within 0.3-30% for more than 90% of samples at or above 10⁻⁶ molar, providing corroboration for chemical speciation models for a wide range of pH and chemistry of YNP thermal waters. Calculated speciation results show that free fluoride, F⁻, and major complexes (, AlF²⁺, and ) account for more than 95% of total fluoride. Occasionally, some complex species like , FeF²⁺, , MgF⁺ and may comprise 1-10% when the concentrations of the appropriate components are high. According to the simulation results by PHREEQC and calculated results, the ratio of main fluoride species to total fluoride varies as a function of pH and the concentrations and ratios of F and Al.     

Theoretical study on the reactivity of sulfate species with hydrocarbons

The abiotic, thermochemically controlled reduction of sulfate to hydrogen sulfide coupled with the oxidation of hydrocarbons, is termed thermochemical sulfate reduction (TSR), and is an important alteration process that affects petroleum accumulations in nature. Although TSR is commonly observed in high-temperature carbonate reservoirs, it has proven difficult to simulate in the laboratory under conditions resembling nature. The present study was designed to evaluate the relative reactivities of various sulfate species in order to provide greater insight into the mechanism of TSR and potentially to fill the gap between laboratory experimental data and geological observations. Accordingly, quantum mechanics density functional theory (DFT) was used to determine the activation energy required to reach a potential transition state for various aqueous systems involving simple hydrocarbons and different sulfate species. The entire reaction process that results in the reduction of sulfate to sulfide is far too complex to be modeled entirely; therefore, we examined what is believed to be the rate limiting step, namely, the reduction of sulfate S(VI) to sulfite S(IV). The results of the study show that water-solvated sulfate anions are very stable due to their symmetrical molecular structure and spherical electronic distributions. Consequently, in the absence of catalysis, the reactivity of is expected to be extremely low. However, both the protonation of sulfate to form bisulfate anions () and the formation of metal-sulfate contact ion-pairs could effectively destabilize the sulfate molecular structure, thereby making it more reactive.Previous reports of experimental simulations of TSR generally have involved the use of acidic solutions that contain elevated concentrations of relative to . However, in formation waters typically encountered in petroleum reservoirs, the concentration of is likely to be significantly lower than the levels used in the laboratory, with most of the dissolved sulfate occurring as , aqueous calcium sulfate ([CaSO₄]_(aq)), and aqueous magnesium sulfate ([MgSO₄]_(aq)). Our calculations indicate that TSR reactions that occur in natural environments are most likely to involve bisulfate ions () and/or magnesium sulfate contact ion-pairs ([MgSO₄]_CIP) rather than ‘free’ sulfate ions () or solvated sulfate ion-pairs, and that water chemistry likely plays a significant role in controlling the rate of TSR.     

Flow and nutrient dynamics in a subterranean estuary (Waquoit Bay, MA, USA): Field data and reactive transport modeling

A two-dimensional (2D) reactive transport model is used to investigate the controls on nutrient (, , PO₄) dynamics in a coastal aquifer. The model couples density-dependent flow to a reaction network which includes oxic degradation of organic matter, denitrification, iron oxide reduction, nitrification, Fe²⁺ oxidation and sorption of PO₄ onto iron oxides. Porewater measurements from a well transect at Waquoit Bay, MA, USA indicate the presence of a reducing plume with high Fe²⁺, , DOC (dissolved organic carbon) and PO₄ concentrations overlying a more oxidizing -rich plume. These two plumes travel nearly conservatively until they start to overlap in the intertidal coastal sediments prior to discharge into the bay. In this zone, the aeration of the surface beach sediments drives nitrification and allows the precipitation of iron oxide, which leads to the removal of PO₄ through sorption. Model simulations suggest that removal of through denitrification is inhibited by the limited overlap between the two freshwater plumes, as well as by the refractory nature of terrestrial DOC. Submarine groundwater discharge is a significant source of to the bay.     

Selenite retention by nanocrystalline magnetite: Role of adsorption, reduction and dissolution/co-precipitation processes

We studied selenite () retention by magnetite () using both surface complexation modeling and X-ray absorption spectroscopy (XAS) to characterize the processes of adsorption, reduction, and dissolution/co-precipitation. The experimental sorption results for magnetite were compared to those of goethite (Fe^IIIOOH) under similar conditions. Selenite sorption was investigated under both oxic and anoxic conditions and as a function of pH, ionic strength, solid-to-liquid ratio and Se concentration. Sorption onto both oxides was independent of ionic strength and decreased as pH increased, as expected for anion sorption; however, the shape of the sorption edges was different. The goethite sorption data could be modeled assuming the formation of an inner-sphere complex with iron oxide surface sites (SOH). In contrast, the magnetite sorption data at low pH could be modeled only when the dissolution of magnetite, the formation of aqueous iron-selenite species, and the subsequent surface complexation of these species were implemented. The precipitation of ferric selenite was the predominant retention process at higher selenite concentrations (>1 × 10⁻⁴ M) and pH < 5, which was in agreement with the XAS results. Sorption behavior onto magnetite was similar under oxic and anoxic conditions. Under anoxic conditions, we did not observe the reduction of selenite. Possible reasons for the absence of reduction are discussed. In conclusion, we show that under acidic reaction conditions, selenite retention by magnetite is largely influenced by dissolution and co-precipitation processes.     

The ionic strength dependence of lead (II) carbonate complexation in perchlorate media

Lead speciation in many aqueous geochemical systems is dominated by carbonate complexation. However, direct observations of Pb²⁺ complexation by carbonate ions are few in number. This work represents the first investigation of the equilibrium over a range of ionic strength. Through spectrophotometric observations of formation at 25 °C in NaHCO₃-NaClO₄ solutions, formation constants of the form were determined between 0.001 and 5.0 molal ionic strength. Formation constant results were well represented by the equation:

     

Kinetics of gypsum crystal growth from high ionic strength solutions: A case study of Dead Sea - seawater mixtures

Gypsum precipitation kinetics were examined from a wide range of chemical compositions , ionic strengths (4.75-10 m) and saturation state with respect to gypsum (1.16-1.74) in seeded batch experiments of mixtures of Ca²⁺-rich Dead Sea brine and -rich seawater. Despite the variability in the experimental solutions, a single general rate law was formulated to describe the heterogeneous precipitation rate of gypsum from these mixtures:

     

Biogenic hydroxysulfate green rust, a potential electron acceptor for SRB activity

Microbiological reduction of a biogenic sulfated green rust , was examined using a sulfate reducing bacterium (Desulfovibrio alaskensis). Experiments investigated whether could serve as a sulfate source for D. alaskensis anaerobic respiration by analyzing mineral transformation. Batch experiments were conducted using lactate as the electron donor and biogenic as the electron acceptor, at circumneutral pH in unbuffered medium. transformation was monitored with time by X-ray diffraction (XRD), Transmission Mössbauer Spectroscopy (TMS), Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS), Transmission Electron Microscopy (TEM) and X-ray Photoelectron Spectroscopy (XPS). The reduction of sulfate anions and the formation of iron sulfur mineral were clearly identified by XPS analyses. TMS showed the formation of additional mineral as green rust (GR) and vivianite. XRD analyses discriminated the type of the newly formed GR as GR1. The formed GR1 was as indicated by DRIFTS analysis. Thus, the results presented in this study indicate that D. alaskensis cells were able to use as an electron acceptor. , vivianite and an iron sulfur compound were formed as a result of reduction by D. alaskensis. Hence, in environments where geochemical conditions promote biogenic formation, this mineral could stimulate the anaerobic respiration of sulfate reducing bacteria.     

Kinetics of organic matter degradation, microbial methane generation, and gas hydrate formation in anoxic marine sediments

Seven sediment cores were taken in the Sea of Okhotsk in a south-north transect along the slope of Sakhalin Island. The retrieved anoxic sediments and pore fluids were analyzed for particulate organic carbon (POC), total nitrogen, total sulfur, dissolved sulfate, sulfide, methane, ammonium, iodide, bromide, calcium, and total alkalinity. A novel method was developed to derive sedimentation rates from a steady-state nitrogen mass balance. Rates of organic matter degradation, sulfate reduction, methane turnover, and carbonate precipitation were derived from the data applying a steady-state transport-reaction model. A good fit to the data set was obtained using the following new rate law for organic matter degradation in anoxic sediments: