Microanalysis of carbon isotope composition in organic matter by secondary ion mass spectrometry

An analytical procedure has been developed for the in situ measurement of carbon isotope composition of organic matter, with a spatial resolution of 20-30 μm, using a Cameca IMS 1270 ion microprobe. Instrumental mass fractionation (IMF) of carbon isotopes was observed to be independent of primary ion beam intensity and sputtering time, but did depend on vacuum conditions and on the chemical composition of the sample. To evaluate such “matrix effects”, a set of 9 standards representative of the natural chemical variability of organic matter was prepared, with H/C atomic ratios and organic carbon contents (C_org) ranging between 0.04 and 1.74 and between 41 and 100 wt.%, respectively. Under the analytical conditions tested, IMF was not found to be influenced by the presence of silicate mineral impurities in the organic matter, but variations in IMF up to 5‰ were observed over the set of standards with the magnitude of IMF negatively correlated to the H/C ratios of samples. Aliphaticity ratios determined using Fourier transform infrared microspectroscopy provided an in situ estimation of H/C ratios with a spatial resolution barely exceeding that of the ion microprobe and permit a correction for matrix effects with a standard error of ± 0.2‰ (1σ). Taking into account all sources of uncertainty, ion microprobe δ¹³C were accurately determined with a ± 0.7‰ (1σ) total uncertainty. The mechanism for the matrix effect of H/C ratios upon IMF is still to be determined but it is likely related to the difference in proportion of atomic vs. molecular carbon ions observed between samples of different H/C ratios.     

New Reference Materials and Assessment of Matrix Effects for SIMS Measurements of Oxygen Isotopes in Garnet

Accurate ion microprobe analysis of oxygen isotope ratios in garnet requires appropriate reference materials to correct for instrumental mass fractionation that partly depends on the garnet chemistry (matrix effect). The matrix effect correlated with grossular, spessartine and andradite components was characterised for the Cameca IMS 1280HR at the SwissSIMS laboratory based on sixteen reference garnet samples. The correlations fit a second‐degree polynomial with maximum bias of ca. 4‰, 2‰ and 8‰, respectively. While the grossular composition range 0–25% is adequately covered by available reference materials, there is a paucity of them for intermediate compositions. We characterise three new garnet reference materials GRS2, GRS‐JH2 and CAP02 with a grossular content of 88.3 ± 1.2% (2s), 83.3 ± 0.8% and 32.5 ± 3.0%, respectively. Their micro scale homogeneity in oxygen isotope composition was evaluated by multiple SIMS sessions. The reference δ¹⁸O value was determined by CO₂ laser fluorination (δ¹⁸O_LF). GRS2 has δ¹⁸O_LF = 8.01 ± 0.10‰ (2s) and repeatability within each SIMS session of 0.30–0.60‰ (2s), GRS‐JH2 has δ¹⁸O_LF = 18.70 ± 0.08‰ and repeatability of 0.24–0.42‰ and CAP02 has δ¹⁸O_LF = 4.64 ± 0.16‰ and repeatability of 0.40–0.46‰.     

Isotopic constraints on ice age fluids in active geothermal systems: Reykjanes, Iceland

The Reykjanes geothermal system is located on the landward extension of the Mid-Atlantic Ridge in southwest Iceland, and provides an on-land proxy to high-temperature hydrothermal systems of oceanic spreading centers. Previous studies of elemental composition and salinity have shown that Reykjanes geothermal fluids are likely hydrothermally modified seawater. However, δD values of these fluids are as low as −23‰, which is indicative of a meteoric water component. Here we constrain the origin of Reykjanes hydrothermal solutions by analysis of hydrogen and oxygen isotope compositions of hydrothermal epidote from geothermal drillholes at depths between 1 and 3 km. δD_EPIDOTE values from wells RN-8, -9, -10 and -17 collectively range from −60 to −78‰, and δ¹⁸O_EPIDOTE in these wells are between −3.0 and 2.3‰. The δD values of epidote generally increase along a NE trend through the geothermal field, whereas δ¹⁸O values generally decrease, suggesting a southwest to northeast migration of the geothermal upflow zone with time that is consistent with present-day temperatures and observed hydrothermal mineral zones. For comparative analysis, the meteoric-water dominated Nesjavellir and Krafla geothermal systems, which have a δD_FLUID of ∼ −79‰ and −89‰, respectively, show δD_EPIDOTE values of −115‰ and −125‰. In contrast, δD_EPIDOTE from the mixed meteoric-seawater Svartsengi geothermal system is −68‰; comparable to δD_EPIDOTE from well RN-10 at Reykjanes.Stable isotope compositions of geothermal fluids in isotopic equilibrium with the epidotes at Reykjanes are computed using published temperature dependent hydrogen and oxygen isotope fractionation curves for epidote-water, measured isotope composition of the epidotes and temperatures approximated from the boiling point curve with depth. Calculated δD and δ¹⁸O of geothermal fluids are less than 0‰, suggesting that fluids of meteoric or glacial origin are a significant component of the geothermal solutions. Additionally, δD_FLUID values in equilibrium with geothermal epidote are lower than those of modern-day fluids, whereas calculated δ¹⁸O_FLUID values are within range of the observed fluid isotope composition. We propose that modern δD_EPIDOTE and δD_FLUID values are the result of diffusional exchange between hydrous alteration minerals that precipitated from glacially-derived fluids early in the evolution of the Reykjanes system and modern seawater-derived geothermal fluids. A simplified model of isotope exchange in the Reykjanes geothermal system, in which the average starting δD_ROCK value is −125‰ and the water to rock mass ratio is 0.25, predicts a δD_FLUID composition within 1‰ of average measured values. This model resolves the discrepancy between fluid salinity and isotope composition of Reykjanes geothermal fluids, explains the observed disequilibrium between modern fluids and hydrothermal epidote, and suggests that rock-fluid interaction is the dominant control over the evolution of fluid isotope composition in the hydrothermal system.     

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.     

Use of B isotopes as a tracer of anthropogenic emissions in the atmosphere of Paris,France

In order to test the potential of B isotopes as a tracer of contamination of the atmosphere, the B isotopic composition of rainwater samples monitored over a year in the centre of Paris, France were determined. Boron concentrations range from 19 nmol/L to 500 nmol/L and δ¹¹B range from 0‰ to +38‰. Mean annual values are 148 nmol/L and +25‰, respectively. The results suggest that variability in B isotopic compositions is mainly caused by mixing of two main sources, although isotopic fractionation during the evaporation–condensation processes may also be important. One source is a marine component, which exhibits a heavy B isotopic composition. The decrease of δ¹¹B in rainwater with increasing NO₃/B and SO₄/B molar ratios suggests that a second source may be anthropogenic emissions. To constrain this end-member, B was determined in urban particulates, which were enriched in the light isotope and the lowest values were consistent with a B-rich fossil fuel composition. These results confirm the great sensitivity of B to anthropogenic sources and the ability of B isotopic ratios to decipher the origin of B in the atmosphere.     

Inter-mineral Fe isotope variations in mantle-derived rocks and implications for the Fe geochemical cycle

Iron isotope and major- and minor-element compositions of coexisting olivine, clinopyroxene, and orthopyroxene from eight spinel peridotite mantle xenoliths; olivine, magnetite, amphibole, and biotite from four andesitic volcanic rocks; and garnet and clinopyroxene from seven garnet peridotite and eclogites have been measured to evaluate if inter-mineral Fe isotope fractionation occurs in high-temperature igneous and metamorphic minerals and if isotopic fractionation is related to equilibrium Fe isotope partitioning or a result of open-system behavior. There is no measurable fractionation between silicate minerals and magnetite in andesitic volcanic rocks, nor between olivine and orthopyroxene in spinel peridotite mantle xenoliths. There are some inter-mineral differences (up to 0.2‰ in ⁵⁶Fe/⁵⁴Fe) in the Fe isotope composition of coexisting olivine and clinopyroxene in spinel peridotites. The Fe isotope fractionation observed between clinopyroxene and olivine appears to be a result of open-system behavior based on a positive correlation between the Δ⁵⁶Fe_{clinopyroxene-olivine} fractionation and the δ⁵⁶Fe value of clinopyroxene and olivine. There is also a significant difference in the isotopic compositions of garnet and clinopyroxene in garnet peridotites and eclogites, where the average Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation is +0.32 ± 0.07‰ for six of the seven samples. The one sample that has a lower Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation of 0.08‰ has a low Ca content in garnet, which may reflect some crystal chemical control on Fe isotope fractionation. The Fe isotope variability in mantle-derived minerals is interpreted to reflect subduction of isotopically variable oceanic crust, followed by transport through metasomatic fluids. Isotopic variability in the mantle might also occur during crystal fractionation of basaltic magmas within the mantle if garnet is a liquidus phase. The isotopic variations in the mantle are apparently homogenized during melting processes, producing homogenous Fe isotope compositions during crust formation.     

Pore water evolution in oilfield sandstones: constraints from oxygen isotope microanalyses of quartz cement

Oxygen isotope microanalyses of authigenic quartz, in combination with temperatures of quartz precipitation constrained by fluid inclusion microthermometry and burial history modelling, are employed to trace the origin and evolution of pore waters in three distinct reservoirs of the Brae Formation in the Miller and Kingfisher Fields (North Sea). Oxygen isotope ratios of quartz cements were measured in situ in nine sandstone thin sections with a Cameca ims-4f ion microprobe. In conjunction with quartz cement paragenesis in the reservoirs, constrained from textural and cathodoluminescence (CL) microscopy studies, pore water evolution was reconstructed from the time of deposition of the sandstones in the Upper Jurassic until the present.CL photomicrographs of quartz overgrowths in the Brae Formation sandstones show three cement zones (A, B and C) which can be related to different oxygen isotope compositions: (1) the earliest, and thinnest, zone A (homogeneous CL pattern with probable δ¹⁸O values between +23‰ and +26‰—direct measurements were not possible) precipitated in the sandstones at temperatures <60 °C; (2) the second zone B (complex CL pattern and directly measured δ¹⁸O values between +15‰ and +18‰) precipitated in the sandstones most likely between 70 and 90 °C; (3) the third zone C (homogeneous CL pattern and directly measured δ¹⁸O values between +16‰ and +22‰) precipitated in the sandstones most likely at temperatures >90 °C. Calculated oxygen isotope compositions of pore waters show that zone A quartz cements, and enclosing concretionary calcite, precipitated from a meteoric-type fluid (∼−7‰) during shallow burial (<1.5 km). Zone B quartz cements precipitated from fluids which evolved in composition from a meteoric-type fluid (δ¹⁸O −7‰) to a more ¹⁸O-enriched fluid (δ¹⁸O −4‰) as burial continued to ∼3.0 km. Data from zone C quartz cements are consistent with further fluid evolution from δ¹⁸O −4‰ to basinal-type fluids with δ¹⁸O similar to the present-day formation water oxygen isotope composition (+0.6‰ at 4.0 km burial). A similar pore water evolution can be derived for all three reservoirs studied, indicating that hydrogeologic evolution was similar across sandstones of the whole Brae Formation.The quartz cement zones observed in the Brae Formation sandstones, and the pore water history derived for the area studied, is analogous to published petrographic and pore water evolution data from the nearby Brent Group reservoirs and from reservoirs located in the Haltenbanken area on the Atlantic margin offshore Norway. Considering quartz cement is a major porosity-occluding phase in many reservoir sandstones, and because pore waters both dissolve quartz and carry the dissolved silica to cementation sites, the data presented are valuable for improving the understanding and prediction of reservoir quality development in sandstones globally.     

Deciphering interaction of regional aquifers in Southern Tunisia using hydrochemistry and isotopic tools

Groundwater is the most important source of water supply in southern Tunisia. Previous hydrogeologic and isotopic studies carried out in this region revealed the existence of two major aquifer systems: the “Complex Terminal” (CT) and the “Continental Intercalaire” (CI). Turonian carbonates constitute one of the major aquifer levels of the CT multilayered aquifer. It extends over most of southern Tunisia, and its hydrodynamic regime is largely influenced by tectonics, lithology and recharge conditions. Forty-eight groundwater samples from the CI and Turonian aquifers were collected between January and April 2004 for chemical and isotopic analyses. Hydrochemistry and isotopic tools were combined to get an insight into the processes controlling chemical composition of groundwater and wide-scale interaction of these two aquifer systems. Analysis of the dissolved constituents revealed that several processes control the observed chemical composition: (i) incongruent dissolution of carbonate minerals, (ii) dissolution of evaporitic minerals, and (iii) cation exchange. Dissolution alone cannot account for the observed high supersaturation states of groundwater with respect to calcite and dolomite. The observed supersaturation is most probably linked to geogenic CO₂ entering water-bearing horizons of the CT and CI aquifers via deep tectonic faults and discontinuities and subsequent degassing in the exploitation wells. Presence of geogenic CO₂ in the investigated region was confirmed by C isotope data of the DIC reservoir. The radiocarbon content of the Turonian samples varied between 9.5 and 43 pmc. For CI samples generally lower values were recorded, between 3.8 and 22.5 pmc. Stable isotope composition of Turonian groundwater samples varied from −8.3 to −5.3‰ for δ¹⁸O and from −60 to −25‰ for δ²H. The corresponding ranges of δ values for the Continental Intercalaire samples were from −8.9‰ to −6.9‰ for δ¹⁸O and from −68.2‰ to −45.7‰ for δ²H. Stable isotope composition of groundwater representing CT and CI aquifers provide strong evidence for regional interaction between both systems.     

Sulfur and carbon isotope records from 1700 to 800 Ma carbonates of the Jixian section, northern China: Implications for secular isotope variations in Proterozoic seawater and relationships to global supercontinental events

This study is a comprehensive, stable isotope survey of the marine carbonate-dominated, upper Paleo- to lower Neoproterozoic stratigraphy of Jixian County, China. Carbonate-associated sulfate (CAS) was extracted and measured for δ³⁴S_CAS using the same samples analyzed for δ¹³C_carbonate. This integrated proxy approach is a step towards a more comprehensive picture of secular variation in the composition of Proterozoic seawater. We specifically sampled marine carbonate intervals from the lower section of the Chuanlinggou Formation, Changcheng Group (ca. 1700 Ma) to the top of the Jingeryu Formation, Qingbaikou Group (ca. 800 Ma). δ¹³C_carbonate values are mostly negative in the upper Paleoproterozoic Changcheng Group, with an ascending trend from −3‰ to 0‰. We observed variation of approximately 0 ± 1‰ in the Mesoproterozoic Jixian Group, and positive values of +2 ± 2‰ characterize the lower Neoproterozoic Qingbaikou Group. Stratigraphic variations in δ³⁴S_CAS are more remarkable in their ranges and magnitudes, including conspicuously high values exceeding +30‰ in the three intervals at ca. 1700 Ma, 1300-1100 Ma, and 1000-900 Ma. In the Changcheng Group, δ³⁴S_CAS values are typically higher than +25‰, with only a few values of less than +15‰. In contrast, most of the data spanning from the Mesoproterozoic Tieling Formation of the Jixian Group to the lower Neoproterozoic Jingeryu Formation of the Qingbaikou Group are highly variable between +10‰ and +25‰, with some values exceeding +25‰.In the late Paleoproterozoic (1700-1600 Ma), a >10‰ decrease in δ³⁴S_CAS and ∼3‰ increase in δ¹³C_carbonate are coincident with, and likely related to, the breakup of Columbia, a supercontinent that predated Rodinia. Carbon and sulfur isotope data from the Mesoproterozoic, when global tectonic activity was comparatively weaker, fall mostly in the ranges of +15 ± 10‰ and 0 ± 1‰, respectively, but fluctuations of >20‰ for δ³⁴S_CAS and >3‰ for the δ¹³C_carbonate at ca. 1450-1400 Ma may reflect subduction and large-scale magmatic activity in island arcs marking the end of Columbia breakup. From the late Mesoproterozoic (ca. 1300-1100 Ma) to the early Neoproterozoic (ca. 800 Ma), the δ¹³C and δ³⁴S of seawater increased gradually with increasing variability. Most impressive areδ³⁴S_CAS values that exceed +30‰ in two intervals at ca. 1300-1100 Ma and ca. 1000-900 Ma, which may reflect the assembly and early breakup of Rodinia. Although gaps in the record remain, and studies of even higher resolution are warranted, our results suggest that changes in paleoceanographic conditions linked to global tectonics strongly influenced the biogeochemical cycles of C and S. Furthermore, periods of the Proterozoic previously noted for their isotopic invariability show clear isotopic expressions of this tectonic activity.     

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.     

Carbon, oxygen, and strontium isotope geochemistry of carbonate rocks of the upper Miocene Kudankulam Formation, southern India: Implications for paleoenvironment and diagenesis

The carbon, oxygen, and strontium isotope compositions of carbonate rocks from the upper Miocene Kudankulam Formation, southern India, were measured to understand palaeoenvironment and carbonate diagenesis of this formation. Both carbon and oxygen isotope ratios of various carbonate phases including whole rocks, ooids, molluscan mold-fill and sparry pore-fill calcite cements are depleted in ¹⁸O and ¹³C compared to those of contemporaneous seawater, indicating that the Kudankulam carbonates underwent extensive meteoric diagenesis. Based on δ¹³C and δ¹⁸O values for sparry calcite cements (pore-fill and molluscan mold-fill) formed in the meteoric diagenetic realm (δ¹³C from −7.8‰ to −6.0‰ and −9.0‰ to −7.0‰; δ¹⁸O from −9.2‰ to −6.5‰ and −9.4‰ to −2.6‰, respectively), it is interpreted that the diagenetic system was open and was proximal to the vadose water recharge zone. The negative δ¹⁸O values of various carbonate components (about −9.4‰ to −4.1‰ for whole rocks; about −8.4‰ to −2.6‰ for ooids) suggest that during the late Miocene the paleoclimate of the study area was humid, unlike today, probably due to the intense Indian monsoon system. The carbon isotope compositions (−7.9‰ to −3.6‰ for whole rocks; −4.9‰ to −1.5‰ for ooids) are consistent with the interpretation that the paleo-ecosystem comprised a significant proportion of C₄ type plants, supporting a scenario of expansion of C₄ plants during the late Miocene in the Indian subcontinent as far south as the southern tip of India. The ⁸⁷Sr/⁸⁶Sr ratios of the Kudankulam carbonates (0.70920 to 0.72130) are much greater than those of the contemporaneous or modern seawater (between 0.7089 and 0.7091) and show a general decrease up-sequence. Such high Sr isotope ratios indicate significant radiogenic ⁸⁷Sr influx to the system from the Archean rocks exposed in the drainage area, implying that the deep-seated Archean rocks were already exposed in southern India by the late Miocene.     

The isotopic composition of particulate organic carbon in mountain rivers of Taiwan

Small rivers draining mountain islands are important in the transfer of terrestrial particulate organic carbon (POC) to the oceans. This input has implications for the geochemical stratigraphic record. We have investigated the stable isotopic composition of POC (δ¹³C_org) in rivers draining the mountains of Taiwan. In 15 rivers, the suspended load has a mean δ¹³C_org that ranges from −28.1±0.8‰ to −22.0±0.2‰ (on average 37 samples per river) over the interval of our study. To investigate this variability we have supplemented suspended load data with measurements of POC in bedrock and river bed materials, and constraints on the composition of the terrestrial biomass. Fossil POC in bedrock has a range in δ¹³C_org from −25.4±1.5‰ to −19.7±2.3‰ between the major geological formations. Using coupled δ¹³C_org and N/C we have found evidence in the suspended load for mixing of fossil POC with non-fossil POC from the biosphere. In two rivers outside the Taiwan Central Range anthropogenic land use appears to influence δ¹³C_org, resulting in more variable and lower values than elsewhere. In all other catchments, we have found that 5‰ variability in δ¹³C_org is not controlled by the variable composition of the biomass, but instead by heterogeneous fossil POC.In order to quantify the fraction of suspended load POC derived from non-fossil sources (F_nf) as well as the isotopic composition of fossil POC (δ¹³C_fossil) carried by rivers, we adapt an end-member mixing model. River suspended sediments and bed sediments indicate that mixing of fossil POC results in a negative trend between N/C and δ¹³C_org that is distinct from the addition of non-fossil POC, collapsing multiple fossil POC end-members onto a single mixing trend. As an independent test of the model, F_nf reproduces the fraction modern (F_mod) in our samples, determined from ¹⁴C measurements, to within 0.09 at the 95% confidence level. Over the sampling period, the mean F_nf of suspended load POC was low (0.29 ± 0.02, n = 459), in agreement with observations from other mountain rivers where physical erosion rates are high and fossil POC enters river channels. The mean δ¹³C_fossil in suspended POC varied between −25.2±0.5‰ and −20.2±0.6‰ from catchment to catchment. This variability is primarily controlled by the distribution of the major geological formations. It also covers entirely the range of δ¹³C_org found in marine sediments which is commonly thought to derive from mixing between marine and terrigenous POC. If land-sourced POC is preserved in marine sediments, then changes in the bulk δ¹³C_org observed offshore Taiwan could instead be explained by changes in the onshore provenance of sediment. The range in δ¹³C_org of fossil organic matter in sedimentary rocks exposed at the surface is large and given the importance of these rocks as a source of clastic sediment to the oceans, care should be taken in accounting for fossil POC in marine deposits supplied by active mountain belts.     

Stable silicon isotopes of groundwater, feldspars, and clay coatings in the Navajo Sandstone aquifer, Black Mesa, Arizona, USA

We present some of the first analyses of the stable isotopic composition of dissolved silicon (Si) in groundwater. The groundwater samples were from the Navajo Sandstone aquifer at Black Mesa, Arizona, USA, and the Si isotope composition of detrital feldspars and secondary clay coatings in the aquifer were also analyzed. Silicon isotope compositions were measured using high-resolution multi-collector inductively coupled mass spectrometry (HR-MC-ICP-MS) (Nu1700 & NuPlasma HR). The quartz dominated bulk rock and feldspar separates have similar δ³⁰Si of −0.09 ± 0.04‰ and −0.15 ± 0.04‰ (±95% SEM), respectively, and clay separates are isotopically lighter by up to 0.4‰ compared to the feldspars. From isotopic mass-balance considerations, co-existing aqueous fluids should have δ³⁰Si values heavier than the primary silicates. Positive δ³⁰Si values were only found in the shallow aquifer, where Si isotopes are most likely fractionated during the dissolution of feldspars and subsequent formation of clay minerals. However, δ³⁰Si decreases along the flow path from 0.56‰ to −1.42‰, representing the most negative dissolved Si isotope composition so far found for natural waters. We speculate that the enrichment in ²⁸Si is due to dissolution of partly secondary clay minerals and low-temperature silcretes in the Navajo Sandstone. The discovery of the large range and systematic shifts of δ³⁰Si values along a groundwater flow path illustrates the potential utility of stable Si isotopes for deciphering the Si cycling in sedimentary basins, tracing fluid flow, and evaluating global Si cycle.     

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.     

Germanium isotopic variations in igneous rocks and marine sediments

A new technique for the precise and accurate determination of Ge stable isotope compositions has been developed and applied to silicate rocks and biogenic opal. The analyses were performed using a continuous flow hydride generation system coupled to a MC-ICPMS. Samples have been purified through anion- and cation-exchange resins to separate Ge from matrix elements and eliminate potential isobaric interferences. Variations of ⁷⁴Ge/⁷⁰Ge ratios are expressed as δ⁷⁴Ge values relative to our internal standard and the long-term external reproducibility of the data is better than 0.2‰ for sample size as low as 15 ng of Ge. Data are presented for igneous and sedimentary rocks, and the overall variation is 2.4‰ in δ⁷⁴Ge, representing 12 times the uncertainty of the measurements and demonstrating that the terrestrial isotopic composition of Ge is not unique. Co-variations of ⁷⁴Ge/⁷⁰Ge, ⁷³Ge/⁷⁰Ge and ⁷²Ge/⁷⁰Ge ratios follow a mass-dependent behaviour and imply natural isotopic fractionation of Ge by physicochemical processes. The range of δ⁷⁴Ge in igneous rocks is only 0.25‰ without systematic differences among continental crust, oceanic crust or mantle material. On this basis, a Bulk Silicate Earth reservoir with a δ⁷⁴Ge of 1.3 ± 0.2‰ can be defined. In contrast, modern biogenic opal such as marine sponges and authigenic glauconite displayed higher δ⁷⁴Ge values between 2.0‰ and 3.0‰. This suggests that biogenic opal may be significantly enriched in light isotopes with respect to seawater and places a lower bound on the δ⁷⁴Ge of the seawater to +3.0‰.This suggests that seawater is isotopically heavy relative to Bulk Silicate Earth and that biogenic opal may be significantly fractionated with respect to seawater. Deep-sea sediments are within the range of the Bulk Silicate Earth while Mesozoic deep-sea cherts (opal and quartz) have δ⁷⁴Ge values ranging from 0.7‰ to 2.0‰. The variable values of the cherts cannot be explained by binary mixing between a biogenic component and a detrital component and are suggestive of enrichment in the light isotope of diagenetic quartz. Further work is now required to determine Ge isotope fractionation by siliceous organisms and to investigate the effect of diagenetic processes during chert lithification.     

Abundance of mass 47 CO2 in urban air, car exhaust, and human breath

Atmospheric carbon dioxide is widely studied using records of CO₂ mixing ratio, δ¹³C and δ¹⁸O. However, the number and variability of sources and sinks prevents these alone from uniquely defining the budget. Carbon dioxide having a mass of 47 u (principally ¹³C¹⁸O¹⁶O) provides an additional constraint. In particular, the mass 47 anomaly (Δ₄₇) can distinguish between CO₂ produced by high temperature combustion processes vs. low temperature respiratory processes. Δ₄₇ is defined as the abundance of mass 47 isotopologues in excess of that expected for a random distribution of isotopes, where random distribution means that the abundance of an isotopologue is the product of abundances of the isotopes it is composed of and is calculated based on the measured ¹³C and ¹⁸O values. In this study, we estimate the δ¹³C (vs. VPDB), δ¹⁸O (vs. VSMOW), δ47, and Δ₄₇ values of CO₂ from car exhaust and from human breath, by constructing ‘Keeling plots’ using samples that are mixtures of ambient air and CO₂ from these sources. δ47 is defined as , where is the R⁴⁷ value for a hypothetical CO₂ whose δ¹³C_VPDB = 0, δ¹⁸O_VSMOW = 0, and Δ₄₇ = 0. Ambient air in Pasadena, CA, where this study was conducted, varied in [CO₂] from 383 to 404 μmol mol⁻¹, in δ¹³C and δ¹⁸O from −9.2 to −10.2‰ and from 40.6 to 41.9‰, respectively, in δ47 from 32.5 to 33.9‰, and in Δ₄₇ from 0.73 to 0.96‰. Air sampled at varying distances from a car exhaust pipe was enriched in a combustion source having a composition, as determined by a ‘Keeling plot’ intercept, of −24.4 ± 0.2‰ for δ¹³C (similar to the δ¹³C of local gasoline), δ¹⁸O of 29.9 ± 0.4‰, δ47 of 6.6 ± 0.6‰, and Δ₄₇ of 0.41 ± 0.03‰. Both δ¹⁸O and Δ₄₇ values of the car exhaust end-member are consistent with that expected for thermodynamic equilibrium at∼200 °C between CO₂ and water generated by combustion of gasoline-air mixtures. Samples of CO₂ from human breath were found to have δ¹³C and δ¹⁸O values broadly similar to those of car exhaust-air mixtures, −22.3 ± 0.2 and 34.3 ± 0.3‰, respectively, and δ47 of 13.4 ± 0.4‰. Δ₄₇ in human breath was 0.76  ± 0.03‰, similar to that of ambient Pasadena air and higher than that of the car exhaust signature.     

Oxygen isotopic studies of the interaction between xenoliths and mafic magma, Voisey’s Bay Intrusion, Labrador, Canada

Sulfide mineralization in the Voisey’s Bay Intrusion, Labrador, Canada, is closely associated with country rock xenoliths that have extensively reacted with basaltic magma. In order to better understand the processes that control the assimilation of country rocks by mafic magma, a detailed study of oxygen isotope systematics related to magma-country rock interaction in the Voisey’s Bay area was undertaken. Protracted interaction of the xenoliths with magma produced refractory mineral assemblages in the xenoliths (2-10 cm in diameter) composed of Ca-rich plagioclase, corundum, hercynite, and minor magnetite. Overgrowth rims of plagioclase and biotite that surround most xenoliths separate the restites from the enclosing igneous matrix. The δ¹⁸O values of minerals from regionally metamorphosed pelitic and quartzofeldspathic protoliths are: plagioclase (8.7-12.3‰), orthoclase (9.5-9.8‰), biotite (5.2-8.7‰), garnet (8.3-10.8‰), pyroxene (8.0-10.1‰), and quartz (9.6-14.0). The δ¹⁸O values of minerals from the hornfels in the contact aureole of the intrusion are consistent with modeling which indicates that as a result of essentially closed system contact metamorphism oxygen isotope values should differ only slightly from those of the protoliths. Hercynite, plagioclase, and corundum separates from the xenoliths have δ¹⁸O values that vary from 2.9‰ to 10.5‰, 5.6‰ to 10.9‰, and 2.0‰ to 6.8‰, respectively. Although a siliceous ¹⁸O-enriched melt has been lost from the xenoliths, corundum, and feldspar δ¹⁸O values are significantly lower than expected through melt loss alone. The relatively low δ¹⁸O values of minerals from the xenoliths may be a function of incomplete isotopic exchange with surrounding mafic magma which had a δ¹⁸O value of ∼5.5‰ to 6.0‰. The high-¹⁸O melt that was released from the xenoliths is partially recorded in the plagioclase overgrowth on the margin of the xenoliths (δ¹⁸O values from 6.2‰ to 10.7‰), and in hercynite that replaced corundum. However, mass balance calculations indicate that a portion of the partial melt must have been transferred to magma that was moving through the conduit system. δ¹⁸O and δD values of biotite surrounding the plagioclase overgrowth range from 5.0‰ to 6.2‰ and −58‰ to −80‰, respectively. These data suggest that the outermost rim associated with many xenoliths has closely approached isotopic equilibrium with uncontaminated mafic magma. The current gabbroic to troctolitic matrix of the xenoliths shows no evidence for contamination by the high-¹⁸O partial melt from the xenoliths. The feldspar and biotite overgrowths on the xenoliths that formed after the motion of the xenoliths relative to the magma had stopped prevented further isotopic exchange between the xenoliths and final magma. The minerals within the xenoliths are not in oxygen isotopic equilibrium with each other, due in part to rapid thermal equilibration, partial melting, and partial exchange with flow through magma.     

Lithium isotope analysis of olivine by SIMS: Calibration of a matrix effect and application to magmatic phenocrysts

High-spatial resolution analysis of light element isotope variations by secondary-ion mass spectrometry (SIMS) has numerous applications in geochemistry and cosmochemistry. Recent attention has focused on ⁷Li/⁶Li variations in magmatic phenocrysts to infer the volatile degassing history of their parent magmas, and on minerals from mantle samples to determine source-region processes and the recycling history of mantle reservoirs. In these studies the effect of mineral composition on the ⁷Li/⁶Li ratio measured by SIMS has been considered secondary, and generally disregarded. We show, using a suite of nine olivines analyzed by MC-ICP-MS or TIMS, that there is a substantial effect of composition on the ⁷Li/⁶Li ratio of olivine measured by SIMS. For magnesian olivine (74 < Fo < 94) the effect is a linear function of composition, with δ⁷Li increasing by 1.3‰ for each mole percent decrease in forsterite component. At higher Fe contents, the relationship ceases to be linear. The composition range over which linear behavior is exhibited appears to depend on instrumental conditions. A calibration of this matrix effect over the linear range is presented, assuming the measurement of ⁷Li/⁶Li relative to an olivine standard of known composition. Application of this calibration to a suite of olivines separated from basaltic lavas from Ko'olau, Hawai'i demonstrates that the matrix effect is responsible for a geologically spurious correlation between δ⁷Li and Mg#. However, after correction, the olivines retain evidence of Li isotope heterogeneity, the degree and nature of which differs in each of the four separates studied. These results emphasize the importance of compositional correction for SIMS measurement of δ⁷Li in olivine, particularly in zoned crystals, and support previous conclusions that Li isotope variability in igneous materials is subject to late-stage disturbance. The significant matrix effect demonstrated for olivine suggests that matrix effects in other minerals require further evaluation.     

Note on the isotopic geochemistry of fossil-lacustrine tufas in carbonate plateau—A study from Dungul region (SW Egypt)

Lithological, chemical, and stable isotope data are used to characterize lacustrine tufas dating back to pre-late Miocene and later unknown times, capping different surfaces of a Tertiary carbonate (Sinn el-Kedab) plateau in Dungul region in the currently hyperarid south-western Egypt. These deposits are composed mostly of calcium carbonate, some magnesium carbonate and clastic particles plus minor amounts of organic matter. They have a wide range of (Mg/Ca)_molar ratios, from 0.03 to 0.3. The bulk-tufa carbonate has characteristic isotope compositions: (δ¹³C_mean = −2.49 ± 0.99‰; δ¹⁸O_mean = −9.43 ± 1.40‰). The δ¹³C values are consistent with a small input from C4 vegetation or thinner soils in the recharge area of the tufa-depositing systems. The δ¹⁸O values are typical of fresh water carbonates. Covariation between δ¹³C and δ¹⁸O values probably is a reflection of climatic conditions such as aridity. The tufas studied are isotopically similar to the underlying diagenetic marine chalks, marls and limestones (δ¹³C_mean = −2.06 ± 0.84‰; δ¹⁸O_mean = −10.06 ± 1.39‰). The similarity has been attributed to common meteoric water signatures. This raises large uncertainties in using tufas (Mg/Ca)_molar, δ¹³C and δ¹⁸O records as proxies of paleoclimatic change and suggests that intrinsic compositional differences in material sources within the plateau may mask climatic changes in the records.