On Dec. 22, 2018, at approximately 20:55-57 local time, Anak Krakatau volcano, located in the Sunda Straits of Indonesia, experienced a major lateral collapse during a period of eruptive activity ...that began in June. The collapse discharged volcaniclastic material into the 250 m deep caldera southwest of the volcano, which generated a tsunami with runups of up to 13 m on the adjacent coasts of Sumatra and Java. The tsunami caused at least 437 fatalities, the greatest number from a volcanically-induced tsunami since the catastrophic explosive eruption of Krakatau in 1883 and the sector collapse of Ritter Island in 1888. For the first time in over 100 years, the 2018 Anak Krakatau event provides an opportunity to study a major volcanically-generated tsunami that caused widespread loss of life and significant damage. Here, we present numerical simulations of the tsunami, with state-of the-art numerical models, based on a combined landslide-source and bathymetric dataset. We constrain the geometry and magnitude of the landslide source through analyses of pre- and post-event satellite images and aerial photography, which demonstrate that the primary landslide scar bisected the Anak Krakatau volcano, cutting behind the central vent and removing 50% of its subaerial extent. Estimated submarine collapse geometries result in a primary landslide volume range of 0.22-0.30 km
, which is used to initialize a tsunami generation and propagation model with two different landslide rheologies (granular and fluid). Observations of a single tsunami, with no subsequent waves, are consistent with our interpretation of landslide failure in a rapid, single phase of movement rather than a more piecemeal process, generating a tsunami which reached nearby coastlines within ~30 minutes. Both modelled rheologies successfully reproduce observed tsunami characteristics from post-event field survey results, tide gauge records, and eyewitness reports, suggesting our estimated landslide volume range is appropriate. This event highlights the significant hazard posed by relatively small-scale lateral volcanic collapses, which can occur en-masse, without any precursory signals, and are an efficient and unpredictable tsunami source. Our successful simulations demonstrate that current numerical models can accurately forecast tsunami hazards from these events. In cases such as Anak Krakatau's, the absence of precursory warning signals together with the short travel time following tsunami initiation present a major challenge for mitigating tsunami coastal impact.
Several lines of evidence have previously been used to suggest that ice retreat after the last glacial maximum (LGM) resulted in regionally-increased levels of volcanic activity. It has been proposed ...that this increase in volcanism was globally significant, forming a substantial component of the post-glacial rise in atmospheric CO2, and thereby contributing to climatic warming. However, as yet there has been no detailed investigation of activity in glaciated volcanic arcs following the LGM. Arc volcanism accounts for 90% of present-day subaerial volcanic eruptions. It is therefore important to constrain the impact of deglaciation on arc volcanoes, to understand fully the nature and magnitude of global-scale relationships between volcanism and glaciation.
The first part of this paper examines the post-glacial explosive eruption history of the Andean southern volcanic zone (SVZ), a typical arc system, with additional data from the Kamchatka and Cascade arcs. In all cases, eruption rates in the early post-glacial period do not exceed those at later times at a statistically significant level. In part, the recognition and quantification of what may be small (i.e. less than a factor of two) increases in eruption rate is hindered by the size of our datasets. These datasets are limited to eruptions larger than 0.1km3, because deviations from power-law magnitude–frequency relationships indicate strong relative under-sampling at smaller eruption volumes. In the southern SVZ, where ice unloading was greatest, eruption frequency in the early post-glacial period is approximately twice that of the mid post-glacial period (although frequency increases again in the late post-glacial). A comparable pattern occurs in Kamchatka, but is not observed in the Cascade arc. The early post-glacial period also coincides with a small number of very large explosive eruptions from the most active volcanoes in the southern and central SVZ, consistent with enhanced ponding of magma during glaciation and release upon deglaciation.
In comparison to non-arc settings, evidence of post-glacial increases in rates of arc volcanism is weak, and there is no need to invoke significantly increased melt production upon ice unloading, as occurred in areas such as Iceland. Non-arc volcanoes may therefore account for a relatively higher proportion of global volcanic emissions in the early post-glacial period than is suggested by the relative contributions of arc and non-arc settings at the present day.
The second part of this paper critically examines global eruption records, in an effort to constrain global-scale changes in volcanic output since the LGM. Accurate interpretation of these records relies on correcting both temporal and spatial variability in eruption recording. In particular, very low recording rates, which also vary spatially by over two orders of magnitude, prevent precise, and possibly even accurate, quantitative analysis. For example, if we assume record completeness for the past century, the number of known eruptions (volcanic explosivity index≥2) from some low-latitude regions, such as Indonesia, is approximately 1 in 20,000 (0.005%) for the period 5–20ka. There is a need for more regional-scale studies of past volcanism in such regions, where current data are extremely sparse. We attempt to correct for recording biases, and suggest a maximum two-fold (but potentially much less) increase in global eruption rates, relative to the present day, between 13 and 7ka. Although volcanism may have been an important source of CO2 in the early Holocene, it is unlikely to have been a dominant control on changes in atmospheric CO2 after the LGM.
Submarine landslides on open continental slopes can be prodigious in scale. They are an important process for global sediment fluxes, and can generate very damaging tsunamis. Submarine landslides are ...far harder to monitor directly than terrestrial landslides, and much greater uncertainty surrounds their preconditioning factors and triggers. Submarine slope failure often occurs on remarkably low (< 2°) gradients that are almost always stable on land, indicating that particularly high excess pore pressures must be involved. Earthquakes trigger some large submarine landslides, but not all major earthquakes cause widespread slope failure. The headwalls of many large submarine landslides appear to be located in water depths that are too deep for triggering by gas hydrate dissociation. The available evidence indicates that landslide occurrence is either weakly (or not) linked to changes in sea level or atmospheric methane abundance, or the available dates for open continental slope landslides are too imprecise to tell. Similarly, available evidence does not strongly support a view that landslides play an important role in methane emissions that cause climatic change. However, the largest and best-dated open continental slope landslide (the Storegga Slide) coincides with a major cooling event 8,200 years ago. This association suggests that caution may be needed when stating that there is no link between large open slope landslides and climate change.
The 1888 Ritter Island volcanic sector collapse triggered a regionally damaging tsunami. Historic eyewitness accounts allow the reconstruction of the arrival time, phase and height of the tsunami ...wave at multiple locations around the coast of New Guinea and New Britain. 3D seismic interpretations and sedimentological analyses indicate that the catastrophic collapse of Ritter Island was preceded by a phase of deep-seated gradual spreading within the volcanic edifice and accompanied by a submarine explosive eruption, as the volcanic conduit was cut beneath sea level. However, the potential impact of the deep-seated deformation and the explosive eruption on tsunami genesis is unclear. For the first time, it is possible to parameterise the different components of the Ritter Island collapse with 3D seismic data, and thereby test their relative contributions to the tsunami. The modelled tsunami arrival times and heights are in good agreement with the historic eyewitness accounts. Our simulations reveal that the tsunami was primarily controlled by the displacement of the water column by the collapsing cone at the subaerial-submarine boundary and that the submerged fraction of the slide mass and its mobility had only a minor effect on tsunami genesis. This indicates that the total slide volume, when incorporating the deep-seated deforming mass, is not directly scalable for the resulting tsunami height. Furthermore, the simulations show that the tsunamigenic impact of the explosive eruption energy during the Ritter Island collapse was only minor. However, this relationship may be different for other volcanogenic tsunami events with smaller slide volumes or larger magnitude eruptions, and should not be neglected in tsunami simulations and hazard assessment.
On May 1st 2008 Mount Chaitén (southern Chile) interrupted a long period of quiescence, generating a sequence of explosive eruptions and causing the evacuation of Chaitén town located a few ...kilometers south of the volcano. The activity was characterized by several explosive events each associated with plumes which reached up to about 19 km above sea level. The products were dispersed across a wide area, with the finest ash reaching the Atlantic coast of Argentina. Our field observations in the proximal-medial area (3–25 km from the vent) indicate that the May 2008 tephra deposit consists of numerous layers, most of which can be correlated with individual eruptive events. These layers vary from extremely fine-grained ash to layers of lapilli and blocks, composed of both juvenile and lithic material. Here we describe the stratigraphy and physical characteristics of the May 2008 deposits, and propose a reconstruction of the timing of the May 2008 events. The deposits are mainly associated with the three main explosive phases which occurred on 1st–2nd May, 3rd–5th May and 6th May, with an estimated bulk tephra volume of 0.5–1.0 km
3
(integration of both exponential and power-law fitting). For the 6th May event, represented by a layer composed mainly of lithic lapilli and blocks (>2 mm), an isopleth map was compiled from which a 19 km plume height was determined, which is in good agreement with satellite observations.
Detailed knowledge of the past history of an active volcano is crucial for the prediction of the timing, frequency and style of future eruptions, and for the identification of potentially at-risk ...areas. Subaerial volcanic stratigraphies are often incomplete, due to a lack of exposure, or burial and erosion from subsequent eruptions. However, many volcanic eruptions produce widely-dispersed explosive products that are frequently deposited as tephra layers in the sea. Cores of marine sediment therefore have the potential to provide more complete volcanic stratigraphies, at least for explosive eruptions. Nevertheless, problems such as bioturbation and dispersal by currents affect the preservation and subsequent detection of marine tephra deposits. Consequently, cryptotephras, in which tephra grains are not sufficiently concentrated to form layers that are visible to the naked eye, may be the only record of many explosive eruptions. Additionally, thin, reworked deposits of volcanic clasts transported by floods and landslides, or during pyroclastic density currents may be incorrectly interpreted as tephra fallout layers, leading to the construction of inaccurate records of volcanism. This work uses samples from the volcanic island of Montserrat as a case study to test different techniques for generating volcanic eruption records from marine sediment cores, with a particular relevance to cores sampled in relatively proximal settings (i.e. tens of kilometres from the volcanic source) where volcaniclastic material may form a pervasive component of the sedimentary sequence. Visible volcaniclastic deposits identified by sedimentological logging were used to test the effectiveness of potential alternative volcaniclastic-deposit detection techniques, including point counting of grain types (component analysis), glass or mineral chemistry, colour spectrophotometry, grain size measurements, XRF core scanning, magnetic susceptibility and X-radiography. This study demonstrates that a set of time-efficient, non-destructive and high-spatial-resolution analyses (e.g. XRF core-scanning and magnetic susceptibility) can be used effectively to detect potential cryptotephra horizons in marine sediment cores. Once these horizons have been sampled, microscope image analysis of volcaniclastic grains can be used successfully to discriminate between tephra fallout deposits and other volcaniclastic deposits, by using specific criteria related to clast morphology and sorting. Standard practice should be employed when analysing marine sediment cores to accurately identify both visible tephra and cryptotephra deposits, and to distinguish fallout deposits from other volcaniclastic deposits.
Global arc magmatism is sustained by a continuous fluid flux that is returned to the mantle in subduction zones. Despite considerable advances in simulations of melting processes, models of arc ...magmatism remain incompletely tested against erupted products. Here, we show that a suite of primitive volcanic rocks from across the southern Chilean arc preserves the signature of a systematic down‐slab gradient in fluid chemistry. The chemical gradient is consistent with predictions from modeling, geothermometry and experiments. We infer that increasing slab‐surface temperatures cause the sub‐arc slab flux to become less water‐rich and increasingly dominated by hydrous melts over a distance of a few kilometers behind the arc front. This change exerts a first‐order control on magma chemistry, and implies discrete melt‐transport pathways through subduction zones. Our results replicate patterns in other arcs, implying common sub‐arc slab‐surface temperature ranges in thermally‐diverse subduction zones.
Key Points
Mafic arc rocks preserve a systematic down‐slab magma chemical gradient
This reflects increasing dominance of hydrous melts in the slab flux
This change, over a few km, is a first‐order control on melt chemistry
Spatial variations in the thickness and grain-size characteristics of tephra fall deposits imply that tephra depositional processes cannot be fully captured by models of single-particle sedimentation ...from the base of the eruption plume. Here, we document a secondary thickness maximum in a ∼9.75 ka tephra fall deposit from Chaitén volcano, Chile (Cha1 eruption). This secondary thickness maximum is notably coarser-grained than documented historical examples, being dominated by medium-grained ash, and an origin via particle aggregation is therefore unlikely. In the region of secondary thickening, we propose that high levels of atmospheric turbulence accelerated particles held within the mid- to lower-troposphere (0 to ∼6 km) towards the ground surface. We suggest that this enhancement in vertical atmospheric mixing was driven by the breaking of lee waves, generated by winds passing over elevated topography beneath the eruption plume. Lower atmospheric circulation patterns may exert a significant control on the dispersal and deposition of tephra from eruption plumes across all spatial scales, particularly in areas of complex topography.
The major explosive eruption of Chaitén volcano, Chile, in May 2008 provided a rare opportunity to track the long‐range dispersal and deposition of fine volcanic ash. The eruption followed ∼10,000 ...years of quiescence, was the largest explosive eruption globally since Hudson, Chile, in 1991, and was the first explosive rhyolitic eruption since Novarupta, Alaska, in 1912. Field examination of distal ashfall indicates that ∼1.6 × 1011 kg of ash (dense rock equivalent volume of ∼0.07 km3) was deposited over ∼2 × 105 km2 of Argentina during the first week of eruption. The minimum eruption magnitude, estimated from the mass of the tephra deposit, is 4.2. Several discrete ashfall units are identifiable from their distribution and grain size characteristics, with more energetic phases showing a bimodal size distribution and evidence of cloud aggregation processes. Ash chemistry was uniform throughout the early stages of eruption and is consistent with magma storage prior to eruption at depths of 3–6 km. Deposition of ash over a continental region allowed the tracking of eruption development and demonstrates the potential complexity of tephra dispersal from a single eruption, which in this case comprised several phases over a week‐long period of intense activity.