O/O and D/H ratios of pedogenic kaolinite in a North American Cenomanian laterite: Paleoclimatic implications

Kaolinite, gibbsite and quartz are the dominant minerals in samples collected from two outcrops of a Cenomanian (∼95 Ma) laterite in southwestern Minnesota. A combination of measured yields and isotope ratios permitted mass balance calculations of the δD and δ¹⁸O values of the kaolinite in these samples. These calculations yielded kaolinite δD values of about −73‰ and δ¹⁸O values of about +18.7‰. The δD and δ¹⁸O values appear to preserve information on the ancient weathering system.If formed in hydrogen and oxygen isotope equilibrium with water characterized by the global meteoric water line (GMWL), the kaolinite δD and δ¹⁸O values indicate a crystallization temperature of 22 (±5) °C. A nominal paleotemperature of 22 °C implies a δ¹⁸O value for the corresponding water of −6.3‰. The combination of temperature and meteoric water δ¹⁸O values is consistent with relatively intense rainfall at that mid-paleolatitude location (∼40°N) on the eastern shore of the North American Western Interior Seaway. The inferred Cenomanian paleosol temperature of ∼22 °C is in general accord with published mid-Cretaceous continental mean annual temperatures (MAT) estimated from leaf margin analyses of fossil plants.When compared with results from a published GCM-based Cenomanian climate simulation which specifies a latitudinal sea surface temperature (SST) gradient that was either near modern or smaller-than-modern, the kaolinite paleotemperature of 22 °C is closer to the GCM-predicted MAT for a smaller equator-to-pole temperature difference in the mid-Cretaceous. Moreover, the warm, kaolinite-derived, mid-paleolatitude temperature of 22 °C is associated with proxy estimates of high concentrations of atmospheric CO₂ in the Cenomanian. The overall similarity of proxy and model results suggests that the general features of Cenomanian continental climate in that North American locale are probably being revealed.     

Coexisting goethite and gibbsite from a high-paleolatitude (55°N) Late Paleocene laterite: Concentration and C/C ratios of occluded CO2 and associated organic matter

A Late Paleocene (∼60 Ma BP) lateritic soil from Northern Ireland (the Antrim paleosol, herein referred to as Nire) contains coexisting goethite, gibbsite, phyllosilicate, and hematite. The Fe(III) oxides exhibit pisolitic and Liesegang-type morphologies that are mutually exclusive in hand specimens. X-ray diffraction (XRD) measurements of Al substituted for Fe in goethite indicate two populations: (1) low-Al, Liesegang-type goethites (∼0 mol% Al) and (2) high-Al, pisolitic goethites (∼9 to ∼24 mol% Al). Selective dissolution and incremental vacuum dehydration-decarbonation were used to determine the concentration and δ¹³C values of CO₂ occluded in the respective structures of the goethites and gibbsites in this complex mixture of Nire lateritic minerals. The Fe(CO₃)OH component in the high-Al goethites appears to retain a proxy carbon isotopic record of vadose zone CO₂ in the ancient soil. The δ¹³C values of CO₂ occluded in coexisting goethites and gibbsites indicate that these minerals did not form in equilibrium with the same environmental CO₂.The measured mole fractions (X) of Fe(CO₃)OH in the high-Al goethites range from 0.0059 (±0.0005) to 0.0077 (±0.0006) and correspond to soil CO₂ concentrations of ∼28,000 to ∼37,000 ppmV. The average values of X and δ¹³C for the four high-Al goethites are 0.0067 ± 0.0007 and −20.1 ± 0.5‰, respectively. The δ¹³C value of the organic matter undergoing oxidation in this midlatitude (∼55°N) Late Paleocene soil appears to have been ∼ −28.2‰. Taken together, these data indicate an atmospheric CO₂ concentration of ∼2400 ppmV (± ∼1200 ppmV) at ∼60 Ma BP. The inferred high concentration of atmospheric CO₂ would have been coincident with the warm global climate of the Late Paleocene and is consistent with the idea that CO₂ plays an important role in climate variation.     

Environmental memory and a possible seasonal bias in the stable isotope composition of (U-Th)/He-dated goethite from the Canadian Arctic

Goethite (Ax-2) from Axel Heiberg Island (∼80°N) on the margin of the Arctic Ocean is the dominant mineral in a sample of “petrified” Eocene wood, but U, Th, and He measurements suggest that the goethite (α-FeOOH) crystallized in the latest Miocene/Pliocene (ca. 5.5 to 2.8 Ma). Measured δD and δ¹⁸O values of Ax-2 are −221 (±6)‰ and −9.6 (±0.5)‰, respectively. The inferred δD and δ¹⁸O values of the ancient water were about −139‰ and −18.6‰, respectively, with a calculated temperature of crystallization of 3 (±5)°C, which compares with the modern summer (J-J-A) temperature of 3 °C and contrasts with a modern MAT of −19 °C. Published results from various biological proxies on nearby Ellesmere Island indicate a Pliocene (∼4 Ma) MAT of either −6 or −0.4 °C and corresponding seasonal amplitudes of about 18 or 13 °C. A conductive heat flow model suggests that a temperature of 3 °C could represent goethite crystallization at depths of ∼100-200 cm (for MAT = −6 °C) or ∼250-450 cm (for MAT = −0.4 °C) over seasonally restricted intervals of time.The δ¹⁸O value of the Ax-2 water (−18.6‰) is more positive than the modern J-J-A precipitation (−22‰). In combination, the paleotemperatures and δ¹⁸O values of ancient waters (from Ax-2 and published results from three Eocene or Pliocene proxy sites on Axel Heiberg and Ellesmere Islands) are consistent with a warm season bias in those isotopic proxies. The results are also consistent with higher proportions of J-J-A precipitation in the annual total. If so, this emphasizes the importance of seasonality at high latitudes even in times of warmer global climates, and suggests that the Arctic hydrologic cycle, as expressed in the seasonal distribution and isotopic composition of precipitation (perhaps modified by a warmer Arctic Ocean), differed from modern.The δ¹³C value of the Fe(CO₃)OH component in the Ax-2 goethite is +6.6‰, which is much more positive than expected if crystallizing goethite incorporated CO₂ derived primarily from oxidation of relict Eocene wood with δ¹³C values of about −24‰. This apparent paradox may be resolved if the goethite is a product of oxidation of ¹³C-rich siderite, which had previously replaced wood in an Eocene methanogenic burial environment. Thus, the goethite retains a carbon isotope “memory” of a diagenetic Eocene event, but a δD and δ¹⁸O record of the latest Miocene/Pliocene Arctic climate.     

Iron weathering products in a CO2 + (H2O or H2O2) atmosphere: Implications for weathering processes on the surface of Mars

Various iron-bearing primary phases and rocks have been weathered experimentally to simulate possible present and past weathering processes occurring on Mars. We used magnetite, monoclinic and hexagonal pyrrhotites, and metallic iron as it is suggested that meteoritic input to the martian surface may account for an important source of reduced iron. The phases were weathered in two different atmospheres: one composed of CO₂ + H₂O, to model the present and primary martian atmosphere, and a CO₂ + H₂O + H₂O₂ atmosphere to simulate the effect of strong oxidizing agents. Experiments were conducted at room temperature and a pressure of 0.75 atm. Magnetite is the only stable phase in the experiments and is thus likely to be released on the surface of Mars from primary rocks during weathering processes. Siderite, elemental sulfur, ferrous sulfates and ferric (oxy)hydroxides (goethite and lepidocrocite) are the main products in a water-bearing atmosphere, depending on the substrate. In the peroxide atmosphere, weathering products are dominated by ferric sulfates and goethite. A kinetic model was then developed for iron weathering in a water atmosphere, using the shrinking core model (SCM). This model includes competition between chemical reaction and diffusion of reactants through porous layers of secondary products. The results indicate that for short time scales, the mechanism is dominated by a chemical reaction with second order kinetics (k = 7.75 × 10⁻⁵ g⁻¹/h), whereas for longer time scales, the mechanism is diffusion-controlled (De_A = 2.71 × 10⁻¹⁰ m²/h). The results indicate that a primary CO₂- and H₂O-rich atmosphere should favour sulfur, ferrous phases such as siderite or Fe²⁺-sulfates, associated with ferric (oxy)hydroxides (goethite and lepidocrocite). Further evolution to more oxidizing conditions may have forced these precursors to evolve into ferric sulfates and goethite/hematite.     

An experimental study of the solubility of baddeleyite (ZrO2) in fluoride-bearing solutions at elevated temperature

The solubility of baddeleyite (ZrO₂) and the speciation of zirconium have been investigated in HF-bearing aqueous solutions at temperatures up to 400 °C and pressures up to 700 bar. The data obtained suggest that in HF-bearing solutions zirconium is transported mainly in the form of the hydroxyfluoride species ZrF(OH)₃° and ZrF₂(OH)₂°. Formation constants determined for these species (Zr⁴⁺ + nF⁻ + mOH⁻ = ZrF_n(OH)_m°) range from 43.7 at 100 °C to 46.41 at 400 °C for ZrF(OH)₃°, and from 37.25 at 100 °C to 43.88 at 400 °C for ZrF₂(OH)₂°.Although the solubility of ZrO₂ is retrograde with respect to temperature, the measured concentrations of Zr are orders of magnitude higher than those predicted from theoretical extrapolations based on simple fluoride species (ZrF³⁺-ZrF₆²⁻). Model calculations performed for zircon show that zirconium can be transported by aqueous fluids in concentrations sufficient to account for the concentration of this metal at conditions commonly encountered in fluoride-rich natural hydrothermal systems.     

Pisolitic laterite geochemistry in the golden grove massive sulphide district, Western Australia

The two known massive Cu-Zn sulphide deposits in the steeply dipping strata of the Golden Grove District of Western Australia — the Gossan Hill Deposit and Scuddles Prospect — both have a multi-element secondary geochemical dispersion halo in the pisolitic layer of laterite of probable Tertiary age. The Gossan Hill geochemical anomaly is about 1.5 × 2 km and has a Bi-Sn-Mo-In-Sb-As-Cu-Au association. Mean compositions of pisolitic laterite samples range from 4 to 18 times background, with maxima for Bi, Sn and Mo reaching 90 to 100 times background. Relief of Gossan Hill, which rises some 80 m above the laterite plain, would have favoured dispersion during lateritization. Nevertheless at Scuddles, where the laterite surface has only a gentle slope across strike, an anomaly in pisolitic laterite measures 1 to 1.5 by 1.5 to 2.0 km. The anomalous association, Bi-Sn-Mo-Sb-As, is similar to that at Gossan Hill; the mean contents in pisolitic laterite samples from the Scuddles anomaly range from 2.3 to 8 times background, with maxima for Bi being 12 times background.At Gossan Hill, individual pisolites and the concretionary skins contain gossan fragments and cassiterite, indicating mechanical dispersion during interitization. However, anomalous Cu, Zn and As in goethites of the concretionary skins of pisolites and nodules from the centre of the anomaly indicate that there has also been hydromorphic dispersion.The geochemical characteristics of pisolitic laterite samples have been represented by empirically derived additive indices based on up to eleven chalcophile elements. One index is suggested as a method to increase anomaly size and simplify interpretation; a more specific index based on Bi, Sn and Se is suggested for defining more closely the actual centres of mineralization.For routine exploration of lateritic terrain where the dips of the host sequence are steep, 320-m spacing of samples on a triangular grid appears suitable for locating anomalies and generally defining their source. However, for reconnaissance exploration, the density of sampling could be significantly reduced for example to 1 km spacing.     

Oxygen isotope exchange between H2O and CO2 at elevated CO2 pressures: Implications for monitoring of geological CO2 storage

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on δ¹³C values of CO₂ to trace the migration of injected CO₂ in the subsurface. More recently the use of δ¹⁸O values of both CO₂ and reservoir fluids has been proposed as a method for quantifying in situ CO₂ reservoir saturations due to O isotope exchange between CO₂ and H₂O and subsequent changes in δ¹⁸O_H2O values in the presence of high concentrations of CO₂. To verify that O isotope exchange between CO₂ and H₂O reaches equilibrium within days, and that δ¹⁸O_H2O values indeed change predictably due to the presence of CO₂, a laboratory study was conducted during which the isotope composition of H₂O, CO₂, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50 °C and up to 19 MPa) and varying CO₂ pressures. Conditions typical for the Pembina Cardium CO₂ Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that δ¹⁸O values of CO₂ were on average 36.4 ± 2.2‰ (1σ, n = 15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5‰ for the experimental temperatures of 50 °C. By using ¹⁸O enriched water for the experiments it was demonstrated that changes in the δ¹⁸O values of water were predictably related to the fraction of O in the system sourced from CO₂ in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO₂ is related to the total volumetric saturation of CO₂ and water as a fraction of the total volume of the system, it is concluded that changes in δ¹⁸O values of reservoir fluids can be used to calculate reservoir saturations of CO₂ in CCS settings given that the δ¹⁸O values of CO₂ and water are sufficiently distinct.     

Incremental vacuum dehydration-decarbonation experiments on a natural gibbsite (α-Al(OH3)): CO2 abundance and δC values

Incremental vacuum dehydration-decarbonation experiments were performed at 190°C on chemically “cleaned” aliquots of a gibbsite-dominated, Eocene-age bauxite sample with evolution of CO₂ and H₂O. “Plateau” F (CO₂/H₂O ratios) and δ¹³C values of the CO₂ derived from gibbsite were attained over the dehydration interval, X_v(H₂) = 0.16 to 0.67 (i.e., 16 to 67% breakdown of gibbsite). The plateau value of F for gibbsite was 0.0043 ± 0.0003, while the corresponding δ¹³C value of evolved CO₂ was −16.0‰±0.4‰. Additional experiments on chemically cleaned aliquots included (1) treatment with a solution of 0.3M Na-Citrate + 0.1M Na-Dithionite and (2) an exchange experiment with 0.1 bar of ¹³C-depleted CO₂ (−46‰) at 105°C for 64.5 h. Neither of these additional treatments resulted in a measurable perturbation of plateau values of F or δ¹³C for CO₂ evolved from gibbsite during dehydroxylation. These results support published work on Holocene samples which suggested that CO₂ occluded in gibbsite may preserve information on δ¹³C values of CO₂ in ancient terrestrial systems. The plateau values of F observed in the Eocene gibbsite indicate that it may be possible to experimentally calibrate a relationship between the concentration of CO₂ occluded in gibbsite and CO₂ in the environment at the time of crystallization. Such a calibration would significantly enhance the value of gibbsite as a source of information on ancient oxidized carbon systems.     

Calcite dissolution kinetics in the system CaCO3-H2O-CO2 at high undersaturation

Dissolution rates of limestone covered by a water film open to a CO₂-containing atmosphere are controlled by the chemical composition of the CaCO₃-H₂O-CO₂ solution at the water-mineral interface. This composition is determined by the Ca²⁺-concentration at this boundary, conversion of CO₂ into H⁺ and in the solution, and by diffusional mass transport of the dissolved species from and towards the water-limestone interface. A system of coupled diffusion-reaction equations for Ca²⁺, , and CO₂ is derived. The Ca²⁺ flux rates at the surface of the mineral are defined by the PWP-empirical rate law. These flux rates by the rules of stoichiometry must be equal to the flux rates of CO₂ across the air-water interface. In the solution, CO₂ is converted into H⁺ and . At low water-film thickness this reaction becomes rate limiting. The time dependent diffusion-reaction equations are solved for free drift dissolution by a finite-difference scheme, to obtain the dissolution rate of calcite as a function of the average calcium concentration in the water film. Dissolution rates are obtained for high undersaturation. The results reveal two regimes of linear dissolution kinetics, which can be described by a rate law F = α_i(m_ic_eq − c), where c is the calcium concentration in the water film, c_eq the equilibrium concentration with respect to calcite. For index i = 0, a fast rate law, which here is reported for the first time, is found with α₀ = 3 × 10⁻⁶ m s⁻¹ and m₀ = 0.3. For c > m₀c_eq, a slow rate law is valid with α₁ = 3 × 10⁻⁷ m  s⁻¹ and m₁ = 1, which confirms earlier work. The numbers given above are valid for film thickness of several tenths of a millimetre and at 20 °C. These rates are proven experimentally, using a flat inclined limestone plate covered by a laminar flowing water film injected at an input point with known flow rate Q and calcium concentration. From the concentration measured after flow distance x the dissolution rates are determined. These experiments have been performed at a carbon-dioxide pressure of 0.00035 atm and also of 0.01 atm. The results are in good agreement to the theoretical predictions.     

CO2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.     

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

Determination of carbon fractionation factor between aqueous carbonate and CO2(g) in two-direction isotope equilibration

The experiments were conducted in the open CO₂ system to find out the equilibrium fractionation between the carbonate ion and CO₂(g). The existence of isotopic equilibrium was checked using the two-direction approach by passing the CO₂−N₂ gases with different δ¹³C compositions (− 1.5‰ and − 23‰) through the carbonate solution with δ¹³C = − 4.2‰. The Δ_{CO3T²⁻−CO2(g)} equilibrium fractionation is given as 6.03 ± 0.17‰ at 25 °C. Discussion is provided about the significance of carbonate complexing in determination of Δ_{CO3T²⁻−CO2(g)} and Δ_{HCO3T⁻−CO2(g)} fractionations. Finally, an isotope numerical model of flow and kinetics of hydration and dehydroxylation is built to predict the isotopic behaviour of the system with time.     

Dissolution and carbonation of a serpentinite: Inferences from acid attack and high P-T experiments performed in aqueous solutions at variable salinity

Dissolution experiments on a serpentinite were performed at 70 °C, 0.1 MPa, in H₂SO₄ solution, in open and closed systems, in order to evaluate the overall dissolution rate of mineral components over different times (4, 9 and 24 h). In addition, the serpentinite powder was reacted with a NaCl-bearing aqueous solution and supercritical CO₂ for 24 h at higher pressures (9-30 MPa) and temperatures (250-300 °C) either in a stirred reactor or in an externally-heated pressure vessel to assess both the dissolution rate of serpentinite minerals and the progress of the carbonation reaction. Results show that, at 0.1 MPa, MgO extraction from serpentinite ranges from 82% to 98% and dissolution rate varies from 8.5 × 10⁻¹⁰ mole m⁻² s⁻¹ to 4.2 × 10⁻⁹ mole m⁻² s⁻¹. Attempts to obtain carbonates from the Mg-rich solutions by increasing their pH failed since Mg- and NH₄- bearing sulfates promptly precipitated. On the other hand, at higher pressures, significant crystallization (5.0-10.4 wt%) of Ca- and Fe-bearing magnesite was accomplished at 30 MPa and 300 °C using 100 g L⁻¹ NaCl aqueous solutions. The corresponding amount of CO₂ sequestered by crystallization of carbonates is 9.4-15.9 mole%. Dissolution rate (from 6.3 × 10⁻¹¹ mole m⁻² s⁻¹ to 1.3 × 10⁻¹⁰ mole m⁻² s⁻¹) is lower than that obtained at 0.1 MPa and 70 °C but it is related to pH values much higher (3.3-4.4) than that (−0.65) calculated for the H₂SO₄ solution.Through a thorough review of previous experimental investigations on the dissolution kinetics of serpentine minerals the authors propose adopting: (i) the log rate [mole m⁻² s⁻¹] value of −12.08 ± 0.16 (1σ), as representative of the neutral dissolution mechanism at 25 °C and (ii) the following relationship for the acidic dissolution mechanism at 25 °C:

log rate=-0.45(±0.09)×pH-10.01(±0.30).     

O/O and D/H in goethite from a North American Oxisol of the Early Eocene climatic optimum

An Early Eocene Oxisol in the Ione Formation of California formed in a coastal continental weathering environment at a paleolatitude of ∼38°N. The dominant minerals in the Oxisol are goethite, quartz, and kaolinite. Material balance calculations were applied to new measurements of chemical composition, D/H, and ¹⁸O/¹⁶O ratios of Oxisol samples to determine the δD (−150 ± 3‰) and δ¹⁸O (−2.4 ± 0.3‰) values of the goethite (α-FeOOH). These data, in combination with the global meteoric water line (MWL), yielded an isotopic temperature of 21(±4) °C. The nominal value of 21 °C contrasts with the modern mean annual temperature (MAT) of 16 °C in that area. The warmer temperature is consistent with formation of the goethite during the Early Eocene climatic optimum. The isotopic composition of the goethite and a temperature of 21 °C imply ancient water with a δD value of −61(±4)‰ and a δ¹⁸O value of −8.9(±0.5)‰. This Early Eocene δ¹⁸O (or δD) value is more negative than values in the range of isotopic scatter observed for modern global precipitation at sites with a MAT of 21 °C.At times of warm global climates, the location of a near-surface atmospheric isotherm would generally shift relative to its location under modern climatic conditions. A simple Rayleigh-type condensation model indicates that, if one “follows the isotherm”, the associated scatter in δD and δ¹⁸O of precipitation in very warm global climates should shift (for a given isotherm) to more negative values that may be detectable in proxy records. The isotopic results from the goethite of the Early Eocene Oxisol appear to add to evidence in support of this idea.     

Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field,Taupo Volcanic Zone,New Zealand

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW_e). Soil CO₂ degassing was quantified with 2663 CO₂ flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m⁻²) using published soil temperature heat flow functions. Both CO₂ flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO₂ gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO₂ emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO₂ emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO₂ emissions due to the loss of >99% of the original CO₂ content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO₂ release to the atmosphere from the high temperature reservoir at Ohaaki. The CO₂ flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO₂ flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m⁻² d⁻¹) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO₂ diffuses through porous media of the soil zone. For high-flux sites (>300 g m⁻² d⁻¹), the δ¹³CO₂ signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO₂ source for each respective upflow. Flux thresholds of <30 g m⁻² d⁻¹ for purely diffusive gas transport, between 30 and 300 g m⁻² d⁻¹ for combined diffusive–advective transport, and ?300 g m⁻² d⁻¹ for purely advective gas transport at Ohaaki were assigned. δ¹³CO₂ values and cumulative probability plots of CO₂ flux data both identified a threshold of ∼15 g m⁻² d⁻¹ by which background (atmospheric and soil respired) CO₂ may be differentiated from hydrothermal CO₂.     

A general algorithm for the O abundance correction to C/C determinations from CO2 isotopologue measurements, including CO2 characterised by ‘mass-independent’ oxygen isotope distributions

It is widely recognised that a significant limitation to the ultimate precision of carbon stable isotope ratio measurements, as obtained from dual-inlet mass spectrometric measurements of CO₂ isotopologue ion abundances at m/z 44, 45, and 46, is the correction for interference from ¹⁷O-bearing molecular ions. Two long-established, alternative procedures for determining the magnitude of this correction are in widespread use (although only one has IAEA approval); their differences lead to small but potentially significant discrepancies in the magnitude of the resulting correction. Furthermore, neither approach was designed to accommodate oxygen three-isotope distributions which do not conform to terrestrial mass-dependent behaviour. Stratospheric CO₂, for example, contains a strongly ‘mass-independent’ oxygen isotope composition. A new strategy for determining the ¹⁷O-bearing ion correction is presented, for application where the oxygen three-isotope characteristics of the analyte CO₂ are accurately known (or assigned) in terms of the slope λ of the three-isotope fractionation line and the ordinate axis intercept 10³ ln(1 + k) on a 10³ ln(1 + δ¹⁷O) versus 10³ ln(1 + δ¹⁸O) plot. At the heart of the approach is the relationship between ¹⁷R, which is the ¹⁷O/¹⁶O ratio of the sample CO₂, and other assigned or empirically determined parameters needed for the δ¹³C evaluation:

     

A model for O/O variations in CO2 evolved from goethite during the solid-state α-FeOOH to α-Fe2O3 phase transition

“Plateau” δ¹⁸O values of CO₂ that evolved from the Fe(CO₃)OH component during isothermal vacuum dehydrations (200-230 °C) of 18 natural goethites range from 8.2 to 28.1‰. In contrast, the measured δ¹⁸O values of the goethite structural oxygen range from −11.3 to 1.7‰. The results of this study indicate that the apparent oxygen isotope fractionation factor (¹⁸α_app) between plateau CO₂ and initial goethite is systematically related to the rate of isothermal vacuum dehydration. The nonlinear correlation and the magnitudes of the ¹⁸α_app values are predicted by a relatively simple mass balance model with the following assumptions: (1) the rate of isothermal vacuum dehydration of goethite (for the interval from 0 to ∼60 to 80% loss of structural hydroxyl hydrogen) can be reasonably well represented by first-order kinetics and (2) isotopic exchange between evolving H₂O vapor and solid occurs only in successive, local transition states. The generally good correspondence between the model predictions and the experimental data seems to validate these assumptions. Thus, the ¹⁸O/¹⁶O ratios of the evolved CO₂ can act as probes into the transient processes operating at the molecular level during the solid-state goethite-to-hematite phase transition. For example, the activation energy for the rate constant associated with the transition state, oxygen isotopic exchange between solid and H₂O vapor, is tentatively estimated as 28 ± 11 KJ/mol. Such knowledge may be of consequence in understanding the significance of ¹⁸O/¹⁶O ratios in hematites from some natural environments (e.g., Mars?).Kinetic data and δ¹⁸O values of CO₂ are routinely obtained in the course of measurements of the abundance and δ¹³C values of the Fe(CO₃)OH in goethite. The observed correlation between ¹⁸α_app and dehydration rates suggests that plateau δ¹⁸O values of evolved CO₂ may provide complementary estimates of the δ¹⁸O values of total goethite structural oxygen (O, OH, CO₂) with an overall precision of about ±1‰. However, because of isotopic exchange during the dehydration process, δ¹⁸O values of the evolved CO₂ do not reflect the original δ¹⁸O values of the CO₂ that was occluded as Fe(CO₃)OH in goethite.     

Experimental tests of the effects of Al substitution on the goethite-water D/H fractionation factor

Twelve goethite samples with different degrees of substitution of Al for Fe were synthesized at 22-48 °C and pH values of 1.5-14 under closed system conditions and used to study the effects of Al substitution on the hydrogen isotopic fractionation between goethite and its ambient water. The syntheses followed two pathways: (1) Fe³⁺ hydrolysis in high pH aqueous solutions; (2) oxidation of Fe²⁺ to Fe³⁺ in mid to low pH solutions. XRD and SEM analyses indicated that, irrespective of temperature and pH, goethite was the predominant product of the syntheses in all of the experiments (with degrees of Al substitution as high as ∼13 mol %). “High temperature nonstoichiometric” (HTN) water is present in all of the samples and rapidly exchanges D/H with ambient vapor at room temperature. Uncertainties in the value of the apparent D/H fractionation factor (α_e-v) between HTN water and ambient exchange water at 22 °C lead to significant uncertainties in determinations of the δD values of structural hydrogen (δD_s) in goethites which contain high proportions of HTN water. As determined for the samples of this study, α_e-v has a nominal value of 0.942 (±0.02). δD_s values determined using an α_e-v value of 0.942 indicate that Al substitution increases the δD value of structural hydrogen in goethite by about 1.4 (±0.4)‰ for each increase in Al of 1 mol %. This dependence on Al is of the same sign as, but somewhat larger in magnitude than, the effect of Al predicted by a published model (∼0.7‰ per mol % Al). The overall uncertainties in the current results suggest that an increase of ∼1‰ per mol % Al, as adopted by previous studies, may be a reasonable estimate with which to adjust δ D_s values of natural goethites to those of the pure FeOOH endmember and could be valid for degrees of Al substitution of up to at least 15 mol %. These synthesis experiments also yield a hydrogen isotopic fractionation factor (^Dα_G-W) between pure goethite (α-FeOOH) and liquid water of 0.900 (±0.006), which is analytically indistinguishable from the published value of 0.905 (±0.004). Thus, use of an ^Dα_G-W value of 0.905 in applications to the FeOOH component of natural goethites is supported by the current study.     

Carbon isotopic fractionation during a mesocosm bloom experiment dominated by Emiliania huxleyi: Effects of CO2 concentration and primary production

We investigated the effect of CO₂ and primary production on the carbon isotopic fractionation of alkenones and particulate organic matter (POC) during a natural phytoplankton bloom dominated by the coccolithophore Emiliania huxleyi. In nine semi-closed mesocosms (∼11 m³ each), three different CO₂ partial pressures (pCO₂) in triplicate represented glacial (∼180 ppmv CO₂), present (∼380 ppmv CO₂), and year 2100 (∼710 ppmv CO₂) CO₂ conditions. The largest shift in alkenone isotopic composition (4-5‰) occurred during the exponential growth phase, regardless of the CO₂ concentration in the respective treatment. Despite the difference of ∼500 ppmv, the influence of pCO₂ on isotopic fractionation was marginal (1-2‰). During the stationary phase, E. huxleyi continued to produce alkenones, accumulating cellular concentrations almost four times higher than those of exponentially dividing cells. Our isotope data indicate that, while alkenone production was maintained, the interaction of carbon source and cellular uptake dynamics by E. huxleyi reached a steady state. During stationary phase, we further observed a remarkable increase in the difference between δ¹³C of bulk organic matter and of alkenones spanning 7-12‰. We suggest that this phenomenon is caused mainly by a combination of extracellular release of ¹³C-enriched polysaccharides and subsequent particle aggregation induced by the production of transparent exopolymer particles (TEP).