Seasonal patterns of carbon chemistry and isotopes in tufa depositing groundwaters of southwestern Japan

Annually laminated carbonates, known as tufas, commonly develop in limestone areas and typically record seasonal patterns of oxygen- and carbon-isotope compositions. δ¹⁸O values are principally controlled by seasonal changes of water temperature, whereas δ¹³C values are the result of complex reactions among the gaseous, liquid, and solid sources of carbon in the system. We examined the processes that cause the seasonal patterns of δ¹³C in groundwater systems at three tufa-depositing sites in southwestern Japan by applying model calculations to geochemical data. Underground inorganic carbon species are exchanged with gaseous CO₂, which is mainly introduced to the underground hydrological system by natural atmospheric ventilation and by diffusion of soil air. These processes control the seasonal pattern of δ¹³C, which is low in summer and high in winter. Among the three sites we investigated, we identified two extreme cases of the degree of carbon exchange between liquid and gaseous phases. For the case with high radiocarbon composition (Δ¹⁴C) and low pCO₂, there was substantial carbon exchange because of a large contribution of atmospheric CO₂ and a small water mass. For the other extreme case, which was characterized by low Δ¹⁴C and high pCO₂, the contribution of atmospheric CO₂ was small and the water mass was relatively large. Our results suggest that at two of the three sites water residence time within the soil profile was longer than 1 year. Our results also suggested a short residence time (less than 1 year) of water in the soil profile at the site with the smallest water mass, which is consistent with large seasonal amplitude of the springwater temperature variations. The Δ¹⁴C value of tufas is closely related to the hydrological conditions in which they are deposited. If the initial Δ¹⁴C value of a tufa-depositing system is stable, ¹⁴C-chronology can be used to date paleo-tufas.     

Experimental constraints on the parental liquid of the Chassigny meteorite: A possible link between the Chassigny meteorite and a Martian Gusev basalt

Most magmas proposed as parental to the Martian SNC meteorites are high in iron and low in alumina. Yet, experiments at low pressures on such liquids have not produced the cumulate or melt-inclusion assemblages seen in the chassignite meteorites. Therefore, elevated pressure experiments under anhydrous and hydrous (water-undersaturated) conditions were conducted on a high-Fe, low-Al liquid proposed to be parental to the Chassigny meteorite. These experiments failed to produce the most magnesian cumulate phases, as well as the olivine hosted kaersutite-bearing melt-inclusion assemblage, of the chassignites. These results suggest that the parental liquid to the chassignite meteorites is both more magnesium and aluminum-rich than the previously considered composition (A^∗; Johnson et al., 1991). The proposed composition is similar to the Martian Adirondack class Gusev basalt Humphrey and suggests a link between the Chassigny meteorite and rocks on the surface of Mars.     

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.     

Cr-spinel/olivine and Cr-spinel/liquid nickel partition coefficients from natural samples

The Nernst partition coefficient of nickel (D^Ni) between Cr-spinel and silicate melt in natural systems has been investigated using mid-ocean ridge basalts (MORB) and other volcanic rocks. The Cr-spinel/olivine D^Ni values in volcanic rocks are between 1.2 and 0.3, indicating that the Cr-spinel/liquid D^Ni values vary from slightly higher to significantly lower than the olivine/liquid D^Ni values in natural systems. The Cr-spinel/liquid D^Ni values from the MORB samples vary between 6 and 11, slightly higher than those from the S-bearing experiments of Satari et al. [Satari P., Brenan J. M., Horn I. and McDonough W. F. (2002) Experimental constraints on the sulfide- and chromite-silicate melt partitioning behavior of rhenium and platinum-group elements. Economic Geology97, 385-398]. The results of the MORB samples and the experiments of Satari et al. (2002) indicate a negative correlation between the Cr-spinel/liquid D^Ni and the X_Cr values in Cr-spinels (Cr cation number on the basis of 3 total cations in the spinel structure). Variations of Cr-spinel/liquid D^Ni values with Cr-spinel compositions can be estimated from an empirical equation based on the results of the MORB samples and the experiments by Satari et al. (2002). The choice of Cr-spinel/liquid D^Ni = 10 for numerical modeling by Righter et al. [Righter K., Leeman W. P. and Hervig R. L. (2006) Partitioning of Ni, Co, and V between spinel-structured oxides and silicate melts: importance of spinel composition. Chemical Geology227, 1-25] is reasonable for basaltic systems. For picritic and komatiitic systems a lower value of ∼5 is more appropriate.     

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.     

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.     

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.