Magma Composition Research Articles

Geochemical Distribution of Platinum Metals, Gold and Silver in Intrusive Rocks of the Norilsk Region

The Norilsk ore district is one of the world leaders in the production of platinum metals. Long-term research focused on the detection of sulfide platinum-copper-nickel ores contributed to the accumulation of a large volume of scientific material on the geology and mineralization of the Norilsk area. Despite this, the issue of the composition of the initial melt for ore-bearing intrusive complexes and its degree of enrichment with noble metals remains open. Intrusive rocks of the Norilsk region are rarely analyzed for their ratio of noble metals. However, the analysis and comparison of geochemical parameters of different types of intrusions allows us to draw important conclusions not only about the composition of the initial magmas of ore-bearing complexes, but also about the formation conditions of the intrusions. This study demonstrates the distribution of platinum metals, gold and silver in the main petrographic differentiates of the Kharaelakh, Talnakh, Vologochan intrusions and Kruglogorsk-type intrusion. The regularities and variations of the distribution of metals depend on the host rocks. There are two series of rocks in the inner structure of the ore-bearing intrusions: 1. Picritic and taxitic gabbro-dolerites enriched in PGE-Au-Ag mineralization which forms disseminated ores at intrusion bottoms (ore-bearing rocks). 2. Olivine-, olivine-bearing, olivine-free gabbro-dolerites and leucogabbro with poor sulfide mineralization at the upper part of the intrusions (ore-free rocks). There is a distinct correlation between PGE, Cu, S and to a lesser extent correlation with Ni in the first rock group, which is a characteristic of sulfide PGE-Cu-Ni deposits. In the second group, correlations are also revealed, but the correlation coefficients are lower. The main element controlling the distribution of platinum metals is copper. The taxitic gabbro-dolerites of the Talnakh intrusion are the most enriched by noble metals. According to noble metal patterns the rocks of the Kharaelakh intrusion show the highest degree of melting of the initial mantle material during the formation of parental magmas chambers. Despite some differences, the geochemical features of the studied rocks indicate the similar characteristics of the accumulation of gold, silver and platinum metals in the intrusions of the Talnakh, Kruglogorsk and Zubovsk types, which allow suggesting the close conditions for the formation of ore mineralization of these intrusions.

Magma Composition and Formation Conditions of Alkaline-Salic Rocks of the Early Mesozoic Adaatsag Bimodal Association of the Kharkhorin Rift Zone, Central Asia: Data from the Study of Melt Inclusions in Minerals

The composition and formation conditions of magmas of comendites of the Early Mesozoic Adaatsag volcanic association are determined using electron and ion microprobe analysis of melt inclusions, and mechanisms favorable for the accumulation of trace and rare earth elements are identified. The melt and fluid inclusions are found in quartz of comendites sampled in various parts of the section of a volcanic sequence. The melt inclusions consist of glass, a gas bubble, and daughter minerals (fluorite, polylithionite, and K-feldspar). The composition of a gas phase in melt inclusions was studied using Raman spectroscopy. The dominant gas phase components include water and hydrogen. The fluid inclusions are aqueous fluids probably corresponding to KF. The thermometric studies of melt inclusions and the analysis of glasses of homogenized melt inclusions in quartz phenocrysts of comendites showed that magmas crystallized from water-saturated rare metal melts with high Li, Zr, F, Rb, Nb, Y, and Th contents at a temperature of 880–930°C and a pressure of 1000 bar at a depth of ~3.5 km accompanied by degassing. The comparison of our data of the study of melt inclusions in phenocrysts of alkaline salic rocks of the coeval Adaatsag, Dzarta-Khuduk, and Sant volcanic associations, which occur within the Kharkhorin Rift Zone, revealed common principles of magma evolution. This allowed us to propose a similar mechanism of their formation, involving the accumulation of many rare and rare-earth elements, as well as volatile components (F, H2O) in the process of crystallization differentiation. Subsequently, a salt melt rich in Li, F and water could be separated from such comendite magmas. The detection of fluoride aqueous inclusions in quartz allows us to suppose the further evolution of the salt melt leading to the appearance of a concentrated aqueous fluid and the possible participation of the latter in metasomatic processes.

A machine learning method for distinguishing detrital zircon provenance

Zircon geochemistry provides a sensitive monitor of its parental magma composition. However, due to the complexity of the uptake of trace elements during zircon growth, identifying source magmas remains challenging, particularly for detrital grains whose petrological context is lost. We use a machine learning-based approach to explore the classifiers for zircon provenance, based on 3794 published, high-quality zircon trace element analyses compiled from I-, S-, and A-type granites. Three supervised machine learning algorithms, namely, Support Vector Machine (SVM), Random Forest (RF), and Multilayer Perceptron (MLP) were used and trained with 11 features, including 7 trace elements (Ce, Eu, Ho, Nb, Ta, Th, and U) and 4 derived trace element ratios (Th/U, U/Yb, Ce/Ce*, and Eu/Eu*). Our results show that all three trained machine learning methods perform very well with accuracy varying from 0.86 to 0.89, and that input–output relationships captured by different ML methods are nearly consistent and can be explained by the known petrological processes. The application of our trained machine learning classifiers to detrital zircon studies will enhance the interpretability of zircon assemblages of different origins. It also helps develop interpretations, approaches, and tools that will benefit, for example, the study of continental crust evolution and mineral exploration.

Explosive volcanism of Piton des Neiges (Reunion Island) and excess age dispersion in sanidine: Insights into magma chamber processes in a hotspot setting

Explosive eruptions of trachytic to trachyandesitic composition produced thick and widespread pyroclastic deposits during the recent volcanic activity of Piton des Neiges (PdN), the largest volcano of La Réunion Island (France). However, without precise geochronology for these eruptions, the future behavior of the dormant volcano is difficult to assess. New high-resolution single crystal 40Ar/39Ar dates on alkali feldspar, together with petrochemical data on juvenile clasts, indicate that instead of multiple and small eruptions, PdN produced a single large-volume event (“Dalle Soudée” eruption) at 198.8 ± 2.5 kyrs in which its entire crustal plumbing system was reactivated by mafic injections in the silicic reservoir. After 60 kyrs of quiescence, PdN activity resumed with several small silicic eruptions. From a dome-forming trachytic event (Belouve eruption) represented by some of the most recent block-and-ash flow deposits near the present-day summit we report in-situ 238U/230Th disequilibrium dates on pristine zircons. Sanidine and zircon crystalized at ca. 875 °C in the trachytic magma. For this event, the corresponding 238U/230Th isochron age of the zircons overlaps with the youngest 40Ar/39Ar dates on single sanidine crystals and yields an age of 42.7 ± 3.8 kyr. The variability in trace elements of zircon argues for a prolonged crystallization during this phase prior to eruption. Similarly, the 40Ar/39Ar sanidine dates show a significant spread in both eruptions, covering ca. 80 kyrs, uncommon for a method with a predicted low closure temperature. Such a pattern of age distribution argues for the presence of inherited argon in sanidines and for long residence times in the magma reservoir. The recent eruptive history of Piton des Neiges, involving evolved magma compositions, long-living silicic reservoirs and millennia-long quiescence in volcanic activity, suggests that future potentially hazardous eruptions from this volcano cannot be excluded.

The nature of the Karoo mantle source region from the perspective of olivine in the Luenha picrites: an oxygen isotopic and chemical study

The Karoo large igneous province has been divided into rift zone and basin-related groups, with picrites from the Luenha river, Mozambique, representing an end-member of the latter. New O isotope, major and trace element data for olivine have been combined with MELTS crystallisation modelling to deconvolve compositional diversity associated with magma differentiation from source-derived heterogeneity. Three olivines populations have been discerned as follows: (1) the ‘main trend’, which records crystallisation from a variety of magma compositions; (2) the ‘low Cu trend’, which is inferred to represent xenocrysts or antecrysts; and (3) the ‘high CaO’ olivines, which record polybaric crystallisation of a primitive, little fractionated magma. The trace element variability in olivine phenocrysts relates partially to sampling of different parts of the same overall magma transport and storage systems, and partly to heterogeneity of parental magmas and their mantle sources. When the measured δ18Oolivine values have been converted into δ18Omelt values, the mean δ18Omelt values for the ‘main trend’ and ‘low Cu’ groups are indistinguishable from each other (5.7 ± 0.1‰, 2σ); however, the mean δ18Omelt value of 6.1 ± 0.1‰ for the ‘high CaO’ group is distinctly enriched. These data record source heterogeneity and suggest contributions from two mantle sources, one with elevated δ18O, and another with more ‘typical’ mantle δ18O. Combining these data with previously reported trace element and Nd and Sr isotope data support derivation from a mantle source similar to non-chondritic bulk silicate earth, but with minor contributions (1–2% for the enriched magmas) from a recycled sedimentary component. This points to the importance of a primitive mantle source for the basin-related successions in Karoo.

Layered Intrusions of Paleoproterozoic Age in the Kola and Karelian Regions

Large-scale layered intrusions of a peridotite–pyroxenite–gabbronorite complex, to which Cr, Ni, Cu, and PGE deposits and ore occurrences are confined, were emplaced into the Baltic paleocontinent 2.50–2.45 Ga. Layered intrusions in the Monchegorsk Ore District, including the Monchepluton and Imandra–Umbarechka Complex, as well as the gabbro-anorthosite complex of the Main Ridge, were analyzed earlier geochemically and isotopically. In the present paper, the authors analyze layered intrusions in the Kola region (Mount Generalskaya) and Karelia (Kivakka, Kovdozero, and the Burakovsky Pluton). The primary composition of mantle magmas for the layered intrusions is assumed to be identical to that of the komatiitic basalts making up the volcanogenic units of the Vetreny Belt and the Imandra–Varzuga zone. A general model for the formation of layered intrusions includes superplume uplift in the early Paleoproterozoic, the generation of mantle magmas and their injection into the lower portion of the earth crust, the formation of deep-seated and intermediate magma chambers, and the intense contamination of the granulite–metamorphic complex followed by the generation of magma chambers provoked by single or multiple injections.

The Fundamental Role of H2O in the Generation of Coeval Sodic and Potassic Granitoids at Continental Arcs: An Example from the Yangtze Block, South China

AbstractThe bulk rock composition of granitoids reflects the composition of their source and the conditions of partial melting, which are functions of the geodynamic setting in which they formed. Granitoids in active continental margins (continental arcs) are dominated by calc-alkaline rocks with subordinate alkaline compositions, although how these different magma compositions formed is not well understood. Neoproterozoic magmatic rocks are widely distributed along the western margin of the Yangtze Block in South China to form the >1000-km long Panxi continental arc system, which is dominated by granitoids with minor mafic–ultramafic and intermediate plutons. The granitoids are subdivided into sodic and potassic variants that occur as belts along the western and eastern sides of the continental arc, respectively. Sodic granitoids from the western part consist of tonalite, granodiorite, and monzogranite with crystallisation ages ranging from 870 Ma to 740 Ma. They have low K2O/Na2O ratios (0.1–1.0) and high Na2O contents (3.5–6.7 wt%), high but variable SiO2 (61–75 wt%) concentrations, and negative to positive whole-rock εNd(t) values (−1.7 to +2.9). Zircon grains from the sodic granitoids have εHf(t) values ranging from +0.3 to +9.6 and δ18O from 3.90‰ to 7.71‰. The potassic granitoids from the eastern side consist of monzogranite and syenogranite with crystallisation ages from 820 Ma to 790 Ma. They have high K2O/Na2O ratios (0.6–2.2), K2O (2.6–5.9 wt%) and SiO2 contents (69–78 wt%), but whole-rock εNd(t) (−0.9 to +2.9) and zircon εHf(t) (+1.8 to +12.9), and δ18O values (2.98‰ to 6.41‰) similar to those of the sodic granitoids. The isotopic compositions of both the sodic and potassic granitoids are similar to those of spatially- and temporally-related mantle-derived (mafic to ultramafic) rocks, and are considered to have been derived from juvenile mafic continental crust. Phase equilibrium modelling shows that the H2O content of the granitoid source rocks played a key role in their petrogenesis, both in lowering solidus temperatures and in controlling the compositions of the derived partial melts. Our results indicate that calc-alkaline sodic granitoids can be formed by water-fluxed melting of juvenile mafic crust at 750–900°C and 9–12 kbar in which the required H2O was derived from the dewatering of underplating mafic arc magmas. By contrast, the potassic granitoids were generated by fluid-absent (H2O-undersaturated) partial melting of a similar juvenile mafic source at 725–900°C and 6–9 kbar. We conclude that the sodic granitoids were derived from partial melting of the newly-formed mafic lower crust in the continental arc, whereas the potassic granitoids were likely generated in the back-arc setting induced by upwelling of asthenospheric mantle.

Petrogenesis of Early Cretaceous volcanic rocks in Santos Basin, offshore SE Brazil: evidence for the Tristan da Cunha plume

ABSTRACT The Santos Basin comprises the most voluminous hydrocarbon accumulations in its ‘Pre-salt’ reservoirs. Magmatic rocks have been recorded in the basin mostly by the tholeiitic basalts (Camboriú Formation) within its rift megasequence at Hauterivian time. Recent studies have shown that magmatic rocks are also found interbedded with Aptian sedimentary rocks in the basin, related to its post-rift (mostly Aptian) megasequence. This paper presents petrological data obtained from sidewall core samples of magmatic rocks interbedded with Early Aptian coquina in a well drilled in the Santos Basin. The well comprises six magmatic sections of alkaline basalts being overlain and underlain by the Barra Velha limestones. The basalts are mostly plagioclase-phyric, amygdaloidal rocks but vary in texture and modal compositions. Magmatic compositions seem to have been slightly affected by alteration processes. Basalts comprise an alkaline, miaskitic series with relatively evolved compositions (MgO = 5.55 ± 1.50 wt%; Ni < 220 ppm; Cr < 400 ppm). The patterns in the Primitive Mantle normalized multielement diagram displays the typical Nb peak (0.3 < La/Nb(N)<0.7) of the OIB. Geochemical modelling showed that basalts within five magmatic sections are not cogenetic by differentiation, representing basaltic flows emitted from different vents. Modal batch partial melting model indicates that it is possible, although unlikely, to generate the differences in La/Yb(N) and Zr/Y ratios of the least evolved samples in the studied well by different partial melting degrees (0.5–11%) from the same mantle source. Binary mixing models using lithogeochemical datapoint to a role for the Tristan da Cunha plume tail in the petrogenesis of the alkaline basalts with no clear evidence for a possible role for the local subcontinental lithospheric mantle. The results of the petrogenetic models are consistent with a shallow MOHO and a high-temperature mantle below part of the Santos Basin in Early Aptian time.

Origin of the Paleoproterozoic basaltic dikes from the central and eastern Dharwar Craton and sills and volcanics from the adjoining Cuddapah Basin, southern India

Reverse fractionation modeling considering energy-constrained assimilation-fractional crystallization is performed to estimate primary magma compositions, degree of crustal contamination, pressure–temperature of equilibrium with mantle, and potential temperatures for the origin of the Paleoproterozoic (~ 2.37–1.88 Ga) basaltic dikes in central and eastern Dharwar Craton and sills and volcanics in the adjoining Cuddapah Basin, southern India. Mineral thermobarometry indicates that the dikes crystallized at upper crustal conditions (~ 1–6 kbar/ ~ 1120–1210 °C). Hence, the reverse fractionation calculations are performed at low pressures by adding olivine + plagioclase + clinopyroxene, olivine + plagioclase and only olivine in equilibrium with melt, and simultaneously subtracting an upper crustal partial melt in small steps until the melt is multiply saturated with lherzolite at a high pressure. The results indicate that the basalts are 5–30% contaminated, and their enriched light rare earth element (REE) patterns can be attributed to upper crustal assimilation. The upper crust was pre-heated to 665–808 °C during dike emplacement. The primary magmas of all basalts were last equilibrated with spinel lherzolite at 10–16.5 kbar/1291–1366 °C, and they resemble pooled polybaric incremental melts generated along a ~ 1450 °C adiabat. The estimated mantle potential temperatures (1293–1515 °C) are similar to Paleoproterozoic ambient mantle temperatures. All basalts and their primary magmas show lower chondrite-normalized DyN/YbN ratios than the plume-derived mid-Proterozoic Mackenzie dikes of Canadian Shield, and the primary magmas show flat REE patterns indicating spinel lherzolite melting. The low estimated potential temperatures, low DyN/YbN ratios, and a spinel-bearing mantle source are at odds with an origin of the basalts from mantle plumes.

An assessment of potential causal links between deglaciation and eruption rates at arc volcanoes

One of the fundamental questions that underpins studies of the interactions between the cryosphere and volcanism is: do causal relationships exist between the ice volume on a volcano and its eruption rate? In particular, it is critical to determine whether the decompression of crustal magma systems via deglaciation has resulted in enhanced eruption rates along volcanic arcs in the middle to high latitudes. Evidence for such a feedback mechanism would indicate that ongoing glacier retreat could lead to future increases in eruptive activity. Archives of eruption frequency, size, and style, which can be used to test whether magma generation and eruption dynamics have been affected by local ice volume fluctuations, exist in the preserved eruptive products of Pleistocene-Holocene volcanoes. For this contribution, we have reviewed time-volume-composition trends for 33 volcanoes and volcanic groups in arc settings affected by glaciation, based on published radiometric ages and erupted volumes and/or compositions of edifice-forming products. Of the 33 volcanic systems examined that have geochronological and volumetric data of sufficient resolution to compare to climatic changes since ∼250 ka, increases in apparent eruption rates during post-glacial periods were identified for 4, with unclear trends identified for a further 12. Limitations in the geochronological and eruption volume datasets of the case studies make it difficult to test whether apparent eruption rates are correlated with ice coverage. Major caveats are: 1) the potential for biased preservation and exposure of eruptive materials within certain periods of a volcano’s lifespan; 2) the relative imprecision of geochronological constraints for volcanic products when compared with high-resolution climate proxy records; 3) the reliance on data only from immediately before and after the Last Glacial Termination (∼18 ka), which are rarely compared with trends throughout the Pleistocene to test the reproducibility of eruptive patterns; and 4) the lack of consideration that eruption rates and magma compositions may be influenced by mantle and crustal processes that operate independently of glacial advance/retreat. Addressing these limitations will lead to improvements in the fields of geochronology, paleoclimatology, and eruption forecasting, which could make valuable contributions to the endeavours of mitigating future climate change and volcanic hazards.

Global volcanic arc magma composition correlates with thickness of both arc crust (Moho depth) and arc lithosphere (LAB depth): A revealing message on arc basement histories and arc magmatism

Past fifty years have seen mounting publications on the genesis of volcanic arc magmas. While details remain debated, it is generally agreed that arc magmas result from slab-dehydration induced mantle wedge melting followed by crustal level differentiation of varying extent and sophistication. Two recent arc magma studies deserve particular attention because they attempt to discuss globally unifying controls on arc magma composition. Both Harvard study (Turner and Langmuir, 2015a,b) and Rice study (Farner and Lee, 2017) show correlations of arc magma composition with crustal thickness and both ascribe the crustal thickness as the principal control on their observed magma compositional variations, yet the physical role of the crustal thickness in their interpretations is markedly different because of (1) the ambiguous use of “crust” and (2) their different magma compositional ranges chosen in discussion. The Harvard study only uses basaltic samples corrected to MgO = 6.0 wt.% to discuss mantle processes and interprets the arc crustal thickness as restricting the mantle wedge melting, i.e., the extent of melting decreases with increasing crustal thickness. The Rice study uses samples of all compositions (basaltic to rhyolitic), whose extent of differentiation increases with increasing crustal thickness, interpreted as Moho-crossing mantle wedge melts travelling greater vertical distance with greater degree of cooling and erupting more evolved compositions above thicker crust than melts erupted above thinner crust without need of invoking mantle wedge processes. We commend these efforts and approve their different approaches but emphasize that the unifying understanding of global arc magmatism requires clearly defined Moho (the base of the crust) and LAB (the lithosphere-asthenosphere boundary) and their intrinsic controls on mantle wedge melting (Harvard Study model) and crustal level magma differentiation (Rice Study model) beneath global arcs. In this study, we use chemical compositions of 36,945 global arc volcanic samples provided by the Rice study together with the literature data on seismic Moho and LAB depths of these sample locations to establish (1) the correlation of crustal thickness interval averaged magma composition with crustal thickness and (2) the correlation of lithosphere thickness interval averaged magma composition with lithosphere thickness. These correlations reaffirm our understanding that the lithospheric mantle must exist beneath volcanic arc crust with a globally averaged LAB/Moho depth ratio of 3.26 ± 0.63, i.e., the arc crust is on average about 31.8% ± 6.1% of the lithosphere thickness. This knowledge forms a solid constraint essential for models of global arc magmatism.

Magma evolution in a complex geodynamic setting, South Harghita volcanic area, East-Central Europe: Constraints from magma compositions and zircon petrochronology

Magma evolution in a complex geodynamic setting, South Harghita volcanic area, East-Central Europe: Constraints from magma compositions and zircon petrochronology

Origin and crustal contamination of magmas of the Chiapanecan Volcanic Arc in Southern Mexico: Constraints from new geochemical, isotopic, and geochronologic data

The processes related to the heterogeneous magma compositions of the Chiapanecan Volcanic Arc (CVA) are still poorly understood. Here we present new field descriptions combined with UPb zircon geochronology and whole-rock major, trace-element, and isotopic (Sr, Nd, Pb) data on poorly studied volcanic structures that surround the active El Chichón volcano in the northwestern part of this arc. These consist of the Santa Fe and La Mina intrusive bodies, volcaniclastic fans (La Catedral, and newly discovered Juárez and El Azufre) as well as Miocene volcanics. The mineralogy of the studied rocks consists of plagioclase + orthopyroxene + clinopyroxene + amphibole + biotite, ± FeTi oxides, ± olivine, ± quartz (this latter one only in the Santa Fe intrusive), and zircon as accessory phase. Young populations of zircons from El Azufre and Juárez lahar deposits were dated by UPb at ~2.1 Ma, being among the oldest yet-found magmatic products of the CVA. Strikingly, the Miocene volcanics have similar chemical characteristics as products of the modern arc. Based on geochemical data, we deduce that magma genesis in the CVA takes place in the mantle wedge, where it is constantly modified by aqueous fluids from the slab, causing relatively high Ba/La and Th/Yb ratios. The 87Sr/86Sr and 143Nd/144Nd ratios, zircon geochronology, and trace elements strongly suggest that crustal contamination was the dominant mechanism of differentiation of the magmas, supported also by Ce/Pb and Nb/U ratios. Furthermore, inherited zircons suggest the presence of Grenvillian and Permian (Chiapas Massif) rocks underneath the CVA. The contamination of the magmas by these basement rocks may be responsible for the K2O-enrichment in all CVA rocks.

Molybdenum isotopes record recycling of subducting sediment in active continental margin, Northeast China

Highly variable molybdenum (Mo) isotope compositions have been reported for mafic arc volcanics above oceanic subduction zones. The heavy Mo isotope signatures of mafic arc magmas are generally attributed to addition of high δ98Mo oceanic crust-derived fluids to the subarc mantle sources. However, subduction zone fluids at subarc depths may occur in the forms of aqueous solution and hydrous melts, and seafloor sediments would also be subducted to subarc depths to produce hydrous melts. If such sediment-derived melts were incorporated into the mantle wedge, it is intriguing what is their effect on the Mo isotope composition of mafic arc magmas. Here we report Mo isotope data for Early Jurassic mafic igneous rocks from Northeast China, which was an active continental margin during subduction of the Paleo-Pacific Plate in the Mesozoic. These rocks have δ98Mo values of −0.26‰ to 0.17‰, which are mostly higher than the normal mantle value of −0.20‰ ± 0.01‰. They also exhibit arc-like trace element compositions, moderate (87Sr/86Sr)i ratios of 0.7042 to 0.7047, and neutral εNd(t) values of −2.8 to 2.2. The high δ98Mo values are associated with low Mo/Ce ratios, enriched SrNd isotopes, and high melt-mobile element contents in the mafic igneous rocks. These geochemical features indicate that the heavy Mo isotope signatures in these rocks were mainly inherited from subducting sediment-derived melts. Subducting Paleo-Pacific oceanic crust-derived aqueous solutions and seafloor sediment-derived hydrous melts would metasomatize the overlying mantle wedge beneath the northeastern China continent, and the sediment-derived melts would have contributed the heavy Mo isotope signatures to the mantle sources of mafic arc magmas in the active continental margin. Therefore, the Mo isotope study of mafic igneous rocks is a powerful means to decipher the recycling of crustal components at convergent plate margins.

Origin and Petrogenesis of Magmatism in Collision-Related Environments: Evidence from the Melikler Volcanics on the Kars Plateau-Turkey in the Turkish-Iranian High Plateau

Abstract The temporal distribution of enriched source components and magmatism in continental collision zones provides critical information about mantle dynamic processes in collision-related environments. This paper presents petrology, mineralogy, K-Ar ages and whole-rock major and trace elements, as well as Sr-Nd-Pb-Hf isotopic compositions of Melikler volcanism in Kars Plateau (KP) in the East Anatolia Collision Zone, NE Turkey, with the aim to understand the role of the subducting slab, the origin of magmatism and the geodynamic evolution in the collision-related environments. Our K-Ar dating results show the Melikler volcanism erupted between 5.29 and 1.7 Ma and allows us to divide it into an early (5.29–2.53 Ma) and a late (2.24–1.7 Ma) stage. Major-trace element abundances, isotopic compositions, EC(R) AFC (energy-constrained recharge, assimilation, and fractional crystallisation) and MELTS model calculations of both stages indicate that the least evolved samples were not affected by significant crustal contamination and fractional crystallisation. More evolved samples of the late stage underwent AFC processes with up to 6–9% crustal assimilation; however, those of the early stage were differentiated from a parental magma composition via AFC (up to 2–7.5% crustal assimilation) and experienced magma replenishment at pressure of 0.5 kbar; thus, both early and late stages have experienced open system conditions. The least evolved samples of both stages across the KP have arc-enriched geochemical and isotopic signatures, characterised by prevalent negative Nb–Ta anomalies and moderately radiogenic Sr, unradiogenic Nd-Hf and highly radiogenic Pb isotopic compositions. These primary melts could be derived from a depleted MORB mantle source metasomatised by sediment melt from the subducting Neotethys oceanic slab. Combined trace elemental and isotopic modelling results suggest that the least evolved samples of the early stage were formed by 2–4% melting of an amphibole-bearing garnet lherzolitic mantle source, which was metasomatised by 0.3–0.5% contribution of subducted slab component with a ratio of sediment melt/AOC (altered oceanic crust) melt about 90:10. A depleted lherzolitic mantle source containing apatite and garnet through inputs of 0.6–0.8% melts derived from the subducted oceanic slab, with 5–10% partial melting degree, could produce the least evolved samples of the late stage. Thermobarometric calculations reveal that the least evolved samples of the late stage are derived from the lithosphere-asthenosphere boundary at a depth of 77–82 km; in contrast, those of the early stage are produced from the lithosphere at a depth of 66–69 km. Literature data and the findings obtained from this study indicate that the onset of the Arabian-Eurasian collision may have occurred in the Oligocene and lithospheric dripping caused by the hard collision that occurred around the Late Miocene-Early Pliocene may produce the Melikler volcanic rocks.

Integrating global geochemical volcano rock composition with eruption history datasets

The major element composition of volcanic rocks carries important information about the source and differentiation processes affecting the magma, the physical properties that allow it to erupt, and its eruptive style. Although global rock geochemical databases exist, these are not linked to volcanic eruption history which hampers our global understanding of the relationship between magma composition and eruption dynamics. Here, we integrate two global databases, the Geochemistry of Rocks of the Oceans and Continents (GEOROC) and the Holocene volcanoes of the world of the Global Volcanism Program (VOTW-GVP). The integration is based on matching the location name, geographic position and eruption time, which is automated by a tool called DashVolcano. The tool is open-source, accessible at https://github.com/feog/DashVolcano, and gives access to the integrated datasets via an interactive dashboard. DashVolcano is based on more than 138,000 volcanic rock samples and provides the basis for the identification of global relationships between eruption styles, volcano types, and rock composition for more than 700 volcanoes and their eruptions for the last 10,000 years. The combined record of the eruptive history and its corresponding geochemical rock composition that DashVolcano provides can be used for characterizing global geochemical differences between volcanoes, and should also prove useful for improved long-term hazard and risk evaluations.

A comprehensive re‐look into Jurassic‐Cretaceous Neotethyan active margin in Iran: Evidence of deep mantle flow into subduction magma factory and formation of transitional basalts

Subduction zones display the most diverse magma compositions from various calc‐alkaline differentiates to alkaline mafic associations, but mechanisms generating Nb‐Enriched Basalts (NEB) and High‐Nb Basalts (HNB) in such settings have received less attention from the researchers. Occurrences of NEB and HNB magmas in the Mesozoic Neotethyan active margin in Iran, known as the Sanandaj‐Sirjan Zone, provide an excellent natural lab for examining this interesting petrologic puzzle. The Middle‐Jurassic typical calc‐alkaline arc plutons exhibit a NW‐SE subduction parallel trend. In the northwestern sector of the zone, the Late Jurassic–Early Cretaceous magmatic rocks (JCMR) define a chain that transects the Middle‐Jurassic arc trend. The two cross‐cutting magmatic chains are chemically and genetically different. The JCMR include three distinct compositional groups: (i) tholeiites that are depleted in high‐field‐strength elements (HFSE) and show subduction‐related geochemical fingerprints, (ii) NEB with higher abundances of Nb and other HFSE, and (iii) HNB with the highest abundances of Nb and other HFSE. The latter two groups do not display typical arc signatures and are transitional between arc and plume‐type magmas. These three groups of rocks show linear major and trace element, and Sr–Nd isotopic variation trends extending from Mid‐Ocean Ridge Basalts‐like tholeiites to Oceanic Island Basalts‐like HNB. Hence, a two‐component mixing process is suggested for explaining the compositional variation of the JCMR. We propose that a vertical window resulted from slab tear created the geodynamic environment for juxtaposition of magmas originated at different mantle depths and formation of the transitional basalts by magma mixing.

Two stages of porphyry Cu mineralization at Jiru in the Tibetan collisional orogen: Insights from zircon, apatite, and magmatic sulfides

Continental collision−related porphyry copper (Cu) deposits provide significant global copper resources, but their genesis remains controversial because it is not clear whether remelting of remnant arc-derived lower-crustal Cu-rich cumulates is critical to their formation. We investigated zircon and apatite compositions of the Jiru porphyry Cu deposit in the Gangdese belt of southern Tibet, which is characterized by weaker early Eocene mineralization and more pronounced Miocene mineralization. Our data demonstrate that apatite hosted in zircon can record the volatile compositions of magma before fluid exsolution. Zircon-hosted apatites from the early Eocene granite that have low XF/XCl ratios (≤3) and are interpreted to have crystallized from volatile-undersaturated magma have Cl contents of 0.96−2.14 wt% (1.86 ± 0.38 wt%; n = 15) and SO3 contents of 0.01−0.14 wt% (0.08 ± 0.04 wt%; n = 15). In contrast, zircon-hosted apatites from the Miocene porphyry that have low XF/XCl ratios (≤10) and are interpreted to have crystallized from volatile-undersaturated magma have Cl contents of 0.40−0.56 wt% (0.47 ± 0.06 wt%; n = 11) and SO3 contents of 0.20−0.92 wt% (0.59 ± 0.29 wt%; n = 11). In addition, the early Eocene magma was relatively reduced (ΔFMQ = 0.86 ± 0.55, where FMQ is fayalite-magnetite-quartz) compared with the oxidized Miocene magma (ΔFMQ = 2.04 ± 0.43). The contrasting magmatic oxidation states and apatite/melt S contents were likely critical in controlling the different scales of Cu mineralization during the early Eocene and Miocene, as more highly oxidized magma can dissolve much larger amounts of sulfur and metals (e.g., Cu). Melt Cl content may not play a critical role in magma fertility, since the Jiru early Eocene granite had higher melt Cl contents than the Miocene porphyry. Magmatic sulfides (pyrrhotite, chalcopyrite, and pyrite) were not recognized in the Miocene zircons but did occur in the early Eocene zircons, and these sulfides were demonstrated to have crystallized before fluid exsolution in the shallow magma reservoir. The Jiru early Eocene magmatic sulfide saturation might not have enhanced the amount of copper available to hydrothermal fluids during the early Eocene. Late-stage sulfide saturation in the previous arc magmas might reduce the potential for collision-related porphyry Cu mineralization that is associated with remelting of the previous arc lower crust.

Paleocene‐Eocene High‐Pressure Carbonation of Western Alps Serpentinites: Positive Feedback Between Deformation and CO2‐CH4 Fluid Ingression Responsible for Slab Slicing?

AbstractMeter‐ to hectometer‐size horizons of carbonate‐bearing talcschists are found along or near the contacts between the different Liguro‐Piemont subduction slices. Through mineral and bulk‐rock geochemistry, fluid inclusion analyses, and titanite U‐Pb geochronology, this study shows that these horizons formed by the transformation of serpentinites, at conditions close to peak burial, due to fluid infiltration sourced from surrounding sediments. These rocks containing large amounts of carbonates and high concentrations of fluid‐mobile elements may be significant for carbon cycling and contribute to the composition of arc magmas once the main constitutive assemblage destabilizes at high‐pressure and high‐temperature conditions. Following fluid infiltration, the formation of talc likely controlled strain localization in these horizons, enhancing in turn the ingression of external fluids and creating a positive feedback loop between deformation and fluid infiltration. We suggest that these rocks may have acted as major rheological weaknesses responsible for material offscraping from the downgoing plate.

A peridotite source for strongly alkalic ultrabasic HIMU lavas of the Oslo Rift, Norway

High alkali (Na2O + K2O) low SiO2 lavas generally have highly fractionated trace element compositions and their petrogenesis is debated. Both pyroxene-rich (pyroxenite) and olivine-rich (peridotite) mantle sources have been proposed for such lavas. HIMU (high 238U/204Pb) lavas usually have such compositions and are often interpreted as being derived from garnet-pyroxenite (eclogite) residue of subducted recycled oceanic crust, possibly including a carbonate component. In some cases, an origin from a garnet peridotite source has also been considered.We present new Mg, Pb, Nd, and Sr isotopic data and elemental data for alkalic ultrabasic lavas (melilitites, nephelinites) from the Oslo Rift, Norway. The magmatic province of the Permo-Carboniferous Oslo Rift is interpreted as part of the much larger Skagerrak-Centered Large Igneous Province (SCLIP) which may be related to a mantle plume source originating from the African Large Low Shear Velocity Province (LLSVP) at the core-mantle boundary. The range of δ26Mg is narrow (−0.32 to −0.28) and is similar to typical mantle values, and therefore is inconsistent with a carbonated source as has been suggested for other similar lavas. The 206Pb/204Pb values of all the lavas are high, which identifies them as HIMU lavas. The εNd of +1 to +2 are substantially lower than for typical HIMU lavas, and in the εNd-εSr diagram they plot on the LoNd array between HIMU and EMI. The lavas are more like archetypal kimberlites (formerly, Group I kimberlites), as these also have εNd slightly higher than the chondritic value and Pb isotopes consistent with a HIMU source. The REE patterns of the lavas are also almost as fractionated as in Group I kimberlites.We used major and trace elements to infer the primary magma composition of these lavas. This primary magma composition is compared with experimental melts of both carbonated peridotite and carbonated eclogite source compositions and shows that their major element compositions are only consistent with either a very high degree of melting of eclogite or a very low degree melting of garnet peridotite. The strongly fractionated REE patterns of the lavas can only be made by low-degree melting of either eclogite or garnet peridotite, which therefore effectively rules out eclogite as a major component of the mantle source of these lavas. Further trace element modeling points to a mantle source composition with >97% garnet peridotite that was lightly carbonated (<0.5 wt% CO2) and melted to a very low degree (∼0.4%) to produce these lavas. The modeled source trace element pattern was constrained by the Nd, Sr, and Pb isotopic compositions of the lavas. The measured high Ti contents of these lavas (3.8–6 wt%) are quantitatively shown to be the result of enrichment by both partial melting and fractional crystallization, using experimentally determined DTi for melting of this type of source. We also show that a low degree of melting of a slightly carbonated source is sufficient to result in the likely very high CO2 content (10–15 wt%) of the primary magma.The rather uniform Pb isotope composition of HIMU magmas suggests that the HIMU source is rather uniform in composition, but the Pb isotope composition requires this source to form late in Earth's history. The similarity of these lavas and kimberlites suggests a very deep source, possibly in the LLSVP that gave origin to the SCLIP. This means that LLSVPs do not represent some primordial layer around the core that escaped homogenization of the mantle by the Moon-forming giant impact. Instead, this layer must have formed late in Earth's history. The constraint that this source contains <3% of garnet-pyroxenite (eclogite) residue of recycled basaltic oceanic crust shows that this layer is primarily not a subducted slab graveyard. The process by which the LLSVPs formed, however, remains obscure.