Magmatic History of Somma-Vesuvius on the Basis of New Geochemical and Isotopic Data from a Deep Borehole (Camaldoli della Torre)

A continuous-coring borehole recently drilled at Camaldoli dellaTorre on the southern slopes of Somma–Vesuvius providesconstraints on the volcanic and magmatic history of the Vesuvianvolcanic area since c. 126 ka BP. The cored sequence includesvolcanic units, defined on stratigraphical, sedimentological,petrological and geochemical grounds, emitted from both localand distal vents. Some of these units are of known age, suchas one Phlegraean pre-Campanian Ignimbrite, Campanian Ignimbrite(39 ka), Neapolitan Yellow Tuff (14· 9 ka) and VesuvianPlinian deposits, which helps to constrain the relative ageof the other units. The main rock types encountered are shoshonite,phonotephrite, latite, trachyte and phonolite. The sequenceincludes, from the base upwards: a thick succession of pyroclasticunits emplaced between 126 and 39 ka, most of them attributedto eruptions that occurred in the Phlegraean area; the CampanianIgnimbrite; the products of a local tuff cone formed between39 ka and the deposition of the products of the earliest activityof the Mt. Somma volcano; the products of the Somma–Vesuviusvolcano, which include from the base upwards a thick sequenceof lavas, pyroclastic rocks and the products of a local spattercone dated between 3· 7 ka and AD 79. The data obtainedfrom the study of the borehole show that, before the CampanianIgnimbrite eruption, low-energy explosive volcanism took placein the Vesuvian area, whereas mostly high-energy explosive eruptionscharacterized the Campi Flegrei activity. In the Vesuvian area,Campanian Ignimbrite deposition was followed by the eruptionof a local tuff cone and a long repose time, which predatedthe formation of the Mt. Somma edifice. Since 18· 3 ka(Pomici di Base eruption) the activity of Somma–Vesuviusbecame mostly explosive with rare lava effusions. The shallowestcored deposits belong to the Camaldoli della Torre cone, formedbetween the Pomici di Avellino and Pomici di Pompei eruptions(3· 7 ka–AD 79). New geochemical and Sr–Nd–Pb–B-isotopicdata on samples from the drilled core, together with those availablefrom the literature, allow us to further distinguish the volcanicrocks as a function of both their provenance (i.e. Phlegraeanor Vesuvian areas) and age, and to identify different magmaticprocesses acting through time in the Vesuvian mantle source(s)and during magma ascent towards the surface. Isotopically distinctmagmas, rising from a mantle source variably contaminated byslab-derived components, stagnated at mid-crustal depths (8–10km below sea level) where magmas differentiated and were probablycontaminated. Contamination occurred either with Hercynian continentalcrust, mostly during the oldest stages of Vesuvian activity(from 39 to 16 ka), or with Mesozoic limestone, mostly duringrecent Vesuvian activity. Energy constrained assimilation andfractional crystallization (EC-AFC) modelling results show thatcontamination with Hercynian crust probably occurred duringdifferentiation from shoshonite to latite. Contamination withlimestone, which is not well constrained with the availabledata, might have occurred only during the transition from shoshoniteto tephrite. From the ‘deep’ reservoir, magmas rosetowards a series of shallow reservoirs, in which they differentiatedfurther, mixed, and fed volcanic activity. KEY WORDS: Somma–Vesuvius; crustal contamination; source heterogeneity; radiogenic and stable isotopes; magmatic system     

Vent conditions for expected eruptions at Vesuvius

Determining consistent sets of vent conditions for next expected eruptions at Vesuvius is crucial for the simulation of the sub-aerial processes originating the volcanic hazard and the eruption impact. Here we refer to the expected eruptive scales and conditions defined in the frame of the EC Exploris project, and simulate the dynamics of magma ascent along the volcanic conduit for sub-steady phases of next eruptions characterized by intensities of the Violent Strombolian (VS), Sub-Plinian 2 (SP2), and Sub-Plinian 1 (SP1) scale. Sets of conditions for the simulations are determined on the basis of the bulk of knowledge on the past history of Vesuvius [Cioni, R., Bertagnini, A., Santacroce, R., Andronico, D., Explosive activity and eruption scenarios at Somma–Vesuvius (Italy): towards a new classification scheme. Journal of Volcanology and Geothermal Research, this issue.]. Volatile contents (H₂O and CO₂) are parameterized in order to account for the uncertainty in their expected amounts for a next eruption. In all cases the flow in the conduit is found to be choked, with velocities at the conduit exit or vent corresponding to the sonic velocity in the two-phase non-equilibrium magmatic mixture. Conduit diameters and vent mixture densities are found to display minimum overlapping between the different eruptive scales, while exit gas and particle velocities, as well as vent pressures, largely overlap. Vent diameters vary from as low as about 5 m for VS eruptions, to 35–55 m for the most violent SP1 eruption scale. Vent pressures can be as low as less than 1 MPa for the lowest volatile content employed of 2 wt.% H₂O and no CO₂, to 7–8 MPa for highest volatile contents of 5 wt.% H₂O and 2 wt.% CO₂ and large eruptive scales. Gas and particle velocities at the vent range from 100–250 m/s, with a tendency to decrease, and to increase the mechanical decoupling between the phases, with increasing eruptive scale. Except for velocities, all relevant vent quantities are more sensitive to the volatile content of the discharged magma for the highest eruptive scales considered.     

3D Gravity Inversion by Growing Bodies and Shaping Layers at Mt. Vesuvius (Southern Italy)

To improve our knowledge of the structural pattern of Mt. Vesuvius and its magmatic system, which represents one of the three volcanoes located in the Neapolitan area (together with Campi Flegrei and Ischia; southern Italy), we analyze here the Bouguer gravity map that is already available through its interpretation by means of 2.5-dimensional modelling. We have carried out a three-dimensional interpretation using a new and original algorithm, known as ‘Layers’, that has been especially processed for this purpose. Layers works in an automatic and non-subjective way, and allows the definition of the structural settings in terms of several layers, each representing a specific geological formation. The same data are also interpreted in terms of isolated and shallow anomalous density bodies using a well tested algorithm known as ‘Growth’. We focus our inversions on the Mt. Vesuvius volcano, while globally analyzing the entire Neapolitan area, in order to investigate the deep structures, and in particular the deep extended ‘sill’ that has been revealed by seismic tomography. The final models generally confirm the global setting of the area as outlined by previous investigations, mainly for the shape and depth of the carbonate basement below Mt. Vesuvius. The presence of lateral density contrasts inside the volcano edifice is also shown, which was only hypothesized in the 2.5-dimensional inversion. Moreover, the models allow us to note a high density body that rises from the top of the carbonate basement and further elongates above sea level. This probably represents an uprising of the same basement, which is just below the volcano and which coincides with the V_P and V_P/V_S anomalies detected under the crater. The three-dimensional results also reveal that the two inversion methods provide very similar models, where the high density isolated body in the Growth model can be associated with the rising high density anomaly in the Layers model. Taking into account the density of these modelled bodies, we would also suggest that they represent solidified magma bodies, as suggested by other studies. Finally, we did not clearly detect any deep anomalous body that can be associated with the sill that was suggested by seismic tomography.     

On the fractal dimension of the fallout deposits: A case study of the 79 A.D. Plinian eruption at Mt. Vesuvius

We investigated the existence of a fractal law (power law) distribution of size pyroclastic fragments erupted during the fallout phase of the 79 A.D. Plinian eruption at Mt. Vesuvius. In particular, we performed a particle size distribution analysis on 18 white and grey pumice samples collected in six sites distributed in the SW sector of Mt. Vesuvius. Our measurements show that the fragmentation of samples in the investigated range (from 32 mm to 850 μm) follows a power law, guaranteeing the scale invariance of the process. The relationship frequency-size distribution of the fragments is verified independently from the nature (i.e., pumices and lithics) and stratigraphic height of the considered samples in the pyroclastic deposit. Therefore, the fractal fragmentation theory can be indicated for evaluating the relationship between the intensity of fragmentation (fractal dimension D) and eruption energy. In this way the apparent chaotic distribution of the particles in the fallout deposits hides a self-organized complexity revealed by the retrieved power law distribution. We further remark that a key aspect of our analysis is the founded evidence that the fractal dimension of the lithics is systematically greater than that of the pumices.     

Shallow and deep reservoirs involved in magma supply of the 1944 eruption of Vesuvius

 During the 1944 eruption of Vesuvius a sudden change occurred in the dynamics of the eruptive events, linked to variations in magma composition. K-phonotephritic magmas were erupted during the effusive phase and the first lava fountain, whereas the emission of strongly porphyritic K-tephrites took place during the more intense fountain. Melt inclusion compositions (major and volatile elements) highlight that the magmas feeding the eruption underwent differentiation at different pressures. The K-tephritic volatile-rich melts (up to 3 wt.% H₂O, 3000 ppm CO₂, and 0.55 wt.% Cl) evolved to reach K-phonotephritic compositions by crystallization of diopside and forsteritic olivine at total fluid pressure higher than 300 MPa. These magmas fed a very shallow reservoir. The low-pressure differentiation of the volatile-poor K-phonotephritic magmas (H₂O<1 wt.%) involved mixing, open-system degassing, and crystallization of leucite, salite, and plagioclase. The eruption was triggered by intrusion of a volatile-rich magma batch that rose from a depth of 11–22 km into the shallow magma chamber. The first phase of the eruption represents the partial emptying of the shallow reservoir, the top of which is within the volcanic edifice. The newly arrived magma mixed with that resident in the shallow reservoir and forced the transition from the effusive to the lava fountain phase of the eruption. Received: 14 September 1998 / Accepted: 10 January 1999     

Stable isotopes of pedogenic carbonates from the Somma–Vesuvius area,southern Italy,over the past 18 kyr: palaeoclimatic implications

Stable isotopes were measured in the carbonate and organic matter of palaeosols in the Somma–Vesuvius area, southern Italy in order to test whether they are suitable proxy records for climatic and ecological changes in this area during the past 18000 yr. The ages of the soils span from ca. 18 to ca. 3 kyr BP. Surprisingly, the Last Glacial to Holocene climate transition was not accompanied by significant change in δ¹⁸O of pedogenic carbonate. This could be explained by changes in evaporation rate and in isotope fractionation between water and precipitated carbonate with temperature, which counterbalanced the expected change in isotope composition of meteoric water. Because of the rise in temperature and humidity and the progressive increase in tree cover during the Holocene, the Holocene soil carbonates closely reflect the isotopic composition of meteoric water. A cooling of about 2°C after the Avellino eruption (3.8 ka) accounts for a sudden decrease of about 1‰ in δ¹⁸O of pedogenic carbonate recorded after this eruption. The δ¹³C values of organic matter and pedogenic carbonate covary, indicating an effective isotope equilibrium between the organic matter, as the source of CO₂, and the pedogenic carbonate. Carbon isotopes suggest prevailing C₃ vegetation and negligible mixing with volcanogenic or atmospheric CO₂. Copyright © 2000 John Wiley & Sons, Ltd.     

Analysis of the magnetic anomaly field of the volcanic district of the Bay of Naples, Italy

We here present and discuss the results of the analysis and qualitative interpretation of two magnetic surveys performed in the Bay of Naples in 1998 and 2000. A map of the Bay of Naples based on the data acquired during these surveys has already been published by the Italian CNR-IAMC Research Institute. We re-processed the same data to produce maps of the pole reduced, analytic signal and horizontal derivative data and correlated them with the bathymetry and the gravimetric data of the area. The analysis shows strong anomalies in the NW and NE volcanic areas of the Bay of Naples, while the central area seems magnetically quiet. In the Phlegrean area the maps clearly show the southern rim of the Phlegrean caldera and demonstrate that while the Magnaghi Canyon is correlated to gravimetric highs and magnetic structures, and can therefore be interpreted as an active lineament, most of Dohrn Canyon is not characterized by volcanic activity and does not correlate to any gravimetric or magnetic structures. An important round-shaped magnetic anomaly is for the first time identified in the central slope of the gulf between the two canyons, probably correlated to a large buried volcanic edifice. In the Vesuvian area some intense circular anomalies, aligned in the NNW–SSE direction, are localized in the Torre del Greco and Torre Annunziata offshore, related to the submerged part of Vesuvius and possibly connected to buried vents.     

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.