Description: Highlights • Masaya caldera is an unusual basaltic caldera in that it formed by voluminous magma extraction during explosive eruptions. We identify the nature, age and volume of these three eruptions of which the first, emplacing the San Antonio tephra, was by far most voluminous. • The by far largest fraction of the 9 km3 DRE erupted volume of this tephra was discharged during a Plinian eruption phase, which was bracketed by phreatomagmatic eruptions. We demonstrate that water contents measured in melt inclusions equilibrated during residence at shallow level shortly before eruption strongly underestimate original water contents during differentiation at higher pressure. We argue that the large fraction of exsolved H2O together with buoyancy pressure from connection to the deeper reservoir drove the eruptive high mass flux needed for the Plinian eruption phase. Masaya is unusual for a basaltic caldera because it formed by piston-subsidence in response to large-volume magma withdrawal by highly explosive eruptions, i.e. in a fashion typical of silicic calderas. The first and most voluminous of the three explosive eruptions formed the 6 ka old basaltic San Antonio Tephra (SAT). This eruption is also unusual in that most of the 9 km3 DRE basaltic magma was discharged by a plinian eruption. The subsequent eruptions of the basaltic Masaya Triple Layer (MTL, 2.1 ka) and the Masaya Tuff/Ticuantepe Lapilli (MT-TIL, 1.9 ka) each discharged 2 km3 DRE magma and enlarged the Masaya caldera. The SAT consists of a lower sequence of alternating scoria lapilli and ash layers, interpreted as an alternation between more or less phreatomagmatically influenced fallout events. These are followed by two prominent well-sorted lapilli layers: the first one formed by a climactic plinian eruption whose column height reached 21–29 km and discharged most of the total erupted mass including about 35 Mt. SO2. The second, lithic-rich lapilli layer probably formed by a phreatoplinian event when partial collapse of the magma chamber roof initiated increasing magma-water interaction which ultimately formed the upper sequence of phreatomagmatic cross-bedded surge deposits, accretionary lapilli-rich tuffs and a final fallout of dense lapilli. Phreatomagmatic activity may have been related to disruption of a hydrothermal system reflected in hydrothermally altered lithics, and/or by the caldera floor subsiding closer to the groundwater table. The bulk-rock chemical composition of the SAT is basaltic but the bimodal glass compositions demonstrate mixing of a basaltic with an andesitic melt probably in the conduit during eruption. The SAT basalt differentiated in a reservoir near the MOHO at 20 km depth by fractional crystallization of olivine, plagioclase, and minor clinopyroxene forming a tholeiitic fractionation trend. Minor intermediate-An plagioclase crystallized from the basaltic melt at H2O concentrations of about 2 wt% as measured by FTIR in melt inclusions. However, a key observation is that the melt inclusions are not in equilibrium with the high-An plagioclases hosting them. Re-equilibration of the inclusions requires initially higher water contents (about 5–6 wt%) which also fits the high Ba/La ~ 80 indicating input from the strongly hydrated subducting slab. Therefore, while the SAT magma evolved under hydrous conditions at depth, it was then stored at shallow level long enough to adjust to the low saturation pressure and to precipitate some intermediate-An plagioclase but still preserving its high temperature (around 1100 °C) and phenocryst-poor composition. Large overpressure due to connection to the deep-seated reservoir and water degassing during ascent limited the storage time at shallow level and drove the unusually intense and voluminous plinian-style eruption that facilitated piston-type collapse of the chamber roof.

Description: Highlights The increase in volcanic activity after the last glacial maximum observed on Iceland has led to one of the most fascinating hypothesis in science in the last decades: that deglaciation may force volcanism. We: - Re-analyzed four longer tephra records with the same statistical method and demonstrated that all contain the ∼41 kyr and ∼100 kyr Millankovitch periodicities. - The frequency spectra of the tephra and δ18O records are significantly correlated supporting the hypothesis that orbital-driven global climate changes interact with the volcanic eruption frequency regionally and globally. - However, the simultaneous analysis of the four best-characterized tephra records shows that correlations and associated time lags suffer from a number of uncertainties including the nature and quality of tephra time series, a wide range in geographic latitudes and geological settings, as well as applied statistical methods Therefore more precise tephra time series (preservation and age optimized) from different regions (glaciated versus non-glaciated) and geological settings (island arcs, continental arcs, intraplate) are needed together with standardized statistical analysis to decipher the impact of these factors on a global perspective of how climate may control volcanism. Abstract The increase in volcanic activity after the last glacial maximum observed on Iceland has led to one of the most fascinating hypothesis in science in the last decades: that deglaciation may force volcanism. Consequently, tephrostratigraphic records of sufficient length that cover multiple glacial cycles have been used to test whether such relationships hold systematically through the Quaternary. Here we review such tephra records that have been linked with climate proxy records such as δ18O in marine sediments, which is a measure of sea-level change and which is thought to be orbitally forced, as it exhibits the characteristic Milankovitch periodicities of precession (∼23 kyr), obliquity (∼41 kyr) and eccentricity (∼100 kyr). Statistical analyses have identified these periodicities also in long tephra records from different latitudes and geotectonic settings, as well as in compiled semi-global records. These studies detect Milankovitch periods in their tephra record, and also a phase shift relative to the δ18O record in such that periods of increased eruption frequencies coincide with the deglaciation period at the glacial/interglacial transition when ice and water loads on the lithosphere change most rapidly. However, there are also disparities in results and interpretations, which may be attributable to the different methods of analysis applied by the studies. We have therefore re-analyzed the four best-characterized tephra records by the same methods. We distinguish between analysis in the frequency domain, a novel approach, and analysis in the time domain, which has been used in previous studies. Analysis in the frequency domain identifies harmonic frequencies that arise from the binary nature of the tephra records and complicate the identification of primary frequencies. However, we show that all four records show spectral density peaks near the main Milankovitch periodicities of 41 and 100 kyr, and that they produce meaningful and significant statistical correlations with each other and the global δ18O record but not with random time series. Although the time-domain correlations with δ18O roughly confirm phase shifts implying peak volcanism during deglaciation, correlation coefficients arising from very noisy records are generally too low for precise constraints on the relative timing. These deficiencies presently hamper the recognition of the physical mechanisms through which global climate changes affect volcanism at both, high-latitude glaciated regions and low-latitude non-glaciated regions.

Description: Highlights • Major eruption of Ilopango volcano, El Salvador occurred in the first half of the 6th century. • Ilopango eruption is consistent with ‘mystery’ eruption of 540 CE that caused global cooling. • Magnitude 7 event ranks as one of the 10 largest on Earth in past 7000 years. • Impacts on the Maya of Central America were severe, including estimated 100,000 + fatalities. Abstract Ilopango volcano (El Salvador) erupted violently during the Maya Classic Period (250–900 CE) in a densely-populated and intensively-cultivated region of the southern Maya realm, causing regional abandonment of an area covering more than 20,000 km2. However, neither the regional nor global impacts of the Tierra Blanca Joven (TBJ) eruption in Mesoamerica have been well appraised due to limitations in available volcanological, chronological, and archaeological observations. Here we present new evidence of the age, magnitude and sulfur release of the TBJ eruption, establishing it as one of the two hitherto unidentified volcanic triggers of a period of stratospheric aerosol loading that profoundly impacted Northern Hemisphere climate and society between circa 536 and 550 CE. Our chronology is derived from 100 new radiocarbon measurements performed on three subfossil tree trunks enveloped in proximal TBJ pyroclastic deposits. We also reassess the eruption magnitude using terrestrial (El Salvador, Guatemala, Honduras) and near-shore marine TBJ tephra deposit thickness measurements. Together, our new constraints on the age, eruption size (43.6 km3 Dense Rock Equivalent of magma, magnitude = 7.0) and sulfur yield (∼9–90 Tg), along with Ilopango's latitude (13.7° N), squarely frame the TBJ as the major climate-forcing eruption of 539 or 540 CE identified in bipolar ice cores and sourced to the tropics. In addition to deepening appreciation of the TBJ eruption's impacts in Mesoamerica, linking it to the major Northern Hemisphere climatic downturn of the mid-6th century CE offers another piece in the puzzle of understanding Eurasian history of the period.

Description: The stratigraphic succession of widespread tephra layers in west-central Nicaragua was emplaced by highly explosive eruptions from mainly three volcanoes: the Chiltepe volcanic complex and the Masaya and Apoyo calderas. Stratigraphic correlations are based on distinct compositions of tephras. The total tephras combine to a total on-shore volume of about 37 km3 produced during the last ∼ 60 ka. The total erupted magma mass, including also distal volumes, of 184 Gt (DRE) distributes to 84% into 9 dacitic to rhyolitic eruptions and to 16% into 4 basaltic to basaltic–andesitic eruptions. The widely dispersed tephra sheets have up to five times the mass of their parental volcanic edifices and thus represent a significant albeit less obvious component of the arc volcanism. Eruption magnitudes (M = log10(m) − 7 with m the mass in kg), range from M = 4.1 to M = 6.3. Most of the eruptions were dominantly plinian, with eruption columns reaching variably high into the stratosphere, but minor phreatomagmatic phases were also involved. Two phreatomagmatic eruptions, one dacitic and one basaltic–andesitic, produced mostly pyroclastic surges but also fallout from high eruption columns. Comparison of fallout tephra dispersal patterns with present-day, seasonally changing height-dependant wind directions suggests that 8 eruptions occurred during the rainy season while 5 took place during the dry season. The tephra succession documents two major phases of erosion. The first phase, 〉 17 ka ago, appears to be related to tectonic activity whereas the second phase may have been caused by wet climatic conditions between 2 to 6 ka ago. The Apoyo caldera had two large plinian, caldera-forming eruptions in rapid succession about 24 ka ago and should be considered a silicic volcano with long repose times. Three highly explosive basaltic eruptions were generated at the Masaya Caldera within the last 6 ka. Since then frequent but small eruptions and lava effusion were largely limited to the caldera interior. The dacitic Chiltepe volcanic complex experienced six plinian eruptions during the last 17 ka and seems to be an accelerating system in which eruption magnitude increased while the degree of differentiation of erupted magma decreased at the same time. We speculate that the Chiltepe system might produce the next large-magnitude silicic eruption in west-central Nicaragua.

Description: Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition. Deposition from these PDCs began only on slopes at or around the material's angle of repose (c. 30°). In unconfined settings, the granular PDCs are interpreted to have been quasi-steady, forming sheets and lobes around the angle of repose. Where flows were confined, sheet-like proximal facies made up around 10% of the deposit volume at the angle of repose, but 90% of the material was deposited from apparently unsteady inertial granular PDCs as a distal levée-and-channel facies on slopes well below the repose angle. Hence, confined PDCs were able to travel up to 50% farther than unconfined flows. In the distal facies the deposit width is inversely correlated to the local slope, and the height of the levées (above the deposit centreline) is positively correlated with slope. Internally the deposits comprise three parts, a coarse-grained fines-free sole layer that laterally connects to levées (Zone I), an ashy matrix-supported central body (Zone II) and an overlying coarse plaster of clasts (Zone II). Trends in grain-size data suggests these zones derive from a continuous un-mixing of coarse particles from the initial bulk material by granular segregation that preferentially drives large particles to the upper free surface of the flow where they are concentrated at the front of flow before being deposited and overrun. By comparison to analogue experiments, we suggest a model of flow and deposition where the temporally and spatially varying mode of deposition is determined by the flow velocity, the local slope, the vertical velocity gradient, the velocity gradient at the free surface and the vertical deposition rate. Using this model, estimated vertical deposition rates of c. 5 cm s− 1 from the Ngauruhoe PDCs agree with those determined in laboratory experiments on inertial granular flows.

Description: The Masaya Caldera Complex has been the site of three highly explosive basaltic eruptions within the last six thousand years. A Plinian eruption ca. 2 ka ago formed the widespread deposits of the Masaya Triple Layer. We distinguish two facies within the Masaya Triple Layer from each other: La Concepción facies to the south and Managua facies to the northwest. These two facies were previously treated as two separated deposits (La Concepción Tephra and the Masaya Triple Layer of Pérez and Freundt, 2006) because of their distinct regional distribution and internal architectures. However, chemical compositions of bulk rock, matrix and inclusion glasses and mineral phases demonstrate that they are the product of a single basaltic magma batch. Additionally, a marker bed containing fluidal-shaped vesicular lapilli allowed us to make a plausible correlation between the two facies, also supported by consistent lateral changes in lithologic structure and composition, thickness and grain size. We distinguish 10 main subunits of the Masaya Triple Layer (I to X), with bulk volumes ranging between 0.02 and 0.22 km3, adding up to 0.86 km3 (0.4 km3 DRE) for the entire deposit. Distal deposits identified in two cores drilled offshore Nicaragua, at a distance of ∼ 170 km from the Masaya Caldera Complex, increase the total tephra volume to 3.4 km3 or ∼ 1.8 km3 DRE of erupted basaltic magma. Isopleth data of five major fallout subunits indicate mass discharges of 106 to 108 kg/s and eruption columns of 21 to 32 km height, affected by wind speeds of 〈 2 m/s to ∼ 20 m/s which increased during the course of the multi-episodic eruption. Magmatic Plinian events alternated with phreatoplinian eruptions and phreatomagmatic explosions generating surges that typically preceded breaks in activity. While single eruptive episodes lasted for few hours, the entire eruption probable lasted weeks to months. This is indicated by changes in atmospheric conditions and ash-layer surfaces that had become modified during the breaks in activity. The Masaya Triple Layer has allowed to reconstruct in detail how a basaltic Plinian eruption develops in terms of duration, episodicity, and variable access of external water to the conduit, with implications for volcanic hazard assessment.

Description: Defining a precise timeline for past eruptions from explosive volcanoes in continental arcs is imperative to forecast future hazards and mitigate volcanic disasters in these often densely populated regions. However, establishing reliable ages for Quaternary eruptions in the Central American Volcanic Arc has been challenging due to the common lack or alteration of suitable K-rich phases for 40Ar/39Ar geochronology, but also from their position in time beyond the reach of 14C dating. This especially holds for the active Amatitlán caldera in Guatemala, from which at least six explosive silicic eruptions have produced tephra blanketing neighboring regions that are today inhabited by millions of people. Zircon, a common datable accessory mineral in Amatitlán caldera magmas, is used here to retrieve eruption ages by applying the novel zircon double-dating method (ZDD) that integrates 238U–230Th disequilibrium dating and (U–Th)/He thermochronology. This approach yielded the first-ever radioisotopic ages of 24 ± 3 ka and 48 ± 6 ka (1σ), respectively, of two of Amatitlán caldera's most recent eruptions (J-tephra and E-tephra). Remarkably, both zircon crystallization and ZDD eruption ages for the older and voluminous T-tephra and L-tephra units significantly post-date existing plagioclase 40Ar/39Ar dates by ca. 26 and 70 kyr, respectively. The ZDD eruption age for T-tephra is 93 ± 4 ka, whereas zircon crystallization ages for L-tephra yield a maximum model eruption age of ca. 124 ka. The strong eruption age divergence between ZDD and plagioclase 40Ar/39Ar dating argues for the presence of inherited or xenocrystic plagioclase in Amatitlán caldera eruptive products. Statistical analysis based on the updated eruptive history suggests a recurrence interval of ca. 17 kyr, which is significantly shorter than previously estimated. The new age data, thus, suggest a more frequent eruptive activity of Amatitlán caldera than formerly thought and underscores the necessity to better understand the current underlying magmatic system and to constrain its past eruptive history more precisely.