Modulation of Ground Deformation and Earthquakes by Rainfall at Vesuvius and Campi Flegrei (Italy)

Volcanoes are complex systems whose dynamics is the result of the interplay between endogenous and exogenous processes. External forcing on volcanic activity by seasonal hydrological variations can influence the evolution of a volcanic system; yet the underlying mechanisms remain poorly understood. In the present study, we analyse ground tilt, seismicity rates and rainfall amount recorded over 6 years (2015–2021) at Vesuvius and Campi Flegrei, two volcanic areas located in the south of Italy. The results indicate that at both volcanoes the ground deformation reflects the seasonality of the hydrological cycles, whereas seismicity shows a seasonal pattern only at Campi Flegrei. A correlation analysis on shorter time scales (days) indicates that at Vesuvius rain and ground tilt are poorly correlated, whereas rain and earthquakes are almost uncorrelated. Instead, at Campi Flegrei precipitations can affect not only ground deformation but also earthquake rate, through the combined action of water loading and diffusion processes in a fractured medium, likely fostered by the interaction with the shallow hydrothermal fluids. Our observations indicate a different behavior between the two volcanic systems: at Vesuvius, rain-induced hydrological variations poorly affect the normal background activity. On the contrary, such variations play a role in modulating the dynamics of those metastable volcanoes with significant hydrothermal system experiencing unrest, like Campi Flegrei.

Introduction

In the recent years an increasing number of observations has been showing how rainfall-induced stress variations play a role in modulating ground deformation and seismic activity. Meteoric water storage and load can produce surface deformation, inducing crustal stresses and modulating seismicity (Craig et al., 2017). In addition, variations of the pore-fluid pressures induced by rainfall can promote both ground deformation (Kümpel et al., 2001; Dal Moro and Zadro, 1998; Westerhaus and Welle, 2002; Lesparre et al., 2017) and earthquake nucleation, often with seasonal patterns well correlated with precipitations (Muço, 1999; Heinzl et al., 2006).

The modulation of ground deformation and seismicity by rainfall assumes a critical role in volcanic environments because it superimposes on the endogenous processes and may influence the volcano dynamics. In particular, for those volcanoes in magmatic/hydrothermal unrest or prone to eruption, meteoric input can act as a trigger for volcanic activity. This interaction occurs on typical time scales going from hours, to days and months. As an example, the effect of meteoric water load and the subsequent interaction with groundwater at Nevado del Ruiz (Colombia) possibly caused deformation events and triggered phreatic explosions and seismicity occurred through 1985–1988 (Banks et al., 1990). The study of the seasonal pattern of eruptions over more than 150 years of available records at Piton de la Fournaise (La Réunion) suggested to explore the possible link with rainfall, by modelling the volumetric ground deformation in response to a water loading/unloading mechanism (Violette et al., 2001). At Soufrière Hills volcano (Montserrat, Lesser Antilles), the response of seismic activity to intense precipitations observed over 3 years (2001–2003) of monitoring has been attributed to both shallow and deeper fluid/solid interactions occurring at different time scales (from hours to days) and related to the percolation of water into cracks (Matthews et al., 2002, 2009). The authors conclude that the rainfall modulates existing, internal processes, rather than generating new events itself. Extreme rainfall is supposed to have triggered the 2018 rift eruption at Kīlauea volcano (Hawaii), by increasing the pore pressure at depths of 1–3 km and promoting the mechanical failure of the volcanic edifice (Farquharson and Amelung, 2020). It is noteworthy that the eruption was not caused by the forceful intrusion of new magma into the rift zone, but by a dyke propagation prompted by meteoric water infiltration. Also the statistical analysis of the occurrence of historical eruptions at Kīlauea strongly suggests a correlation between rainfall and volcanic activity (Farquharson and Amelung, 2020).

In the present paper we investigate and compare the possible links among rainfall, ground deformation and seismicity at two of volcanoes currently characterized by different activity: Vesuvius, which is in a quiescent state, and Campi Flegrei, a caldera experiencing unrest.

Vesuvius and Campi Flegrei are potentially dangerous volcanoes located in the area of Naples (Southern Italy). Vesuvius is a stratovolcano which experienced its last eruption in 1944 (Sbrana et al., 2020). It is characterized by a great variability in eruptive style: highly explosive, sub-Plinian and Plinian eruptions alternating with periods of Strombolian and effusive open-conduit activity, that typically followed long periods of quiescence (Cioni et al., 2008 and references therein). After the last subplinian eruption, in 1631, Vesuvius was marked by open-conduit activity, which ended with the eruption of 1944. The alternating activity of Vesuvius resulted in different deposits including lavas and pyroclastic products. The ground deformation inferred from more than 20 years of tiltmeter signals shows a common subsidence pattern of the southern part of the volcanic edifice, related to joint effects of gravitational sliding and extensional tectonic stress (Ricco et al., 2013, 2021). This background tilt showed gradual or abrupt interruptions in both trend and amplitude during phases of high seismic activity (Ricco et al., 2021), resulting in more complex deformation patterns and indicating a close link between ground tilt and seismicity. The most significant episodes related to variations of internal dynamics occurred in 2017–2019 (Ricco et al., 2021). Vesuvius seismicity consists mainly of volcano-tectonic (VT) earthquakes of low-to-moderate magnitudes (up to 3.6) and mainly located along the crater axis, up to a depth of 6 km below sea level (bsl) (Madonia et al., 2008). The source mechanisms are related to the interaction of the regional and local stress fields, as well as to fluid-driven rock fracturing triggered by pressure variations in the shallow hydrothermal system. Since 2003, low-frequency (LF) and long-period (LP) earthquakes have been observed and associated with mechanisms of brittle slow failure of dry rocks, and resonance of pre-existing fluid-filled cracks, respectively (Bianco et al., 2005; Petrosino et al., 2020a). At Vesuvius a correlation between rainfall and seismicity has not yet been observed; modulation of seismic activity on short (daily) time scales has been interpreted to be linked to a number of possible mechanisms related with the cooling/warming diurnal cycle of the volcanic edifice and/or with the daily oscillation of the geomagnetic field (Mazzarella and Scafetta, 2016). On long time scales (decades) it has been hypothesized that sea level changes could modulate the eruptive cycles (Bragato, 2015).

Campi Flegrei is a caldera generated by at least two major collapses, related to the Campanian Ignimbrite (∼39 ka), and the Neapolitan Yellow Tuff eruptions (∼15 ka). The caldera experiences slow subsidence alternating with fast ground uplift accompanied by seismicity. After the last eruption of Monte Nuovo in 1538, Campi Flegrei experienced a phase of subsidence, followed by episodes of uplift between 1950–1952 (about 0.75 m), 1969–1972 (1.77 m), and 1982–1984 (1.79 m). Since 1985, the area has undergone a further phase of subsidence, with only brief uplift events recorded in 1989, 1994, and 2000. A new episode of ground uplift started in 2005 and it is still ongoing with accelerated rate during 2011–2013 and subsequently since September 2020; to date the overall vertical deformation is about 0.8 m. Uplift episodes are likely related to the pressurization of hot fluids (CO₂ and H₂O) exsolving from a deep magmatic body towards a shallow hydrothermal system where they mix with meteoric components (Gresse et al., 2017). Hydrothermal fluids and gases are released at the surface of Solfatara crater and Pisciarelli fumarolic field, the most thermally active sites of the caldera (Caputo et al., 2020; Cusano et al., 2021). Uplifts are accompanied by seismicity consisting of both sequences of VT earthquakes and LP events. The most remarkable LP swarm occurred in October 2006 (Saccorotti et al., 2007; Falanga and Petrosino, 2012; De Lauro et al., 2012) and it was located at depths of about 500 m bsl beneath the Solfatara crater. The LP generation mechanism is related to the acoustic resonance of a crack filled by a water-gas mixture of hydrothermal origin (Cusano et al., 2008). The VT activity is characterized by the alternation of swarms lasting a few hours, and phases of low seismicity rate (Petrosino et al., 2018; Bellucci Sessa et al., 2021). Since 2018 the activity is increasing and the most relevant earthquakes (Md 3.1 and 3.3) occurred on December 6th, 2019 and April 26th, 2020. The earthquakes are mainly located beneath the Solfatara-Pozzuoli area, at depths up to 4 km bsl, with the highest concentration between 1–2 km bsl (Bellucci Sessa et al., 2021). The source mechanism is a brittle shear failure induced by the pressurization of the hydrothermal system and the consequent fluid flow towards the surface (Saccorotti et al., 2007). As shown in recent works, there is evidence that the volcanic system responds the to the rainfall both in term of ground deformation and seismic activity. In particular, diffusive processes into the highly fractured rocks, pore pressure and permeability variations have been advocated to explain the correlation between rain and ground tilt (Ricco et al., 2019; Petrosino et al., 2020b), as well as rain and seismicity (Petrosino et al., 2018, Scafetta and Mazzarella 2021).

In continuity with the previous studies, in the present work we perform a statistical correlation analysis of the time series of ground tilt, seismicity rates and rainfall amount recorded over 6 years (2015–2021). The results indicate that, although at both Vesuvius and Campi Flegrei ground deformation is sensitive to the input of meteoric water, the two volcanic systems respond in a different way to rain-induced hydrological variations. In particular, at a stressed caldera like Campi Flegrei the rainfall can modulate not only ground deformation on different time scales but also earthquake rate. These observations suggest that the Campi Flegrei dynamics is the result of a complex interplay between endogenous and exogenous processes.

Data and Methods

The borehole tiltmeter network at Vesuvius is currently composed of four instruments (TRC, IMB, CMG and CMT) which started to operate contemporaneously since September 2016 (Figure 1). At Campi Flegrei the network consists of three borehole tiltmeters (CMP, ECO and HDM) working since 2015 (Figure 1). Both networks are managed by Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Napoli–Osservatorio Vesuviano (INGV-OV). The tiltmeter stations consist of digital sensors model “Lily Self-Leveling Borehole Tiltmeter” equipped with a self-leveling bubble electrolyte with a dynamic range of ±330 μradians and a resolution less than 5 nradians. The tiltmeter package includes a magnetic compass and a temperature sensor. Ground tilt variations are measured along two orthogonal directions NS and EW; they are recorded with a sampling rate of 1 sample per minute. The depths of the bore hole installations range from 20 to 28 m (Ricco et al., 2018).

FIGURE 1

FIGURE 1. Map of Vesuvius (A) and Campi Flegrei (B) with the location of the borehole tiltmeters. The base maps were created by the INGV-OV Laboratory of Geomatics and Thermography, ^© LGT INGV OV, 2015.

We used tiltmeter data from October 1st 2016 to March 31st 2021 for Vesuvius, and from April 1st 2015 to March 31st 2021 for Campi Flegrei. The number of VT earthquakes for the same time intervals was obtained by the INGV-OV seismic catalogue, available at the url (https://www.ov.ingv.it/index.php/cataloghi-sismici-vulcani-napoletani?category[0]=73&category_children=1&own=0, last access 2021/04/30). Rainfall amount data sampled at 1 day for Pozzuoli and Torre del Greco towns were downloaded by the historical archive of the website https://www.3bmeteo.com/meteo/ora/storico (last access 2021/04/30).

Borehole tilt signals contain several components: trend (which is often related to the dynamics of the area and consists of net movement), seasonal and tidal deformation, and a residual which is related to isolated events (such as earthquakes), rainfall-induced deformation and short period (less than a few days) fluctuations (Kümpel et al., 2001; Garcia et al., 2010; De Lauro et al., 2018). The typical representation of tilt data (cumulative deformation over a certain time interval; Supplementary Figure S1) highlights the long period trend, minimizing short period noise which is often considered as a disturbance to discard. An equivalent way to represent data is by differencing the time series (Chatfield, 2013): in this case data will reveal short term variations, with the trend being removed. Thus, considering the time derivative of the tiltmeter signals is a useful approach to compare the ground deformation with rainfall amount (Breitenberger, 1999; Dal Moro and Zadro, 1998).

For our analyses, we used the differenced ground tilt and calculated the root mean square (RMS) according to:

$x_{rms}=\sqrt{\frac{1}{N}{\displaystyle \sum }_{i=1}^N{x}_i^2}$

where x represents the time series and N is the number of samples in the selected time window. After having tried several durations between 1 and 30 days, we chose a 1-day time interval, long enough to stabilize fast ground tilt fluctuations/oscillations, nevertheless still suitable to show evidence of ground deformation at short time scales. This is a good compromise between too much detail and a too generalized overview. Thus, the obtained RMS represents an estimate of the daily average ground tilt.

We estimated the cross-correlation among RMS tilt, rainfall amount and earthquake number in order to obtain an indication of the similarity of the time series as a function of the time lag (Chatfield, 2013). The normalized cross-correlation function is defined as:

$C_{xy}\left(k\right)=\frac{1}{N\sqrt{C_{xx}\left(0\right)}\sqrt{C_{yy}\left(0\right)}}{\displaystyle \sum }_{i=1}^{N-k}\left(x_i-{\mu }_x\right)\left(y_{i+k}-{\mu }_y\right)$

where N is the number of samples, (μ_x, μ_y) and $\left(\sqrt{C_{xx}\left(0\right)},\sqrt{C_{yy}\left(0\right)}\right)$ denote the mean and the zero lag autocorrelation of the time series, and k < N are positive integer values. The maxima (or minima if the signals are negatively correlated) of the cross-correlation function correspond to the points in time where the signals are best aligned. It is important to remark that correlation tests a relationship between two variables, but it does not imply causation.

Results

A generalized overview of the temporal pattern of ground deformation, earthquakes and rainfall throughout the entire time interval of analysis is shown in Figure 2, in which we compare the daily RMS of both the components of the ground tilt, and the monthly distributions of VTs and rainfall amount. A simple visual inspection shows the seasonality of the RMS tilt at both Vesuvius and Campi Flegrei areas, reflecting the pattern of the hydrological year (Petrosino et al., 2018). The comparison among the time series indicates a possible correlation among the three physical observables (especially between tilt and rainfall), which appear more evident at Campi Flegrei. The patterns of the RMS tilt (nearly overlapped in Figure 2) are almost independent of the site, as confirmed by the estimates of the cross-correlation matrix among the different tilt components and sensors (Supplementary Figure S2). Such analysis provides correlation coefficients greater than 0.75 for the tiltmeters installed at both areas (except for TRC site, which has a slight lower value), suggesting a common origin of the observed daily average tilt.

FIGURE 2

FIGURE 2. Daily RMS tilt and monthly distribution of rainfall (purple line) and earthquakes (yellow bars) at Vesuvius (A) and Campi Flegrei (B). The light-grey boxes in the background mark the wet season of the hydrological year. For Campi Flegrei, the RMS tilt calculated over a 30-day time window and averaged over the components of all instruments is also shown (black line).

A further observation regards the RMS values of the tilt at Campi Flegrei which are higher than those recorded at Vesuvius: with the exception of few peaks, RMS is far below <0.05 μrad/day in the dry season and <0.1 μrad/day at Vesuvius, while at Campi Flegrei these thresholds are often exceeded. This behavior suggests that the daily induced ground response to rainfall is more prominent at Campi Flegrei.

The cross-correlation analysis of the daily sampled time series provides further insights. In order to have a mean spatial estimate of the cross-correlation functions, at Campi Flegrei we averaged the RMS of the tilt over all the components and instruments. At Vesuvius we excluded the time series recorded at TRC site from the average estimate, considering its lower correlation with the signals acquired at the other three instruments. Moreover, at Campi Flegrei we extracted from the seismic catalogue only earthquakes with inter-event times less than 1 day because it has been inferred that this clustered seismicity (often grouped in intense seismic swarms) is closely related with exogenous processes, while the background seismicity (inter-event times >1 day) accounts for endogenous dynamics (Petrosino et al., 2018).

The results are shown in Figure 3. In general, at Vesuvius the correlations among the time series are low (C(k) < 0.2). Maxima values of correlation between tilt and rain occur at lags [0–4] and [19–24]. Tilt and VTs as well as rain and VTs are almost uncorrelated at least at the daily time scale (Figure 3A). A similar analysis performed only on the ground tilt recorded at TRC site confirms the low degree of correlation among tilt, rain and VTs (Supplementary Figure S3).

FIGURE 3

FIGURE 3. (A) On the left, daily values of RMS tilt, rainfall and earthquake number for Vesuvius. On the right, cross-correlation functions among the three time series (lag = 1 corresponds to 1 day). The blue horizontal lines mark the confidence limits (95%), above which the correlation between the time-series is considered significant. (B) On the left, daily values of RMS tilt, rainfall and earthquake number for Campi Flegrei. On the right, cross-correlation functions among the three time series.

The cross-correlation among rainfall, ground tilt and earthquakes at Campi Flegrei shows a different pattern. From Figure 3B, the maximum value (0.53) of cross-correlation between RMS tilt and the rainfall daily time series is at zero lag (within the 1-day uncertainty which corresponds to the sampling rate of the available rainfall data) suggesting a nearly instantaneous response of the ground deformation to meteoric input. Then, the cross-correlation function smoothly decreases during the first 10 days, indicating a persistent relationship between the two observables. Moreover, although seismicity and tilt, as well as seismicity and rainfall, appear poorly correlated, the corresponding cross-correlation functions have both a maximum value (above the 95% significance level) at a time lag of 11 days.

In order to take into account the eventual delay related to meteoric water diffusion into the ground, at Campi Flegrei data were further analysed by considering the RMS tilt, rainfall amount and VT number calculated over a time window of 30 days. Then we estimated the cross-correlation coefficients from the linear regression between rain and tilt as well as rain and earthquakes (Figure 4A). We observed an increase in both the cross-correlation values (0.77 rain and tilt, 0.36 rain and VTs) on 30 days, compared with the daily values. The obtained regression relationship between tilt and rain is:

$y=0.026+0.000086x$

FIGURE 4

FIGURE 4. Results of the regression analysis at Campi Flegrei. (A) On the left, linear regression between rainfall amount and RMS tilt. On the right, linear regression between rainfall amount and VT number. The red dotted lines mark the 95% confidence interval. All the variables have been estimated over a 30-day moving window. (B) On the left, linear regression between rainfall amount and RMS tilt, considering only rainfall values exceeding the 60th percentile of the 6-years distribution. On the right, linear regression between rainfall amount and VT number using the same rain threshold. The red dotted lines mark the 95% confidence interval. All the variables have been estimated over 30-day moving windows.

The equation provides an intercept value (0.026 μrad) which is consistent with the 30-day RMS tilt observed during dry seasons (black line in Figure 2B). In addition, referring to the time interval of 30 days, the regression law predicts an increase of 0.01 μrad/mm of the ground tilt for each increase of 120 mm of rain (which is the median values of the 6-year distribution of the 30-day precipitation amounts; Supplementary Figure S4). Such an increase is about two orders of magnitude lower that the monthly variations of the tilt rates related to the volcanic uplift (Supplementary Figure S1; https://www.ov.ingv.it/index.php/monitoraggio-e-infrastrutture/bollettini-tutti/bollett-mensili-cf).

We also investigated the possibility of a threshold effect related to heavy and/or long-lasting precipitations. Because the average rainfall amount lies between the 45th and 55th percentiles of the total distribution (Knapp et al., 2015), we extracted only those RMS tilt and earthquake data corresponding to rainfall exceeding the 60th percentile (143 mm in 30 days) of the distribution calculated over 6 years (Supplementary Figure S4). In this case, the regression of the 30-day-average data provides correlation coefficients of 0.64 and 0.63 for rain and tilt and rain and VTs, respectively (Figure 4B). While the rain and VT cross-correlation increases considering the rainfall threshold, the rain and tilt cross-correlation is slightly higher when no threshold is applied. It is likely that the threshold introduction acts as a cut off for the nearly-instantaneous induced deformation effect and accounts mainly for the contribution related to a diffusive process.

Discussion

Rainfall-induced ground deformations have been observed in many areas, especially for very shallow borehole instrument installations (Goulty, N. R., et al., 1979; Sakata and Sato, 1986; Roeloffs et al., 1989; Westerhaus and Welle, 2002; Kümpel et al., 2001; Lesparre et al., 2017). In our study cases, we notice a different behavior of the two volcanic areas of Vesuvius and Campi Flegrei. At Vesuvius a relationship between rainfall and ground deformation appears on seasonal time scale, as the RMS tilt shows amplitude variations which reflect the pattern of wet and dry seasons of the hydrological year (Figure 2A). On the daily time scale, the correlation between rainfall and ground tilt seems low. Moreover, no correlation with earthquakes is further observed at any of the investigated time scales.

At Campi Flegrei both the seasonal patterns of ground tilt (Figure 2B) and earthquake rates (Figure 2B and Supplementary Figure S5) have a striking similarity with the hydrological cycles, as the RMS tilt and VT number are systematically higher during the wet season. Here the seasonal tilt variations are more pronounced compared with those of Vesuvius. In addition, an evident correlation among rainfall, ground tilt and earthquakes also arises on the time scale of days. In particular, seismicity correlates with both rainfall and tilt (Figure 3B), with a delayed response which is focused at best when considering cumulative 30-days values and precipitation threshold (Figure 4). This suggests that a complex interplay among the three observables exists and that rainfall could have an impact on the dynamics of the system, as already pointed out in recent papers (Petrosino et al., 2018; Ricco et al., 2019; Petrosino et al., 2020b).

Several mechanisms have been proposed to take into account the possible hydromechanical coupling between rain and ground deformation observed in both tectonic and volcanic environments: among them, the most relevant are: 1) water loading essentially due to the ground compression by the added water mass, and 2) infiltration and diffusion processes which are able to modify the pore pressure of the medium (Kümpel et al., 2001; Westerhaus and Welle, 2002; Lesparre et al., 2017). The theory of poroelasticity (Wang, 2017; Wang and Kumpel, 2003) is able to explain the physical process at the origin of ground deformation related to water infiltration. According to this model, the tilt amplitude is proportional to the pore pressure gradient and depends on the rock poroelastic parameters. In the framework of the poroelasticity theory, variations of the pore-fluid pressure also play an important role in the generation of earthquakes, especially in hydrothermal/volcanic environments (Madonia et al., 2008; Savage et al., 2015; Ricco et al., 2021). Triggering of fluid-induced seismicity can be activated by diffusive processes related to pore pressure and permeability changes in porous and/or fractured saturated rocks (Shapiro et al., 2003; Andajani et al., 2020).

At Campi Flegrei such a mechanism has been invoked to explain the rainfall triggering of VT earthquakes (Scafetta and Mazzarella, 2021). The authors suggest that rainfall water percolates into the highly fractured and water-saturated rocks, mixing with deeper hot hydrothermal fluids, thus inducing pressure and permeability changes. By using a time-delayed rain function to model the diffusion process inside the soil, a delay of some days to 1–2 weeks is estimated so that water percolation affects seismicity in the depth range of 0–2.5 km. In addition, statistical results indicate that time-delayed models, compared to rapid seismic response to water infiltration, provide the highest correlations between the rainfall and seismic events at Campi Flegrei and thus represent the most probable scenario (Scafetta and Mazzarella, 2021).

From our analyses of geophysical time series recorded at Campi Flegrei, both a nearly prompt response to rainfall and a delayed one arise. As inferred in past studies the effect of the rain on ground tilt can be both immediate and delayed (Wyatt and Berger, 1980; Westernhaus and Welle, 2002; Garcia et al., 2010; Lesparre et al., 2017). Recently, Meurers et al. (2021) show that water accumulation on the terrain surface causes short-term (a few hours) tilt and gravity changes, while long-term (>a few days/weeks) variations occur frequently after long-lasting and/or heavy rain. On the other hand, seismic activity often appears as a lagged feedback to precipitations (Muço, 1999; Heinzl et al., 2006; Craig et al., 2017; Scafetta and Mazzarella, 2021). Therefore, on the basis of the existing mechanisms of hydrological response, we hypothesize a combined action of water loading (which is almost instantaneous) mainly effective on the ground deformation, and a diffusion process (on the time scale of 10–30 days), which also affect (besides the tilt) the seismicity. Indeed, the delay of 11 days of earthquake occurrence with respect to rain and ground deformation could likely be related to the time of the percolation of the shallow meteoric fluid in the fracture system of the medium. This is supported by the increase of the cross-correlation values cumulated on 30 days. We also have to consider an additional possible effect of the ground deformation which under the water loading and diffusive process generates itself elastic stress on the rocks, in turn favoring the occurrence of VTs.

Interestingly, the introduction of a threshold for the precipitation amount further increases the correlation between rain and VT number, suggesting that below a certain rainfall rate the diffusion process is only weakly effective in triggering earthquakes. This is consistent with the results of Scafetta and Mazzarella (2021) who found that seismic swarms are more likely to occur during wet days rather than in dry ones (see also Supplementary Figure S5). On the other hand, the correlation between rain and tilt slightly decreases when considering a minimum value for the precipitation amount, likely because the threshold could act as a cut off for the nearly instantaneous water loading effect of light rain on the ground deformation.

The results of the present paper provide the ground for systematical investigation of the link among tilt, seismicity and rain, identifying the RMS of instantaneous tilt on short/medium time scale (days, weeks) as a suitable observable which can better highlight (at least at the first-order) ground response to rainfall. Moreover, we have shown how the RMS pattern can enhance the seasonal behavior of the ground tilt time series; therefore it can be used as a quick indicator of seasonality in ground deformation data.

Interesting inferences can be obtained by cross-correlating the rainfall, RMS tilt and earthquake rates. Despite the assumption of a simple linear model to describe the ground response to rain (Garcia et al., 2010), we were able to highlight basic differences between Vesuvius and Campi Flegrei. The clearer response observed at Campi Flegrei, compared to Vesuvius, could be likely related to differences in the fracture density of the medium as well as in rock porosity and permeability (Zollo et al., 2006; Andajani et al., 2020; Ricco et al., 2021). Moreover, at Campi Flegrei the presence of a shallower aquifer, water saturated rocks and geothermal fluids (Gresse et al., 2017) can amplify the rain-induced ground response (Kümpel et al., 2001). The temperatures associated with the hydrothermal system (de Lorenzo et al., 2001; Zollo et al., 2006) and the presence of pressurized gases discharging at Solfatara affect the infiltration and propagation of meteoric water through porous rock and fractures, in turn influencing the poroelastic response of the volcanic system. As a result, Campi Flegrei caldera, differently from Vesuvius, is prone to respond to the stress variations related to pore-fluid pressure changes.

Finally, besides the medium properties and the presence of an active shallow hydrothermal system, a further factor which can control the response to rain-induced hydrological variations is the state of the volcano. Vesuvius is characterized by a low-level background activity, poorly influenced by external forcing. On the contrary, at Campi Flegrei endogenous processes related to the ongoing unrest likely drive the system towards metastable conditions. Although the rain-induced deformation is two orders of magnitude lower than that related to the volcanic uplift, we can reasonably infer that rainfall can have an impact in modulating the dynamics of such a metastable system, which is more sensitive to exogenous triggering (Petrosino et al., 2018).

The present outcomes are the first steps for the future development of a finer modeling of water load and diffusive processes related to rainfall phenomena. Additional factors such as duration and intensity of the precipitations, pre-rainfall moisture conditions of the soil, rock physical parameters (fracture density, porosity, permeability) and threshold non linear mechanisms (Yamauchi 1987; Dal Moro and Zadro, 1998) should be eventually taken into account with the ultimate goal to determine a more accurate as possible transfer function between rainfall and ground response and seismicity.

Data Availability Statement

The dataset of the daily RMS of the ground tilt at Vesuvius and Campi Flegrei can be found in the Zenodo repository https://doi.org/10.5281/zenodo.5196688. The INGV-OV seismic catalogue is available at the https://www.ov.ingv.it/index.php/cataloghi-sismici-vulcani-napoletani?category[0]=73&category_children=1&own=0. Rainfall amount data were downloaded by the historical archive of the website https://www.3bmeteo.com/meteo/ora/storico.

Author Contributions

SP and CR conceived the research study, analyzed data and interpreted the results. IA handled the processing of raw tiltmeter data and drew Figure 1. SP wrote the first draft; all the authors read, corrected and approved the final article.

Funding

This work is part of the INGV research project LAST: Linked Analysis of Seismic and Tiltmeter observables (Ricerca Libera 2019 - Funding n. 53/2020). Publication fee has been supported by the “Institutional Funding for Editorial Activity” of INGV, Sezione di Napoli - Osservatorio Vesuviano.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We acknowledged the INGV-OV Seismic Laboratory for having compiled the earthquake catalogue, and the INGV-OV Laboratory of Geomatics and Thermography (LGT) for having provided the base maps of Figure 1. We wish to thank Vincenzo Augusti for providing regular maintenance of the tiltmeter networks. Two reviewers are acknowledged for their suggestions, which improved the clarity of the article.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.758602/full#supplementary-material

References

Andajani, R. D., Tsuji, T., Snieder, R., and Ikeda, T. (2020). Spatial and Temporal Influence of Rainfall on Crustal Pore Pressure Based on Seismic Velocity Monitoring. Earth, Planets and Space 72 (1), 1–17. doi:10.1186/s40623-020-01311-1

CrossRef Full Text | Google Scholar

A. Zollo, P. Capuano, and M. Corciulo (Editors) (2006). Geophysical Exploration of the Campi Flegrei (Southern Italy) Caldera'Interiors: Data, Methods and Results (Napoli: GNV (National Group of Vulcanology), Doppiavoce Editore).

Google Scholar

Banks, N. G., Carvajal, C., Mora, H., and Tryggvason, E. (1990). Deformation Monitoring at Nevado del Ruiz, Colombia-October 1985-March 1988. J. Volcanol. Geotherm. Res. 41 (1–4), 269–295. doi:10.1016/0377-0273(90)90092-T

CrossRef Full Text | Google Scholar

Bellucci Sessa, E., Castellano, M., and Ricciolino, P. (2021). GIS Applications in Volcano Monitoring: the Study of Seismic Swarms at the Campi Flegrei Volcanic Complex, Italy. Adv. Geosci. 52, 131–144. doi:10.5194/adgeo-52-131-2021

CrossRef Full Text | Google Scholar

Bianco, F., Cusano, P., Petrosino, S., Castellano, M., Buonocunto, C., Capello, M., et al. (2005). Small-aperture Array for Seismic Monitoring of Mt. Vesuvius. Seismological Res. Lett. 76 (3), 344–355. doi:10.1785/gssrl.76.3.344

CrossRef Full Text | Google Scholar

Bragato, P. L. (2015). Italian Seismicity and Vesuvius' Eruptions Synchronize on a Quasi 60 year Oscillation. Earth Space Sci. 2 (5), 134–143. doi:10.1002/2014EA000030

CrossRef Full Text | Google Scholar

Braitenberg, C. (1999). The Hydrologic Induced Strain-Tilt Signal: a Review. Bull. D’inf. Marées Terrestres 131 (10), 171–210.

Google Scholar

Caputo, T., Cusano, P., Petrosino, S., Sansivero, F., and Vilardo, G. (2020). Spectral Analysis of Ground thermal Image Temperatures: what We Are Learning at Solfatara Volcano (Italy). Adv. Geosci. 52, 55–65. doi:10.5194/adgeo-52-55-2020

CrossRef Full Text | Google Scholar

Chatfield, C. (2013). The Analysis of Time Series: Theory and Practice. New York: Springer.

Google Scholar

Cioni, R., Bertagnini, A., Santacroce, R., and Andronico, D. (2008). Explosive Activity and Eruption Scenarios at Somma-Vesuvius (Italy): towards a New Classification Scheme. J. Volcanology Geothermal Res. 178 (3), 331–346. doi:10.1016/j.jvolgeores.2008.04.024

CrossRef Full Text | Google Scholar

Craig, T. J., Chanard, K., and Calais, E. (2017). Hydrologically-Driven Crustal Stresses and Seismicity in the New Madrid Seismic Zone. Nat. Commun. 8 (1), 1–11. doi:10.1038/s41467-017-01696-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Cusano, P., Caputo, T., De Lauro, E., Falanga, M., Petrosino, S., Sansivero, F., et al. (2021). Tracking the Endogenous Dynamics of the Solfatara Volcano (Campi Flegrei, Italy) through the Analysis of Ground Thermal Image Temperatures. Atmosphere 12 (8), 940. doi:10.3390/atmos12080940

CrossRef Full Text | Google Scholar

Cusano, P., Petrosino, S., and Saccorotti, G. (2008). Hydrothermal Origin for Sustained Long-Period (LP) Activity at Campi Flegrei Volcanic Complex, Italy. J. Volcanology Geothermal Res. 177 (4), 1035–1044. doi:10.1016/j.jvolgeores.2008.07.019

CrossRef Full Text | Google Scholar

Dal Moro, G., and Zadro, M. (1998). Subsurface Deformations Induced by Rainfall and Atmospheric Pressure: Tilt/strain Measurements in the NE-Italy Seismic Area. Earth Planet. Sci. Lett. 164 (1-2), 193–203. doi:10.1016/s0012-821x(98)00203-9

CrossRef Full Text | Google Scholar

De Lauro, E., Falanga, M., and Petrosino, S. (2012). Study on the Long-Period Source Mechanism at Campi Flegrei (Italy) by a Multi-Parametric Analysis. Phys. Earth Planet. Interiors 206-207, 16–30. doi:10.1016/j.pepi.2012.06.006

CrossRef Full Text | Google Scholar

De Lauro, E., Petrosino, S., Ricco, C., Aquino, I., and Falanga, M. (2018). Medium and Long Period Ground Oscillatory Pattern Inferred by Borehole Tiltmetric Data: New Perspectives for the Campi Flegrei Caldera Crustal Dynamics. Earth Planet. Sci. Lett. 504, 21–29. doi:10.1016/j.epsl.2018.09.039

CrossRef Full Text | Google Scholar

de Lorenzo, S., Gasparini, P., Mongelli, F., and Zollo, A. (2001). Thermal State of the Campi Flegrei Caldera Inferred from Seismic Attenuation Tomography. J. geodynamics 32 (4-5), 467–486. doi:10.1016/s0264-3707(01)00044-8

CrossRef Full Text | Google Scholar

Falanga, M., and Petrosino, S. (2012). Inferences on the Source of Long-Period Seismicity at Campi Flegrei from Polarization Analysis and Reconstruction of the Asymptotic Dynamics. Bull. Volcanol 74 (6), 1537–1551. doi:10.1007/s00445-012-0612-2

CrossRef Full Text | Google Scholar

Farquharson, J. I., and Amelung, F. (2020). Extreme Rainfall Triggered the 2018 Rift Eruption at KĪlauea Volcano. Nature 580 (7804), 491–495. doi:10.1038/s41586-020-2172-5

PubMed Abstract | CrossRef Full Text | Google Scholar

García, A., Hördt, A., and Fabian, M. (2010). Landslide Monitoring with High Resolution Tilt Measurements at the Dollendorfer Hardt Landslide, Germany. Geomorphology 120 (1-2), 16–25. doi:10.1016/j.geomorph.2009.09.011

CrossRef Full Text | Google Scholar

Goulty, N. R., Davis, P. M., Gilman, R., and Motta, N. (1979). Meteorological Noise in Wire Strainmeter Data from Parkfield, California. Bull. Seismological Soc. America 69 (6), 1983–1988. doi:10.1785/bssa0690061983

CrossRef Full Text | Google Scholar

Gresse, M., Vandemeulebrouck, J., Byrdina, S., Chiodini, G., Revil, A., Johnson, T. C., et al. (2017). Three-Dimensional Electrical Resistivity Tomography of the Solfatara Crater (Italy): Implication for the Multiphase Flow Structure of the Shallow Hydrothermal System. J. Geophys. Res. Solid Earth 122 (11), 8749–8768. doi:10.1002/2017jb014389

CrossRef Full Text | Google Scholar

Hainzl, S., Kraft, T., Wassermann, J., Igel, H., and Schmedes, E. (2006). Evidence for Rainfall‐triggered Earthquake Activity. Geophys. Res. Lett. 33 (19). doi:10.1029/2006gl027642

CrossRef Full Text | Google Scholar

Knapp, A. K., Hoover, D. L., Wilcox, K. R., Avolio, M. L., Koerner, S. E., La Pierre, K. J., et al. (2015). Characterizing Differences in Precipitation Regimes of Extreme Wet and Dry Years: Implications for Climate Change Experiments. Glob. Change Biol. 21 (7), 2624–2633. doi:10.1111/gcb.12888

PubMed Abstract | CrossRef Full Text | Google Scholar

Kümpel, H.-J., Lehmann, K., Fabian, M., and Mentes, G. (2001). Point Stability at Shallow Depths: Experience from Tilt Measurements in the Lower Rhine Embayment, Germany, and Implications for High-Resolution GPS and Gravity Recordings. Geophys. J. Int. 146 (3), 699–713. doi:10.1046/j.1365-246x.2001.00494.x

CrossRef Full Text | Google Scholar

Lesparre, N., Boudin, F., Champollion, C., Chéry, J., Danquigny, C., Seat, H. C., et al. (2017). New insights on fractures deformation from tiltmeter data measured inside the Fontaine de Vaucluse karst system. Geophys. J. Int. 208 (3), 1389–1402. doi:10.1093/gji/ggw446

CrossRef Full Text | Google Scholar

Madonia, P., Federico, C., Cusano, P., Petrosino, S., Aiuppa, A., and Gurrieri, S. (2008). Crustal Dynamics of Mount Vesuvius from 1998 to 2005: Effects on Seismicity and Fluid Circulation. J. Geophys. Res. Solid Earth 113 (B5). doi:10.1029/2007jb005210

CrossRef Full Text | Google Scholar

Matthews, A. J., Barclay, J., Carn, S., Thompson, G., Alexander, J., Herd, R., et al. (2002). Rainfall‐induced Volcanic Activity on Montserrat. Geophys. Res. Lett. 29 (13), 22–31. doi:10.1029/2002gl014863

CrossRef Full Text | Google Scholar

Matthews, A. J., Barclay, J., and Johnstone, J. E. (2009). The Fast Response of Volcano-Seismic Activity to Intense Precipitation: Triggering of Primary Volcanic Activity by Rainfall at Soufrière Hills Volcano, Montserrat. J. volcanology geothermal Res. 184 (3-4), 405–415. doi:10.1016/j.jvolgeores.2009.05.010

CrossRef Full Text | Google Scholar

Mazzarella, A., and Scafetta, N. (2016). Evidences for Higher Nocturnal Seismic Activity at the Mt. Vesuvius. J. Volcanology Geothermal Res. 321, 102–113. doi:10.1016/j.jvolgeores.2016.04.026

CrossRef Full Text | Google Scholar

Meurers, B., Papp, G., Ruotsalainen, H., Benedek, J., and Leonhardt, R. (2021). Hydrological Signals in Tilt and Gravity Residuals at Conrad Observatory (Austria). Hydrol. Earth Syst. Sci. 25 (1), 217–236. doi:10.5194/hess-25-217-2021

CrossRef Full Text | Google Scholar

Muço, B. (1999). Statistical Investigation on Possible Seasonality of Seismic Activity and Rainfall-Induced Earthquakes in Balkan Area. Phys. earth Planet. interiors 114 (3-4), 119–127. doi:10.1016/s0031-9201(99)00051-5

CrossRef Full Text | Google Scholar

Petrosino, S., Cusano, P., and Madonia, P. (2018). Tidal and Hydrological Periodicities of Seismicity Reveal New Risk Scenarios at Campi Flegrei Caldera. Sci. Rep. 8 (1), 13808–13812. doi:10.1038/s41598-018-31760-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrosino, S., and Cusano, P. (2020a). Low Frequency Seismic Source Investigation in Volcanic Environment: the Mt. Vesuvius Atypical Case. Adv. Geosci. 52, 29–39. doi:10.5194/adgeo-52-29-2020

CrossRef Full Text | Google Scholar

Petrosino, S., Ricco, C., De Lauro, E., Aquino, I., and Falanga, M. (2020b). Time Evolution of Medium and Long-Period Ground Tilting at Campi Flegrei Caldera. Adv. Geosci. 52, 9–17. doi:10.5194/adgeo-52-9-2020

CrossRef Full Text | Google Scholar

Ricco, C., Petrosino, S., Aquino, I., Cusano, P., and Madonia, P. (2021). Tracking the Recent Dynamics of Mt. Vesuvius from Joint Investigations of Ground Deformation, Seismicity and Geofluid Circulation. Sci. Rep. 11 (1), 965–1014. doi:10.1038/s41598-020-79636-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricco, C., Aquino, I., Augusti, V., D'auria, L., Del Gaudio, C., and Scarpato, G. (2018). Improvement and Development of the Tiltmetric Monitoring Networks of Neapolitan Volcanoes. Ann. Geophys. 61 (1). doi:10.4401/ag-7496

CrossRef Full Text | Google Scholar

Ricco, C., Aquino, I., Borgstrom, S. E. P., and Del Gaudio, C. (2013). 19 Years of Tilt Data on Mt. Vesuvius: State of the Art and Future Perspectives. Ann. Geophys. 56 (4), S0453. doi:10.4401/ag-6459

CrossRef Full Text | Google Scholar

Ricco, C., Petrosino, S., Aquino, I., Del Gaudio, C., and Falanga, M. (2019). Some Investigations on a Possible Relationship between Ground Deformation and Seismic Activity at Campi Flegrei and Ischia Volcanic Areas (Southern Italy). Geosciences 9 (5), 222. doi:10.3390/geosciences9050222

CrossRef Full Text | Google Scholar

Roeloffs, E. A., Burford, S. S., Riley, F. S., and Records, A. W. (1989). Hydrologic Effects on Water Level Changes Associated with Episodic Fault Creep Near Parkfield, California. J. Geophys. Res. 94 (B9), 12387–12402. doi:10.1029/jb094ib09p12387

CrossRef Full Text | Google Scholar

Saccorotti, G., Petrosino, S., Bianco, F., Castellano, M., Galluzzo, D., La Rocca, M., et al. (2007). Seismicity Associated with the 2004–2006 Renewed Ground Uplift at Campi Flegrei Caldera, Italy. Phys. Earth Planet. Interiors 165 (1-2), 14–24. doi:10.1016/j.pepi.2007.07.006

CrossRef Full Text | Google Scholar

Sakata, S., and Sato, H. (1986). Borehole-type Tiltmeter and Three-Component Strainmeter for Earthquake Prediction. J,Phys,Earth 34 (Suppl. ment), S129–S140. doi:10.4294/jpe1952.34.supplement_s129

CrossRef Full Text | Google Scholar

Savage, M. K., Ferrazzini, V., Peltier, A., Rivemale, E., Mayor, J., Schmid, A., Massin, F., Got, J.-L., Battaglia, J., DiMuro, A., Staudacher, T., Rivet, D., Taisne, B., and Shelley, A. (2015). Seismic anisotropy and its precursory change before eruptions at Piton de la Fournaise volcano, La Réunion. J. Geophys. Res. Solid Earth 120 (5), 3430–3458. doi:10.1002/2014jb011665

CrossRef Full Text | Google Scholar

Sbrana, A., Cioni, R., Marianelli, P., Sulpizio, R., Andronico, D., and Pasquini, G. (2020). Volcanic Evolution of the Somma-Vesuvius Complex (Italy). J. Maps 16 (2), 137–147. doi:10.1080/17445647.2019.1706653

CrossRef Full Text | Google Scholar

Scafetta, N., and Mazzarella, A. (2021). On the Rainfall Triggering of Phlegraean Fields Volcanic Tremors. Water 13 (2), 154. doi:10.3390/w13020154

CrossRef Full Text | Google Scholar

Shapiro, S. A., Patzig, R., Rothert, E., and Rindschwentner, J. (2003). Triggering of Seismicity by Pore-Pressure Perturbations: Permeability-Related Signatures of the Phenomenon. in” Thermo-Hydro-Mechanical Coupling in Fractured Rock. Basel: Birkhäuser, 1051–1066. doi:10.1007/978-3-0348-8083-1_16

CrossRef Full Text | Google Scholar

Violette, S., de Marsily, G., Carbonnel, J. P., Goblet, P., Ledoux, E., Tijani, S. M., et al. (2001). Can rainfall trigger volcanic eruptions? A mechanical stress model of an active volcano: 'Piton de la Fournaise', Reunion Island. Terra Nova 13 (1), 18–24. doi:10.1046/j.1365-3121.2001.00297.x

CrossRef Full Text | Google Scholar

Wang, H. F. (2017). Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology. Princeton: Princeton University Press.

Google Scholar

Wang, R., and Kümpel, H. J. (2003). Poroelasticity: Efficient Modeling of Strongly Coupled, Slow Deformation Processes in a Multilayered Half-Space. Geophysics 68 (2), 705–717. doi:10.1190/1.1567241

CrossRef Full Text | Google Scholar

Westerhaus, M., and Welle, W. (2002). Environmental Effects on Tilt Measurements at Merapi Volcano. Bull. Inf. Marees Terr 137, 10917–10926.

Google Scholar

Wyatt, F., and Berger, J. (1980). Investigations of Tilt Measurements Using Shallow Borehole Tiltmeters. J. Geophys. Res. 85 (B8), 4351–4362. doi:10.1029/jb085ib08p04351

CrossRef Full Text | Google Scholar

Yamauchi, T. (1987). Anomalous Strain Response to Rainfall in Relation to Earthquake Occurrence in the Tokai Area, Japan. J,Phys,Earth 35 (1), 19–36. doi:10.4294/jpe1952.35.19

CrossRef Full Text | Google Scholar