Search for: All records

Stromboli volcano (Italy), always active with low energy explosive activity, is a very attractive place for visitors, scientists, and inhabitants of the island. Nevertheless, occasional more intense eruptions can present a serious danger. This study focuses on the modeling and estimation of their inter-event time and temporal rate. With this aim we constructed a new historical catalog of major explosions and paroxysms through a detailed review of scientific literature of the last ca. 140 years. The catalog includes the calendar date and phenomena descriptions for 180 explosive events, of which 36 were paroxysms. We evaluated the impact of the main sources of uncertainty affecting the historical catalog. In particular, we categorized as uncertain 45 major explosions that reportedly occurred before 1985 and tested the effect of excluding these events from our analysis. Moreover, after analyzing the entire record in the period [1879, 2020], we separately considered, as sequences, events in [1879, 1960] and in [1985, 2020] because of possible under recording issues in the period [1960, 1985]. Our new models quantify the temporal rate of major explosions and paroxysms as a function of time passed since the last event occurred. Recurrence hazard levels are found to be significantly elevated in the weeks and months following a major explosion or paroxysm, and then gradually decrease over longer periods. Computed hazard functions are also used to illustrate a methodology for estimating order-of-magnitude individual risk of fatality under certain basis conditions. This study represents a first quantitatively formal advance in determining long-term hazard levels at Stromboli.

 

Effective volcanic hazard management in regions where populations live in close proximity to persistent volcanic activity involves understanding the dynamic nature of hazards, and associated risk. Emphasis until now has been placed on identification and forecasting of the escalation phase of activity, in order to provide adequate warning of what might be to come. However, understanding eruption hiatus and post-eruption unrest hazards, or how to quantify residual hazard after the end of an eruption, is also important and often key to timely post-eruption recovery. Unfortunately, in many cases when the level of activity lessens, the hazards, although reduced, do not necessarily cease altogether. This is due to both the imprecise nature of determination of the “end” of an eruptive phase as well as to the possibility that post-eruption hazardous processes may continue to occur. An example of the latter is continued dome collapse hazard from lava domes which have ceased to grow, or sector collapse of parts of volcanic edifices, including lava dome complexes. We present a new probabilistic model for forecasting pyroclastic density currents (PDCs) from lava dome collapse that takes into account the heavy-tailed distribution of the lengths of eruptive phases, the periods of quiescence, and the forecast window of interest. In the hazard analysis, we also consider probabilistic scenario models describing the flow’s volume and initial direction. Further, with the use of statistical emulators, we combine these models with physics-based simulations of PDCs at Soufrière Hills Volcano to produce a series of probabilistic hazard maps for flow inundation over 5, 10, and 20 year periods. The development and application of this assessment approach is the first of its kind for the quantification of periods of diminished volcanic activity. As such, it offers evidence-based guidance for dome collapse hazards that can be used to inform decision-making around provisions of access and reoccupation in areas around volcanoes that are becoming less active over time. 

Statistical emulators are a key tool for rapidly producing probabilistic hazard analysis of geophysical processes. Given output data computed for a relatively small number of parameter inputs, an emulator interpolates the data, providing the expected value of the output at untried inputs and an estimate of error at that point. In this work, we propose to fit Gaussian Process emulators to the output from a volcanic ash transport model, Ash3d. Our goal is to predict the simulated volcanic ash thickness from Ash3d at a location of interest using the emulator. Our approach is motivated by two challenges to fitting emulators—characterizing the input wind field and interactions between that wind field and variable grain sizes. We resolve these challenges by using physical knowledge on tephra dispersal. We propose new physically motivated variables as inputs and use normalized output as the response for fitting the emulator. Subsetting based on the initial conditions is also critical in our emulator construction. Simulation studies characterize the accuracy and efficiency of our emulator construction and also reveal its current limitations. Our work represents the first emulator construction for volcanic ash transport models with considerations of the simulated physical process. 

Abstract This study presents a new method, called the Radial Interpolation Method, to interpolate data characterized by an approximately radial pattern around a relatively constrained central zone, such as the ground deformation patterns shown in many active volcanic areas. The method enables the fast production of short-term deformation maps on the base of spatially sparse ground deformation measurements and can provide uncertainty quantification on the interpolated values, fundamental for hazard assessment purposes and deformation source reconstruction. The presented approach is not dependent on a priori assumptions about the geometry, location and physical properties of the source, except for the requirement of a locally radial pattern, i.e., allowing multiple centers of symmetry. We test the new method on a synthetic point source example, and then, we apply the method to selected time intervals of real geodetic data collected at the Campi Flegrei caldera during the last 39 years, including examples of leveling, Geodetic Precise Traversing measurements and Global Positioning System. The maps of horizontal displacement, calculated inland, show maximum values lying along a semicircular annular region with a radius of about 2–3 km in size. This semi-annular area is marked by mesoscale structures such as faults, sand dikes and fractures. The maps of vertical displacement describe a linear relation between the maximum vertical uplift measured and the volume variation. The multiplicative factor in the linear relation is about 0.3 × 10 6  m 3 /cm if we estimate the proportion of the Δ V that is captured by the GPS network onland and we use this to estimate the full Δ V . In this case, the 95% confidence interval on K because of linear regression is ± 5%. Finally, we briefly discuss how the new method could be used for the production of short-term vent opening maps on the base of real-time geodetic measurements of the horizontal and vertical displacements. 

The behaviour of low-viscosity, pressure-driven compressible pore fluid flows in viscously deformable porous media is studied here with specific application to gas flow in lava domes. The combined flow of gas and lava is shown to be governed by a two-equation set of nonlinear mixed hyperbolic–parabolic type partial differential equations describing the evolution of gas pore pressure and lava porosity. Steady state solution of this system is achieved when the gas pore pressure is magmastatic and the porosity profile accommodates the magmastatic pressure condition by increased compaction of the medium with depth. A one-dimensional (vertical) numerical linear stability analysis (LSA) is presented here. As a consequence of the pore-fluid compressibility and the presence of gravitation compaction, the gradients present in the steady-state solution cause variable coefficients in the linearized equations which generate instability in the LSA despite the diffusion-like and dissipative terms in the original system. The onset of this instability is shown to be strongly controlled by the thickness of the flow and the maximum porosity, itself a function of the mass flow rate of gas. Numerical solutions of the fully nonlinear system are also presented and exhibit nonlinear wave propagation features such as shock formation. As applied to gas flow within lava domes, the details of this dynamics help explain observations of cyclic lava dome extrusion and explosion episodes. Because the instability is stronger in thicker flows, the continued extrusion and thickening of a lava dome constitutes an increasing likelihood of instability onset, pressure wave growth and ultimately explosion.