Crustal thickness control on Sr/Y signatures of recent arc magmas: an Earth scale perspective

High Mg# of the continental crust explained by calc-alkaline differentiation

We used compiled geochemical data to investigate the mechanisms that control Mg# (molar ratio of Mg/(Mg + Fe_T)) in andesitic arc lavas. We find that andesites from mature continental arcs with crustal thickness of >45 km have systematically higher Mg# than those from oceanic arcs with crustal thickness of <30 km. The elevated Mg# in continental arc lavas results from strong Fe depletion during high-pressure differentiation favored in thick crusts. This proposal is reinforced by our compiled melting/crystallization experiment data. We show that the Mg# characteristics of continental arc lavas match that of the continental crust. These findings suggest that the formation of many high-Mg# andesites and the continental crust may not require slab-melt/peridotite interactions. Instead, the high Mg# of the continental crust can be explained by intracrustal calc-alkaline differentiation processes in magmatic orogens.

A rapid change in magma plumbing taps porphyry copper deposit-forming magmas

Porphyry-type deposits are a vital source of green technology metals such as copper and molybdenum. They typically form in subduction-related settings from large, long-lived magmatic systems. The most widely accepted model for their formation requires that mantle-derived magmas undergo an increase in volatiles and ore-forming constituents in mid- to lower crustal reservoirs over millions of years, however, this is mostly based on observations from shallow, sporadically exposed parts of porphyry systems. To examine this paradigm, we have evaluated the timeframe and geochemical signatures of magmatism in a ~ 8 km palaeodepth cross-section through plutonic and volcanic rocks of the classic Yerington magmatic system, Nevada. We show that the magmas in the upper parts of the system (< 8 km) underwent a major and rapid change in chemistry over a period of < 200 kyrs that is coincident with the initiation of ore formation. We attribute this change to a shift from extraction of quartz monzodiorite and quartz monzonite magmas evolving in mid-crustal reservoirs, and that had relatively poor ore-forming potential, to extraction of volatile-rich granitic magmas from greater (~ 30 km) depths. As the granites crystallised, late stage melts were intruded through the carapace as aplite dykes which contain traceable expressions of the porphyry deposit-forming fluids. The rapid nature of the shift in ore-forming potential narrows the temporal-geochemical footprint of magmas associated with porphyry mineralisation and provides new constraints for exploration models.

Chemical Mohometry: Assessing Crustal Thickness of Ancient Orogens Using Geochemical and Isotopic Data

Convergent plate boundaries are key sites for continental crustal formation and recycling. Quantifying the evolution of crustal thickness and paleoelevation along ancient convergent margins represents a major goal in orogenic system analyses. Chemical and in some cases isotopic compositions of igneous rocks formed in modern supra-subduction arcs and collisional belts are sensitive to Moho depths at the location of magmatism, implying that igneous suites from fossil orogens carry information about crustal thickness from the time they formed. Several whole-rock chemical parameters correlate with crustal thickness, some of which were calibrated to serve as "mohometers," that is, quantitative proxies of paleo-Moho depths. Based on mineral-melt partition coefficients, this concept has been extended to detrital zircons, such that combined chemical and geochronological information extracted from these minerals allows us to reconstruct the crustal thickness evolution using the detrital archive. We discuss here the mohometric potential of a variety of chemical and isotopic parameters and show that their combined usage improves paleocrustal thickness estimates. Using a MATLAB® app developed for the underlying computations, we present examples from the modern and the deeper time geologic record to illustrate the promises and pitfalls of the technique. Since arcs are in isostatic equilibrium, mohometers are useful in reconstructing orogenic paleoelevation as well. Our analysis suggests that many global-scale correlations between magma composition and crustal thickness used in mohometry originate in the sub-arc mantle; additional effects resulting from intracrustal igneous differentiation depend on the compatible or incompatible behavior of the involved parameters.

Zinc systematics quantify crustal thickness control on fractionating assemblages of arc magmas

Understanding the processes leading to the broad chemical variability of arc magmas is an essential, yet not fully elucidated, issue in Earth Sciences. Here, I show that Zn–MgO–SiO₂ systematics of magmatic arc rocks correlate significantly with arc thickness. Because Zn–MgO–SiO₂ systematics are mostly controlled by fractionation of different mineral phases, this suggests a systematic change in the proportions of fractionating mineral assemblages depending on arc thickness. Using a mass balance model with a Monte Carlo approach, I show that Zn–MgO–SiO₂ systematics can be quantitatively explained by a continuous transition from plagioclase-dominated fractionating assemblages in thin arcs to amphibole-garnet-magnetite-dominated assemblages in increasingly thicker arcs. Most likely, such a systematic change results from the increase of average depth of magma differentiation that is ultimately controlled by arc thickness. Results presented have implications on the causes of different geochemical trends in arcs, the role of arcs as H₂O filters, and their association with porphyry deposits.

Data on the arc magmatism developed in the Antarctic Peninsula and Patagonia during the Late Triassic-Jurassic: A compilation of new and previous geochronology, geochemistry and isotopic tracing results

We present the results of U-Pb zircon dating conducted using laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS), isotopic tracing analyses of Hf in zircon and Sr-Nd in whole-rock and whole-rock major oxides, and trace element abundances of 12 plutonic and volcanic rocks present on the Antarctic Peninsula. The dataset is presented in combination with the results of previous studies conducted in both Patagonia and the Antarctic Peninsula. These results were filtered for concordant 206Pb-238U zircon ages and topology of the 40Ar/39Ar age spectra. These results may be useful for researchers studying the geological evolution of southern Gondwana, West Antarctica or Patagonia. The interpretation of this dataset is found in the co-submitted paper by Bastias, et al. (2021a) titled 'A revised interpretation of the Chon Aike magmatic province: active margin origin and implications for the opening of the Weddell Sea'.

Record of Crustal Thickening and Synconvergent Extension from the Dajiamang Tso Rift, Southern Tibet

North-trending rifts throughout south-central Tibet provide an opportunity to study the dynamics of synconvergent extension in contractional orogenic belts. In this study, we present new data from the Dajiamang Tso rift, including quantitative crustal thickness estimates calculated from trace/rare earth element zircon data, U-Pb geochronology, and zircon-He thermochronology. These data constrain the timing and rates of exhumation in the Dajiamang Tso rift and provide a basis for evaluating dynamic models of synconvergent extension. Our results also provide a semi-continuous record of Mid-Cretaceous to Miocene evolution of the Himalayan-Tibetan orogenic belt along the India-Asia suture zone. We report igneous zircon U-Pb ages of ~103 Ma and 70–42 Ma for samples collected from the Xigaze forearc basin and Gangdese Batholith/Linzizong Formation, respectively. Zircon-He cooling ages of forearc rocks in the hanging wall of the Great Counter thrust are ~28 Ma, while Gangdese arc samples in the footwalls of the Dajiamang Tso rift are 16–8 Ma. These data reveal the approximate timing of the switch from contraction to extension along the India-Asia suture zone (minimum 16 Ma). Crustal-thickness trends from zircon geochemistry reveal possible crustal thinning (to ~40 km) immediately prior to India-Eurasia collision onset (58 Ma). Following initial collision, crustal thickness increases to 50 km by 40 Ma with continued thickening until the early Miocene supported by regional data from the Tibetan Magmatism Database. Current crustal thickness estimates based on geophysical observations show no evidence for crustal thinning following the onset of E–W extension (~16 Ma), suggesting that modern crustal thickness is likely facilitated by an underthrusting Indian lithosphere balanced by upper plate extension.

Crustal thickening and endogenic oxidation of magmatic sulfur

Porphyry ore deposits, Earth's most important resources of copper, molybdenum, and rhenium, are strongly associated with felsic magmas showing signs of high-pressure differentiation and are usually found in places with thickened crust (>45 kilometers). This pattern is well-known, but unexplained, and remains an outstanding problem in our understanding of porphyry ore deposit formation. We approach this problem by investigating the oxidation state of magmatic sulfur, which controls the behavior of ore-forming metals during magma differentiation and magmatic-hydrothermal transition. We use sulfur in apatite to reconstruct the sulfur oxidation state in the Gangdese batholith, southern Tibet. We find that magma sulfate content increased abruptly after India-Eurasia collision. Apatite sulfur content and the calculated magma S6+/ΣS ratio correlate with whole-rock dysprosium/ytterbium ratio, suggesting that residual garnet, favored in thickened crust, exerts a first-order control on sulfur oxidation in magmatic orogens. Our findings link sulfur oxidation to internal petrogenic processes and imply an intrinsic relationship of magma oxidation with synmagmatic crustal thickening.

Gold endowments of porphyry deposits controlled by precipitation efficiency

Porphyry deposits are natural suppliers of most copper and significant gold to our society. Whereas the Cu-richest (Au-poor) porphyries are related to Andean-type subduction and typical calc-alkaline magmatism, the Au-richest porphyries are associated with high-K calc-alkaline to alkaline magmatism in late to post-subduction or post-collision and extensional settings, and subordinately with calc-alkaline magmatism. The reasons behind these associations and the large variations in metal endowments of porphyry Cu–Au deposits remain obscure. Here, I show that porphyry Cu–Au deposits define two distinct trends in Au vs. Cu tonnage plots (Cu-rich and Au-rich). Metal endowments for both trends grow larger the longer the mineralization process. However, Au is precipitated at much higher rates in Au-rich than in Cu-rich porphyry deposits. Using Monte Carlo simulations of petrologic processes, I show that whereas Cu-rich porphyries require large amounts of magma and water to be formed, Au-rich porphyries are the result of a better efficiency of Au precipitation.

Porphyry copper and gold deposits are the dominant natural suppliers of these metals to our society, yet the large variations in metal endowments of porphyry Cu–Au deposits remain obscure. Here, the author shows that Cu-rich porphyries require large amounts of magma and water to be formed, while Au-rich porphyries are the result of a better efficiency of Au precipitation.

Nb/Ta systematics in arc magma differentiation and the role of arclogites in continent formation

The surfaces of rocky planets are mostly covered by basaltic crust, but Earth is unique in that it also has extensive regions of felsic crust, manifested in the form of continents. Exactly how felsic crust forms when basaltic magmas are the dominant products of melting the mantles of rocky planets is unclear. A fundamental part of the debate is centered on the low Nb/Ta of Earth’s continental crust (11–13) compared to basalts (15–16). Here, we show that during arc magma differentiation, the extent of Nb/Ta fractionation varies with crustal thickness with the lowest Nb/Ta seen in continental arc magmas. Deep arc cumulates (arclogites) are found to have high Nb/Ta (average ~19) due to the presence of high Nb/Ta magmatic rutiles. We show that the crustal thickness control of Nb/Ta can be explained by rutile saturation being favored at higher pressures. Deep-seated magmatic differentiation, such as in continental arcs and other magmatic orogens, is thus necessary for making continents.

It is unclear why the crust among many arcs is so silica rich when there is a missing mafic component. Here the authors compile a global evaluation of Nb/Ta and examine xenoliths from Arizona to derive a model of why Nb/Ta fractionation is favored in the crust.

The redox "filter" beneath magmatic orogens and the formation of continental crust

The two most important magmatic differentiation series on Earth are the Fe-enriching tholeiitic series, which dominates the oceanic crust and island arcs, and the Fe-depleting calc-alkaline series, which dominates the continental crust and continental arcs. It is well known that calc-alkaline magmas are more oxidized when they erupt and are preferentially found in regions of thick crust, but why these quantities should be related remains unexplained. We use the redox-sensitive behavior of europium (Eu) in deep-seated, plagioclase-free arc cumulates to directly constrain the redox evolution of arc magmas at depth. Primitive arc cumulates have negative Eu anomalies, which, in the absence of plagioclase, can only be explained by Eu being partly reduced. We show that primitive arc magmas begin with low oxygen fugacities, similar to that of mid-ocean ridge basalts, but increase in oxygen fugacity by over two orders of magnitude during magmatic differentiation. This intracrustal oxidation is attended by Fe depletion coupled with fractionation of Fe-rich garnet. We conclude that garnet fractionation, owing to its preference for ferrous over ferric iron, results in simultaneous oxidation and Fe depletion of the magma. Favored at high pressure and water content, garnet fractionation explains the correlation between crustal thickness, oxygen fugacity, and the calc-alkaline character of arc magmas.

Andean surface uplift constrained by radiogenic isotopes of arc lavas

Climate and tectonics have complex feedback systems which are difficult to resolve and remain controversial. Here we propose a new climate-independent approach to constrain regional Andean surface uplift. ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd ratios of Quaternary frontal-arc lavas from the Andean Plateau are distinctly crustal (>0.705 and <0.5125, respectively) compared to non-plateau arc lavas, which we identify as a plateau discriminant. Strong linear correlations exist between smoothed elevation and ⁸⁷Sr/⁸⁶Sr (R² = 0.858, n = 17) and ¹⁴³Nd/¹⁴⁴Nd (R² = 0.919, n = 16) ratios of non-plateau arc lavas. These relationships are used to constrain 200 Myr of surface uplift history for the Western Cordillera (present elevation 4200 ± 516 m). Between 16 and 26°S, Miocene to recent arc lavas have comparable isotopic signatures, which we infer indicates that current elevations were attained in the Western Cordillera from 23 Ma. From 23–10 Ma, surface uplift gradually propagated southwards by ~400 km.

Multiple complex tectonic and climatic processes have formed the Andes, which today provides a unique ecological niche. Here, Scott et al. investigate how the chemical composition of lavas from stratovolcanoes can be used to give insight on the uplift of the Andes over the last 200 million years.

Evaluating the importance of metamorphism in the foundering of continental crust

The metamorphic conditions and mechanisms required to induce foundering in deep arc crust are assessed using an example of representative lower crust in SW New Zealand. Composite plutons of Cretaceous monzodiorite and gabbro were emplaced at ~1.2 and 1.8 GPa are parts of the Western Fiordland Orthogneiss (WFO); examples of the plutons are tectonically juxtaposed along a structure that excised ~25 km of crust. The 1.8 GPa Breaksea Orthogneiss includes suitably dense minor components (e.g. eclogite) capable of foundering at peak conditions. As the eclogite facies boundary has a positive dP/dT, cooling from supra-solidus conditions (T > 950 ºC) at high-P should be accompanied by omphacite and garnet growth. However, a high monzodioritic proportion and inefficient metamorphism in the Breaksea Orthogneiss resulted in its positive buoyancy and preservation. Metamorphic inefficiency and compositional relationships in the 1.2 GPa Malaspina Pluton meant it was never likely to have developed densities sufficiently high to founder. These relationships suggest that the deep arc crust must have primarily involved significant igneous accumulation of garnet–clinopyroxene (in proportions >75%). Crustal dismemberment with or without the development of extensional shear zones is proposed to have induced foundering of excised cumulate material at P > 1.2 GPa.

Quantifying Crustal Thickness in Continental Collisional Belts: Global Perspective and a Geologic Application

We present compiled geochemical data of young (mostly Pliocene-present) intermediate magmatic rocks from continental collisional belts and correlations between their whole-rock Sr/Y and La/Yb ratios and modern crustal thickness. These correlations, which are similar to those obtained from subduction-related magmatic arcs, confirm that geochemistry can be used to track changes of crustal thickness changes in ancient collisional belts. Using these results, we investigate temporal variations of crustal thickness in the Qinling Orogenic Belt in mainland China. Our results suggest that crustal thickness remained constant in the North Qinling Belt (~45–55 km) during the Triassic to Jurassic but fluctuates in the South Qinling Belt, corresponding to independently determined tectonic changes. In the South Qinling Belt, crustal thickening began at ~240 Ma and culminated with 60–70-km-thick crust at ~215 Ma. Then crustal thickness decreased to ~45 km at ~200 Ma and remained the same to the present. We propose that coupled use of Sr/Y and La/Yb is a feasible method for reconstructing crustal thickness through time in continental collisional belts. The combination of the empirical relationship in this study with that from subduction-related arcs can provide the crustal thickness evolution of an orogen from oceanic subduction to continental collision.

Cordillera Zealandia: A Mesozoic arc flare-up on the palaeo-Pacific Gondwana Margin

Two geochemically and temporally distinct components of the Mesozoic Zealandia Cordilleran arc indicate a shift from low to high Sr/Y whole rock ratios at c. 130 Ma. Recent mapping and a reappraisal of published Sr-Nd data combined with new in-situ zircon Hf isotope analyses supports a genetic relationship between the two arc components. A reappraisal of geophysical, geochemical and P-T estimates demonstrates a doubling in thickness of the arc to at least 80 km at c. 130 Ma. Contemporaneously, magmatic addition rates shifted from ~14 km³/my per km of arc to a flare-up involving ~100 km³/my per km of arc. Excursions in Sr-Nd-Hf isotopic ratios of flare-up rocks highlight the importance of crust-dominated sources. This pattern mimics Cordilleran arcs of the Americas and highlights the importance of processes occurring in the upper continental plates of subduction systems that are incompletely reconciled with secular models for continental crustal growth.

Quantifying crustal thickness over time in magmatic arcs

We present global and regional correlations between whole-rock values of Sr/Y and La/Yb and crustal thickness for intermediate rocks from modern subduction-related magmatic arcs formed around the Pacific. These correlations bolster earlier ideas that various geochemical parameters can be used to track changes of crustal thickness through time in ancient subduction systems. Inferred crustal thicknesses using our proposed empirical fits are consistent with independent geologic constraints for the Cenozoic evolution of the central Andes, as well as various Mesozoic magmatic arc segments currently exposed in the Coast Mountains, British Columbia, and the Sierra Nevada and Mojave-Transverse Range regions of California. We propose that these geochemical parameters can be used, when averaged over the typical lifetimes and spatial footprints of composite volcanoes and their intrusive equivalents to infer crustal thickness changes over time in ancient orogens.