A model for the numerical simulation of tephra fall deposits

A simple semianalytical model to simulate ash dispersion and deposition produced by sustained Plinian and sub-Plinian eruption columns based on the 2D advection–dispersion equation was applied. The eruption column acts as a vertical line source with a given mass distribution and neglects the complex dynamics within the eruption column. Thus, the use of the model is limited to areas far from the vent where the dynamics of the eruption column play a minor role. Vertical wind and diffusion components are considered negligible with respect to the horizontal ones. The dispersion and deposition of particles in the model is only governed by gravitational settling, horizontal eddy diffusion, and wind advection. The model accounts for different types and size classes of a user-defined number of particle classes and changing settling velocity with altitude. In as much as wind profiles are considered constant on the entire domain, the model validity is limited to medium-range distances (about 30–200 km away from the source).The model was used to reconstruct the tephra fall deposit from the documented Plinian eruption of Mt. Vesuvius, Italy, in 79 A.D. In this case, the model was able to broadly reproduce the characteristic medium-range tephra deposit. The results support the validity of the model, which has the advantage of being simple and fast to compute. It has the potential to serve as a simple tool for predicting the distribution of ash fall of hypothetical or real eruptions of a given magnitude and a given wind profile. Using a statistical set of frequent wind profiles, it also was used to construct air fall hazard maps of the most likely affected areas around active volcanoes where a large eruption is expected to occur.     

Complex changes in eruption dynamics during the 79 AD eruption of Vesuvius

The 79 AD eruption of Vesuvius included 8 eruption units (EU1–8) and several complex transitions in eruptive style. This study focuses on two important transitions: (1) the abrupt change from white to gray pumice during the Plinian phase of the eruption (EU2 to EU3) and (2) the shift from sustained Plinian activity to the onset of caldera collapse (EU3 to EU4). Quantification of the textural features within individual pumice clasts reveals important changes in both the vesicles and groundmass crystals across each transition boundary. Clasts from the white Plinian fall deposit (EU2) present a simple story of decompression-driven crystallization followed by continuous bubble nucleation, growth and coalescence in the eruptive conduit. In contrast, pumices from the overlying gray Plinian fall deposit (EU3) are heterogeneous and show a wide range in both bubble and crystal textures. Extensive bubble growth, coalescence, and the onset of bubble collapse in pumices at the base of EU3 suggest that the early EU3 magma experienced protracted vesiculation that began during eruption of the EU2 phase and was modified by the physical effects of syn-eruptive mingling-mixing. Pumice clasts from higher in EU3 show higher bubble and crystal number densities and less evidence of bubble collapse, textural features that are interpreted to reflect more thorough mixing of two magmas by this stage of the eruption, with consequent increases in both vesiculation and crystallization. Pumice clasts from a short-lived, high column at the onset of caldera collapse (EU4) continue the trend of increasing crystallization (enhanced by mixing) but, unexpectedly, the melt in these clasts is more vesicular than in EU3 and, in the extreme, can be classified as reticulite. We suggest that the high melt vesicularity of EU4 reflects strong decompression following the partial collapse of the magma chamber.Editorial responsibility: D.B. Dingwell     

Digital topography of volcanoes from radar interferometry: an example from Mt Vesuvius,Italy

A new airborne radar technique can generate digital topographic data for volcanoes at a scale of 10 m spatial and 1–5 m vertical, with a swath width of 6.4 km. Called TOPSAR, the intrument is an interferometric radar flown on the NASA DC-8 aircraft. TOPSAR data permit the quantification of volcano slopes, volumes, and heights, and as such will be valuable for the analysis of lava flows, domes, and lahar channels. This instrument will be flown over several volcanoes in the near future, providing volcanologists with valuable data sets for the analysis of high-resolution topography. We briefly illustrate the potential use of TOPSAR data through examples from Mt Somma and Vesuvius, Italy.     

Shallow structure of the Somma–Vesuvius volcano from 3D inversion of gravity data

A gravity investigation was carried out in the Somma–Vesuvius complex area (Campania, Italy) based on a dataset recently enlarged with new measurements. These cover the volcanic top and fill some other important spatial gaps in previous surveys. Besides the new gravity map of the Vesuvius, we also present the results of a 3D inverse modelling, carried out by using constraints from deep well exploration and seismic reflection surveys. The resulting density model provides a complete reconstruction of the top of the carbonate basement. This is relevant mostly on the western side of the survey area, where no significant information was previously available. Other new information regards the Somma–Vesuvius structure. It consists of an annular volume of rocks around the volcanic vent and that extends down to the carbonate basement. It results to be denser with respect to the surrounding sedimentary cover of the Campanian Plain and to the material located just along the central axis of the volcanic structure. The coherence between these features and other geophysical evidences from previous studies, will be discussed together with the other results of this research.     

Volcaniclastic debris-flow occurrences in the Campania region (Southern Italy) and their relation to Holocene–Late Pleistocene pyroclastic fall deposits: implications for large-scale hazard mapping

The Campania Region (southern Italy) is characterized by the frequent occurrence of volcaniclastic debris flows that damage property and loss of life (more than 170 deaths between 1996 and 1999). Historical investigation allowed the identification of more than 500 events during the last four centuries; in particular, more than half of these occurred in the last 100 years, causing hundreds of deaths. The aim of this paper is to quantify debris-flow hazard potential in the Campania Region. To this end, we compared several elements such as the thickness distribution of pyroclastic fall deposits from the last 18 ka of the Vesuvius and Phlegrean Fields volcanoes, the slopes of relieves, and the historical record of volcaniclastic debris flows from A.D. 1500 to the present. Results show that flow occurrence is not only a function of the cumulative thickness of past pyroclastic fall deposits but also depends on the age of emplacement. Deposits younger than 10 ka (Holocene eruptions) apparently increase the risk of debris flows, while those older than 10 ka (Late Pleistocene eruptions) seem to play a less prominent role, which is probably due to different climatic conditions, and therefore different rates of erosion of pyroclastic falls between the Holocene and the Late Pleistocene. Based on the above considerations, we compiled a large-scale debris-flow hazard map of the study area in which five main hazard zones are identified: very low, low, moderate, high, and very high.     

Luminescence dating of scoria fall and lahar deposits from Somma–Vesuvius,Italy

Luminescence dating has been applied to scoria and lahar deposits from Somma–Vesuvius, Italy. Samples include scoria from the AD472 and 512 (or 536) eruptions and lahar deposits. In order to find a stable luminescence signal which is less affected by anomalous fading, infrared stimulated luminescence (IRSL) signals at elevated temperatures after bleaching with IR at 50 °C (termed post-IR IRSL; pIRIR) were tested at different preheat and elevated stimulation temperatures. The fading rates of both IRSL and pIRIR signals reduced dramatically with increasing preheat and pIRIR stimulation temperatures. A pIRIR signal measured at 290 °C after a preheat at 320 °C (60 s) and an IR stimulation at 50 °C (100 s) was selected to calculate the equivalent dose (D_e). The gamma spectrometry results indicate that the U-series nuclides are not in equilibrium and there is a large ²²⁶Ra excess. The dose rates and ages were calculated by assuming a ²²⁶Ra excess (over its parent ²³⁰Th) at deposition, and that this unsupported excess then decayed to the present level. The resulting luminescence ages of the two scoria samples agreed with the expected ages, and the ages of the lahar deposits indicate that they are associated with the AD1631 eruption.     

Behaviour of chlorine prior and during the 79 A.D. Plinian eruption of Vesuvius (southern Italy) as inferred from the present distribution in glassy mesostases and whole-pumices

This study deals with the distribution of chlorine in glassy mesostases and whole-pumices from the 79 A.D. Plinian eruption (Somma–Vesuvius volcanic complex, Italy). This explosive event produced a prominent Plinian fall deposit followed by flows and surges. The fall deposit can be divided into two sub-units on the basis of an abrupt change in colour at approximately mid-level: a phonolitic white pumice layer at the base and a tephriphonolitic grey pumice layer at the top. Due to its hybrid nature (a mixture of k-tephritic, tephriphonolitic and phonolitic magmas), information on chlorine behaviour in tephriphonolitic magma has been inferred only by means of mass balance calculations. In the white pumices, chlorine concentrations show constant whole-pumice values, whereas glassy mesostases display significant compositional variations. These variations have been linked to the cryptocrystallisation of leucite in glassy mesostases, which affected the original melt compositions just before and during the eruptive event. In this framework, whole-pumice appears to better represent the pre-eruptive melt compositions. Using chlorine concentration in whole-pumices, a three-stage model of chlorine behaviour prior and during the eruptive event is predicted: (i) free variation during Rayleigh fractionation, according to a system with variance greater than zero; (ii) exsolution of volatile chlorine compounds (e.g., metal chlorides), when chlorine reaches its solubility limit in silicate melt, in coexisting hyper-saline and in dilute immiscible fluids; the variance of the system is zero at a given temperature and pressure; (iii) residual syn-eruptive variable enrichment of chlorine in the melt due to cryptocrystallization of leucite, suggesting a very minor loss of chlorine in the gas phase during the eruption. Although chlorine does not behave as a volatile element during the eruption, it is present in the volcanic plume. This is due to the postulated ‘excess' fluid phase containing chlorine that formed in the magma chamber prior to the eruption. The homogeneous distribution of chlorine in whole pumices, in contrast with a well-established chemical and isotopic layering in Vesuvian magmas prior to Pompei eruption, suggests that the trace element zonation is not directly linked to chlorine distribution in silicic melts.     

Physical Modelling and Human Survival in Pyroclastic Flows

Volcanic eruptions increasingly present catastrophic natural risks with hundreds of millions of people now living in areas of active volcanism and major conurbations around active eruptive centres. Interdisciplinary studies in disaster reduction have an important role in volcanic emergency management through advancing our understanding of the physical impacts of eruptive phenomena and the causes of death and injury in explosive eruptions. Numerical modelling of pyroclastic flows, amongst the most destructive of eruptive phenomena, provides new opportunities to improve the evaluation of the potential destructiveness of volcanic events and their human impacts in densely populated areas. In this work, the results of numerical modelling of pyroclastic flow propagation at Vesuvius have been analysed in terms of the physical parameters (temperature, ash in air concentration, and dynamic pressure) that are most critical for human survival. Our numerical simulations of eruptions of Vesuvius indicate that a large area exists where total destruction may not be inevitable in small to medium scale events, a finding that has prompted us to explore further the implications for human survival as part of an interdisciplinary approach to disaster reduction. The lessons of modelling at Vesuvius should be integrated into civil protection plans for other urban centres threatened by volcanoes.     

Three Decades of Seismic Activity at Mt. Vesuvius: 1972–2000

We analyse the seismic catalogue of the local earthquakes which occurred at Somma-Vesuvius volcano in the past three decades (1972–2000). The seismicity in this period can be described as composed of a background level, characterised by a low and rather uniform rate of energy release and by sporadic periods of increased seismic activity. Such relatively intense seismicity periods are characterised by energy rates and magnitudes progressively increasing in the critical periods. The analysis of the b value in the whole period evidences a well-defined pattern, with values of b progressively decreasing, from about 1.8 at the beginning of the considered period, to about 1.0 at present. This steady variation indicates an increasing dynamics in the volcanic system. Within this general trend it is possible to identify a substructure in the time sequence of the seismic events, formed by the alternating episodes of quiescence and activity. The analysis of the source moment tensor of the largest earthquakes shows that the processes at the seismic source are generally not consistent with simple double-couples, but that they are compatible with isotropic components, mostly indicating volumetric expansion. These components are shown to be statistically significant for most of the analysed events. Such focal mechanisms can be interpreted as the effect of explosion phenomena, possibly related to volatile exsolution from the crystallising magma. The availability of a reduced amount of high quality data necessary for the inversion of the source moment tensor, the still limited period of systematic observation of Vesuvius micro-earthquakes and, above all, the absence of eruptive events during such interval of time, cannot obviously permit the outlining of any formal premonitory signal. Nevertheless, the analysis reported in this paper indicates a progressively evolving dynamics, characterised by a generally increasing trend in the seismic activity in the volcanic system and by a significant volumetric component of recent major events, thus posing serious concern for a future evolution towards eruptive activity.