1. Introduction
More than one hundred years of observations of active volcanoes show that pre-eruptive stages are usually accompanied by seismicity, geochemical variations, and ground deformation—i.e., unrest. Unrest can also end without an eruption; the possibility of having eruptive or non-eruptive unrest may depend on the response to the passage of seismic waves and deformation induced by earthquakes and/or the activity of adjacent volcanoes, which may interact through stress variations caused by pressure and/or permeability changes and/or magmatic connections, e.g., because of dykes and/or a common magmatic system at depth and/or direct fluid transfer (e.g., [
1]).
Short- and long-term interactions among neighbouring volcanoes have only recently been considered, with some examples reported in the literature. In this section, we summarise three well-known examples of short- and long-term volcanic interactions; then, we briefly describe Campi Flegrei and Vesuvio volcanoes; lastly, we report the main results of this work.
On a short timescale, temporally correlated activity occurred at Mauna Loa and Kilauea volcanoes (Hawaii) in 2002. This correlation has been ascribed to a shared pulse of magma or, alternatively, the arrival of a pulse of magma that entered Mauna Loa’s plumbing system. According to the latter ascription, the magma pulse could have slightly perturbed Kilauea’s already-overpressurised adjacent shallow magma system, increasing its pressure and triggering the effusive event [
2]. However, Mauna Loa and Kilauea had anti-correlated activities over the past century, possibly due to competition for the magma supply [
3].
Both correlated and anti-correlated volcanic activities were found in Japan at Aira—associated with Sakurajima—and Kirishima—about 20 km from Aira—between 2009 and 2013, around the 2011 Kirishima eruption [
4]. Both volcanoes inflated before the lava extrusion phase; Aira deflated during the eruption and stopped deflating almost concomitantly with the Kirishima replenishment phase. Interaction between the two volcanoes has been ascribed to a hydraulic connection between their plumbing systems, fed by a common deep reservoir located in the mantle [
4].
In the longer term, anti-correlated eruptive histories occurred at the nearby Bolsena and Latera calderas (Italy). Much of Bolsena’s subsidence broadly coincided with activity at Latera (west rim of Bolsena); moreover, the similar compositions of the erupted magmas suggest the existence of a common magmatic reservoir [
5].
Despite this encouraging evidence, reported volcano interactions are still limited, and a general framework for understanding these processes is lacking. However, as these interactions may determine the eruption of at least one of the volcanoes, their definition and interpretation may affect eruption forecasting. The knowledge of volcano interactions may significantly benefit from additional observations, such as those reported in this work related to Campi Flegrei and Vesuvio, Italy. Campi Flegrei and Vesuvio are an iconic pair of active and high-risk volcanoes. They are located approximately 25 km apart inside a highly populated metropolitan area, which includes the city of Naples (
Figure 1a,b). Thus, it is of great importance to identify signs of any possible interaction between these two volcanoes and recognise indications of changes in hazards associated with even one of them.
Campi Flegrei’s and Vesuvio’s volcanic histories are intriguing, complex, and much studied, and there is a large amount of literature on the anatomies of the two volcanic systems. To frame the analyses performed in this work and our results, it is essential to summarise, albeit partially, the main characteristics of the behaviours of Campi Flegrei and Vesuvio during recent decades and the related available data, paying special attention to deformation data. This is especially important as regards Vesuvio, since Campi Flegrei caldera has recently received much more attention because of its ongoing uplift (e.g., [
6] and references therein).
In brief, Campi Flegrei, west of Naples, consists of a caldera formed before 15 ka (e.g., [
7,
8]). Several tens of eruptions (at least 60) and the resurgence of the caldera centre— i.e., long-term floor uplift—have occurred since then; the caldera-wide uplift has been largest near Rione Terra (R in
Figure 1a). The last eruption occurred in 1538 AD, after about 3 ka of quiescence, with more than 1 ka of subsidence since Roman times, and after a period of increasing seismicity and uplift from about 1400 to 1536; the eruption was followed by deflation [
9]. Ground deformation episodes that occurred during the last century have been measured by precise levelling surveys initiated in 1905 [
10]. More recently, 3-D ground deformation has been obtained by continuous Global Positioning System (cGPS) observations provided by the Neapolitan Volcanoes Continuous GPS (NeVoCGPS) network; the oldest stations started operating in 2000 [
11]. Since the 1990s, differential synthetic aperture radar interferometry (DInSar) data from ERS/Envisat and other missions’ imagery have also been available. Four main unrest episodes occurred between 1950 and 1952, between 1969 and 1972, between 1982 and 1984, and between 2005 and the present, with maximum uplifts of about 0.7 m, 1.7 m, 1.8 m, and more than 1 m at the time of writing, respectively (e.g., [
10,
12]); 1985–2002 was characterised by long-lasting decreasing-rate subsidence overlapped by a few short (several months) uplift phases—usually known as “mini uplifts”—also culminating in the caldera centre.
Figure 2 shows ground level changes at Rione Terra—located in the area of maximum vertical displacement and indicated by R in
Figure 1a—from 1980 to 2022.
Figure 1.
Map of Campi Flegrei and Vesuvio volcanoes, UTM WGS84 33N coordinates. (
a) Topography; colours indicate altitudes [
13]. Black dots and labels: L, LICO cGPS station, here used as a reference for comparing DInSAR and cGPS displacements; R, Rione Terra, in the area of maximum vertical movement of Campi Flegrei; S, Solfatara crater; P, Pisciarelli fumaroles; A, Agnano plain. (
b) Satellite image (map data: Google, ©2022 Landsat/Copernicus) showing the densely inhabited urban area of Naples. Cyan and yellow discontinuous lines with triangles indicate the borders of the Campanian Ignimbrite (CI, about 39 ka) and Neapolitan Yellow Tuff (NYT, about 15 ka) caldera collapses, respectively (adapted from [
14]). The white curve indicates the Mt Somma caldera rim. The red polygon encloses the area that is used here as reference for the analysis of DInSAR data.
Figure 1.
Map of Campi Flegrei and Vesuvio volcanoes, UTM WGS84 33N coordinates. (
a) Topography; colours indicate altitudes [
13]. Black dots and labels: L, LICO cGPS station, here used as a reference for comparing DInSAR and cGPS displacements; R, Rione Terra, in the area of maximum vertical movement of Campi Flegrei; S, Solfatara crater; P, Pisciarelli fumaroles; A, Agnano plain. (
b) Satellite image (map data: Google, ©2022 Landsat/Copernicus) showing the densely inhabited urban area of Naples. Cyan and yellow discontinuous lines with triangles indicate the borders of the Campanian Ignimbrite (CI, about 39 ka) and Neapolitan Yellow Tuff (NYT, about 15 ka) caldera collapses, respectively (adapted from [
14]). The white curve indicates the Mt Somma caldera rim. The red polygon encloses the area that is used here as reference for the analysis of DInSAR data.
There is wide agreement that the 1982–1984 unrest episode resulted from the emplacement of a sill-like magmatic intrusion at about 3.5 km depth offshore Rione Terra, although the hydrothermal system may have also played an important role (e.g., [
15,
16,
17,
18,
19,
20]). The causative process (intrusion of magma and/or magmatic fluids or instability of the hydrothermal system) for post-1985 unrest episodes is still debated (e.g., [
20,
21]). The dynamics of the sill-like source centred offshore Rione Terra is consistent with the majority of the caldera-scale deformation during both inflation and deflation periods since 1980; almost all the residual deformation is confined to the area of the Solfatara fumarolic field (S in
Figure 1a; e.g., [
6,
15]) and connected to the very shallow hydrothermal system. The feeding system of magma and/or magmatic fluids consists of a main storage reservoir at about 8–9 km depth containing 80–90% melt, interacting with shallow (mainly 2–4 km deep) reservoirs and recharged by a deeper source (e.g., [
22,
23,
24,
25]).
Ground deformation at Campi Flegrei has been—and still is—accompanied by seismicity. As for the main unrest episode that occurred between 1982 and 1984, seismicity started in 1982, reached significant magnitudes in 1983, was very intense during 1984, and waned at the end of the same year. Epicentres were scattered within the whole caldera, although two clusters could be identified, with the major one in an area between Rione Terra and the Agnano plain (A in
Figure 1a) and the minor one at sea; the hypocentre depth distribution exhibits a maximum between 1 and 3 km, with only a few events occurring down to 7 km [
26]. Few events and swarms occurred between 1985 and 2000, coevally with mini-uplifts. After 2000, the seismic rate slightly increased; however, earthquakes were still rare and occurred in low-magnitude swarms until 2014, when the events became more frequent and were often higher in magnitude (e.g., [
26,
27]). Until 2011, the earthquake depth distribution was similar to that between 1982 and 1984 [
26], but with no significant events deeper than 4 km. After 2014, seismicity became shallower (usually < 2 km), except for the main swarms [
27]. The onset of microseismicity is, at least in part, ascribed to hot fluid injections from depths (e.g., [
19]), regulated by a rigid caprock [
28] that acts as a barrier to magmatic fluids coming from the shallow reservoir and/or the deep magma chamber.
Vesuvio (
Figure 1)—or, more precisely, Somma-Vesuvio—is a stratovolcano with a summit caldera (Mt Somma) and a recent cone (Vesuvio) resulting from several Plinian eruptions, with the last one occurring in 79 AD (e.g., [
29]). Sub-Plinian eruptions occurred in 472 and 1631 AD. The latter was followed by semi-persistent activity that lasted until 1944; the quiescence duration correlates with the eruption size at the beginning of an eruptive period [
30]. An approximately 4 km deep magma reservoir fed eruptions in the last 2000 years (e.g., [
31,
32]), although there is currently no evidence of stored melt [
33].
Levelling surveys at Somma-Vesuvio were initiated at the beginning of the 1900s, but unfortunately, most of the few benchmarks of the early network were lost. The levelling network was subsequently expanded; however, until 1985, it still consisted of a few benchmarks placed in a poor configuration, which prevented geodetic surveyors from performing a rigorous adjustment process and robust quantitative analyses [
34]. The analysis of 1986–2009 data shows strong subsidence—up to −1.5 cm/year—inside the caldera rim, probably because of the sliding and compaction of young incoherent materials [
34].
As for DInSAR data, we are not aware of any studies comparing deformation patterns related to different time periods. Published works discuss the mean ground displacement rates (velocities) during the periods 1992–2000 [
35], 1992–2003 [
36], and 1992–2010 [
37]. The approximate vertical displacement (see
Section 2.1) time series from 2003 to 2013—two measurement points inside the caldera rim and a third point around 5 km west of the Vesuvio summit—are shown in [
38]. DInSAR displacement rates confirmed subsidence within the caldera rim and showed subsidence spots at a discontinuous “girdle” about 10 km from the summit [
35,
36,
37]. The mean 1992–2003 velocities also show the relative subsidence of the E and SE sectors of the Vesuvian area relative to the NW and W sectors and an eastward component ascribed to large-scale extension [
36], while the mean 1992–2010 velocities show the subsidence of the S-SW flank as well as uplift W-NW and a few kilometres SE of the volcanic edifice [
37]. Velocities are on the order of a few mm/year, apart from inside the caldera rim. Subsidence spots were ascribed to the loading of volcanic and sedimentary units at the contact between volcanic and non-volcanic rocks [
35]. The overall deformation pattern of the Vesuvian area and the presence of an offshore arc-shaped ridge, formed by the folding and cutting of recent strata by minor thrust faults, were ascribed to the volcano’s gravitational spreading on its ductile sedimentary substratum [
36]. However, it was recently suggested that the deformation of Somma-Vesuvio is also partly caused by sagging—i.e., the flexure of a brittle layer into an underlying low-viscosity ductile layer—based on the mean 1992–2010 ground velocity, analogue experiments, and numerical modelling [
37].
Since 2001, 3-D ground deformation has also been obtained by cGPS data provided by the NeVoCGPS network—nine stations, three of which started operating in 2001; station displacements are given with respect to six stations located outside the Neapolitan volcanic area [
11]. cGPS observations confirm the strong subsidence inside the caldera rim and indicate the general subsidence of the Vesuvian area with no clear pattern; horizontal velocities are not consistent with the spreading of the volcano edifice (see also [
39]).
Seismicity has been recorded with modern instrumentation since the beginning of the 1970s; the network has been progressively improved in the subsequent decades. Since the last eruption in 1944, intense seismicity was recorded during the periods 1978–1980, 1989–1990, 1995–1996, and 2018; a major seismic crisis occurred between 1999 and 2000, lasting until the beginning of 2002 (e.g., [
27,
40,
41,
42,
43]). Seismicity is concentrated close to the Vesuvio crater axis, within a radial distance of less then 2 km. Two main seismogenic volumes have been identified: one with hypocentres located at depth below the crater and clustered mostly between 1 and 6 km below sea level, and another one with hypocentres clustered above sea level and located within the volcanic edifice. Deep magmatic fluids were probably involved in at least the 1999–2000 crisis (e.g., [
27] and references therein). Deep microseismicity and magmatic fluid injections are most likely due to the dynamics of a large sill-like magma chamber at about 8 km depth—presently containing about 10% melt [
44,
45,
46]—which also supplied magma to the approximately 4 km deep magma reservoir feeding eruptions over the past two millennia.
Since 2002, deep seismicity significantly decreased, and its absolute minimum occurred in 2005 (e.g., [
43,
47,
48]). From 2008, after a slight increase, deep seismicity continued to decrease, whereas shallow earthquakes were increasingly recorded [
27]. The time–depth distribution from 1999 to 2012 shows that seismic strain was mostly released by the deeper volume until 2002 and within the shallower volume since then (e.g., [
43]). More recent seismicity consists of many hundreds of low-magnitude, mostly shallow earthquakes per year, mainly located in the crater area and sometimes occurring as low-energy swarms; since 2018, an increase in the number of earthquakes per year has been recorded (e.g., [
27,
48]).
Magmatic volatiles are actively transported by groundwater flowing along the main faults and fractures. Groundwater samples collected in the southern and western sectors and in the eastern and northern sectors of the volcano are chemically different: in the southern and western sectors, the contributions of deep fluids are more evident, while in the eastern and northern ones, shallow meteoric fluids dominate (e.g., [
47]). A fumarolic field covers the rim and the inner flanks of the Vesuvio crater, and diffuse soil CO
2 emissions come from some sectors of the cone [
49]. The temperature of the groundwater and fumarolic field emissions suggest that the activity of the hydrothermal system decreased from 1998 to 2010, with the exception of a few periods [
47]. The CO
2 soil flux, equilibrium P
, and CO
2/CH
4 fumarolic ratio culminated in 2002 at the end of the 1999–2001 seismic crisis, indicating a coeval increase in the relative amount of magmatic fluids and in the Vesuvio activity (e.g., [
49]); other changes have been ascribed to the effects of permeability variations accompanying seismogenic (mainly tectonic) processes (e.g., [
48]).
The existence of possible links between Campi Flegrei and Vesuvio is suggested by the presence of a wide low-velocity layer—interpreted as an extended magmatic body recharged by a deeper source—at about 8 km depth beneath Campi Flegrei and Vesuvio (e.g., [
23,
44]) and the similarity in the petrological features of the deposits from the two volcanoes (e.g., [
50,
51,
52]). The analysis of 1992–2008 DInSAR time series related to Rione Terra (
Figure 1a) and the Vesuvio summit cone area suggested correlated activity during the 2001 and 2006 mini-uplifts, superimposed on prominent long-term subsidence [
53]. However, uplift subsequently affecting Campi Flegrei was not accompanied by correlated changes at Vesuvio, and synchronicity appears to be in opposition to the historical records showing that pronounced subsidence periods at Campi Flegrei are associated with eruptions of Vesuvio and vice versa [
53].
Campi Flegrei’s and Vesuvio’s deformation histories in recent decades are of great interest; the turn of the 2000s was quite important for the Campi Flegrei dynamics, so in this work, we searched for changes in Vesuvio’s behaviour between the late 1990s and early 2000s. We analysed ground displacement in the Vesuvian area and its surroundings during the period 1993–2010 using ERS/ENVISAT imagery. Empirical Orthogonal Function (EOF) analysis revealed a sudden large-scale displacement rate change around 2001; the pattern of the second EOF mode and the differences between pre-2001 and post-2002 velocity maps are fully consistent with the deflation of a deep pressurised source. The coincidence between the transition from deflation to inflation at Campi Flegrei and the onset of deflation at Vesuvio may suggest the possible transfer of magma and/or magmatic fluids between the two plumbing systems. Our results help define Vesuvio’s present deformation style and may be of primary importance in identifying and properly assessing any possible reactivation processes.
3. Results
Figure 7a,b show maps of the LOS (ascending and descending orbits) and approximate upward and eastward mean displacement rates (velocities) in the Vesuvian area for the period 1993–2010. The LOS velocity maps are similar to
Figure 3 in [
37], apart from a uniform velocity translation due to a different reference area. As already pointed out by others, the most evident features are the occurrence of subsidence at the summit and S-SW flank of the volcanic edifice [
37], subsidence spots along a discontinuous “girdle” about 10 km from the Vesuvio summit, and the occurrence of uplift W-NW and S-SE of the volcanic edifice [
35,
36,
37].
A closer look at the two most prominent subsidence spots is provided in
Figure 8. The subsidence rate decreased over time (
Figure 8b); intriguingly, the local urbanisation rate—here defined as the number of new buildings per year in the municipalities where the spots are located—also decreased from the early 1990s compared to previous decades (
Figure 8c). In this respect, a third deformation spot—which is located about 5 km north of the “girdle”—is particularly instructive: the subsidence rate is not constant (
Figure 8b), and satellite images acquired in December 1985 (
Figure 8c) and August 2002 (
Figure 8d) show the construction of an industrial area during that time interval.
To identify possible common modes in the ground displacement field, we applied an EOF analysis to the 1992–2010 DInSAR displacement time series inside two circular regions 15 km and 9 km in radius, centred on the magenta star in
Figure 7; the radii of the two circular regions were chosen in order to exclude the westernmost part of the region shown in
Figure 7, where Campi Flegrei dynamics may contribute to the horizontal movements, and include (15 km radius) or exclude (9 km radius) the “girdle” of discontinuous strong subsidence. In both regions, the data are arranged into two different matrices:
, containing ascending- and descending-orbit LOS displacements, and
, containing approximate upward and eastward displacements. Thus, the total number of EOF analyses is four. The ECs and eigenvalues of the first modes are shown in
Figure 9. There is practically no difference between the results obtained by analysing
and
: only the first two modes are statistically significant; the EC of the first mode follows a quasi-linear trend, and the EC of the second mode exhibits a sudden slope change between pre-2001 and post-2002 data.
Figure 10a–d show the spatial patterns of the second EOF mode when the radius of the analysed circular region is 9 km. When the radius of the analysed circular region is 15 km, the spatial patterns of the second EOF mode—not shown here for brevity—are very similar to those in
Figure 10a–d inside the first 9 km from the centre. Bluish colours imply a mean velocity change from positive to negative values around 2002, because the displacement time series related to a selected EOF mode is given by the EOF coefficient at the pixel times the mode EC.
The spatial patterns of the second EOF mode closely resemble the displacement patterns for a deep (several kilometres) deflating pressurised cavity (
Figure 4a,b,e,f). We have inverted the EOF coefficients in
Figure 10a,b (ascending- and descending-orbit LOS displacements) for the best-fit parameters of a very small (point) pressurised spheroidal cavity embedded in a homogeneous elastic half-space (
Section 2.3). Pixels around the Vesuvio cone have been excluded because of the local sliding and compaction of young incoherent materials [
34].
Figure 10e,f show modelled ascending- and descending-orbit LOS displacements when the source depth is 7300 m and the polar-to-equatorial radius is 1.1. The position of the source surface projection and polar-to-equatorial radius are well constrained by the inversion; the depth uncertainty is a few hundred metres.
Figure 9 and
Figure 10 indicate that the years before 2001 and those after 2002 may be characterised by different ground dynamics. To test this hypothesis, we simply computed and compared pre-2001 and post-2002 pixel-wise mean velocities. Before 2001 (
Figure 11a,c and
Figure 12a,c), no noticeable pattern is apparent in the Vesuvian area as in the approximate upward velocity (
Figure 12a), apart from strong subsidence close to the summital crater and in various spots along a quasi-circular ring around Vesuvio—i.e., the “girdle” in [
35,
36]. The approximate eastward velocity is characterised by westward motion west of Vesuvio and eastward motion east of Vesuvio (
Figure 12c); the two side lobes might be related to the spreading of the volcanic edifice on its ductile sedimentary substratum, although the west lobe extends farther than expected to the west, possibly because of the strong subsidence of the Campi Flegrei caldera. We do not ascribe the eastward velocity pattern to regional extension [
36] just because of the lobed shape.
Post-2002 velocity maps (
Figure 11b,d and
Figure 12b,d) show very different features, which are easily interpretable when looking at the approximate upward and eastward velocity components. After 2002, the strong subsidence spots are less evident than before 2001, and the whole Vesuvian area was clearly subsiding (
Figure 12b); the two side lobes in
Figure 12d are consistent with the deflation of a pressurised cavity. Post-2002, upward and horizontal cGPS velocities are also consistent with the deflation of a pressurised cavity (
Figure 12b,d).
The difference between post-2002 and pre-2001 velocities is shown in
Figure 13. The “girdle” subsidence spots are characterised by positive values (
Figure 13c), thus confirming the subsidence rate decrease (
Figure 8). Maps related to the approximate upward and eastward velocities (
Figure 13c,d) are consistent with subsidence due to deep deflation (see
Figure 4), although the west lobe extends farther than expected to the west, as in the case before 2001 and probably for the same reasons.
4. Discussion
In this work, we used and analysed multi-year DInSAR ground displacement time series related to the Vesuvian area. The SAR phase may be affected by tropospheric effects; however, while seasonal oscillations were successfully removed from the time series of DInSAR displacements from ENVISAT/RADARSAT imagery at Vesuvio by utilising surface weather and radiosonde data, the trends were practically unaffected by the atmospheric correction [
38].
The second EOF mode shows a clear change between 2001 and 2002; its spatial pattern is consistent with a quasi-spherical deflating cavity embedded in a homogeneous elastic half-space at about 7300 m depth beneath the Vesuvio summit (
Figure 10a,b,e,f). The actual source depth is probably larger because of the “focusing effect” related to shallow soft materials [
15].
The main features of the pre-2001 ground velocity field (
Figure 11a,c and
Figure 12a,c) were already pointed out in [
35,
36]: strong subsidence occurs inside the caldera ring close to the summital crater and in various spots along a quasi-circular “girdle” at about 10 km around Vesuvio; the approximate eastward velocity is consistent with the spreading of the volcanic edifice, with some possible contributions from regional extensional tectonics and/or the strong Campi Flegrei subsidence. The summital crater subsidence is probably caused by the sliding and compaction of young incoherent materials [
34]. As for the subsidence spots, a possible link to water table changes has been rejected, and they have been ascribed to the loading of volcanic and sedimentary units at the contact between volcanic and non-volcanic rocks [
35]. However, the spot subsidence rate decreased over time (
Figure 8b), and spots are located in areas where faster and faster urbanisation occurred from the 1960s to the 1990s, when the rate of the construction of new buildings rapidly decreased (
Figure 8c). Thus, we believe that the spot subsidence might be ascribed to the loading of buildings at the contact between volcanic and non-volcanic rocks. Groundwater could actually contribute to the differences between the pre-2001 upward velocities in the southern and western sectors in comparison with the eastern and northern sectors (
Figure 12c), because groundwater circulation differs in the two areas [
47].
If the pre-2001 ground velocity field is mainly due to the spreading of the volcanic edifice onto its sedimentary substratum and/or to the regional tectonics, we are thinking of long-lasting processes—thousands of years. Thus, the post-2002 ground velocity field, which differs substantially from the pre-2001 one, probably includes the previously acting processes and a new phenomenon that began around 2001–2002. Taking the difference between the post-2002 and pre-2001 ground velocity fields could highlight the characteristics of this new phenomenon (
Figure 13), except for the areas of decreasing local subsidence induced by human activity (reddish areas in
Figure 13c) and the possible transient contribution of the Campi Flegrei dynamics to the horizontal movements in the westernmost region. The ground deformation around Vesuvio is consistent with the deflation of a several-kilometres-deep pressurised cavity (
Figure 4) but also with volcano sagging into an underlying low-viscosity ductile layer, to which part of the 1992–2010 ground velocity field was attributed by [
37]. Analogue experiments show a possible transition from spreading-dominated to sagging-dominated behaviour after cone emplacement [
39]. However, sagging is also a long-lasting process, and the sudden onset of the post-2002 phenomenon makes sagging unlikely.
The hypothesis of the deflation of a several-kilometres-deep pressurised cavity is also supported by the significant decrease in Vesuvio seismicity at depths larger than about 3 km since around 2002 [
27,
43]. Deep magmatic fluids were probably involved in at least the 1999–2000 seismic crisis (e.g., [
27] and references therein). A coeval increase in the ascent of magmatic fluids is suggested by the time series of CO
2 soil flux, equilibrium P
, and the CO
2/CH
4 fumarolic ratio [
49]. In addition, the groundwater temperature and fumarolic emissions suggest that the activity of the hydrothermal system generally decreased from 1998 to 2010 [
47,
49]. As already mentioned in
Section 1, deep microseismicity and magmatic fluid injections are most likely due to the dynamics of a large sill-like magma reservoir at about 8 km depth, presently containing about 10% melt [
44,
45,
46], which also supplied magma to the approximately 4 km deep magma reservoir feeding eruptions over the past two millennia.
The coeval onset of deflation below Vesuvio (this work) and inflation below Campi Flegrei [
6] seems consistent with the observation that historically pronounced subsidence periods at Campi Flegrei are associated with eruptions of Vesuvio and vice versa [
53]. However, 1992–2008 DInSAR time series show two mini-uplifts at the Vesuvio summital area in 2001 and 2006 [
53], i.e., when two mini-uplifts were also recorded at Campi Flegrei (
Figure 2). This apparent inconsistency—previously identified in [
53]—could be linked to the somewhat different nature of the synchronous mini-uplifts compared to the post-2002 uplift of Campi Flegrei and the subsidence of the Vesuvian area. The ground deformation pattern of the 2001 and 2006 Campi Flegrei mini-uplifts was the same as that of the post-1985 subsidence [
15], and both the subsidence and the mini-uplifts were due to an approximately 3.5 km deep source. The coeval Vesuvio mini-uplifts affected the summital area and were probably caused by a shallow—around 2 km deep—source [
53]. Thus, the Campi Flegrei and Vesuvio deep reservoirs were not involved in the mini-uplifts. On the contrary, the ongoing Campi Flegrei uplift is probably driven by the inflation of a 8 km deep reservoir [
6]. As for the post-2002 deflation below Vesuvio, a several-kilometres-deep deformation source is suggested by the pattern of the second EOF mode, as well as the large extension of the affected area and the apparent southward shift of the approximate upward velocity change pattern compared to the approximate eastward velocity change pattern (see
Figure 4a,b,e,f and
Figure 13c,d). In view of the above, an anti-correlation between the deformation signals of Campi Flegrei and Vesuvio seems to exist only when the observed deformation in both volcanoes is connected to the deep reservoirs, suggesting a connection/interaction between the deep reservoirs of Campi Flegrei and Vesuvio.
Unfortunately, in the Vesuvian area, the deformation rate is small (a few mm/yr), the central part of the deforming area is poorly covered by data, and three-dimensional rheological parameters are poorly known, particularly those at depth. A reliable inversion of the pattern of the second EOF mode and/or the difference between post-2002 and pre-2001 ground velocities for the best-fit deformation source is out of the scope of this work. Our simple modelling of the second EOF mode for ascending- and descending-orbit LOS displacements suggest depressurisation deeper than 7000 m. Ground displacement time series from GPS do not show clear changes in velocity over time [
11], suggesting that deflation below Vesuvio continues. In the near future, we aim to extend our analysis to DInSAR ground displacement times series after 2010 using Sentinel-1 and—if made available—TerraSar-X imagery. Such an extension could improve the signal-to-noise ratio and enlarge the area covered by data, hopefully clarifying whether there has been an increase in the Vesuvio deflation rate corresponding to the Campi Flegrei uplift acceleration since 2010.
5. Conclusions
The active Campi Flegrei and Vesuvio volcanoes are located about 25 km apart inside a highly populated metropolitan area, which includes the city of Naples; the most recent eruptions occurred in 1538 and 1944, respectively. After a large uplift, Campi Flegrei was generally subsiding from 1985 to 2001, while it has been uplifting at an accelerating rate since then, probably driven by deep magma inflation. To search for possible interactions between Campi Flegrei and Vesuvio, we analysed Vesuvio dynamics during the period 1992–2010 using DInSAR data from ERS-ENVISAT imagery.
Previous works computed and analysed the mean ground velocity field in the Vesuvian area during the periods 1992–2000 [
35], 1992–2003 [
36], and 1992–2010 [
37], without looking for temporal variations. Here, we show that the second mode of an Empirical Orthogonal Function analysis of the displacement time series indicates a sudden change in trend—i.e., mean velocity—between pre-2001 and post-2002 data. Its spatial pattern is in agreement with what is expected from a quasi-spherical deflating cavity embedded in a homogeneous elastic half-space at approximately 7300 m depth beneath the Vesuvio summit. The actual source depth is probably larger because of the “focusing effect” related to shallow soft materials [
15].
The pre-2001 ground velocity field confirms strong subsidence—up to 5 mm/y—occurring inside the caldera ring close to the summital crater and in various spots along a quasi-circular “girdle” at about 10 km around Vesuvio. The subsidence rate in the spots decreased over time between 1992 and 2010. The approximate eastward velocity is characterised by eastward displacement in a lobe east of Vesuvio and by westward displacement in a lobe west of Vesuvio; this is consistent with the spreading of the volcanic edifice on its ductile sedimentary substratum and/or extensional tectonics—as suggested by [
36] based on 1992–2003 mean ground velocities—with some possible contribution in the westernmost part from the strong Campi Flegrei subsidence. In contrast to prior works, we ascribe the subsidence spots to previous rapid urbanisation, not a load of volcanic and sedimentary units at the contact between volcanic and non-volcanic rocks—as suggested by [
35] based on 1992–2000 mean ground velocities—because the subsidence rate decreased over time between 1992 and 2010; moreover, we do not ascribe the eastward velocity field to regional extension—as suggested by [
36]—because of the lobed pattern.
Post-2002 velocity maps show general subsidence in the whole Vesuvian area, westward displacement in a lobe east of Vesuvio, and eastward displacement in a lobe west of Vesuvio. Thus, it probably includes both spreading and a new phenomenon that began around 2001–2002. The second EOF mode and the difference between post-2002 and pre-2001 ground velocities demonstrate a ground deformation field that is fully consistent with the deflation of a several-kilometres-deep pressurised cavity. Deep deflation is also in agreement with Vesuvio seismic and geochemical data. We do not invoke volcano sagging into a low-viscosity underlying ductile layer—as suggested by [
37] based on mean 1992–2010 ground velocities—because of the sudden onset of post-2002 subsidence.
The coincidence between the transition from deflation to inflation at Campi Flegrei and the onset of deflation at Vesuvio may suggest the possible transfer of magma and/or magmatic fluids between the two plumbing systems. What drives the inflation of CF and how does that relate to the coincidental deflation of Vesuvio in terms of fluid mechanics? Maybe a degassing–depressurisation mechanism (e.g., [
64]) is involved; however, we think that answering this question requires more data, not only geodetic: data from experimental petrology, geochemistry, seismology, etc., covering a long period of time, would also be needed. In the near future, we aim to extend our analysis to ground displacement times series after 2010 using Sentinel-1 and, if made available, TerraSar-X imagery, obtain more reliable characteristics of the deflating source, and hopefully clarify whether there has been an increase in the Vesuvio deflation rate corresponding to the Campi Flegrei uplift acceleration occurring since 2010.