Dawsonite synthesis and reevaluation of its thermodynamic properties from solubility measurements: Implications for mineral trapping of CO2

Over the last decade, a significant research effort has focused on determining the feasibility of sequestering large amounts of CO₂ in deep, permeable geologic formations to reduce carbon dioxide emissions to the atmosphere. Most models indicate that injection of CO₂ into deep sedimentary formations will lead to the formation of various carbonate minerals, including the common phases calcite (CaCO₃), dolomite (CaMg(CO₃)₂), magnesite (MgCO₃), siderite (FeCO₃), as well as the far less common mineral, dawsonite (NaAlCO₃(OH)₂). Nevertheless, the equilibrium and kinetics that control the precipitation of stable carbonate minerals are poorly understood and few experiments have been performed to validate computer codes that model CO₂ sequestration.In order to reduce this uncertainty we measured the solubility of synthetic dawsonite according to the equilibrium: , from under- and oversaturated solutions at 50-200 °C in basic media at 1.0 mol · kg⁻¹ NaCl. The solubility products (Q_s) obtained were extrapolated to infinite dilution to obtain the solubility constants (. Combining the fit of these values and fixing  at 25 °C, which was derived from the calorimetric data of Ferrante et al. [Ferrante, M.J., Stuve, J.M., and Richardson, D.W., 1976. Thermodynamic data for synthetic dawsonite. U.S. Bureau of Mines Report Investigation, 8129, Washington, D.C., 13p.], the following thermodynamic parameters for the dissolution of dawsonite were calculated at 25 °C: , and . Subsequently, we were able to derive values for the Gibbs energy of formation (, enthalpy of formation ( and entropy ( of dawsonite. These results are within the combined experimental uncertainties of the values reported by Ferrante et al. (1976). Predominance diagrams are presented for the dawsonite/boehmite and dawsonite/bayerite equilibria at 100 °C in the presence of a saline solution with and without silica-containing minerals.     

Sorption of phosphate onto calcite; results from batch experiments and surface complexation modeling

The adsorption of phosphate onto calcite was studied in a series of batch experiments. To avoid the precipitation of phosphate-containing minerals the experiments were conducted using a short reaction time (3 h) and low concentrations of phosphate (?50 μM). Sorption of phosphate on calcite was studied in 11 different calcite-equilibrated solutions that varied in pH, P_CO2, ionic strength and activity of Ca²⁺, and . Our results show strong sorption of phosphate onto calcite. The kinetics of phosphate sorption onto calcite are fast; adsorption is complete within 2-3 h while desorption is complete in less than 0.5 h. The reversibility of the sorption process indicates that phosphate is not incorporated into the calcite crystal lattice under our experimental conditions. Precipitation of phosphate-containing phases does not seem to take place in systems with ?50 μM total phosphate, in spite of a high degree of super-saturation with respect to hydroxyapatite (SI_HAP ? 7.83). The amount of phosphate adsorbed varied with the solution composition, in particular, adsorption increases as the activity decreases (at constant pH) and as pH increases (at constant activity). The primary effect of ionic strength on phosphate sorption onto calcite is its influence on the activity of the different aqueous phosphate species. The experimental results were modeled satisfactorily using the constant capacitance model with >CaPO₄Ca⁰ and either >CaHPO₄Ca⁺ or > as the adsorbed surface species. Generally the model captures the variation in phosphate adsorption onto calcite as a function of solution composition, though it was necessary to include two types of sorption sites (strong and weak) in the model to reproduce the convex shape of the sorption isotherms.     

Co-precipitation of radium in high ionic strength systems: 2. Kinetic and ionic strength effects

High concentrations of naturally occurring radium pose environmental and health concerns in natural and industrial systems. The adsorption of Ra²⁺ in saline water is limited compared to its adsorption in fresh water, but the process of co-precipitation may be effective in decreasing its concentration. However, despite its importance, Ra co-precipitation has rarely been studied in high ionic strength environments such as those in evaporitic systems.The fate of Ra in the reject brine of a desalination plant was studied via evaporation batch experiments at ionic strengths (I) ranging from 0.7 to 7.0 mol kg⁻¹. Precipitation sequences revealed that Ra co-precipitated with barite, even though the latter was a trace mineral compared to the precipitated gypsum. The concentration-based effective partition coefficient, , for the co-precipitation reaction was 1.04 ± 0.01. This value of is significantly lower than the value for relatively diluted solutions (1.8 ± 0.1). This low value of is mainly the result of a kinetic effect but is also slightly affected by the ionic strength.Both effects are quantitatively examined in the present paper. It is suggested that a kinetic effect influences the nucleation of (Ra,Ba)SO₄, reducing the value of the partition coefficient. This kinetic effect is caused by the favorable nucleation of a more soluble phase (i.e., a phase with a higher BaSO₄ fraction). An additional decrease in the partition coefficient results from the ionic strength effect. Considering the activity of Ra²⁺ and Ba²⁺ in the solution (rather than their concentration) makes it possible to determine the activity-based partition coefficient (), which accounts for the ionic strength effect. was calculated empirically from the experiments and theoretically via a kinetic model. The two derived values are consistent with one another and indicate the combined effect of ionic strength and precipitation kinetics.Finally, the common assumption that γ_Ra²⁺/γ_Ba²⁺=1 was re-examined using a numerical model to predict the experimental results. As the ionic strength increases, this assumption becomes less appropriate for predicting the change in as calculated in the experiments. Understanding the co-precipitation of Ra in such systems is crucial for risk assessments in which both Ra concentration and ionic strength are relatively high.     

Arsenate uptake by calcite: Macroscopic and spectroscopic characterization of adsorption and incorporation mechanisms

Batch uptake experiments and X-ray element mapping and spectroscopic techniques were used to investigate As(V) (arsenate) uptake mechanisms by calcite, including adsorption and coprecipitation. Batch sorption experiments in calcite-equilibrated suspensions (pH 8.3; PCO₂ = 10^−3.5 atm) reveal rapid initial sorption to calcite, with sorption rate gradually decreasing with time as available sorption sites decrease. An As(V)-calcite sorption isotherm determined after 24 h equilibration exhibits Langmuir-like behavior up to As concentrations of 300 μM. Maximum distribution coefficient values (K_d), derived from a best fit to a Langmuir model, are ∼190 L kg⁻¹.Calcite single crystals grown in the presence of As(V) show well-developed rhombohedral morphology with characteristic growth hillocks on surfaces at low As(V) concentrations (?5 μM), but habit modification is evident at As(V) concentrations ?30 μM in the form of macrostep development preferentially on the − vicinal surfaces of growth hillocks. Micro-X-ray fluorescence element mapping of surfaces shows preferential incorporation of As in the − vicinal faces relative to + vicinals. EXAFS fit results for both adsorption and coprecipitation samples confirm that As occurs in the 5+ oxidation state in tetrahedral coordination with oxygen, i.e., as arsenate. For adsorption samples, As(V) forms inner-sphere surface complexes via corner-sharing with Ca octahedra. As(V) coprecipitated with calcite substitutes in carbonate sites but with As off-centered, as indicated by two Ca shells, and with likely disruption of local structure. The results indicate that As(V) interacts strongly with the calcite surface, similar to often-cited analog phosphate, and uptake can occur via both adsorption and coprecipitation reactions. Therefore, calcite may be effective for partial removal of dissolved arsenate from aquatic and soil systems.     

Iron isotopes constrain biologic and abiologic processes in banded iron formation genesis

The voluminous 2.5 Ga banded iron formations (BIFs) from the Hamersley Basin (Australia) and Transvaal Craton (South Africa) record an extensive period of Fe redox cycling. The major Fe-bearing minerals in the Hamersley-Transvaal BIFs, magnetite and siderite, did not form in Fe isotope equilibrium, but instead reflect distinct formation pathways. The near-zero average δ⁵⁶Fe values for magnetite record a strong inheritance from Fe³⁺ oxide/hydroxide precursors that formed in the upper water column through complete or near-complete oxidation. Transformation of the Fe³⁺ oxide/hydroxide precursors to magnetite occurred through several diagenetic processes that produced a range of δ⁵⁶Fe values: (1) addition of marine hydrothermal , (2) complete reduction by bacterial dissimilatory iron reduction (DIR), and (3) interaction with excess that had low δ⁵⁶Fe values and was produced by DIR. Most siderite has slightly negative δ⁵⁶Fe values of ∼ −0.5‰ that indicate equilibrium with Late Archean seawater, although some very negative δ⁵⁶Fe values may record DIR. Support for an important role of DIR in siderite formation in BIFs comes from previously published C isotope data on siderite, which may be explained as a mixture of C from bacterial and seawater sources.Several factors likely contributed to the important role that DIR played in BIF formation, including high rates of ferric oxide/hydroxide formation in the upper water column, delivery of organic carbon produced by photosynthesis, and low clastic input. We infer that DIR-driven Fe redox cycling was much more important at this time than in modern marine systems. The low pyrite contents of magnetite- and siderite-facies BIFs suggests that bacterial sulfate reduction was minor, at least in the environments of BIF formation, and the absence of sulfide was important in preserving magnetite and siderite in the BIFs, minerals that are poorly preserved in the modern marine record. The paucity of negative δ⁵⁶Fe values in older (Early Archean) and younger (Early Proterozoic) BIFs suggests that the extensive 2.5 Ga Hamersley-Transvaal BIFs may record a period of maximum expansion of DIR in Earth’s history.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

Kinetic analysis of cation exchange in birnessite using time-resolved synchrotron X-ray diffraction

In this study, we applied time-resolved synchrotron X-ray diffraction (TRXRD) to develop kinetic models that test a proposed two-stage reaction pathway for cation exchange in birnessite. These represent the first rate equations calculated for cation exchange in layered manganates. Our previous work has shown that the substitution of K, Cs, and Ba for interlayer Na in synthetic triclinic birnessite induces measurable changes in unit-cell parameters. New kinetic modeling of this crystallographic data supports our previously postulated two-stage reaction pathway for cation exchange, and we can correlate the kinetic steps with changes in crystal structure. In addition, the initial rates of cation exchange, R (ų min⁻¹), were determined from changes in unit-cell volume to follow these rate laws: R = 1.75[]^0.56, R = 41.1[]^1.10, R = 1.15[]^0.50. Thus, the exchange rates for Na in triclinic birnessite decreased in the order: Cs ? K > Ba. These results are likely a function of hydration energy differences of the cations and the preference of the solution phase for the more readily hydrated cation.     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

An experimental study of magnesite dissolution rates at neutral to alkaline conditions and 150 and 200 °C as a function of pH, total dissolved carbonate concentration, and chemical affinity

Steady-state magnesite dissolution rates were measured in mixed-flow reactors at 150 and 200 °C and 4.6 < pH < 8.4, as a function of ionic strength (0.001 M ? I ? 1 M), total dissolved carbonate concentration (10⁻⁴ M < ΣCO₂ < 0.1 M), and distance from equilibrium. Rates were found to increase with increasing ionic strength, but decrease with increasing temperature from 150 to 200 °C, pH, and aqueous CO₃²⁻ activity. Measured rates were interpreted using the surface complexation model developed by Pokrovsky et al. (1999a) in conjunction with transition state theory (Eyring, 1935). Within this formalism, magnesite dissolution rates are found to be consistent with

     

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.     

Determination of Sr and Ba partition coefficients between apatite from fish (Sparus aurata) and seawater: The influence of temperature

The Sr/Ca and Ba/Ca ratios in inorganic apatite are strongly dependent on the temperature of the aqueous medium during precipitation. If valid in biogenic apatite, these thermometers would offer the advantage of being more resistant to diagenesis than those calibrated on biogenic calcite and aragonite. We have reared seabreams (Sparus aurata) in tanks with controlled conditions during experiments lasting for more than 2 years at 13, 17, 23 and 27 °C, in order to determine the variations in Sr and Ba partitioning relative to Ca (D^Sr and D^Ba, respectively) between seawater and fish apatitic hard tissues (i.e. teeth and bones), as a function of temperature. The sensitivity of the Sr and Ba thermometers (i.e. ∂D^Sr/∂T and ∂D^Ba/∂T, respectively), are similar in bone (/∂T = 0.0036 ± 0.0003 and /∂T = 0.0134 ± 0.0026, respectively) and enamel (/∂T = 0.0037 ± 0.0005 and /∂T = 0.0107 ± 0.0026, respectively). The positive values of ∂D^Sr/∂T and ∂D^Ba/∂T in bone and enamel indicate that D^Sr and D^Ba increase with increasing temperature, a pattern opposite to that observed for inorganic apatite. This distinct thermodependent trace element partitioning between inorganic and organic apatite and water highlights the contradictory effects of the crystal-chemical and biological controls on the partitioning of Ca, Sr and Ba in vertebrate organisms. Taking into account the diet Sr/Ca and Ba/Ca values, it is shown that the bone Ba/Ca signature of fish can be explained by Ca-biopurification and inorganic apatite precipitation, whereas both of these processes fail to predict the bone Sr/Ca values. Therefore, the metabolism of Ca as a function of temperature still needs to be fully understood. However, the biogenic Sr thermometer is used to calculate an average seawater temperature of 30.6 °C using the Sr/Ca compositions of fossil shark teeth at the Cretaceous/Tertiary boundary, and a typical seawater Sr/Ca ratio of 0.02. Finally, while the present work should be completed with data obtained in natural contexts, it is clear that Sr/Ca and Ba/Ca ratios in fossil biogenic apatite already constitute attractive thermometers for marine paleoenvironments.     

Sorption of yttrium and rare earth elements by amorphous ferric hydroxide: Influence of solution complexation with carbonate

The influence of solution complexation on the sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide was investigated at 25 °C over a range of pH (4.0-7.1) and carbonate concentrations . Distribution coefficients, defined as , where [MS_i]_T is the total concentration of sorbed YREE, M_T is the total YREE concentration in solution, and [S_i] is the concentration of amorphous ferric hydroxide, initially increased in magnitude with increasing carbonate concentration, and then decreased. The initial increase of is due to sorption of YREE carbonate complexes , in addition to sorption of free YREE ions (M³⁺). The subsequent decrease of , which is more extensive for the heavy REEs, is due to the increasing intensity of YREE solution complexation by carbonate ions. The competition for YREEs between solution complexation and surface complexation was modeled via the equation:

     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.     

Mechanisms of iron oxide transformations in hydrothermal systems

Coexistence of magnetite and hematite in hydrothermal systems has often been used to constrain the redox potential of fluids, assuming that the redox equilibrium is attained among all minerals and aqueous species. However, as temperature decreases, disequilibrium mineral assemblages may occur due to the slow kinetics of reaction involving the minerals and fluids. In this study, we conducted a series of experiments in which hematite or magnetite was reacted with an acidic solution under H₂-rich hydrothermal conditions (T = 100-250 °C, ) to investigate the kinetics of redox and non-redox transformations between hematite and magnetite, and the mechanisms of iron oxide transformation under hydrothermal conditions. The formation of euhedral crystals of hematite in 150 and 200 °C experiments, in which magnetite was used as the starting material, indicates that non-redox transformation of magnetite to hematite occurred within 24 h. The chemical composition of the experimental solutions was controlled by the non-redox transformation between magnetite and hematite throughout the experiments. While solution compositions were controlled by the non-redox transformation in the first 3 days in a 250 °C experiment, reductive dissolution of magnetite became important after 5 days and affected the solution chemistry. At 100 °C, the presence of maghemite was indicated in the first 7 days. Based on these results, equilibrium constants of non-redox transformation between magnetite and hematite and those of non-redox transformation between magnetite and maghemite were calculated. Our results suggest that the redox transformation of hematite to magnetite occurs in the following steps: (1) reductive dissolution of hematite to and (2) non-redox transformation of hematite and to magnetite.     

Thermodynamic properties of anhydrous smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1 as determined by calorimetric methods. Part I: Heat capacities, heat contents and entropies

The heat capacities of the anhydrous international reference clay minerals, smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1, were measured by low temperature adiabatic calorimetry and differential scanning calorimetry, from 6 to 520 K (at 1 bar). The samples were chemically purified and Na-saturated. Dehydrated clay fractions <2 μm were studied. The structural formulae of the corresponding clay minerals, obtained after subtracting the remaining impurities, are K_0.026Na_0.435Ca_0.010(Si_3.612Al_0.388) (Al_1.593Mg_0.228Ti_0.011)O₁₀(OH)₂ for smectite MX-80, K_0.762Na_0.044(Si_3.387Al_0.613) (Al_1.427Mg_0.241O₁₀(OH)₂ for illite IMt-2 and K_0.530Na_0.135(Si_3.565Al_0.435)(Al_1.709Mg_0.218Ti_0.005)O₁₀(OH)₂for mixed-layer ISCz-1. From the heat capacity values, we determined the molar entropies, standard entropies of formation and heat contents of these minerals. The following values were obtained at 298.15 K and 1 bar:

	(J mol−1 K−1)	S0 (J mol−1 K−1)
Smectite MX-80	326.13 ± 0.10	280.56 ± 0.16
Illite IMt-2	328.21 ± 0.10	295.05 ± 0.17
Mixed-layer ISCz-1	320.79 ± 0.10	281.62 ± 0.15

Full-size table

     

Surface complexation of Cu on birnessite (δ-MnO2): Controls on Cu in the deep ocean

Hexagonal birnessite (δ-MnO₂) is a close analogue to the dominant phase in hydrogenetic marine ferromanganese crusts and nodules. These deposits contain ∼0.25 wt.% Cu which is believed to be scavenged from the overlying water column where Cu concentrations are near 0.1 μg/L. Here, we measured the sorption of Cu on δ-MnO₂ as a function of pH and surface loading. We characterized the nature of the Cu sorption complex at pH 4 and 8 using EXAFS spectroscopy and find that, at pH 4, Cu sorbs to birnessite by inner-sphere complexation on the {0 0 1} surface at sites above Mn vacancies to give a three to fourfold coordinated complex with 6 Mn neighbors at ∼3.4 Å. At pH 8, however, we find that some Cu has become structurally incorporated into the MnO₂ layer by occupying the vacancy sites to give 6 Mn neighbors at ∼2.91 Å. Density functional calculations on and clusters predict a threefold coordinated surface complex and show that the change from surface complexation to structural incorporation is a response to protonation of oxygens surrounding the vacancy site. Consequently, we propose that the transformation between sorption via surface complex and vacancy site occupancy should be reversible. By fitting the Cu sorption as a function of surface loading and pH to the formation of the observed and predicted surface complex, we developed a surface complexation model (in the basic Stern approximation) for the sorption of Cu onto birnessite. Using this model, we demonstrate that the concentration of inorganic Cu in the deep ocean should be several orders of magnitude lower than the observed total dissolved Cu. We propose that the observed total dissolved Cu concentration in the oceans reflects solubilization of Cu by microbially generated ligands.