Volcanic–plutonic parity and the differentiation of the continental crust

A zircon case for super-wet arc magmas

Arc magmas have higher water contents (2-6 wt.% H2O) than magmas generated in other tectonic environments, with a growing body of evidence suggesting that some deep arc magmas may be 'super-wet' (>6 wt.% H2O). Here, we use thermodynamic modelling to show that the behaviour of zirconium during magmatic differentiation is strongly sensitive to melt water contents. We demonstrate that super-wet magmas crystallise zircon with low, homogeneous titanium concentrations (75th percentile <10 ppm) due to a decrease in zircon saturation temperatures with increasing melt H2O. We find that zircon titanium concentrations record a transition to super-wet magmatism in Central Chile immediately before the formation of the world's largest porphyry copper deposit cluster at Río Blanco-Los Bronces. Broader analysis shows that low, homogeneous zircon titanium concentrations are present in many magmatic systems. Our study suggests that super-wet magmas are more common than previously envisaged and are fundamental to porphyry copper deposit mineralisation.

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.

Extensive crystal fractionation of high-silica magmas revealed by K isotopes

Fractional crystallization plays a critical role in generating the differentiated continental crust on Earth. However, whether efficient crystal-melt separation can occur in viscous felsic magmas remains a long-standing debate because of the difficulty in discriminating between differentiated melts and complementary cumulates. Here, we found large (~1 per mil) potassium isotopic variation in 54 strongly peraluminous high-silica (silicon dioxide >70 weight %) leucogranites from the Himalayan orogen, with potassium isotopes correlated with trace elemental proxies (e.g., strontium, rubidium/strontium, and europium anomaly) for plagioclase crystallization. Quantitative modeling requires up to ~60 to 90% fractional crystallization to account for the progressively light potassium isotopic composition of the fractionated leucogranites, while plagioclase accumulation results in enrichment of heavy potassium isotopes in cumulate leucogranites. Our findings strongly support fractional crystallization of high-silica magmas and highlight the great potential of potassium isotopes in studying felsic magma differentiation.

Interplay between oceanic subduction and continental collision in building continental crust

Generation of continental crust in collision zones reflect the interplay between oceanic subduction and continental collision. The Gangdese continental crust in southern Tibet developed during subduction of the Neo-Tethyan oceanic slab in the Mesozoic prior to reworking during the India-Asia collision in the Cenozoic. Here we show that continental arc magmatism started with fractional crystallization to form cumulates and associated medium-K calc-alkaline suites. This was followed by a period commencing at ~70 Ma dominated by remelting of pre-existing lower crust, producing more potassic compositions. The increased importance of remelting coincides with an acceleration in the convergence rate between India and Asia leading to higher basaltic flow into the Asian lithosphere, followed by convergence deceleration due to slab breakoff, enabling high heat flow and melting of the base of the arc. This two-stage process of accumulation and remelting leads to the chemical maturation of juvenile continental crust in collision zones, strengthening crustal stratification.

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.

Episodic growth of felsic continents in the past 3.7 Ga

Continents form the most accessible parts of Earth, but their complex compositions make their origin difficult to investigate. A novel approach based on a comprehensive compilation of samarium-neodymium isotopic compositions of detrital sedimentary rocks is here used to unravel continental growth through time. This record reveals that continents were as felsic as today in the past 3.7 Ga (billion years) and that their growth was not continuous but episodic. Reworking of preexisting crust was a ubiquitous process during most of Earth history, but at least six periods of continental growth can be identified every 500 to 700 Ma (million years) in the past 3.7 Ga. This recurrence could be accounted for by changes in tectonic plate velocities favoring periods of rapid subduction and enhanced production of juvenile felsic crust.

Lithium systematics in global arc magmas and the importance of crustal thickening for lithium enrichment

Much of the world's Li deposits occurs as basinal brines in magmatic orogens, particularly in continental volcanic arcs. However, the exact origin of Li enrichment in arc magmatic systems is not clear. Here, we show that, globally, primitive arc magmas have Li contents and Li/Y ratios similar to mid-ocean ridge basalts, indicating that the subducting slab has limited contribution to Li enrichment in arc magmas. Instead, we find that Li enrichment is enhanced by lower degrees of sub-arc mantle melting and higher extents of intracrustal differentiation. These enrichment effects are favored in arcs with thick crust, which explains why magmatism and differentiation in continental arcs, like the Andes, reach greater Li contents than their island arc counterparts. Weathering of these enriched source rocks mobilizes and transports such Li into the hydrologic system, ultimately developing Li brines with the combination of arid climate and the presence of landlocked extensional basins in thickened orogenic settings.

Constraining crustal silica on ancient Earth

Accurately quantifying the composition of continental crust on Hadean and Archean Earth is critical to our understanding of the physiography, tectonics, and climate of our planet at the dawn of life. One longstanding paradigm involves the growth of a relatively mafic planetary crust over the first 1 to 2 billion years of Earth history, implying a lack of modern plate tectonics and a paucity of subaerial crust, and consequently lacking an efficient mechanism to regulate climate. Others have proposed a more uniformitarian view in which Archean and Hadean continents were only slightly more mafic than at present. Apart from complications in assessing early crustal composition introduced by crustal preservation and sampling biases, effects such as the secular cooling of Earth's mantle and the biologically driven oxidation of Earth's atmosphere have not been fully investigated. We find that the former complicates efforts to infer crustal silica from compatible or incompatible element abundances, while the latter undermines estimates of crustal silica content inferred from terrigenous sediments. Accounting for these complications, we find that the data are most parsimoniously explained by a model with nearly constant crustal silica since at least the early Archean.

Geologic evidence for an icehouse Earth before the Sturtian global glaciation

Snowball Earth episodes, times when the planet was covered in ice, represent the most extreme climate events in Earth's history. Yet, the mechanisms that drive their initiation remain poorly constrained. Current climate models require a cool Earth to enter a Snowball state. However, existing geologic evidence suggests that Earth had a stable, warm, and ice-free climate before the Neoproterozoic Sturtian global glaciation [ca. 717 million years (Ma) ago]. Here, we present eruption ages for three felsic volcanic units interbedded with glaciolacustrine sedimentary rocks from southwest Virginia, USA, that demonstrate that glacially influenced sedimentation occurred at tropical latitudes ca. 751 Ma ago. Our findings are the first geologic evidence of a cool climate teetering on the edge of global glaciation several million years before the Sturtian Snowball Earth.

The causes of spatiotemporal variations in erupted fluxes and compositions along a volcanic arc

Decades of study on volcanic arcs have provided insight into the overarching processes that control magmatism, and how these processes manifest at individual volcanoes. However, the causes of ubiquitous and dramatic intra-arc variations in volcanic flux and composition remain largely unresolved. Investigating such arc-scale issues requires greater quantitative comparison of geophysical and geochemical data, linked through sets of common intensive variables. To work towards these goals, we use observed lava compositions to estimate the heat budget associated with Quaternary volcanism in the Cascades Arc and compare this to the heat required to produce the observed geophysical properties of the crust. Here we show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily related to variations in the flux of basalt into the crust, rather than variations in their crustal storage history. This approach shows promise for studying other large-scale frontier geologic problems in volcanic arcs.

The primary causes of dramatic variations in volcanic flux and composition along strike in subduction zones remain largely unknown. Here we use a promising new approach to show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily due to variations in the flux of basalt into the base of the crust, rather than crustal magma storage.

U-Pb constraints on pulsed eruption of the Deccan Traps across the end-Cretaceous mass extinction

Temporal correlation between some continental flood basalt eruptions and mass extinctions has been proposed to indicate causality, with eruptive volatile release driving environmental degradation and extinction. We tested this model for the Deccan Traps flood basalt province, which, along with the Chicxulub bolide impact, is implicated in the Cretaceous-Paleogene (K-Pg) extinction approximately 66 million years ago. We estimated Deccan eruption rates with uranium-lead (U-Pb) zircon geochronology and resolved four high-volume eruptive periods. According to this model, maximum eruption rates occurred before and after the K-Pg extinction, with one such pulse initiating tens of thousands of years prior to both the bolide impact and extinction. These findings support extinction models that incorporate both catastrophic events as drivers of environmental deterioration associated with the K-Pg extinction and its aftermath.

Neoproterozoic glacial origin of the Great Unconformity

The Great Unconformity, a profound gap in Earth's stratigraphic record often evident below the base of the Cambrian system, has remained among the most enigmatic field observations in Earth science for over a century. While long associated directly or indirectly with the occurrence of the earliest complex animal fossils, a conclusive explanation for the formation and global extent of the Great Unconformity has remained elusive. Here we show that the Great Unconformity is associated with a set of large global oxygen and hafnium isotope excursions in magmatic zircon that suggest a late Neoproterozoic crustal erosion and sediment subduction event of unprecedented scale. These excursions, the Great Unconformity, preservational irregularities in the terrestrial bolide impact record, and the first-order pattern of Phanerozoic sedimentation can together be explained by spatially heterogeneous Neoproterozoic glacial erosion totaling a global average of 3-5 vertical kilometers, along with the subsequent thermal and isostatic consequences of this erosion for global continental freeboard.

Chemical differentiation, cold storage and remobilization of magma in the Earth's crust

The formation, storage and chemical differentiation of magma in the Earth's crust is of fundamental importance in igneous geology and volcanology. Recent data are challenging the high-melt-fraction 'magma chamber' paradigm that has underpinned models of crustal magmatism for over a century, suggesting instead that magma is normally stored in low-melt-fraction 'mush reservoirs'^1-9. A mush reservoir comprises a porous and permeable framework of closely packed crystals with melt present in the pore space^1,10. However, many common features of crustal magmatism have not yet been explained by either the 'chamber' or 'mush reservoir' concepts^1,11. Here we show that reactive melt flow is a critical, but hitherto neglected, process in crustal mush reservoirs, caused by buoyant melt percolating upwards through, and reacting with, the crystals¹⁰. Reactive melt flow in mush reservoirs produces the low-crystallinity, chemically differentiated (silicic) magmas that ascend to form shallower intrusions or erupt to the surface^11-13. These magmas can host much older crystals, stored at low and even sub-solidus temperatures, consistent with crystal chemistry data^6-9. Changes in local bulk composition caused by reactive melt flow, rather than large increases in temperature, produce the rapid increase in melt fraction that remobilizes these cool- or cold-stored crystals. Reactive flow can also produce bimodality in magma compositions sourced from mid- to lower-crustal reservoirs^14,15. Trace-element profiles generated by reactive flow are similar to those observed in a well studied reservoir now exposed at the surface¹⁶. We propose that magma storage and differentiation primarily occurs by reactive melt flow in long-lived mush reservoirs, rather than by the commonly invoked process of fractional crystallization in magma chambers¹⁴.

How deceptive are microstructures in granitic rocks? Answers from integrated physical theory, phase equilibrium, and direct observations

In this contribution, we address the vexed question of the extent to which microstructures in granitic rocks reflect their igneous histories or have been masked by later events. The previous works have tended to address the problem either using theoretical or modelling considerations, or by interpretation of observed microstructures. Here, we use an approach that integrates the theory of microstructural development and the results of experimental phase-equilibrium studies with direct observation of natural examples on a variety of scales. We show that the predictions of the theoretical and experimental approaches agree perfectly with the mesoscopic and microscopic evidence from granitic rocks themselves. Our conclusion is that although, in many cases, granitic rock microstructures have been modified by near-solidus reactions and crystallisation, in the absence of tectonic deformation the fundamental elements of their igneous heritage remain intact. This means that it is perfectly in order to infer aspects of crystallisation sequences, magmatic reactions, and magma flow through careful microstructural observations. Thus, our answer to the question of how deceptive granitic textures are is, in most instances, 'not very'. However, some undeformed plutons have undergone fluid-driven alteration, and others have been affected by contact metamorphism. Thus, each case should be examined on its own merits.

Potassic, high-silica Hadean crust

Understanding Hadean (>4 Ga) Earth requires knowledge of its crust. The composition of the crust and volatiles migrating through it directly influence the makeup of the atmosphere, the composition of seawater, and nutrient availability. Despite its importance, there is little known and less agreed upon regarding the nature of the Hadean crust. By analyzing the ⁸⁷Sr/⁸⁶Sr ratio of apatite inclusions in Archean zircons from Nuvvuagittuq, Canada, we show that its protolith had formed a high (>1) Rb/Sr ratio reservoir by at least 4.2 Ga. This result implies that the early crust had a broad range of igneous rocks, extending from mafic to highly silicic compositions.

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.

A record of deep-ocean dissolved O2 from the oxidation state of iron in submarine basalts

The oxygenation of the deep ocean in the geological past has been associated with a rise in the partial pressure of atmospheric molecular oxygen (O₂) to near-present levels and the emergence of modern marine biogeochemical cycles. It has also been linked to the origination and diversification of early animals. It is generally thought that the deep ocean was largely anoxic from about 2,500 to 800 million years ago, with estimates of the occurrence of deep-ocean oxygenation and the linked increase in the partial pressure of atmospheric oxygen to levels sufficient for this oxygenation ranging from about 800 to 400 million years ago. Deep-ocean dissolved oxygen concentrations over this interval are typically estimated using geochemical signatures preserved in ancient continental shelf or slope sediments, which only indirectly reflect the geochemical state of the deep ocean. Here we present a record that more directly reflects deep-ocean oxygen concentrations, based on the ratio of Fe³⁺ to total Fe in hydrothermally altered basalts formed in ocean basins. Our data allow for quantitative estimates of deep-ocean dissolved oxygen concentrations from 3.5 billion years ago to 14 million years ago and suggest that deep-ocean oxygenation occurred in the Phanerozoic (541 million years ago to the present) and potentially not until the late Palaeozoic (less than 420 million years ago).

Seismic signature of active intrusions in mountain chains

Intrusions are a ubiquitous component of mountain chains and testify to the emplacement of magma at depth. Understanding the emplacement and growth mechanisms of intrusions, such as diapiric or dike-like ascent, is critical to constrain the evolution and structure of the crust. Petrological and geological data allow us to reconstruct magma pathways and long-term magma differentiation and assembly processes. However, our ability to detect and reconstruct the short-term dynamics related to active intrusive episodes in mountain chains is embryonic, lacking recognized geophysical signals. We analyze an anomalously deep seismic sequence (maximum magnitude 5) characterized by low-frequency bursts of earthquakes that occurred in 2013 in the Apennine chain in Italy. We provide seismic evidences of fluid involvement in the earthquake nucleation process and identify a thermal anomaly in aquifers where CO₂ of magmatic origin dissolves. We show that the intrusion of dike-like bodies in mountain chains may trigger earthquakes with magnitudes that may be relevant to seismic hazard assessment. These findings provide a new perspective on the emplacement mechanisms of intrusive bodies and the interpretation of the seismicity in mountain chains.

Vertically extensive and unstable magmatic systems: A unified view of igneous processes

Volcanoes are an expression of their underlying magmatic systems. Over the past three decades, the classical focus on upper crustal magma chambers has expanded to consider magmatic processes throughout the crust. A transcrustal perspective must balance slow (plate tectonic) rates of melt generation and segregation in the lower crust with new evidence for rapid melt accumulation in the upper crust before many volcanic eruptions. Reconciling these observations is engendering active debate about the physical state, spatial distribution, and longevity of melt in the crust. Here we review evidence for transcrustal magmatic systems and highlight physical processes that might affect the growth and stability of melt-rich layers, focusing particularly on conditions that cause them to destabilize, ascend, and accumulate in voluminous but ephemeral shallow magma chambers.

Long-term Geochemical Evolution of Lithogenic Versus Anthropogenic Distribution of Macro and Trace Elements in Household Attic Dust

Attic dusts were examined as historical archives of anthropogenic emissions, with the goal of elucidating the enrichment pathways associated with hydrothermal exploitation of Cu, Pb, and Zn minerals in the Bregalnica River basin in the eastern part of the Republic of Macedonia. Dust samples were collected from 84 settlements. Atomic emission spectrometry and mass spectrometry with inductively coupled plasma were applied as analytical techniques for the determination of 69 element contents. Multivariate analysis was applied for the extraction of dominant geochemical markers. The lithogenic distribution was simplified to six dominant geochemical markers: F1: Ga-Nb-Ta-Y-(La-Gd)-(Eu-Lu); F2: Be-Cr-Li-Mg-Ni; F3: Ag-Bi-Cd-Cu-In-Mn-Pb-Sb-Te-W-Zn; F4: Ba-Cs-Hf-Pd-Rb-Sr-Tl-Zr; F5: As-Co-Ge-V; and F6: К-Na-Sc-Ti. The anthropogenic effects on the air pollution were marked by a dominance of F3 and secondary dominance of F5. The fifth factor also was determined as a lithogenic marker for the occurrence of the very old Rifeous shales. The first factor also presents a very unique association that despite the heterogeneity relays on natural phenomena of tracking the deposition in areas of Proterosoic gneisses; related to the distribution of fine particles was associated with carbonate-silicate volcanic rocks. Intensive poly-metallic dust depositions were recorded only in the surroundings of localities where the hydrothermal extractions are implemented. Long-term deposition can be considered as pollution indexes for these hot spots. This mainly affects the Cd, Pb, and Zn deposition that is as high as 25, 3900, and 3200 mg/kg, respectively.

Multivariate extraction of dominant geochemical markers for deposition of 69 elements in the Bregalnica River basin, Republic of Macedonia (moss biomonitoring)

Atmospheric deposition was investigated using the terrestrial moss species Hypnum cupressiforme (Hedw.) and Homolothecium lutescens (Hedw.) in the Bregalnica River basin, Republic of Macedonia. Long-term emission occurs in this area due to the hydrothermal exploitation of Pb-Zn deposits (Sasa and Zletovo mines) and copper ore exploitation and floatation (Bučim mine). Determination of the chemical elements was conducted using atomic emission spectrometry with inductively coupled plasma (ICP-AES) and mass spectrometry with inductively coupled plasma (ICP-MS). A combination of multivariate techniques (PCA, FA and CA) was applied for data processing and identification of element association with lithogenic/anthropogenic origin. Seven dominant factors were extracted from the total of 69 analysed elements. Spatial distribution maps were constructed for the determination and localisation of smaller areas with higher contents of certain anthropogenic elements. In this way, the influences of selected human activities on local air pollution can be determined. The summarised data show quantification of the element distributions. This not only allows the determination of the distribution of hazardous elements but also presents complete characterisation of element deposition in the environs of mines.

Archean upper crust transition from mafic to felsic marks the onset of plate tectonics

The Archean Eon witnessed the production of early continental crust, the emergence of life, and fundamental changes to the atmosphere. The nature of the first continental crust, which was the interface between the surface and deep Earth, has been obscured by the weathering, erosion, and tectonism that followed its formation. We used Ni/Co and Cr/Zn ratios in Archean terrigenous sedimentary rocks and Archean igneous/metaigneous rocks to track the bulk MgO composition of the Archean upper continental crust. This crust evolved from a highly mafic bulk composition before 3.0 billion years ago to a felsic bulk composition by 2.5 billion years ago. This compositional change was attended by a fivefold increase in the mass of the upper continental crust due to addition of granitic rocks, suggesting the onset of global plate tectonics at ~3.0 billion years ago.