Campi Flegrei Research Articles

Compositional variation and zoning of epidote supergroup minerals in the Campi Flegrei geothermal field, Naples, Italy

Abstract. Authigenic epidote supergroups are an abundant accessory mineral in the calcium–aluminum silicate and thermometamorphic hydrothermal zones of the Campi Flegrei (Phlegraean Fields) geothermal field located west of Naples, Italy. Geothermal exploration for high-enthalpy fluid produced drill core and cuttings to ∼ 3 km depth in the Mofete (MF1, MF2, MF5) and San Vito (SV1, SV3) wells, where measured down-hole temperatures of epidote-bearing samples range from 270–350 ∘C and from 285–390 ∘C for the Mofete and San Vito areas, respectively. Two epidote group (epidote, clinozoisite), some rare earth element (REE)-bearing, and two allanite group (allanite-(Ce), ferriallanite-(Ce)) minerals were identified by electron microprobe. The allanite group is light rare earth element (LREE, La–Gd) enriched, Ce dominant, with REE + Y that varies from 30.59 wt %–14.32 wt %. Complex compositional variation such as oscillatory, sector, and complex (mixed) zoning is a ubiquitous feature observed in the epidote group, which occurs as veins, in vugs, as various size masses, and as isolated single crystals. Compositional zoning is caused by variable Fe ↔ Al3+ substitution and XFe [(Fe3+) / (Fe3++ Al)] ranges from 0.06–0.33 (Fe3+=0.185–0.967 apfu). XFe tends to decrease with increasing temperature in the Mofete wells, but its distribution is more complex in the San Vito wells, which records recent fault displacement. The variety and complexity of the epidote supergroup zoning suggest rapid fluid composition and/or intensive parameter fluctuations in the local hydrothermal system.

Bouguer Anomaly Re-Reduction and Interpretative Remarks of the Phlegraean Fields Caldera Structures (Southern Italy)

Phlegraean Fields is a large, active caldera located in the densely populated westernmost sector of Naples’s Bay (Southern Italy). Several Bouguer anomaly surveys are available for this area with different resolution and accuracy; gravity data derive from the integration of stations placed below and above the sea level as the caldera develops both onshore and offshore. The comparison of these maps with the Digital Elevation Model shows a still remaining Terrain Effect hiding the shallower and deep caldera structure’s signal. This effect has an impact on the modelling of the gravity source’s depth and geometry. In this research, we apply a geologically constrained terrain correction method to the higher resolution Free Air dataset available for the study area to enhance the complete Bouguer reduction. The correlation analysis between the residual and the topography allows us to assess the quality of the outcomes. The results represent an improvement in the anomalies’ isolation and clearly show a continuous circular-like clustering of maxima related to the geometry of the caldera rim. The minima are associated with volcano-tectonic depression filled with pyroclastic and sediment. Furthermore, features alignments overlap the fault systems, along which the volcanic activity occurred.

The Effectiveness of Geotrails to Support Sustainable Development in the Campi Flegrei Active Volcanic Area

The Effectiveness of Geotrails to Support Sustainable Development in the Campi Flegrei Active Volcanic Area

Space-time analysis of informational properties of GPS time series recorded at the Campi Flegrei caldera (Italy)

In this work the time dynamics of GPS time series recorded at the Campi Flegrei caldera from 2000 to 2019 was investigated by using the Fisher Information Measure (FIM) and the Shannon entropy power (SEP), two informational methods that allow the detection of changes in the dynamical behavior of a complex system, like the volcanic one, quantifying respectively the order/organization and the disorder/uncertainty. By jointly analyzing the FIM and the SEP through the Fisher-Shannon Information Plane (FSIP) it was shown that the SEP discriminates quite well the vertical displacement time series, while the FIM discriminates better the horizontal ones, which are featured by a more complex behavior. Globally, the most striking change in the SEP is found in 2012, while at local scale abrupt loss of order (indicated by low values of FIM) well correlate with the occurrence of seismic swarms. A further informational quantity, the complexity C (given by the product between SEP and FIM), has been revealed to be sensitive to the baseline deformation level increases and to the small-amplitude modulation of deformation signals within the caldera.

Numerical evaluation of thermal and hydrodynamic effects caused by heat production well on geothermal Phlegraean Fields

This study describes the geothermal response of the Phlegraean Fields as well as the impact of changes in its thermal and hydrodynamic properties brought on by a deep borehole heat exchanger (DBHE). For this purpose, we have developed a specialized model based on the Galerkin Method (GM) and the iterative Newton–Raphson algorithm to perform a transient simulation of heat transfer with fluid flow in porous media by solving the related system of coupled non-linear differential equations. A two-dimensional domain characterized with an anisotropic saturated porous media and a non-uniform grid is simulated. Extreme characteristics, such as non-uniformity in the distribution of the thermal source, are implemented as well as the fluid flow boundary conditions. While simulating the undisturbed geothermal reservoir and reaching the steady temperature, stream function, and velocity components, a DBHE is placed into the domain to evaluate its impact on the thermal and fluid flow fields. This research aims to identify and investigate the variables involved in the Phlegraean Fields and provide a numerical approach to accurately simulate the thermodynamic and hydrodynamic effects induced in a reservoir by a DBHE. The results show a maximum temperature change of 107.3°C in 200 years of service in the study area and a 65-year time limit is set for sustainable geothermal energy production.

3D structure of the Campi Flegrei caldera central sector reconstructed through short-period magnetotelluric imaging

The Campi Flegrei caldera experienced an unrest phase dating to 2005, which primary expression is the impressive ground uplift, accompanied by increasing degassing and seismic activities. Such last two phenomena developed mainly in the caldera central sector, including the Solfatara–Pisciarelli complex. However, the inner structure of such an area is still not defined, and this originates a poor understanding of the ongoing unrest. This paper describes the results of a new magnetotelluric survey performed in the Campi Flegrei caldera central sector. Through the inversion of data collected in 47 independent soundings, a 3D model of the electrical resistivity has been retrieved, which evidenced a partition of the investigated structure. The Agnano–Astroni area seems to be associated with a liquid-dominated geothermal reservoir, whereas the Solfatara–Pisciarelli area seems to be characterized by a single mixed liquid and gasses-dominated geothermal reservoir, which supplies the main caldera fumaroles. The proposed reconstruction of the geometrical characteristics of the hydrothermal system and the primary fluid rising pathways gives substantial clues about the significance of the detected structures in the evolution of the caldera unrest.

Structure and Present State of the Astroni Volcano in the Campi Flegrei Caldera in Italy Based on Multidisciplinary Investigations

AbstractDespite its known reconstructed volcanic history, the structural setting and present state of the Astroni Volcano of the Campi Flegrei caldera in Italy are still poorly defined. Through structural, geophysical, and geochemical investigations, we elucidate the structure and present volcanic activity of the Astroni Volcano, which hosts tuff cones, scoriae cones, lava domes, and lakes on the crater floor. A volcano‐tectonic analysis focused on the entire volcano edifice, coupled with electrical resistivity tomography of the shallower part of the Astroni crater, revealed the main rock formations, faults, and possible fluid patterns within the first 150 m depth. Two main NE–SW and NW–SE trending fault sets were imaged using electrical resistivity modeling and measurements along the wall of the volcanic edifice; they likely delimit a maar‐like structure resulting from the highest energetic subplinian Astroni 6 eruption event and acted as magma pathways during the late eruptive activity stage. A 3D view of the reconstructed resistivity model revealed both deep root‐conduit‐like structures and shallower dome‐like shapes for volcanic edifices on the crater floor. Gas and carbon compositions in the NNE sectors of the Astroni Lago Grande are similar to those of the Solfatara fumarole fluids, suggesting common hydrothermal origin and a possible link with a deep hydrothermal reservoir. This fluid‐emission area along the border of the younger volcanic structure exhibits a +40°C maximum soil‐temperature anomaly. The proposed volcano‐tectonic architecture should improve the unrest scenarios in case of reactivation in this Campi Flegrei caldera sector and the monitoring strategies for the Astroni Volcano.

Thirty years of volcano geodesy from space at Campi Flegrei caldera (Italy)

This work provides the mean ground deformation rates and ground displacement time series of the Campi Flegrei caldera (Italy) retrieved by satellite remote sensing data analysis from 1992 to 2021. Synthetic Aperture Radar (SAR) images acquired by ERS 1–2 (1992–2002), Envisat (2003–2011) and Cosmo-SkyMed (2011–2021) are processed by multi-temporal SAR Interferometry (InSAR) approach using the same technique, parameters and reference system, to obtain for the first time a homogeneous and time-continuous dataset. The validation of the InSAR products is carried out by comparison with the measurements provided by precise levelling lines and cGNSS stations. The produced outcomes offer an overview on the temporal behaviour of ground deformation at Campi Flegrei along an unprecedented time window of about 30 years and can be exploited by the scientific community for supporting and improving the knowledge of the dynamics of the caldera.

Potentially toxic elements in the soils of Campi Flegrei (south Italy) and the immediate surroundings: Spatial distribution, origin and probabilistic human health risk

Campi Flegrei is an active volcanic field in south Italy where the potentially toxic elements (PTEs) are of growing concern because the intensive anthropogenic and volcanic activities might pose adverse human health effects. In this article, 394 topsoils (0–15 cm) are collected for instrumental analysis of the <2 mm fraction. The geochemical maps indicate that higher concentrations of Pb, Zn, Cd, Cr, Hg, Ni and Sb are related to the urban area, but greater levels of As, Tl, Co, Cu, Se and V are observed in the other parts. A robust principal component analysis detected: (1) the Pb–Zn–Hg–Cd–Sb–Cr–Ni association that probably highlights anthropogenic activities such as heavy traffic load and fossil fuel combustion in the urbanized area; (2) the Al–Fe–Mn–Ti–Tl–V–Co–As–U–Th association that mostly reveals the contribution of pyroclastic deposits; and (3) the Na–K–B association that feasibly indicates the weathering degree. The probabilistic health risk modeling for the children under 6 years old shows negligible Pb and Zn non-carcinogenic risk and unexpected Pb carcinogenic risk for exposure through soil ingestion. However, for the inhalation pathway, the children aged <1 year old have the highest chance (90%) of acceptable (i.e. between 1E-6 and 1E-4) Pb carcinogenic health risk. This should not be overlooked because Naples is under high environmental pressure and previous studies reported the increased Pb and Zn quantities in soil over a 26-year timespan. Overall, the results of geostatistical interpolation, compositional data analysis and probabilistic health risk modeling potentially uncover the link between soil geochemistry and human health.

Clues of Ongoing Deep Magma Inflation at Campi Flegrei Caldera (Italy) from Empirical Orthogonal Function Analysis of SAR Data

Understanding the characteristics of a volcanic system is always important and becomes crucial when the volcano is in pluri-decadal unrest and located in a densely populated area, such as Campi Flegrei caldera, Italy. Ground deformation is a powerful source of information for studying the activity of magmatic sources and hydrothermal systems, even when difficult to detect otherwise. Here, we use ground displacement from ERS-ENVISAT and Sentinel-1A SAR images to investigate the 1993–2000 subsidence and part (2015–2022) of the present unrest. For each of the two time periods, we combine the line-of-sight displacements to obtain vertical and eastward displacements, and apply the empirical orthogonal function analysis to these latter time series—treated as a single data set—to decompose space-time fields into separated modes, consisting of uncorrelated spatial patterns and associated temporal evolutions. We only retain the first mode, since it captures the main deformation during both investigated periods, is the sole mode related to long-lasting (years) processes, and is less affected by noise than original data. Our analyses: (i) Confirm that most of the deformation is related to the activity of a 3–4 km deep sill-like source, which is inflated by magma and/or magmatic fluids during periods of unrest and deflates during periods of subsidence; (ii) Evidence ongoing deformation linked to local fluid migration in the Solfatara area; (iii) Identify persistent deformation features where peculiar fluid migration processes occurred during the 1982–1984 unrest; (iv) Most importantly, provide direct evidence of deep magma inflation at least since 2015, thus giving a strong warning of increasing risk at Campi Flegrei. Results demonstrate the capability of our approach to disclose hardly detectable processes and suggest a tool to monitor the activity of the deep magmatic source. Our approach can be useful also in other volcanic systems.

Data analysis of the unsteadily accelerating GPS and seismic records at Campi Flegrei caldera from 2000 to 2020

Ongoing resurgence affects Campi Flegrei caldera (Italy) via bradyseism, i.e. a series of ground deformation episodes accompanied by increases in shallow seismicity. In this study, we perform a mathematical analysis of the GPS and seismic data in the instrumental catalogs from 2000 to 2020, and a comparison of them to the preceding data from 1983 to 1999. We clearly identify and characterize two overlying trends, i.e. a decennial-like acceleration and cyclic oscillations with various periods. In particular, we show that all the signals have been accelerating since 2005, and 90–97% of their increase has occurred since 2011, 40–80% since 2018. Nevertheless, the seismic and ground deformation signals evolved differently—the seismic count increased faster than the GPS data since 2011, and even more so since 2015, growing faster than an exponential function The ground deformation has a linearized rate slope, i.e. acceleration, of 0.6 cm/yr2 and 0.3 cm/yr2 from 2000 to 2020, respectively for the vertical (RITE GPS) and the horizontal (ACAE GPS) components. In addition, all annual rates show alternating speed-ups and slow-downs, consistent between the signals. We find seven major rate maxima since 2000, one every 2.8–3.5 years, with secondary maxima at fractions of the intervals. A cycle with longer period of 6.5–9 years is also identified. Finally, we apply the probabilistic failure forecast method, a nonlinear regression that calculates the theoretical time limit of the signals going to infinity (interpreted here as a critical state potentially reached by the volcano), conditional on the continuation of the observed nonlinear accelerations. Since 2000, we perform a retrospective analysis of the temporal evolution of these forecasts which highlight the periods of more intense acceleration. The failure forecast method applied on the seismic count from 2001 to 2020 produces upper time limits of [0, 3, 11] years (corresponding to the 5th, 50th and 95th percentiles, respectively), significantly shorter than those based on the GPS data, e.g. [0, 6, 21] years. Such estimates, only valid under the model assumption of continuation of the ongoing decennial-like acceleration, warn to keep the guard up on the future evolution of Campi Flegrei caldera.

Ancient Roman Mortars from Anfiteatro Flavio (Pozzuoli, Southern Italy): A Mineralogical, Petrographic and Chemical Study

This research is based on the characterization of ancient mortars from the Anfiteatro Flavio (Pozzuoli) dating back to the 1st and 2nd century CE through a multi-analytical approach involving macroscopic, mineralogical, petrographic, and chemical investigations. The goal that has been set is to deepen knowledge about mortar mix design, the provenance of used raw materials, and secondary minerogenetic processes that have occurred within ancient Roman mortars. Results show that: (i) raw materials for mortar preparation have a local provenance, i.e., Phlegraean Fields area; (ii) mortars can be considered as hydraulic; (iii) calcite presence could be due to a non-complete calcination process, an improper slaking or to exposition of materials to the subaerial environment; (iv) gypsum is due to calcite sulfation process; (v) halite presence is due to marine aerosol exposition. The achieved information testifies that, for at least two centuries, Roman builders considered the identified mortar mix as optimal for their buildings, but also contributes to the understanding of their technical skills and represents an important first step to planning future restoration operations.

A simple two-state model interprets temporal modulations in eruptive activity and enhances multivolcano hazard quantification.

Volcanic activity typically switches between high-activity states with many eruptions and low-activity states with few or no eruptions. We present a simple two-regime physics-informed statistical model that allows interpreting temporal modulations in eruptive activity. The model enhances comprehension and comparison of different volcanic systems and enables homogeneous integration into multivolcano hazard assessments that account for potential changes in volcanic regimes. The model satisfactorily fits the eruptive history of the three active volcanoes in the Neapolitan area, Italy (Mt. Vesuvius, Campi Flegrei, and Ischia) which encompass a wide range of volcanic behaviors. We find that these volcanoes have appreciably different processes for triggering and ending high-activity periods connected to different dominant volcanic processes controlling their eruptive activity, with different characteristic times and activity rates (expressed as number of eruptions per time interval). Presently, all three volcanoes are judged to be in a low-activity state, with decreasing probability of eruptions for Mt. Vesuvius, Ischia, and Campi Flegrei, respectively.

The petrogenesis of chemically zoned, phonolitic, Plinian and sub-Plinian eruptions of Somma-Vesuvius, Italy: Role of accessory phase removal, independently filled magma reservoirs with time, and transition from slightly to highly silica undersaturated magmatic series in an ultrapotassic stratovolcano

A mineralogical, major, LA-ICP-MS trace element mineral chemistry and bulk-rock geochemical study of juvenile samples of the Mercato, Avellino, Pompeii and Pollena eruptions, collected in stratigraphically and volcanologically well-characterized sections of the Somma-Vesuvius stratovolcano (Roman Magmatic Province), along with reference data on the 1944 CE and the Pomici di Base eruptions, highlights the compositional variability of bulk-rock and glass from leucite phonotephrite to garnet-bearing phonolite . The latter products have extreme fractionation of trace elements (e.g., La/Yb n = 126, Zr/Y = 89, Zr/Hf =78, Nb/Ta = 40; Th/U = 2.3), very low Sc, V, Y, HREE and very high As, Tl, Cs, Pb, Th and U. The Pomici di Base products, older than the eruptions described above, range from leucite-bearing shoshonites to trachytes , are devoid of garnet and belong to an independent liquid-line-of-descent, having also different fractionation between trace elements (e.g., La/Yb n = 15; Zr/Y = 12.4, Zr/Hf = 50, Nb/Ta = 15.6; Th/U = 3 in the Pomici di Base trachytes). A marked chemical variability of the observed phases is found. The geochemistry of garnet, amphibole , clinopyroxene and other phases shows wide variations of concentrations and elemental ratios (e.g., La/Yb n up to 520 in the sadanagaite coexisting with garnet). The magmatic evolution is dominated by low-pressure, oxidized, nearly closed-system fractional crystallization of clinopyroxene, plagioclase , leucite, ±magnetite, ±biotite, ±olivine and apatite in the transition tephrite-phonotephrite, and of potassic sanidine (±hyalophane), Fe-clinopyroxene, melanite garnet±Fe-amphibole in the transition tephriphonolite-phonolite. Mineralogic and geochemical evidence and modelling points out the existence of independent, zoned magma batches throughout the activity of the stratovolcano, which possibly started to crystallize at similar depths. The Somma-Vesuvius magmas thus evolved in shallow independent reservoirs with respect to those of the neighbouring volcanic complexes (Campi Flegrei, Ischia) as shown, for instance, by the contrasting compositional trends of clinopyroxene and amphibole, and have very limited evidence of crustal contamination (and/or carbonate assimilation). • Pollena, Pompeii, Avellino and Mercato are zoned eruptions with phonolites • The magmas which formed the various eruptions evolved independently • Different garnets dominate HREE, Y and HFSE; amphibole is strongly HREE-poor • The phonolites and their mineral compositions have no tendency to become peralkaline • F, Cl and H 2 O dominate the volatile budget

Experimental investigation of trace element partitioning between amphibole and alkali basaltic melt: Toward a more general partitioning model with implications for amphibole fractionation at deep crustal levels

Abstract Time-series experiments were carried out in a piston-cylinder apparatus at 0.8 GPa and 1030–1080 °C using a hydrous K-basalt melt as the starting material to determine the element partition coefficients between amphibole and silicate glass. Major, minor, and trace element compositions of amphibole and glass were determined with a combination of electron microprobe and laser ablation inductively coupled plasma mass spectrometry. Results show that the main mineral phase is calcic amphibole, and the coexisting glass compositions range from basaltic trachyandesite to andesite. We estimated the ideal radius, the maximum partition coefficient and the apparent Young’s modulus of the A, M1-M2-M3, and M4-M4′ sites of amphibole. The influence of melt and amphibole composition, temperature, and pressure on the partition coefficients between amphiboles and glasses has also been investigated by comparing our data with a literature data set spanning a wide range of pressures (0.6–2.5 GPa), temperatures (780–1100 °C), and compositions (from basanite to rhyolite). Finally, we modeled a deep fractional crystallization process using the amphibole-melt partition coefficients determined in this study, observing that significant amounts of amphibole crystallization (&amp;gt;30 wt%) well reproduce the composition of an andesitic melt similar to that of the calc-alkaline volcanic products found in Parete and Castelvolturno boreholes (NW of Campi Flegrei, Italy).

Submarine Stratigraphy of the Eastern Bay of Naples: New Seismo-Stratigraphic Data and Implications for the Somma-Vesuvius and Campi Flegrei Volcanic Activity

The submarine stratigraphy of the eastern Bay of Naples is studied through seismo-stratigraphic data correlated with borehole data. Multichannel seismic profiles are interpreted in order to reconstruct the stratigraphic relationships between the Quaternary marine seismic units and the volcanic acoustic substratum. Seven seismic units are recognized based on the geological interpretation of seismic profiles and using seismo-stratigraphic criteria. The top of the lowest seismic unit was correlated with the Campanian Ignimbrite (Southern Campania Volcanic Zone). The stratigraphic setting of the eastern Bay of Naples is characterized by NE-SW trending seismic structures, probably corresponding with tuff rings. These tuff rings can be compared with the Porto Miseno, the Archiaverno and Averno, and the Astroni tuff rings (Campi Flegrei). Offshore, the Somma-Vesuvius a seismic unit was interpreted as the fallout deposits representing the base of the AD 79 eruption. However, since a branch of the isopach of 5 m of the “Pomici di Avellino” pyroclastic deposits is very close to the Tyrrhenian coastline and near our GRNA01 and GRNA03 seismic profiles, we cannot exclude that the seismic unit could be also correlated with the deposits of this eruption.

The Fucino 250–170 ka tephra record: New insights on peri-Tyrrhenian explosive volcanism, central mediterranean tephrochronology, and timing of the MIS 8-6 climate variability

The Fucino Basin, central Italy, with its long and continuous history of Quaternary sediment accumulation, is one of the richest Mediterranean Middle Pleistocene tephra records. Here, we present a new detailed investigation of tephra layers (tephras) of the 250–170 thousand years before present (ka) interval, corresponding to the entire Marine Isotope Stage (MIS) 7 and parts of the MIS 8 and MIS 6. The investigated tephras have been characterised in terms of major, minor and trace elements, Sr–Nd isotopic compositions and 40Ar/39Ar ages. For correlation purposes, glass compositions and several new 40Ar/39Ar ages of selected proximal pyroclastic units spanning the same temporal interval from Vulsini (Latera Volcanic Complex), Sabatini, and Vico volcanic systems, central Italy, were measured. The late MIS 8-early MIS 6 Fucino tephras were backtracked to their corresponding volcanic sources, which include the Vulsini, Vico, Sabatini, Roccamonfina, Ischia and Campi Flegrei volcanic systems. While some of these tephras have been correlated to specific eruption units, other layers are currently not documented or described in near-vent sections, thus highlighting previously unrecognised events generated by these volcanic systems. Furthermore, the new high precision 40Ar/39Ar ages provide improved temporal constraints for Fucino making it one of the most detailed and chronologically best constrained tephra records for central Mediterranean MIS 7 tephrochronology. The Fucino record thus provides new integrative information for reconstructing the explosive history of Italian volcanoes during the investigated time interval. Furthermore, the geochronological constraints provide the basis for future paleoclimatic investigations at local and regional scale.

On the feasibility and usefulness of high-performance computing in probabilistic volcanic hazard assessment: An application to tephra hazard from Campi Flegrei

For active volcanoes, knowledge about probabilities of eruption and impacted areas becomes valuable information for decision-makers to develop short- and long-term emergency plans, for which probabilistic volcanic hazard assessment (PVHA) is needed. High-resolution or spatially extended PVHA requires extreme-scale high-performance computing systems. Within the framework of ChEESE (Center of Excellence for Exascale in Solid Earth; www.cheese-coe.eu), an effort was made to generate exascale-suitable codes and workflows to collect and process in some hours the large amount of data that a quality PVHA requires. To this end, we created an optimized HPC-based workflow coined PVHA_HPC-WF to develop PVHA for a volcano. This tool uses the Bayesian event tree methodology to calculate eruption probabilities, vent-opening location(s), and eruptive source parameters (ESPs) based on volcano history, monitoring system data, and meteorological conditions. Then, the tool interacts with the chosen hazard model, performing a simulation for each ESP set or volcanic scenario (VS). Finally, the resulting information is processed by proof-of-concept-subjected high-performance data analytics (HPDA) scripts, producing the hazard maps which describe the probability over time of exceeding critical thresholds at each location in the investigated geographical domain. Although PVHA_HPC-WF can be adapted to other hazards, we focus here on tephra (i.e., lapilli and ash) transport and deposition. As an application, we performed PVHA for Campi Flegrei (CF), Italy, an active volcano located in one of the most densely inhabited areas in Europe and under busy air traffic routes. CF is currently in unrest, classified as being in an attention level by the Italian Civil Protection. We consider an approximate 2,000 × 2,000 × 40 km computational domain with 2 km grid resolution in the horizontal and 40 vertical levels, centered in CF. To explore the natural variability and uncertainty of the eruptive conditions, we consider a large number of VSs allowing us to include those of low probability but high impact, and simulations of tephra dispersal are performed for each of them using the FALL3D model. Results show the potential of HPC to timely execute a vast range of simulations of complex numerical models in large high-resolution computational domains and analyze great volumes of data to obtain quality hazard maps.

Fault systems in the offshore sector of the Campi Flegrei caldera (southern Italy): Implications for nested caldera structure, resurgent dome, and volcano-tectonic evolution

The structure of a caldera may influence its activity, making its understanding crucial for hazard assessment. Here, we analysed high-resolution seismic profiles in the Campi Flegrei (southern Italy) offshore sector. We recognised two main fault systems, including those associated with the formation of the caldera and those affecting the resurgent dome. The former system comprises three broadly concentric fault zones (inner, medial and outer ring fault zones) depicting a nested caldera geometry. Considering the relations between faults and seismic units that represent the marine and volcaniclastic successions filling the caldera, all ring faults were formed during the Campanian Ignimbrite eruption (40 ka) and subsequently reactivated during the Neapolitan Yellow Tuff eruption (15 ka). In this last caldera-forming event, the inner and medial fault zones accommodated most of the collapse and were episodically reactivated during the younger volcano-tectonic activity. The second fault system occurs in the apical zone of the resurgent dome and comprises dominantly high-angle normal faults that are mainly related to the volcano-tectonic collapse that followed the Agnano-Monte Spina Plinian eruption (4.55 ka). Finally, we provide a volcano-tectonic evolutionary model of the last 40 kyr, considering the interplay among ring and dome faults activity, volcaniclastic sedimentation, ground deformation and sea-level changes.

Raman focal point on Roman Egyptian blue elucidates disordered cuprorivaite, green glass phase and trace compounds

The discussed comparative analyses of Roman Imperial pigment balls and fragmentary murals unearthed in the ancient cities of Aventicum and Augusta Raurica (Switzerland) by means of Raman microspectroscopy pertain to a predecessor study on trace compounds in Early Medieval Egyptian blue (St. Peter, Gratsch, South Tyrol, Northern Italy). The plethora of newly detected associated minerals of the raw materials surviving the synthesis procedure validate the use of quartz sand matching the composition of sediments transported by the Volturno river into the Gulf of Gaeta (Campania, Southern Italy) with a roasted sulphidic copper ore and a mixed-alkaline plant ash as fluxing agent. Thus, the results corroborate a monopolised pigment production site located in the northern Phlegrean Fields persisting over the first centuries A.D., this in line with statements of the antique Roman writers Vitruvius and Pliny the Elder and recent archaeological evidences. Beyond that, Raman spectra reveal through gradual peak shifts and changes of band width locally divergent process conditions and compositional inhomogeneities provoking crystal lattice disorder in the chromophoric cuprorivaite as well as the formation of a copper-bearing green glass phase, the latter probably in dependency of the concentration of alkali flux, notwithstanding that otherwise solid-state reactions predominate the synthesis.