Ignimbrite Flare-up Research Articles

The Neogene–Quaternary volcanism of the Carpathian–Pannonian region: from initial plate tectonic models to quantitative petrogenetic explanations

A wide range of volcanic rocks formed in the Carpathian–Pannonian region during the last 20 Myr, closely associated with the tectonic evolution of the area: (1) high-SiO 2 rhyolites and rhyodacites; (2) calc-alkaline-type volcanic suites (basalt–andesite–dacite–rhyolite and adakitic rocks); (3) potassic–ultrapotassic volcanic rocks (shoshonites, latites, leucitites and lamproites); and (4) alkali sodic volcanic rocks (basalts, trachyandesites and trachytes). During the 18.1–14.4 Ma silicic ignimbrite flare-up volcanism, explosive eruptions produced large volume of tephras, which covered extended areas of Europe. Calc-alkaline-type volcanism took place in syn-extensional and post-collisional settings. There is no classic volcanic arc in the Carpathian–Pannonian Region, and so the term ‘back-arc basin’ should be used for the Pannonian Basin only in a structural geology sense. Alkali basalt volcanism occurred mostly in the peripheral areas of the Pannonian Basin primarily owing to edge-driven convective upwelling of heterogeneous asthenospheric mantle. The youngest eruptions took place 100–30 ka ago, mostly in the geodynamically still active southeast Carpathians. After 50 years of intense research, there are still many open questions to be solved. Interdisciplinary understanding of the results of different disciplines is crucial to better constrain the evolution of this tectonically complex area, to reveal its geoenergy and to assess potential natural hazards.

Permian-Triassic volcanic and plutonic records of the Argentine Frontal Cordillera: A review with new U–Pb and Hf-isotope zircon data

The Permian-Triassic magmatism of western Argentina and Chile represents one of the most outstanding silicic magmatic events of the southwestern Gondwana margin, notably marked by the development of the Choiyoi Magmatic Province (CMP). We provide a comprehensive review of its volcanic and plutonic record in the Argentine Frontal Cordillera. The volcanic rocks form three distinct sequences. The oldest is depicted by the Las Lozas volcanic sequence of the northern Frontal Cordillera for which new U–Pb zircon data (288 ± 2 Ma and εHft values ranging from −3.97 to +0.73) reassigns these outcrops to the early Cisuralian, aligning with volcanic records of northern Chile (ca. 297-288 Ma). The middle sequence, deposited during the late Cisuralian-late Guadalupian interval, is ascribed to the Choiyoi Group, which is characterised by a transition from andesitic (ca. 280-270 Ma) to rhyolitic compositions (ca. 270-262 Ma), including a remarkable mid-Guadalupian ignimbrite flare-up event (ca. 265 Ma). The upper sequence, composed of andesite-dacite-rhyolite, is associated with the Guanaco Sonso basin situated in the westernmost region, deposited from the late Lopingian to the middle Triassic period (ca. 254-240 Ma). The volcanic successions were developed in an extensional/transtensional setting, with facies and thickness variations controlled by normal faults, some of them active during the eruption of caldera-forming ignimbrites. Regarding the plutonic component, it comprises over sixty granitoid bodies forming the Colangüil and El Portillo batholiths and scattered stocks throughout the Frontal Cordillera. The early-stage plutons (ca. 285-272 Ma) exhibit calc-alkaline tonalite-diorite to granodiorite-monzogranite compositions and overlap in age with the lower andesitic section of the Choiyoi Group, while the late-stage plutons (ca. 265-252 Ma) display syenodiorite and alkali granite compositions and are in most cases younger than the Choiyoi Group succession. Compiled U–Pb zircon geochronological data reveal a distinctive Permian-age phase characterized by a rapid expansion of magmatism from the Gondwana margin towards its interior, followed by a slower westward shift of the main magmatic belt, predominantly recorded in the Frontal Cordillera of Argentina.

Onset and tempo of ignimbrite flare-up volcanism in the eastern and central Mogollon-Datil volcanic field, southern New Mexico, USA

The Cenozoic ignimbrite flare-up (40–18 Ma) generated multiple volcanic fields in the southwestern United States and northern Mexico resulting from asthenospheric mantle upwelling after removal of the Farallon slab. The correlation of tuffs to one another and to source calderas within these volcanic fields is essential for determining spatiotemporal patterns in volcanism and magma geochemistry, which have been used to deduce migration of the Farallon slab at depth and associated mantle melting. However, the correlation of Eocene–Oligocene tuffs in the southwestern U.S. is difficult because of post-emplacement erosion and faulting. This study focuses on spatiotemporal patterns of the initial episode of ignimbrite flare-up activity (ca. 36.5–33.8 Ma) in the Mogollon-Datil volcanic field in south-central New Mexico, USA. We show that alkali feldspar major and trace element geochemistry is an effective tool for correlating tuffs when combined with high-precision, single-crystal 40Ar/39Ar geochronology and bulk-rock geochemistry. Using these data, we correlate several tuff units and differentiate other tuffs that have the same eruption age but very different geochemistry, and we conclude that there was a broadly northwestward migration in volcanism over time. The new tuff correlations are used to investigate spatiotemporal variations in magma geochemistry, erupted volumes, and recurrence intervals during the initial episode of Mogollon-Datil volcanic field volcanism. Early-erupted tuffs restricted to the eastern Mogollon-Datil volcanic field share similarities with western U.S. topaz rhyolites, which suggests that the silicic magmas were generated by partial melting of mafic lower crustal rocks. We also find differences in the compositions, crystallinities, and mineral assemblages between the early- and late-erupted tuffs. The early-erupted tuffs tend to have single-feldspar mineralogies, lower feldspar Or contents, large negative Eu anomalies, and low-whole–rock Ba concentrations. Conversely, late-erupted tuffs have two feldspar plus quartz assemblages, lesser Eu anomalies, higher whole-rock Ba concentrations, and feldspars have higher Or contents. Thus, we suggest that for some of the early eruptions, after magmas underwent crystal fractionation in the crust, the silicic melt largely separated from the crystalline mush prior to eruption, whereas late-erupted tuff magmas underwent crystal fractionation at near the eutectic minimum and were remobilized and erupted with a larger proportion of their crystalline mush. Using our new ages, correlations, and previously published data, we find that the initial phase of Mogollon-Datil volcanic field volcanism produced at least 15 eruptions between 36.5 Ma and 33.8 Ma, with a minimum total erupted volume of ~1350 km3 and an average recurrence interval of 170 k.y. However, eruptions were generally smaller in volume (most <15 km3) than in other coeval fields, and most eruptions (n = 11) occurred in the first 1.2 m.y. (ca. 36.5–35.3 Ma) of activity. Altogether, our work sheds new light on variations in the composition, timing, and migration of volcanism during the initial phase of Mogollon-Datil volcanic field activity and highlights the utility of feldspar geochemistry in both “fingerprinting” tuffs and elucidating magma evolution.

High-precision U[sbnd]Pb geochronology of an ignimbrite flare-up in a silicic large igneous province (Chon Aike, Patagonia)

Silicic large igneous provinces (SLIPs) are periods of particularly voluminous felsic volcanism in the geologic record. Previous work has suggested an overall long lifespan for SLIPs of 20 to 40 Myr, commonly punctuated by shorter-lived ‘flare-ups’ of higher volcanic productivity (1–5 Myr), but detailed studies of individual flare-up events are lacking. The Jurassic Chon Aike SLIP (CASP) is the product of an exceptional geological event wherein voluminous felsic volcanism (ca. 219,000 km3) was generated predominately via crustal anatexis over 45 Myr. We focus on the Late Jurassic El Quemado Complex (EQC), which marked the final stages of felsic volcanism for the CASP. In-situ U-Pb ages for the EQC previously suggested a duration of ∼5 Myr; however, new high-precision CA-ID-TIMS ages indicate the total durations of ignimbrite successions were shorter than 350 kyr. A compilation of zircon U-Pb ages for the eight CASP formations in Patagonia and the Antarctic Peninsula reveals changes in volcanic duration between formations deposited in the intraplate relative to the continental margin, suggesting a spatial control over the magmatic lifespans of geographically restricted systems. We suggest that the rate of magmatism in the CASP was primarily controlled by heterogeneities in the crust, largely between Proterozoic igneous crust and younger, metasedimentary crust that was recently accreted relative to the timing of Jurassic volcanism. The observed duration of volcanism was modulated by these differences in crustal properties from the injection of mafic magmas into the lower crust during extension and mantle upwelling, following rollback of the subducting oceanic slab towards the Paleo-Pacific margin of Gondwana. The short duration in the EQC is the product of unique overlapping conditions that favored crustal melting and resulted in the voluminous ignimbrite flareup (104 km3) of some of the Earth's highest δ18O magmas measured.

Precursors to a continental-arc ignimbrite flare-up: Early central volcanoes of the San Juan Mountains, Colorado, USA

Abstract Our newly acquired and recently published map, geochronologic, and compositional data for early intermediate-composition central volcanoes in the northeastern San Juan Mountains provide insights about the broad magmatic precursors to the large continental-arc ignimbrite flare-up in the mid-Cenozoic Southern Rocky Mountain volcanic field (SRMVF). Initial volcanism migrated from central Colorado to northern New Mexico ca. 40–29 Ma, as part of a more regional trend of southward-progressing mid-Cenozoic magmatism in the U.S. segment of the North American Cordillera. Within the San Juan locus, which represents the largest preserved erosional remnant of the SRMVF and site of most intense eruptive activity, new 40Ar/39Ar and U-Pb zircon ages show that eruptions at many individual edifices began nearly concurrently, at ca. 35 Ma, with peak activity at 34–32 Ma. Broadly similar precursor effusive volcanism characterizes other major loci of continental-arc ignimbrite magmatism along the western American cordilleras, but none of these sites records early volcanism as voluminous, spatially widespread, well exposed, or compositionally diverse as the San Juan locus in Colorado. Early San Juan volcanism was larger in volume than the later ignimbrites, constituting about two thirds the total erupted. Early lava and breccias are as much as 700–900 m thick where exposed along eroded flanks of the San Juan Mountains; drill holes have penetrated sections as thick as 2600 m. The early volcanoes were dominantly andesitic, with lesser dacite and minor rhyolite. Such volcanism is widely interpreted as initiated by basaltic magma from the mantle, but mid-Cenozoic basalt is almost nonexistent at the San Juan locus—an absence inferred to be due to extensive lower-crustal assimilation and fractionation. The early volcanic rocks are calc-alkaline and typical of high-K continental-arc volcanism; they become modestly more alkalic and enriched in trace elements such as light rare earth elements, Zr, Nb, and Th from the San Juan locus northeastward into central Colorado. Such variations may reflect synmagmatic crustal thickening and deeper levels of primary magma generation concurrent with assembly of upper-crustal magma bodies that could support large ignimbrite eruptions. Substantial uncertainties remain for growth histories of the early volcanoes, however, because of unexposed lower parts of edifices, eroded upper parts, and limited availability of mineral phases that could be dated reliably. Although the early volcanoes are widely distributed within the SRMVF, many are clustered at sites of subsequent ignimbrite calderas. The precursor edifices are inferred to record incubation stages in construction of overall translithospheric batholithic-scale magmatic systems. Prolonged processes of incremental magma generation, accumulation, fractionation, and solidification intermittently generated sufficient liquid to erupt large ignimbrites. Maturation of focused eruptions and intrusions was prolonged, 5 m.y. or more, prior to the culminating ignimbrite at some centers in the San Juan Mountains. Some large-volume ignimbrites and related calderas, including the ~5000 km3 Fish Canyon Tuff and associated La Garita caldera, formed as much as several million years later than peak growth of associated precursor volcanoes, recording a sustained interval of diminished eruptive activity as the magma reservoir increased in volume and evolved to more silicic compositions capable of supporting a subsequent large ignimbrite eruption. Dike configurations at early volcanoes that were active nearly concurrently in the SRMVF vary from symmetrically radial to more parallel trends. The contrasting dike geometries are inferred to record possible multiple fluctuations from compressive to weakly extensional regional stress, concurrent with destabilization of the prior flat-slab plate configuration that triggered mid-Cenozoic ignimbrite flare-ups along the Cordilleran margin of the North American plate. These apparent fluctuations in regional stress preceded development of substantial extensional strain in the Southern Rocky Mountain region; outflow ignimbrite sheets of the SRMVF spread across subsequent horst-and-graben structures of the Rio Grande rift without complementary thickness variations.

Connecting the Dots: the Lava Domes' Perspective of Magmatism Related to an Ignimbrite Flare-Up

Abstract Dome eruptions associated with rhyolitic calderas offer an important insight into how extremely large (> > 10 km3), rhyolitic magma systems are constructed through time. We focus on rhyolitic calderas in the central Taupō Volcanic Zone leading to, during, and immediately following the 350- to 240-ka ignimbrite flare-up. We identified 103 dome eruptions that are dated between ca. 650 and 150 ka and collated 239 literature whole-rock compositions from these domes. For each composition, we modelled the pressure of magma extraction from the magma mush and the mineral assemblage of the mush using the rhyolite-MELTS geobarometer. We calculated extraction temperatures using zircon saturation geothermometry. We show that magmas are extracted from typically quartz-bearing magma mush at a wide range of depths (~50–425 MPa, ~2–16 km) and temperatures (~750°C to ~850°C). Throughout the central TVZ, there are two dominant extraction pressure modes at 1) 150–175 MPa and 2) 250–325 MPa, consistent with 1) the depth of the brittle–ductile transition (~6 km) and just below typical pre-eruptive storage depths of other TVZ magmas (100–150 MPa, ~4–6 km); and 2) partial melt regions imaged below ~8 km by previous geophysical studies. In some regions, there is a clear correlation between crustal structures, the depth of magma extraction, and the composition of the magmas. In the Whakamaru caldera, the domes erupted inside the caldera following caldera collapse are extracted from ~225 to ~350 MPa at ~810°C and have orthopyroxene-bearing compositions dissimilar to the caldera-forming eruption. These domes are aligned along normal faults, suggesting that rifting creates pathways for magma extraction from a deeper mush rejuvenated by recharge. The domes erupted along the structural margins of the Whakamaru caldera have very evolved, hornblende-bearing compositions, similar to the caldera-forming eruption and shallow, colder extraction from ~100 to ~200 MPa at ~770°C, suggesting the mush feeding these domes is a remnant of the older caldera-forming magma system mobilized along the caldera-bounding faults. Two structural levels of magma extraction at ~6 km and 9 to 12 km are persistent throughout the flare-up period and across the central TVZ region, demonstrating the need for further investigation into the factors controlling the depth of mush development.

Flare Up of Hot-Dry-Reduced Ignimbrites Related to Extension in the Cascades Arc: The Deschutes Formation, Central Oregon

Abstract Ignimbrite flare-ups are rare periods of intense silicic volcanism during which the pyroclastic volume and eruptive frequency is more than an order of magnitude higher than background activity. Investigating the compositional differences between flare-up and steady-state magmas provides critical constraints on the petrogenetic causes for the event and can offer unique opportunities to investigate the role of large-scale tectonic or geodynamic processes in arc magmatism. In this study, we focus on the bimodal Deschutes Formation ignimbrite flare-up of Central Oregon, which erupted unusually high volumes of pyroclastic material 6.25–5.45 Ma from a new axis of volcanism in the Cascades arc. This episode is marked by increased eruption rates and eruption of more silicic compositions relative to the Quaternary Cascade arc, which rarely erupts rhyolites. Ignimbrites are crystal-poor (<10%) dacite to rhyolites (mostly 65–77 wt.% SiO2) with anhydrous mineral assemblages and higher FeO/MgO, Y, Eu/Eu*, MREE and Zr/Sr, indicating drier magmatic evolution compared to the Quaternary arc, and are more similar to those from the rear-arc High Lava Plains (HLP) province that lies to the east. Magnetite-ilmenite oxybarometry indicates that Deschutes Formation felsic magmas tend to be hotter and more reduced (NNO-1 to NNO) than the Quaternary arc (NNO to NNO + 1.5). Rhyolite-MELTS geobarometry suggests complex storage of diverse Deschutes Formation magmas within the shallow crust (50–250 MPa), and the common co-eruption of multiple plagioclase populations, pumice compositions, and compositionally banded pumice suggest variable degrees of mixing and mingling of distinct magmas. Deschutes magmas also have low δ18Oplagioclase values that indicate partial melting and assimilation of hydrothermally altered shallow crust. Trace element systematics and rhyolite-MELTS modeling suggests that felsic pumice cannot be produced by simple fractionation of co-erupted mafic pumice or basaltic lavas, and requires a crustal melting origin, and trace elements and Pb isotopes suggest that young mafic crust may have been the primary protolith. We suggest that partial melting produced low-Si rhyolite melt (~72 wt.%) that acted as both a parent for the most evolved rhyolites, and as a mixing endmember to create the dacite to rhyodacite magmas with heterogenous plagioclase populations. Unlike the predominantly calc-alkaline basalts erupted in the Quaternary Cascade arc, Deschutes Formation primary basalts are mostly low-K tholeiites, indicative of decompression melting. These are similar to the compositions erupted during a contemporaneous pulse of low-K tholeiite volcanism across the whole HLP that reached into the Cascades rear-arc. We suggest that intra-arc extension focused decompression melts from the back-arc into the arc and that tensional stresses allowed this high flux of hot-dry-reduced basalt throughout the crustal column, causing partial melting of mafic protoliths and the production of hot-dry-reduced rhyolite melts. Depletion of incompatible elements in successive rhyolites implies progressive depletion in fertility of the protolith. Extension also allowed for the establishment of a robust hydrothermal system, and assimilation of hydrothermally-altered rocks by magmas residing in a shallow, complex storage network lead to low δ18O melts. Our findings suggest the integral role that extensional tectonics played in producing an unusual ignimbrite flare-up of hot-dry-reduced rhyolite magmas that are atypical of the Cascades arc and may be an important contributor to flare-ups at arcs worldwide.

Geochemical constraints on Eocene ignimbrite flare-up in the Urumieh-Dokhtar magmatic arc, northwestern Iran

Geochemical constraints on Eocene ignimbrite flare-up in the Urumieh-Dokhtar magmatic arc, northwestern Iran

Following magma: The pathway of silicic magmas from extraction to storage during an ignimbrite flare-up, Taupō Volcanic Zone, New Zealand

Very large silicic magma systems can generate tens to thousands of cubic kilometres of very evolved, rhyolitic magmas that erupt in caldera-forming events. These magma systems may incorporate multiple bodies of magma mush, from which eruptible, crystal-poor magma is extracted and stored. We examine the pressure of magma extraction from the mush and storage of eruptible magma bodies from seven caldera-forming eruptions in the central Taupō Volcanic Zone, New Zealand, during a 350 ka to 240 ka ignimbrite flare-up. For the first time, we use rhyolite-MELTS geobarometry to calculate extraction and storage pressures for the same pumice clasts (representing magma parcels), revealing that each eruption was fed by multiple, compositionally distinct, mush and eruptible magma bodies. Some of the eruptible bodies were stored contiguously with the magma mush; but, on average, our calculations reveal that individual parcels of eruptible magma were stored ∼6 km shallower than their extraction depths, demonstrating that magmas migrate an appreciable vertical distance in the build-up to eruption. The flare-up includes three different geometries of magmatic system: 1) a “mature” type in which multiple eruptible magma bodies sit adjacent to each other in the shallow crust (∼100 MPa) and are in contact with the top of a vertically extensive, heterogeneous mush system; 2) a “mush model” type in which a shallow eruptible magma body is directly in contact with a shallow mushy cumulate; and 3) a “vertically discretised” type in which eruptible magma bodies are stored at different depths and may also be extracted from separate mush bodies. The flare-up began with a “mature” type system, with the very large (>2200 km3 magma) Whakamaru Group eruptions. Following Whakamaru, the eruptions were smaller (erupting 50 – 150 km3 magma), and initially had deep, “vertically discretised” magma systems. The subsequent eruptions reveal the vertical expansion of a new system, until the last eruption (Ohakuri) had a “mature”-type organisation. This sequence of flare-up eruptions hints at how very large silicic magma systems can spatially and vertically mature over a geologically short timeframe (∼100 kyrs) from a patchwork of vertically discretised magma bodies to a mature magma system with a shallow storage zone and compositionally heterogeneous magma mush beneath.

Andesites and evolution of the continental crust: Perspectives from the Central Volcanic Zone of the Andes

Named for the Andes, andesites (53%–63% SiO2) are the archetypal magma erupted at magmatic arcs. They have been established as the average composition of continental crust and as such are integral to the growth and evolution of the continental crust. However, andesites are quite variable in trace element and isotopic composition reflecting disparate paths of origin. Herein we return to the original site of their identification, the Central Andes, and use a comprehensive dataset of published and unpublished trace elements and isotopes to show that during the past 6 Myr two distinct types of andesite have erupted in the Central Volcanic Zone (CVZ), which correspond with different geodynamic conditions. Consistent with previous work, we confirm that major composite cones and minor centers of the steady state (low magmatic flux) Quaternary CVZ arc have trace element and isotopic characteristics consistent with magma generation/fractionation in the lower crust. Within the Quaternary arc centers, there are also significant latitudinal variations that correspond with the age, composition, and P-T conditions of the lower crust. However, in contrast to this prevailing model, in the 21–24°S segment 6–1 Ma andesites from ignimbrites and lava domes associated with the peak of the regional Neogene ignimbrite flare-up have compositions that indicate these andesites are hybrids between mantle-derived basalts and upper crustal lithologies. Since ∼1 Ma, andesites in young silicic lava domes associated with the regional flare-up are compositionally indistinguishable from proximal Quaternary arc centers, indicating a return to steady-state magmatism and lower crustal production of andesites. We propose that the transition from upper crustal to lower crustal andesite production results from a decrease in mantle heat input and subsequent relaxation of the regional geotherm during the waning of the flare-up event. The two modes of andesite production have significant implications for the production and evolution of the CVZ arc crust. During the flare-up, prodigious amounts of basalt were emplaced into the mid-crust, resulting in the production of large volumes of hybrid intermediate magmas in the mid and upper crust. In contrast, the lower crustal differentiation recorded in the Quaternary steady state arc andesites would result in the formation of a dense crystalline residue in the lower crust and an overall densification of the lower crust. Over time, gravity instabilities associated with this densification may ultimately aid in the delamination of the dense lower crustal root, triggering flare-ups. These differences in andesite production may help explain the cyclicity (flare-up cycles) observed in mature continental arcs and emphasizes that andesite is not a monotonous composition and can vary with depth-dependent intra-crustal differentiation related to magmatic flux.

Tracing widespread Early Miocene ignimbrite eruptions and petrogenesis at the onset of the Carpathian-Pannonian Region silicic volcanism

The Carpathian-Pannonian Region (CPR) hosted some of the largest silicic volcanic eruptions in Europe during the Early and Middle Miocene, contemporaneously with major lithospheric thinning of the Pannonian Basin. This was recorded as an ignimbrite flare-up event from approximately 18.1–14.4 Ma. To gain in-depth perspectives on the eruption chronology, tephrostratigraphy, and petrogenesis at the onset of CPR silicic volcanism, we applied a multi-proxy approach to Lower Miocene rhyolitic ignimbrites and pyroclastic fall deposits from the northern CPR to the Dinaride Lake System. High-precision zircon U-Pb geochronology distinguished two Lower Miocene groups of volcaniclastic rocks at ∼ 18.1 Ma and ∼ 17.3 Ma. Based on combined tephrostratigraphic signatures we propose that the ∼ 18.1 Ma Kalnik and ∼ 17.3 Ma Eger eruptions produced widespread (intermediate to) large caldera-forming massive rhyolitic ignimbrites, deposited across northern and southwestern regions of the CPR. Due to easterly winds that carried volcanic ash hundreds of kilometers to the southwest, Eger eruption products also reached distal intra-montane Dinaride lacustrine basins, recorded as pyroclastic fall deposits. Heterogeneous major and trace elemental compositions of ∼ 18.1 Ma volcanic glass shards suggest that the Kalnik eruption was sourced from complex silicic magmatic systems, with simultaneous tapping of two discrete melt bodies during the eruption. The homogeneous geochemical composition of ∼ 17.3 Ma glasses is distinct from the older glasses. Integrated zircon and bulk glass Nd-Hf isotope compositions have a positive correlation, defining a regional mantle array, and are more radiogenic in the younger phase of volcanism. The recorded systematic isotopic change, moving from older more crustal signatures to younger more juvenile compositions, imply that during the period of lithospheric thinning of the Pannonian Basin the region underwent more complex variations in the interaction between metasomatized lithospheric mantle-derived magmas and various crustal components than previously recognized.

Fluorine-rich mafic lower crust in the southern Rocky Mountains: The role of pre-enrichment in generating fluorine-rich silicic magmas and porphyry Mo deposits

AbstractFluorine-rich granites and rhyolites occur throughout the southern Rocky Mountains, but the origin of F-enrichment has remained unclear. We test if F-enrichment could be inherited from ancient mafic lower crust by: (1) measuring amphibole compositions, including F and Cl contents, of lower crustal mafic granulite xenoliths from northern Colorado to determine if they are unusually enriched in halogens; (2) analyzing whole-rock elemental and Sr, Nd, and Pb isotopic compositions for upper crustal Cretaceous to Oligocene igneous rocks in Colorado to evaluate their sources; and (3) comparing batch melting models of mafic lower crustal source rocks to melt F and Cl abundances derived from biotite data from the F-rich silicic Never Summer batholith. This approach allows us to better determine if the mafic lower crust was pre-enriched in F, if it is concentrated enough to generate F-rich anatectic melts, and if geochemical data support an ancient lower crustal origin for the F-rich rocks in the southern Rocky Mountains.Electron microprobe analyses of amphibole in lower crustal mafic granulite xenoliths show they contain 0.56–1.38 wt% F and 0.45–0.73 wt% Cl. Titanium in calcium amphibole thermometry indicates that the amphiboles equilibrated at high to ultrahigh temperature conditions (805 to 940 °C), and semiquantitative amphibole thermobarometry indicates the amphiboles equilibrated at 0.5 to 1.0 GPa prior to entrainment in magmas during the Devonian. Mass balance calculations, based on these new measurements, indicate parts of the mafic lower crust in Colorado are at least 3.5 times more enriched in F than average mafic lower crust. Intrusions coeval with the Laramide Orogeny (75 to 38 Ma) pre-date F-rich magmatism in Colorado and have Sr and Nd isotopic compositions consistent with mafic lower crust ± mantle sources, but many of these intrusions contain elevated Sr/Y ratios (>40) that suggest amphibole was a stable phase during magma generation. The F-rich igneous rocks from the Never Summer igneous complex and Colorado Mineral Belt also have Sr and Nd isotopic compositions that overlap with the lower crustal mafic granulite xenoliths, but they have lower Sr/Y, higher Nb and Y abundances, and distinctly less radiogenic 206Pb/204Pbi compositions than preceding Laramide magmatism. Batch melt modeling indicates low-degree partial melts derived from rocks similar to the mafic lower crustal xenoliths we analyzed can yield silicic melts with >2000 ppm F, similar to estimated F melt concentrations for silicic melts that are interpreted to be parental to evolved leucogranites.We suggest that F-rich silicic melts in the southern Rocky Mountains were sourced from garnet-free mafic lower crust, and that fluid-absent breakdown of amphibole in ultrahigh temperature metamorphic rocks was a key process in their generation. Based on the composition of high-F amphibole measured from lower crustal xenoliths, the temperature of amphibole breakdown and melt generation for these F-enriched source rocks is likely >100 °C higher than similar lower crust with low or average F abundances. As such, these source rocks only melted during periods of unusually high heat flow into the lower crust, such as during an influx of mantle-derived magmas related to rifting or the post-Laramide ignimbrite flare-up in the region. These data have direct implications for the genesis of porphyry Mo mineralization, because they indicate that pre-enrichment of F in the deep crust could be a necessary condition for later anatexis and generation of F-rich magmas.

Geochronology and glass geochemistry of major Pleistocene eruptions in the Main Ethiopian Rift: Towards a regional tephrostratigraphy

The Main Ethiopian Rift (MER) is renowned as a focus of investigations into human origins. It is also the site of many large volcanic calderas, whose eruptions have spanned the timeframe of speciation, cultural innovation, and dispersal of our species. Yet, despite their significance for dating human fossils and cultural materials, the timing and geochemical signatures of some of the largest eruptions have remained poorly constrained at best. Here, through a programme of field surveys, geochemical analysis and 40Ar/39Ar dating, we report the ages of MER ignimbrites and link them to widespread tephra layers found in sequences of archaeological and paleoenvironmental significance. We date major eruptions of Fentale (76 ± 18 ka), Shala (ca. 145–155 ka), Kone (184 ± 42 ka and ca. 200 ± 12 ka) and Gedemsa (251 ± 47 ka) volcanoes, and correlate a suite of regionally important tephra horizons. Geochemical analysis highlights the predominantly peralkaline rhyolitic melt compositions (7.5–12 wt% Na2O + K2O, 70–76 wt% SiO2) across the central MER and remarkable similarity in incompatible trace element ratios, limiting the correlation of deposits via glass composition alone. However, by integrating stratigraphic and geochronological evidence from proximal deposits, lake sediment cores and distal outcrops at archaeological sites, we have traced ash layers associated with the ca. 177 ka Corbetti, ca. 145–155 ka Shala and ca. 108 ka Bora-Baricha-Tullu-Moye eruptions across southern Ethiopia. In addition to strengthening the tephrochronological framework that supports paleoenvironmental and archaeological work in the region, our findings have wider implications for evaluating the hypothesis of a middle Pleistocene ‘ignimbrite flare-up’ in the MER, and for evaluating the impacts of these great eruptions on landscapes, hydrology, and human ecology.

Large-magnitude (VEI ≥ 7) ‘wet’ explosive silicic eruption preserved a Lower Miocene habitat at the Ipolytarnóc Fossil Site, North Hungary

During Earth’s history, geosphere-biosphere interactions were often determined by momentary, catastrophic changes such as large explosive volcanic eruptions. The Miocene ignimbrite flare-up in the Pannonian Basin, which is located along a complex convergent plate boundary between Europe and Africa, provides a superb example of this interaction. In North Hungary, the famous Ipolytarnóc Fossil Site, often referred to as “ancient Pompeii”, records a snapshot of rich Early Miocene life buried under thick ignimbrite cover. Here, we use a multi-technique approach to constrain the successive phases of a catastrophic silicic eruption (VEI ≥ 7) dated at 17.2 Ma. An event-scale reconstruction shows that the initial PDC phase was phreatomagmatic, affecting ≥ 1500 km2 and causing the destruction of an interfingering terrestrial–intertidal environment at Ipolytarnóc. This was followed by pumice fall, and finally the emplacement of up to 40 m-thick ignimbrite that completely buried the site. However, unlike the seemingly similar AD 79 Vesuvius eruption that buried Pompeii by hot pyroclastic density currents, the presence of fallen but uncharred tree trunks, branches, and intact leaves in the basal pyroclastic deposits at Ipolytarnóc as well as rock paleomagnetic properties indicate a low-temperature pyroclastic event, that superbly preserved the coastal habitat, including unique fossil tracks.

Timescales and thermal evolution of large silicic magma reservoirs during an ignimbrite flare-up: perspectives from zircon

Four voluminous ignimbrites (150–500 km3) erupted in rapid succession at 27 Ma in the central San Juan caldera cluster, Colorado. To reconstruct the timescales and thermal evolution of these magma reservoirs, we used zircon ID-TIMS U–Pb geochronology, zircon LA-ICP-MS geochemistry, thermal modeling, and zircon age and crystallization modeling. Zircon geochronology reveals dispersed zircon age spectra in all ignimbrites, with decreasing age dispersion through time that we term a ‘chimney sweeping’ event. Zircon whole-grain age modeling suggests that 2σ zircon age spans represent approximately one-quarter of total zircon crystallization timescales due to the averaging effect of whole-grain, individual zircon ages, resulting in zircon crystallization timescales of 0.8–2.7 m.y. Thermal and zircon crystallization modeling combined with Ti-in-zircon temperatures indicates that magma reservoirs were built over millions of years at relatively low magmatic vertical accretion rates (VARs) of 2–5 × 10–3 m y−1 (2–5 × 10–6 km3 y−1 km−2), and we suggest that such low VARs were characteristic of the assembly of the greater San Juan magmatic body. Though we cannot unequivocally discern between dispersed zircon age spectra caused by inheritance (xenocrystic or antecrystic) versus prolonged crystallization from the same magma reservoir (autocrystic), our findings suggest that long-term magma input at relatively low VARs produced thermally mature upper crustal magma reservoirs resulting in protracted zircon crystallization timescales. Compiling all U–Pb ID-TIMS zircon ages of large ignimbrites, we interpret the longer timescales of subduction-related ignimbrites as a result of longer term, lower flux magmatism, and the shorter timescales of Snake River Plain ignimbrites as a result of shorter term, higher flux magmatism.

Magmatism and Extension in the Anaconda Metamorphic Core Complex of Western Montana and Relation to Regional Tectonics

AbstractMetamorphic core complexes (MCCs) are a product of crustal extension, but their dynamics are still debated. Early research suggests that the formation of MCCs in the western United States was due to gravitational collapse of crust that had been thickened during Cordilleran orogenesis. However, the instability of overthickened crust alone cannot explain the diachronous formation of core complexes with a strong spatial dependency, as there was relatively uniform crustal thickness along strike of the Cordillera. For this reason, there is an interest in what role other lithospheric processes (such as subducted slab removal) play in the evolution of MCCs. We investigate the role of such processes by determining the temporal relation between magmatism and extension in the Anaconda MCC (AMCC) of western Montana. Geologic mapping, zircon U‐Pb geochronology, and zircon (U‐Th)/He thermochronology reveal that the initiation of extension in the AMCC in the Eocene (∼53 Ma) began at least 3 Myr after the emplacement of voluminous Paleocene two‐mica plutons. We interpret that the AMCC is an example of a core complex that was primed for extension by magmatic thermal weakening and suggest that foundering of the Farallon flat slab and the onset of the ignimbrite flareup in western Montana was responsible for the initiation of AMCC extension. An updated compilation of MCC cooling ages and Cenozoic volcanic activity across the western United States supports previous interpretations that the removal of Farallon oceanic lithosphere likely initiated MCC exhumation in some regions.

Tephrochronology of North America's most recent arc-sourced ignimbrite flare-up: The Deschutes Formation of the Central Oregon Cascades

The dangerous juxtaposition of human populations and active volcanic arcs worldwide makes elucidating the history of explosive eruptions of these arcs a crucial endeavor. The well-preserved Deschutes Formation, which forms part of the Deschutes Basin of central Oregon, USA, contains a distal record of a period of unusually frequent explosive eruptions between 5.45 and 6.25 Ma, which is atypical of the Cascades arc as a whole. However, no previous study has established the chronology, frequency, or total volume of all of the ignimbrites and tephra-fall units exposed within the Formation in order to compare to baseline activity of the arc. We use multivariate statistics on major and trace element compositions of pumice glass to correlate pyroclastic units across the Deschutes Basin, and combine these with recently published 40Ar/39Ar dates to establish, for the first time, a detailed and comprehensive tephrochronology of the Deschutes Formation. Our results suggest that at least 78 statistically distinct explosive eruptions occurred between 6.25 and 5.45 Ma. Eruption frequency is not constant throughout the 800 kyr of pyroclastic volcanism; a maximum rate of 1.7 to 3.0 explosive eruptions per 10 kyr occurred between 5.76 and 5.68 Ma and decreased for 340 kyr prior to graben formation. We estimate that the volume for 51 ignimbrite eruptions is 290 to 565 km3 (170 to 280 km3 DRE), assuming an equal westward flow component and a fall:flow ratio of 0.5:1 to 1.9:1. Using a single isopach method, we also estimate that a minimum volume of 110 km3 (45 km3 DRE) was deposited as 27 pumice-fall units which have no correlated ignimbrites. Thus, we estimate that a cumulative pyroclastic volume of 400 to 675 km3 (210 to 330 km3 DRE) was erupted in the central Oregon Cascades during this time. The volumetric rate (2.7 to 4.1 km3/10 kyr, DRE) and average frequency of ignimbrite eruptions (1.0 to 1.5 eruptions/ 10kyr) are more than a factor of ten greater than the Quaternary Cascades arc. Thus, the Deschutes Formation records North America's most recent arc-sourced ignimbrite flare-up. Although smaller in volume than most well-recognized flare-up events, this represents an important category of flare-ups that may be more common in arcs than previously recognized. We hypothesize that a heightened flux of basalt, possibly induced by slab-rollback, was focused beneath the arc and into the shallow crust by minor amounts of crustal extension. This extension allowed for the high flux of basalt to be stored at shallow levels beneath a new arc locus within fertile crust, resulting in an unusually silicic and explosive period of Cascade arc history.

Magmatic Evolution of Zoned and Unzoned Ignimbrites: Evidence for a Complex Crustal Architecture Feeding four Rapid-sequence, Caldera-forming Eruptions in the San Juan Mountains, Colorado

Abstract The last four caldera-forming ignimbrites in the central San Juan caldera cluster, Colorado, erupted 1400 km3 in ≤80 kyr and alternated between zoned crystal-poor rhyolite to crystal-rich dacite and unzoned, crystal-rich dacite. The zoned 150 km3 Rat Creek Tuff (26·91 Ma), unzoned 250 km3 Cebolla Creek Tuff, and zoned 500 km3 Nelson Mountain Tuff (26·90 Ma) formed the nested San Luis caldera complex with slightly offset calderas, and the unzoned 500 km3 Snowshoe Mountain Tuff (26·87 Ma) formed the Creede caldera to the south. The Rat Creek Tuff, Nelson Mountain Tuff, and Snowshoe Mountain Tuff have similar mineral assemblages of plagioclase, sanidine, quartz, biotite, hornblende, clinopyroxene, Fe–Ti oxides, and accessory zircon, titanite, and apatite. The Cebolla Creek Tuff differs from the other three ignimbrites with more abundant hornblende and a lack of quartz and sanidine. Trace element compositions of interstitial glass are unique to each ignimbrite, correlating with mineral assemblages and inferred crystallization depths. Glass, feldspar, hornblende, and clinopyroxene thermobarometry calculations provide evidence for vertically extensive crustal magma reservoirs with a common magmatic zone at ∼435–470 MPa (∼16–17 km) showing a transition into shallow pre-eruptive reservoirs between ∼110 and 340 MPa (∼4–13 km), similar to the estimated magma reservoir architecture of the Altiplano Puna Volcanic Complex. The upper portions of the eruptible parts of the magma reservoirs of the Rat Creek Tuff (215 ± 50 MPa, ∼810–820 °C), Cebolla Creek Tuff (340 ± 20 MPa, ∼860–880 °C), Nelson Mountain Tuff (215 ± 20 MPa, ∼745–800 °C), and Snowshoe Mountain Tuff (110 ± 40 MPa, 825 ± 10 °C) occupied shallow levels in the crust similar to other magma reservoirs of the central San Juan caldera cluster. Trace element modelling correlates with a deep crystallization signature in the unzoned Cebolla Creek Tuff and Snowshoe Mountain Tuff, typified by a flat trend in Ba versus Sr whole-rock and glass chemistry. The zoned Rat Creek Tuff and Nelson Mountain Tuff are typified by a steep trend in Ba versus Sr chemistry interpreted as a shallower crystallization signature. Similarly, the unzoned Cebolla Creek Tuff and Snowshoe Mountain Tuff have flatter slopes in FeO versus An space of plagioclase chemistry interpreted as more abundant deep plagioclase crystallization and a difficulty in physically mixing with Fe-rich mafic recharge magma owing to higher viscosity. The zoned Rat Creek Tuff and Nelson Mountain Tuff have higher slopes in FeO versus An space of plagioclase chemistry interpreted as more abundant shallow plagioclase crystallization and more feasible mixing with Fe-rich mafic recharge magma owing to lower viscosity. The eruption of the Rat Creek Tuff was probably triggered by mafic injection, but the other three ignimbrites lack mingling textures in pumice, suggesting that other mechanisms were important in causing such large eruptions. After a prolonged period of mantle-derived magma injection and crustal heating (∼25 000 km3 Conejos Formation erupted during ∼35–29 Ma), the San Juan magmatic body became a robust and versatile producer of diverse eruptible magmas in short time periods during its Oligocene ignimbrite flare-up.

Spatio-Temporal Shifts in Magmatism and Mineralization in Northern Colorado Beginning in the Late Eocene

AbstractMagmatism in northern Colorado beginning in the late Eocene is associated with the formation of Pb-Zn-Ag carbonate-replacement and polymetallic vein deposits, the onset of caldera-forming magmatism, and eventually, the formation of rift-related, F-rich Mo porphyries (“Climax-type” intrusions). We use high-precision U/Pb zircon geochronology to better evaluate the temporal framework of magmatism and mineralization in the region. Our results demonstrate that mineralization in the Leadville area occurred between 43.5 and 39.7 Ma and was followed by mesothermal mineralization in the Montezuma area at approximately 38.7 Ma. Mineralization is associated with a suite of approximately 43 to 39 Ma intermediate magmatic centers that extended from Twin Lakes through Montezuma. The oldest porphyries associated with F-rich Mo prospects and deposits (Middle Mountain; 36.45 Ma) intruded 900 kyr after the start of the ignimbrite flare-up in the region. Spatiotemporal analyses reveal that the pattern of magmatism shifted in orientation between 40 and 35 Ma. We propose a model wherein magmatism before 39 Ma was the result of fluids evolved from the subducted Farallon slab being focused through weak zones in the lithospheric mantle and into the lower crust. This was followed by a more diffuse and higher power melting event that corresponds to a distinct change in the spatial patterns of magmatism. Our data suggest that low-grade Mo porphyry deposits can form close in time to calderas. We hypothesize that the transition from subduction to extensional tectonics in the region was responsible for this more widespread melting and a distinct shift in the style of magmatic-hydrothermal mineralization.

Zircon geochronology and O-Hf isotopes of Cappadocian ignimbrites: New insights into continental crustal architecture underneath the Central Anatolian Volcanic Province, Turkey

During the Mio–Pliocene between ca. 10 and 5 Ma, the Central Anatolian Volcanic Province (CAVP) formed as the focus of a high eruptive flux episode with voluminous high-silica, caldera-forming eruptions. Magmatic compositions of ignimbrites emplaced during this episode were uniformly rhyodacitic–rhyolitic, reflecting melt generation in an intracrustal hot zone heated by intrusions of basalt that were extracted from fertile asthenospheric mantle undergoing decompression. The clear imprint of crustal rocks on CAVP magmas provides unique opportunity to constrain the deep crustal architecture of central Anatolia. Here, we present a detailed study of zircon from key CAVP ignimbrites including U-Pb ages and the first coupled O- and Hf-isotopic compositions for this silicic igneous province, along with complementary whole-rock Nd isotope data. Ignimbrites with higher δ18O and lower εHf values (Zelve, Cemilköy, Gördeles) are typically compositionally more evolved compared to those with intermediate (Kavak, Kızılkaya) and more mantle-like compositions (Sarımadentepe, Sofular, Tahar), which is consistent with magmas acquiring their isotopic compositions through crustal assimilation coupled with fractional crystallization. Intriguingly, O-Hf-zircon compositions from Kızılkaya ignimbrite, one of the largest and youngest CAVP eruptions, cover almost the same isotopic range than all CAVP ignimbrites collectively. Despite evidence for a significant crustal influence revealed by O-Nd-Hf isotopic compositions of CAVP ignimbrites, inherited zircon is absent in all units. This is in stark contrast with previous studies of Kızılkaya ignimbrite for which abundant Proterozoic–Archean zircon populations have been reported. Isotopic compositions of evolved magmas in the CAVP can be explained by interaction of mantle-derived basalt with post-Archean crust, with minor contributions of hydrothermally altered rocks. No tangible evidence exists for a magmatic contribution of Archean crustal components to the Mio–Pliocene ignimbrite flare-up in Central Anatolia, suggesting that regional basement rocks largely represent Gondwanan crust that was originally consolidated and subsequently became reprocessed during post-Archean times.