As demonstrated at Anak Krakatau on December 22
, 2018, tsunamis generated by volcanic flank collapse are incompletely understood and can be devastating. Here, we present the first high-resolution ...characterisation of both subaerial and submarine components of the collapse. Combined Synthetic Aperture Radar data and aerial photographs reveal an extensive subaerial failure that bounds pre-event deformation and volcanic products. To the southwest of the volcano, bathymetric and seismic reflection data reveal a blocky landslide deposit (0.214 ± 0.036 km
) emplaced over 1.5 km into the adjacent basin. Our findings are consistent with en-masse lateral collapse with a volume ≥0.175 km
, resolving several ambiguities in previous reconstructions. Post-collapse eruptions produced an additional ~0.3 km
of tephra, burying the scar and landslide deposit. The event provides a model for lateral collapse scenarios at other arc-volcanic islands showing that rapid island growth can lead to large-scale failure and that even faster rebuilding can obscure pre-existing collapse.
Understanding the location, intensity, and likely duration of volcanic hazards is key to reducing risk from volcanic eruptions. Here, we use a novel near-real-time dataset comprising Volcanic Ash ...Advisories (VAAs) issued over 10 years to investigate global rates and durations of explosive volcanic activity. The VAAs were collected from the nine Volcanic Ash Advisory Centres (VAACs) worldwide. Information extracted allowed analysis of the frequency and type of explosive behaviour, including analysis of key eruption source parameters (ESPs) such as volcanic cloud height and duration. The results reflect changes in the VAA reporting process, data sources, and volcanic activity through time. The data show an increase in the number of VAAs issued since 2015 that cannot be directly correlated to an increase in volcanic activity. Instead, many represent increased observations, including improved capability to detect low- to mid-level volcanic clouds (FL101–FL200, 3–6 km asl), by higher temporal, spatial, and spectral resolution satellite sensors. Comparison of ESP data extracted from the VAAs with the Mastin et al. (J Volcanol Geotherm Res 186:10–21,
2009a
) database shows that traditional assumptions used in the classification of volcanoes could be much simplified for operational use. The analysis highlights the VAA data as an exceptional resource documenting global volcanic activity on timescales that complement more widely used eruption datasets.
This study compares and evaluates one-dimensional (1D) and three-dimensional (3D) numerical models of volcanic eruption columns in a set of different inter-comparison exercises. The exercises were ...designed as a blind test in which a set of common input parameters was given for two reference eruptions, representing a strong and a weak eruption column under different meteorological conditions. Comparing the results of the different models allows us to evaluate their capabilities and target areas for future improvement. Despite their different formulations, the 1D and 3D models provide reasonably consistent predictions of some of the key global descriptors of the volcanic plumes. Variability in plume height, estimated from the standard deviation of model predictions, is within ~20% for the weak plume and ~10% for the strong plume. Predictions of neutral buoyancy level are also in reasonably good agreement among the different models, with a standard deviation ranging from 9 to 19% (the latter for the weak plume in a windy atmosphere). Overall, these discrepancies are in the range of observational uncertainty of column height. However, there are important differences amongst models in terms of local properties along the plume axis, particularly for the strong plume. Our analysis suggests that the simplified treatment of entrainment in 1D models is adequate to resolve the general behaviour of the weak plume. However, it is inadequate to capture complex features of the strong plume, such as large vortices, partial column collapse, or gravitational fountaining that strongly enhance entrainment in the lower atmosphere. We conclude that there is a need to more accurately quantify entrainment rates, improve the representation of plume radius, and incorporate the effects of column instability in future versions of 1D volcanic plume models.
•We present the main results of an eruptive column model inter-comparison exercise.•Simulations with standard inputs for strong and weak eruptive plumes were performed.•We compare results of empirical, one-dimensional, and three-dimensional models.•Results allowed for evaluating model capabilities and areas for model improvement.
The uncertainties associated with tephra thickness measurements are calculated and implications for volume estimates are presented. Statistical methods are used to analyse the large dataset of Walker ...and Croasdale J Geol Soc 127:17-55, 1971 of the Fogo A plinian deposit, São Miguel, Azores. Dirichlet tessellation demonstrates that Walker and Croasdale’s measurements are highly clustered spatially and the area represented by a single measurement ranges between 0.5 and 10 km
2
. K-means cluster analysis shows that lower thickness uncertainties are associated with closely spaced measurements. Re-examination and analysis of Fogo A fall deposits show thickness uncertainties are about 9% for measured thickness while uncertainty associated with natural variance ranges, between 10 and 40%, with an average error of 30%. Correlations between measurement uncertainties and natural variance are complex and depend on a unit’s thickness, position within a succession and distance from source. Normative error increases as tephra thickness decreases. The degree to which thickness measurement error impacts on volume uncertainty depends on the number of measurements within a given dataset and their associated uncertainty. The uncertainty in volume associated with thickness uncertainty calculated herein for Fogo A is 1.3%, equivalent to a volume of 0.02 km
3
. However uncertainties associated with smaller datasets can be much larger; for example typically exceeding 10% for less than 20 data points.
Tephra layers in marine sediment cores from scientific ocean drilling largely record high‐magnitude silicic explosive eruptions in the Japan arc for up to the last 20 million years. Analysis of the ...thickness variation with distance of 180 tephra layers from a global data set suggests that the majority of the visible tephra layers used in this study are the products of caldera‐forming eruptions with magnitude (M) > 6, considering their distances at the respective drilling sites to their likely volcanic sources. Frequency of visible tephra layers in cores indicates a marked increase in rates of large magnitude explosive eruptions at ∼8 Ma, 6–4 Ma, and further increase after ∼2 Ma. These changes are attributed to major changes in tectonic plate interactions. Lower rates of large magnitude explosive volcanism in the Miocene are related to a strike‐slip‐dominated boundary (and temporary cessation or deceleration of subduction) between the Philippine Sea Plate and southwest Japan, combined with the possibility that much of the arc in northern Japan was submerged beneath sea level partly due to previous tectonic extension of northern Honshu related to formation of the Sea of Japan. Changes in plate motions and subduction dynamics during the ∼8 Ma to present period led to (1) increased arc‐normal subduction in southwest Japan (and resumption of arc volcanism) and (2) shift from extension to compression of the upper plate in northeast Japan, leading to uplift, crustal thickening and favorable conditions for accumulation of the large volumes of silicic magma needed for explosive caldera‐forming eruptions.
Key Points
Increase in volcanic activity in Japan at ∼8, 4–6, and ∼2 Ma
Tephra layers in distal sediment cores represent magnitude 6 or greater eruptions
Tectonic changes driving rates of explosive volcanism in Japan
Characterization of explosive volcanic eruptive processes from interpretation of deposits is a key for assessing volcanic hazard and risk, particularly for infrequent large explosive eruptions and ...those whose deposits are transient in the geological record. While eruption size—determined by measurement and interpretation of tephra fall deposits—is of particular importance, uncertainties for such measurements and volume estimates are rarely presented. Here, tephra volume estimates are derived from isopach maps produced by modeling raw thickness data as cubic B-spline curves under tension. Isopachs are objectively determined in relation to original data and enable limitations in volume estimates from published maps to be investigated. The eruption volumes derived using spline isopachs differ from selected published estimates by 15–40 %, reflecting uncertainties in the volume estimation process. The formalized analysis enables identification of sources of uncertainty; eruptive volume uncertainties (>30 %) are much greater than thickness measurement uncertainties (~10 %). The number of measurements is a key factor in volume estimate uncertainty, regardless of method utilized for isopach production. Deposits processed using the cubic B-spline method are well described by 60 measurements distributed across each deposit; however, this figure is deposit and distribution dependent, increasing for geometrically complex deposits, such as those exhibiting bilobate dispersion.
In this study, we combine reconstructions of volcanological data sets and inputs from a structured expert judgment to produce a first long‐term probability map for vent opening location for the next ...Plinian or sub‐Plinian eruption of Somma‐Vesuvio. In the past, the volcano has exhibited significant spatial variability in vent location; this can exert a significant control on where hazards materialize (particularly of pyroclastic density currents). The new vent opening probability mapping has been performed through (i) development of spatial probability density maps with Gaussian kernel functions for different data sets and (ii) weighted linear combination of these spatial density maps. The epistemic uncertainties affecting these data sets were quantified explicitly with expert judgments and implemented following a doubly stochastic approach. Various elicitation pooling metrics and subgroupings of experts and target questions were tested to evaluate the robustness of outcomes. Our findings indicate that (a) Somma‐Vesuvio vent opening probabilities are distributed inside the whole caldera, with a peak corresponding to the area of the present crater, but with more than 50% probability that the next vent could open elsewhere within the caldera; (b) there is a mean probability of about 30% that the next vent will open west of the present edifice; (c) there is a mean probability of about 9.5% that the next medium‐large eruption will enlarge the present Somma‐Vesuvio caldera, and (d) there is a nonnegligible probability (mean value of 6–10%) that the next Plinian or sub‐Plinian eruption will have its initial vent opening outside the present Somma‐Vesuvio caldera.
Key Points
New vent opening probability maps for a future Plinian/sub‐Plinian eruption of Somma‐Vesuvio volcano are developed
The maps incorporate uncertainty estimates and uniquely include the special case of Somma‐Vesuvio caldera enlargement
New volcanological/structural data sets (from the companion paper) are combined to create spatial probability density functions
Coignimbrite plumes provide a common and effective mechanism by which large volumes of fine‐grained ash are injected into the atmosphere. Nevertheless, controls on formation of these plumes as a ...function of eruptive conditions are still poorly constrained. Herein, two 1‐D axysymmetric steady state models were coupled, the first describing the parent pyroclastic density current and the second describing plume rise. Global sensitivity analysis is applied to investigate controls on coignimbrite plume formation and describe coignimbrite source and the maximum plume height attained. For a range of initial mass flow rates between 108 and 1010 kg/s, modeled liftoff distance (the distance at which neutral buoyancy is attained), assuming radial supercritical flow, is controlled by the initial flow radius, gas mass fraction, flow thickness, and temperature. The predicted decrease in median grain size between flow initiation and plume liftoff is negligible. Calculated initial plume vertical velocities, assuming uniform liftoff velocity over the pyroclastic density current invasion area, are much greater (several tens of m/s) than those previously used in modeling coignimbrite plumes (1 m/s). Such velocities are inconsistent with the fine grain size of particles lofted into coignimbrite plumes, highlighting an unavailability of large clasts, possibly due to particle segregation within the flow, prior to plume formation. Source radius and initial vertical velocity have the largest effect on maximum plume height, closely followed by initial temperature. Modeled plume heights are between 25 and 47 km, comparable with Plinian eruption columns, highlighting the potential of such events for distributing fine‐grained ash over significant areas.
Key Points
Coignimbrite plumes act as an efficient mechanism of ash injection into the atmosphere
Sensitivity analysis enables investigation of relationship between model inputs and plume height
Initial flow mass flux controls both location of plume formation and maximum plume height
We present developments to the physical model and the open-source numerical code IMEX_SfloW2D (de' Michieli Vitturi et al., 2019). These developments consist of a generalization of the depth-averaged ...(shallow-water) fluid equations to describe a polydisperse fluid–solid mixture, including terms for sedimentation and entrainment, transport equations for solid particles of different sizes, transport equations for different components of the carrier phase, and an equation for temperature/energy. Of relevance for the simulation of volcanic mass flows, vaporization and entrainment of water are implemented in the new model. The model can be easily adapted to simulate a wide range of volcanic mass flows (pyroclastic avalanches, lahars, pyroclastic surges), and here we present its application to transient dilute pyroclastic density currents (PDCs). The numerical algorithm and the code have been improved to allow for simulation of sub- to supercritical regimes and to simplify the setting of initial and boundary conditions. The code is open-source. The results of synthetic numerical benchmarks demonstrate the robustness of the numerical code in simulating transcritical flows interacting with the topography. Moreover, they highlight the importance of simulating transient in comparison to steady-state flows and flows in 2D versus 1D. Finally, we demonstrate the model capabilities to simulate a complex natural case involving the propagation of PDCs over the sea surface and across topographic obstacles, through application to Krakatau volcano, showing the relevance, at a large scale, of non-linear fluid dynamic features, such as hydraulic jumps and von Kármán vortices, to flow conditions such as velocity and runout.
Tephra layers are frequently used to reconstruct past volcanic activity. Inferences made from tephra layers rely on the assumption that the preserved tephra layer is representative of the initial ...deposit. However, a great deal can happen to tephra after it is deposited; thus, tephra layer taphonomy is a crucial but poorly understood process. The overall goal of this research was to gain greater insight into the taphonomy of terrestrial tephra layers, specifically the extent to which deposit thickness is altered over time, with implications for tephra volume estimation. We approached this by a)conducting a new survey of the tephra layer from the recent, well-studied eruption of Mount St Helens on May 18th, 1980 (MSH1980); b)modelling the tephra layer thickness using a mathematical technique and c)comparing our results with an equivalent model based on measurements taken immediately after the eruption. In this way, we aimed to quantify any losses and transformations that have occurred. During our study, we collected measurements of tephra layer thickness from 86 locations ranging from <20 to >600 km from the vent. Geochemical analysis was used to identify tephra of uncertain origin. Our results indicated that the extant tephra layer at undisturbed sites was representative of the original deposit: overall, preservation in these locations (in terms of thickness, stratigraphy and geochemistry) had been remarkably good. However, the isopach maps generated from our measurements diverged from isopachs produced in the same way, but derived from the original survey data. Furthermore, our estimates of the quantity of tephra produced during the eruption greatly exceeded previous estimates of the fallout volume. In our study, inaccuracies in the modelled fallout arose from issues of sampling strategy, rather than taphonomy. Our results demonstrate the sensitivity of volcanological reconstructions to measurement location, and the great importance of reliably measured low/zero values in reconstructing tephra deposits.
•The tephra layer produced by the 1980 eruption of Mount St Helens was surveyed.•Isopach maps derived from a)our survey & b)measurements made in 1980 were compared.•Preservation of the tephra layer was good in undisturbed areas.•Despite good preservation, our model overestimated the volume of the tephra deposit.•Our results demonstrate the sensitivity of reconstructions to measurement location.