Large Eruptions Research Articles

A re-appraisal of the ENSO response to volcanism with paleoclimate data assimilation

The potential for explosive volcanism to affect the El Niño-Southern Oscillation (ENSO) has been debated since the 1980s. Several observational studies, based largely on tree-ring proxies, have since found support for a positive ENSO phase in the year following large eruptions. In contrast, recent coral data from the heart of the tropical Pacific suggest no uniform ENSO response to explosive volcanism over the last millennium. Here we leverage paleoclimate data assimilation to integrate both tree-ring and coral proxies into a reconstruction of ENSO state, and re-appraise this relationship. We find only a weak statistical association between volcanism and ENSO, and identify the selection of volcanic events as a key variable to the conclusion. We discuss the difficulties of conclusively establishing a volcanic influence on ENSO by empirical means, given the myriad factors affecting the response, including the spatiotemporal details of the forcing and ENSO phase preconditioning.

Limitations of assuming internal mixing between different aerosol species: a case study with sulfate geoengineering simulations

Abstract. Simulating the complex aerosol microphysical processes in a comprehensive Earth system model can be very computationally intensive; therefore many models utilize a modal approach, where aerosol size distributions are represented by observation-derived lognormal functions, and internal mixing between different aerosol species within an aerosol mode is often assumed. This approach has been shown to yield satisfactory results across a large array of applications, but there may be cases where the simplification in this approach may produce some shortcomings. In this work we show specific conditions under which the current approximations used in some modal approaches might yield incorrect answers. Using results from the Community Earth System Model v1 (CESM1) Geoengineering Large Ensemble (GLENS) project, we analyze the effects in the troposphere of a continuous increasing load of sulfate aerosols in the stratosphere, with the aim of counteracting the surface warming produced by non-mitigated increasing greenhouse gas (GHG) concentrations between 2020–2100. We show that the simulated results pertaining to the evolution of sea salt and dust aerosols in the upper troposphere are not realistic due to internal mixing assumptions in the modal aerosol treatment, which in this case reduces the size, and thus the settling velocities, of those particles and ultimately changes their mixing ratio below the tropopause. The unnatural increase of these aerosol species affects, in turn, the simulation of upper tropospheric ice formation, resulting in an increase in ice clouds that is not due to any meaningful physical mechanisms. While we show that this does not significantly affect the overall results of the simulations, we point to some areas where results should be interpreted with care in modeling simulations using similar approximations: in particular, in the evolution of upper tropospheric clouds when large amounts of sulfate are present in the stratosphere, as after a large explosive volcanic eruption or in similar stratospheric aerosol injection cases. Finally, we suggest that this can be avoided if sulfate aerosols in the coarse mode, the predominant species in these situations, are treated separately from other aerosol species in the model.

Post-caldera volcanism reveals shallow priming of an intra-ocean arc andesitic caldera: Hunga volcano, Tonga, SW Pacific

Intra-oceanic arcs are typically associated with intermediate (andesitic) cone volcanoes. However, caldera volcanoes may also form in these settings from very large eruptions, resulting in sudden changes to the magma reservoir. These reservoirs can then produce either semi-continuous or intermittent low-intensity volcanism between major caldera-producing or caldera-deepening eruptions, providing insights into the post-caldera evolution of the system. Hunga volcano (Kingdom of Tonga, Southwest Pacific) is a large mainly submarine edifice that produced a series of caldera-forming eruptions ~900 years ago. Since then, numerous smaller-scale subaerial and submarine eruptions occurred, the most recent forming new islands in 2009 and 2014/15. Pyroclastic deposits associated with these latest eruptions have identical (range ~ 0.1 wt% of all major oxides) andesitic composition that overlap with the primitive end of the slightly wider compositional range of the caldera-forming episodes. Texturally simple plagioclase, clinopyroxene and orthopyroxene phenocrysts in pre-, syn- and post-caldera pyroclasts point to a single shallow storage reservoir at 5–8 km depth. Lack of complex zonation indicates that this reservoir is constantly resupplied by low-flux inputs of basaltic andesite magma and is large enough that convective mixing rapidly homogenises new inputs. The reservoir feeds intermittent, low-intensity, post-caldera volcanism with constant andesite composition, driven possibly by magmatic overpressure and “leakage” of gas-rich magma pockets around the edges of the caldera. More primitive and compositionally variable basaltic andesites formed a lava-dominated edifice prior to the caldera-forming event. This suggests a causal link between magma supply dynamics and caldera priming relating to the maturing of the plumbing system and formation of a sustained subvolcanic andesite magma reservoir.

Clear and Present Dangers: The Multiple Health Hazards of Volcanic Eruptions.

Clear and Present Dangers: The Multiple Health Hazards of Volcanic Eruptions.

Development characteristics and petroleum geological significance of Permian pyroclastic flow volcanic rocks in Western Sichuan Basin, SW China

By examining field outcrops, drilling cores and seismic data, it is concluded that the Middle and Late Permian “Emeishan basalts” in Western Sichuan Basin were developed in two large eruption cycles, and the two sets of igneous rocks are in unconformable contact. The lower cycle is dominated by overflow volcanic rocks; while the upper cycle made up of pyroclastic flow volcanic breccia and pyroclastic lava is typical explosive facies accumulation. With high-quality micro-dissolution pores and ultra-fine dissolution pores, the upper cycle is a set of high-quality porous reservoir. Based on strong heterogeneity and great differences of pyroclastic flow subfacies from surrounding rocks in lithology and physical properties, the volcanic facies and volcanic edifices in Western Sichuan were effectively predicted and characterized by using seismic attribute analysis method and instantaneous amplitude and instantaneous frequency coherence analysis. The pyroclastic flow volcanic rocks are widely distributed in the Jianyang area. Centering around wells YT1, TF2 and TF8, the volcanic rocks in Jianyang area had 3 edifice groups and an area of about 500 km2 which is the most favorable area for oil and gas exploration in volcanic rocks.

Nutrition in Abrupt Sunlight Reduction Scenarios: Envisioning Feasible Balanced Diets on Resilient Foods.

Abrupt sunlight reduction scenarios (ASRS) following catastrophic events, such as a nuclear war, a large volcanic eruption or an asteroid strike, could prompt global agricultural collapse. There are low-cost foods that could be made available in an ASRS: resilient foods. Nutritionally adequate combinations of these resilient foods are investigated for different stages of a scenario with an effective response, based on existing technology. While macro- and micronutrient requirements were overall met, some—potentially chronic—deficiencies were identified (e.g., vitamins D, E and K). Resilient sources of micronutrients for mitigating these and other potential deficiencies are presented. The results of this analysis suggest that no life-threatening micronutrient deficiencies or excesses would necessarily be present given preparation to deploy resilient foods and an effective response. Careful preparedness and planning—such as stock management and resilient food production ramp-up—is indispensable for an effective response that not only allows for fulfilling people’s energy requirements, but also prevents severe malnutrition.

Estimating melt fraction in silicic systems using Bayesian inversion of magnetotelluric data

The location, volume and physical states of magma reservoirs are primary controls on the eruptive behavior of volcanic systems. Fundamental to understanding and monitoring these systems is the ability to identify reservoir size and physical properties, in particular melt fraction which plays an important role in the rheology and stability of a magmatic system. Large silicic volcanic eruptions in the geological record suggest that extensive pockets of melt-rich silicic magma must exist in the subsurface but such melt pockets have not been detected by geophysics. This has led to the question of whether the reservoirs that feed large volcanic eruptions are only melt-rich for a short time and thus would only be detected by geophysics shortly prior to an eruption. Magnetotelluric data measure the electrical resistivity of the subsurface and are sensitive to subsurface fluids and partial melts making it a powerful tool for imaging subvolcanic magma reservoirs. This study examines the ability for magnetotelluric data to accurately estimate melt fraction using both stochastic Bayesian inversion and deterministic regularized inversion. Results from synthetic modelling indicate that magnetotelluric data are best able to predict the melt fraction for the thick melt-rich layer using both inversion methods, though both methods under-estimate the true amount of melt. In addition, magnetotelluric data can accurately detect changes in melt fraction from crystal-rich mush (0.1 melt fraction) to melt-rich magma (0.9 melt fraction) for thick layers. Thickness is a key parameter which provides a method to assess the total volume of melt present, but it is difficult to estimate using smooth regularized inversions. Fixed-dimension Bayesian inversions provide estimates of layer thickness and their uncertainties and it is shown that estimates of total volume of melt are more accurate for both thin (0.2 km) and thick (1 km) layers than individual estimates of melt fraction or thickness alone. The melt fraction of the thin melt-rich layer is under-estimated using both the deterministic inversion method and the most probable solution of the Bayesian inversion. The Bayesian distribution of solutions has long tails which includes the true solution and better represents the uncertainty of the modelled parameters. As we consider past large melt-rich caldera eruptions at Taupo (New Zealand), Toba (Indonesia), Yellowstone and Long Valley (USA), our analysis suggests that thin melt-rich, potentially eruptible zones involved in large (e.g. >500 km3) eruptions could be systematically misinterpreted as thicker non-eruptible crystal-rich mush zones. Testing the sensitivity of MT to detect melt accumulation in active magmatic systems has important implications for volcanic hazard assessment and the use of Bayesian inversions allows for a better understanding of the range of possible interpretations.

Integrated on‐land‐offshore stratigraphy of the Campi Flegrei caldera: New insights into the volcano‐tectonic evolution in the last 15 kyr

AbstractSilicic calderas are volcanic systems whose unrest evolution is more unpredictable than other volcano types because they often do not culminate in an eruption. Their complex structure strongly influences the post‐collapse volcano‐tectonic evolution, usually coupling volcanism and ground deformation. Among such volcanoes, the Campi Flegrei caldera (southern Italy) is one of the most studied. Significant long‐ and short‐term ground deformations characterise this restless volcano. Several studies performed on the marine‐continental succession exposed in the central sector of the Campi Flegrei caldera provided a reconstruction of ground deformation during the last 15 kyr. However, considering that over one‐third of the caldera is presently submerged beneath the Pozzuoli Gulf, a comprehensive stratigraphic on‐land‐offshore framework is still lacking. This study aims at reconstructing the offshore succession through analysis of high‐resolution single and multichannel reflection seismic profiles and correlates the resulting seismic stratigraphic framework with the stratigraphy reconstructed on‐land. Results provide new clues on the causative relations between the intra‐caldera marine and volcaniclastic sedimentation and the alternating phases of marine transgressions and regressions originated by the interplay between ground deformation and sea‐level rise. The volcano‐tectonic reconstruction, provided in this work, connects the major caldera floor movements to the large Plinian eruptions of Pomici Principali (12 ka) and Agnano Monte Spina (4.55 ka), with the onset of the first post‐caldera doming at ca. 10.5 ka. We emphasise that ground deformation is usually coupled with volcanic activity, which shows a self‐similar pattern, regardless of its scale. Thus, characterising the long‐term deformation history becomes of particular interest and relevance for hazard assessment and definition of future unrest scenarios.

The 1815 Tambora eruption: Its significance to the understanding of large-explosion caldera formations

Volcanic calderas, plentiful on the Earth and the moon, have been of much interest to volcanologists because of their large dimensions and extensive volumes of ejecta. Here, we consider the dynamics of caldera-forming by major explosive eruptions, examining how the breakdown of the earth's surface is caused by violent igneous activity. This leads to the definition of “typical explosion caldera”, which is a prototype of several newly-formed calderas in the historical timescale. There are three examples of such calderas: Tambora (Sumbawa), Krakatau (Sunda Straits), and Novarupta (Alaska). Tam- bora Caldera is the best example of a well-documented, recently formed typical explosion caldera, with no significant subsequent eruptions occurring after its formation. The subsurface structure of Tambora Caldera is discussed and compared to the 1883 eruption of Krakatau, the second largest eruption in historical times. Then, contrasting with the typically basaltic “collapse-type” calderas, a “Tambora-caldera type” is defined as a large “explosion-type” caldera, that may reach up to 10 km in diameter. The Tambora- type caldera concept is useful to qualify and understand the structure and components of other major calderas in the world. Fully developed larger explosion calderas such as Aso and Aira Calderas in Kyushu, Japan are discussed and explained as composite calderas based on geophysical data. Those calderas have repeatedly ejected massive pyroclastic products causing their original structures to grow wider than 10 km.

The role of large volcanic eruptions in stratospheric ozone depletion and degradation of coniferous forests

The role of large volcanic eruptions in stratospheric ozone depletion and degradation of coniferous forests

Береговые дюны острова Уруп (Курильские острова, северо-западная Пацифика): архив изменений палеоклимата и природной среды

The Late Holocene phases of aeolian processes have been reconstructed on the basis of dune fields construction. The palaeoenvironmental studies were targeted to following problems: 1) to determine the periods of higher and lower activity of aeolian processes; 2) to establish the age of aeolian deposits and buried soils; 3) to retrace the deve­lopment of coastal landscapes. The age was determined by radiocarbon dating of paleosols and tephrostratigraphy. The correlation of tephra was performed using data on the volcanic glass chemical composition. The dunes formed during the cooling accompanied by a minor regression. Six buried soils found in the dunes reflect stabilization and overgrowing of dune fields. The longest period of dune stabilization began after a cold event 2800–2600 cal yr BP and lasted until the Little Ice Age. Paleosols contain the tephra of large volcanic eruptions on Urup (Kolokol volcano), Simushir (Zavaritsky volcano) and Iturup (tephra). Pollen analysis allows us to retrace the development of coastal landscapes. Thickets of dwarf pine developed during cooling, birch forests spread in the Medieval Warm Period, and herb meadows were widely represented on the dunes. Human impact on the coastal palaeovegetation was found. Aeolian sedimentation was high during the Little Ice Age. One of the factors of dune reactivation during the Little Ice Age was increased winter storminess associated with the East Asian winter monsoon. Evidence of active cyclogenesis is the increasing proportion of allochthonous pollen. The modern reactivation of aeolian processes is associated with human activity and storm erosion of dunes.

Evolution of the Popocatépetl Volcanic Complex: constraints on periodic edifice construction and destruction by sector collapse

Popocatépetl is one of the most active volcanoes in North America. Its current predominantly mild activity is contrasted by a history of large effusive and explosive eruptions and sector collapse events, first summarized by Espinasa-Pereña and Martín-Del Pozzo (2006). Since then, a wealth of new radiometric, geophysical and volcanological data have been published, requiring a re-evaluation of the evolution of the Popocatépetl Volcanic Complex (PVC). We combined existing literature data with new field observations, aerial imagery and digital elevation model interpretations to produce an updated and improved reconstruction of the growth and evolution of the PVC throughout its history. This will be fundamental for the assessment and mitigation of risks associated with potential future high-magnitude activity of the PVC. The PVC consists of four successive volcanic edifices separated by three sector collapse events producing avalanche deposits: Tlamacas (>538–>330 ka, described here for the first time); Nexpayantla ( c. 330 to >98 ka); Ventorrillo ( c. 98–23.5 ka); and Popocatépetl (<23.5 ka). The newly described Tlamacas collapse propagated towards the ENE, forming part of the Mayorazgo avalanche deposit. Supplementary material: Supplementary files 1–7 are available at https://doi.org/10.6084/m9.figshare.c.5709190

Large silicic magma bodies and very large magnitude explosive eruptions

Over the last 20 years, new concepts have emerged into understanding the processes that lead to build up to large silicic explosive eruptions based on integration of geophysical, geochemical, petrological, geochronological and dynamical modelling. Silicic melts are generated within magma systems extending throughout the crust by segregation from mushy zones. Segregated melt layers become unstable and can assemble into ephemeral upper crustal magma chambers rapidly prior to eruption. In the next 10 years, we can expect major advances in dynamical models as well as in analytical and geophysical methods, which need to be underpinned in field research.

Nature, timing and magnitude of buried Late Cretaceous magmatism on the central West Iberian Margin

AbstractThe magma‐poor West Iberian Margin (WIM), as part of the Peri‐Atlantic alkaline province, records multiple evidence for intra‐plate post‐rift magmatism. Based on high‐resolution multichannel seismic data, this work discusses the presence of large volcanic and intrusive features in the Estremadura Spur, providing evidence for important magmatic activity during the drifting of the continental margin. Our observations reveal distinct voluminous fissure‐fed effusive sequences and the details of the 2800 m‐high Fontanelas compound volcano, including its external and internal architecture, secondary vents and associated lava flows, all of which were probably extruded at intermediate water depths. Numerous and morphologically diverse sills and sill complexes are also described, attesting to the presence of a Late Cretaceous shallow magmatic plumbing system in the area. Magmatism in this region is interpreted as having occurred during two main pulses and types of activity: (1) Coniacian to lower Campanian(?) age, characterised by fissural and fault‐controlled volcanism, which mostly extruded massive lobate/sheet lava flows; and (2) a second voluminous intrusive and extrusive event of mid to late Campanian age, which includes the intrusion of the Estremadura Spur laccolith and the prominent Fontanelas compound volcano with associated dendritic lava flows. The inferred volumes of the first fissure‐fed effusive event suggest a large eruption magnitude, comparable to some of the largest historical effusive eruptions. The second magmatic pulse led to the emplacement of discrete clusters of sills and sill complexes, as well as the construction of the ca. 2.8‐km‐high Fontanelas volcano, suggesting a syn‐rift structural inheritance that controlled the location of the Estremadura Spur Intrusion and the Fontanelas volcanic area. Altogether, a total volume of rock exceeding 1.452 km3 is estimated to have been emplaced or extruded in this region in a relatively short period, attesting to the prominence of the magmatism in this sector of the WIM.

Magnetic Expression of Hydrothermal Systems Hosted by Submarine Calderas in Subduction Settings: Examples from the Palinuro and Brothers Volcanoes

Volcanism is the most widespread expression of cyclic processes of formation and/or destruction that shape the Earth’s surface. Calderas are morphological depressions resulting from the collapse of a magma chamber following large eruptions and are commonly found in subduction-related tectono-magmatic regimes, such as arc and back-arc settings. Some of the most impressive examples of seafloor hydrothermal venting occur within submarine calderas. Here, we show the results of magnetic investigations at two hydrothermally active submarine calderas, i.e., Palinuro Seamount in the Southern Tyrrhenian Sea, Italy, and Brothers volcano of the Kermadec arc, New Zealand. These volcanoes occur in different geodynamic settings but show similarities in the development of their hydrothermal systems, both of which are hosted within calderas. We present a new integrated model based on morphological, geological and magnetic data for the Palinuro caldera, and we compare this with the well-established model of Brothers caldera, highlighting the differences and common features in the geophysical expressions of both hydrothermal systems. For consistency with the results at Brothers volcano, we build a model of demagnetised areas associated with hydrothermal alteration derived from 3D inversion of magnetic data. Both these models for Brothers and Palinuro show that hydrothermal up-flow zones are strongly controlled by caldera structures which provide large-scale permeability pathways, favouring circulation of the hydrothermal fluids at depth.

Implications of Multiple Disequilibrium Textures in Quartz-Hosted Embayments

The faces of volcanic phenocrysts may be marked by imperfections occurring as holes that penetrate the crystal interior. When filled with glass these features, called embayments or reentrants, have been used to petrologically constrain magmatic ascent rate. Embayment ascent speedometry relies on the record of disequilibrium preserved as diffusion-limited volatile concentration gradients in the embayment glass. Clear, glassy embayments are carefully selected for speedometry studies. The use and subsequent descriptions of pristine embayments overrepresent their actual abundance. Here, we provide a textural analysis of the number, morphology, and filling characteristics of quartz-hosted embayments. We target a collection of large (i.e., >20 km3 erupted volume) silicic eruptions, including the Bishop Tuff, Tuff of Bluff Point, Bandelier Tuff, Mesa Falls Tuff, and Huckleberry Ridge Tuff in the United States, Oruanui Tuff in New Zealand, Younger Toba Tuff in Indonesia, the Kos Plateau Tuff in Greece, and the Giant Pumice from La Primavera caldera in Mexico. For each unit, hundreds of quartz crystals were picked and the total number of embayment-hosting crystals were counted and categorized into classifications based on the vesicularity and morphology. We observed significant variability in embayment abundance, form, and vesicularity across different eruptions. Simple, cylindrical forms are the most common, as are dense glassy embayments. Increasingly complex shapes and a range of bubble textures are also common. Embayments may crosscut or deflect prominent internal cathodoluminescence banding in the host quartz, indicating that embayments form by both dissolution and growth. We propose potential additional timescales recorded by embayment disequilibrium textures, namely, faceting, bubbles, and the lack thereof. Embayment formation likely occurs tens to hundreds of years before eruption because embayment surfaces are rounded instead of faceted. Bubble textures in embayments are far from those predicted by equilibrium solubility. Homogenous nucleation conditions likely allow preservation of pressures much greater than magmastatic inside embayments. Our textural observations lend insight into embayment occurrence and formation and guide further embayment studies.

Deciphering Cryptic Multi-Stage Crystal-Melt Separation during Construction of the Tonglu Volcanic–Plutonic Complex, SE China

Abstract Revealing the origin of explosive eruptions of silica-rich magma is of paramount importance for understanding the evolution of continental crust and volcanic hazards. However, it remains controversial as to how the erupted magmas form and how they connect with plutonic realms, partly owing to the rarity and the obscurity of the ‘cumulate’ complementary to such eruptions of silica-rich magmas. Here the issues are explored by comparing the volcanic rocks (rhyodacite + rhyolite) and their associated subvolcanic intrusions (monzodiorite + monzonite + quartz monzonite) within Tonglu volcanic basin (SE, China). The Tonglu plutonic and volcanic units are consistent with each other in age (~130 Ma), space and source materials (e.g. Sr-Nd-Hf isotopes), strongly suggesting that they are cogenetic. Mineral mapping demonstrates that abundant plagioclase clusters (and chains in monzodiorite) occur in plutonic units, implying the processes of crystal gathering and/or accumulation. Rhyolite-MELTS modeling, and geochemical studies coupled with textural observations suggest that the Tonglu volcanic and plutonic rocks represent the residual melts and the complementary cumulate residues, respectively. The compositional and mineral variations in the plutonic rocks can be explained by two-stage, low-pressure crystal-melt separation of a dacitic magma. The monzodiorite represents the first-stage cumulate that was unsaturated in zircon and biotite/K-feldspar on the basis of low Zr and Ba concentrations and the occurrence of these two phases in the interstices between plagioclase and amphibole. The monzonite and quartz monzonite are the second-stage cumulates after saturation of zircon and biotite/K-feldspar as indicated by abrupt increases in Zr and Ba concentrations and zircon inclusions within euhedral biotite. Mass balance calculation and textural estimation indicate that the cumulates actually are a mixture of crystals and melt, containing ~40 vol% interstitial melt. Interstitial zircons from monzodiorite, largely crystallized from trapped melt, show contrasting trace-element trends (e.g. Ti, Zr/Hf, Eu/Eu*, Gd/Yb) to those of the other rock types (i.e. monzonite, quartz monzonite, rhyodacite and rhyolite), consistent with fractional crystallization (FC). We interpret these divergences to heating-induced partial dissolution of a basal crystal framework of monzodiorite due to recharges of hot mafic magmas, but the limited dissolution fails to rejuvenate the crystal mush. The Tonglu volcanic–plutonic system demonstrates that compositional distillation via crystal–liquid separation within the upper crust is an effective mechanism by which a potentially eruptible rhyolitic cap can be generated. Recharges of hot mafic magmas at shallow crustal levels may prolong the lifespan of granitic magma reservoirs. A combined study including texture, geochemistry, thermodynamic modeling and mass-balance calculations can help us identify the fingerprints of cumulates in felsic magma systems and thus track the processes responsible for producing large eruptions of silica-rich magmas.

Seasonal Changes in Arctic Cooling After Single Mega Volcanic Eruption

To investigate the pure long-term influence of single mega volcanic eruption (SMVE) of universal significance on Arctic temperature changes in summer and winter, the Samalas eruption in Indonesia which is the largest eruption over the past millennium is selected as an ideal eruption for simulation study based on Community Earth System Model. The significant Arctic cooling lasts for 16 years after the Samalas eruption. The obvious Arctic cooling shifts from summer to winter, and this seasonal change of cooling after the SMVE only exists in the high-latitude Arctic region. The cooling range in Arctic summer is larger than that in winter during the first 2 years, due to the strong weakening effect of volcanic aerosol on summer incident solar radiation and the snow-ice positive feedback caused by the rapid expansion of summer sea ice, while the winter sea ice in the same period doesn’t increase obviously. Starting from the third year, the Arctic winter cooling is more intense and lasting than summer cooling. The direct weakening effect of aerosol on solar radiation, which is the main heat source in Arctic summer, is greatly weakened during this period, making summer cooling difficult to sustain. However, as the main heat source in Arctic winter, the sea surface upward longwave radiation, sensible heat, and latent heat transport still maintain a large decrease. Furthermore, sea ice expansion and albedo increase result in the decrease in solar radiation and heat absorbed and stored by the ocean in summer. And the isolation effect of sea ice expansion on air-sea heat transfer in winter during this period makes the heat transfer from the ocean to the atmosphere correspondingly reduce in winter, thus intensifying the Arctic winter cooling. Additionally, the Arctic Oscillation (AO) changes from the negative phase to the positive phase in summer after the SMVE (such as Samalas), while it is reversed in winter. This phase change of AO is also one of the reasons for the seasonal changes in Arctic cooling.

Magma Chamber Formation by Dike Accretion and Crustal Melting: 2D Thermo‐Compositional Model With Emphasis on Eruptions and Implication for Zircon Records

AbstractWe present a 2D thermo‐elastic model of a magma body formation in a granitic crust by injection of rhyolitic or basaltic dikes and sills. Elastic analytical solutions enable the computation of rock displacement in response to magma intrusion and evacuation during volcanic eruptions. Phase diagrams for magma and rocks govern melting/crystallization behavior and temperature evolution in 16 million computational cells. Calculated temperature histories are used to predict zircon crystallization/dissolution, their ages, and isotopic ratios within individual batches of magma and rocks. Incremental dike injection naturally generates magma batches of melt that form kilometer‐wide clusters appearing in different parts of the domain with diverse shapes, horizontal and vertical interconnectivity. The volume of eruptive melt strongly depends on magma influx rates Q, the width of the injection region W, and eruptions. For example, rhyolitic dike injection with Q = 0.125 m3/s with W = 5 km during 100 ka generates ∼50 km3 of eruptive melt while no significant melt forms if W = 10 km. Injection of basaltic dikes into the granitic crust generates comparable amounts of rhyolitic melt from dike‐residual and country rocks. High injection intensities produce magma reservoirs capable of large eruptions, while repetitive eruptions lead to shrinkage of magma bodies. High heat input causes host rock zircons to lose a significant portion of their old cores. Magmatic zircons in the periphery form quickly due to rapid cooling of intruded dikes while in the central part crystals can grow during several hundreds of ka. The model predicts highly heterogeneous O and Hf isotopic ratios recorded in zircons.

A latest Pleistocene and Holocene composite tephrostratigraphic framework for northeastern North America

Lakes and bogs in northeastern North America preserve tephra deposits sourced from multiple volcanic systems in the Northern Hemisphere. However, most studies of these deposits focus on specific Holocene intervals and the latest Pleistocene, providing snapshots rather than a full picture. We combine new data with previous work, supplemented by a broad review of the characteristics and ages of potential source regions and volcanoes, to develop the first composite tephrostratigraphic framework covering the last ∼14,000 years for this region. We report new cryptotephra records from three ombrotrophic peat bogs—Irwin Smith (Michigan), Bloomingdale (New York), and Sidney Bog (Maine)—as well as new analyses and age models from previously reported sites, Nordan's Pond Bog (Newfoundland) and Thin-Ice Pond (Nova Scotia). A new tephra (Iliinsky) from the NGRIP and GRIP ice cores is also presented as it can be correlated to new data from these terrestrial records and helps validate radiocarbon age models. We identify 21 new tephra in addition to the 15 already known, several of which cover the entire region – the White River Ash east, Newberry Pumice, Ruppert (NDN-230), and Mazama. For the first time we find Mount St. Helens Yn (ca. 3660 cal yr BP) and a set P tephra (∼3000–2550 cal yr BP), and confirm the presence of Jala Pumice from Volcan Ceboruco, Mexico, and KS1 from Ksudach volcano, Kamchatka. We describe new “ultra-distal” tephra, including the early Holocene KS2 eruption, and propose correlations to volcanoes Iliinsky and Shiveluch of Kamchatka, and Ushishir of the Kurile Islands. Not all of these tephra represent large eruptions, with several plausible correlations to sub-Plinian events. Using Bayesian age-modeling, we present new age estimates for the newly described tephra, for tephra with previously poor age control, and for several proximal correlatives. Overall, we demonstrate northeastern North America's importance for providing transcontinental linkages between paleoenvironmental records and providing insights into ash distribution from different styles and sizes of eruptions.