Evolution Of Arc Magmas Research Articles

Two distinct early Paleozoic subduction zones in the eastern Central Asian Orogenic Belt: Evidence of subduction recycling and arc evolution

This paper presents zircon U–Pb dating and Hf isotopic, and whole-rock geochemical and Sr–Nd isotopic data for early Paleozoic igneous rocks as well as zircon U–Pb dating for the sedimentary rocks of the Songnen Massif within the eastern Central Asian Orogenic Belt (CAOB), aiming to reveal arc magma evolution process. The zircon U–Pb dating results of the siltstone of the Xiaojingou Formation yielded 450 Ma, 471 Ma, 502 Ma, 773 Ma, 818 Ma and 949 Ma age peaks, as well as multi-episodic Paleo-proterozoic ages (2486-1612 Ma), implying its deposition time is younger than ∼ 450 Ma. The zircon U–Pb dating results reveal that the Songnen Massif exists three episodes of volcanic rocks of the middle Cambrian (502 Ma), Middle Ordovician (453-449 Ma) and the Silurian (435-420 Ma). The middle Cambrian basalts are low-K tholeiitic to calc-alkaline series, weakly depleted εHf(t) and εNd(t) values. They were derived from the weakly depleted mantle metasomatized by slab fluid in an immature island arc setting. The Middle Ordovician to Silurian volcanic rocks show SiO2 of 49.4–58.1 wt%, plotting into high-K calc-alkaline series, which are also featured by higher Th/Nd and (Hf/Sm)PM ratios, together with moderate depleted εNd(t) values, suggesting their derivation of the depleted mantle wedge metasomatized by sediment-derived melt. In contrast, coeval early Paleozoic volcanic rocks in the eastern Xing’an Massif with the MORB-like εNd(t) values are characterized by an origin of the depleted mantle wedge and the distinct interaction of oceanic crust melt and mantle peridotite. Combined with the variation of calculated crust thickness, we suggested that there was a continental arc in the eastern Songnen Massif, which developed from immaturity to maturity during middle Cambrian to Silurian, while the early Paleozoic island arc in the eastern Xing’an Massif is characterized by ocean ridge subduction process.

Spatial and temporal evolution of magmatic arcs: Geochemical diversity and tectonic controls

We review the main petrologic and geochemical features of magmatic rocks from continental arcs compared to their intra-oceanic equivalents. Both types share a number of characteristics, including their typical HFSE-depleted incompatible element patterns and isotopic Sr, Nd, Pb and Hf features, intermediate between those of N-MORB and continental crust. They are thought to derive from the partial melting of the mantle wedge metasomatized by fluids emanating from the downgoing plate. They are dominantly hydrous fluids provided by the breakdown of amphibolites and antigorite serpentinites, but a significant role is attributed to felsic hydrous melts derived from oceanic metabasalts and associated sediments. The more enriched character of magmas from continental arcs with respect to those from intra-oceanic arcs is ascribed to their genesis from a less depleted mantle wedge in a milder thermal regime, to the enrichment of their source via subduction erosion of the continental margin, and finally to their contamination during ascent and storage within a thicker crust.

Lower crustal control in the iron isotope systematics of plutonic xenoliths from Adak Island, Central Aleutians, with implications for arc magma geochemistry

We present bulk-rock and mineral Fe isotope data of ultramafic to mafic xenoliths and basaltic to andesitic lavas from Adagdak Volcano (Adak Island, Central Aleutians) to study the effects of early differentiation on the Fe isotopic evolution of island arc basalts and their crystallization products. The Fe isotope composition of ultramafic cumulate xenoliths increases from dunite (δ56Fe = –0.09 to –0.02 ‰) to clinopyroxenite (δ56Fe = +0.06 to +0.09 ‰), consistent with higher modal proportions of clinopyroxene (δ56Fe = –0.05 to +0.11 ‰) relative to olivine (δ56Fe = –0.10 to +0.06 ‰) in the latter. Mid-crustal cumulate amphibole gabbro and hornblendite cumulates also record heavier Fe isotope compositions (δ56Fe = +0.04 to +0.08 ‰) due to the abundance of isotopically heavy amphibole (δ56Fe = +0.07 to +0.09 ‰) and magnetite (δ56Fe = +0.11 to +0.13 ‰) in these rocks. High inter-mineral fractionations observed in spinel-olivine and spinel-clinopyroxene pairs (Δ56Fespl-ol = +0.12 to +0.28 and Δ56Fespl-cpx = +0.06 to +0.19) suggest that spinel is not recording equilibrium crystallization conditions for the ultramafic assemblages, likely due to subsolidus Fe-Mg exchange. Our data also include Fe isotope measurements of one mantle dunite (δ56Fe = +0.03 ± 0.05 ‰). Five Adagdak lavas, spanning from basalts to andesites, yield a narrow range of δ56Fe between +0.03 and +0.06 ‰. Our results highlight the potential of amphibole in driving the Fe isotope depletion trends observed in many erupted arc lavas, as amphibole hosts 28–99 % of the FeOT budget in the amphibole gabbro and hornblendite cumulates. This is also supported by single-crystal synchrotron Mössbauer spectroscopy of two amphibole grains, the first from an amphibole gabbro and the second from a hornblendite, which yield Fe3+/ΣFe ratios of 0.55 ± 0.06 and 0.58 ± 0.02, respectively. Water content and hydrogen isotope compositions determined by secondary-ion mass spectrometry from the same amphibole grains indicate partial dehydrogenation. Using Rayleigh fractionation modeling to account for oxidation during post-crystallization dehydrogenation, we calculate magmatic Fe3+/ΣFe ratios of 0.41 ± 0.04 for the amphibole gabbro and 0.30 ± 0.05 for the hornblendite. These data are then used to estimate an appropriate Fe force constant for Adagdak amphibole and quantitatively evaluate the effects of amphibole fractionation. Through a fractional crystallization model, we show how arc melts may experience periods of increasing δ56Fe during olivine-dominated fractionation, followed by decreasing δ56Fe once magnetite and amphibole saturate as cumulate phases. Notably, this dichotomy between fractionation of isotopically light versus heavy cumulate assemblages and its effects on the Fe isotope evolution of arc magmas is not captured by the Adagdak lava record, highlighting the utility of cumulates in chronicling the early isotopic evolution of magmatic systems.

Tectonic and Crustal Processes Drive Multi-Million Year Arc Magma Evolution Leading up to Porphyry Copper Deposit Formation in Central Chile

Abstract Subduction zone magmatism is a major control of volcanism, the generation of modern continental crust and the formation of economically important porphyry Cu–(Mo–Au) deposits. Reading the magmatic record of individual arc segments and constraining the rates of magmatic changes are critical in order to fully understand and quantify the processes that drive magma evolution in subduction settings during arc growth. This study focuses on the San Francisco Batholith and the Rio Blanco-Los Bronces porphyry deposit cluster in central Chile, which provides an igneous rock record over ~13.5 Myr of arc evolution. We use whole-rock geochemistry, zircon geochronology and Hf isotope geochemistry to track changes in the crustal magmatic system of this arc segment during crustal thickening and porphyry Cu deposit formation. By combining the analytical dataset with Monte Carlo fractional crystallisation and assimilation fractional crystallisation modelling, we test a model for significant crustal involvement during magma evolution. Systematic and continuous increases in Dy/Yb, La/Yb, V/Sc and Sr/Y in the magmas over time indicate a transition in the main fractionation assemblage from plagioclase-dominated to amphibole-dominated that reflects deeper crystallisation and/or a higher meltwater content. Concomitant decreases in εHf and Th/La as well as increasing Ba/Th are best explained by assimilation of progressively deeper crustal lithologies from low (Chilenia) to high Ba/Th (Cuyania) basement terranes. Our study highlights that an increasingly hydrous magma and a deepening locus of crustal magma differentiation and assimilation, driven by crustal thickening contemporaneous with increased tectonic convergence and ingression of the aseismic Juan Fernandez ridge, can account for all investigated aspects of the multi-Myr magmatic evolution leading up to the formation of the Rio Blanco-Los Bronces porphyry Cu deposits. Our findings corroborate the importance of high-pressure differentiation of hydrous magma for the formation of Andean-style porphyry deposits. Once magmas favourable for porphyry Cu mineralisation were generated in the lower crust, multiple episodes of efficient magma migration into the upper crust fed several, discrete, shallow magmatic-hydrothermal systems over ~3.5 Myr to form the world’s largest known Cu resource at Rio Blanco-Los Bronces.

High-Mg andesites and Nb-enriched basalts in the Dulate arc, East Junggar (NW China): Evidence of ocean ridge subduction

Abstract The Dulate arc, located in East Junggar (NW China) in the southern Central Asian orogenic belt, records a Devonian magmatic arc evolution, offering a window to understanding the orogenic processes of the Central Asian orogenic belt. Here we present new geochemical and isotopic data for Late Devonian high-Mg andesite (HMA) and Nb-enriched basalt (NEB) suites from the Qiakuerte area, East Junggar. The HMA samples are typical subduction-related volcanic rocks. They have SiO2 contents ranging from 53.30 to 54.59 wt%, high MgO (5.0–5.26 wt%), and high Mg# values (~55) and show enrichments in large ion lithophile elements (LILEs) and depletions in high field strength elements (HFSEs). The HMA samples have high (La/Yb)N ratios and Sr/Y (~6.5 and 50–59, respectively) with no Eu anomalies. The HMA samples have high Na2O (~3.3 wt%) and low K2O (~2.5 wt%) and Th (~2.4 ppm) contents, combined with positive εNd(t) and low (87Sr/86Sr)i values. These characteristics suggest that the samples were formed mainly through interactions between subducted oceanic melts and mantle peridotites. Compared to normal arc basalts, the NEB samples have higher concentrations of Nb (~20 ppm), higher primitive mantle–normalized Nb/La (0.50–0.58), and higher ratios of Nb/U (9.4–14.6). The NEB samples also have positive εNd(t) and low (87Sr/86Sr)i values, indicating that their source was mantle wedge that had been metasomatized by slab melt. Considering the widespread presence of A-type granites, the abnormally high heat flow, and the tectonic characteristics of East Junggar, we conclude that a slab window created by the subduction of an ocean ridge was responsible for the melting of slab and the formation of the NEB-HMA suites. These processes may have also played a key role in the tectonic evolution processes of East Junggar during the Late Devonian.

Detrital isotopic record of a retreating accretionary orogen: An example from the Patagonian Andes

Abstract U-Pb zircon geochronology and isotopic records have played an influential role in our understanding of convergent margin dynamics. Orogenic cyclicity models link tectonic regimes with magmatic isotopic signatures in advancing orogens, relating compressional regimes with evolved signatures and extension with juvenile signatures; however, such frameworks may not apply for retreating orogens, which commonly produce substantial crustal heterogeneities during backarc rifting and ocean spreading. We explore the Mesozoic to Cenozoic Patagonian Andes tectonic evolution, combining U-Pb zircon ages, bulk rock εNd, and new detrital zircon εHf from the retroarc basin to understand the associated magmatic arc evolution during retreat and advance of the margin. Our results reveal a protracted phase of isotopically juvenile magmatism between 150 and 80 Ma, which began during backarc extension and persisted long after the margin switched to a contractional regime. We propose that the prolonged juvenile isotopic trend started mainly due to trenchward migration of the arc during backarc extension (150–120 Ma) and persisted due to partial melting of underthrusted juvenile attenuated and oceanic crust during backarc basin closure (120–80 Ma). This interpretation implies that tectonic stress alone does not predict isotopic trends, and factors like assimilation or the composition of underthrusted crust are important controls on magmatic isotopic composition, especially in retreating and transitional orogens.

Alkaline magmas in shallow arc plutonic roots: a field and experimental investigation of hydrous cumulate melting in the southern Adamello batholith

Despite the first-order importance of crystallisation–differentiation for arc magma evolution, several other processes contribute to their compositional diversity. Among them is the remelting of partly crystallised magmas, also known as cumulate melting or ‘petrological cannibalism’. The impact of this process on the plutonic record is poorly constrained. We investigate a nepheline-normative dyke suite close to the Blumone gabbros, a large amphibole-gabbro unit of the Tertiary Southern Alpine Adamello igneous complex. The compositions of the studied dykes are characterised by low SiO2 (43–46 wt. %), MgO (5.0–7.2 wt. %), Ni (18–40 μg/g), and high Al2O3 (20.2–22.0 wt. %) contents. Phenocrystic plagioclase in these dykes exhibits major, trace, and Sr isotope compositions similar to Blumone cumulate plagioclase, suggesting a genetic link between the nepheline-normative dykes and the amphibole-gabbro cumulates. We tested this hypothesis by performing saturation experiments on a nepheline-normative dyke composition in an externally heated pressure vessel at 200 MPa between 975 and 1100 °C at fO2 conditions close to the Ni–NiO buffer. Plagioclase and spinel are near-liquidus phases at and above 1050 °C, contrasting with the typical near-liquidus olivine ± spinel assemblage in hydrous calc-alkaline basalts. The alkaline nature of the dykes results from the abundance of amphibole in the protolith, consistent with melting of amphibole-gabbro cumulates. We modelled the heat budget from the repeated injection of basaltic andesite into a partly crystallised amphibole-gabbro cumulate. The results of this model show that no more than 7% of the cumulate pile reaches temperatures high enough to produce nepheline-normative melts. We propose that such nepheline-normative dykes are a hallmark of hydrous cumulate melting in subvolcanic plumbing systems. Therefore, ne-normative dykes in arc batholiths may indicate periods with high magma fluxes.

Magmatic arc evolution during the tectonic closure of the Rocas Verdes basin: insights from Cretaceous–earliest Paleocene intrusive rocks of Navarino Island (55°S), Fuegian Andes

Navarino Island is located in the southernmost part of the Fuegian Andes, south of the Beagle Channel. Its geological record documents the complex tectonic history of Tierra del Fuego, which includes the opening and closure of the Rocas Verdes basin, Cordillera arc collision and subsequent subduction processes. The geology of the island is mostly composed of the Cretaceous Yahgán Formation, a marine metasedimentary sequence, which is intruded by diverse plutons that are mostly exposed on the northwestern tip of the island. We herein present a new dataset that shows the presence of three Cretaceous–earliest Paleocene magmatic suites of active margin magmatism emplaced during the early stage of the Fuegian Andes, which are referred to as (1) the Dientes de Navarino Microdioritic Sills, a suite of pre-tectonic microdioritic sills that formed at c. 101–97 Ma; (2) the Castores Plutonic Complex, a series of pre- to syntectonic gabbroic to tonalitic plutons emplaced at c. 90–87 Ma, and (3) the Samantha Monzonites, a suite of isolated monzonitic to monzodioritic post-tectonic plutons that formed at c. 66–65 Ma. These distinct magmatic episodes are recognized by field observations, geological mapping, petrography and whole-rock geochemistry integrated with amphibole and biotite 40 Ar/ 39 Ar and U–Pb LA-ICP-MS. The geochemical compositions of these rocks are consistent with a continental arc setting that existed during the interval c. 101–65 Ma. Whereas the three pulses spatially overlap in Navarino Island, the arc magmatism shows a migration (or expansion) throughout the Late Cretaceous. The locus of the arc then migrates at c. 68–66 Ma towards the SW. We suggest that this trenchward migration at c. 68–66 Ma may be associated with a change in the subduction angle. The three Cretaceous–earliest Paleocene plutonic pulses recorded in Navarino Island formed during the early stages of development of the Fuegian Andes, and are pre-, syn- and post-tectonic with respect to a major compressional event caused by the collision and obduction of the back-arc Rocas Verdes oceanic floor. Supplementary material: A detailed dataset is available at https://doi.org/10.6084/m9.figshare.c.6675416

Phosphorus deficit in continental crust induced by recycling of apatite-bearing cumulates

Abstract Trace-element anomalies in the continental crust (e.g., Nb-Ta-Pb) are useful to decipher the formation and evolution of continents. We found that the trace-element patterns of upper, middle, and lower crust are all marked by strongly negative phosphorus anomalies, which were not likely to have been generated directly by mantle melting. In order to address this “crust composition paradox,” we compiled major- and trace-element data from global arc magmas, igneous minerals, and deep arc cumulates and found that the phosphorus deficit in continental crust is coupled with an elevated Th/La ratio. This feature can be best explained by the fractional crystallization of apatite during arc magma evolution rather than the influence of a subducted slab. Accumulation of apatite-bearing mafic cumulates in deep arcs followed by foundering into the upper mantle causes the mismatch in the phosphorus systematics between the mantle and the continental crust. This process plays a critical role in shaping the chemical composition of continental crust.

A crystal mush origin for the high Sr/Y plutons: Insights from the Nyemo Igneous Complex in the Gangdese arc, southern Tibet

High Sr/Y plutons generated in continental arc resemble geochemically the Archean trondhjemite-tonalite-granodiorite suite and are widely used to interpret processes of intermediate-felsic crust generation during the Archean. However, this geochemical signature could be overinterpreted without detailed petrographic observation. Here we carry out a study of detailed field and petrographic observations, geochronology , and geochemistry for the Nyemo Igneous Complex (NIC) in the central Gangdese arc, southern Tibet . Our results show that the NIC is composed of gabbros , monzogabbros, monzodiorites, monzonites , and monzogranites that were formed at 91–89 Ma. Petrographic observations and geochemical data suggest that the diverse lithologies of the NIC can be ascribed to a cogenetic magma fractionation sequence defined by olivine → orthopyroxene + clinopyroxene + apatite + Fe-Ti oxides → plagioclase + biotite → K-feldspar → amphibole . The early crystallization of biotite and relatively late crystallization of amphibole indicates that the primary melts of the NIC have low water contents. The monzodiorites and monzogranites from the NIC have high Sr (412–1026 ppm), low Y (4.30–25.4 ppm), and thus high Sr/Y ratios (16–140). The orthocumulate texture defined by euhedral phenocrysts of normally zoned plagioclase and subhedral biotite indicates these high Sr/Y rocks are unlikely pure melts and instead represent fossil plagioclase-rich magma mush. We emphasize that the high Sr/Y signature does not necessarily reflect the pressures of magma generation involved of garnet/amphibole as residual or fractional phases. Late Cretaceous arc magmatic rocks in the eastern and central Gangdese arc show different crystallization sequences, requiring the presence of contemporaneous primary melts with distinct H 2 O contents along the Gangdese arc, which probably were caused by variable water input from the Neo-Tethyan oceanic lithosphere . • Nyemo Igneous Complex (NIC) in the Gangdese arc represents a fractionation sequence. • The primary melts of the Nyemo Igneous Complex have low water contents. • High Sr/Y plutons in the NIC represent fossil plagioclase-rich magma mush. • Arc magma evolution is controlled by the water input from the oceanic lithosphere.

A Role for Crustal Assimilation in the Formation of Copper-Rich Reservoirs at the Base of Continental Arcs

Abstract Understanding the behavior of chalcophile elements during the evolution of arc magmas is critical to refining models for the formation and distribution of porphyry copper deposits used in mineral exploration. Because magmas in continental arcs undergo copper depletion during their early differentiation, a widely held hypothesis posits that the removed copper is locked at the base of the crust in copper-rich cumulates that form due to early sulfide saturation. Testing this hypothesis requires direct evidence for such copper-rich reservoirs and a comprehensive understanding of the mechanisms driving sulfide saturation. Interaction between oxidized magmas and reducing crustal material in island arcs has been shown to be an efficient process causing sulfide saturation. However, the extent to which crustal assimilation impacts the flux of chalcophile elements during magmatism in thick continental arcs remains to be established. Here, we provide a deep perspective into these problems by studying a suite of subarc cumulate rocks from the Acadian orogen, New England (USA). These cumulates record the imprint of subduction zone magmatism and represent the residues left behind during the genesis of intermediate to evolved Acadian magmas (ca. 410 Ma). We find that the most primitive Acadian cumulates are enriched in copper (up to ~730 µg g–1) hosted by sulfide phases, providing direct evidence for the formation of lower crustal copper-rich reservoirs. The Acadian cumulates reveal a wide range of δ34S values, from –4.9‰ in the ultramafic rocks to 8‰ in the most evolved mafic rocks. The negative δ34S values observed in the most primitive and copper-rich cumulates (avg –3‰) reflect the assimilation of isotopically light sulfur from surrounding sulfidic and graphite-bearing metasedimentary rocks (δ34S of –19 to –12‰), whereas the more evolved cumulates with positive δ34S signatures may have formed from different magma batches that experienced less sediment assimilation. The assimilation of these reducing metasedimentary rocks caused a critical drop in oxygen fugacity (~DFMQ –2.5 to –1.9; FMQ = fayalite-quartz-magnetite buffer) in the evolving magmas, ultimately leading to extensive sulfide saturation and the consequent formation of copper-rich subarc cumulates. Assimilation-driven sulfide saturation may be a common process at the root of thickened arc crusts that triggers the formation of lower crustal copper-rich reservoirs, which play a pivotal role in the fate of copper during arc magmatism. Thus, deeply buried reducing metasedimentary crustal material at the base of continental arcs can act as a barrier to the magmatic flux of chalcophile elements and may play a crucial role in the genesis and distribution of porphyry copper deposits.

Amphibole fractionation and its potential redox effect on arc crust: Evidence from the Kohistan arc cumulates

Abstract Arc magmas, a major contributor to continental crust growth, are thought to be more oxidized than mid-ocean ridge basalts as reflected by an enrichment in ferric iron relative to ferrous iron. But how arc magmas become oxidized is hotly debated. It is acknowledged that the fractionation of common Fe-rich phases (e.g., amphibole) may change the Fe valence of the derivative melt. Amphibole has Fe-rich compositions commonly found in arc systems. We present high-precision (±0.01%) Fe valence data of amphibole and cumulates from the Kohistan arc determined by Mössbauer spectroscopy and bulk cumulate Fe3+/ΣFe ratios by wet chemistry. We evaluate the Fe3+/ΣFe trend of Fe-rich amphibole during arc magma fractionation. Our results show that bulk Fe3+/ΣFe ratio of cumulates in mature island arc settings is mainly controlled by amphibole due to its abundance and decreases (from 0.4 to 0.2) with decreasing Mg# as Fe3+/ΣFe ratios decrease in amphibole (from 0.35 to 0.2). Our modeling suggests that amphibole fractionation from parental arc magmas to Fe3+/ΣFe above 0.2 may lead to an increase in Fe3+/ΣFe ratios of middle–late stage residual melts by 0.1–0.3, and the partitioning of ferrous Fe into amphibole increases with the evolution of arc magmas. Our findings highlight the importance of petrological processes in the magma that contribute to the production of fertile arc crust.

Cumulate Mush Hybridization by Melt Invasion: Evidence from Compositionally Diverse Amphiboles in Ultramafic–Mafic Arc Cumulates within the Eastern Gangdese Batholith, Southern Tibet

Abstract Amphibole plays an important role in the petrogenesis and evolution of arc magmas, but its role is not completely understood yet. Here, a field, petrological, geochronological and geochemical study is carried out on ultramafic–mafic arc cumulates with textural and chemical heterogeneities and on associated host diorites from the eastern Gangdese Batholith, southern Tibet to explore the problem. The cumulates occur as a large body in diorite host-rocks. The core of the body consists of coarse-grained Cpx hornblendite with a porphyritic texture. Towards the contact with the host diorite, the coarse-grained Cpx hornblendite grades to relatively homogeneous fine-grained melagabbro. Zircon U–Pb dating indicates that they all crystallized at 200 ± 1 Ma. Textural features and whole-rock and mineral chemical data reveal that both the Cpx hornblendite and the melagabbro are mixtures of two different mineral assemblages that are not in equilibrium: (1) brown amphibole and its clinopyroxene inclusions; (2) matrix clinopyroxene + green amphibole + plagioclase + quartz + accessory phases. Clinopyroxene and brown amphibole from the first assemblage are enriched in middle rare earth elements (MREE) relative to light REE (LREE) and heavy REE (HREE), and are weakly depleted in Ti, whereas clinopyroxene and green amphibole from the second assemblage are characterized by LREE enrichment over MREE–HREE and more marked Sr and Ti depletion. The higher Mg#, MgO and Cr of the late-formed green amphibole than the early formed brown amphibole suggest that the two assemblages are not on the same liquid line of descent. Given the close relations of the three rock types in the exposed crustal section, the cumulates are interpreted to have formed in an open system, in which an ultramafic cumulate body consisting of the first assemblage reacted with the host dioritic melt to form new clinopyroxene and amphibole of the second assemblage. The melt calculated to be in equilibrium with the first mineral assemblage resembles an average continental arc basalt, which is less evolved than the host dioritic melt, responsible for the second mineral assemblage. On the basis of whole-rock Sr–Nd–Hf isotopic similarity of the cumulates and a host diorite sample, we argue that the host diorites were formed through crystal fractionation from the parent melt of the first assemblage. Results of least-squares mass-balance calculations suggest that the quantities of the host dioritic melts, involved in the generation of these modified cumulates, vary from ~25 to ~44 %. The presence of magmatic epidote in the host diorites and Al-in-Hb geobarometry indicate that the reaction that occurred when the dioritic melts percolated through the cumulate body was at ~6 kbar. Both the brown and green amphiboles are enriched in MREE relative to HREE, and can impart to residual melts a strong geochemical signature of amphibole fractionation (low Dy/Yb). Thus, we conclude that fractional crystallization and melt–rock reaction are two mechanisms by which amphibole controls arc magma petrogenesis and evolution.

Evidence of late Palaeozoic and Middle Triassic magmatism in the Sakar-Strandzha Zone, SE Bulgaria: Regional geodynamic implications

ABSTRACT Late Palaeozoic granitoids, meta-granitoids and meta-volcanic rocks predominate in the metamorphic basement of the northern and western parts of the Sakar-Strandzha Zone (SASTZ) in southeast Bulgaria, together with the subsidiary Triassic meta-granitoids and meta-volcanic rocks. Generally, igneous minerals and textures are preserved, except in the meta-granitoids and meta-volcanic rocks that experienced a low- to high-grade metamorphic overprint. The volcanic rocks have a peraluminous and high-K calc-alkaline composition, and the granitoids range between I- to S-type compositions, typical of volcanic arcs and syn-collisional settings. LILE and LREE-enrichment and Nb-Ta anomalies characterize the intrusive and extrusive rock suites. U-Pb zircon geochronology has yielded crystallization ages between 245 and 237 Ma for the majority of the studied igneous rocks, and between 297 and 295 Ma for a small group of igneous rocks. Early Permian and Middle Triassic igneous suites of the northern and western SASTZ have similar compositions and a similar tectonic setting when compared to Late Carboniferous-Early Permian intrusive and extrusive suites of the adjacent Sakar unit of the SASTZ, confirming a common regional late Palaeozoic-early Mesozoic tectono-magmatic event. As the Late Carboniferous-Permian to Middle Triassic magmatic arc components extend across the SASTZ, they trace the time-correspondent active continental margin along the Eurasian plate during subduction of the Paleotethys oceanic lithosphere. The late Palaeozoic Eurasian active continental margin magmatic arc evolution of the SASTZ can be linked with the Serbo-Macedonian-Rhodope zones to the southwest, where coeval meta-granitoids document the same geodynamic context. By contrast, the Triassic igneous suite of the SASTZ is unrelated to the Serbo-Macedonian-Rhodope zones, where Triassic meta-ophiolite and meta-granitoids record Neotethys rifting.

Genesis of Intermediate and Silicic Arc Magmas Constrained by Nb/Ta Fractionation

AbstractResolving the geochemical discrepancies in Nb/Ta between the bulk continental crust and its parental mantle melts can provide insights into the processes by which basaltic magmas differentiated into the felsic crust. Deep‐seated garnet pyroxenite cumulates (i.e., arclogite) have high Nb/Ta values due to the presence of high Nb/Ta‐rutile. For this reason, rutile‐bearing arclogites having Nb‐rich composition were proposed to complement the low Nb/Ta nature of the continental crust and evolved arc magmas, that is to say, continental arc magma evolution dominantly occurred in the stability field of garnet. However, this suggestion is challenged by the incomplete understanding of the role of rutile in the Nb/Ta fractionation during arc magma evolution due to insufficient partitioning data at low temperatures (T ≤ 1,000°C). Here, we determined Nb and Ta partitioning ratios in rutile‐felsic melt systems at 2.0 GPa, 850–1,000°C and ∼4–20 wt.% H2O. Our results confirm that equilibrium crystallization of rutile with respect to Nb and Ta will not impart lower Nb/Ta values to any derivative liquids, whereas diffusive fractionation of Nb/Ta in siliceous melts during the growth of rutile can reduce Nb/Ta ratios in melts. However, a positive correlation between Nb/Ta and Dy/Yb for mature continental arc magmas, incipient continental arc magmas, island arc magmas, and global arc magmas indicates that amphibole rather than rutile‐bearing arclogite was responsible for causing low Nb/Ta ratios of many evolved intermediate rocks. Accordingly, the chemical differentiation from primary basaltic to average andesitic compositions of the continental crust mainly occurred in the stability field of amphibole.

Mafic Microgranular Enclaves Formed by Gas-driven Filter Pressing During Rapid Cooling: an Example from the Gangdese Batholith in Southern Tibet

AbstractMafic microgranular enclaves (MMEs), widespread in intermediate to felsic arc plutons, carry significant information on the genesis and evolution of arc magmas, yet their origin remains debatable. Here, we examine MME-host diorite pairs from the c.200 Ma Cuijiu Igneous Complex in the eastern Gangdese Batholith, southern Tibet, to constrain the petrogenesis of MMEs and the evolution of arc magmas. Within the complex, MMEs are essentially similar to their host diorites with similar emplacement ages (∼200 Ma), mineral assemblages and mineral compositions, as well as whole-rock Sr–Nd–Hf and zircon Hf isotopic compositions. However, MMEs have higher modal contents of hornblende and biotite, and are enriched in compatible elements and depleted in incompatible elements. Zircons from some MME samples are characterized by dark cathodoluminescence (CL) cores overgrown by light-CL rims of varying thickness. The dark-CL cores show higher Th, U and rare earth elements (REE) abundances than the light-CL rims. Based on comparison with co-genetic mafic melts and mass-balance calculations, we propose that the MMEs were early-crystallized cumulates (autoliths) related to their host diorites. The chilled textures, flow microstructures and pillow shapes suggest that the MMEs experienced rapid cooling before being captured by the host magmas. The rapid cooling may result from contact between ascending diorite magmas and cooler wall rocks. As the magmas quickly crystallized, they reached second boiling and vesiculation, and separated into fine-grained crystal-rich margins and melt-rich centres. Gradients in crystallinity and pressure expelled interstitial melts from the crystal-rich margins to the crystal-poor centres, leading to crystal-liquid separation (gas-driven filter pressing). The dark-CL zircon cores with high Th and U abundances may crystallize from highly evolved interstitial melts within the crystal-rich margins. The fine-grained crystal-rich margins were subsequently captured and dragged as MMEs before their complete crystallization by later ascending host magmas. This differentiation process could have occurred over several kilometres of magma ascent, and have played an important role in the polybaric fractional crystallization of the Cuijiu Igneous Complex, feeding more differentiated andesitic magmas to upper crustal mushes.

Subduction zone processes and crustal growth mechanisms at Pacific Rim convergent margins: modern and ancient analogues

AbstractContinents grow mainly through magmatism, relamination, accretionary prism development, sediment underplating, tectonic accretion of seamounts, oceanic plateaus and oceanic lithosphere, and collisions of island arcs at convergent margins. The modern Pacific–Rim subduction zone environments present a natural laboratory to examine the nature of these processes. The papers in this special issue focus on the: (1) modern and ancient accretionary margins of Japan; (2) arc–continent collision zone in the Taiwan orogenic belt; (3) accreting versus non-accreting convergent margins of the Americas; and (4) several examples of ancient convergent margins of East Asia. Subduction erosion and sediment underplating are important processes, affecting the melt evolution of arc magmas by giving them special crustal isotopic characteristics. Oblique arc–continent collisions cause strong deformation partitioning that results in orogen-parallel extension, crustal exhumation and wrench faulting in the hinterland, and thrust faulting–folding in the foreland. Trench-parallel widths of subducting slabs exert major control on slab geometries, the degree of coupling–decoupling between the lower and upper plates, and subduction velocity partitioning. An initially large width of the subducting Palaeo-Pacific Plate against East Asia caused flat subduction and resistance to slab rollback during the Triassic Period. These conditions resulted in shortening across SE China. Foundering and delamination of the flat slab during the Early Jurassic Epoch led to slab segmentation and reduced slab widths, followed by slab steepening and rollback. This pull-away tectonics induced lithospheric extension and magmatism in SE China during Late Jurassic – Cretaceous time. Melting of subducted carbonaceous sediments commonly produces networks of silicate veins in CLM that may subsequently undergo partial melting, producing ultrapotassic magmas.

Noble gas magmatic signature of the Andean Northern Volcanic Zone from fluid inclusions in minerals

Trace volatile elements like He are key for understanding the mantle source signature of magmas and to better constrain the relative roles of subduction and crustal processes to the variability of along-arc chemical and isotopic signatures of magmatic fluids. Here we report on noble gas abundances and isotopic data of Fluid Inclusions (FIs) in eruptive products and/or fumarolic gases from the Colombia-Ecuador segment of Andean Northern Volcanic Zone (NVZ). FIs in olivine phenocrysts from Ecuador (El Reventador, Cotopaxi and Tungurahua) yield air-normalized corrected 3He/4He ratios of 7.0–7.4 RA, within the MORB range (8 ± 1 RA). With exception of the Cotopaxi lavas (opx < <oliv.), these are indistinguishable of those obtained for their cogenetic orthopyroxene pairs and of gas emissions previously reported in literature. Olivine phenocrysts from Nevado del Ruiz fissure lavas also yield the highest 3He/4He (8.5 ± 0.3 RA) for this volcanic system, which is in the range of fumarolic gases for Galeras (previously reported as high as 8.8 RA and here measured to a maximum of 8.3 ± 0.1 RA). Our dataset highlights disparities between isotope signatures of eruptive products from Ecuador (avg. ~7.2 RA) and those reported for the Colombian portion of the NVZ (avg. ~8.5 RA). Previous studies on the geochemistry of erupted products put in evidence significant along-arc variations ascribed either to the involvement of different slab components, or to variable depths of evolution of arc magmas within the continental crust. However, the same variation is not discernible in the signature of noble gases, especially helium, from FIs and gas emissions analyzed in this study, with little inter-variation between Cotopaxi, Reventador and Tungurahua (all within 0.2 RA from the Ecuador average of 7.2) and Galeras and Nevado del Ruiz, whose maximum values differ by ~0.3 RA. We therefore suggest a homogenous MORB-like 3He/4He signature for the mantle wedge beneath this arc segment, whereby along-arc variations in crustal thickness (from <35 km at the northernmost part of the segment to ≥50 km at the Ecuadorian arc segment) may factor largely into the variability recorded on our data set. The first CO2/3He ratios obtained in FIs from Andean rocks support the hypothesis of increasing crustal contamination from Colombia to Ecuador, concomitant with increasing crustal thicknesses under the respective arc regions.

Redox processes in subduction zones: Progress and prospect

Oxygen fugacity (fO2) is an intensive variable that describes the redox state of a system. By controlling the valence state of multivalent elements, fO2 affects the stability of iron-bearing minerals, dominants the species of volatile elements (e.g., carbon and sulfur), and controls the partitioning behaviors of multivalent elements (e.g., iron, vanadium, cerium, europium). Thus, fO2 plays a key role in understanding the generation and differentiation of arc magmas, the formation of magmatic-hydrothermal deposits, and the nature of magmatic volatiles. Subduction zones are an important site for arc magmatism and fluid action, and the study of redox processes is indispensable in subduction zone geochemistry. In this paper, we first introduce the concept, expression, and estimation methods of fO2. Then we retrospect the history and progress about the oxidation state of the metasomatized mantle wedge, summarize the redox property of slab-derived fluids, and review the latest progress on redox evolution of arc magmas during magma generation and differentiation. The main conclusions include: (1) despite its wide variation range, fO2 of the mantle wedge is generally higher than that of the oceanic mantle; (2) the redox property of the subducting slab-derived fluids is still controversial and the mechanism for the oxidization of the mantle wedge remains unclear; (3) how the fO2 varies during the generation and differentiation of the arc magmas is debated. We propose that the crux in deciphering the oxidization mechanism of the mantle wedge is to determine the mobility of iron, carbon and sulfur in subducting slab-derived fluids (especially solute-rich fluid or supercritical fluid); the key in understanding the redox evolution during arc magma generation and differentiation is to determine the partition coefficients of Fe3+ and Fe2+ between ferromagnesian minerals and silicate melts.

Pre-eruptive storage, evolution, and ascent timescales of a high-Mg basaltic andesite in the southern Cascade Arc

Crustal storage and transit play a critical role in the compositional evolution of arc magmas; however, the enigmatic nature of lower crustal magma storage and early differentiation limit our understanding of the connections between the physical processes of subduction zones and the architecture of the arc crust. We present new geochemical data and applications of existing barometric and chronometric tools to interrogate the mantle source compositions, crustal storage depths, and ascent timescales of a primitive, high-Mg, clinopyroxene-bearing cinder cone, the basaltic andesite of Box Canyon, located in the vicinity of Lassen Volcanic National Park, CA, in the southern Cascades. Petrographic examination in addition to bulk and in situ geochemical analyses (XRF, LA-ICP-MS, and EPMA) of tephra and lava-derived samples reveals co-crystallization of clinopyroxene and olivine as phenocrysts and glomerocrysts with ~ Mg# 80 and > Fo85, respectively. Phase equilibria experiments of analogous Cascade Arc magma compositions estimate crustal storage of the observed phase assemblage at pressures in the lower crust > 700 MPa. Reverse zonation in olivine and clinopyroxene phenocryst interiors from core-rim analytical profiles record a lower crustal mafic mixing event. Results from one-dimensional, multi-elemental olivine and clinopyroxene diffusion models fit to these interior mixing zones provide an assessment of available trace element (Ti, La, Yb, Ce, and Nd) diffusion chronometers in clinopyroxene by considering multiple element profiles within two phases that experienced the same pre-eruptive conditions. While multi-elemental and multi-phase diffusion timescales span two to three orders of magnitude, Ni in olivine profiles provide the most-robust estimate of 19.1 ± 8.6 years from mixing-to-eruption. Our results provide new constraints on arc crustal differentiation processes, indicate rapid crustal transit timescales in agreement with a growing database of diffusion-based ascent timescales of basaltic magmas, and demonstrate significant, systematic deviation of diffusion timescales calculated from natural zonation of rare earth elements in clinopyroxene and Ti in clinopyroxene and olivine.