Proportions, Timing, and Re-equilibration Progress During the 1959 Summit Eruption of Kīlauea: an Example of Magma Mixing Processes Operating During OIB Petrogenesis

Petrographic and chemical analysis of scoria samples collected during the 1959 Kīlauea summit eruption illustrate the progress of thermal and chemical homogenization of the melts, and the gradual growth and/or re-equilibration of olivine phenocrysts, over the course of the eruption. Glass compositions show that thermal equilibration was largely complete within the span of the eruption, while chemical homogenization was a work in progress. The olivine phenocryst population, known to contain conspicuous antecrystic components, is also hybrid within the euhedral population. The bulk of the olivine reached the level of the erupting magma on November 18-19, 1959. Zoning patterns in olivine phenocrysts show that initially unzoned grains developed normal zoning by the end of the eruption. Reverse zoning in relatively Fe-rich olivine phenocrysts (interpreted as cognate to the stored magma) was progressively eliminated from November 21 to December 19, 1959, by diffusive re-equilibration between crystals and melt. Toward the end of the eruption, the only olivine composition in direct contact with the melt was Fo84-86, with the original rim compositional heterogeneity gone in 4-5 weeks’ time. Activity in December 1959 differed from that in November, as high fountaining events were more closely spaced and almost all samples were picritic, with bulk MgO ≥16.5 wt %. Three different levels were in play during the 1959 eruption: a deep source for high-MgO melts and forsteritic (Fo87-89) olivines, an intermediate source for the bulk of the stored magma, and a shallower source for the most differentiated magma. This model is consistent with geophysical, petrologic and chemical observations. Comparison of the 1959 eruption with results from older explosive deposits suggest that stored and recharge melts and olivine from the deeper parts of Kīlauea’s plumbing are similar in composition to those observed or inferred in the 1959 eruption, so they behave similarly during extrusive and explosive periods alike.

Establishing genetic relationships between the Takidani pluton and two large silicic eruptions in the Northern Japan Alps

The Takidani pluton (1.1-1.6 Ma) represents a shallow magmatic reservoir at the base of an exhumed caldera floor. The deposits of two large caldera-forming eruptions including the Nyukawa Pyroclastic Flow Deposit (1.76 Ma; crystal-rich dacite) and the Chayano Tuff and Ebisutoge Pyroclastic Deposits (1.75 Ma; a sequence of crystal-poor rhyolite) are distributed concentrically around the pluton. We use major and trace element chemistry of whole-rock, glass and minerals to show (1) that the crystal-rich dacite (>400 km3 DRE; dense rock equivalent) is the erupted portion of a shallow mush zone constituting the Takidani pluton and (2) that the crystal-poor rhyolite (>100 km3 DRE) was extracted from a deeper part of this vertically extended magmatic plumbing system. Whole-rock geochemistry indicates that the Nyukawa and Takidani compositions were produced dominantly through crystal fractionation of amphibole, pyroxene and plagioclase in the mid-to-lower crust and subsequently emplaced in the upper crust prior to eruption and solidification, respectively. The crystal-poor Chayano-Ebisutoge rhyolite (>100 km3 DRE) is compositionally distinct from the Nyukawa and Takidani magmas and its generation is associated with a substantial contribution of crustal melts. Yet, plagioclase and orthopyroxene textures and chemistry provide strong evidence that the ascending rhyolite percolated through the upper Takidani-Nyukawa mush zone prior to eruption. Overgrowth of “rhyolitic plagioclase” on “xenocrystic dacitic plagioclase” typical of the Takidani-Nyukawa magmas indicates that the extraction and accumulation of the rhyolitic melts could have occurred in less than 10 kyr (i.e. time between eruptions) prior to eruption providing maximum timescales for pre-eruption storage. Overall, our findings show a progressive growth and thermal maturation of a vertically extended magmatic plumbing system over hundreds of thousands of years and imply that large volcanic eruptions can occur in relatively short succession without dramatic changes in the plumbing system, thus, complicating the identification of signs of an impending large eruption.

Multi-stage metamorphism of the South Altyn ultrahigh-pressure metamorphic belt, West China: insights into tectonic evolution from continental subduction to arc–backarc extension

The South Altyn ultrahigh-pressure (UHP) metamorphic belt is claimed to host the deepest subducted continental crust based on the discovery of former stishovite, and thus can provide unique insights into the tectonic evolution from deep continental subduction and exhumation to arc–backarc extension. In this paper, we present detailed studies of petrography, mineral chemistry, phase equilibria modelling and zircon U-Pb dating for three representative samples involving garnet amphibolite (A1531 & A1533) and associated garnet-biotite gneiss (A1534) from the UHP belt. Three phases of metamorphism are inferred for the rocks. The first phase high pressure (HP)–UHP-type eclogite facies is represented by the mineral assemblages of garnet and phengite inclusions in zircon and garnet cores with the high grossular (XGrs = 0.33–0.34). The Si contents of 3.40–3.53 and 3.24–3.25 p.f.u. in phengite inclusions yield pressure conditions of >1.7–2.3 GPa for A1533 and 2.5–2.55 GPa for A1534 at a fixed temperature of 770 °C. The second phase medium-pressure (MP)-type overprinting of garnet amphibolite facies shows P–T conditions of 0.8–1.2 GPa/750–785 °C based on the stability fields of corresponding mineral assemblages, the measured isopleths of Ti contents in biotite and amphibole cores, and XGrs in garnet. The third phase low-pressure (LP) type overprinting includes early-stage heating to peak granulite facies followed by cooling towards a late-stage amphibolite facies. The peak granulite facies is represented by the high Ti amphibole mantle, high Zr titanite and the intergrowths of clinopyroxene + ilmenite in A1533 & A1531, with P–T conditions of 800–875 °C/0.80–0.95 GPa. The late-stage is defined by the solidus assemblages, giving P–T conditions of 0.5–0.7 GPa/720–805 °C. U-Pb geochronology on metamorphic zircons from A1533 and A1534 gives three ages of c. 500 Ma, c. 482 Ma and c. 460 Ma. They are interpreted to represent the HP–UHP, MP and LP types of metamorphism respectively, based on cathodoluminescence images, mineral inclusions and trace element patterns. Combining the regional geology and metamorphic evolution from the Altyn Orogen, a tectonic model is inferred, including the following tectonic scenarios. The small Altyn Microcontinent was subducted to great mantle depths with dragging of the surrounding vast oceanic lithosphere to undergo the HP–UHP eclogite facies metamorphism during the early subduction stage (c. 500 Ma) of the Proto-Tethys Ocean. Then, the subducted slabs were exhumed to a thickened crust region to be overprinted by the MP-type assemblages at c. 482 Ma. Finally, an arc–backarc extension was operated within the thickened crust region due to the retreat of subduction zones. It caused evident heating and the LP-type metamorphic overprinting at c. 460 Ma, with a fairly long interval of 30–40 Myr after the HP–UHP metamorphism, distinct from the short interval of <5–10 Myr in the Bohemian Massif.

Variable modes of formation for tonalite-trondhjemite-granodiorite-diorite (TTG)-related porphyry-type Cu ± Au deposits in the Neoarchean southern Abitibi subprovince (Canada): Evidence from petrochronology and oxybarometry

Most known porphyry Cu ± Au deposits are associated with moderately oxidized and sulfur-rich, calc-alkaline to mildly alkalic arc-related magmas in the Phanerozoic. In contrast, sodium-enriched tonalite-trondhjemite-granodiorite-diorite (TTG) magmas predominant in the Archean are hypothesized to be unoxidized and sulfur-poor, which together preclude porphyry Cu deposit formation. Here, we test this hypothesis by interrogating the causative magmas for the ~2.7 Ga TTG-related Côté Gold, St-Jude, and Clifford porphyry-type Cu ± Au deposit settings in the Neoarchean southern Abitibi subprovince. New and previously published geochronological results constrain the age of emplacement of the causative magmas at ~2.74 Ga, ~2.70 Ga, and ~2.69 Ga, respectively. The dioritic and trondhjemitic magmas associated with Côté Gold and St-Jude evolved along a plagioclase-dominated fractionation trend, in contrast to amphibole-dominated fractionation for tonalitic magma at Clifford. Analyses of zircon grains from the Côté Gold, St-Jude, and Clifford igneous rocks yielded εHf(t) ± SD values of 4.5 ± 0.3, 4.2 ± 0.6, and 4.3 ± 0.4, and δ18O ± SD values of 5.40 ± 0.11 ‰, 3.91 ± 0.13 ‰, and 4.83 ± 0.12 ‰, respectively. These isotopic signatures indicate that although these magmas are mantle-sourced with minimal crustal contamination, for the St-Jude and Clifford settings the magmas or their sources may have undergone variable alteration by heated seawater or meteoric fluids. Primary barometric minerals (i.e., zircon, amphibole, apatite, and magnetite-ilmenite) that survived variable alteration and metamorphism (up to greenschist facies) were used for estimating fO2 of the causative magmas. Estimation of magmatic fO2 values, reported relative to the fayalite-magnetite-quartz buffer as ΔFMQ, using zircon geochemistry indicate that the fO2 values of the St-Jude, Côté Gold, and Clifford magmas increase from ΔFMQ -0.3 ± 0.6, ΔFMQ +0.8 ± 0.4, to ΔFMQ +1.2 ± 0.4, respectively. In contrast, amphibole chemistry yielded systematically higher fO2 values of ΔFMQ +1.6 ± 0.3 and ΔFMQ +2.6 ± 0.1 for Côté Gold and Clifford, respectively, which are consistent with previous studies that indicate amphibole may overestimate the fO2 of intrusive rocks by up to one log unit. Micro X-ray absorption near edge structure (μ-XANES) spectrometric determination of sulfur (i.e., S6+/ΣS) in primary apatite yielded ≥ΔFMQ -0.3 and ΔFMQ +1.4–1.8 for the St-Jude and Clifford, respectively. The magnetite-ilmenite mineral pairs from the Clifford tonalite yielded ΔFMQ +3.3 ± 1.3 at equilibrium temperatures of 634 ± 21 °C, recording the redox state of the late stage of magma crystallization. Electron probe microanalyses revealed that apatite grains from Clifford are enriched in S (up to 0.1 wt. %) relative to those of Côté Gold and St-Jude (below the detection limit), which is attributed to either relatively oxidized or sulfur-rich features of the Clifford tonalite. We interpret these results to indicate the deposits at Côté Gold and Clifford formed from mildly (~ΔFMQ +0.8 ± 0.4) to moderately (~ΔFMQ +1.5) oxidized magmas where voluminous early sulfide saturation was probably limited, whereas the St-Jude deposit represents a rare case whereby the ingress of externally derived hydrothermal fluids facilitated metal fertility in a relatively reduced magma chamber (~ΔFMQ +0). Furthermore, we conclude that variable modes of formation for these deposits and, in addition, the apparent rarity of porphyry-type Cu-Au deposits in the Archean may be attributed to either local restriction of favorable metallogenic conditions, and/or preservation, or an exploration bias.

Late Eocene two-pyroxene trachydacites from the southern Qiangtang Terrane, central Tibetan Plateau: High-temperature melting of overthickened and dehydrated lower crust

Orthopyroxene-bearing granitic rock (e.g., charnockite) is relatively rare but provides an excellent opportunity to probe the thermal and tectonic evolution of deep orogenic crust because of its distinct mineral assemblage. Here we present petrological, mineralogical, elemental, and Sr–Nd–Hf–O isotopic data for late Eocene (ca. 36 Ma; zircon U–Pb ages) volcanic rocks exposed in the Ejiu region in the southern Qiangtang Terrane to investigate how the central Tibetan crust evolved to its modern thickness and thermal state. The Ejiu volcanic rocks (EVRs) are trachydacites with anhydrous mineral assemblages (i.e., two pyroxenes, sanidine, plagioclase, and ilmenite, without amphibole and biotite) and geochemical characteristics (e.g., high P2O5 and TiO2) that resemble those of charnockite-type magmatic rocks. Mineral and whole-rock thermometry and hygrometry suggests that the parent magma crystallized under hot (~1000 °C) and dry (H2O < 2 wt.%) condition. Besides, the EVRs display adakitic affinities according to their high SiO2 and Al2O3 contents, high Sr/Y, La/Yb, and Gd/Yb ratios, and low Y and Yb contents, without marked negative Eu anomalies. The calculated melts in equilibrium with pyroxenes also display adakitic compositions (e.g., high Sr/Y and La/Yb ratios), indicating that the adakitic compositions of the EVRs did not result from late-stage magmatic evolution. In addition, the melts of the EVRs were saturated in TiO2, as inferred from the high TiO2 contents of these rocks and the presence of ilmenite. An integrated analysis of the geochemical, petrological, and mineralogical data suggests that the EVRs were neither evolutional products nor partial melts of hydrous mafic materials at normal crustal pressures, but were formed by fusion of an eclogitized mafic protolith with residue containing garnet and rutile but lacking amphibole and plagioclase. The whole-rock Sr–Nd and zircon Hf isotope compositions of the EVRs [(87Sr/86Sr)i = 0.7053 to 0.7066; εNd(t) = −1.40 to −0.99; zircon εHf(t) = +1.08 to +5.31] indicate that the parental protolith was relatively juvenile in nature, but also contained some supracrustal materials given the high zircon δ18O values [zircon δ18O = +8.21‰ to +11.00‰]. The above arguments lead us to propose that of partial melting of a previously dehydrated—but chemically undepleted—mafic lower continental crust at high pressure (>1.5 GPa) and high temperature (>1000 °C) generated the EVRs. Based on a synthesis of independent geological and geophysical data, we further suggest that the southern Qiangtang Terrane crust of the central Tibetan Plateau was thick, dry, and elevated during the Late Cretaceous to early Eocene time, and that it became abnormally hot owing to the ascending asthenosphere after lithospheric foundering during the middle Eocene.

Melt Percolation, Melt-Rock Reaction and Oxygen Fugacity in Supra-Subduction Zone Mantle and Lower Crust from the Leka Ophiolite Complex, Norway

Samples of peridotites and pyroxenites from the mantle and lower crustal sections of the Leka Ophiolite Complex (LOC; Norway) are examined to investigate the effects of melt-rock reaction and oxygen fugacity variations in the sub-arc oceanic lithosphere. The LOC is considered to represent supra-subduction zone (SSZ) oceanic lithosphere, but also preserves evidence of pre-SSZ magmatic processes. Here we combine field and microstructural observations with mineral chemical and structural analyses of different minerals from the major lithologies of the LOC. Wehrlite and websterite bodies in both the mantle and lower crust contain clinopyroxene likely formed at a pre-SSZ stage, characterised by high Al, high Cr, low Mg crystal cores. These clinopyroxenes also exhibit low Al, low Cr, high Mg outer rims and intracrystalline dissolution surfaces, indicative of reactive melt percolation during intrusion and disruption of these lithologies by later, SSZ-related, dunite-forming magmas. Chromian-spinel compositional variations correlate with lithology; dunite-chromitite Cr-spinels are characterised by relatively uniform and high TiO2 and Al2O3, indicating formation by melt-rock reaction associated with SSZ processes. Harzburgite Cr-spinel compositions are more variable but preserve a relatively high Al2O3, low TiO2 endmember that may reflect crystallisation in a pre-SSZ oceanic spreading centre setting. An important finding of this study is that the LOC potentially preserves the petrological signature of a transition between oceanic spreading centre processes and subsequent supra-subduction zone magmatism. Single crystal Cr-spinel Fe3+/ΣFe ratios calculated on the basis of stoichiometry (from electron microprobe [EPMA] and crystal structural [X-ray diffraction; XRD] measurements) correlate variably with those calculated by point-source (single crystal) Mössbauer spectroscopy. Average sample EPMA Fe3+/ΣFe ratios overestimate or underestimate the Mössbauer-derived values for harzburgites, and always overestimate the Mössbauer Fe3+/ΣFe ratios for dunites and chromitites. The highest Fe3+/ΣFe ratios, irrespective of method of measurement, are therefore generally associated with dunites and chromitites, and yield calculated log(fO2)FMQ values of up to ~+1.8. While this lends support to the formation of the dunites and chromitites during SSZ-related melt percolation in the lower part of the LOC, it also suggests that these melts were not highly oxidised, compared to typical arc basalts (fO2FMQ of >+2). This may in turn reflect the early (forearc) stage of subduction zone activity preserved by the LOC and implies that some of the arc tholeiitic and boninitic lava compositions preserved in the upper portion of the ophiolite are not genetically related to the mantle and lower crustal rocks, against which they exhibit tectonic contacts. Our new data also have implications for the use of ophiolite chromitites as recorders of mantle oxidation state through time; a global comparison suggests that the Fe3+/ΣFe signatures of ophiolite chromitites are likely to have more to do with local environmental petrogenetic conditions in sub-arc systems than large length-scale mantle chemical evolution.

Phase equilibrium and trace-element modeling of the partial melting of basaltic rocks under low pressure: Implications for plagiogranite generation

Phase equilibria and trace-element modeling using two previously reported basaltic bulk-rock compositions (samples D11 and 104-16), were carried out in this study, in order to better understand mechanism of low-pressure (LP) partial melting of mafic rocks and associated melt compositions. The T–MH2O pseudosections for both samples at three pressures (i.e. 0.5, 1.0 and 2.0 kbar) display that the H2O-stability field gradually increased with decreasing pressure within the T–MH2O range of 600–1100 °C and 0–12 mol.%. The H2O contents of 10, 5.0, and 0.5 mol.% were selected on the basis of the T–MH2O pseudosections to calculate P–T pseudosections over a P–T window of 0.1–3 kbar and 600–1100 °C, so that the reactions of both the H2O-fluxed and -absent meltings at LP conditions can be investigated. The solidus displays a negative or near-vertical P–T slope, and occurs between 710 and 900 °C at pressure between 0.1 and 3.0 kbar. LP melting of metabasites is attributed to the reactions of the hydrous mineral (hornblende and/or biotite) melting and anhydrous mineral (plagioclase, orthopyroxene, and augite) melting. The hydrous mineral melting is gradually replaced by anhydrous mineral melting as pressure decreasing, as the stability of hornblende decreases with falling pressure. With increasing temperature at a given pressure, the modeled melt compositions are expressed as progressions of the granite-granodiorite-gabbroic diorite fields for sample D11and granite-quartz monzonite-monzonite-gabbroic diorite fields for sample 104-16 on the total alkali–silica diagram. The modeled melts produced through the H2O-fluxed melting display higher Al2O3, CaO, MgO, and lower SiO2 and K2O than those formed by H2O-absent melting at the same P–T conditions. Furthermore, the modeled melts formed by H2O-absent melting, become richer in Al2O3, CaO, MgO, FeO, Na2O, but poorer in SiO2 and K2O as increasing water content. The results of trace-element modeling suggests that the nearly flat REE patterns of modeled bulk-rock composition are inherited by all the modeled melts, and the negative Eu anomalies and Sr depletion of the modeled melts gradually decrease as melting degree increases. Combined with the geochemical characteristics of natural oceanic plagiogranites, which have low K2O contents and flat or slightly LREE-depleted REE patterns, our results imply that a bulk-rock composition with low K2O (<0.17 wt.%) and slightly LREEs depletion is the most likely protolith composition (e.g. basalt D11) for plagiogranites, and the compositions of modeled melts formed by LP H2O-absent partial melting of the basalt D11 at relatively high temperatures (1000–1025 °C) are coincident with those of 1256D tonalites.

The Tethyan Himalaya igneous province: Early melting products of the Kerguelen mantle plume

The Kerguelen large igneous province (LIP) has been related to mantle plume activity since at least 120 Ma. There are some older (147–130 Ma) magmatic provinces on circum-eastern Gondwana, but the relationship between these provinces and the Kerguelen mantle plume remains controversial. Here we present petrological, geochronological, geochemical, and Sr–Nd–Hf–Pb–Os isotopic data for high-Ti mafic rocks from two localities (Cuona and Jiangzi) in the eastern Tethyan Himalaya igneous province (147–130 Ma). Zircon grains from these two localities yielded concordant weighted mean 206Pb/238U ages of 137.25 ± 0.98 and 131.28 ± 0.78 Ma (2σ), respectively. The analyzed mafic rocks are enriched in high field strength elements and have positive Nb–Ta anomalies relative to Th and La, which have ocean island basalt-like characteristics. The Cuona basalts were generated by low degrees of melting (3–5%) of garnet lherzolites (3–5 vol.% garnet), and elsewhere the Jiangzi diabases were formed by relatively lower degrees of melting (1–3%) of garnet lherzolite (1–5 vol.% garnet). The highly radiogenic Os and Pb isotopic compositions of the Jiangzi diabases were produced by crustal contamination, but the Cuona basalts experienced the least crustal contamination given their relatively low γOs(t), 206Pb/204Pbi, 207Pb/204Pbi, and 208Pb/204Pbi values. Major and trace element geochemical and Sr–Nd–Hf–Pb–Os isotope data for the Cuona basalts are similar to products of the Kerguelen mantle plume head. Together with high mantle potential temperatures (>1500°C), this suggests that the eastern Tethyan Himalaya igneous province (147–130 Ma) was an early magmatic product of the Kerguelen plume. A mantle plume initiation model can explain the temporal and spatial evolution of the Kerguelen LIP, and pre-continental break-up played a role in the breakup of eastern Gondwana, given the >10 Myr between initial mantle plume activity (147–130 Ma) and continental break-up (132–130 Ma). Like studies of Re-Os isotopes in other LIPs, the increasing amount of crustal assimilation with distance from the plume stem can explain the variations in radiogenic Os.