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.

Pre-eruptive uplift in the Emeishan? [Correspondence]

In their correspondence, He and colleagues question our conclusion of little or no uplift preceding Emeishan volcanism that we reported in our letter1. Debate concerns the nature of the contact between the Maokou limestone and Emeishan volcanics, the depositional environment and volumetric significance of mafic hydromagmatic deposits (MHDs), and evidence for symmetrical domal thinning. MHDs in the Daqiao section are separated from the Maokou limestone by 100 m of subaerial basaltic lavas, but elsewhere MHDs — previously interpreted as basal conglomerates2, 3 — directly overlie the Maokou2, 3. MHDs thus feature strongly in basal sections of the Emeishan lava succession, as also recently shown4 elsewhere in the Emeishan. An irregular surface at the top of the Maokou limestone has been interpreted as an erosional unconformity2, 3, but clastic deposits presented as evidence of this erosion2, 3 are MHDs produced by explosive magma–water interaction1. A clear demonstration that this irregular top surface is an erosional truncation of limestone reef facies (slope/rim, flat, lagoonal) is currently lacking, but is critical because reefs and carbonate platforms show considerable natural relief of tens of metres. The persistent hot, wet climate since the Oligocene has produced well-developed weathering profiles on exposed Palaeozoic marine sedimentary sequences5, but weathering and karst relief of the uppermost Maokou limestone underlying the flood basalts have not been properly documented, nor shown to be of middle Permian age and immediately preceding emplacement of the large igneous province.

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.

The El Chichon stratospheric cloud : solid particulates and settling rates

Sampling of the El Chichón stratospheric cloud in early May and in late July, 1982, showed that a significant proportion of the cloud consisted of solid particles between 2 μm and 40 μm size. In addition, many particles may have been part of larger aggregates or clusters that ranged in size from < 10 μm to > 50 μm. The majority of individual grains were angular aluminosilicate glass shards with various amounts of smaller, adhering particles. Surface features on individual grains include sulfuric acid droplets and larger (0.5 μm to 1 μm) sulfate gel droplets with various amounts of Na, Mg, Ca and Fe. The sulfate gels probably formed by the interaction of sulfur-rich gases and solid particles within the cloud soon after eruption. Ca-sulfate laths may have formed by condensation within the plume during eruption, or alternatively, at a later stage by the reaction of sulfuric acid aerosols with ash fragments within the stratospheric cloud. A Wilson-Huang formulation for the settling rate of individual particles qualitatively agrees with the observed particle-size distribution for a period at least four months after injection of material into the stratosphere. This result emphasizes the importance of particle shape in controlling the settling rate of volcanic ash from the stratosphere.

Were intercalated komatiites and dacites at the Black Swan nickel sulphide mine, Yilgarn Craton, Western Australia, emplaced as extrusive lavas or intrusive bodies? The significance of breccia textures and contact relationships

Intercalated Archean komatiites and dacites sit above a thick footwall dacite unit in the host rock succession at the Black Swan Nickel Mine, north of Kalgoorlie in the Yilgarn Craton, Western Australia. Both lithofacies occur in units that vary in scale from laterally extensive at the scale of the mine lease to localized, thin, irregular bodies, from > 100 m thick to only centimetres thick. Some dacites are only slightly altered and deformed, and are interpreted to post-date major deformation and alteration (late porphyries). However, the majority of the dacites display evidence of deformation, especially at contacts, and metamorphism, varying from silicification and chlorite alteration at contacts to pervasive low grade regional metamorphic alteration represented by common assemblages of chlorite, sericite and albite. Texturally, the dacites vary from entirely massive and coherent to partially brecciated to totally brecciated. Strangely, some dacites are coherent at the margins and brecciated internally. Breccia textures vary from cryptically defined, to blocky, closely packed, in situ jig-saw fit textures with secondary minerals in fractures between clasts, to more apparent matrix rich textures with round clast forms, giving apparent conglomerate textures. Some clast zones have multi-coloured clasts, giving the impression of varied provenance. Strangely however, all these textural variants have gradational relationships with each other, and no bedding or depositional structures are present. This indicates that all textures have an in situ origin. The komatiites are generally altered and pervasively carbonate veined. Preservation of original textures is patchy and local, but includes coarse adcumulate, mesocumulate, orthocumulate, crescumulate-harrisite and occasionally spinifex textures. Where original contacts between komatiites and dacites are preserved intact (i.e. not sheared or overprinted by alteration), the komatiites have chilled margins, whereas the dacites do not. The margins of the dacites are commonly silicified, and inclusions of dacite occur in komatiite, even at the top contacts of komatiite units, but komatiite clasts do not occur in the dacites. The komatiites therefore were emplaced as sills into the dacites, and the intercalated relationships are interpreted as intrusive. The brecciation and alteration in the dacites are interpreted as being largely due to hydraulic fracturing and alteration induced by contact metamorphic effects and hydrothermal alteration deriving from the intrusion of komatiites into the felsic pile. The absence of autobreccia and hyaloclastite textures in the dacites suggest that they were emplaced as an earlier intrusive (sill?) complex at a high level in the crust.

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.

Late Oligocene to middle Miocene rifting and syn-extensional magmatism in the southwestern Sierra Madre occidental, Mexico: the beginning of the Gulf of California rift

Although Basin and Range style extension affected several areas of western Mexico since the Late Eocene, extension in the Gulf of California region (the Gulf Extensional Province GEP) is thought to have started as subduction waned and ended at ~14 12.5 Ma. A general consensus also exists in considering the mid Miocene Comondú group as a suprasubduction volcanic arc. Our new integration of the geology of the south east Gulf region, backed by 43 new Ar Ar and U Pb mineral ages and geochemical studies, document a widespread phase of extension in the southern GEP between latest Oligocene and Early Miocene that subsequently focused in the region of the future Gulf in the Middle Miocene. Upper Oligocene to Lower Miocene rocks across the southern Sierra Madre Occidental (SMO)(northern Nayarit and southern Sinaloa) were affected by major ~N S to NNW striking normal faults prior to ~21 Ma. Then, between ~21 and 11 Ma, a system of NNW-SSE high angle extensional faults continued extending the southwestern side of the SMO. Rhyolitic domes, shallow intrusive bodies, and lesser basalts were emplaced along this extensional belt at 20-17 Ma. In northern Sinaloa, large grabens were floored by huge dome complexes at ~21-17 Ma and filled by continental sediments with interlayered basalts dated at 15-14 Ma, a setting and timing very similar to Sonora. Early to Middle Miocene volcanism, including the largely volcaniclastic Comondú strata in Baja California Sur, was thus emplaced in rift basins and was likely associated to decompression melting of upper mantle (inducing crustal partial melting) rather than to fluxing by fluids from the young subducting plate. Along the Nayarit and Sinaloa coast, flatlying basaltic lava flows dated at 11-10 Ma are exposed just above the present sea level. Here, crustal thickness is almost half that in the unextended core of the SMO, implying significant lithosphere stretching before ~11 Ma. Our study shows that rifting began much earlier than Late Miocene and provided a fundamental control on the style and composition of volcanism from at least 30 Ma. We envision a sustained period of lithospheric stretching and magmatism during which the pace and breadth of extension changed at ~20-18 Ma to be narrower and likely more rapid, and again at ~12.5 Ma, when the kinematics of rifting became more oblique.

Early-mid Cretaceous tectonic evolution of eastern Gondwana : from silicic LIP magmatism to continental rupture

The Early–mid Cretaceous marks the confluence of three major continental-scale events in eastern Gondwana: (1) the emplacement of a Silicic Large Igneous Province (LIP) near the continental margin; (2) the volcaniclastic fill, transgression and regression of a major epicontinental seaway developed over at least a quarter of the Australian continent; and (3) epeirogenic uplift, exhumation and continental rupturing culminating in the opening of the Tasman Basin c. 84 Ma. The Whitsunday Silicic LIP event had widespread impact, producing both substantial extrusive volumes of dominantly silicic pyroclastic material and coeval first-cycle volcanogenic sediment that accumulated within many eastern Australian sedimentary basins, and principally in the Great Australian Basin system (>2 Mkm3 combined volume). The final pulse of volcanism and volcanogenic sedimentation at c. 105–95 Ma coincided with epicontinental seaway regression, which shows a lack of correspondence with the global sea-level curve, and alternatively records a wider, continental-scale effect of volcanism and rift tectonism. Widespread igneous underplating related to this LIP event is evident from high paleogeothermal gradients and regional hydrothermal fluid flow detectable in the shallow crust and over a broad region. Enhanced CO2 fluxing through sedimentary basins also records indirectly, large-scale, LIP-related mafic underplating. A discrete episode of rapid crustal cooling and exhumation began c. 100–90 Ma along the length of the eastern Australian margin, related to an enhanced phase of continental rifting that was largely amagmatic, and probably a switch from wide–more narrow rift modes. Along-margin variations in detachment fault architecture produced narrow (SE Australia) and wide continental margins with marginal, submerged continental plateaux (NE Australia). Long-lived NE-trending cross-orogen lineaments controlled the switch from narrow to wide continental margin geometries.

Rapid, long-distance dispersal by pumice rafting

Pumice is an extremely effective rafting agent that can dramatically increase the dispersal range of a variety of marine organisms and connect isolated shallow marine and coastal ecosystems. Here we report on a significant recent pumice rafting and long-distance dispersal event that occurred across the southwest Pacific following the 2006 explosive eruption of Home Reef Volcano in Tonga. We have constrained the trajectory, and rate, biomass and biodiversity of transfer, discovering more than 80 species and a substantial biomass underwent a .5000 km journey in 7–8 months. Differing microenvironmental conditions on the pumice, caused by relative stability of clasts at the sea surface, promoted diversity in biotic recruitment. Our findings emphasise pumice rafting as an important process facilitating the distribution of marine life, which have implications for colonisation processes and success, the management of sensitive marine environments, and invasive pest species.

Post-orogenic Mesozoic basins and magmatism

This chapter describes the later Mesozoic history of Queensland, when the broad epicratonic basins that underlie most of the state west of the Great Dividing Range received the greater part of their sediment infill after ~210 Ma(middle Norian). The final major orogenic event—the Hunter Bowen Orogeny—had abated in the Tasmanides. These basins preserve relatively thin sedimentary successions that extend over about two-thirds of the area of the state...

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.