Source of coesite inclusions within inherited magmatic zircon from Sulu UHP rocks, eastern China, and their bearing for fluid–rock interaction and SHRIMP dating

The 5-km deep Chinese Continental Scientific Drilling Main Hole penetrated a sequence of ultrahigh pressure (UHP)-metamorphic rocks consisting mainly of eclogite, gneiss and garnet-peridotite with minor schist and quartzite. Zircon separates taken from thin layers of schist and gneiss within eclogite were investigated. Cathodoluminescence images of zircon grains show that they have oscillatory zoned magmatic cores and unzoned to patchy zoned metamorphic rims. Zircon rims contain rare coesite and calcite inclusions whereas cores contain inclusions of both low- P minerals (e.g. feldspar, biotite and quartz) and coesite and other eclogite-facies minerals such as phengite and jadeite. The zircon cores give highly variable ²⁰⁶Pb/²³⁸U ages ranging from 760 to 431 Ma for schist and from 698 to 285 Ma for gneiss, and relatively high but variable Th/U ratios (0.16–1.91). We suggest that the coesite and other eclogite facies mineral inclusions in zircon cores were not magmatic but formed through metasomatic processes caused by fluids during UHP metamorphism, and that the fluids contain components of SiO₂, Al₂O₃, K₂O, FeO, MgO, Na₂O and H₂O. Metasomatism of the Sulu UHP rocks during continental subduction to mantle depths has partly altered magmatic zircon cores and reset isotopic systems. This study provides key evidence that mineral inclusions within magmatic zircon domains are not unequivocal indicators of the formation conditions of the respective domain. This finding leads us to conclude that the routine procedure for dating of metamorphic events solely based on the occurrence of mineral inclusions in zoned zircon could be misleading and the data should be treated with caution.     

Oxygen isotope exchange and closure temperatures in cooling rocks

Retrograde exchange of oxygen isotopes between minerals in igneous and metamorphic rocks by means of diffusion is explored using a finite difference computer model, which predicts both the zonation profile of δ¹⁸O within grains, and the bulk δ¹⁸O value of each mineral in the rock. Apparent oxygen isotope equilibrium temperatures that would be observed in these rocks are calculated from the δ¹⁸O values of each mineral pair within the rock. In systems which cool linearly from a sufficiently high temperature or at a low enough cooling rate, such that the final oxygen isotope values are not dependent upon the initial oxygen isotope values ('slow cooling'), the apparent oxygen isotope temperature derived for a rock composed of a single mineral pair can be shown to be simply related to the Dodson closure temperatures ( T _c) for the two phases and the mode of the rock. Adding a third phase into a system which undergoes 'slow' cooling will cause the apparent temperature derived for the two minerals already present to differ from the simple relationship for a two-phase system. In some systems oxygen isotope reversals can be developed. If cooling is not 'slow', then the mineral δ¹⁸O values resulting from cooling will be partly dependent upon the initial temperature of the system concerned. The model successfully simulates the mineral δ¹⁸O values that are often observed in granitic rocks. Application of the model will help in assessing the validity of oxygen isotope thermometry in different geological settings, and allows quantitative prediction of the oxygen isotope fractionations that are developed in cooling closed systems.     

Pressure–temperature conditions and retrograde paths of eclogites, garnet–glaucophane rocks and schists from South Sulawesi, Indonesia

High-pressure metamorphic rocks exposed in the Bantimala area, c . 40  km north-east of Ujung Pandang, were formed as a Cretaceous subduction complex with fault-bounded slices of melange, chert, basalt, turbidite, shallow-marine sedimentary rocks and ultrabasic rocks. Eclogites, garnet–glaucophane rocks and schists of the Bantimala complex have estimated peak temperatures of T  =580–630 °C at 18  kbar and T  =590–640 °C at 24  kbar, using the garnet–clinopyroxene geothermometer. The garnet–omphacite–phengite equilibrium is used to estimate pressures. The distribution coefficient K _D1=[( X _pyr)³( X _grs)⁶/( X _di)⁶]/[(Al/Mg)_M2,wm (Al/Si)_T2,wm]³ among omphacite, garnet and phengite is a good index for metamorphic pressures. The K _D1values of the Bantimala eclogites were compared with those of eclogites with reliable P–T  estimates. This comparison suggests that peak pressures of the Bantimala eclogites were P =18–24  kbar at T  =580–640 °C. These results are consistent with the P–T  range calculated using garnet–rutile–epidote–quartz and lawsonite–omphacite–glaucophane–epidote equilibria.     

Controls on the first appearance of jadeitic pyroxene, northern Diablo Range, California

A reaction producing jadeitic pyroxene in metagreywackes of the northern Diablo Range has been identified on the basis of mineral distribution, isograd patterns and composition of coexisting minerals. The appearance of jadeitic pyroxene (∼Jd₈₀) is closely followed by the disappearance of pumpellyite, which indicates that pumpellyite plays a major role in the pyroxene-producing reaction. A new projection from hematite, lawsonite, chlorite, quartz and H₂O on to the NaAlO₂-FeO-MgO ternary confirms the role of pumpellyite in pyroxene production and suggests a reaction of the form: 1.00 pumpellyite + 0.31 chlorite + 8.71 albite + 0.70 hematite + 2.00 H₂O = 8.54 jadeite + 0.57 glaucophane + 3.09 lawsonite + 5.26 quartz. Metagreywackes of the northern Diablo Range were metamorphosed under conditions of P _H2O= P _total at 200-300 °C and 7.5-10.0 kbar. Despite the low temperatures attained during metamorphism, the assumption of equilibrium yields results consistent with field observations and phase relations.     

Fluid-mineral equilibria in prehnite-pumpellyite to greenschist facies metabasites near Flin Flon, Manitoba, Canada: implications for petrogenetic grids

A sequence of regional metamorphic isograds indicating a range from prehnite-pumpellyite to lower amphibolite facies was mapped in metabasites near Flin Flon, Manitoba. The lowest grade rocks contain prehnite + pumpellyite and are cut by younger brittle faults containing epidote + chlorite + calcite. Isobaric temperature- X _CO2 and pressure-temperature (constant X _CO2) diagrams were calculated to quantify the effects of CO₂ in the metamorphic fluid on the stability of prehnite-pumpellyite facies minerals in metabasites containing excess quartz and chlorite. Prehnite and, to a lesser extent, pumpellyite are stable only in fluids with X _co2 <0.002. For X _co2>0.002, epidote + chlorite + calcite assemblages are stable. Our calculated phase relations are consistent with regional metamorphism in the Flin Flon area in the presence of an H₂O-rich fluid and a more CO₂-rich fluid in the later fault zones. We believe that the potential effects of small amounts of CO₂ in the metamorphic fluid should be assessed when considering the pressure-temperature implications of mineral assemblages in low-grade metabasites.     

Shallow-level metamorphic fluids in a high uplift rat metamorphic belt; Alpine Schist, New Zealand

Abstract Fluids, some of which are CO₂-rich (up to 40 mol.% CO₂) and some of which are highly saline (up to 18 wt% NaCl equivalent), are trapped as fluid inclusions in quartz-calcite (∼ metallic minerals) veins which cross-cut the pumpellyite-actinolite to amphibolite facies rocks of the Alpine Schist. Fluids were commonly trapped as immiscible liquid-vapour mixes in quartz and calcite showing open-space growth textures. Fluid entrapment occurred at fluid pressures near 500 bars (possibly as low as 150 bars) at temperatures ranging from 260 to 330° C. Saline fluids may have formed by partitioning of dissolved salts into an aqueous phase on segregation of immiscible fluids from a low-density CO₂-rich fluid. Calcite deposited by these fluids has δ¹³C ranging from – 8.4 to – 11.5 and δ¹⁸O from + 4 to + 13. Isotopic data, fluid compositions and mode of occurrence suggest that the fluids are derived from high-grade metamorphic rocks. Fluid interaction with wall-rock has caused biotite crystallization and/or recrystallization in some rocks and retrogression of biotite to chlorite in other rocks.
Fluid penetration through the rock is almost pervasive in many areas where permeability, probably related to Alpine Fault activity, has focussed fluids on a regional scale into fractured rocks. The fluid flow process is made possible by high uplift-rates (in excess of 10 mm/year) bringing hot rocks near to the surface.     

Stable isotopic and fluid inclusion evidence for meteoric fluid penetration into an active mountain belt; Alpine Schist, New Zealand

Calcite and quartz veins have formed, and are forming, in steeply dipping fissures in the actively rising Alpine Schist metamorphic belt of New Zealand. The fluids that deposited these minerals were mostly under hydrostatic pressure almost down to the brittle-ductile transition, which has been raised to 5-6 km depth by rapid uplift. Some fluids were trapped under lithostatic pressures. Fluids in the fissure veins were immiscible H₂O + NaCl-CO₂ mixtures at 200-350 C. Bulk fluid composition is 15-20 mol% CO₂ and <4.3 total mol CH₄+ N₂+ Ar/100mol H₂O. Water hydrogen isotopic ratio δD_H2O in the fissure veins spans -29 to -68‰, δ¹⁸O_H2O -0.7 to 8.5‰, and bulk carbon isotopic ratio δ¹³C ranges from -3.7 to -11.7‰. The oxygen and hydrogen isotopic data suggest that the water has a predominantly meteoric source, and has undergone an oxygen isotope shift as a result of interaction with the host metamorphic rock. Similar fluids were present during cooling and uplift. Dissolved carbon is not wholly derived from residual metamorphic fluids; part may be generated by oxidation of graphite.     

Genesis of Kanggur Gold Deposit in Eastern Tianshan Orogenic Belt, NW China: Fluid Inclusion and Oxygen and Hydrogen Isotope Constraints

Abstract. Primary fluid inclusions in quartz and carbonates from the Kanggur gold deposit are dominated by aqueous inclusions, with subsidiary CO₂-H₂O inclusions that have a constant range in CO₂ content (10–20 vol %). Microthermometric results indicate that total homogenization temperatures have a wide but similar range for both aqueous inclusions (120 to 310C) and CO₂-H₂O inclusions (140 to 340C). Estimates of fluid salinity for CO₂-H₂O inclusions are quite restricted (5.9∼10.3 equiv. wt% NaCl), whereas aqueous inclusions show much wider salinity ranging from 2.2 to 15.6 equivalent wt %NaCl.
The 6D values of fluid inclusions in carbonates vary from -45 to -61 %, in well accord with the published δD values of fluid inclusions in quartz (-46 to -66 %). Most of the δ¹⁸O and δD values of the ore-forming fluids can be achieved by exchanged meteoric water after isotopic equilibration with wall rock by fluid/rock interaction at a low water/rock ratio. However, the exchanged meteoric water alone cannot explain the full range of δ¹⁸O and δD values, magmatic and/or meta-morphic water should also be involved. The wide salinity in aqueous inclusions may also result from mixing of meteoric water and magmatic and/or metamorphic water.     

The (Na–Ca)amphibole–albite–chlorite–epidote–quartz geothermobarometer in the system S–A–F–M–C–N–H2O. 1. An empirical calibration

The proposed geothermobarometer is based on an empirical calibration which takes account of two equilibria involving the tremolite, edenite, pargasite and hastingsite components in amphiboles. It has applications to assemblages found in metabasic rocks of widely different chemical compositions (magnesian to Fe-rich metabasalts), and for metamorphism ranging from lower greenschist to highest amphibolite facies. Knowing the Si_(T1), Al^iv, Al^vi, Fe³⁺, Fe²⁺, Mg, Ca, Na_M4, Na_A and A vacancy in an amphibole, and the Al³⁺ and X _Mg in coexisting epidote and chlorite, it is possible to calculate two values of In K _d for this assemblage. These equilibria involve edenite-tremolite and (pargasite/hastingsite)-tremolite end-members in amphibole (the calculation program is given). For these equilibria, the isopleths (iso-values of K _d) have been calculated for 0.27 < X _Mg < 0.75 and 0 < X _Fe3+= Fe³⁺/(Fe³⁺+ Al^vi) < 0.8. It is then possible to determine pressure and temperature directly when X _Mg, X _Fe3+, In K _d for tremoliteedenite and In K _d for (pargasite/hastingsite)-tremolite are known. Application of this geothermobarometer is limited to Ca-free plagioclase assemblages, and complete P–T paths can be drawn only if all the minerals are considered together. Phase relations at successive stages of crystallization can be constrained by studying the relationships between the coexisting minerals, their zoning and the metamorphic fabrics.     

High Precision Determination of 48Ca/42Ca Ratio by TIMS for Ca Isotope Fractionation Studies

For studies of mass-dependent fractionation of calcium isotopes in natural materials, the ⁴⁸Ca/⁴²Ca ratio is a superior choice to the conventionally measured ⁴⁴Ca/⁴⁰Ca ratio for two important reasons. These are (1) mass-dependent fractionation can be determined free from the effects of inherited or ingrown radiogenic ⁴⁰Ca and (2) this ratio increases the spread of measured isotopic masses by 50%, resulting in statistically better resolution of fractionation, assuming similar precision. A third, though strictly technical, advantage is the inherent ability of a mass spectrometer to measure ratios close to unity (⁴⁸Ca/⁴²Ca) more precisely than very small or large ratios (⁴⁴Ca/⁴⁰Ca). However, because of the very low natural abundance of both ⁴⁸Ca and ⁴²Ca, their ratio has been very difficult to measure, the only attempt so far being on a high mass resolution MC-ICP-MS with a precision of 0.33%. We report here determination of the ⁴⁸Ca/⁴²Ca ratio by the more commonly available and user-friendly multi-collector TIMS using a ⁴³Ca-⁴⁶Ca double-spike, with a significantly better precision of 0.18% (2s). The ⁴⁸Ca/⁴⁰Ca or ⁴⁴Ca/⁴⁰Ca ratio can also be measured in the same mass spectrometer run to provide complementary information on any radiogenic component.     

Stable isotope study of carbonate sediments from the Coorong Area, South Australia

ABSTRACT
The mineralogy and isotope geochemistry of carbonate minerals in the Coorong area are determined by the water chemistry of different depositional environments ranging from seawater to evaporitically modified continental water. The different isotopic compositions of coexisting calcite and dolomite suggest that each of the above two minerals was formed from water of composition and origin unique to that specific mineral. In addition, the dolomite was not formed by simple solid state cation exchange.
The occurrence of two types of dolomite was shown by isotope analysis and SEM observations. The dolomite, which is isotopically light (δ¹³C = -1 to -2% 0 ; δ¹⁸O=+3 to +5%₀) and of fine grain size (˜ 0·5 μm) probably precipitated under the influence of evaporitically modified continental water. Coarser grained dolomite (up to 4 μm) is isotopically heavier (δ¹³C=+3 to +4%₀; δ¹⁸O=+5 to + 6%₀) contains Mg in excess of Ca and was formed in or close to equilibrium with atmospheric CO² probably by the dolomitization of aragonite.     

Mass Transfer, Oxygen Isotopic Variation and Gold Precipitation in Epithermal System: A Case Study of the Hishikari Deposit, Southern Kyushu, Japan

Abstract: Transportation of various kinds of elements occurred in wall rocks (Quaternary andesites) during the hydrothermal alteration accompanied by the Hishikari epithermal gold mineralization. For example, K₂O and MgO contents of wall rocks decrease away from the gold-quartz veins, while (CaO+Na₂O) content increases, and SiO₂ content is variable near the veins. Hydrothermal alteration zoning and bulk compositional variations in wall rocks suggest that the mixing of hydrothermal solution and acidic groundwater took place an important role as the cause for the hydrothermal alteration and bulk compositional variations. The relationship between dissolved silica concentration and temperature of hydrothermal solution mixed with groundwater is obtained based on precipitation kinetics-fluid flow–mixing model, and the computed results are compared with the distribution of SiO₂ minerals (quartz and cristobalite) in the hydrothermal alteration zones. This comparison suggests that the most reasonable flow rate of fluids migrating through hydrothermal alteration zones, and A/M (A: surface area of rocks interacting with fluid, M: mass of fluid) are estimated to be ca. 10^-4.2 m/sec, and ca. 0.10 m²/kg, respectively. The mixing of two fluids (hydrothermal solution and acidic groundwater) can also explain δ¹⁸O zoning in the altered country rocks, hydrothermal alteration zoning from K-feldspar through K-mica to kaolinite from the center (veins) to margin, and deposition of gold.     

The evolution of fluorine-rich felsic magmas: source dichotomy, magmatic convergence and the origins of topaz granite

Topaz granite is alkali-feldspar granite that contains essential albite, quartz, K-feldspar, lithium-mica, and topaz. As a group topaz granites are characterized by their extreme enrichment in F (up to 3 wt%) and a wide variety of lithophile elements. They can be subdivided into a 'low-P₂O₅ subtype' (P₂O₅ < 0.1 wt%, Al₂O₃ < 14.5 wt%, SiO₂ > 73 wt%) and a 'high-P₂O₅ subtype' (P₂O₅ > 0.4 wt%, Al₂O₃ > 14.5 wt%, SiO₂ < 73 wt%), the δ¹⁸O values of which indicate a dichotomy of source rock: the low-P₂O₅ subtype (δ¹⁸O < 10‰) having a meta-igneous protolith and the high-P₂O₅ subtype (δ¹⁸O > 10 ‰) a source with a significant component of pelitic material. The unusually high F contents enhance the efficacy of melt segregation and crystal-melt fractionation and so facilitate extreme differentiation in topaz granite magmas. Very low melt volumes restrict the bulk composition of the partial melts regardless of the nature of the source; and extreme fractionation forces them along a path of magmatic convergence, to produce a group of granitic rocks with near-minimum compositions so enriched in a variety of lithophile elements (Li, Nb, Ta, Sn) that economic mineralization often results.     

Hydrogen, Oxygen and Sulfur Isotope Studies of Seafloor Hydrothermal System at the Desmos Caldera, Manus Back-arc Basin, Papua New Guinea: An Analogue of Terrestrial Acid Hot Crater-lake

Abstract. The Onsen site is an active submarine hydrothermal system hosted by the Desmos caldera in the Eastern Manus Basin, Papua New Guinea. The hydrothermal fluid is very acidic (pH=1.5) and abundant native sulfur is deposited around the vent. The δ³⁴S values of native sulfur range from -6.5 to -9.3 %o. δ³⁴S values of H₂S and SO₄ in the hydrothermal fluid are -4.3 to -9.9 %o and +18.6 to +20.0 %o, respectively. These δ³⁴S values are significantly lower than those of the other hydrothermal systems so far reported. These low δ³⁴S values and the acidic nature of the vent fluids suggest that volcanic SO₂ gas plays an important role on the sulfur isotope systematic of the Onsen hydrothermal system. Relationship among the δ³⁴S values of S-bearing species can be successively explained by the model based on the disproportionation reaction starting from the volcanic SO₂ gas. The predicted δ³⁴S values of SO₂ agree with the measured whole rock δ³⁴S values. δD and δ¹⁸O values of clay minerals separated from the altered rock samples also suggest the contribution of the magmatic fluid to the hydrothermal system. Present stable isotopic study strongly suggests that the Onsen hydrothermal site in the Desmos caldera is a magmatic submarine hydrothermal system.     

Using estimated thermodynamic properties to model accessory phases: the case of tourmaline

Accessory phases and minor components in minerals are commonly ignored in thermodynamic modelling. Such an approach seems unwarranted, as accessory phases can represent a significant element reservoir and minor components can substantially change their host mineral's stability field. However, a lack of thermodynamic data prohibits assessment of these effects. In this contribution, the polyhedron method is used to estimate the thermodynamic properties of tourmaline, a common and widespread accessory phase, stable over a range of P–T–X conditions. The polyhedron method allows Δ H , S , V , C _P and V _m ( T , P ) properties to be estimated from a linear stoichiometric summation over the fractional properties of its polyhedron constituents. To allow for estimates of tourmaline, fractional thermodynamic properties for B^III and B^IV polyhedra were derived. Mixing contributions to molar volume were evaluated and symmetrical mixing parameters derived for Al-Mg, Al-Fe and Al-Li interaction on tourmaline's Y-site and T-site Al-Si interaction. Evaluation of the estimated properties using experimental and natural equilibria between tourmaline and melts, minerals and hydrothermal fluids, shows that reliable semi-quantitative results are obtained. The boron contents in fluids coexisting with tourmaline are calculated to within an order of magnitude of measured content, and where anchor-points are available, agreement improves to within a factor of 2. Including tourmaline in petrogenetic modelling of metamorphic rocks indicates that its presence leads to disappearance of staurolite and garnet, among others, and modifies the X _Mg of coexisting phases, in line with observations on natural rocks.     

Refinements in BaSO4 to CO2 Preparation and δ18 O Calibration of the Sulfate Reference Materials NBS-127, IAEA SO-5 and IAEA SO-6

Refinements have been made to achieve over 99% yield in the conversion of CO to CO₂ in order to improve the reproducibility and accuracy of δ¹⁸ O measurements in sulfates. BaSO₄ (10-15 mg) was mixed with an identical amount of spectrographic-grade graphite and loaded into a Pt boat. The mixture was gradually heated to 1100 °C to reduce sulfate to CO and CO₂; the former gas was simultaneously converted to CO₂ by a glow discharge between Pt electrodes immersed in a magnetic field (produced by a pair of external neodymium magnets). A small memory effect was noticed during the analysis (less than 0.3‰ per 10‰ difference in δ¹⁸ O between two subsequently analysed samples). The memory effect, however, was suppressed by repetitive preparation of the same specimen. CO₂ produced in this way from sulfate reference samples was analysed on a dual inlet and triple collector mass spectrometer along with CO₂ equilibrated with VSMOW, GISP and SLAP water reference samples. To avoid large departures of measured isotope ratios from ¹⁸O/¹⁶O of the working calibrator we used CO₂ gas prepared from ocean water sulfate for this purpose. The calibrated δ¹⁸ O values (in ‰) obtained in this way for NBS-127, IAEA SO-5 and IAEA SO-6 reference materials were 8.73 ± 0.05, 12.20 ± 0.07 and -10.43 ± 0.12, respectively.     

P-T grids for silica-undersaturated granulites in the systems MAS (n+ 4) and FMAS (n+ 3)-tools for the derivation of P-T paths of metamorphism

Abstract Considering the minerals cordierite (Cd), sapphirine (Sa), hypersthene (Hy), garnet (Ga), spinel (Sp), sillimanite (Si) and corundum (Co) in the system FeO-MgO-Al₂O₃-SiO₂ (FMAS), the stable invariant points are [Co], [Ga], [Cd] and [Sa]. Constraints imposed by experimental data for the system MAS indicate that under low P _H2o conditions the invariant points occur at high temperature (> 900° C) and intermediate pressure (7-10 kbar). This temperature is higher than that commonly advocated for granulite facies metamorphism. In granulites Fe-Mg exchange geothermometers may yield temperatures of 100–150° C below peak metamorphic conditions and evidence for peak temperatures is best preserved by relict high-temperature assemblages and by Al-rich cores in orthopyroxene. Application of the FMAS grid to some well-documented granulite occurrences introduces important constraints on their P-T histories. Rocks of different bulk compositions, occurring in close proximity in the field, may record distinct segments of their P-T paths. This applies particularly to rocks with evidence for reaction in the form of coronas, symplectites and zoned minerals. Consideration of curved reaction boundaries and _XMs isopleths may explain apparently contradictory results for the stability of cordierite obtained from low-temperature experiments and thermochemical calculations on the one hand and hightemperature experimental data on the other.     

Reaction textures in calc-silicate granulites from the Bolingen Islands, Prydz Bay, East Antarctica: implications for the retrograde P–T path

Calc-silicate granulites from the Bolingen Islands, Prydz Bay, East Antarctica, exhibit a sequence of reaction textures that have been used to elucidate their retrograde P–T path. The highest temperature recorded in the calc-silicates is represented by the wollastonite- and scapolite-bearing assemblages which yield at least 760°C at 6 kbar based on experimental results. The calc-silicates have partially re-equilibrated at lower temperatures (down to 450°C) as evidenced by the successive reactions: (1) wollastonite + scapolite + calcite = garnet + CO₂, (2) wollastonite + CO₂= calcite + quartz, (3) wollastonite + plagioclase = garnet + quartz, (4) scapolite = plagioclase + calcite + quartz, (5) garnet + CO₂+ H₂O = epidote + calcite + quartz, and (6) clinopyroxene + CO₂+ H₂O = tremolite + calcite + quartz.
The reaction sequence observed indicates that a CO₂ was relatively low in the wollastonite-bearing rocks during peak metamorphic conditions, and may have been further lowered by local infiltration of H₂O from the surrounding migmatitic gneisses on cooling. Fluid activities in the Bolingen calc-silicates were probably locally variable during the granulite facies metamorphism, and large-scale CO₂ advection did not occur.
A retrograde P–T path, from the sillimanite stability field ( c. 760°C at 6 kbar) into the andalusite stability field ( c. 450°C at <3 kbar), is suggested by the occurrence of secondary andalusite in an adjacent cordierite–sillimanite gneiss in which sillimanite occurs as inclusions in cordierite.     

Determination of the Absolute Hydrogen Isotopic Ratio of V-SMOW and SLAP

By mixing ¹H₂O and ²H₂O, both with accurately known purity, samples were prepared with ²H/¹H ratios close to those of the international isotopic water standards: V-SMOW and SLAP. A mass spectrometrical comparison of these calibration samples with the actual water standards revealed:
²H/¹H of V-SMOW = (155.95 ± 0.08) × 10^−6
2H/¹H of SLAP = (89.12 ± 0.07) × 10⁻⁶
δ²H_V-SMOW(SLAP) =−428.5 ± 0.5 %