The structure of the Campanian Plain and the activity of the Neapolitan volcanoes (Italy)

The central Campanian Plain is dominated by the structural depression of Acerra whose origin is tectonic, but may have been enlarged and further depressed after the eruption of the Campanian Ignimbrite (42-25 ka). The deposits of the Campanian Ignimbrite are possibly the results of multiple eruptions with huge pyroclastic deposits that covered all the Campanian Plain.The more recent activity of Vesuvius, Campi Flegrei and Procida occurred on the borders of Acerra depression and resulted from a reactivation of regional faults after the Campanian Ignimbrite cycle. The activity of Vesuvius produced the building of a stratovolcano mostly by effusive and plinian explosive eruptions. The Campi Flegrei area, on the contrary, was dominated by the eruption of the Neapolitan Yellow Tuff at 12 ka that produced a caldera collapse of the Gulf of Pozzuoli. The caldera formation controlled the emplacement of the recent activity of Campi Flegrei and the new volcanoes were formed only within the caldera or along its rim.     

Evaluating fracture patterns within a resurgent caldera: Campi Flegrei, Italy

Understanding deformation of active calderas allows their dynamics to be defined and their hazard mitigated. The Campi Flegrei resurgent caldera (Italy) is one of the most active and hazardous volcanoes in the world, characterized by post-collapse resurgence, eruptions, ground deformation, and seismicity. An original structural analysis provides an overview of the main fracture zones. NW-SE and NE-SW fractures (normal or transtensive faults and extensional fractures) predominate along the rim and within the caldera, suggesting a regional control, both during and after the collapses. While the NE-SW fractures are ubiquitous in the deposits of the last ∼37 ka, NW-SE fractures predominate in the last 4.5 ka, during resurgence. The most recently (<4.5 ka) strained area lies in the caldera center (Solfatara area), where the faults, with an overall ∼ENE-WSW extension direction, appear to be associated with the bending due to resurgence. Solfatara lies immediately to the east of the most uplifted part of the caldera (Pozzuoli area), where domes form and culminate both on the long-term (resurgence, accompanied by volcanic activity) and short-term deformation (1982–1984 bradyseism, accompanied by seismic and hydrothermal activity). Similar volcano-tectonic behavior characterizes the short- and long-term uplifts, and only the intensity of the tectonic and volcanic activity varies, being related to varying amounts of uplift. Seismicity and hydrothermal manifestations occur during the bradyseisms, with moderate uplift, while surface faulting and eruptions occur during resurgence, with higher uplift. The features observed at Campi Flegrei are found at other major calderas, suggesting consistent behavior of large magmatic systems.     

History of earthquakes and vertical ground movement in Campi Flegrei caldera, Southern Italy: comparison of precursory events to the A.D. 1538 eruption of Monte Nuovo and of activity since 1968

The record of felt earthquakes around Naples Bay in southern Italy is probably complete since the mid-15th century. According to this record, intense earthquake swarms originating beneath Campi Flegrei, an explosive caldera located along the north coast of Naples Bay, have occurred only twice: (1) before the only historical eruption in Campi Flegrei in 1538; and (2) from mid-1983 to December 1984. Earthquake activity during the earlier period, which began at least a few years, and possibly as many as 30 years, before the 1538 eruption, damaged many buildings in the city of Pozzuoli, located near the center of Campi Flegrei. Minor seismic activity, which consisted of only a few felt earthquakes, occurred from 1970 to 1971. The second period of intense earthquake swarms lasted from mid-1983 to 1984, again damaging many buildings in Pozzuoli. Two periods of uplift along the shoreline within Campi Flegrei have also been noted since the mid-15th century: (1) during the few decades before the 1538 eruption; and (2) as two distinct episodes since 1968. Uplift of a few meters probably occurred a few decades before the 1538 eruption; uplift of as much as 3.0 m has occurred in Pozzuoli since 1968.These similarities strongly suggest that, for the first time in 440 years, the same process that caused intense local earthquake swarms and uplift in the early 1500's and led to an eruption in 1538, has again occurred beneath Campi Flegrei. Though no major seismicity or uplift has occurred since December 1984, because of the large amount of extensional strain accumulated during the past two decades, if a third episode of seismicity and rapid uplift occurs, it may lead to an eruption within several months after the resumption of activity.     

Vertical ground movements in the Campi Flegrei caldera as a chaotic dynamic phenomenon

The ground level in the Campi Flegrei caldera has never been stationary in the last 2,000 years. Historical data, and a nearly continuous tide-gauge record 20 years long, show that uplift and sinking have taken place on a variety of different time scales. In addition, the Campi Flegrei volcanic system appears to be sensitive to weak external forces such as tidal forces. We infer from these elements that the Campi Flegrei system is far from thermodynamic equilibrium, and suggest that its dynamics may be chaotic. We analyze the short-term variations of the ground level, and find that they can be described in a low-dimensional phase space. The dynamics of the Campi Flegrei system seems to have been phase-locked with tidal forces in the period following the 1970–1972 climax, and to have undergone a transition to chaos in some moment that preceded the presently continuing sinking phase.     

<Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of tuff vents in the Campi Flegrei caldera (southern Italy): toward a new chronostratigraphic reconstruction of the Holocene volcanic activity

The Campi Flegrei hosts numerous monogenetic vents inferred to be younger than the 15 ka Neapolitan Yellow Tuff. Sanidine crystals from the three young Campi Flegrei vents of Fondi di Baia, Bacoli and Nisida were dated using ⁴⁰Ar/³⁹Ar geochronology. These vents, together with several other young edifices, occur roughly along the inner border of the Campi Flegrei caldera, suggesting that the volcanic conduits are controlled by caldera-bounding faults. Plateau ages of ∼9.6 ka (Fondi di Baia), ∼8.6 ka (Bacoli) and ∼3.9 ka (Nisida) indicate eruptive activity during intervals previously interpreted as quiescent. A critical revision, involving calendar age correction of literature ¹⁴C data and available ⁴⁰Ar/³⁹Ar age data, is presented. A new reference chronostratigraphic framework for Holocene Phlegrean activity, which significantly differs from the previously adopted ones, is proposed. This has important implications for understanding the Campi Flegrei eruptive history and, ultimately, for the evaluation of related volcanic risk and hazard, for which the inferred history of its recent activity is generally taken into account.     

Procida volcanic history: new insights into the evolution of the Phlegraean Volcanic District (Campania region,Italy)

New stratigraphic, major- and trace-element, and Sr-, Nd- and Pb- isotopic data on volcanic deposits older than 14 ka from the island of Procida, Italy, are presented and compared with published analyses from the rest of the Phlegraean Volcanic District (PVD). Procida rocks range in composition from basalt to shoshonite and trachyte and show ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb ratios varying within the ranges 0.70523–0.70678, 0.512539–0.512630, 18.99–19.29, 15.67–15.69 and 39.10–39.39, respectively. The presence of a compositional gap in the range SiO₂=54–59 wt % is evidence of magma bimodality, and suggests that the feeding magmatic system was formed by at least two different reservoirs located at different depths. Geochemical and isotopic variations with increasing differentiation can be explained by fractional crystallization mechanisms, that in some cases were associated with crustal contamination that occurred in both deeper and shallower reservoirs; the most evolved magmas formed in the shallower one. Mixing/mingling processes also occurred. The variation of isotopic composition through time observed both for Procida and for Campi Flegrei and Ischia rocks is evidence of strong affinity between magmas that erupted on the entire PVD until about 40 ka. This indicates that they share a common origin and a common plumbing system. Most of the PVD eruptive centers active until about 40 ka lie within a NE-SW-oriented volcano-tectonic belt extending from the southeastern part of Ischia, through Procida and Torregaveta volcano to the northeastern sector of the present Campi Flegrei caldera. This not only indicates the existence of a link between regional structures and volcanism in the area, but also suggests that PVD magma genesis and evolution were strongly regulated by extensional tectonics. In the last 40 ka the mafic rocks erupted along this extensional structure – from Torregaveta and the islands of Ischia and Procida – indicate that it still represents an important crustal discontinuity that focuses mantle-derived magmas. Procida trachybasalts are enriched in large ion lithophile elements (LILE) and light and middle rare earth elements (LREE and MREE), and show slight negative anomalies in the high field strength elements (HFSE) relative to average MORB. A slight depletion in HREE is present. Trace element and isotope systematics can be referred to a lithospheric mantle source. The lithospheric mantle source carries intra-plate and slab-derived components, the latter probably inherited from a previous subduction event.     

Volcanic CO2 flux measurement at Campi Flegrei by tunable diode laser absorption spectroscopy

Near-infrared room temperature tunable diode lasers (TDL) have recently found increased usage in atmospheric chemistry and air monitoring research, but applications in volcanology are still limited to a few examples. Here, we explored the potential of a commercial infrared laser unit (GasFinder 2.0 from Boreal Laser Ltd) for measurement of volcanic CO₂ mixing ratios, and ultimately for estimating the volcanic CO₂ flux. Our field tests were conducted at Campi Flegrei near Pozzuoli, Southern Italy, where the GasFinder was used during three campaigns in October 2012, January 2013 and May 2013 to repeatedly measure the path-integrated mixing ratios of CO₂ along cross sections of the atmospheric plumes of two major fumarolic fields (Solfatara and Pisciarelli). By using a tomographic post-processing routine, we resolved, for each of the two fields, the contour maps of CO₂ mixing ratios in the atmosphere, from the integration of which (and after multiplication by the plumes’ transport speeds) the CO₂ fluxes were finally obtained. We evaluate a total CO₂ output from the Campi Flegrei fumaroles of ～490 Mg/day, in line with independent estimates based on in situ (Multi-GAS) observations. We conclude that TDL technique may enable CO₂ flux quantification at other volcanoes worldwide.     

The Neapolitan Yellow Tuff — A large volume multiphase eruption from Campi Flegrei,Southern Italy

the Neapolitan Yellow Tuff (NYT) (12 ka BP) is considered to be the product of a single eruption. Two different members (A and B) have been identified and can be correlated around the whole of Campi Flegrei. Member A is made up of at least 6 fall units including both ash and lapilli horizons. The basal stratified ash unit (A1) is interpreted to be a phreatoplinian fall deposit, since it shows a widespread dispersal (>1000 km²) and a constant thickness over considerable topography. The absence of many lapilli fall units in proximal and medial areas testifies to the erosive power of the intervening pyroclastic surges. The overlying member B was formed by many pyroclastic flows, radially distributed around Campi Flegrei, that varied widely in their eruptive and emplacement mechanisms. In some of the most proximal exposures coarse scoria and lithic-rich deposits, sometimes welded, have been identified at the base of member B. Isopach and isopleth maps of fall-units, combined with the distribution of the coarse proximal facies, indicate that the eruptive vent was located in the NE area of Campi Flegrei. It is considered that the NYT eruption produced collapse of a caldera approximately 10 km diameter within Campi Flegrei. The caldera rim, located by geological and borehole evidence, is now largely buried by the products of more recent eruptions. Initiation of caldera collapse may have been contemporaneous with the start of the second phase (member B). It is suggested that there was a single vent throughout the eruption rather than the development of multiple or ring vents. Chemical data indicate that different levels of a zoned trachyte-phonolite magma chamber were tapped during the eruption. The minimum volume of the NYT is calculated to be about 50 km³ (DRE), of which 35 km³ (70%) occurs within the caldera.     

Seismic images and rock properties of the very shallow structure of Campi Flegrei caldera (southern Italy)

In September 2001, an extensive active-seismic investigation (Serapis experiment) was carried out in the Gulfs of Naples and Pozzuoli, with the aim of investigating and reconstructing the shallow crustal structure of the Campi Flegrei caldera, and possibly identifying its feeding system at depth. The present study provides a joint analysis of the very shallow seismic reflection data and tomographic images based on the Serapis dataset. This is achieved by reflection seismic sections obtained by the 3D data gathering and through refined P-velocity images of the shallowest layer of Pozzuoli Gulf (z < 1,000 m). From the refined Vp model, the overall picture of the velocity distribution confirms the presence of a complex arc-shaped anomaly that borders the bay offshore. The deeper part of the anomaly (beneath 700 m, with Vp > 3,500 m/s) correlates with units made up of agglomerate tuffs and interbedded lava, which form the southern edge of the caldera, which was probably formed following the two large ignimbritic eruptions that marked the evolutionary history of the area under study. The upper part of the anomaly that tends to split into two parallel arcs is correlated with dikes, volcanic mounds and hydrothermal alteration zones noted in previous shallow reflection seismic analyses. The depth of the transition between the upper and lower parts of the anomaly is characterized by an abrupt Vp increase on the one-dimensional (1D) profiles extracted from the 3D tomographic model and by the presence of a strong reflector located at about 0.6/0.7 s Two Way Time (TWT) on Common Mid Point gathers. The move-out velocity analysis and stack of the P–P and P–S reflections at the layer bottom allowed to estimate relatively high Vp/Vs values (3.7 ± 0.9). This hypothesis has been tested by a theoretical rock physical modeling of the Vp/Vs ratio as a function of porosity suggesting that the shallow layer is likely formed by incoherent, water saturated, volcanic and marine sediments that filled Pozzuoli Bay during the post-caldera activity.     

Stress pattern at Campi Flegrei from focal mechanisms of the 1982–1984 earthquakes (southern Italy)

Campi Flegrei is a Holocene volcanic area in the Campanian Plain (Southern Italy) within the Apennine Chain, a neogenic thrust belt built up since the Middle Miocene. The volcanic complex consists of a c. 12-km-diameter caldera containing several monogenetic volcanoes, the youngest of which (Monte Nuovo) erupted in 1538. Since at least Roman times, the area has also been affected by slow vertical movements (bradiseismic activity). Between 1982 and 1985, this slow motion was interrupted by a period of more rapid displacement which caused a maximum uplift of 180 cm in the town of Pozzuoli. To model the local stress field associated with the uplift, the Angelier inversion technique has been applied using the focal mechanisms of forty-nine earthquakes which occurred between April 1982 and December 1984. The results show that doming coupled with a regional extensional strain can account for the seismic phenomena that affected the area.     

Hazard and risk assessment in a complex multi-source volcanic area: the example of the Campania Region, Italy

In order to zone the territory of Campania Region (southern Italy) with regard to the hazard related to future explosive activity of Somma-Vesuvio, Campi Flegrei, and Ischia Island, we drew a multi-source hazard map for tephra and pyroclastic flows. This map, which merges the areas possibly endangered by the three volcanic sources, takes into account a large set of tephra fall and pyroclastic flow events that have occurred in the last 10 ka. In detail, for fall products at Campi Flegrei and Somma-Vesuvio we used the dispersal of past eruption products as deduced by field surveys and their recurrence over the whole area. For pyroclastic flows, the field data were integrated with VEI = 4 simulated events; about 100 simulations sourcing from different points of the area were performed, considering the different probability of vent opening. The spatial recurrence of products of both past eruptions and simulated events was used to assign a weight to the area endangered by the single volcanic sources. The sum of these weights in the areas exposed to the activity of two sources and/or to different kinds of products was used to draw a hazard map, which highlights the spatial trend and the extent of the single equivalent classes at a regional scale. A multi-source risk map was developed for the same areas as the graphic result of the product of volcanic hazard and exposure, assessed in detail from a dasymetric map. The resulting multi-source hazard and risk maps are essential tools for communication among scientists, local authorities, and the public, and may prove highly practical for long-term regional-scale mitigation planning.     

A methodology for the evaluation of long-term volcanic risk from pyroclastic flows in Campi Flegrei (Italy)

The volcanological history of Campi Flegrei suggests that the most frequent eruptions are characterized by the emplacement of pyroclastic flow and surge deposits erupted from different vents scattered over a 150-km² caldera. The evaluation of volcanic risk in volcanic fields is complex because of the lack of a central vent. To approach this problem, we subdivided the entire area of Campi Flegrei into a regular grid and evaluated the relative spatial probability of opening of vents based on geological, geophysical and geochemical data. We evaluated the volcanic risk caused by pyroclastic flows based on the formula proposed by UNESCO (1972), R=H×V×Va, where H is the hazard, V is the vulnerability and Va is the value of the elements at risk. The product H×V was obtained by performing simulations of type eruptions centered in each cell of the grid. The simulation is based on the energy cone scheme proposed by Sheridan and Malin [J. Volcanol. Geotherm. Res. 17 (1983) 187–202], hypothesizing a column collapse height of 100 m for eruptions of VEI=3 and 300 m for eruptions of VEI=4 with a slope angle of 6°. Each simulation has been given the relative probability value associated with the corresponding cell. We made use of the GIS software ArcView 3.2 to evaluate the intersection between the energy cone and the topography. The superposition of the areas invaded by pyroclastic flows (124 simulations for VEI=3 and 37 for VEI=4) was used to obtain the relative hazard map of the area. The relative volcanic risk map is obtained by superimposing the urbanization maps.     

Vent area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (Campi Flegrei, Italy): new insights from directional fabric through image analysis

In order to provide new information about the source area and depositional mechanisms of the Upper Member of the Neapolitan Yellow Tuff (NYT), a prominent pyroclastic deposit of the Campi Flegrei Volcanic District (southern Italy), statistics on directional fabric, by means of computer-assisted image analysis on 32 rock samples, were compiled. Seventeen samples were collected along vertical direction on two selected exposures and fifteen were taken from outcrops widely distributed all around the Campi Flegrei Volcanic District. Fabric measurements within the investigated successions reveal a vertically homogeneous direction of the mean particle iso-orientation, with considerable variability in the strength of particle iso-orientation even at cm-scale. The existence of particle iso-orientation can be related to continuous sedimentation from a concentrated bedload region beneath suspension currents, producing massive or inversely graded beds by traction carpet sedimentation. The considerable vertical variability in the strength of iso-orientation is the result of very unstable flow regimes, up to the extreme condition of discrete depositional events, with a variable combination of traction carpet and/or direct suspension sedimentation. The vertical homogeneity in the mean orientation values, found in the investigated sections, may derive from the sequential deposition of laminae to thin beds, whose relatively flat upper surfaces were unable to significantly deflect the depositional system of the following currents. According to the observed homogeneous mean particle orientation values along the investigated vertical profiles, samples collected through areal distribution are considered representative of the local paleo-flow directions of the whole deposit. The mean directions of the samples collected areally show two different coherent patterns which point to the existence of two different source areas. The first, which includes all samples from the northern outcrops, appears to converge in a narrow area about 2 km NE of the town of Pozzuoli, largely in coincidence with the inferred area on the basis of the pumice fall distribution. The second, which includes samples from Capo Miseno and Posillipo areas, points to the central part of the Pozzuoli Bay, about 4 km offshore the town of Pozzuoli.     

Statistically coherent approach involving log-ratio transformation of geochemical data enabled tephra correlations of two late Pleistocene tephra from the eastern Adriatic shelf

Tephrostratigraphic correlations commonly rely on geochemical composition supported by additional constraints (e.g., multiple stratigraphically ordered tephra, geochronological-stratigraphical constraints, and isotopic determinations), which provide key clues to restrict the number of possible candidates and disambiguate the correlation of a specific tephra among compositionally similar volcanic sources/tephra. However, such additional data may not be available or acquirable, leaving the geochemical data as sole, but challenging viable approach. In this study, two geochronologically poorly constrained late Pleistocene tephra from the eastern Adriatic –from a sand profile on Mljet Island (M-53/2) and from a marine sediment core from Pirovac Bay (PROS 721)– were correlated to known eruptions using only geochemical data (major and trace elements of glass shards), which were treated using both log ratio transformed and raw data. After the statistical treatment of the geochemical data using bivariate plots, linear discrimination analyses and selbal algorithm, the tephra M-53/2 and PROS 721 were suitably correlated with the widespread tephra generated during the Campi Flegrei eruptions of Massereia del Monte (Y-3 marine tephra, 29.0 ± 0.8 ka) and Neapolitan Yellow Tuff (14.5 ± 0.4 ka), respectively. This study showed that the correlation was hardly tenable when using the raw data, as opposed to compositional approach, which yielded satisfactory results. As a consequence, the distribution of Massereia del Monte/Y-3 tephra extended far toward the northeast, while a better chronological model, for reconstructing the paleoenvironmental changes at the Pirovac Bay location and the Holocene sea-level dynamics, could be obtained.     

Volcanic hazard assessment at the restless Campi Flegrei caldera

Eruption forecasting and hazard assessments at the restless Campi Flegrei caldera, within the Neapolitan volcanic area, have been performed using stratigraphical, volcanological, structural and petrological data.On the basis of the reconstructed variation of eruption magnitude through time, we hypothesize that the most probable maximum expected event is a medium-magnitude explosive eruption, fed by trachytic magma. Such an eruption could likely occur in the north-eastern sector of the caldera floor that is under a tensile stress regime, when the ongoing deformation will generate mechanical failure of the rocks. A vent could open also in the western sector, at the intersection of two fault systems contemporaneously activated, as happened in the last eruption at Monte Nuovo. The eruption could likely be preceded by precursors apparent to the population, such as ground deformation, seismicity and increase in gas emissions. It will probably alternate between magmatic and phreatomagmatic phases with the generation of tephra fallout, and dilute and turbulent pyroclastic currents. During and/or after the eruption, the re-mobilization of ash by likely heavy rains, could probably generate mud flows.In order to perform a zoning of the territory in relation to the expected volcanic hazards, we have constructed a comprehensive hazard map. On this map are delimited (I) areas of variable probability of opening of a new vent, (II) areas which could be affected by variable load of fallout deposits, and (III) areas over which pyroclastic currents could flow. The areas in which a vent could likely open have been defined on the basis of the dynamics of the ongoing deformation of the caldera floor. To construct the fallout hazard map we have used the frequency of deposition of fallout beds thicker than 10 cm, the frequency of load on the ground by tephra fallout and the direction of dispersal axes of the deposits of the last 5 ka, and the limit load of collapse for the variable types of roof construction. The pyroclastic-current hazard map is based on the areal distribution and frequency of pyroclastic-current deposits of the last 5 ka.Editorial Responsibility: T. Druitt     

Probability hazard map for future vent opening at the Campi Flegrei caldera,Italy

The Campi Flegrei caldera is a restless structure affected by general subsidence and ongoing resurgence of its central part. The persistent activity of the system and the explosive character of the volcanism lead to a very high volcanic hazard that, combined with intense urbanization, corresponds to a very high volcanic risk. One of the largest sources of uncertainty in volcanic hazard/risk assessment for Campi Flegrei is the spatial location of the future volcanic activity. This paper presents and discusses a long-term probability hazard map for vent opening in case of renewal of volcanism at the Campi Flegrei caldera, which shows the spatial conditional probability for the next vent opening, given that an eruption occurs. The map has been constructed by building a Bayesian inference scheme merging prior information and past data. The method allows both aleatory and epistemic uncertainties to be evaluated. The probability map of vent opening shows that two areas of relatively high probability are present within the active portion of the caldera, with a probability approximately doubled with respect to the rest of the caldera. The map has an immediate use in evaluating the areas of the caldera prone to the highest volcanic hazard. Furthermore, it represents an important ingredient in addressing the more general problem of quantitative volcanic hazards assessment at the Campi Flegrei caldera.     

Hazard assessment at volcanic fields: the Campi Flegrei case history

Volcanological analysis of the 10 000 yr –1538 explosive activity at Campi Flegrei shows that the most common explosive eruptions are characterized by the emplacement of flow or surge deposits, originating from the interaction between magma and shallow and/or sea water. The minimum volumes of pyroclastic products range between 0.04 and 0.7 km³; the proximal areas covered by these products range from 3–4 to 40–50 km². The pyroclastic flow and surge deposits occurring inside the caldera have been strongly controlled by pre-existent morphology; because of this, the area of present Napoli city was blanketed by approximately 5 m of pyroclastic deposits, during the last 5000 yr.Previous analysis suggests that the presence of even very low topographic obstacles may influence pyroclastic density current run out such that future eruptive deposits would mainly be confined inside the caldera rim. We suggest that a future eruption at Campi Flegrei would not seriously involve the urbanized area of Napoli city located on the hills. On the contrary, the plains located on the eastern side of the caldera (Fuorigrotta, Bagnoli) would be the most damaged area.     

Magma chamber evolution prior to the Campanian Ignimbrite and Neapolitan Yellow Tuff eruptions (Campi Flegrei,Italy)

The Campi Flegrei (Campanian Region, Italy) experienced two cataclysmic caldera-forming eruptions which produced the Campanian Ignimbrite (39 ka, CI) and the Neapolitan Yellow Tuff (15 ka, NYT). We studied the minor eruptions before both these large events to understand magma chamber evolution leading towards such catastrophic eruptions. Major, trace element, and Sr and Nd isotope compositions of pre-Campanian Ignimbrite and pre-Neapolitan Yellow Tuff products define distinct geochemical groups, which are here interpreted as distinct magma batches. These batches do not show any transitional trend towards the CI and NYT eruptions. The CI and NYT systems are decoupled geochemically and isotopically. At least one of the pre-CI and one of the pre-NYT erupted magma batches qualifies as mixing endmembers for the large CI and NYT eruptions, and thus, must have been stored in reservoirs for some time to remain available for the CI and NYT eruptions. The least evolved, isotopically distinct magma compositions that are typical of the last phases of the NYT and CI eruptions did not occur before caldera-forming events. Based on the new data, we propose the following scenario: Multiple magma chambers with distinct compositions existed below the Campi Flegrei before the CI and NYT eruptions and remained generally separated for some time unless new magma was recharged. In each case, one of the residing magma reservoirs was recharged by a new large-volume magma input of intermediate composition from a deeper differentiating magma reservoir. This may have triggered the coalescence of the previously separated reservoirs into one large chamber which fed the cataclysmic caldera-forming eruption. Large magma chambers in the Campi Flegrei may therefore be ephemeral features, interrupted by periods of evolution in individual, separated magma reservoirs.     

Radon variation in groundwaters in the Campi Flegrei Caldera (southern Italy) during and after the 1982–1984 bradyseismic crisis

Groundwater Rn concentration has been monitored at several sites in the Campi Flegrei Caldera since the time of the last episode of bradyseismic activity. The results concerning the period 1983–88 indicate that the Rn concentration is spatially variable due to the different chemico-physical characteristics of the waters or structural features. Rn concentration is generally influenced by the long-term changes of meteorological factors. Rainfalls appear to affect the Rn content of the most superficial aquifers. Several short-term spike-like Rn anomalies have been detected both during and after the period of the bradyseismic crisis. They do not show correlation with temporal changes of the pH and/or temperature of waters. A systematic correlation between Rn anomalies and episodic strong rainfall events is also absent. Furthermore, a clear relationship between Rn spikes and both uplift velocity changes and seismic activity has not been found.The overall evidence makes it very difficult to understand the temporal changes of groundwater Rn concentration that have occurred during and after the 1982–84 bradyseismic crisis. Nevertheless, a clear causal relationship between Rn concentration and the dynamics of the bradyseism is observed in some monitored sites. Its appearance may have been favoured by the low Rn concentration in the groundwaters during periods of quiescent activity.     

The Breccia Museo formation, Campi Flegrei, southern Italy: geochronology, chemostratigraphy and relationship with the Campanian Ignimbrite eruption

The Breccia Museo is one of the most debated volcanic formations of the Campi Flegrei volcanic district. The deposit, made up of six distinctive stratigraphic units, has been interpreted by some as the proximal facies of the major caldera-forming Campanian Ignimbrite eruption, and by others as the product of several, more recent, independent and localized events. New geochemical and chemostratigraphical data and Ar–Ar age determinations for several units of the Breccia Museo deposits (~39 ka), correlate well with the Campanian Ignimbrite-forming eruption. The chemical zoning of the Breccia Museo deposits is interpreted here to be a consequence of a three-stage event that tapped a vertically zoned trachytic magma chamber. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.