An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

Kinetics of Cd sorption, desorption and fixation by calcite: A long-term radiotracer study

Time-dependent sorption and desorption of Cd on calcite was studied over 210 days utilizing ¹⁰⁹Cd as a tracer to distinguish between ‘labile’ and ‘non-labile’ forms of sorbed Cd. Stabilizing the calcite suspensions for 12 months under atmospheric P_CO2 and controlled temperature was necessary to reliably follow Cd dynamics following initial sorption. Results revealed time-dependant Cd sorption and marked desorption hysteresis by calcite under environmentally relevant conditions. Data obtained were fitted to a first-order kinetic model and a concentric shell diffusion model. Both models described the progressive transfer of Cd²⁺ to a less reactive form within calcite and subsequent desorption of Cd subject to different initial contact times. The kinetic model provided a better fit to the combined sorption and desorption data (R² = 0.992). It differentiates between two ‘pools’ of sorbed Cd²⁺ on calcite, ‘labile’ and ‘non-labile’, in which labile sorbed Cd is in immediate equilibrium with the free Cd²⁺ ion activity in solution whereas non-labile Cd is kinetically restricted. For the diffusion model (R² = 0.959), the rate constants describing Cd dynamics in calcite produced a half-life for Cd desorption of ∼175 d, for release to a ‘zero-sink’ solution. Results from this study allow comment on the likely mechanisms occurring at the calcite surface following long-term Cd sorption.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

The age of the martian meteorite Northwest Africa 1195 and the differentiation history of the shergottites

Samarium-neodymium isotopic analyses of unleached and acid-leached mineral fractions from the recently identified olivine-bearing shergottite Northwest Africa 1195 yield a crystallization age of 347 ± 13 Ma and an value of +40.1 ± 0.9. Maskelynite fractions do not lie on the Sm-Nd isochron and appear to contain a martian surface component with low ¹⁴⁷Sm/¹⁴⁴Nd and ¹⁴³Nd/¹⁴⁴Nd ratios that was added during shock. The Rb-Sr system is disturbed and does not yield an isochron. Terrestrial Sr appears to have affected all of the mineral fractions, although a maximum initial ⁸⁷Sr/⁸⁶Sr ratio of 0.7016 is estimated by passing a 347 Ma reference line through the maskelynite fraction that is least affected by contamination. The high initial value and the low initial ⁸⁷Sr/⁸⁶Sr ratio, combined with the geologically young crystallization age, indicate that Northwest Africa 1195 is derived from a source region characterized by a long-term incompatible-element depletion.The age and initial Sr and Nd isotopic compositions of Northwest Africa 1195 are very similar to those of Queen Alexandra Range 94201, indicating these samples were derived from source regions with similar Sr-Nd isotopic systematics. These similarities suggest that these two meteorites share a close petrogenetic relationship and might have been erupted from a common volcano. The meteorites Yamato 980459, Dar al Gani 476, Sayh al Uhaymir 005/008, and Dhofar 019 also have relatively old ages between 474 and 575 Ma and trace element and/or isotopic systematics that are indicative of derivation from incompatible-element-depleted sources. This suggests that the oldest group of meteorites is more closely related to one another than they are to the younger meteorites that are derived from less incompatible-element-depleted sources. Closed-system fractional crystallization of this suite of meteorites is modeled with the MELTS algorithm using the bulk composition of Yamato 980459 as a parent. These models reproduce many of the major element and mineralogical variations observed in the suite. In addition, the rare earth element systematics of these meteorites are reproduced by fractional crystallization using the proportions of phases and extents of crystallization that are calculated by MELTS. Other shergottites that demonstrate enrichments in incompatible-elements and have evolved Sr and Nd isotopic systematics have some geochemical systematics that are similar to those observed in the depleted group. Most notably, although they exhibit a very limited range of incompatible trace element and isotopic compositions, they have highly variable major element compositions. This is also consistent with evolution from a common mantle source region by variable amounts of fractional crystallization. If this scenario is correct, it suggests that the combined effects of source composition and fractional crystallization are likely to account for the major element, trace element, and isotopic diversity of all shergottites.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

Association of dissolved aluminum with silica: Connecting molecular structure to surface reactivity using NMR

We studied uptake mechanisms for dissolved Al on amorphous silica by combining bulk-solution chemistry experiments with solid-state Nuclear Magnetic Resonance techniques (²⁷Al magic-angle spinning (MAS) NMR, ²⁷Al{¹H} cross-polarization (CP) MAS NMR and ²⁹Si{¹H} CP-MAS NMR). We find that reaction of Al (1 mM) with amorphous silica consists of at least three reaction pathways; (1) adsorption of Al to surface silanol sites, (2) surface-enhanced precipitation of an aluminum hydroxide, and (3) bulk precipitation of an aluminosilicate phase. From the NMR speciation and water chemistry data, we calculate that 0.20 (±0.04) tetrahedral Al atoms nm⁻² sorb to the silica surface. Once the surface has sorbed roughly half of the total dissolved Al (∼8% site coverage), aluminum hydroxides and aluminosilicates precipitate from solution. These precipitation reactions are dependent upon solution pH and total dissolved silica concentration. We find that the Si:Al stoichiometry of the aluminosilicate precipitate is roughly 1:1 and suggest a chemical formula of NaAlSiO₄ in which Na⁺ acts as the charge compensating cation. For the adsorption of Al, we propose a surface-controlled reaction mechanism where Al sorbs as an inner-sphere coordination complex at the silica surface. Analogous to the hydrolysis of , we suggest that rapid deprotonation by surface hydroxyls followed by dehydration of ligated waters results in four-coordinate (>SiOH)₂Al(OH)₂ sites at the surface of amorphous silica.     

Low-oxygen and chemical kinetic constraints on the geochemical niche of neutrophilic iron(II) oxidizing microorganisms

Neutrophilic iron oxidizing bacteria (FeOB) must actively compete with rapid abiotic processes governing Fe(II) oxidation and as a result have adapted to primarily inhabit low-O₂ environments where they can more successfully compete with abiotic Fe(II) oxidation. The spatial distribution of these microorganisms can be observed through the chemical gradients they affect, as measured using in situ voltammetric analysis for dissolved Fe(II), Fe(III), O₂, and FeS_(aq). Field and laboratory determination of the chemical environments inhabited by the FeOB were coupled with detailed kinetic competition studies for abiotic and biotic oxidation processes using a pure culture of FeOB to quantify the geochemical niche these organisms inhabit. In gradient culture tubes, the maximum oxygen levels, which were associated with growth bands of Sideroxydans lithotrophicus (ES-1, a novel FeOB), were 15-50 μM. Kinetic measurements made on S. lithotrophicus compared biotic/abiotic (killed control) Fe oxidation rates. The biotic rate can be a significant and measurable fraction of the total Fe oxidation rate below O₂ concentrations of approximately 50 μM, but biotic Fe(II) oxidation (via the biotic/abiotic rate comparison) becomes difficult to detect at higher O₂ levels. These results are further supported by observations of conditions supporting FeOB communities in field settings. Variablity in cell densities and cellular activity as well as variations in hydrous ferrous oxide mineral quantities significantly affect the laboratory kinetic rates. The microbial habitat (or geochemical niche) where FeOB are active is thus largely controlled by the competition between abiotic and biotic kinetics, which are dependent on Fe(II) concentration, P_O2, temperature and pH in addition to the surface area of hydrous ferric oxide minerals and the cell density/activity of FeOB. Additional field and lab culture observations suggest a potentially important role for the iron-sulfide aqueous molecular cluster, FeS_(aq), in the overall cycling of iron associated with the environments these microorganisms inhabit.     

Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): II. Reaction-transport modeling

A steady-state reaction-transport model is applied to sediments retrieved by gravity core from two stations (S10 and S13) in the Skagerrak to determine the main kinetic and thermodynamic controls on anaerobic oxidation of methane (AOM). The model considers an extended biomass-implicit reaction network for organic carbon degradation, which includes extracellular hydrolysis of macromolecular organic matter, fermentation, sulfate reduction, methanogenesis, AOM, acetogenesis and acetotrophy. Catabolic reaction rates are determined using a modified Monod rate expression that explicitly accounts for limitation by the in situ catabolic energy yields. The fraction of total sulfate reduction due to AOM in the sulfate-methane transition zone (SMTZ) at each site is calculated. The model provides an explanation for the methane tailing phenomenon which is observed here and in other marine sediments, whereby methane diffuses up from the SMTZ to the top of the core without being consumed. The tailing is due to bioenergetic limitation of AOM in the sulfate reduction zone, because the methane concentration is too low to engender favorable thermodynamic drive. AOM is also bioenergetically inhibited below the SMTZ at both sites because of high hydrogen concentrations (∼3-6 nM). The model results imply there is no straightforward relationship between pore water concentrations and the minimum catabolic energy needed to support life because of the highly coupled nature of the reaction network. Best model fits are obtained with a minimum energy for AOM of ∼11 kJ mol⁻¹, which is within the range reported in the literature for anaerobic processes.     

Seasonal δS variations in two high elevation snow pits measured by S-S double spike thermal ionization mass spectrometry

δ³⁴S and sulfate concentrations were determined in snow pit samples using a thermal ionization mass spectrometric technique capable of 0.2‰ accuracy and requires ≈5 μg (0.16 μmol) natural S. The technique utilizes a ³³S-³⁶S double spike for instrumental mass fractionation correction, and has been applied to snow pit samples collected from the Inilchek Glacier, Kyrgyzstan and from Summit, Greenland. These δ³⁴S determinations provide the first high-resolution seasonal data for these sites, and are used to estimate seasonal sulfate sources. Deuterium (δD) and oxygen (δ¹⁸O) isotope data show that the Inilchek and Summit snow pit samples represent precipitation over ≈20 months.The δ³⁴S values for the Inilchek ranged from +2.6 ± 0.4‰ to +7.6 ± 0.4‰ on sample sizes ranging from 0.3 to 1.8 μmol S. δ³⁴S values for Greenland ranged from +3.6 ± 0.7‰ to +13.3 ± 5‰ for sample sizes ranging from 0.05 to 0.29 μmol S. The concentration ranged from 92.6 ± 0.4 to 1049 ± 4 ng/g for the Inilchek and 18 ± 9 to 93 ± 6 ng/g for the Greenland snow pit. Anthropogenic sulfate dominates throughout the sampled time interval for both sites based on mass balance considerations. Additionally, both sites exhibit a seasonal signature in both δ³⁴S and concentration. The thermal ionization mass spectrometric technique has three advantages compared to gas source isotopic methods: (1) sample size requirements of this technique are 10-fold less permitting access to the higher resolution S isotope record of low concentration snow and ice, (2) the double spike technique permits δ³⁴S and S concentration to be determined simultaneously, and (3) the double spike is an internal standard.