Contributions to Mineralogy and Petrology最新文献

The unexpected discovery of felsic magma by the Iceland Deep Drilling Project-1 (IDDP-1) in the Krafla volcanic system (KVS) presents a unique opportunity to investigate pre-eruptive lithium (Li) dynamics and establish a more direct connection between magma reservoirs and volcanic deposits. Our study provides new insights into Li abundances and isotope compositions in bulk-rock, minerals, and groundmass glass from rhyolitic lavas at KVS, encompassing various stages of groundmass crystallisation. Additionally, we examined felsic cuttings retrieved from the IDDP-1 well, comprising crystal-poor obsidian and crystal-bearing to -rich ‘felsite’ particles. Groundmass glasses from surface lavas show limited variability in K/Na, indicating limited secondary hydration of the glasses and that their Li contents seem to not be affected by this post-eruptive process. Lithium inventories in groundmass glasses and minerals within lavas exhibit variations consistent with the cooling history of the deposit, resembling patterns seen in Snake River Plain ignimbrites. Lithium contents of glassy rhyolitic lavas, whether bulk-rock (avg. 27.2 ± 3.1 μg/g) or groundmass glass (average 28.4 ± 4.7 μg/g), and their bulk isotopic compositions (avg. δ⁷Li =+ 4.4 ± 0.2‰) overlap with those observed in IDDP-1 obsidian cuts (avg. 24.9 μg/g Li in bulk, 28.6 ± 1.5 μg/g in groundmass glass, and δ⁷Li = 4.5 ± 0.2‰). Glassy lavas lacking spherulites may potentially preserve pristine magmatic Li element and isotope compositions, while areas with extensive groundmass crystallisation reveal Li enrichments in phenocrysts. Plagioclases in slowly cooled parts of the deposit record a two-fold increase in Li contents compared to plagioclase found in glassy counterparts, along with evidence of open-system degassing marked by heavier bulk Li isotope compositions and lower bulk Li contents of the crystallised lava portions (avg. δ⁷Li = +7.2 ± 0.1‰ and 7 ± 0.8 μg/g Li) relative to bulk glassy lithologies (avg. δ⁷Li = +4.1 ± 0.1‰ and 28 ± 2 μg/g Li). Partition coefficients derived from IDDP-1 cuts successfully predict Li inventories in vitrophyres of rhyolites on the surface of the KVS. Lithium isotope compositions of the crystal-rich IDDP-1 cuts are significantly heavier (avg. δ⁷Li = +7.2 ± 0.2‰) than lavas and IDDP-1 obsidian cuts, casting doubt on the notion that the IDDP-1 rhyolitic magma could result from the melting of felsite lenses in the KVS. Lithium contents in groundmass glasses within IDDP-1 crystal-rich cuts show higher Li contents (avg. 55.1–60.7 μg/g), correlating with the higher crystal content and an increase in other incompatible elements (avg. 250 μg/g Rb) relative to obsidian cuttings (avg. 75 μg/g Rb).

If primitive arc magmas are primarily oxidized or if they acquire their oxidizing character during crustal evolution remains debated. Mineral-melt V partitioning is extremely redox sensitive and has the potential to resolve this debate, but critical low-temperature, hydrous experimental partitioning data are limited. We present new experimental V partitioning data for olivine-melt, spinel-melt, and clinopyroxene-melt based on the ≤ 1200 °C, hydrous basalt phase equilibrium experiments of Pichavant et al. (Geochim Cosmochim Acta 66:2193–2209, 2002) and Pichavant and Macdonald (Contrib Mineral Petrol 154:535–558, 2007). Combined with published experimental data, we use our olivine-melt V partitioning coefficients to show that—contrary to previous conclusions—hydrous melt composition, i.e. melt H₂O concentration in addition to the concentration of other melt components, affects V partitioning and thus calculated fO₂, systematically overestimating fO₂ for olivine-melt equilibrated at high melt H₂O concentration. False absolute fO₂ values and false oxidation or reduction trends may be inferred. Based on these findings, we have calibrated a set of new olivine-melt oxybarometers applicable to hydrous arc magmas partially crystallized at ~ 1025–1350 °C. In a case study on a high-Mg basalt from St. Vincent, Lesser Antilles, we show that olivine-melt V oxybarometry records oxidizing near-liquidus conditions (~ QFM + 1.8) and possibly limited oxidation (~ 0.3 log units) during intermittent crustal magma storage.

We have measured the dependence of the Fe–Mg interdiffusion coefficient, D_Fe-Mg, on the ferrosilite component in orthopyroxene, which so far has not been experimentally calibrated. Diffusion couples, consisting of approximately 1 µm thin-films were deposited by pulsed laser ablation on orthopyroxene crystals of En₉₁Fs₉ composition. Diffusion experiments were carried out at atmospheric pressure in vertical gas mixing furnaces (CO-CO₂) at temperatures between 950 and 1100 °C at constant fO₂ = 10^–7 Pa. Using a focused ion beam-scanning electron microscope (FIB-SEM), FIB-foils were cut from diffusion couples before and after annealing. Diffusion profiles were extracted by using combined backscattered electron (BSE) imaging and energy dispersive X-ray (EDXS) mapping on FIB-foils which allowed to resolve concentration gradients within 20 nm. The microstructure of the diffusion experiments was investigated using transmission electron microscopy (TEM). Using this method, we obtained the first experimentally determined data on the dependence of D_Fe-Mg on the ferrosilite content in orthopyroxene at different temperatures, appearing to increase with increasing a temperature. For the temperature range (950 – 1100 °C), log fO₂ = -7 Pa, along [001] in the composition Fs9, the log D_Fe-Mg yields the following Arrhenius equation:

The effect of X_Fe on D_Fe-Mg, given by (D(X_{Fe} ) = D(X_{Fe} = 0.09) cdot {10}^{m (X_{Fe} - 0.09)},) can be calculated by the following parameterization where m follows a linear regression of m on temperature:

By considering the D_Fe-Mg dependence on X_Fe, the timescales of natural processes obtained from modelling the compositional zoning of natural crystal may considerably differ from previously estimated timescales.

Rock textures observed via thin section are skewed from their true 3D nature. This is due to various cut effects—artifacts introduced due to the lower dimensional nature of the thin section relative to the rock. These cut effects can be corrected, and several methods have been developed to invert crystal shape and crystal size, but with each process performed separately and sequentially. With the ongoing adoption of electron backscatter diffraction (EBSD) by petrologists, an additional data stream has now become available: the 3D orientation of 2D grain sections. For EBSD analysis, no stereological corrections are typically applied for interpreting the data. This study tests whether this orientational information is skewed due to a fabric cut effect. We test this by numerically generating synthetic crystal datasets representative of several crystal shapes and population sizes. We find that EBSD orientational data has a fabric cut effect since crystals oriented with long axes perpendicular to the thin section are more likely to be sampled compared to those with long axes oriented parallel to it. This effect must be accounted for to interpret the true 3D fabric accurately. Towards this end, we develop two new tools for working with EBSD-derived fabric: (1) a simple first-order test for determining if a measured fabric exceeds that of the fabric cut effect, and (2) a method of inverting cut fabrics that provides robust error estimations. We demonstrate the applicability and accuracy of these methods using a range of synthetic examples and a natural sample. With these newly developed tools, there is clear potential for a new textural toolkit framework, to further our ability to correct for the various cut effects while also providing accurate uncertainty estimates.

We applied elastic thermobarometry on garnet-bearing migmatites along two transects through the crustal section at Sierra Valle Fértil-La Huerta, Argentina. We performed quartz-in-garnet barometry and zircon-in-garnet thermometry on metapelites from different paleo-depths across the crustal section. Our work recovers entrapment pressures ranging from 240 to 1330 MPa and entrapment temperatures between 691 and 1574 °C. The entrapment conditions are broadly consistent with anticipated pressures and temperatures along the crustal section derived previously using conventional, thermodynamic thermobarometers. The quartz-in-garnet barometer reproduces those conventionally established entrapment conditions when samples only experienced conditions within the alpha-quartz stability field. Raman-derived pressures for samples that experienced beta-quartz reference conditions are commonly much higher than those established by conventional barometry. Samples that preserve compressive (positive) residual pressures best reproduce reference entrapment pressures. Entrapment temperatures show high variability and overestimation of temperature conditions compared to conventional results. These results indicate elastic thermobarometry over- or under-estimates crystallization conditions in rocks crystallized at high temperatures, as is common in the Famatinian Arc deep-crust. We suggest that modeling quartz behavior across the alpha–beta transition may present challenges, as does shape maturation, viscous deformation, and radiation damage in zircon.

The origin of the wide range of sulfur isotope compositions (i.e., δ³⁴S) measured in arc rocks remains debated. While the observed δ³⁴S variability has been attributed to slab-related fluids that flux the sub-arc mantle, others have argued that it primarily reflects crustal-derived processes by some combination of magmatic differentiation, country rock assimilation, and/or degassing. Here, we present new whole rock sulfur isotopes for the Late Cretaceous Bear Valley Intrusive Suite (BVIS) that represents a continuous arc crustal section in the southern Sierra Nevada Batholith, exposing lower crustal mafic cumulates and cogenetic mid-upper crustal tonalites. Our data reveal a range of δ³⁴S-depleted values (–1.2 to − 5.1‰) for the BVIS with overlapping δ³⁴S between mafic cumulates and tonalites. Complementary δ³⁴S measurements of structurally concordant metasedimentary pendants indicate δ³⁴S-depleted values (–11.5 to − 5.2‰) for deep metasedimentary rocks compared to δ³⁴S-enriched values (+ 1.6 to + 6.4‰) for shallower ones. Quantitative mixing models suggest that assimilation of crustal-derived sulfur from metasedimentary rocks in the lower crust can account for the δ³⁴S-depleted values in the BVIS, whereas assimilation of shallower ones is unlikely. Sulfur degassing modelling indicates that the range of δ³⁴S-depleted values observed within mid-upper crustal tonalites can be reproduced by degassing  ~60–80% of the initial melt sulfur at fO₂ ≤ FMQ + 1 with initial H₂O content of 10–12 wt%. Finally, the identical ranges of δ³⁴S values within the tonalites and mafic cumulates argue for limited sulfur isotope fractionation related to magmatic sulfide immiscibility. Although assimilation, magma degassing and sulfide immiscibility are not mutually exclusive during crustal magmatic processes, field, thermal and geochemical evidence favor lower crustal-derived sulfur assimilation as the primary mechanism to explain the range of δ³⁴S- depleted values within the mafic cumulates, which are ultimately inherited by the derivative tonalitic melts. Overall, this study emphasizes that deep crustal magmatic processes can severely influence the early δ³⁴S evolution of arc magmas.

We performed ¹H and ²⁹Si NMR and infrared measurements, and first-principles calculations to clarify the nature of OH defects in MgSiO₃ orthoenstatite. An orthoenstatite sample synthesized at 7 GPa and 1200 °C from a composition of MgSiO₃ + 0.1 wt% H₂O yielded two ¹H MAS NMR peaks near 5.9 and 7.6 ppm that are correlated in 2D NMR spectra, and two infrared bands near 3361 and 3066 cm^− 1 that correspond to the previously reported A3 and A4 bands. The first-principles calculations confirmed that they are due to a pair of protons in a Mg (M2) vacancy. The previously reported A1 and A2 infrared bands near 3687 and 3592 cm^− 1 for orthoenstatite synthesized at low silica activities were confirmed to arise from four protons in a SiB vacancy. The latter is predicted to give two additional OH stretching bands associated with two strongly hydrogen-bonded O3b-H bonds with frequencies below the spectral range reported thus far. The previously reported infrared absorption coefficients were thus revised to account for the undetected bands. ¹H NMR may be used to quantitatively detect all four protons (expected at 1–12 ppm). Other mantle minerals should also be examined for potentially overlooked OH defects with strong hydrogen bonding.

This work investigates the relationships between partial melting, melt extraction, pluton growth and silicic volcanism in garnet-bearing felsic volcanic rocks that were extruded in the Iberian Pyrite Belt, at ca. 345 Ma. The garnets are of peritectic origin, displaying textural and chemical features of disequilibrium crystallization during partial melting reactions involving biotite at high temperatures (up to 870 °C) in the middle-lower crust. Major element composition suggests compositional equilibrium with the entrained and pinitized peritectic cordierite, but reveals some subsequent homogenization by diffusion. Trace element maps and spot analyses of garnet show, nonetheless, significant trace element variations, reflecting biotite and Y-REE-P-rich accessory phase breakdown during partial melting reactions. Peritectic garnet and cordierite growth resulted in the preservation of Th- and Y-rich prograde suprasolidus monazite, which constrains the timing of partial melting of the metapelitic protolith at ca. 356.8 ± 2.4 Ma. The zircon cargo further shows that a significant amount of zircon crystals from previously crystallized felsic melts were also remobilized and erupted. These were likely stored in an upper crustal pluton that grew episodically since ca. 390 Ma during voluminous melt generation periods within the middle to lower crust, which also resulted in voluminous volcanism. The geochemical trends of the felsic volcanic rocks reflect the entrainment of xenoliths of peritectic garnet, cordierite and feldspar, and as such, the garnet-bearing felsic volcanic rocks represent an erupted mixture of a lower-temperature (ca. 770 °C) silicic melt and autocrysts, and peritectic phases and zircon crystals from previously crystallized and stored felsic melts.

Orogenic peridotites in the crystalline basement of the northwestern Bohemian Massif contain multiphase solid inclusions (MSI), which are interpreted to be crystallisation products of trapped former carbonate–silicate melts metasomatizing their host rocks. We applied conventional thermobarometry and forward thermodynamic modelling to constrain the P–T evolution ranging from the peak metamorphic conditions of the investigated harzburgite and lherzolite, through entrapment of the melts in the outer parts of garnets, to the (re)-equilibration of the MSI assemblages. The peak conditions of c. 1100 °C/4.5–5.5 GPa are recorded by garnet cores and large pyroxene porphyroclasts. The melt entrapment, during which garnet outer parts grew, was associated with influx of the metasomatizing liquids and probably took place during the early stage of the exhumation. Thermodynamic model of amphibole-free MSI assemblage comprising kinoshitalite/Ba-rich phlogopite (approximated by phlogopite in the model), dolomite, magnesite, clinopyroxene, orthopyroxene, garnet and chromite provided robust estimate of P and T of its (re)-equilibration, c. 900–1000 °C, 1.8–2.2 GPa. Furthermore, the lack of olivine reflects co-existence of COH fluid with high X(CO₂) = CO₂/(CO₂ + H₂O) ≥ 0.6. Models employing identical P–T–X(CO₂) parameters successfully reproduced the other two amphibole-bearing assemblages observed. The modelled stability fields show perfect alignment with a characteristic isobaric segment of the solidus curve of carbonated peridotite. This co-incidence implies that the (re)-equilibration corresponds to the termination of the melt crystallisation once the near-isothermal exhumation path intersected the solidus. Decreased solubility of silicates at the carbonated peridotite “solidus ledge”, inferred from the published experimental data, as well as concentric textures of some MSI indicates sequential crystallisation from the early silicates to late dolomite. The carbonated “solidus ledge” is a relatively narrow boundary in the lithospheric mantle capable of an abrupt immobilisation of fluxing or transported carbonated melts. The investigated rocks are estimated to store approximately 0.02 kg C/m³ (or 6 ppm C) occurring as carbonates in the MSI.

Understanding the compositional heterogeneity of the deep mantle requires identifying the nature of recycled crustal materials in the sources of mantle-plume-related magmas. However, it is still unclear whether or not the deep mantle contains recycled carbonates from the Earth’s surface. In this study, we present comprehensive data on whole-rock high-precision zinc isotopes, as well as major- and trace-element geochemistry, and Sr–Nd–Pb isotopes of basalts on Hainan Island to examine the influence of recycled materials (particularly carbonates) on the mantle source heterogeneity of the Hainan mantle plume. The basalts have highly variable δ⁶⁶Zn values ranging from 0.21‰ to 0.42‰. These variable Zn isotopic compositions cannot be accounted for by processes such as post-magmatic alteration and crustal contamination, or by fractional crystallization and partial melting; instead, they reflect mantle heterogeneity. Comparisons of the major- and trace-element compositions (e.g., CaO and TiO₂ contents and Zn/Fe and Fe/Mn ratios), FC3MS and FCKANTMS peridotite and pyroxenite melting parameters, as well as pseudo-ternary projections of the primary Hainan basaltic magmas with experimental data suggest that the primary magmas were partial melts of silica-deficient pyroxenitic lithologies with peridotite residue. The heterogeneous geochemical and lithological compositions of the Hainan basalts indicate that recycled sedimentary carbonates and siliceous rocks were important constituents in their mantle source. Quantitative modeling reveals that the addition of 1–10% subducted sediments into the source of the Hainan basalts closely reproduces their Zn–Sr–Nd–Pb isotopic values. The source component with the heaviest Zn isotopic composition measured for the Hainan basalt samples could have contained more than 9% recycled carbonate. Our findings provide insights into the role of subducted materials to mantle heterogeneity and highlight the contribution of subducted sedimentary carbonates in the deep recycling of oceanic slabs, including—in the case of the Hainan mantle plume—recycled deep mantle (i.e., the mantle transition zone and lower mantle).