Large fractionation of calcium isotopes during cave-analogue calcium carbonate growth

We have measured δ^44/42Ca of laboratory-precipitated calcite grown in an experimental setup that closely replicates stalagmite formation. Calcium solutions were dripped onto two different substrates in tightly-controlled conditions and calcite precipitated due to rapid CO₂ degassing. With seeded glass slides as the substrate, we observe a Ca isotope ratio in the calcite which is ∼0.5‰ per amu lower than that in the growth solution. This fractionation is generally almost twice that observed in previously published calcite growth experiments and indicates a large kinetic effect on Ca isotopes in the stalagmite growth environment. The precipitate forming near the spot where the drip lands shows slightly greater solution-to-precipitate fractionation than calcite further from the drip reflecting a decrease in this kinetic fractionation as precipitation continues. We interpret these results in the context of the model of Fantle and DePaolo (2007) which involves surface entrapment of light Ca isotopes to decrease calcite δ^44/42Ca, and depletion of Ca from the solution in the direct vicinity of the growing calcite to increase calcite δ^44/42Ca. In the stalagmite setting, the second of these effects is minimized so that calcite Ca isotope ratios are unusually light. This interpretation suggests that stalagmite Ca isotope ratios should decrease with the saturation state of the drip water (i.e. with the growth rate of calcite). Ca isotopes might therefore allow reconstruction of surface entrapment of trace metals and isotopes more generally and might, for instance, allow an assessment of the appropriate relationship between oxygen isotope fractionation and temperature for periods of past growth in stalagmites.     

Kinetic isotopic fractionation during diffusion of ionic species in water

Experiments specifically designed to measure the ratio of the diffusivities of ions dissolved in water were used to determine . The measured ratio of the diffusion coefficients for Li and K in water (D_Li/D_K = 0.6) is in good agreement with published data, providing evidence that the experimental design being used resolves the relative mobility of ions with adequate precision to also be used for determining the fractionation of isotopes by diffusion in water. In the case of Li, we found measurable isotopic fractionation associated with the diffusion of dissolved LiCl (D_⁷Li/D_⁶Li=0.99772±0.00026). This difference in the diffusion coefficient of ⁷Li compared to ⁶Li is significantly less than that reported in an earlier study, a difference we attribute to the fact that in the earlier study Li diffused through a membrane separating the water reservoirs. Our experiments involving Mg diffusing in water found no measurable isotopic fractionation (D_²⁵Mg/D_²⁴Mg=1.00003±0.00006). Cl isotopes were fractionated during diffusion in water (D_³⁷Cl/D_³⁵Cl=0.99857±0.00080) whether or not the co-diffuser (Li or Mg) was isotopically fractionated. The isotopic fractionation associated with the diffusion of ions in water is much smaller than values we found previously for the isotopic fractionation of Li and Ca isotopes by diffusion in molten silicate liquids. A major distinction between water and silicate liquids is that water surrounds dissolved ions with hydration shells, which very likely play an important but still poorly understood role in limiting the isotopic fractionation associated with diffusion.     

Li/Ca in multiple species of benthic and planktonic foraminifera: thermocline, latitudinal, and glacial-interglacial variation

Li/Ca ratios were measured in planktonic and benthic foraminifera from a variety of hydrographic settings to investigate the factors influencing lithium incorporation into foraminiferal tests including temperature, dissolution, pressure, and interspecies differences. Down-core measurements of planktonic (Orbulina universa, Globigerinoides ruber, and Globigerinoides sacculifer) and benthic foraminifera (calcitic Cibicides wuellerstorfi and aragonitic Hoeglandina elegans) show a systematic variation in Li/Ca with δ¹⁸O through the last glacial-interglacial transition. All species examined exhibit an increase in Li/Ca between 14 to 50% from the Holocene to the last glacial maximum. Li/Ca generally increases with decreasing temperature as seen in a latitudinal transect of planktonic O. universa and down-slope benthic species along the Bahama Bank margins. Postdepositional dissolution possibly causes a decrease in planktonic foraminiferal Li/Ca along the Sierra Leone Rise, and increased water depth causes a decrease in benthic foraminiferal Li/Ca in the deep Caribbean. However, none of these effects are sufficient to account for the observed glacial-interglacial changes. Physiological factors such as calcification rate may affect the Li/Ca content of foraminiferal calcite. The calcification rate in turn may be a function of carbonate ion concentration of ambient ocean water. This work shows that incorporation of lithium by foraminifera appears to be influenced by factors other than seawater composition and does not appear to be dominated by changes in temperature, dissolution, or pressure. We hypothesize that the consistent increase in foraminiferal Li/Ca during the last glacial maximum may be linked to changes in seawater carbonate ion concentration. Important parameters to be tested include calcification rate and foraminiferal test size and weight. If foraminiferal Li/Ca is dominantly controlled by calcification rate as a function of seawater carbonate ion concentration, then Li/Ca may act as a proxy of past atmospheric CO₂.     

Lithium isotopes in global mid-ocean ridge basalts

The lithium isotope compositions of 30 well-characterized samples of glassy lavas from the three major mid-ocean ridge segments of the world, spanning a wide range in radiogenic isotope and elemental content and sea floor physical parameters, have been measured. The overall data set shows a significant range in δ⁷Li (+1.6 to +5.6), with no global correlation between Li isotopes and other geochemical or tectonic parameters. The samples with the greatest lithophile element depletion (N-MORB: K₂O/TiO₂ < 0.09) display an isotopic range consistent with the extant database. Samples with greater trace element enrichment display a greater degree of isotopic variability and trend toward heavier compositions (δ⁷Li = +2.4 to +5.6), but are not distinct on average from N-MORB. Together with published data, N-MORB is estimated to have mean δ⁷Li = +3.4 ± 1.4‰ (2σ), consistent with the estimate for an uncontaminated MORB source based on pristine peridotite xenoliths. Locally, where sampling density permits, sources of Li isotope heterogeneity may be evaluated. Sample sets from the East Pacific Rise show correlation of δ⁷Li with halogen concentration ratios. This is interpreted at 15.5°N latitude to represent incorporation of <5 weight percent recycled subduction-modified mantle in the MORB source. At 9.5°N latitude the data are more consistent with shallow level magma chamber contamination by seawater-derived components (<0.5 wt.%).     

Experimental determination of the role of diffusion on Li isotope fractionation during basaltic glass weathering

In order to use lithium isotopes as tracers of silicate weathering, it is of primary importance to determine the processes responsible for Li isotope fractionation and to constrain the isotope fractionation factors caused by each process as a function of environmental parameters (e.g. temperature, pH). The aim of this study is to assess Li isotope fractionation during the dissolution of basalt and particularly during leaching of Li into solution by diffusion or ion exchange. To this end, we performed dissolution experiments on a Li-enriched synthetic basaltic glass at low ratios of mineral surface area/volume of solution (S/V), over short timescales, at various temperatures (50 and 90 °C) and pH (3, 7, and 10). Analyses of the Li isotope composition of the resulting solutions show that the leachates are enriched in ⁶Li (δ⁷Li = +4.9 to +10.5‰) compared to the fresh basaltic glass (δ⁷Li = +10.3 ± 0.4‰). The δ⁷Li value of the leachate is lower during the early stages of the leaching process, increasing to values close to the fresh basaltic glass as leaching progresses. These low δ⁷Li values can be explained in terms of diffusion-driven isotope fractionation. In order to quantify the fractionation caused by diffusion, we have developed a model that couples Li diffusion with dissolution of the glassy silicate network. This model calculates the ratio of the diffusion coefficients of both isotopes (a = D₇/D₆), as well as its dependence on temperature, pH, and S/V. a is mainly dependent on temperature, which can be explained by a small difference in activation energy (0.10 ± 0.02 kJ/mol) between ⁶Li⁺ and ⁷Li⁺. This temperature dependence reveals that Li isotope fractionation during diffusion is low at low temperatures (T < 20 °C), but can be significant at high temperatures. However, concerning hydrothermal fluids (T > 120 °C), the dissolution rate of basaltic glass is also high and masks the effects of diffusion. These results indicate that the high δ⁷Li values of river waters, in particular in basaltic catchments, and the fractionated values of hydrothermal fluids are mainly controlled by precipitation of secondary phases.     

B and Li isotopes as intrinsic tracers for injection tests in aquifer storage and recovery systems

Boron and Li isotopes have been tested as environmental tracers of treated sewage injected into a sandy aquifer (Shafdan reclamation project, Israel). During a 38 days injection test in a newly dug injection well, a conservative artificial tracer (Br⁻) was monitored together with δ¹¹B and δ⁷Li in the injectate, in the unsaturated soil zone (porous cup) and an observation well in the aquifer. In spite of B and Li concentrations in the injectate close to background values, significant shifts of the isotope signatures could be observed over the duration of the injection test. Boron isotope ratios show a breakthrough curve delayed with respect to Br⁻ breakthrough due to some reversible sorption on the aquifer material. No isotope fractionation was observed in the unsaturated or the saturated zone so that B isotopes can be considered as conservative in the investigated part of the aquifer system. Lithium isotopes are strongly fractionated, probably due to sorption processes. Lithium concentrations point to a Li sink in the system, δ⁷Li values vary strongly with a tendency of ⁷Li depletion in the liquid phase over the duration of the experiment. This is opposite to the expected preferential sorption of ⁶Li onto clay minerals. Boron isotopes reveals a valuable tracer of artificial recharge of freshwaters derived from treated sewage, both for short term tracer tests and for long-term monitoring of artificial recharge, even if in aquifers with higher clay contents, sorption-linked isotope fractionation cannot be excluded. More data are needed on Li isotope fractionation in natural groundwater systems to assess the potential of this tracer as monitoring tool.     

Effects of seawater carbonate ion concentration and temperature on shell U, Mg, and Sr in cultured planktonic foraminifera

We investigate the sensitivity of U/Ca, Mg/Ca, and Sr/Ca to changes in seawater [CO₃²⁻] and temperature in calcite produced by the two planktonic foraminifera species, Orbulina universa and Globigerina bulloides, in laboratory culture experiments. Our results demonstrate that at constant temperature, U/Ca in O. universa decreases by 25 ± 7% per 100 μmol [CO₃²⁻] kg⁻¹, as seawater [CO₃²⁻] increases from 110 to 470 μmol kg⁻¹. Results from G. bulloides suggest a similar relationship, but U/Ca is consistently offset by ∼+40% at the same environmental [CO₃²⁻]. In O. universa, U/Ca is insensitive to temperature between 15°C and 25°C. Applying the O. universa relationship to three U/Ca records from a related species, Globigerinoides sacculifer, we estimate that Caribbean and tropical Atlantic [CO₃²⁻] was 110 ± 70 μmol kg⁻¹ and 80 ± 40 μmol kg⁻¹ higher, respectively, during the last glacial period relative to the Holocene. This result is consistent with estimates of the glacial-interglacial change in surface water [CO₃²⁻] based on both modeling and on boron isotope pH estimates. In settings where the addition of U by diagenetic processes is not a factor, down-core records of foraminiferal U/Ca have potential to provide information about changes in the ocean’s carbonate concentration.Below ambient pH (pH < 8.2), Mg/Ca decreased by 7 ± 5% (O. universa) to 16 ± 6% (G. bulloides) per 0.1 unit increase in pH. Above ambient pH, the change in Mg/Ca was not significant for either species. This result suggests that Mg/Ca-based paleotemperature estimates for the Quaternary, during which surface-ocean pH has been at or above modern levels, have not been biased by variations in surface-water pH. Sr/Ca increased linearly by 1.6 ± 0.4% per 0.1 unit increase in pH. Shell Mg/Ca increased exponentially with temperature in O. universa, where Mg/Ca = 0.85 exp (0.096*T), whereas the change in Sr/Ca with temperature was within the reproducibility of replicate measurements.     

Seasonal water-column dynamics of dissolved inorganic carbon stable isotopic compositions (δCDIC) in small hardwater lakes in Minnesota and Montana

The carbon stable isotopic value of dissolved inorganic carbon (δ¹³C_DIC) was measured over several years at different depths in the water column in six carbonate-precipitating temperate lakes. δ¹³C_DIC behavior in three of these lakes departed from the conventional model wherein epilimnetic waters are seasonally enriched relative to all hypolimnetic waters, and in general δ¹³C_DIC values in the water column were not readily correlated to parameters such as lake stratification, algal productivity, hydraulic residence time, or water chemistry. Additionally, the processes implicated in generating the δ¹³C_DIC values of individual lakes differ between lakes with similar δ¹³C_DIC compositions. Each lake thus initially appears idiosyncratic, but when the effects of carbonate mineral equilibria, microbial activity, and lake residence time are viewed in terms of the magnitude of distinct DIC pools and fluxes in stratified lakes, generalizations can be made that allow lakes to be grouped by δ¹³C_DIC behavior. We recognize three modes in the relationship between δ¹³C_DIC values and DIC concentration ([DIC]) of individual lakes: (A) δ¹³C_DIC values decreasing with increasing [DIC]; (B) δ¹³C_DIC values increasing with increasing [DIC]; (C) δ¹³C_DIC values decreasing with increasing [DIC] but increasing again at the highest [DIC]. This approach is useful both in understanding δ¹³C_DIC dynamics in modern hardwater lakes and in reconstructing the environmental changes recorded by sedimentary δ¹³C components in the lacustrine paleorecord.     

A large epeiric methanogenic Bambuí sea in the core of Gondwana supercontinent?

Carbon isotope compositions of both sedimentary carbonate and organic matter can be used as key proxies of the global carbon cycle and of its evolution through time,as long as they are acquired from waters where the dissolved inorganic carbon(DIC)is in isotope equilibrium with the atmospheric CO2.However,in shallow water platforms and epeiric settings,the influence of local to regional parameters on carbon cycling may lead to DIG isotope variations unrelated to the global carbon cycle.This may be especially true for the terminal Neoproterozoic,when Gondwana assembly isolated waters masses from the global ocean,and extreme positive and negative carbon isotope excursions are recorded,potentially decoupled from global signals.To improve our understanding on the type of information recorded by these excursions,we investigate the pairedδ^13Ccarb andδ^13Corg evolution for an increasingly restricted late Ediacaran-Cambrian foreland system in the West Gondwana interior:the basal Bambui Group.This succession represents a 1~(st)-order sedimentary sequence and records two majorδ^13Ccarb excursions in its two lowermost lower-rank sequences.The basal cap carbonate interval at the base of the first sequence,deposited when the basin was connected to the ocean,hosts antithetical negative and positive excursions forδ^13Ccarb andδ^13Corg,respectively,resulting inΔ^13C values lower than 25‰.From the top of the basal sequence upwards,an extremely positiveδ^13Ccarb excursion is coupled toδ^13Corg,reaching values of+14‰and-14‰,respectively.This positive excursion represents a remarkable basin-wide carbon isotope feature of the Bambui Group that occurs with only minor changes inΔ^13C values,suggesting change in the DIC isotope composition.We argue that this regional isotopic excursion is related to a disconnection between the intrabasinal and the global carbon cycles.This extreme carbon isotope excursion may have been a product of a disequilibria between the basin DIC and atmospheric CO2 induced by an active methanogenesis,favored by the basin restriction.The drawdown of sulfate reservoir by microbial sulfate reduction in a poorly ventilated and dominantly anoxic basin would have triggered methanogenesis and ultimately methane escape to the atmosphere,resulting in a^13C-enriched DIC influenced by methanogenic CO2.Isolated basins in the interior of the Gondwana supercontinent may have represented a significant source of methane inputs to the atmosphere,potentially affecting both the global carbon cycle and the climate.     

Stable isotope profiles (Ca,O, C) through modern brachiopod shells of T. septentrionalis and G. vitreus: Implications for calcium isotope paleo-ocean chemistry

Brachiopod shells are widely used as an archive to reconstruct elemental and isotopic composition of seawater. Studies, focused on oxygen and carbon isotopes over the last decades, are increasingly extending to the emerging calcium isotope system. To date, only little attention has been paid to test the reliability of fossil brachiopods on their modern counterparts.In this context, the present study investigates two modern brachiopods, Terebratulina septentrionalis (eastern Canada, 5–30 m depth, 7.1 °C seasonal temperature variation, two-layer shell) and Gryphus vitreus (northern Mediterranean, 200 m depth, constant all-year round temperature, three-layer shell). Both species were sampled along the ontogenetic growth direction and calcium, oxygen, and carbon isotopes as well as elemental concentration were measured. Calcium isotopes were analyzed on TIMS. The elemental composition was analyzed by LA-ICP-MS and ICP-AES.The results indicate an intra-specimen δ^44/40Ca variation ranging from 0.16 to 0.33‰, pointing to a fairly homogenous distribution of calcium isotopes in brachiopod shells. However, in the light of the suggested 0.7‰ increase in calcium isotopes over the Phanerozoic such intra-specimen variations constrain ocean reconstruction. δ^44/40Ca values of T. septentrionalis do not seem to be affected by growth rate. Calcium isotopic values of G. vitreus are heavy in the central part of the shell and trend towards lighter values in peripheral areas approaching the maximum isotopic composition of T. septentrionalis. The maximum inter-species δ^44/40Ca difference of 0.62‰ between T. septentrionalis and G. vitreus indicates that care should be taken when using different taxa, species with different strontium content or brachiopods with specialized shell structure, such as G. vitreus, for ocean water reconstruction in terms of Ca isotopic composition. T. septentrionalis may record Ca isotopic fractionation related to seasonal seawater temperature variations in its shell but this is difficult to resolve at the current analytical precision. Average δ¹⁸O-derived temperatures of the two investigated species are close to on-site measured temperatures.     

Interaction between shallow and deep aquifers in the Tivoli Plain (Central Italy) enhanced by groundwater extraction: A multi-isotope approach and geochemical modeling

In the Tivoli Plain (Rome, Central Italy) the interaction between shallow and deep groundwater flow systems enhanced by groundwater extraction has been investigated using isotopic and chemical tracers. A conceptual model of the groundwater flowpaths has been developed and verified by geochemical modeling. A combined hydrogeochemical and isotopic investigation using ion relationships such as DIC/Cl⁻, Ca/(Ca + Mg)/SO₄/(SO₄ + HCO₃), and environmental isotopes (δ¹⁸O, δ²H, ⁸⁷Sr/⁸⁶Sr, δ³⁴S and δ¹³C) was carried out in order to determine the sources of recharge of the aquifer, the origin of solutes and the mixing processes in groundwater of Tivoli Plain. Multivariate statistical methods such as principal component analysis and Cluster analyses have confirmed the existence of different geochemical facies and the role of mixing in the chemical composition of the groundwater.     

Behavior of stable isotopes of dissolved oxygen,dissolved inorganic carbon and nitrate in groundwater at a former wood treatment facility containing hydrocarbon contamination

The speciation and mobility of a variety of chemical species in groundwater are strongly influenced by redox and pH conditions. Dissolved O₂ (DO) and dissolved inorganic C (DIC) concentrations are significant controls of these conditions, respectively. It is not always clear what the major processes are that influence changes in the concentration of DO and DIC across a groundwater flowpath. The combined use of the stable isotope compositions of DO (δ¹⁸O–DO) and DIC (δ¹³C–DIC) has the potential to help investigators discriminate between sources and sinks of DO and DIC in groundwater systems.     

Carbon, sulfur, and oxygen isotope evidence for a strong depth gradient and oceanic oxidation after the Ediacaran Hankalchough glaciation

In order to understand spatial variations of stable isotope geochemistry in the Quruqtagh basin (northwestern China) in the aftermath of an Ediacaran glaciation, we analyzed carbonate carbon isotopes (δ¹³C_carb), carbonate oxygen isotopes (δ¹⁸O_carb), carbonate associated sulfate sulfur (δ³⁴S_CAS) and oxygen isotopes (δ¹⁸O_CAS), and pyrite sulfur isotopes (δ³⁴S_py) of a cap dolostone atop the Ediacaran Hankalchough glacial diamictite at four sections. The four studied sections (YKG, MK, H and ZBS) represent an onshore-offshore transect in the Quruqtagh basin. Our data show a strong paleobathymetry-dependent isotopic gradient. From the onshore to offshore sections, δ¹³C_carb values decrease from −2‰ to −16‰ (VPDB), whereas δ¹⁸O_carb values increase from −4‰ to −1‰ (VPDB). Both δ³⁴S_CAS and δ³⁴S_py show stratigraphic variations in the two onshore sections (MK and YKG), but are more stable in the two offshore sections (H and ZBS). δ¹⁸O_CAS values of onshore samples are consistent with terrestrial oxidative weathering of pyrite. We propose that following the Hankalchough glaciation seawater in the Quruqtagh basin was characterized by a strong isotopic gradient. The isotopic data may be interpreted using a three-component mixing model that involves three reservoirs: deep-basin water, surface water, and terrestrial weathering input. In this model, the negative δ¹³C_carb values in the offshore sections are related to the upwelling of deep-basin water (where anaerobic oxidation of dissolved organic carbon resulted in ¹³C-depleted DIC), whereas sulfur isotope variations are strongly controlled by surface water sulfate and terrestrial weathering input derived from oxidative weathering of pyrite. The new data provide evidence for the oceanic oxidation following the Hankalchough glaciation.     

Calcium isotopes in a proglacial weathering environment: Damma glacier, Switzerland

The biogeochemical cycling and isotopic fractionation of calcium during the initial stages of weathering were investigated in an alpine soil chronosequence (Damma glacier, Switzerland). This site has a homogeneous silicate lithology and minimal biological impacts due to sparse vegetation cover. Calcium isotopic compositions, obtained by TIMS using a ⁴³Ca-⁴⁶Ca double spike, were measured in the main Ca pools. During this very early stage of weathering, the young soils which have formed (δ^44/42Ca=+0.44‰) were indistinguishable to the rocks from which they were derived (δ^44/42Ca=+0.44‰) and stream water (δ^44/42Ca=+0.48‰) was also within error of the average rock. This lack of variation indicates that the dissolution of the bulk silicate rock does not strongly fractionate Ca isotopes. The only Ca pool which was strongly fractionated from bulk rock was vegetation, which exhibited an enrichment of light Ca isotopes. Significant Ca isotope fractionation between bulk rock and the dissolved flux of Ca is likely to only occur where the Ca biogeochemical cycle is dominated by secondary processes such as biological cycling, adsorption and secondary mineral precipitation.     

Controls on calcium isotope fractionation in cultured planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera

Specimens of two species of planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera, were grown under controlled laboratory conditions at a range of temperatures (18-31 °C), salinities (32-44 psu) and pH levels (7.9-8.4). The shells were examined for their calcium isotope compositions (δ^44/40Ca) and strontium to calcium ratios (Sr/Ca) using Thermal Ionization Mass Spectrometry and Inductively Coupled Plasma Mass Spectrometry. Although the total variation in δ^44/40Ca (∼0.3‰) in the studied species is on the same order as the external reproducibility, the data set reveals some apparent trends that are controlled by more than one environmental parameter. There is a well-defined inverse linear relationship between δ^44/40Ca and Sr/Ca in all experiments, suggesting similar controls on these proxies in foraminiferal calcite independent of species. Analogous to recent results from inorganically precipitated calcite, we suggest that Ca isotope fractionation and Sr partitioning in planktic foraminifera are mainly controlled by precipitation kinetics. This postulation provides us with a unique tool to calculate precipitation rates and draws support from the observation that Sr/Ca ratios are positively correlated with average growth rates. At 25 °C water temperature, precipitation rates in G. siphonifera and G. ruber are calculated to be on the order of 2000 and 3000 μmol/m²/h, respectively. The lower δ^44/40Ca observed at ?29 °C in both species is consistent with increased precipitation rates at high water temperatures. Salinity response of δ^44/40Ca (and Sr/Ca) in G. siphonifera implies that this species has the highest precipitation rates at the salinity of its natural habitat, whereas increasing salinities appear to trigger higher precipitation rates in G. ruber. Isotope effects that cannot be explained by precipitation rate in planktic foraminifera can be explained by a biological control, related to a vacuolar pathway for supply of ions during biomineralization and a pH regulation mechanism in these vacuoles. In case of an additional pathway via cross-membrane transport, supplying light Ca for calcification, the δ^44/40Ca of the reservoir is constrained as −0.2‰ relative to seawater. Using a Rayleigh distillation model, we calculate that calcification occurs in a semi-open system, where less than half of the Ca supplied by vacuolization is utilized for calcite precipitation. Our findings are relevant for interpreting paleo-proxy data on δ^44/40Ca and Sr/Ca in foraminifera as well as understanding their biomineralization processes.     

δB, Sr, Mg and B in a modern Porites coral: the relationship between calcification site pH and skeletal chemistry

Sr/Ca, B/Ca, Mg/Ca and δ¹¹B were determined at high spatial resolution across ∼1 year of a modern Hawaiian Porites lobata coral by secondary ion mass spectrometry (SIMS). We observe significant variations in B/Ca, Mg/Ca, Sr/Ca and δ¹¹B over short skeletal distances (nominally equivalent to periods of <20 days). This heterogeneity probably reflects variations in the composition of the extracellular calcifying fluid (ECF) from which the skeleton precipitates. Calcification site pH (total scale) was estimated from skeletal δ¹¹B and ranged from 8.3 to 8.8 (± ∼0.1) with a mean of ∼8.6. Sr/Ca and B/Ca heterogeneity is not simply correlated with calcification site pH, as might be expected if Ca-ATPase activity increases the pH and decreases the Sr/Ca and B(OH)₄⁻/CO₃²⁻ ratios of the ECF. We produced a simple model of the ECF composition and the skeleton deposited from it, over a range of calcium transport and carbonate scenarios, which can account for these observed geochemical variations. The relationship between the pH and Sr/Ca of the ECF is dependent on the concentration of DIC at the calcification site. At higher DIC concentrations the ECF has a high capacity to buffer the [H⁺] changes induced by Ca-ATPase pumping. Conversely, at low DIC concentrations, this buffering capacity is reduced and ECF pH changes more rapidly in response to Ca-ATPase pumping. The absence of a simple correlation between ECF pH and skeletal Sr/Ca implies that calcification occurred under a range of DIC concentrations, reflecting variations in the respiration and photosynthesis of the coral and symbiotic zooxanthellate in the overlying coral tissues. Our observations have important implications for the use of coral skeletons as indicators of palaeo-ocean pH.     

Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Lithium elemental and isotopic variations in rock-melt interaction

Despite the occurrence of highly variable lithium (Li) elemental distribution and isotopic fractionation in mantle mineral, the mechanism of Li heterogeneity and fractionation remains a controversial issue. We measured Li contents and isotopic compositions of olivine and clinopyroxene xenocrysts and phenocrysts from kamafugite host lavas, as well as minerals in melt pockets occurring as metasomatic products in peridotite xenoliths from the Western Qinling, central China. The olivine xenocrysts in the kamafugites show compositional zonation. The cores have high Mg# (100 × Mg/(Mg+Fe); 91.0–92.2) and Li abundances (5.63–21.7 ppm), low CaO contents (≤0.12 wt%) and low δ⁷Li values (−39.6 to −6.76‰), which overlap with the compositional ranges of the olivines in the melt pockets as well as those in peridotite xenoliths. The rims of the olivine xenocrysts display relatively low Mg# (85.9–88.2), high CaO contents (0.19–0.38 wt%) and high δ⁷Li values (18.3–26.9‰), which are comparable to the olivine phenocrysts (Mg#: 86.4–87.1; CaO: 0.20–0.28 wt%; Li: 12.4–36.8 ppm; δ⁷Li: 18.1–26.0‰) and the silicate-melt metasomatized olivines. The clinopyroxene phenocrysts and clinopyroxenes in the melt pockets have no distinct characteristics with respect to the Li abundances and δ⁷Li values, but show higher and lower CaO contents, respectively, than the clinopyroxenes from silicate and carbonatite metasomatized samples. These features indicate that Li concentration and isotopic signatures of the cores of the xenocrysts recorded carbonatite melt-peridotite reaction (carbonatite metasomatism) at mantle depth, and the variations in the rims probably resulted from xenocryst–host magma interaction during ascent. Our results reveal that the interaction with carbonatite and silicate melts gave rise to an increase in Li abundance in minerals of peridotite xenoliths at mantle depth or during transportation. In terms of δ⁷Li, the carbonatite and silicate melts produced remarkably contrasting δ⁷Li variations in olivine. Based on the systematic variations of Li abundances and Li isotopes in olivines, we suggest that the δ⁷Li value of olivine is a more important indicator than that of clinopyroxene in discriminating carbonatite and silicate melt interaction agents with peridotites.     

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

Evaluation of petroleum hydrocarbon biodegradation in shallow groundwater by hydrogeochemical indicators and C, S-isotopes

Biodegradation is one of the main natural attenuation processes in groundwater contaminated with petroleum hydrocarbons. In this work, preliminary studies have been carried out by analyzing the concentrations of total petroleum hydrocarbons (TPH), dissolved inorganic carbon (DIC), dominant terminal electron accepters or donors, as well as δ ¹³C_DIC and δ ³⁴S_SO4, to reveal the biodegradation mechanism of petroleum hydrocarbons in a contaminated site. The results show that along groundwater flow in the central line of the plume, the concentrations of electron acceptors, pH, and E _h increased but TPH and DIC decreased. The δ ¹³C_DIC values of the contaminated groundwater were in the range of ?14.02 to ?22.28 ‰_PDB and ?7.71 to 8.36 ‰_PDB, which reflected a significant depletion and enrichment of ¹³C, respectively. The increase of DIC is believed to result from the non-methanogenic and methanogenic biodegradation of petroleum hydrocarbon in groundwater. Meanwhile, from the contaminated source to the downgradient of the plume, the ³⁴S in the contaminated groundwater became more depleted. The Rayleigh model calculation confirmed the occurrence of bacterial sulfate reduction as a biodegradation pathway of the petroleum hydrocarbon in the contaminated aquifers. It was concluded that stable isotope measurements, combined with other biogeochemical measurements, can be a useful tool to prove the occurrence of the biodegradation process and to identify the dominant terminal electron-accepting process in contaminated aquifers.