Paired Measurements of Foraminiferal d18O and Mg/Ca Ratios of Indian Monsoons Reconstructed from Holocene to Last Glacial Record

The effect of seasonally reversing monsoons in the northern Indian Ocean is to impart significant changes in surface salinity (SS). Here, we report SS changes during the last 32 kyr in the Lakshadweep Sea (southeastern Arabian Sea) estimated from paired measurements of d18O and sea surface temperature (SST) using Globigerinoides sacculifer, an upper mixed layer dwelling foraminifera. The heaviest d18OG.sacculifer (–0.07±0.08‰) is recorded between 23 and 15 ka, which could be defined as the last glacial maximum (LGM). The d18OG.sacculifer shift between the LGM and Holocene is 2.07‰. The SST shows an overall warming of 2°C from the LGM to Holocene (28°C to 30°C). However, coldest SSTs are observed prior to LGM, i.e., ~27 ka. The SS was higher (~38 psu) throughout most of the recorded last glacial period (32.5–15 ka). This high salinity together with generally lower SSTs suggests a period of sustained weaker summer or stronger winter monsoons. The deglacial warming is associated with rapid reorganization of monsoons and is reflected in decreased salinity to a modern level of ~ 36.5 psu, within a period of ~5 kyr. This indicates intensification of summer monsoons during cold to warm climate transition.     

Climatology of the Eastern Arabian Sea during the last glacial cycle reconstructed from paired measurement of foraminiferal δO and Mg/Ca

Paired measurements of Mg/Ca and δ¹⁸O of Globigerenoides sacculifer from an Eastern Arabian Sea (EAS) sediment core indicate that sea-surface temperature (SST) varied within 2°C and sea-surface salinity within 2 psu during the last 100 ka. SST was coldest (∼ 27°C) during Marine Isotope Stage (MIS) 4 and 2. Sea-surface salinity was highest (∼ 37.5 psu) during most of the last glacial period (∼ 60–18 ka), concurrent with increased δ¹⁸O_G.sacculifer and C/N ratios of organic matter and indicative of sustained intense winter monsoons. SST time series are influenced by both Greenland and Antarctic climates. However, the sea-surface salinity time series and the deglacial warming in the SST record (beginning at ∼ 18 ka) compare well with the LR04 benthic δ¹⁸O-stack and Antarctic temperatures. This suggests a teleconnection between the climate in the Southern Hemisphere and the EAS. Therefore, the last 100-ka variability in EAS climatology appears to have evolved in response to a combination of global climatic forcings and regional monsoons. The most intense summer monsoons within the Holocene occurred at ∼ 8 ka and are marked by SST cooling of ∼ 1°C, sea-surface salinity decrease of 0.5 psu, and δ¹⁸O_G.sacculifer decrease of 0.2‰.     

Oceanography in northwestern Europe during the last interglacial from intrashell δO ranges in Littorina littorea gastropods

Coastal sea-surface temperature (SST) and sea-surface salinity (SSS), including seasonality, in northwest (NW) Europe during the early phase of the Eemian interglacial ca. 125 ka ago were reconstructed from Littorina littorea (common periwinkle) gastropods. The results were based on intra-annual δ¹⁸O analyses in recent and fossil shells, mainly originating from the sea of Kattegat (Sweden) and the English Channel (United Kingdom), and confined to intertidal settings. The Eemian L. littorea shells indicated annual SSTs in the range 8–18°C for the English Channel and 8–26°C for Kattegat. All specimens from the Eemian sites experienced summer SSTs of ca. 1–3°C above recent conditions. The estimated winter SST in the English Channel during the Eemian was comparable to modern measurements of ca. 8°C. However, the Kattegat region displayed Eemian winter SST approximately 8°C warmer than today, and similar to conditions in the western English Channel. The recent-fossil isotope analogue approach indicated high SSS above 35 practical salinity units (psu) for a channel south of England in full contact with the North Atlantic Ocean during the last interglacial. In addition, the Kattegat shells indicated a SSS of ca. 29 psu, which points out a North Sea affinity for this region during the Eemian.     

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.     

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

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

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.     

Continental-scale patterns in modern wood cellulose δO: Implications for interpreting paleo-wood cellulose δO

The isotopic composition of ancient wood may be a useful archive of some climatic or geochemical conditions of the past, but presently there are many uncertainties that constrain such interpretations. We sampled naturally growing, predominantly native trees in forested regions of North America and the Caribbean to evaluate the strength of the relationships among cellulose δ¹⁸O (δ¹⁸O_cel), relative humidity (RH), precipitation δ¹⁸O (δ¹⁸O_ppt), and mean annual temperature (MAT) at the continental scale, and the general range of variability in δ¹⁸O_cel associated with site hydrologic conditions and species differences. We found up to 4‰ differences among different species growing at the same site, that conifer cellulose at a site is more enriched than angiosperm cellulose by 1.5‰ (p < 0.001), and that differences in landscape position, reflecting differing access to the water table, produced differences of <1‰ in δ¹⁸O_cel. At the continental scale, δ¹⁸O_cel was strongly influenced by modeled δ¹⁸O_ppt (R² = 0.80, p < 0.001). Average summer minimum RH (RH_min) combined with δ¹⁸O_ppt explained more of the variability (R² = 0.93, p < 0.001) in δ¹⁸O_cel across North American and Caribbean forests. MAT and δ¹⁸O_cel were also strongly correlated across North America (R = 0.91 and 0.95, p < 0.001, for angiosperms and conifers, respectively). The difference between δ¹⁸O_ppt and δ¹⁸O_cel is not constant (varying from 35-44‰) and is inversely correlated with δ¹⁸O_ppt. The relationships among δ¹⁸O_ppt, RH_min, δ¹⁸O_cel, and MAT established for North America and the Caribbean applied reasonably well when δ¹⁸O_cel was used to estimate MAT and δ¹⁸O_ppt in Asia, Europe, and South America, but there were important exceptions. The most accurate predictions of MAT and δ¹⁸O_ppt from δ¹⁸O_cel require RH_min. Predictions of δ¹⁸O_ppt and MAT made from δ¹⁸O_cel alone produced errors of up to 8‰ and 16 °C, respectively.     

Oxygen and iron isotope constraints on near-surface fractionation effects and the composition of lunar mare basalt source regions

Oxygen and iron isotope analyses of low-Ti and high-Ti mare basalts are presented to constrain their petrogenesis and to assess stable isotope variations within lunar mantle sources. An internally-consistent dataset of oxygen isotope compositions of mare basalts encompasses five types of low-Ti basalts from the Apollo 12 and 15 missions and eight types of high-Ti basalts from the Apollo 11 and 17 missions. High-precision whole-rock δ¹⁸O values (referenced to VSMOW) of low-Ti and high-Ti basalts correlate with major-element compositions (Mg#, TiO₂, Al₂O₃). The observed oxygen isotope variations within low-Ti and high-Ti basalts are consistent with crystal fractionation and match the results of mass-balance models assuming equilibrium crystallization. Whole-rock δ⁵⁶Fe values (referenced to IRMM-014) of high-Ti and low-Ti basalts range from 0.134‰ to 0.217‰ and 0.038‰ to 0.104‰, respectively. Iron isotope compositions of both low-Ti and high-Ti basalts do not correlate with indices of crystal fractionation, possibly owing to small mineral-melt iron fractionation factors anticipated under lunar reducing conditions.The δ¹⁸O and δ⁵⁶Fe values of low-Ti and the least differentiated high-Ti mare basalts are negatively correlated, which reflects their different mantle source characteristics (e.g., the presence or absence of ilmenite). The average δ⁵⁶Fe values of low-Ti basalts (0.073 ± 0.018‰, n = 8) and high-Ti basalts (0.191 ± 0.020‰, n = 7) may directly record that of their parent mantle sources. Oxygen isotope compositions of mantle sources of low-Ti and high-Ti basalts are calculated using existing models of lunar magma ocean crystallization and mixing, the estimated equilibrium mantle olivine δ¹⁸O value, and equilibrium oxygen-fractionation between olivine and other mineral phases. The differences between the calculated whole-rock δ¹⁸O values for source regions, 5.57‰ for low-Ti and 5.30‰ for high-Ti mare basalt mantle source regions, are solely a function of the assumed source mineralogy. The oxygen and iron isotope compositions of lunar upper mantle can be approximated using these mantle source values. The δ¹⁸O and δ⁵⁶Fe values of the lunar upper mantle are estimated to be 5.5 ± 0.2‰ (2σ) and 0.085 ± 0.040‰ (2σ), respectively. The oxygen isotope composition of lunar upper mantle is identical to the current estimate of Earth’s upper mantle (5.5 ± 0.2‰), and the iron isotope composition of the lunar upper mantle overlaps within uncertainty of estimates for the terrestrial upper mantle (0.044 ± 0.030‰).     

SIMS analyses of silicon and oxygen isotope ratios for quartz from Archean and Paleoproterozoic banded iron formations

Banded iron formations (BIFs) are chemical marine sediments dominantly composed of alternating iron-rich (oxide, carbonate, sulfide) and silicon-rich (chert, jasper) layers. Isotope ratios of iron, carbon, and sulfur in BIF iron-bearing minerals are biosignatures that reflect microbial cycling for these elements in BIFs. While much attention has focused on iron, banded iron formations are equally banded silica formations. Thus, silicon isotope ratios for quartz can provide insight on the sources and cycling of silicon in BIFs. BIFs are banded by definition, and microlaminae, or sub-mm banding, are characteristic of many BIFs. In situ microanalysis including secondary ion mass spectrometry is well-suited for analyzing such small features. In this study we used a CAMECA IMS-1280 ion microprobe to obtain highly accurate (±0.3‰) and spatially resolved (∼10 μm spot size) analyses of silicon and oxygen isotope ratios for quartz from several well known BIFs: Isua, southwest Greenland (∼3.8 Ga); Hamersley Group, Western Australia (∼2.5 Ga); Transvaal Group, South Africa (∼2.5 Ga); and Biwabik Iron Formation, Minnesota, USA (∼1.9 Ga). Values of δ¹⁸O range from +7.9‰ to +27.5‰ and include the highest reported δ¹⁸O values for BIF quartz. Values of δ³⁰Si have a range of ∼5‰ from −3.7‰ to +1.2‰ and extend to the lowest δ³⁰Si values for Precambrian cherts. Isua BIF samples are homogeneous in δ¹⁸O to ±0.3‰ at mm- to cm-scale, but are heterogeneous in δ³⁰Si up to 3‰, similar to the range in δ³⁰Si found in BIFs that have not experienced high temperature metamorphism (up to 300 °C). Values of δ³⁰Si for quartz are homogeneous to ±0.3‰ in individual sub-mm laminae, but vary by up to 3‰ between multiple laminae over mm-to-cm of vertical banding. The scale of exchange for Si in quartz in BIFs is thus limited to the size of microlaminae, or less than ∼1 mm. We interpret differences in δ³⁰Si between microlaminae as preserved from primary deposition. Silicon in BIF quartz is mostly of marine hydrothermal origin (δ³⁰Si < −0.5‰) but silicon from continental weathering (δ³⁰Si ∼ 1‰) was an important source as early as 3.8 Ga.     

Silicon isotope and trace element constraints on the origin of ∼3.5 Ga cherts: Implications for Early Archaean marine environments

Silicon (Si) isotope variability in Precambrian chert deposits is significant, but proposed explanations for the observed heterogeneity are incomplete in terms of silica provenance and fractionation mechanisms involved. To address these issues we investigated Si isotope systematics, in conjunction with geochemical and mineralogical data, in three well-characterised and approximately contemporaneous, ∼3.5 Ga chert units from the Pilbara greenstone terrane (Western Australia).We show that Si isotope variation in these cherts is large (−2.4‰ to +1.3‰) and was induced by near-surface processes that were controlled by ambient conditions. Cherts that formed by chemical precipitation of silica show the largest spread in δ³⁰Si (−2.4‰ to +0.6‰) and are characterised by positive Eu, La and Y anomalies and overall depletions in lithophile trace elements. Silicon isotope systematics in these orthochemical deposits are explained by (1) mixing between hydrothermal fluids and seawater, and/or (2) fractionation of hydrothermal fluids by subsurface losses of silica due to conductive cooling. Rayleigh-type fractionation of hydrothermal fluids was largely controlled by temperature differences between these fluids and seawater. Lamina-scale Si isotope heterogeneity within individual chemical chert samples up to 2.2‰ is considered to reflect the dynamic nature of hydrothermal activity. Silicified volcanogenic sediments lack diagnostic REE+Y anomalies, are enriched in lithophile elements, and exhibit a much more restricted range of positive δ³⁰Si (+0.1‰ to +1.1‰), which points to seawater as the dominant source of silica.The proposed model for Si isotope variability in the Early Archaean implies that chemical cherts with the most negative δ³⁰Si formed from pristine hydrothermal fluids, whereas silicified or chemical sediments with positive δ³⁰Si are closest to pure seawater deposits. Taking the most positive value found in this study (+1.3‰), and assuming that the Si isotope composition of seawater is governed by input of fractionated hydrothermal fluids, we infer that the temperature of ∼3.5 Ga seawater was below ∼55 °C.     

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.     

Cl,N, C isotopic analysis of common agro-chemicals for identifying non-point source agricultural contaminants

The isotopic compositions of commercially available herbicides were analyzed to determine their respective ¹⁵N, ¹³C and ³⁷Cl signatures for the purposes of developing a discrete tool for tracing and identifying non-point source contaminants in agricultural watersheds. Findings demonstrate that of the agrochemicals evaluated, chlorine stable isotopes signatures range between δ³⁷Cl = −4.55‰ and +3.40‰, whereas most naturally occurring chlorine stable isotopes signatures, including those of road salt, sewage sludge and fertilizers, vary in a narrow range about the Standard Mean Ocean Chloride (SMOC) between −2.00‰ and +1.00‰. Nitrogen stable isotope values varied widely from δ¹⁵N = −10.86‰ to +1.44‰ and carbon stable isotope analysis gave an observed range between δ¹³C = −37.13‰ and −21.35‰ for the entire suite of agro-chemicals analyzed. When nitrogen, carbon and chlorine stable isotope analyses were compared in a cross-correlation analysis, statistically independent isotopic signatures exist suggesting a new potential tracer tool for identifying herbicides in the environment.     

Soils at the hyperarid margin: The isotopic composition of soil carbonate from the Atacama Desert, Northern Chile

We evaluate the impact of exceptionally sparse plant cover (0-20%) and rainfall (2-114 mm/yr) on the stable carbon and oxygen composition of soil carbonate along elevation transects in what is among the driest places on the planet, the Atacama Desert in northern Chile. δ¹³C and δ¹⁸O values of carbonates from the Atacama are the highest of any desert in the world. δ¹³C (VPDB) values from soil carbonate range from −8.2‰ at the wettest sites to +7.9‰ at the driest. We measured plant composition and modeled respiration rates required to form these carbonate isotopic values using a modified version of the soil diffusion model of [Cerling (1984) Earth Planet. Sci. Lett.71, 229-240], in which we assumed an exponential form of the soil CO₂ production function, and relatively shallow (20-30 cm) average production depths. Overall, we find that respiration rates are the main predictor of the δ¹³C value of soil carbonate in the Atacama, whereas the fraction C₃ to C₄ biomass at individual sites has a subordinate influence. The high average δ¹³C value (+4.1‰) of carbonate from the driest study sites indicates it formed—perhaps abiotically—in the presence of pure atmospheric CO₂.δ¹⁸O (VPDB) values from soil carbonate range from −5.9‰ at the wettest sites to +7.3‰ at the driest and show much less regular variation with elevation change than δ¹³C values. δ¹⁸O values for soil carbonate predicted from local temperature and δ¹⁸O values of rainfall values suggest that extreme (>80% in some cases) soil dewatering by evaporation occurs at most sites prior to carbonate formation. The effects of evaporation compromise the use of δ¹⁸O values from ancient soil carbonate to reconstruct paleoelevation in such arid settings.     

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.     

O and S isotopic composition of dissolved and attached oxidation products of pyrite by Acidithiobacillus ferrooxidans: Comparison with abiotic oxidations

The acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans, plays a part in the pyrite oxidation process and has been widely studied in order to determine the kinetics of the reactions and the isotopic composition of dissolved product sulphates, but the details of the oxidation processes at the surface of pyrite are still poorly known. In this study, oxygen and sulphur isotopic compositions (δ¹⁸O and δ³⁴S) were analyzed for dissolved sulphates and water from experimental aerobic acidic (pH < 2) pyrite oxidation by A. ferrooxidans. The oxidation products attached to the pyrite surfaces were studied for their morphology (SEM), their chemistry (Raman spectroscopy) and for their δ¹⁸O (ion microprobe). They were compared to abiotically (Fe³⁺, H₂O₂, O₂) oxidized pyrite surface compounds in order to constrain the oxidation pathways and to look for the existence of potential biosignatures for this system.The pyrite dissolution evolved from non-stoichiometric (during the first days) to stoichiometric (with increasing time) resulting in dissolved sulphates having distinct δ¹⁸O (e.g. +11.0‰ and −2.0‰, respectively) and δ³⁴S (+4.5‰ and +2.8‰, respectively) values. The “oxidation layer” at the surface of pyrite is complex and made of iron oxides, sulphate, polysulphide, elemental sulphur and polythionates. Bio- and Fe³⁺-oxidation favour the development of monophased micrometric bumps made of hematite or sulphate while other abiotic oxidation processes result in more variable oxidation products. The δ¹⁸O of these oxidation products at the surface of oxidized pyrites are strongly variable (from ≈−40‰ to ≈+30‰) for all experiments.Isotopic fractionation between sulphates and pyrite, Δ³⁴S_SO4-pyrite, is equal to −1.3‰ and +0.4‰ for sulphates formed by stoichiometric and non-stoichiometric processes, respectively. These two values likely reflect either a S-S or a Fe-S bond breaking process. The Δ¹⁸O_SO4-H2O and Δ¹⁸O_SO4-O2 are estimated to be ≈+16‰ and ≈−25‰, respectively. These values are higher than previously published data and may reflect biological effects. The large δ¹⁸O heterogeneity measured at the surfaces of oxidized pyrites, whatever the oxidant, may be related (i) to the existence of local surface environments isolated from the solution in which the oxidation processes are different and (ii) to the stabilization at the pyrite surface of reaction intermediates that are not in isotopic equilibrium with the solution. Though the oxygen isotopic composition of surface oxidation products cannot be taken as a direct biosignature, the combined morphological, chemical and isotopic characterization of the surfaces of oxidized pyrites may furnish clues about a biological activity on a mineral surface.     

Behaviour of Li and its isotopes during metasomatism of French Massif Central lherzolites

Li behaviour and distribution in the mantle were investigated by ion microprobe in situ measurements on co-existing olivine (ol), orthopyroxene (opx), clinopyroxene (cpx) and amphibole (amp) in xenoliths from the French Massif Central. The fertile spinel lherzolites of this study record increasing degrees of mantle metasomatism, from unmetasomatised anhydrous samples through cryptically metasomatised samples to highly metasomatised amphibole-rich samples. In anhydrous lherzolites, Li is preferentially incorporated into olivine (1.1-1.4 ppm, average values) compared to pyroxenes (0.2-0.9 ppm). The hydrous samples clearly show enrichment of Li in ol (1.5-5.0 ppm), opx (1.1-2.4 ppm) and cpx (2.4-5.4 ppm), while amphibole incorporates less Li than the co-existing phases (0.8-1.3 ppm). Average δ⁷Li values range from +7.6 to +14.5‰ in ol, from 5.1 to +13.7‰ in opx and from 8.8 to +10.3‰ in cpx from the anhydrous lherzolites. A layered peridotite sample (Sdi) shows higher Li content in all phases, with lighter isotopic composition in opx and cpx (−0.6 and −2‰ average δ⁷Li values, respectively). In the hydrous lherzolites average δ⁷Li values both overlap and extend beyond these ranges in ol (up to 17.5 ‰) and in opx (up to 22.9‰), and vary widely in cpx (−2.7 to +9.7‰). Low δ⁷Li values are observed in some opx (−10.4‰) and cpx (−13‰) from sample Sdi, and in cpx from three hydrous samples (from −9.7 to −5.3‰). The different anhydrous phases from the hydrous samples show large intra-grain variations in Li isotopic ratios (e.g., up to 18‰) compared to the same phases from the anhydrous samples (mostly less than 6‰), excepting sample Sdi which has up to 20.4‰ variation in cpx. Similar to the anhydrous silicates, amphiboles show a wide variation of δ⁷Li values on the intra-grain scale (2-27‰). These variations are interpreted to result from fractionation processes during metasomatism by a silicate melt undergoing compositional changes as it percolates through and reacts with the peridotite phases. Thus Li abundances and isotopic in situ measurements are useful for tracing metasomatic processes but the heterogeneities observed in the samples preclude any identification of a specific mantle source by its Li signature.     

Metamorphic growth and recrystallization of zircons in extremely O-depleted rocks during eclogite-facies metamorphism: Evidence from U-Pb ages, trace elements, and O-Hf isotopes

A combined in situ SIMS and LA-(MC)-ICPMS study of U-Pb ages, trace elements, O and Lu-Hf isotopes was conducted for zircon from eclogite-facies metamorphic rocks in the Sulu orogen. The two microbeam techniques sampled various depths of zircon domains, revealing different element and isotope relationships between residual magmatic cores and new metamorphic rims and thus the geochemical architecture of metamorphic zircons which otherwise cannot be recognized by the single microbeam technique. This enables discrimination of metamorphic growth from different subtypes of metamorphic recrystallization. Magmatic cores with U-Pb ages of 769 ± 9 Ma have positive δ¹⁸O values of 0.1-10.1‰, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high REE contents, and steep MREE-HREE patterns with negative Eu anomalies. They are interpreted as crystallizing from positive δ¹⁸O magmas during protolith emplacement. In contrast, newly grown domains have concordant U-Pb ages of 204 ± 4 to 252 ± 7 Ma and show negative δ¹⁸O values of −10.0‰ to −2.2‰, low Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, low REE contents, and flat HREE patterns with weak to no Eu anomalies. They are interpreted as growing from negative δ¹⁸O fluids that were produced by metamorphic dehydration of high-T glacial-hydrothermally altered rocks during continental subduction-zone metamorphism. Differences in δ¹⁸O between different domains within single grains vary from 0.8‰ to 12.5‰, suggesting different degrees of O isotope exchange between the positive δ¹⁸O magmatic core and the negative δ¹⁸O metamorphic fluid during the metamorphism. The magmatic zircons underwent three subtypes of metamorphic recrystallization, depending on their accessibility to negative δ¹⁸O fluids. The zircons recrystallized in solid-state maintained positive δ¹⁸O values, and REE and Lu-Hf isotopes of protolith zircon, but their U-Pb ages are lowered. The zircons recrystallized through dissolution exhibit negative δ¹⁸O values similar to the metamorphic growths, almost completely reset U-Pb ages, and partially reset REE systems. The zircons recrystallized through replacement show variably negative δ¹⁸O values, and partially reset REE, and U-Pb and Lu-Hf isotopic systems. Therefore, this study places robust constraints on the origin of metamorphic zircons in eclogite-facies rocks and provides a methodological framework for linking the different types of metamorphic zircons to petrological processes during continental collision.     

High precision SIMS oxygen three isotope study of chondrules in LL3 chondrites: Role of ambient gas during chondrule formation

We report high precision SIMS oxygen three isotope analyses of 36 chondrules from some of the least equilibrated LL3 chondrites, and find systematic variations in oxygen isotope ratios with chondrule types. FeO-poor (type I) chondrules generally plot along a mass dependent fractionation line (Δ¹⁷O ∼ 0.7‰), with δ¹⁸O values lower in olivine-rich (IA) than pyroxene-rich (IB) chondrules. Data from FeO-rich (type II) chondrules show a limited range of δ¹⁸O and δ¹⁷O values at δ¹⁸O = 4.5‰, δ¹⁷O = 2.9‰, and Δ¹⁷O = 0.5‰, which is slightly ¹⁶O-enriched relative to bulk LL chondrites (Δ¹⁷O ∼ 1.3‰). Data from four chondrules show ¹⁶O-rich oxygen isotope ratios that plot near the CCAM (Carbonaceous Chondrite Anhydrous Mineral) line. Glass analyses in selected chondrules are systematically higher than co-existing minerals in both δ¹⁸O and Δ¹⁷O values, whereas high-Ca pyroxene data in the same chondrule are similar to those in olivine and pyroxene phenocrysts.Our results suggest that the LL chondrite chondrule-forming region contained two kinds of solid precursors, (1) ¹⁶O-poor precursors with Δ¹⁷O > 1.6‰ and (2) ¹⁶O-rich solid precursors derived from the same oxygen isotope reservoir as carbonaceous chondrites. Oxygen isotopes exhibited open system behavior during chondrule formation, and the interaction between the solid and ambient gas might occur as described in the following model. Significant evaporation and recondensation of solid precursors caused a large mass-dependent fractionation due to either kinetic or equilibrium isotope exchange between gas and solid to form type IA chondrules with higher bulk Mg/Si ratios. Type II chondrules formed under elevated dust/gas ratios and with water ice in the precursors, in which the ambient H₂O gas homogenized chondrule melts by isotope exchange. Low temperature oxygen isotope exchange may have occurred between chondrule glasses and aqueous fluids with high Δ¹⁷O (∼5‰) in LL the parent body. According to our model, oxygen isotope ratios of chondrules were strongly influenced by the local solid precursors in the proto-planetary disk and the ambient gas during chondrule melting events.     

Stable isotope geochemistry of ground and surface waters associated with undisturbed massive sulfide deposits; constraints on origin of waters and water-rock reactions

Stable isotopes (H, O, C) were determined for ground and surface waters collected from two relatively undisturbed massive sulfide deposits (Halfmile Lake and Restigouche) in the Bathurst Mining Camp (BMC), New Brunswick, Canada. Additional waters from active and inactive mines in the BMC were also collected. Oxygen and hydrogen isotopes of surface and shallow groundwaters from both the Halfmile Lake and Restigouche deposits are remarkably uniform (− 13 to − 14‰ and − 85 to − 95‰ for δ¹⁸O_VSMOW and δ²H_VSMOW, respectively). These values are lighter than predicted for northern New Brunswick and, combined with elevated deuterium excess values, suggest that recharge waters are dominated by winter precipitation, recharged during spring melting. Deeper groundwaters from the Restigouche deposit, and from active and inactive mines have heavier δ¹⁸O_VSMOW ratios (up to − 10.8‰) than shallow groundwaters suggesting recharge under warmer climate or mixing with Shield-type brines. Some of the co-variation in Cl concentrations and δ¹⁸O_VSMOW ratios can be explained by mixing between saline and shallow recharge water end-members. Carbon isotopic compositions of dissolved inorganic carbon (DIC) are variable, ranging from − 15 to − 5‰ δ¹³C_VPDB for most ground and surface waters. Much of the variation in the carbon isotopes is consistent with closed system groundwater evolution involving soil zone CO₂ and fracture zone carbonate minerals (calcite, dolomite and siderite; average = − 6.5‰ δ¹³C_VPDB). The DIC of saline Restigouche deposit groundwater is isotopically heavy (∼+ 12‰ δ¹³C_VPDB), indicating carbon isotopic fractionation from methanogenesis via CO₂ reduction, consistent with the lack of dissolved sulfate in these waters and the observation of CH₄-degassing during sampling.     

Boron isotopic composition of atmospheric precipitations and liquid-vapour fractionations

Boron isotope compositions (δ¹¹B) and B concentrations of rains and snows were studied in order to characterize the sources and fractionation processes during the boron atmospheric cycle. The ¹¹B/¹⁰B ratios of instantaneous and cumulative rains and snows from coastal and continental sites show a large range of variations, from −1.5 ± 0.4 to +26.0 ± 0.5‰ and from −10.2 ± 0.5 to +34.4 ± 0.2‰, respectively. Boron concentrations in rains and snows vary between 0.1 and 3.0 ppb. All these precipitation samples are enriched in ¹⁰B compared to the ocean value (δ¹¹B = +39.5‰). An empirical rain-vapour isotopic fractionation of +31‰ is estimated from three largely independent methods. The deduced seawater-vapour fractionation is +25.5‰, with the difference between the rain and seawater fractionations principally reflecting changes in the speciation of boron in the liquid with ∼100% B(OH)₃ present in precipitations. A boron meteoric water line, δD = 2.6δ¹¹B − 133, is proposed which describes the relationship between δD and δ¹¹B in many, but not all, precipitations. Boron isotopic compositions of precipitations can be related to that of the seawater reservoir by the seawater-vapour fractionation and one or more of (1) the rain-vapour isotopic fractionation, (2) evolution of the δ¹¹B value of the atmospheric vapour reservoir via condensation-precipitation processes (Rayleigh distillation process), (3) any contribution of vapour from the evaporation of seawater aerosols, and (4) any contribution from particulate matter, principally sea salt, continental dust and, perhaps more regionally, anthropogenic sources (burning of biomass and fossil fuels). From the δ¹¹B values of continental precipitations, a sea salt contribution cannot be more than a percent or so of the total B in precipitation over these areas.