The largest volcanic eruptions on Earth

Large igneous provinces (LIPs) are sites of the most frequently recurring, largest volume basaltic and silicic eruptions in Earth history. These large-volume (N1000 km3 dense rock equivalent) and large-magnitude (NM8) eruptions produce areally extensive (104–105 km2) basaltic lava flow fields and silicic ignimbrites that are the main building blocks of LIPs. Available information on the largest eruptive units are primarily from the Columbia River and Deccan provinces for the dimensions of flood basalt eruptions, and the Paraná–Etendeka and Afro-Arabian provinces for the silicic ignimbrite eruptions. In addition, three large-volume (675– 2000 km3) silicic lava flows have also been mapped out in the Proterozoic Gawler Range province (Australia), an interpreted LIP remnant. Magma volumes of N1000 km3 have also been emplaced as high-level basaltic and rhyolitic sills in LIPs. The data sets indicate comparable eruption magnitudes between the basaltic and silicic eruptions, but due to considerable volumes residing as co-ignimbrite ash deposits, the current volume constraints for the silicic ignimbrite eruptions may be considerably underestimated. Magma composition thus appears to be no barrier to the volume of magma emitted during an individual eruption. Despite this general similarity in magnitude, flood basaltic and silicic eruptions are very different in terms of eruption style, duration, intensity, vent configuration, and emplacement style. Flood basaltic eruptions are dominantly effusive and Hawaiian–Strombolian in style, with magma discharge rates of ~106–108 kg s−1 and eruption durations estimated at years to tens of years that emplace dominantly compound pahoehoe lava flow fields. Effusive and fissural eruptions have also emplaced some large-volume silicic lavas, but discharge rates are unknown, and may be up to an order of magnitude greater than those of flood basalt lava eruptions for emplacement to be on realistic time scales (b10 years). Most silicic eruptions, however, are moderately to highly explosive, producing co-current pyroclastic fountains (rarely Plinian) with discharge rates of 109– 1011 kg s−1 that emplace welded to rheomorphic ignimbrites. At present, durations for the large-magnitude silicic eruptions are unconstrained; at discharge rates of 109 kg s−1, equivalent to the peak of the 1991 Mt Pinatubo eruption, the largest silicic eruptions would take many months to evacuate N5000 km3 of magma. The generally simple deposit structure is more suggestive of short-duration (hours to days) and high intensity (~1011 kg s−1) eruptions, perhaps with hiatuses in some cases. These extreme discharge rates would be facilitated by multiple point, fissure and/or ring fracture venting of magma. Eruption frequencies are much elevated for large-magnitude eruptions of both magma types during LIP-forming episodes. However, in basaltdominated provinces (continental and ocean basin flood basalt provinces, oceanic plateaus, volcanic rifted margins), large magnitude (NM8) basaltic eruptions have much shorter recurrence intervals of 103–104 years, whereas similar magnitude silicic eruptions may have recurrence intervals of up to 105 years. The Paraná– Etendeka province was the site of at least nine NM8 silicic eruptions over an ~1 Myr period at ~132 Ma; a similar eruption frequency, although with a fewer number of silicic eruptions is also observed for the Afro- Arabian Province. The huge volumes of basaltic and silicic magma erupted in quick succession during LIP events raises several unresolved issues in terms of locus of magma generation and storage (if any) in the crust prior to eruption, and paths and rates of ascent from magma reservoirs to the surface.

New insights into Crustal Contributions to Large-volume Rhyolite Generation in the Mid-Tertiary Sierra Madre Occidental Province, Mexico, revealed by U Pb Geochronology

Voluminous (≥3·9 × 105 km3), prolonged (∼18 Myr) explosive silicic volcanism makes the mid-Tertiary Sierra Madre Occidental province of Mexico one of the largest intact silicic volcanic provinces known. Previous models have proposed an assimilation–fractional crystallization origin for the rhyolites involving closed-system fractional crystallization from crustally contaminated andesitic parental magmas, with <20% crustal contributions. The lack of isotopic variation among the lower crustal xenoliths inferred to represent the crustal contaminants and coeval Sierra Madre Occidental rhyolite and basaltic andesite to andesite volcanic rocks has constrained interpretations for larger crustal contributions. Here, we use zircon age populations as probes to assess crustal involvement in Sierra Madre Occidental silicic magmatism. Laser ablation-inductively coupled plasma-mass spectrometry analyses of zircons from rhyolitic ignimbrites from the northeastern and southwestern sectors of the province yield U–Pb ages that show significant age discrepancies of 1–4 Myr compared with previously determined K/Ar and 40Ar/39Ar ages from the same ignimbrites; the age differences are greater than the errors attributable to analytical uncertainty. Zircon xenocrysts with new overgrowths in the Late Eocene to earliest Oligocene rhyolite ignimbrites from the northeastern sector provide direct evidence for some involvement of Proterozoic crustal materials, and, potentially more importantly, the derivation of zircon from Mesozoic and Eocene age, isotopically primitive, subduction-related igneous basement. The youngest rhyolitic ignimbrites from the southwestern sector show even stronger evidence for inheritance in the age spectra, but lack old inherited zircon (i.e. Eocene or older). Instead, these Early Miocene ignimbrites are dominated by antecrystic zircons, representing >33 to ∼100% of the dated population; most antecrysts range in age between ∼20 and 32 Ma. A sub-population of the antecrystic zircons is chemically distinct in terms of their high U (>1000 ppm to 1·3 wt %) and heavy REE contents; these are not present in the Oligocene ignimbrites in the northeastern sector of the Sierra Madre Occidental. The combination of antecryst zircon U–Pb ages and chemistry suggests that much of the zircon in the youngest rhyolites was derived by remelting of partially molten to solidified igneous rocks formed during preceding phases of Sierra Madre Occidental volcanism. Strong Zr undersaturation, and estimations for very rapid dissolution rates of entrained zircons, preclude coeval mafic magmas being parental to the rhyolite magmas by a process of lower crustal assimilation followed by closed-system crystal fractionation as interpreted in previous studies of the Sierra Madre Occidental rhyolites. Mafic magmas were more probably important in providing a long-lived heat and material flux into the crust, resulting in the remelting and recycling of older crust and newly formed igneous materials related to Sierra Madre Occidental magmatism.

Petrology and geochemistry of the Quaternary Caldera-forming, Phonolitic Granadilla Eruption, Tenerife (Canary Islands)

The Granadilla eruption at 600 ka was one of the largest phonolitic explosive eruptions from the Las Cañadas volcano on Tenerife, producing a classical plinian eruptive sequence of a widespread pumice fall deposit overlain by an ignimbrite. The eruption resulted in a major phase of caldera collapse that probably destroyed the shallow-level magma chamber system. Granadilla pumices contain a diverse phenocryst assemblage of alkali feldspar + biotite + sodian diopside to aegirine–augite + titanomagnetite + ilmenite + nosean/haüyne + titanite + apatite; alkali feldspar is the dominant phenocryst and biotite is the main ferromagnesian phase. Kaersutite and partially resorbed plagioclase (oligoclase to sodic andesine) are present in some eruptive units, particularly in pumice erupted during the early plinian phase, and in the Granadilla ignimbrite at the top of the sequence. Associated with the kaersutite and plagioclase are small clots of microlitic plagioclase and kaersutite interpreted as quenched blebs of tephriphonolitic magma within the phonolite pumice. The Granadilla Member has previously been recognized as an example of reverse-then-normal compositional zonation, where the zonation is primarily expressed in terms of substantial variations in trace element abundances with limited major element variation (cryptic zonation). Evidence for cryptic zonation is also provided by the chemistry of the phenocryst phases, and corresponding changes in intensive parameters (e.g. T, f O2, f  H2O). Geothermometry estimates indicate that the main body of phonolite magma had a temperature gradient from 860 °C to ∼790 °C, with hotter magma (≥900 °C) tapped at the onset and terminal phases of the eruption. The reverse-then-normal chemical and thermal zonation reflects the initial tapping of a partially hybridized magma (mixing of phonolite and tephriphonolite), followed by the more sequential tapping of a zoned and relatively large body of highly evolved phonolite at a new vent and during the main plinian phase. This suggests that the different magma types within the main holding chamber could have been laterally juxtaposed, as well as in a density-stratified arrangement. Correlations between the presence of mixed phenocryst populations (i.e. presence of plagioclase and kaersutite) and coarser pumice fall layers suggest that increased eruption vigour led to the tapping of hybridized and/or less evolved magma probably from greater depths in the chamber. New oxygen isotope data for glass and mineral separates preclude syn-eruptive interaction between the vesiculating magma and hydrothermal fluids as the cause of the Sr isotope disequilibrium identified previously for the deposit. Enrichment in radiogenic Sr in the pumice glass has more likely been due to low-temperature exchange with meteoric water that was enriched in 87Sr by sea spray, which may be a common process affecting porous and glassy pyroclastic deposits on oceanic islands.

Pulling apart the mid to late Cenozoic magmatic record of the Gulf of California : no room for the Comondú Arc

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. Significantly, incorrect tectonic interpretations can be made when the spatial-temporal-compositional trends of, and source contributions to magmatism are not properly considered. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

25 million years to break a continent : early to middle Miocene rifting and syn-extensional magmatism in the southern Gulf of California

Cenozoic extension in western Mexico has been divided into two episodes separated by the change from convergence to oblique divergence at the plate boundary. The Gulf Extensional Province is thought to have started once subduction ended at ~12.5 Ma whereas early extension is classified as Basin and Range. Mid-Miocene volcanism of the Comondú group has been considered as a subduction-related arc, whereas post ~12.5 Ma volcanism would be extension-related. Our new integration of the continental onshore and offshore geology of the south-east Gulf region, backed by tens of Ar-Ar and U-Pb ages and geochemical studies, document an early-mid Miocene rifting and extension-related bimodal to andesitic magmatism prior to subduction termination. Between ~21 and 11 Ma a system of NNW-SSE high-angle extensional faults rifted the western side of the Sierra Madre Occidental (SMO) ignimbrite plateau. In Nayarit, rhyolitic domes and some basalts were emplaced along this extensional belt at 18-17 Ma. These rocks show strong antecrystic inheritance but an absence of Mesozoic and older xenocrysts, suggesting a genesis in the mid-upper crust triggered by extension-induced basaltic influx. In Sinaloa, large grabens were floored by huge dome complexes at ~21-17 Ma and filled by continental sediments with interlayered basalts dated at 15 Ma. Mid-Miocene volcanism, including the largely volcaniclastic Comondú strata in Baja California, was thus emplaced in rift basins and appears associated to decompression melting rather than subduction. Along the coast, flat-lying basaltic lava flows dated at 11-10 Ma are exposed just above the present sea level. Here crustal thickness is 25-20 Km, almost half that in the core of the SMO, implying significant lithosphere stretching before ~11 Ma. This mafic pulse, with relatively high Ti but still clear Nb-Ta negative spikes, may be related to the detachment of the lower part of the subducted slab, allowing asthenosphere to flow into parts of the mantle previously fluxed by subduction fluids. Very uniform OIB-like lavas appear in late Pliocene and Pleistocene, only 18 m.y. after the onset of rifting and ~9 m.y. after the end of subduction. Our study shows that rifting began much earlier than Late Miocene and progressively overwhelmed subduction in generating magmatism.

Pulling apart the mid to late Cenozoic magmatic record of the Gulf of California : is there a Comondú arc?

The composition of the lithosphere can be fundamentally altered by long-lived subduction processes such that subduction-modified lithosphere can survive for 100's Myrs. Incorrect petrotectonic interpretations result when spatial-temporal-compositional trends of, and source contributions to, magmatism are not properly considered. Western Mexico has had protracted Cenozoic magmatism developed mostly in-board of active oceanic plate subduction beneath western North America. A broad range of igneous compositions from basalt to high-silica rhyolite were erupted with intermediate to silicic compositions in particular, showing calc-alkaline and other typical subduction-related geochemical signatures. A major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a bimodal volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Extension became more focussed ~18 Ma resulting in localised volcanic activity along the future site of the Gulf of California. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. Past tectonic interpretations of Comondú-age volcanism may have been incorrect as these regional temporal-compositional changes are alternatively interpreted as a result of increased mixing of mantle-derived basaltic and crust-derived rhyolitic magmas in an active rift environment rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate.

Consequences of extensional tectonics on volcanic eruption style, compositions and source regions : new insights from the southern Sierra Madre Occidental and Gulf of California regions, western Mexico

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

The largest eruptions in Earth's history

Large igneous provinces (LIPs) host the most frequently recurring, largest volume basaltic & silicic eruptions on Earth. The largest volume (>1000 km^3 DRE) and magnitude (>M8) eruptions produce areally extensive (10^4-10^5 km^2) basaltic flow fields and sills, and silicic ignimbrites that are the main LIP building blocks. Basaltic and silicic eruptions have comparable magnitudes, but silicic ignimbrite volumes may be significantly underestimated due to unrecognized and correlated, but voluminous co-ignimbrite ash deposits. Magma composition is no barrier to individual eruption volume. Despite similar magnitudes, flood basaltic and silicic eruptions are very different in eruption mechanism, duration, intensity, vent configuration, and emplacement style. Flood basalts are dominantly effusive Hawaiian-Strombolian, with magma discharge rates of ~10^7-10^8 kg s^-1, and produce dominantly compound pahoehoe flow fields over eruption durations most likely >10 yrs. Most silicic eruptions are moderately to highly explosive, producing cocurrent pyroclastic fountains (rarely Plinian) and suggested to be of short-duration (hours to days) and high intensity (~10^11 kg s^-1). Eruption frequencies are elevated for largemagnitude eruptions of both magma types during LIP formation. In basalt-dominated provinces, large magnitude (>M8) eruptions have much shorter recurrence intervals (10^3-10^4 years) than similar magnitude silicic eruptions (~10^5 years). The huge volumes of magma erupted rapidly in LIPs raises several unresolved issues in terms of locus of magma generation and storage (if any) in the crust prior to eruption, the paths and rates of ascent from magma reservoirs to the surface, and relative aerosol contributions to the stratosphere from the flood basaltic and rhyolitic eruptions.