Li-Sr-Nd isotope signatures of the plume and cratonic lithospheric mantle beneath the margin of the rifted Tanzanian craton (Labait)

Lithium concentrations and isotopic compositions of olivine and ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd of coexisting clinopyroxene from peridotite xenoliths from the Quaternary Labait volcano, Tanzania, document the influence of rift-related metasomatism on the ancient cratonic mantle. Olivines show negative correlations between Fo content and both δ⁷Li and Li concentrations. Olivines in iron-rich peridotites (Fo_85–87) have high Li concentrations (3.2–4.8 ppm) and heavy δ⁷Li (+5.2 to +6.6). In contrast, olivines in ancient, refractory peridotites have lower Li concentrations (∼2 ppm) and relatively light δ⁷Li (+2.6 to +3.5). This reflects mixing between ancient, refractory cratonic lithosphere and asthenosphere-derived rift magmas. A uniquely fertile, deformed, high-temperature garnet lherzolite, interpreted to be from the base of the lithosphere, has a ⁸⁷Sr/⁸⁶Sr of 0.7029 and ¹⁴³Nd/¹⁴⁴Nd of 0.51286, similar to HIMU oceanic basalts. It provides the best estimate of the Sr–Nd isotope composition of the upwelling mantle (i.e., plume, sensu lato) underlying this portion of the East African Rift, and is slightly less radiogenic compared to previous estimates of the plume that were based on rift basalts. Although elevated δ⁷Li are not exclusive to HIMU source regions, the data collectively indicate that the plume beneath Labait has HIMU characteristics in Sr, Nd and Li isotope composition. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Influence of temperature,composition, silica activity and oxygen fugacity on the H<Subscript>2</Subscript>O storage capacity of olivine at 8 GPa

Olivine crystals were grown in the presence of a hydrous silicate fluid during multi-anvil experiments at 8 GPa and 1,000–1,600°C. Experiments were conducted both in a simple system (FeO–MgO–SiO₂–H₂O) and in a more complex system containing additional elements (CaO–Na₂O–Al₂O₃–Cr₂O₃–TiO₂–FeO–MgO–SiO₂–H₂O). Silica activity was buffered by the presence of either pyroxene (high a _SiO2) or ferropericlase (low a _SiO2), and was buffered by the presence of Ni + NiO or Fe + FeO, or constrained by the presence of Fe₂O₃. Raman spectroscopy was used to identify pyroxene polymorphs in the run products. Clinoenstatite was present in the 1,000°C experiment, and enstatite in experiments at 1,400–1,520°C. The H₂O content of olivine was measured using secondary ion mass spectroscopy, and infrared spectroscopy was used to investigate the nature of hydrous defects. The H₂O storage capacity of olivine decreases with increasing temperature at 8 GPa. In contrast to previous experimental results at ≤2 GPa, no significant effect of varying oxygen fugacity is evident, but H₂O storage capacity is enhanced under conditions of low silica activity. No significant growth of low wavenumber (<3,400 cm⁻¹) peaks, generally associated with high at low pressure, was observed in the FTIR spectra of olivine from the high experiments. Our experiments show that previous high pressure H₂O storage capacity measurements for olivine synthesized under more oxidizing conditions than the Earth’s mantle are not likely to be compromised by the of the experiments. However, the considerable effect of temperature on H₂O storage capacity in olivine must be taken into account to avoid overestimation of the bulk upper mantle H₂O storage capacity.     

Chemical and boron-isotope variations in tourmalines from an S-type granite and its source rocks: the Erongo granite and tourmalinites in the Damara Belt,Namibia

Tourmaline is widespread in metapelites and pegmatites from the Neoproterozoic Damara Belt, which form the basement and potential source rocks of the Cretaceous Erongo granite. This study traces the B-isotope variations in tourmalines from the basement, from the Erongo granite and from its hydrothermal stage. Tourmalines from the basement are alkali-deficient schorl-dravites, with B-isotope ratios typical for continental crust (δ¹¹B average −8.4‰ ± 1.4, n = 11; one sample at −13‰, n = 2). Virtually all tourmaline in the Erongo granite occurs in distinctive tourmaline-quartz orbicules. This “main-stage” tourmaline is alkali-deficient schorl (20–30% X-site vacancy, Fe/(Fe + Mg) 0.8–1), with uniform B-isotope compositions (δ¹¹B −8.7‰ ± 1.5, n = 49) that are indistinguishable from the basement average, suggesting that boron was derived from anatexis of the local basement rocks with no significant shift in isotopic composition. Secondary, hydrothermal tourmaline in the granite has a bimodal B-isotope distribution with one peak at about −9‰, like the main-stage tourmaline, and a second at −2‰. We propose that the tourmaline-rich orbicules formed late in the crystallization history from an immiscible Na–B–Fe-rich hydrous melt. The massive precipitation of orbicular tourmaline nearly exhausted the melt in boron and the shift of δ¹¹B to −2‰ in secondary tourmaline can be explained by Rayleigh fractionation after about 90% B-depletion in the residual fluid. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Low-pressure,water-assisted anatexis of basic dykes in a contact metamorphic aureole,Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite–gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative δ¹⁸O whole-rock values (+0.2 to −3.4‰), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest δ¹⁸O values (−2.8, −4.2, −4.4 and −7.6‰, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic δ¹⁸O values up to the anatectic zone, where the values are lower and more heterogeneous. Low δ¹⁸O values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Oxygen isotopic compositions of Allende Type C CAIs: Evidence for isotopic exchange during nebular melting and asteroidal metamorphism

In situ oxygen isotopic measurements of primary and secondary minerals in Type C CAIs from the Allende CV3 chondrite reveal that the pattern of relative enrichments and depletions of ¹⁶O in the primary minerals within each individual CAI are similar to the patterns observed in Types A and B CAIs from the same meteorite. Spinel is consistently the most ¹⁶O-rich (Δ¹⁷O = −25‰ to −15‰), followed by Al,Ti-dioside (Δ¹⁷O = −20‰ to −5‰) and anorthite (Δ¹⁷O = −15‰ to 0‰). Melilite is the most ¹⁶O-depleted primary mineral (Δ¹⁷O = −5‰ to −3‰). We conclude that the original melting event that formed Type C CAIs occurred in a ¹⁶O-rich (Δ¹⁷O  −20‰) nebular gas and they subsequently experienced oxygen isotopic exchange in a ¹⁶O-poor reservoir. At least three of these (ABC, TS26F1 and 93) experienced remelting at the time and place where chondrules were forming, trapping and partially assimilating ¹⁶O-poor chondrule fragments. The observation that the pyroxene is ¹⁶O-rich relative to the feldspar, even though the feldspar preceded it in the igneous crystallization sequence, disproves the class of CAI isotopic exchange models in which partial melting of a ¹⁶O-rich solid in a ¹⁶O-poor gas is followed by slow crystallization in that gas. For the typical (not associated with chondrule materials) Type C CAIs as well for as the Types A and B CAIs, the exchange that produced internal isotopic heterogeneity within each CAI must have occurred largely in the solid state. The secondary phases grossular, monticellite and forsterite commonly have similar oxygen isotopic compositions to the melilite and anorthite they replace, but in one case (CAI 160) grossular is ¹⁶O-enriched (Δ¹⁷O = −10‰ to −6‰) relative to melilite (Δ¹⁷O = −5‰ to −3‰), meaning that the melilite and anorthite must have exchanged its oxygen subsequent to secondary alteration. This isotopic exchange in melilite and anorthite likely occurred on the CV parent asteroid, possibly during fluid-assisted thermal metamorphism.     

Metallographic cooling rates and origin of IVA iron meteorites

We have determined the metallographic cooling rates for 13 IVA irons using the most recent and most accurate metallographic cooling rate model. Group IVA irons have cooling rates that vary from 6600 °C/Myr at the low-Ni end of the group to 100 °C/Myr at the high-Ni end of the group. This large cooling rate range is totally incompatible with cooling in a mantled core which should have a uniform cooling rate. Thermal and fractional crystallization models have been used to describe the cooling and solidification of the IVA asteroid. The thermal model indicates that a metallic body of 150 ± 50 km in radius with less than 1 km of silicate on the outside of the body has a range of cooling rates that match the metallographic cooling rates in IVA irons in the temperature range 700-400 °C where the Widmanstätten pattern formed. The fractional crystallization model for Ni with initial S contents between 3 and 9 wt% is consistent with the measured variation of cooling rate with bulk Ni and the thermal model. New models for impacts in the early solar system and the evolution of the primordial asteroid belt allow us to propose that the IVA irons crystallized and cooled in a metallic body that was derived from a differentiated protoplanet during a grazing impact. Other large magmatic iron groups, IIAB, IIIAB, and IVB, also show significant cooling rate ranges and are very likely to share a similar history.