Carbon mass-balance modelling and carbon isotope exchange processes in dynamic caves

Diverse interpretations have been made of carbon isotope time series in speleothems, reflecting multiple potential controls. Here we study the dynamics of ¹³C and ¹²C cycling in a particularly well-constrained site to improve our understanding of processes affecting speleothem δ¹³C values. The small, tubular Grotta di Ernesto cave (NE Italy) hosts annually-laminated speleothem archives of climatic and environmental changes. Temperature, air pressure, pCO₂, dissolved inorganic carbon (DIC) and their C isotopic compositions were monitored for up to five years in soil water and gas, cave dripwater and cave air. Mass-balance models were constructed for CO₂ concentrations and tested against the carbon isotope data. Air advection forces winter pCO₂ to drop in the cave air to ca. 500 ppm from a summer peak of ca. 1500 ppm, with a rate of air exchange between cave and free atmosphere of approximately 0.4 days. The process of cave ventilation forces degassing of CO₂ from the dripwater, prior to any calcite precipitation onto the stalagmites. This phase of degassing causes kinetic isotope fractionation, i.e. ¹³C-enrichment of dripwater whose δ¹³C_DIC values are already higher (by about 1‰) than those of soil water due to dissolution of the carbonate rock. A subsequent systematic shift to even higher δ¹³C values, from −11.5‰ in the cave drips to about −8‰ calculated for the solution film on top of stalagmites, is related to degassing on the stalagmite top and equilibration with the cave air. Mass-balance modelling of C fluxes reveals that a very small percentage of isotopically depleted cave air CO₂ evolves from the first phase of dripwater degassing, and shifts the winter cave air composition toward slightly more depleted values than those calculated for equilibrium. The systematic ¹³C-enrichment from the soil to the stalagmites at Grotta di Ernesto is independent of drip rate, and forced by the difference in pCO₂ between cave water and cave air. This implies that speleothem δ¹³C values may not be simply interpreted either in terms of hydrology or soil processes.     

Speleothem calcite farmed in situ: Modern calibration of δO and δC paleoclimate proxies in a continuously-monitored natural cave system

Understanding the relationships between speleothem stable isotopes (δ¹³C δ¹⁸O) and in situ cave forcing mechanisms is important to interpreting ancient stalagmite paleoclimate records. Cave studies have demonstrated that the δ¹⁸O of inorganically precipitated (low temperature) speleothem calcite is systematically heavier than the δ¹⁸O of laboratory-grown calcite for a given temperature. To understand this apparent offset, rainwater, cave drip water, groundwater, and modern naturally precipitated calcite (farmed in situ) were grown at multiple locations inside Hollow Ridge Cave in Marianna, Florida. High resolution micrometeorological, air chemistry time series and ventilation regimes were also monitored continuously at two locations inside the cave, supplemented with periodic bi-monthly air gas grab sample transects throughout the cave.Cave air chemistry and isotope monitoring reveal density-driven airflow pathways through Hollow Ridge Cave at velocities of up to 1.2 m s⁻¹ in winter and 0.4 m s⁻¹ in summer. Hollow Ridge Cave displays a strong ventilation gradient in the front of the cave near the entrances, resulting in cave air that is a mixture of soil gas and atmospheric CO₂. A clear relationship is found between calcite δ¹³C and cave air ventilation rates estimated by proxies pCO₂ and ²²²Rn. Calcite δ¹³C decreased linearly with distance from the front entrance to the interior of the cave during all seasons, with a maximum entrance-to-interior gradient of Δδ¹³C_CaCO3 = −7‰. A whole-cave “Hendy test” at multiple contemporaneous farming sites reveals that ventilation induces a +1.9 ± 0.96‰ δ¹³C offset between calcite precipitated in a ventilation flow path and calcite precipitated on the edge or out of flow paths. This interpretation of the “Hendy test” has implications for interpreting δ¹³C records in ancient speleothems. Calcite δ¹³C_CaCO3 may be a proxy not only for atmospheric CO₂ or overlying vegetation shifts but also for changes in cave ventilation due to dissolution fissures and ceiling collapse creating and plugging ventilation windows.Farmed calcite δ¹⁸O was found to exhibit a +0.82 ± 0.24‰ offset from values predicted by both theoretical calculations and laboratory-grown inorganic calcite. Unlike δ¹³C_CaCO3, oxygen isotopes showed no ventilation effects, i.e. Δδ¹⁸O_CaCO3 appears to be a function of growth temperature only although we cannot rule out a small effect of (unmeasured) gradients in relative humidity (evaporation) accompanying ventilation. Our results support the findings of other cave investigators that water-calcite fractionation factors observed in speleothem calcite are higher that those measured in laboratory experiments. Cave and laboratory calcite precipitates may differ mainly in the complex effects of kinetic isotope fractionation. Combining our data with other recent speleothem studies, we find a new empirical relationship for cave-specific water-calcite oxygen isotope fractionation across a range of temperatures and cave environments:

1000lnα=16.1(10³T^-1)-24.6     

Controls on dissolved inorganic carbon and δC in cave waters from DeSoto Caverns: Implications for speleothem δC assessments

Unraveling the factors controlling the carbon chemistry and transport of carbon within extant karst systems has important implications concerning the assessment of time-series δ¹³C records of speleothems. Here we report the results of a 3-year study of total dissolved inorganic carbon [DIC] and δ¹³C_DIC from cave waters at DeSoto Caverns (Southeastern USA) that offer valuable insight on carbon transport and the accompanied isotope fractionations from end-member sources to speleothems.[DIC] and δ¹³C_DIC values of cave waters range from 0.2 to 6.0 mM and 2.7 to −12.9 (‰ VPDB), respectively. [DIC] and δ¹³C_DIC of “seasonal drips” show seasonal, albeit noisy, variability and are inversely related (δ¹³C_DIC = −2.49[DIC] + 0.64, r² = 0.84). A shallow pool fed by multiple drips shows a bimodal δ¹³C_DIC distribution with an isotopically heavier mode during winter (−4‰ to −5‰ VPDB) relative to summer months (−9‰ to −10‰ VPDB). A multi-year trend of decreasing water availability during the study period is not reflected in a response of cave water carbon chemistry suggesting that rainfall amount may not be a significant controlling factor of the carbon chemistry. Coupled cave air winter ventilation/summer stagnation and varying CO₂ fluxes through the soil horizon and epikarst exert the strongest influence on seasonal [DIC] and δ¹³C_DIC variability. Measured values of high [DIC] and low δ¹³C_DIC from cave waters collected during the summer/early fall closely approximate isotopic equilibrium conditions. Conversely, low [DIC] and high δ¹³C_DIC values during winter/early months indicate kinetically enhanced isotopic fractionations within the cave waters. The kinetically enhanced isotopic fractionation of partitioned between degassed CO₂ and precipitated CaCO₃(1000lnα_{[(CO2-HCO3)+(CaCO3(AR)-HCO3)]/2}) is greater by about a factor of two (−6.7 ± 0.3‰) relative to the same isotopic fractionation under equilibrium conditions (−3.1‰).On the basis of ¹⁴C mass balance and paired ¹⁴C-U/Th measurements we estimate that on average about ∼23% of C delivered annually by the drips to the aragonite stalagmites is derived from ¹⁴C-dead dolomite cap while the remainder of ∼77% is derived from ¹⁴C-live biomass. δ¹³C measurements of aragonite (n = 12) sampled from the tips of active speleothems during the summer months are consistent with theoretical aragonite δ¹³C values calculated using the shallow pool summer/early fall data thus confirming the δ¹³C seasonality in both drips and coeval aragonite. δ¹³C values of an active stalagmite section spanning the last 200 years show a normal distribution with a mean of −7.1 ± 1.2‰ (n = 81) and a mode of −7‰ to −8‰ that are statistically indistinguishable from the annual mean and mode of all dripwaters. Thus secular time-series δ¹³C records of stalagmites at DeSoto Caverns with resolving power >10⁻¹ year will likely carry the imprints of drip annual means that record climate-driven δ¹³C seasonal biases.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

Magnesium stable isotope fractionation in marine biogenic calcite and aragonite

This survey of magnesium stable isotope compositions in marine biogenic aragonite and calcite includes samples from corals, sclerosponges, benthic porcelaneous and planktonic perforate foraminifera, coccolith oozes, red algae, and an echinoid and brachiopod test. The analyses were carried out using MC-ICP-MS with an external repeatability of ±0.22‰ (2SD for δ²⁶Mg; n = 37), obtained from a coral reference sample (JCp-1).Magnesium isotope fractionation in calcitic corals and sclerosponges agrees with published data for calcitic speleothems with an average Δ²⁶Mg_{calcite-seawater} = −2.6 ± 0.3‰ that appears to be weakly related to temperature. With one exception (Vaceletia spp.), aragonitic corals and sclerosponges also display uniform Mg isotope fractionations relative to seawater with Δ²⁶Mg_{biogenic aragonite-seawater} = −0.9 ± 0.2.Magnesium isotopes in high-Mg calcites from red algae, echinoids and perhaps some porcelaneous foraminifera as well as in all low-Mg calcites (perforate foraminifera, coccoliths and brachiopods) display significant biological influences. For planktonic foraminifera, the Mg isotope data is consistent with the fixation of Mg by organic material under equilibrium conditions, but appears to be inconsistent with Mg removal from vacuoles. Our preferred model, however, suggests that planktonic foraminifera synthesize biomolecules that increase the energetic barrier for Mg incorporation. In this model, the need to remove large quantities of Mg from vacuole solutions is avoided. For the high-Mg calcites from echinoids, the precipitation of amorphous calcium carbonate may be responsible for their weaker Mg isotope fractionation.Disregarding superimposed biological effects, it appears that cation light isotope enrichments in CaCO₃ principally result from a chemical kinetic isotope effect, related to the incorporation of cations at kink sites. In this model, the systematics of cation isotope fractionations in CaCO₃ relate to the activation energy required for cation incorporation, which probably reflects the dehydration of the cation and the crystal surface and bond formation at the incorporation site. This kinetic incorporation model predicts (i) no intrinsic dependence on growth rate, unless significant back reaction upon slow growth reduces the isotope fractionation towards that characteristic for equilibrium isotope partitioning (this may be observed for Ca isotopes in calcites), (ii) a small decrease of isotope fractionation with increasing temperature that may be amplified if higher temperatures promote back reaction and (iii) a sensitivity to changes in the activation barrier caused by additives such as anions or biomolecules or by the initial formation of amorphous CaCO₃.     

Seasonal dripwater Mg/Ca and Sr/Ca variations driven by cave ventilation: Implications for and modeling of speleothem paleoclimate records

A 4-year study in a central Texas cave quantifies multiple mechanisms that control dripwater composition and how these mechanisms vary at different drip sites. We monitored cave-air compositions, in situ calcite growth, dripwater composition and drip rate every 4-6 weeks. Three groups of drip sites are delineated (Groups 1-3) based on geochemical variations in dripwater composition. Quantitative modeling of mineral-solution reactions within the host carbonate rock and cave environments is used to identify mechanisms that can account for variations in dripwater compositions. The covariation of Mg/Ca (and Sr/Ca) and Sr isotopes is key in delineating whether Mg/Ca and Sr/Ca variations are dictated by water-rock interaction (i.e., calcite or dolomite recrystallization) or prior calcite precipitation (PCP). Group 1 dripwater compositions reflects a narrow range of the extent of water-rock interaction followed by varying amounts of prior calcite precipitation (PCP). Group 2 dripwater compositions are controlled by varying amounts of water-rock interaction with little to no PCP influence. Group 3 dripwater compositions are dictated by variable extents of both water-rock interaction and PCP. Group 1 drip sites show seasonal variations in dripwater Mg/Ca and Sr/Ca, whereas the other drip sites do not. In contrast to the findings of most previous dripwater Mg/Ca-Sr/Ca studies, these seasonal variations (at Group 1 drip sites) are independent of changes in water flux (i.e., rainfall and/or drip rate), and instead significantly correlate with changes in cave-air CO₂ concentrations. These results are consistent with lower cave-air CO₂, related to cool season ventilation of the cave atmosphere, enhancing calcite precipitation and leading to dripwater geochemical evolution via PCP. Group 1 dripwater Mg/Ca and Sr/Ca seasonality and evidence for PCP as a mechanism that can account for that seasonality, have two implications for many other regions where seasonal ventilation of caves is likely: (1) speleothem trace-element records may provide seasonal signals, and (2) such records may be biased toward recording climate conditions during the season when calcite is depositing. Additionally, we use our results to construct a forward model that illustrates the types of speleothem Mg/Ca and Sr/Ca variations that would result from varying controls on dripwater compositions. The model provides a basis for interpreting paleo-dripwater controls from high frequency Mg/Ca and Sr/Ca variations for speleothems from caves at which long term monitoring studies are not feasible.     

Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths

We present whole rock Li and Mg isotope analyses of 33 ultramafic xenoliths from the terrestrial mantle, which we compare with analyses of 30 (mostly chondritic) meteorites. The accuracy of our new Mg isotope ratio measurement protocol is substantiated by a combination of standard addition experiments, the absence of mass independent effects in terrestrial samples and our obtaining identical values for rock standards using two different separation chemistries and three different mass-spectrometric introduction systems. Carbonaceous, ordinary and enstatite chondrites have irresolvable mean stable Mg isotopic compositions (δ²⁵Mg = −0.14 ± 0.06; δ²⁶Mg = −0.27 ± 0.12‰, 2SD), but our enstatite chondrite samples have lighter δ⁷Li (by up to ∼3‰) than our mean carbonaceous and ordinary chondrites (3.0 ± 1.5‰, 2SD), possibly as a result of spallation in the early solar system. Measurements of equilibrated, fertile peridotites give mean values of δ⁷Li = 3.5 ± 0.5‰, δ²⁵Mg = −0.10 ± 0.03‰ and δ²⁶Mg = −0.21 ± 0.07‰. We believe these values provide a useful estimate of the primitive mantle and they are within error of our average of bulk carbonaceous and ordinary chondrites. A fuller range of fresh, terrestrial, ultramafic samples, covering a variety of geological histories, show a broad positive correlation between bulk δ⁷Li and δ²⁶Mg, which vary from −3.7‰ to +14.5‰, and −0.36‰ to + 0.06‰, respectively. Values of δ⁷Li and δ²⁶Mg lower than our estimate of primitive mantle are strongly linked to kinetic isotope fractionation, occurring during transport of the mantle xenoliths. We suggest Mg and Li diffusion into the xenoliths is coupled to H loss from nominally anhydrous minerals following degassing. Diffusion models suggest that the co-variation of Mg and Li isotopes requires comparable diffusivities of Li and Mg in olivine. The isotopically lightest samples require ∼5-10 years of diffusive ingress, which we interpret as a time since volatile loss in the host magma. Xenoliths erupted in pyroclastic flows appear to have retained their mantle isotope ratios, likely as a result of little prior degassing in these explosive events. High δ⁷Li, coupled with high [Li], in rapidly cooled arc peridotites may indicate that these samples represent fragments of mantle wedge that has been metasomatised by heavy, slab-derived fluids. If such material is typically stirred back into the convecting mantle, it may account for the heavy δ⁷Li seen in some oceanic basalts.     

Oxygen and carbon isotopic systematics of aragonite speleothems and water in Furong Cave, Chongqing, China

To understand oxygen and carbon stable isotopic characteristics of aragonite stalagmites and evaluate their applicability to paleoclimate, the isotopic compositions of active and fossil aragonite speleothems and water samples from an in situ multi-year (October 2005-July 2010) monitoring program in Furong Cave located in Chongqing of China have been examined. The observations during October 2005-June 2007 show that the meteoric water is well mixed in the overlying 300-500-m bedrock aquifer, reflected by relatively constant δ¹⁸O, ±0.11-0.14‰ (1σ), of drip waters in the cave, which represents the annual status of rainfall water. Active cave aragonite speleothems are at oxygen isotopic equilibrium with drip water and their δ¹⁸O values capture the surface-water oxygen isotopic signal. Aragonite-to-calcite transformation since the last glaciation is not noticeable in Furong stalagmites. Our multi-year field experiment approves that aragonite stalagmite δ¹⁸O records in this cave are suitable for paleoclimate reconstruction. With high U, 0.5-7.2 ppm, and low Th, 20-1270 ppt, the Furong aragonite stalagmites provide very precise chronology (as good as ±20s yrs (2σ)) of the climatic variations since the last deglaciation. The synchroneity of Chinese stalagmite δ¹⁸O records at the transition into the Bølling-Allerød (t-BA) and the Younger Dryas from Furong, Hulu and Dongge Caves supports the fidelity of the reconstructed East Asian monsoon evolution. However, the Furong record shows that the cold Older Dryas (OD) occurred at 14.0 thousand years ago, agreeing with Greenland ice core δ¹⁸O records but ∼200 yrs younger than that in the Hulu record. The OD age discrepancy between Chinese caves can be attributable to different regionally climatic/environmental conditions or chronological uncertainty of stalagmite proxy records, which is limited by changes in growth rate and subsampling intervals in absolute dating. Seasonal dissolved inorganic carbon δ¹³C variations of 2-3‰ in the drip water and 5-7‰ in the pool and spring waters are likely attributed to variable degrees of CO₂ degassing in winter and summer. The variable δ¹³C values of active deposits from −11‰ to 0‰ could be caused by kinetically mediated CO₂ degassing processes. The complicated nature of pre-deposition kinetic isotopic fractionation processes for carbon isotopes in speleothems at Furong Cave require further study before they can be interpreted in a paleoclimatic or paleoenvironmental context.     

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO₄) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (δ^44/40Ca = −2.01 ± 0.15‰) but are different from hydrothermal and cold seep barite samples (δ^44/40Ca = −4.13 to −2.72‰). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, Δ^44/40Ca = −3.42 to −2.40‰. Temperature, saturation state, , and aCa²⁺/aBa²⁺ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower ⁴⁴Ca/⁴⁰Ca than calcite, by −9‰ at 0 °C and −8‰ at 25 °C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower δ^44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca²⁺ substitution for Ba²⁺ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca²⁺ and other elements in inorganic and biogenic minerals.     

Mg isotope constraints on soil pore-fluid chemistry: Evidence from Santa Cruz, California

Mg isotope ratios (²⁶Mg/²⁴Mg) are reported in soil pore-fluids, rain and seawater, grass and smectite from a 90 kyr old soil, developed on an uplifted marine terrace from Santa Cruz, California. Rain water has an invariant ²⁶Mg/²⁴Mg ratio (expressed as δ²⁶Mg) at −0.79 ± 0.05‰, identical to seawater δ²⁶Mg. Detrital smectite (from the base of the soil profile, and therefore unweathered) has a δ²⁶Mg value of 0.11‰, potentially enriched in ²⁶Mg by up to 0.3‰ compared to the bulk silicate Earth Mg isotope composition (although within the range of all terrestrial silicates). The soil pore-waters show a continuous profile with depth for δ²⁶Mg, ranging from −0.99‰ near the surface to −0.43‰ at the base of the profile. Shallow pore-waters (<1 m) have δ²⁶Mg values that are similar to, or slightly lower than the rain waters. This implies that the degree of biological cycling of Mg in the pore-waters is relatively small and is quantified as <32%, calculated using the average Mg isotope enrichment factor between grass and rain (δ²⁶Mg_grass-δ²⁶Mg_rain) of 0.21‰. The deep pore-waters (1-15 m deep) have δ²⁶Mg values that are intermediate between the smectite and rain, ranging from −0.76‰ to −0.43‰, and show a similar trend with depth compared to Sr isotope ratios. The similarity between Sr and Mg isotope ratios confirms that the Mg in the pore-waters can be explained by a mixture between rain and smectite derived Mg, despite the fact that Mg and Sr concentrations may be buffered by the exchangeable reservoir. However, whilst Sr isotope ratios in the pore-waters span almost the complete range between mineral and rain inputs, Mg isotopes compositions are much closer to the rain inputs. If Mg and Sr isotope ratios are controlled uniquely by a mixture, the data can be used to estimate the mineral weathering inputs to the pore-waters, by correcting for the rain inputs. This isotopic correction is compared to the commonly used chloride correction for precipitation inputs. A consistent interpretation is only possible if Mg isotope ratios are fractionated either by the precipitation of a secondary Mg bearing phase, not detected by conventional methods, or selective leaching of ²⁴Mg from smectite. There is therefore dual control on the Mg isotopic composition of the pore-waters, mixing of two inputs with distinct isotopic compositions, modified by fractionation. The data provide (1) further evidence for Mg isotope fractionation at the surface of the Earth and (2) the first field evidence of Mg isotope fractionation during uptake by natural plants. The coherent behaviour of Mg isotope ratios in soil environments is encouraging for the development of Mg isotope ratios as a quantitative tracer of both weathering inputs of Mg to waters, and the physicochemical processes that cycle Mg, a major cation linked to the carbon cycle, during continental weathering.     

Processes controlling dripwater hydrochemistry variations in Xueyu Cave,SW China: implications for speleothem palaeoclimate signal interpretations

Cave dripwater hydrochemistry responds to environmental changes, both within and outside of the cave, and thereby conveys this information to any stalagmite fed by the drips. As stalagmites are important archives of climate proxy information, understanding how dripwater hydrochemistry responds to environmental forcing is critical. However, despite the large number of speleothems in SW China, the response of dripwater to regional climate variability is not yet adequately understood. A 3‐year study of three drip sites in Xueyu Cave, Chongqing Municipality, SW China, revealed the most important mechanisms controlling dripwater chemical variability. The principal chemical indices (pH, specific conductivity, Ca²⁺, Mg²⁺, Sr²⁺ and ) in collected dripwaters and the local climate data were analysed in this study. The principal controls on the hydrochemistry were found to be the external climate and its changes, groundwater residence time, cave ventilation and prior calcite precipitation (PCP) processes. Dripwater hydrochemistry showed strongly coherent seasonal patterns despite the fact that all sites are Ca–HCO₃ type waters and supersaturated with calcite. Seasonal changes in dripwater hydrochemistry were influenced by the soil and vadose zone CO₂ content as well as groundwater residence time in the upper karst zone. Cave‐air CO₂ seasonal variations were consistent with changes in dripwater PCO₂ and cave ventilation. Trace element ratios (Mg/Ca and Sr/Ca) of dripwater were controlled by PCP processes. Seasonal variations in dripwater Mg/Ca and Sr/Ca ratios in Xueyu Cave showed inverse changes with the Asia Monsoon Index during the monitoring period, reflecting the seasonal climate changes that may have been recorded in the speleothems. Based on a linear regression of PCO₂ and the Ca²⁺ data in the cold–dry winter season, a 130‐ppm shift in cave‐air PCO₂ results in a 1‐ppm shift in dripwater Ca²⁺ concentration in Xueyu Cave. This study illustrates the importance of understanding factors controlling the changes in the composition of dripwater before it reaches the speleothem.     

Response of the coccolithophores Emiliania huxleyi and Coccolithus braarudii to changing seawater Mg and Ca concentrations: Mg/Ca, Sr/Ca ratios and δCa, δMg of coccolith calcite

Calcium and magnesium concentrations in seawater have varied over geological time scales. On short time scales, variations in the major ion composition of seawater influences coccolithophorid physiology and the chemistry of biogenically produced coccoliths. Validation of those changes via controlled laboratory experiments is a crucial step in applying coccolithophorid based paleoproxies for the reconstruction of past environmental conditions. Therefore, we examined the response of two species of coccolithophores, Emiliania huxleyi and Coccolithus braarudii, to changes in the seawater Mg/Ca ratio (≈0.5 to 10 mol/mol) by either manipulating the magnesium or calcium concentration under controlled laboratory conditions. Concurrently, seawater Sr/Ca ratios were also modified (≈2 to 40 mmol/mol), while keeping salinity constant at 35. The physiological response was monitored by measurements of the cell growth rate as well as the production rates of particulate inorganic and organic carbon, and chlorophyll a. Additionally, coccolithophorid calcite was analyzed for its elemental composition (Sr/Ca and Mg/Ca) as well as isotope fractionation of calcium and magnesium (Δ^44/40Ca and Δ^26/24Mg). Our results reveal that physiological rates were substantially influenced by changes in seawater calcium rather than magnesium concentration within the range estimated to have occurred over the past 250 million years when coccolithophores appear in the fossil record. All physiological rates of E. huxleyi decreased at a calcium concentration above 25 mmol L⁻¹, whereas C. braarudii displayed a higher tolerance to increased seawater calcium concentrations. Partition coefficient of Sr was calculated as 0.36 ± 0.04 (±2σ) independent of species. Partition coefficient of Mg²⁺ increased with increasing seawater Ca²⁺ concentrations in both coccolithophore species. Calcium isotope fractionation was constant at 1.1 ± 0.1‰ (±2σ) and not altered by changes in seawater Mg/Ca ratio. There is a well-defined inverse linear relationship between calcium isotope fractionation and partition coefficient of Sr²⁺ in all experiments, suggesting similar controls on both proxies in the investigated species. Magnesium isotope ratios were relatively stable for seawater Mg/Ca ratios ranging from 1 to 5, with a higher degree of fractionation in Emiliania huxleyi (by ≈0.2‰ in Δ^26/24Mg). Although Mg/Ca ratios in the calcite of coccolithophores and foraminifera are similar, the former have considerably higher Δ^26/24Mg (by >+3‰), presumably due to differences in calcification mechanisms between the two taxa. These observations suggest, a physiological control over magnesium elemental and isotopic fractionation during the process of calcification in coccolithophores.     

Magnesium and iron isotopes in 2.7 Ga Alexo komatiites: Mantle signatures, no evidence for Soret diffusion, and identification of diffusive transport in zoned olivine

Komatiites from Alexo, Canada, are well preserved and represent high-degree partial mantle melts (∼50%). They are thus well suited for investigating the Mg and Fe isotopic compositions of the Archean mantle and the conditions of magmatic differentiation in komatiitic lavas. High precision Mg and Fe isotopic analyses of 22 samples taken along a 15-m depth profile in a komatiite flow are reported. The δ²⁵Mg and δ²⁶Mg values of the bulk flow are −0.138 ± 0.021‰ and −0.275 ± 0.042‰, respectively. These values are indistinguishable from those measured in mantle peridotites and chondrites, and represent the best estimate of the composition of the silicate Earth from analysis of volcanic rocks. Excluding the samples affected by secondary Fe mobilization, the δ⁵⁶Fe and δ⁵⁷Fe values of the bulk flow are +0.044 ± 0.030‰, and +0.059 ± 0.044‰, respectively. These values are consistent with a near-chondritic Fe isotopic composition of the silicate Earth and minor fractionation during komatiite magma genesis. In order to explain the early crystallization of pigeonite relative to augite in slowly cooled spinifex lavas, it was suggested that magmas trapped in the crystal mush during spinifex growth differentiated by Soret effect, which should be associated with large and coupled variations in the isotopic compositions of Mg and Fe. The lack of variations in Mg and Fe isotopic ratios either rules out the Soret effect in the komatiite flow or the effect is effaced as the solidification front migrates downward through the flow crust. Olivine separated from a cumulate sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg values relative to the bulk flow, which modeling shows can be explained by kinetic isotope fractionation associated with Fe-Mg inter-diffusion in olivine. Such variations can be used to identify diffusive processes involved in the formation of zoned minerals.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.     

Carbon and hydrogen isotope fractionation associated with the aerobic microbial oxidation of methane, ethane, propane and butane

Carbon isotope fractionation factors associated with the aerobic consumption of methane (C1), ethane (C2), propane (C3), and n-butane (C4) were determined from incubations of marine sediment collected from the Coal Oil Point hydrocarbon seep field, located offshore Santa Barbara, CA. Hydrogen isotope fractionation factors for C1, C2 and C3 were determined concurrently. Fresh sediment samples from two seep areas were each slurried with sea water and treated with C1, C2, C3 or C4, or with mixtures of all four gases. Triplicate samples were incubated aerobically at 15 °C, and the stable isotope composition and headspace levels of C1-C4 were monitored over the course of the experiment. Oxidation was observed for all C1-C4 gases, with an apparent preference for C3 and C4 over C1 and C2 in the mixed-gas treatments. Fractionation factors were calculated using a Rayleigh model by comparing the δ¹³C and δD of the residual C1-C4 gases to their headspace levels. Carbon isotope fractionation factors (reported in ε or (α-1) × 1000 notation) were consistent between seep areas and were −26.5‰ ± 3.9 for C1, −8.0‰ ± 1.7 for C2, −4.8‰ ± 0.9 for C3 and −2.9‰ ± 0.9 for C4. Fractionation factors determined from mixed gas incubations were similar to those determined from individual gas incubations, though greater variability was observed during C1 consumption. In the case of C1 and C3 consumption, carbon isotope fractionation appears to decrease as substrate becomes limiting. Hydrogen isotope fractionation factors determined from the two seep areas differed for C1 oxidation but were similar for C2 and C3. Hydrogen isotope fractionation factors ranged from −319.9‰ to −156.4‰ for C1 incubations, and averaged −61.9‰ ± 8.3 for C2 incubations and −15.1‰ ± 1.9 for C3 incubations. The fractionation factors presented here may be applied to estimate the extent of C1-C4 oxidation in natural gas samples, and should prove useful in further studying the microbial oxidation of these compounds in the natural environment.     

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ⁶⁵Cu values (based on ⁶⁵Cu/⁶³Cu) of enargite (δ⁶⁵Cu = −0.01 ± 0.10‰; 2σ) and chalcopyrite (δ⁶⁵Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ⁶⁵Cu < 1‰). These mineral samples show lower δ⁶⁵Cu values than stream waters (1.38‰ ? δ⁶⁵Cu ? 1.69‰). The average isotopic fractionation (Δ_aq-min = δ⁶⁵Cu_aq − δ⁶⁵Cu_min, where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in ⁶⁵Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite (Δ_aq-min^o=-0.57±0.14‰, where min^o refers to the starting mineral) and no apparent fractionation for enargite (Δ_aq-min^o=0.14±0.14‰). Abiotic fractionation is attributed to preferential oxidation of ⁶⁵Cu⁺ at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of ⁶⁵Cu_aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ⁶⁵Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.     

Stable isotope variations in modern tropical speleothems: Evaluating equilibrium vs. kinetic isotope effects

Applications of speleothem calcite geochemistry in climate change studies require the evaluation of the accuracy and sensitivity of speleothem proxies to correctly infer paleoclimatic information. The present study of Harrison’s Cave, Barbados, uses the analysis of the modern climatology and groundwater system to evaluate controls on the C and O isotopic composition of modern speleothems. This new approach directly compares the δ¹⁸O and δ¹³C values of modern speleothems with the values for their corresponding drip waters in order to assess the degree to which isotopic equilibrium is achieved during calcite precipitation. If modern speleothems can be demonstrated to precipitate in isotopic equilibrium, then ancient speleothems, suitable for paleoclimatic studies, from the same cave environment may also have been precipitated in isotopic equilibrium. If modern speleothems are precipitated out of isotopic equilibrium, then the magnitude and direction of the C and O isotopic offsets may allow specific kinetic and/or equilibrium isotopic fractionation mechanisms to be identified.Carbon isotope values for the majority of modern speleothem samples from Harrison’s Cave fall within the range of equilibrium values predicted from the combined use of (1) calcite-water fractionation factors from the literature, (2) measured temperatures, and (3) measured δ¹³C values of the dissolved inorganic carbon of drip waters. Calcite samples range from ∼0.8‰ higher to ∼1.1‰ lower than predicted values. The ¹³C depletions are likely caused by kinetically driven departures in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions, caused by rapid calcite growth. ¹³C enrichments can be accounted for by Rayleigh distillation of the HCO₃⁻ (aq) reservoir during degassing of ¹³C-depleted CO₂.Modern speleothems from Harrison’s Cave are not in O isotopic equilibrium with their corresponding drip waters and are 0.2‰ to 2.3‰ enriched in ¹⁸O relative to equilibrium values. δ¹⁸O variations in modern calcite are likely controlled by kinetically driven changes in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions to nonequilibrium conditions, consistent with rapid calcite growth. In contrast to δ¹³C, δ¹⁸O values of modern calcite may not be affected by Rayleigh distillation during degassing because CO₂ hydration and hydroxylation reactions will buffer the O isotopic composition of the HCO₃⁻ (aq) reservoir. If the effects of Rayleigh distillation manifest themselves in the O isotopic system, they will result in ¹⁸O enrichment in the HCO₃⁻ (aq) reservoir and ultimately in the precipitated CaCO₃.     

Controls on the D/H ratios of plant leaf waxes in an arid ecosystem

The extent to which leaf water D-enrichment (transpiration) and soil water D-enrichment (evaporation) affect the D/H ratio of plant leaf waxes remains a contentious issue, with important implications for paleohydrologic reconstructions. In this study we measure δD values of precipitation (δD_p), groundwater (δD_gw), plant xylem water (δD_xw) and leaf water (δD_lw) to understand their impact on the δD values of plant leaf wax n-alkanes (δD_wax) in an arid ecosystem. Our survey includes multiple species at four sites across an aridity gradient (80-30% relative humidity) in southern California.We find that many species take up groundwater or precipitation without significant fractionation. D-enriched soil water is a minor source even in species known to perform and utilize waters from hydraulic lift, such as Larrea tridentata (+10‰). Measurements of leaf water isotopic composition demonstrate that transpiration is an important mechanism for D-enrichment of leaf waters (+74 ± 20‰, 1σ), resulting in the smallest net fractionation yet reported between source water and leaf waxes (L. tridentata −41‰; multi-species mean value is −94 ± 21‰, 1σ). We find little change in leaf water D-enrichment or net fractionation across the climatic gradient sampled by our study, suggesting that a net fractionation of ca. −90‰ may be appropriate for paleohydrologic reconstructions in semi-arid to arid environments. Large interspecies offsets in net fractionations (1σ = 21‰) are potentially troublesome, given the observed floristic diversity and the likelihood of species assemblage changes with climate shifts.     

An experimental study of magnesium-isotope fractionation in chlorophyll-a photosynthesis

Measurements are presented of the magnesium isotopic composition of chlorophyll-a, extracted from cyanobacteria, relative to the isotopic composition of the culture medium in which the cyanobacteria were grown. Yields of 50-93% chlorophyll-a were achieved from the pigment extracts of Synechococcus elongatus, a unicellular cyanobacteria. This material was then digested using concentrated nitric acid to extract magnesium. Separation was accomplished using columns of cation-exchange resin, which achieved a 103 ± 10% yield of magnesium from chlorophyll-a. This procedure ensured accurate measurement of the magnesium-isotopic ratios without isobaric interferences using a multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). We find a slight depletion in the heavier isotopes of magnesium in chlorophyll-a relative to culture medium, early growth phase: Δ²⁶Mg = −0.71(±0.35)‰ and Δ²⁵Mg = −0.37(±0.18)‰; late growth phase: Δ²⁶Mg = −0.53(±0.20)‰ and Δ²⁵Mg = −0.26(±0.11)‰, due to an apparent mass-dependent fractionation. We suggest that the small fractionation results from chelation during intracellular processes. A likely candidate for this chelation step involves the magnesium-chelatase enzyme, which mediates the insertion of magnesium to the tetrapyrrole ring during chlorophyll-a biosynthesis. Proof of this hypothesis can be tested with biological controls whereby steps in the enzymatic pathways of chlorophyll synthesis are selectively suppressed.