A record of changes in the Gran Sasso groundwater before, during and after the 2016 Amatrice earthquake, central Italy
Introduction
We performed continuous recordings (May 2015 - January 2017) of hydraulic pressure and electrical conductivity of groundwater in the 190 m-long horizontal S13 borehole drilled next to the deep underground laboratories of Gran Sasso (LNGS-INFN), located in the core of the Gran Sasso carbonate aquifer (central Italy) at a distance of about 39 km south-eastward from the 24 August 2016 Amatrice earthquake (6.0 Mw) epicenter. Using a 3-channel, 24-bit ADC we achieved a sampling rate of groundwater physical properties up to 50 Hz for each channel. We focused on the analysis of data recorded before, during and after the Amatrice earthquake, describing and discussing in detail the evidence for significant hydraulic pressure and electrical conductivity anomalies recorded before the main shock. We identified unambiguous signals in the hydraulic pressure data starting on 19 August, i.e. five days before the 24 August mainshock. A more careful analysis allowed us to detect the inception of a weak change up to 40 days before the Amatrice earthquake and a significant variation in the electrical conductivity data about 60 days before. The data revealed highly dynamic aquifer behaviour associated with the uprising of geogas probably related to the preparation stage of the Amatrice earthquake.
Earthquake prediction is a widely recognized goal as well as a challenging scientific problem because of its complexity and its social impact. It is easy to find many papers, reviews, debates, discussions and comments in the scientific literature1-8. Chemical and physical groundwater parameters have been monitored in seismogenic areas worldwide, principally by measuring and analysing water level and hydrochemistry changes in wells, springs and streams to find possible correlations with local and regional seismicity and between their spatial and temporal variations and strain processes9-28. Moreover, in the past few decades many studies have suggested a connection between fault mechanics and underground fluid dynamics, and a large quantity of experimental results and models have appeared in the literature29-37. Among these studies, those concerning the fault frictional properties in the presence of fluid can provide unique insights into fluid-rock interactions34.
In this framework, we decided to monitor some physical parameters of the deep groundwater of the large Gran Sasso carbonate aquifer (Abruzzi region, central Apennines, Fig. 1) through a horizontal borehole (named S13) placed in the highway tunnel very close to the underground laboratories of Gran Sasso (LNGS-INFN). The Gran Sasso chain is placed in a seismically active area of central Apennines, as demonstrated by the 6 April 2009, Mw 6.3 L'Aquila earthquake38,39 and by the recent 24 August 2016, Mw 6.0 Amatrice earthquake40,41 (Fig. 1), and provides a unique opportunity to study an inner, unexploited portion of a very large carbonate aquifer (about 1,000 km2). This area of the central Apennines is monitored by a rather dense regional seismic network42 along with the national seismic network of INGV-ONT (Istituto Nazionale di Geofisica e Vulcanologia - Osservatorio Nazionale Terremoti)43.Main earthquakes location map, from Iside data-base53, of August 2016-January 2017 period (http://cnt.rm.ingv.it/iside); brown squares represent principal towns in the area, black lines for regional boundaries, green star for LNGS (Laboratori Nazionali del Gran Sasso) of INFN (Istituto Nazionale di Fisica Nucleare) close to the S13 horizontal borehole (distance about 200 m); from northern: red stars for 6.5 Mw of 30 October 2016, 6.0 Mw of 24 August 2016, 5.1 Mw of 18 January 2017 (09:25:40 UT), 5.5 Mw of 18 January 2017 (10:14:10 UT), 5.4 Mw of 18 January 2017 (10:25:24 UT) and 5.0 Mw of 18 January 2017 (13:33:37 UT). We also show the L'Aquila earthquake of 6 April 2009 (6.3 Mw - 01:32:40 UT).
PHOTO (COLOR)
Convincing evidence exists in the Gran Sasso chain and its neighbouring areas that high fluid pressure contributed to the rupture during the Mw 6.3 L'Aquila earthquake of 6 April 200944-47. Temporal variations of seismic velocity, anisotropy, slow-slip and ground deformation were also observed before the main event48-50. Specifically, some investigators observed clear variations in the seismic wave propagation properties about a week before the L'Aquila mainshock and inferred that fluids played a key role in the fault failure process44-48. Similarly, during the 24 August 2016, Mw 6.0 Amatrice earthquake there was evidence of hydrogeochemical and hydrogeological changes before and during the seismic sequence51,52, that included about 80,000 events53. Several geochemical anomalies were interpreted as reliable seismic precursors for an extensional tectonic setting51.
We made our observations from a tunnel at 965 m a.s.l. (Fig. 2), in the core of the Gran Sasso aquifer beneath a 1,400 m-thick layer of rocks. The horizontal S13 borehole (Figs 2 and 3) was drilled around the end of the 1980s in a hall excavated in Meso-Cenozoic limestones and dolostones close to the underground laboratories of Gran Sasso (LNGS-INFN). It has a horizontal length of slightly over 190 m and dips gently upward by ~5° (Fig. 3). The S13 borehole intercepts a thrust fault near its end (Figs 2 and 3). The initial 175 m of S13 borehole are tubed with a well casing, while the last 10 m drain inside the upper Triassic dolomite54,55.(a) Scheme (not in scale) of the Gran Sasso aquifer transversal to the highway tunnels and passing through the underground laboratories of Gran Sasso (LNGS), borehole hall and S13 area77. UZ: Unsaturated Zone; SZ: Saturated Zone; KH: Karst Horizon; RA: Regional Aquiclude; T: permeability boundary (regional Thrust); UT: local thrust named Upper Thrust; WT: Water Table; HT: Highway Tunnels; UL LNGS: Underground Laboratories; BH: Borehole Hall; CA: Calderone glacier (high elevation water reservoir - preferential recharge area); 1: overflow spring (CP: Capo Pescara spring); 2: preferential groundwater flowpath area; 3: preferential groundwater flowing toward the UL; PR: Preferential Recharge; DR: Diffuse Recharge; S13: monitored horizontal borehole. The hydrogeological relationships in the square are showed into details in (b). (b) and (c) Detailed hydrogeological relationships between Calderone glacier acting as a water reservoir for the carbonate aquifer down below; i: hydraulic gradient; k: hydraulic conductivity (kh: karst horizon; lim: limestone; dol: dolomite). The hydrogeological relationships in the square are showed in detail in (c).
PHOTO (COLOR)
Top: geological map of underground laboratories (LNGS-INFN), right side, and borehole hall at left (redrawn from [54]). The figure shows the characteristic of the 6 horizontal boreholes (S13, S14, S15, S16, S17 and S18). The S13 borehole was monitored in this study. Bottom: geological logs of S13 and S14 boreholes (top scheme and Fig. 2 for location): stratified limestone (Scaglia Fm); detrital massive limestone (detrital Fucoidi Fm); stratified limestone with cherty layers (Terratta and Limestone with radiolarians Fms); stratified dolomite (Dolomia principale Fm); UT: upper thrust; UL: Underground Laboratories. The hydraulic pressure (MPa) measured during the boring are reported54. The black bold lines of S13 borehole represent the iron tube (150 mm of section).
PHOTO (COLOR)
Hydrogeological Setting
Most of the Abruzzi region is a seismic area characterized by normal faulting earthquakes reaching maximum intensity values up to XI MCS, corresponding to magnitudes close to 7.0; the most destructive known events occurred in 1349, 1461, 1703, 1915 and 200956. The seismic activity of the Abruzzi region is quite like that of most of the central-southern Apennines, with predominant normal faulting earthquakes generated by predominant NW-SE trending fault systems of the Apenninic chain.57-62. The Gran Sasso chain is in the northern part of the Abruzzi region and corresponds to the main topographic and geological feature of central Apennines: its culmination (Corno Grande, 2,912 m a.s.l.) is the highest peak in peninsular Italy. The LNGS-INFN underground laboratories, with a rock overburden of 1,400 m, are placed in this area together with the S13 borehole site (Fig. 2). No significant historical earthquakes are known in the Gran Sasso chain, but paleoseismological studies reported multiple surface-faulting events; the magnitude of the two most recent strong events (5th-3rd cent. BC and 6th-5th millennium BC) have been estimated at about 7.063.
During the 6 April 2009, Mw 6.3 L'Aquila earthquake and the related seismic sequence, several investigators studied groundwater post-seismic changes (spring discharge, water table and water content of main ions, 222Rn, CO2, water isotopes and Uranium) of the Gran Sasso aquifer64-68. Later, the S13 borehole was specifically selected with the goal of continuous monitoring of hydraulic pressure, electrical conductivity and temperature at a very high sampling frequency. This choice was guided by the observation that S13 borehole clearly recorded water table changes induced by the 2009 L'Aquila earthquake64,65.
The 1,000 km2-wide fractured and fault-partitioned Gran Sasso carbonate aquifer is hosting the nuclear physics underground laboratories of Gran Sasso (LNGS-INFN) and six horizontal perforations of the borehole hall, among which the S13 borehole (Fig. 3 for details)54. This aquifer is a representative carbonate aquifer of the Mediterranean domain and over the past few decades has been the subject of investigations in the framework of several projects concerning hydrogeology, hydrogeochemistry and isotope hydrology69-78. The Gran Sasso aquifer is composed of Meso-Cenozoic basin-to-slope and reef-platform carbonate rocks, which are arranged in a thrust belt geometry formed during the Upper Miocene Apenninc orogeny69,70. The thrust belt is subsequenty displaced by Plio-Quaternay normal faults generating by the NE-migration of the Tyrrhenian post-orogenic extensional front54,55,79. The Gran Sasso aquifer is hydraulically confined to the north and to the east by the Upper Miocene terrigenous lithologies as aquiclude. The permeability boundary is represented by the regional north-verging overthrust causing the overlaying of the Gran Sasso aquifer (in the hangingwall) onto terrigenous aquiclude (in the footwall) (Figs 2 and 3)70,77.
The underground laboratories (LNGS-INFN) and the borehole area intersect the upper thrust (UT in Figs 2 and 3). Its hangingwall and footwall are formed by lower Liassic limestones-upper Triassic dolostones and upper Cretaceous-upper Jurassic mud-supported detrital and cherty limestones, respectively54,55. Groundwater from the core aquifer is partially drained by the two highway tunnels and by the underground laboratories (UL LNGS in Fig. 2); it is important to stress that there are no pumping systems in the whole area, and groundwater is only passively drained from the two sides of the Gran Sasso chain78. Since the 1980s, groundwater flowing in the tunnels has been exploited as drinking water. The decrease of spring discharge measured in the period 1979-1990 was interpreted as a transient phenomenon caused by the tunnel excavation. Once the aquifer reached a new steady condition, however, progressive stabilization of discharge was expected70. Hydraulic pressure measured at the end of the 150-200 m-long horizontal boreholes falls in the range of 0.5-0.7 MPa (piezometric head of 50-70 m), except for S13 borehole, which exhibits a much higher pressure, about 2.0-2.5 MPa corresponding to a piezometric head of 200-250 m (Figs 2 and 3)54.
Instrumentation
Our experimental equipment includes a 3-channels 24-bit ADC (mod. SL06 by Sara Electronic Instruments - http://www.sara.pg.it/) set up for continuous local recording. We started data acquisition on 1 May 2015 by continuous high-frequency sampling (10-50 Hz for each channel) of groundwater hydraulic pressure (Gems sensor, 3500 series, 0-4 MPa range), temperature (Pt1000) and electrical conductivity (http://www.bc-electronics.it, mod. Sl-311). From the end of 2015 the sampling rate was fixed at 20 Hz for each channel80. Figure 4 shows the linear behaviour and fit of hydraulic pressure vs its flow of S13 borehole, Fig. 5 shows a simple hydraulic scheme of the apparatus located at the S13 borehole wellhead.Hydraulic pressure (MPa) vs flow rate (l/min) of the S13 monitored borehole. The vertical and horizontal error bars are + /− 5% of value. The linear fit leads to the following equation: P(MPa) = −0.25*flow_rate (l/min) + 2.7.
PHOTO (COLOR)
Hydraulic scheme (not in scale) of the experimental apparatus. The length of the S13 horizontal borehole is 190 m. (a) old analogic manometer; (b) hydraulic valve that is always open during the data acquisition periods; (c) hydraulic pressure sensor; (d) hydraulic valve not completely close to enable the measurement of temperature and electrical conductivity in the box (g); (e) temperature sensor; (f) electrical conductivity sensor; (g) transparent plastic container housing the temperature and electrical conductivity sensors; water is expelled when reaching about three quarters of the volume of the container.
PHOTO (COLOR)
The hydrological characteristics of the S13 borehole showed a specific peculiarity: high hydraulic pressure with respect to very low flow rate (Fig. 4). To measure temperature and electrical conductivity we needed to spill some water, but to avoid disturbance in hydraulic pressure measurement we decided to limit as much as possible the spilled flow. To this end, we installed the hydraulic system shown in Fig. 5; a manual hydraulic valve (d) in Fig. 5 permitted to regulate the water flow. The hydraulic pressure sensor was connected directly to the S13 horizontal borehole output, (c) in Fig. 5, while the other two sensors (temperature (e) and electrical conductivity (f) devices in Fig. 5) are immersed in a small box fed, (g) in Fig. 5, by a continuous supply of water at atmospheric pressure80.
We decided to work with a flow rate of 0.85 l/min (from May 2015) and then we increased it on 5 April 2016 to 1.97 l/min. The increased flow rate carried out on 5 April 2016 resulted in a pressure drop of about 0.3 MPa (Fig. 6). Under these conditions of high pressures and outlet of water from an extremely small opening, the flow rate was not in a laminar regime but in a fully developed turbulence.Hydraulic pressure (MPa) plots. Data from 03 August 2015 (doy 215) to 5 September 2016 (doy 249), about 400 days of data. At the begin of April 2016, we increased flow rate and the hydraulic pressure decreased from 2.35 MPa to 2.08 MPa. The red plot is only the translation of the bottom black line.
PHOTO (COLOR)
A seismic station of the national network53 of Istituto Nazionale di Geofisica e Vulcanologia (INGV) is located about 250 m from the S13 borehole inside the underground laboratory of Gran Sasso (LNGS-INFN). The seismic station (international code: GIGS) was deployed in the framework of the GINGER experiment81-83 and is equipped with two broadband seismometers, including a Nanometrics Trillium 240 s (see http://iside.rm.ingv.it/iside/standard/info%5fstazione.jsp?page=sta&sta=2571 for further details) and a Guralp CMG 3 T 360 s. This instrumentation is used both for continuous microseismic monitoring of the Gran Sasso chain and for recording global seismicity.
Data Discussion and Interpretation
During the months of continuous monitoring (May 2015-Jan 2017) the hydraulic pressure signal of the S13 horizontal borehole showed the recharge and depletion cycle of the Gran Sasso aquifer (Fig. 6) and displayed also the sun/moon tides (Fig. 7). The temperature followed both seasonal and daily variations, the former linked always to the hydrogeological cycle, the latter due to the closeness of the highway tunnel with the S13 hall, while the electrical conductivity did not show any significant periodical variations. Apart from these expected trends, at the end of 2015 the hydraulic pressure data showed coseismic effects related to two large earthquakes occurring during the monitored period, respectively located at global and regional distance: the 16 Sept 2015, Illapel, central Chile earthquake (Mw 8.3, distance about of 12,000 km), and the 17 Nov 2015, Lefkada, western Greece earthquake (Mw 6.5, distance about 700 km). Details of the waveforms recorded can be found in [80]; the arrival of Love and Rayleigh waves from the Chile earthquake and P- and S- waves from the Greece earthquake are very clear80. The signal recorded in the pressure channel during the earthquake represents the response induced by wave arrivals exciting the whole aquifer.Tidal signal in the hydraulic pressure data of S13 borehole. (a) data from 5 to 15 March 2016, about 6 months before the Amatrice earthquake; (b) data from 16 to 26 August 2017, about one year after the Amatrice earthquake.
PHOTO (COLOR)
A few days after the 24 August 2016 (01:36:32 UT, doy 237) Amatrice earthquake (Mw 6.0), we recovered our data to search for any changes possibly induced by the earthquake. We expected the signals to be large due the relative closeness of the Amatrice event (about 39 km to the monitored site) and its notable magnitude. As we observed an unexpected behaviour only in the hydraulic pressure and electrical conductivity channels, in the following we will refer only to these observations.
Figure 8 shows unprocessed data (at a 20 Hz sampling rate, red line) of the hydraulic pressure from 15 August (doy 229) to 5 September 2016 (doy 249) using a 60 s moving average (blue line). The presence of the average values reveals the structural change of the signal (asymmetry) that will be further illustrated in Figs 9-11 and 13. In the first days (15 August-18 August 2016) we observed the aquifer recharge superimposed on a small tidal modulation; a similar behaviour was seen again only around the end of August. Between 19 August - i.e. a week before the earthquake - and the end of August 2016 we observed a different regime of amplitude fluctuations of the hydraulic pressure. To emphasize this effect, we performed a detrend analysis in hydraulic pressure data for the May 2015 - Sept 2016 period. Figure 9a shows the detrended of hydraulic pressure of S13 borehole from April to August 2016 (about 150 days) while the Fig. 9b shows the detrend analysis only from 15 August (doy 228) to 29 August 2016 (doy 242). From the start of data collection (May 2015), detrended data fell inside the +/−0.001 MPa band (signalled by blue lines in Fig. 9a,b). This was true until 19 August 2016 (five days before Amatrice earthquake) when large and asymmetric fluctuations appeared, lasting till the end of August 2016 (Figs 8 and 9b). In the period from May 2015 to July 2016 we did not observe any variations. In Fig. 8 and in a minor way also in Fig. 6, moreover, we observe a variety of steps. In most cases the steps are up or down with durations of less than 1 minute. These are probably related to pore-pressure phenomena. Coseismic steps are comparatively slower (about ten minutes)9,13,16,18-21,23-25.Plot of the hydraulic pressure (in MPa) of S13 borehole, red line, at 20 Hz sampling rate (from 15 August, doy 228, to 4 September, doy 248, 2016). The blue line represents the 60 s moving average.
PHOTO (COLOR)
(a) Detrended of hydraulic pressure of S13 borehole from April to August 2016. (b) Detrended of hydraulic pressure of S13 borehole from 15 (doy 228) to 29 August 2016 (doy 242), shaded rectangle in (a).The blue arrow in the top of (b) represents the Amatrice earthquake occurrence. In the period from May 2015 to March 2016 we did not observe any variations.
PHOTO (COLOR)
(a) Plot of the hydraulic pressure (in MPa) at 20 Hz sampling rate recorded on 19 August 2016 (doy 232) from 10:00 to 24:00 UT (red line). The blue line represents the 60 s moving average. (b) A 30 s section of recorded data: after 13:00 UT (lower line) and after 22:00 UT (upper line) of 19 August 2016. (c) Histograms of the data distributions: the asymmetric distribution refers to one hour of recording from 15:00 to 16:00 UT of 20 August 2016 (doy 233), while the symmetric distribution refers to one hour from 12:00 to 13:00 UT of 19 August 2016 (doy 232).
PHOTO (COLOR)
Number of negative micropulses (red line) each hour from 19 August (doy 232) to 2 September 2016 (doy 246), while the blue line shows the integral trend in arbitrary scale, black vertical line represents the Amatrice earthquake occurrence and green horizontal line is the zero level. In the period from May 2015 to July 2016 we did not observe any negative micropulses.
PHOTO (COLOR)
We next focus on the regime change that occurred on 19 August 2016 (doy 232) and on the resulting anomalies in the hydraulic pressure signal characteristics from S13 borehole. Figure 10a shows the onset of this anomalous behaviour on 19 August 2016. A transition from one regime to a different regime was recorded over a period of about 6 hours (from 14:00 to 20:00 UT). This change occurred not progressively but following a series of "flip-flop" cycles. This change is marked by two distinct features: an increase in the average hydraulic pressure (about 0.004 MPa), and a more dramatic change in the fine structure of the signal. The latter effect can be better appreciated by looking at the two panels (Fig. 10b,c) representing the behaviour of hydraulic pressure data before and after the transition. The presence of asymmetric fluctuations (negative micropulses) is evident following the transition (Fig. 10b). This behaviour is shown more quantitatively in the histograms displayed in Fig. 10c. From 19 August 2016 (doy 232) onwards, this behaviour is observed systematically, although to a varying degree. Figure 11 shows the hourly number of negative micropulses (red line) from 19 August to 2 September 2016, while the blue line shows the integral trend on an arbitrary scale. It is of interest to note that, starting from about 36 hours before the 24 August mainshock there is a large increase of negative micropulses, culminating at the time of occurrence of the Amatrice earthquake. Moreover, in the period from May 2015 to July 2016 we did not observe any negative micropulses.
To explain the presence of these negative micropulses we go back to the details of the hydraulic setup. The flux of water from the S13 borehole is determined by the balance of two effects: the high impedance of the source associated with the high level of fracturing of the rocks crossed by the S13 borehole, and the small area of the valve orifice (O(1) mm2). The density and viscosity of the water limit the flow rate to about 2 l/min with an inlet/outlet pressure ratio of about 25 (and flow speed of about 10 m/s). Any increase in flow rate implies a sudden decrease of hydraulic pressure as showed by the trend of Fig. 4. Then when a macroscopic bubble (O(1) ml at 1 bar) of geogas reaches the orifice of the valve an instantaneous increase of flow rate occurs due to the sudden decrease of density and viscosity. This results in a rapid decrease of the hydraulic pressure, as measured above the valve, producing a negative micropulse in our recording (on the order of 2 × 10−3 MPa). Therefore, the counting rate of negative micropulses is a proxy for the abundance of gas bubbles in groundwater.
Figure 12 focuses on 24 August 2016 from 01:36:30 to 01:38:00 UT, when the Amatrice earthquake occurred. Data from the 3-components GIGS broadband seismic station are shown (red lines) next to the hydraulic pressure signal (black line) from S13 borehole as a reference. Figure 12 shows a main feature: the remarkable variation of hydraulic pressure (0.2 MPa pp) following the arrival of the S-waves; afterwards, a slow increase of about 0.02 MPa (corresponding to a piezometric head change of about 2 m) were recorded over a period of about 5 hours following the end of the earthquake signal (Fig. 8). Although similar observations have already been described in the literature9,84, the results we obtained exhibit a remarkable magnitude and a great level of detail80. The novelty of our time series of hydraulic pressure lies also in the presence of unambiguous signals starting several days before the mainshock of 24 August 2016.Waveforms corresponding to the Amatrice earthquake, 24 August 2016, 01:36:32 UT, doy 237, (Mw 6.0); the red lines represent the vertical, north-south and east-west components of broadband seismic station (GIGS) equipped with a Trillium 240 s seismometer; the black line is the hydraulic pressure expressed in MPa. To note the remarkable variation (0.2 Mpa pp) during the S-waves arrival. The time scale begins at 01:36:30 (UT) on 24 August 2016.
PHOTO (COLOR)
To investigate these rather unexpected phenomena in a quantitative fashion, we performed an elementary statistical analysis of the hydraulic pressure time series looking at the first three moments of the complete time-series. More specifically, after subtracting the (60 s) moving average we evaluated the standard deviation (m2)1/2, the skewness m3/(m2)3/2 and the kurtosis (m4/m22- 3) of the resulting distribution. Figures 13a,b show the results of the statistical analysis of hydraulic pressure data from S13 borehole. The standard deviation (Fig. 13a) is computed from 27 July (doy 209) to 4 September 2016 (doy 248). The Kurtosis and Skewness (Fig. 13b) values are computed from 15 June (doy 167) to 4 September 2016 (doy 248). The evidence of a signal change was very clear about five days before the main shock, although the change was somewhat detectable starting 40 days before the earthquake when a slow decrease in kurtosis is particularly evident. We checked the relevance of this signal by looking at the kurtosis of the full dataset (1 May 2015 to 5 September 2016) and we found no significant deviation from the zero value until 40 days before the main shock. It is worth to note that the 24 August mainshock was not preceded by any sizable foreshock; particularly, in the circular area of radius 10 km from the Mw 6.0 epicenter we counted only one event with Mw greater to 1.8 (15 Aug 2016, doy 228, at 19:00:41 UT - Mw 1.9) in the period 1 May - 24 August 201653.Statistical analysis: (a) standard deviation from 27 July (doy 209) to 4 September 2016 (doy 248) from hydraulic pressure data of S13 borehole; (b) kurtosis (blue) and skewness (red) values from 15 June (doy 167) to 4 September 2016 (doy 248) from hydraulic pressure data of S13 borehole. The statistical indicators are evaluated each 60 s. The plotted values are one-hour averages. Data missing (12-19 July 2016, doy 194-201) is due to instrumentation maintenance. The Amatrice earthquake (Mw 6.0) occurs on 24 August 2016 (doy 237) at 01:36:32 UT.
PHOTO (COLOR)
We next address the electrical conductivity data. Figure 14a shows the electrical conductivity signal averaged over 1 hour, starting 140 days before the Amatrice earthquake (occurring at the 0 in Fig. 14a) and lasting 20 days after. A significant change was noted starting 60 days before the mainshock, which can be correlated with the kurtosis analysis of hydraulic pressure (blue line in Fig. 13b). In the period from May 2015 to May 2016 we did not observe any significant variation of electrical conductivity; note again that no significant foreshocks or swarms were present in the area. Figure 14b shows a zoom of three hours (from 00:00 to 03:00 UT of 24 August 2016) of electrical conductivity data averaged every minute; a change in conductivity occurred at 01:51 UT, and the Amatrice earthquake strikes at 01:36. There were 15 minutes of lag time, but considering that the horizontal S13 tube contains about 2,000 l of water flowing at 2 l per minute, this was traced back, to the upper thrust zone, to around 17 hours before the mainshock.Electrical conductivity data of S13 borehole from 140 days before the Amatrice earthquake (24 August 2016 (doy 237) at 01:36:32 UT). The anomaly starts about 60 days before the mainshock in accordance with weak deviation from zero of Kurtosis of hydraulic pressure data in Fig. 13 (blue line); in the period from May 2015 to May 2016 we did not observe any significant variation. (b) Electrical conductivity data of Amatrice earthquake of 24 August 2016 (01:36:32 UT) from 00:00 to 03:00 UT.
PHOTO (COLOR)
During the period of data analysis and interpretation (between September and the middle of October 2016), which was taken a few months, we made further modifications and tests of our hydraulic system. At the end of October, the seismic sequence restarted with two events of magnitude greater than 5.0 (26 October 2016: 5.4 Mw at 17:10:36 UT and 5.9 Mw at 19:18:07 UT) culminating on October 30, 2016 with a 6.5 Mw (06:40:17 UT)53. During this period, only the pressure sensor was continuously recording, but with a different hydraulic port setup, so a direct comparison with older data is difficult. Figure 15a shows the 30 October earthquake (6.5 Mw at 06:40:17 UT) and Fig. 15b the sequence of 18 January 201753 (data from 06:00 to 18:00 UT); the main events are: 09:25:40 UT (Mw 5.1), 10:14:10 UT (Mw 5.5), 10:25:24 UT (Mw 5.4) and 13:33:37 UT (Mw 5.0). The amplitudes of four events with Mw > 5 are not directly correlated with magnitude because the distance between the earthquake epicentre and the S13 site decreases (about 26.6, 25.5, 25.1 and 24.6 km respectively); in addition, it is easy to notice the presence of a dozen earthquakes of smaller magnitude.(a) Plot of the hydraulic pressure (in MPa) of S13 borehole at 20 Hz of sampling rate for 30 October 2016 earthquake, Mw 6.5, at 06:40:17 UT. (b) Plot of the hydraulic pressure (in MPa) of S13 borehole at 20 Hz of sampling rate for the sequence of 18 January 2017 (data from 08:00 to 20:00 UT); the main events are: 09:25:40 UT (Mw 5.1), 10:14:10 UT (Mw 5.5), 10:25:24 UT (Mw 5.4) and 13:33:37 UT (Mw 5.0)53. It is also possible to note a small upward coseismic offset of about 0.002 MPa, as well as a dozen earthquakes of minor magnitude.
PHOTO (COLOR)
The only realistic comparison that can be made between these two periods concerns the coseismic effect following the two mainshocks of 24 August (6.0 Mw) and 30 October 2016 (6.5 Mw). The upward coseismic offset related to the two earthquakes is clearly shown by the superimposed plots of the hydraulic pressure (in MPa), with moving average at 1 minute (the average smoothes out the elastic part of earthquake signals, principally P-, S- and Coda- waves) (Fig. 16). The positive coseismic offset is almost the same (the data are aligned and superimposed forcing at (0,0) the start of two earthquakes) (Fig. 16).Superimposed plots of the hydraulic pressure of S13 borehole (in MPa), with moving average at 1 minute, of 24 August 2016, Mw 6.0 (red line) and 30 October 2016, Mw 6.5 (blue line) earthquakes (the average smooths the elastic part of earthquake signals). The upward coseismic offset is the same (data are shifted only in time and hydraulic pressure level).
PHOTO (COLOR)
Conclusions
The presence of the S13 horizontal borehole near the underground laboratories of Gran Sasso (LNGS-INFN) in central Italy, gave us the unique opportunity to investigate the pristine groundwater reservoir in the core of the Gran Sasso aquifer located in the very seismically active area of central Apennines (Fig. 1). The S13 horizontal borehole was drilled at the end of the 1980s during excavation works for a tunnel across the Gran Sasso chain, located about 39 km south-eastward from the Amatrice earthquake epicenter (Mw 6.0) (Fig. 3).
Starting in May 2015, we used the S13 borehole to measure continuously hydraulic pressure, groundwater temperature and electrical conductivity with a high sampling rate (10 Hz to 50 Hz) for each channel. Monitoring was still ongoing when the 24 August 2016 (01:36:32 UT, doy 237), Amatrice earthquake (Mw 6.0) hit about 39 km from the study site (Figs 1-3). The acquired record showed distinct and interesting signals in the hydraulic pressure and electrical conductivity data of the Gran Sasso groundwater before, during and after this earthquake. We recognized anomalies that can be classified both as long-term seismic precursors (anomalies of hydraulic pressure signals and electrical conductivity about 40-60 days before the Amatrice earthquake - Figs 13 and 14) and as short-term seismic precursors (about 5 days in the hydraulic pressure data - Figs 8, 9, 11 and 13).
For many years, scientists have observed changes in the amount, direction, and discharge of water flow in streams and in the water level of wells. Furthermore, changes in fluid pressure in the subsurface have been documented following large earthquakes by observing fluctuations as large as several meters in water levels resulting from mid-size regional earthquakes or even from very large earthquakes located thousands of kilometers away. The authors of ref.9 for example, have supplied a review of water level variations in surface wells coincident with the arrival of teleseismic waves. Our data, however, represent the first observations of earthquake-related signals in hydraulic pressure measured with continuous high-frequency sample rate (20 Hz) monitoring of the S13 borehole placed deep within a Gran Sasso carbonate aquifer; S13 also intersects a thrust fault near its end (Figs 2 and 3).
Groundwater variations are believed to respond to crustal deformation processes well before the occurrence of significant earthquakes20,21. Several mechanisms have been proposed to explain such variations during the seismic cycle, including the dilatancy theory which assumes that, during the pre-seismic phase, the rock volume increase, prior to rupture, is favoured by the micro-fracturing within the earthquake source volume. Micro-fracturing induces also rock permeability increase favouring groundwater variations and mobilization or uprising of geogas9,50,85.
Groundwater monitoring at the S13 borehole revealed the existence of highly dynamic aquifer behaviour. The data collected during these twenty-one months turned out to be quite a useful tool to: (i) characterize the background in groundwater monitoring, and (ii) identify more clearly the observed anomalies and their possible link with phenomena of tectonic interest. In this respect, the anomalies observed in the hydraulic pressure 5 days before the 24 August 2016, Amatrice earthquake are obvious.
We interpret the presence of negative micropulses in hydraulic pressure data (Fig. 10b) as the result of a possible uprising of deep geogas, mainly CO2, consistent with recent studies suggesting a link between fault mechanics, seismicity, underground fluids dynamics and CO2 uprising86-94. More precisely, during the L'Aquila earthquake (April 6th, 2009, 6.3 Mw), the evolution of seismicity was being driven in part by the poro-elastic response of trapped reservoirs of high-pressure fluid, presumably CO2, and postseismic fluid flow initiated by the main shock44.
Further investigation of the relationship between earthquakes and changes in groundwater parameters near large seismogenic faults, are needed for a full understanding of preseismic, coseismic and postseismic processes.
Acknowledgements
The authors thank Prof. Stefano Ragazzi, Director of the Laboratorio Nazionale del Gran Sasso of Istituto Nazionale di Fisica Nucleare (LNGS-INFN), for providing facilities and assistence during the field experiment, Dr. Gianluca Valensise (INGV) and Dr. Bernard Chouet for help in improving the manuscript. Prof. C. Doglioni, President of Istituto Nazionale di Geofisica e Vulcanologia (INGV) and Prof. F. Ferroni, President of Istituto Nazionale di Fisica Nucleare (INFN), are also acknowledged for the continuous support and encouragements. GDC acknowledges enlightening discussions with D. Benedetto and A. Galante. GDL acknowledges M. Mariotti for technical help on the instrumentation. Many thanks to Kerry Jane Rhoden for the careful revision of the English language. Finally, a strong thanks to the Editor for encouraging and the anonymous referees which have considerably improved the manuscript with their suggestions.
Author Contributions
All authors had the idea within the G.D.L.'s PhD project (supervisor: M.T.); G.D.L. and G.D.C. performed the experimental work and the analyses; M.T. created Figures 2 and 3; G.D.L. and G.D.C. created Figures 1 and 4-16; M.T. wrote the hydrogeological setting of the study area; G.D.L. and G.D.C. wrote the main manuscript text; all authors discussed the results.
Competing Interests
The authors declare no competing interests.
References
[
References
1
Cicerone RD, Ebel JE, Britton J, A systematic compilation of earthquake precursors, Tectonophysics, 2009, 476, 371, 3962009Tectp.476..371C, 10.1016/j.tecto.2009.06.008
2
Wyss M, Cannot earthquakes be predicted?, Science, 1997, 278, 487, 4881997Sci...278..487W, 10.1126/science.278.5337.487
3
Scholz CH, Sykes LR, Aggarwal YP, Earthquake Prediction: A Physical Basis, Science, 1973, 181, 803, 8101973Sci...181..803S, 10.1126/science.181.4102.803
4
Wyss M, Booth DC, The IASPEI procedure for the evaluation of earthquake precursors, Geophys. J. Int., 1997, 131, 423, 4241997GeoJI.131..423W, 10.1111/j.1365-246X.1997.tb06587.x
5
Geller RJ, Earthquake prediction: a critical review, Geophys. J. Int., 1997, 131, 425, 4501997GeoJI.131..425G, 10.1111/j.1365-246X.1997.tb06588.x
6
Kagan YY, Are earthquakes predictable?, Geophys. J. Int., 1997, 131, 505, 5251997GeoJI.131..505K, 10.1111/j.1365-246X.1997.tb06595.x
7
Nature debates, http://www.nature.com/nature/debates/earthquake/index.html.
8
Geller RJ, Jackson DD, Kagan YY, Mulargia F, Earthquakes cannot be predicted, Science, 1997, 275, 1616, 1617, 10.1126/science.275.5306.1616
9
Manga M., Wang C.-Y., Earthquake Hydrology, Treatise on Geophysics, 2015, 305, 328
Roeloffs EA, Hydrologic precursors to earthquakes: A review, Pure and Applied Geophysics, 1988, 126, 177, 2091988PApGe.126..177R, 10.1007/BF00878996
Kissin IG, Grinevsky AO, Main features of hydrogeodynamic earthquake precursors, Tectonophysics, 1990, 178, 277, 2861990Tectp.178..277K, 10.1016/0040-1951(90)90154-Z
Tsunogai U, Wakita H, Precursory Chemical Changes in Ground Water: Kobe Earthquake, Japan, Science, 1995, 269, 61, 631995Sci...269...61T, 10.1126/science.269.5220.61
Shi Z, Wang G, Manga M, Wang C-Y, Continental-scale water-level response to a large earthquake, Geofluids, 2015, 15, 310, 320, 10.1111/gfl.12099
Ingebritsen SE, Manga M, Earthquakes: Hydrogeochemical precursors, Nature Geoscience, 2014, 7, 697, 6982014NatGe...7..697I, 10.1038/ngeo2261
Claesson L, Hydrogeochemical changes before and after a major earthquake, Geology, 2004, 32, 641, 6442004Geo....32..641C, 10.1130/G20542.1
Chen C-H, Anomalous frequency characteristics of groundwater level before major earthquakes in Taiwan, Hydrology and Earth System Sciences, 2013, 17, 1693, 17032013HESS...17.1693C, 10.5194/hess-17-1693-2013
Skelton A, Changes in groundwater chemistry before two consecutive earthquakes in Iceland, Nature Geoscience, 2014, 7, 752, 7562014NatGe...7..752S, 10.1038/NGEO2250
Roeloffs E, Quilty E, Case 21: Water level and strain changes preceding and following the August 4, 1985 Kettleman Hills, California, earthquake, Pure and Applied Geophysics, 1997, 149, 21, 601997PApGe.149...21R, 10.1007/BF00945160
Roeloffs EA, Persistent water level changes in a well near Parkfield, California, due to local and distant earthquakes, J. Geophys. Res., 1998, 103, 869, 8891998JGR...103..869R, 10.1029/97JB02335
King CY, In search of earthquake precursors in the water-level data of 16 closely clustered wells at Tono, Japan, Geophys. J. Int., 2000, 143, 469, 4772000GeoJI.143..469K, 10.1046/j.1365-246X.2000.01272.x
King CY, Earthquake related water-level changes at 16 closely clustered wells in Tono, Central Japan, J. Geophys. Res., 1999, 104, 13073, 130821999JGR...10413073K, 10.1029/1999JB900080
King CY, Koizumi N, Kitagawa Y, Hydrogeochemical anomalies and the 1995 Kobe earthquake, Science, 1995, 269, 38, 391995Sci...269...38K, 10.1126/science.269.5220.38
Chia Y, Implications of coseismic groundwater level changes observed at multiple-well monitoring stations, Geophys. J. Int., 2008, 172, 293, 3012008GeoJI.172..293C, 10.1111/j.1365-246X.2007.03628.x
Chia Y, Chiu JJ, Chiang Y-H, Lee T-P, Liu C-W, Spatial and temporal changes of groundwater level induced by thrust faulting, Pure and Applied Geophysics, 2008, 165, 5, 162008PApGe.165....5C, 10.1007/s00024-007-0293-5
Montgomery DR, Manga M, Stream flow and water well responses to earthquakes, Science, 2003, 300, 2047, 20492003Sci...300.2047M, 10.1126/science.1082980
Muir-Wood R, King GCP, Hydrological Signatures of Earthquake Strain, J. Geophys. Res., 1993, 98, B12, 22035, 220681993JGR....9822035M, 10.1029/93JB02219
Brodsky EE, Roeloffs E, Woodcock D, Gall I, Manga M, A mechanism for sustained groundwater pressure changes induced by distant earthquakes, J. Geophys. Res., 2003, 108, 23902003JGRB..108.2390B, 10.1029/2002JB002321
Wakita H, Water wells as possible indicators of tectonic strain, Science, 1975, 189, 553, 5551975Sci...189..553W, 10.1126/science.189.4202.553
Doglioni C, Barba S, Carminati E, Riguzzi F, Fault on-off versus coseismic fluids reaction, Geoscience Frontiers, 2014, 5, 767, 780, 10.1016/j.gsf.2013.08.004
Chiarabba C, Chiodini G, Continental delamination and mantle dynamics drive topography, extension and fluid discharge in the Apennines, Geology, 2013, 41, 6, 715, 7182013Geo....41..715C, 10.1130/G33992.1
McCaig A, Fluid flow through fault zones, Nature, 1989, 340, 6001989Natur.340..600M, 10.1038/340600a0
Barnhoorn A, Cox SF, Robinson DJ, Senden T, Stress- and fluid-driven failure during fracture array growth: Implications for coupled deformation and fluid flow in the crust, Geology, 2010, 38, 779, 7822010Geo....38..779B, 10.1130/G31010.1
Di Toro G, Fault lubrication during earthquakes, Nature, 2011, 471, 494, 4992011Natur.471..494D, 10.1038/nature09838
Violay M, Effect of water on the frictional behaviour of cohesive rocks during earthquakes, Geology, 2013, 42, 27, 302014Geo....42...27V, 10.1130/G34916.1
Hirth G, Beeler NM, The role of fluid pressure on frictional behaviour at the base of the seismogenic zone, Geology, 2015, 43, 223, 2262015Geo....43..223H, 10.1130/G36361.1
Sibson RH, Fluid involvement in normal faulting, J. Geodynamics, 2000, 29, 469, 4992000JGeo...29..469S, 10.1016/S0264-3707(99)00042-3
Scuderia MM, Collettinia C, Marone C, Frictional stability and earthquake triggering during fluid pressure stimulation of an experimental fault, Earth Planetary Science Letters, 2017, 477, 84, 962017E&PSL.477...84S, 10.1016/j.epsl.2017.08.009
Chiarabba C, The2009 L'Aquila (central Italy) Mw 6.3 earthquake: Main shock and aftershocks, Geophys. Res. Lett., 2009, 36, L183082009GeoRL..3618308C, 10.1029/2009GL039627
Valoroso L, Radiography of a normal fault system by 64,000 high-precision earthquake locations: The 2009 L'Aquila (central Italy) case study, J. Geophys. Res., 2010, 118, 1156, 1176, 10.1002/jgrb.50130
Anzidei, M. & Pondrelli, S. (eds). The Amatrice seismic sequence: preliminary data and results., Annals of Geophysics, Special Issue 59, Fast Track5, http://www.annalsofgeophysics.eu/index.php/annals/issue/view/515 (2016).
Chiaraluce L, The2016 Central Italy Seismic Sequence: A First Look at the Mainshocks, Aftershocks, and Source Models, Seismological Research Letters, 2017, 88, 757, 771, 10.1785/0220160221
De Luca G, La rete sismica regionale Abruzzo e sua integrazione con la RSN. Riassunti estesi del I° Workshop Tecnico Monitoraggio sismico del territorio nazionale: stato dell'arte e sviluppo delle reti di monitoraggio sismico. AA. VV. a cura di Cattaneo M. & Moretti M. - Miscellanea INGV, Roma 20-21 dicembre, 2011, 2010, 22, 23
Delladio A, Monitoraggio sismico del territorio nazionale. Riassunti estesi del I° Workshop Tecnico Monitoraggio sismico del territorio nazionale: stato dell'arte e sviluppo delle reti di monitoraggio sismico. AA. VV. a cura di Cattaneo M. & Moretti M. - Miscellanea INGV, Roma 20-21 dicembre, 2011, 2010, 11, 16
Terakawa T, Zoporowski A, Galvan B, Miller SA, High-pressure fluid at hypocentral depths in the L'Aquila region inferred from earthquake focal mechanisms, Geology, 2010, 38, 995, 9982010Geo....38..995T, 10.1130/G31457.1
Savage MK, The role of fluids in earthquake generation in the 2009 M, w 6.3 L'Aquila, Italy, earthquake and its foreshocks, Geology, 2010, 38, 1055, 10562010Geo....38.1055S, 10.1130/focus112010.1
Di Luccio F, Ventura G, Di Giovambattista R, Piscini A, Cinti FR, Normal faults and thrusts reactivated by deep fluids: The 6 April 2009 M, w 6.3 L'Aquila earthquake, central Italy, J. Geophys. Res., 2010, 115, B063152010JGRB..115.6315D, 10.1029/2009JB007190
Malagnini L, Lucente FP, De Gori P, Akinci A, Munafo I, Control of pore fluid pressure diffusion on fault failure mode: Insights from the 2009 L'Aquila seismic sequence, J. Geophys. Res., 2012, 117, B053022012JGRB..117.5302M, 10.1029/2011JB008911
Lucente FP, Temporal variation of seismic velocity and anisotropy before the 2009 M, w 6.3 L'Aquila earthquake, Italy, Geology, 2010, 38, 1015, 10182010Geo....38.1015P, 10.1130/G31463.1
Borghi A, Aoudia A, Javed F, Barzaghi R, Precursory slow-slip loaded the 2009 L'Aquila earthquake sequence, Geophys. J. Int., 2016, 205, 776, 7842016GeoJI.205..776B, 10.1093/gji/ggw046
Moro M, New insights into earthquake precursors from InSAR, Scientific Reports, 2017, 72017NatSR...712035M, 10.1038/s41598-017-12058-35607348
Barberio MD, Barbieri M, Billi A, Doglioni C, Petitta M, Hydrogeochemical changes before and during the 2016 Amatrice-Norcia seismic sequence (central Italy), Scientific Reports, 2017, 72017NatSR...711735B, 10.1038/s41598-017-11990-85601465
Petitta M, Water table and discharge changes associated with the 2016-2017 seismic sequence in central Italy: hydrogeological data and conceptual model for fractured carbonate aquifers, Hydrogeology Journal, 2017, 26, 1009, 10262018HydJ...26.1009P, 10.1007/s10040-017-1717-7
ISIDe working group - version 1.0, http://cnt.rm.ingv.it/iside, 10.13127/ISIDe (2016).
Catalano, P. G. & Salza, R. Laboratori sotterranei. Relazione idrogeologica e sistemi di drenaggio e canalizzazione. [INFN underground laboratories. Hydrogeology and the groundwater drainage network report]., INFN report. 180 pp (2003).
Catalano PG, Cavinato G, Salvini F, Tozzi M, Analisi strutturale nei Laboratori dell'INFN del Gran Sasso d'Italia, Mem. Soc. Geol. It., 1986, 35, 647, 655
Rovida, A., Locati, M., Camassi, R., Lolli, B. & Gasperini, P. (eds). CPTI15, the 2015 version of the Parametric Catalogue of Italian Earthquakes. Istituto Nazionale di Geofisica e Vulcanologia; doi: http://doi.org/10.6092/INGV.IT-CPTI15 (2016).
Valensise G, Pantosti D, The investigation of potential earthquake sources in peninsular Italy: A review, J. Seismol., 2001, 5, 287, 3062001JSeis...5..287V, 10.1023/A:1011463223440
Chiarabba C, Jovane L, Di Stefano R, A new view of Italian seismicity using 20 years of instrumental recordings, Tectonophysics, 2005, 395, 251, 2682005Tectp.395..251C, 10.1016/j.tecto.2004.09.013
De Luca G, Cattaneo M, Monachesi G, Amato A, Seismicity in the Umbria-Marche region integrating the Italian national and regional networks, Tectonophysics, 2009, 476, 121, 1352009Tectp.476..121D, 10.1016/j.tecto.2008.11.032
De Luca G, A detailed analysis of two seismic sequences in Abruzzo, Central Apennines, Italy, J. Seismol., 2000, 4, 1, 21, 10.1023/A:1009870327083
Frepoli A, Seismic sequences and swarms in the Latium-Abruzzo-Molise Apennines (central Italy): New observations and analysis from a dense monitoring of the recent activity, Tectonophysics, 2017, 712-713, 312, 3292017Tectp.712..312F, 10.1016/j.tecto.2017.05.026
Bagh S, Background seismicity in the central Apennines of Italy: The Abruzzo region case study, Tectonophysics, 2007, 444, 80, 922007Tectp.444...80B, 10.1016/j.tecto.2007.08.009
Galli P, Galadini F, Moro M, Giraudi C, New paleoseismological data from the Gran Sasso d'Italia area (central Apennines), Geophys. Res. Lett., 2002, 29, 7, 1134, 11372002GeoRL..29.1134G, 10.1029/2001GL013292
Amoruso A, Crescentini L, Petitta M, Rusi S, Tallini M, Impact of the April 6, 2009 L'Aquila earthquake on groundwater flow in the Gran Sasso carbonate aquifer, Central Italy, Hydrol. Process., 2011, 25, 1754, 17642011HyPr...25.1754A, 10.1002/hyp.7933
Adinolfi Falcone R, Changes on groundwater flow and hydrochemistry of the Gran Sasso carbonate aquifer due to the 2009 L'Aquila earthquake, Ital. J. Geosci. (Boll. Soc. Geol. It.), 2012, 131, 459, 474, 10.3301/IJG.2011.34
Chiodini G, Geochemical evidence for and characterization of CO, 2 rich gas sources in the epicentral area of the Abruzzo 2009 earthquakes, Earth and Planetary Science Letters, 2011, 304, 389, 3982011E&PSL.304..389C, 10.1016/j.epsl.2011.02.016
Cigolini C, Laiolo M, Coppola D, The LVD signals during the early-mid stages of the L'Aquila seismic sequence and the radon signature of some aftershocks of moderate magnitude, J. Environ. Radioactivity, 2015, 139, 56, 65, 10.1016/j.jenvrad.2014.09.017
Plastino W, Uranium groundwater anomalies and L'Aquila earthquake, 6, th April 2009 (Italy), J. Environ. Radioactivity, 2010, 101, 45, 50, 10.1016/j.jenvrad.2009.08.009
Barbieri M, Boschetti T, Petitta M, Tallini M, Stable isotopes (, 2H, 18O and, 87Sr/, 86Sr) and hydrochemistry monitoring for groundwater hydrodinamics analysis in a karst aquifer (Gran Sasso, Central Italy), Appl. Geochem., 2005, 20, 2063, 2081, 10.1016/j.apgeochem.2005.07.008
Petitta M, Tallini M, Idrodinamica sotterranea del massiccio del Gran Sasso (Abruzzo): indagini idrologiche, idrogeologiche e idrochimiche (1994-2001), Boll. Soc. Geol. D'It, 2002, 121, 3, 343, 363
Amoruso A, Crescentini L, Petitta M, Tallini M, Parsimonious recharge/discharge modeling in carbonate fractured aquifers: the groundwater flow in the Gran Sasso aquifer (Central Italy), Journal of Hydrology, 2012, 476, 136, 1462013JHyd..476..136A, 10.1016/j.jhydrol.2012.10.026
Amoruso A, Crescentini L, Martino S, Petitta M, Tallini M, Correlation between groundwater flow and deformation in the fractured carbonate Gran Sasso aquifer (INFN underground laboratories, Central Italy), Water Resour. Res., 2014, 50, 1, 19, 10.1002/2013WR014491
Tallini M, Parisse B, Petitta M, Spizzico M, Long-term spatio-temporal hydrochemical and, 222Rn tracing to investigate groundwater flow and water-rock interaction in the Gran Sasso (central Italy) carbonate aquifer, Hydrogeology Journal, 2013, 21, 7, 1447, 14672013HydJ...21.1447T, 10.1007/s10040-013-1023-y
Tallini M, Isotope hydrology and geochemical modeling: new insights into recharge process and water-rock interaction in the Gran Sasso fissured carbonate aquifer (Central Italy), Environ. Earth Sci., 2014, 72, 12, 4957, 4971, 10.1007/s12665-014-3364-9
Fiorillo F, Long term trend and fluctuations of karst spring discharge in a Mediterranean area (Central-Southern Italy), Environ. Earth Sci., 2015, 74, 153, 172, 10.1007/s12665-014-3946-6
Scozzafava M, Tallini M, Report: Net infiltration in the Gran Sasso Massif of central Italy using the Thornthwaite water budget and curve-number method, Hydrogeology Journal, 2001, 9, 461, 4752001HydJ....9..461S, 10.1007/s100400100151
Adinolfi Falcone R, Falgiani A, Petitta M, Tallini M, Characteristics of Gran Sasso INFN laboratory groundwater (inferred from 1996-1998 spot sampling data) to fine-tune the conceptual model of water-rock interaction in carbonate aquifer, Internal Report INFN, LNGS/GEO, 2006, 01, 42 pp
Adinolfi Falcone R, Chemical and isotopic (δ, 18O‰, δ, 2H‰, δ, 13C‰, 222Rn) multi-tracing for groundwater conceptual model of carbonate aquifer (Gran Sasso INFN underground laboratory - central Italy), Journal of Hydrology, 2008, 357, 368, 3882008JHyd..357..368A, 10.1016/j.jhydrol.2008.05.016
D'Agostino N, Chamot-Rooke N, Funiciello R, Jolivet L, Speranza F, The role of pre-existing thrust faults and topography on the styles of extension in the Gran Sasso range (central Italy), Tectonophysics, 1998, 292, 229, 2541998Tectp.292..229D, 10.1016/S0040-1951(98)00070-5
De Luca, G., Di Carlo, G. & Tallini, M. Hydraulic pressure variations of groundwater in the Gran Sasso underground laboratory during the Amatrice earthquake of August 24, 2016. Annals of Geophysics 59, Fast Track 5; doi: 10.4401/AG-7200 (2016).
Di Virgilio ADV, GINGER: A feasibility study, Eur. Phys. J. Plus, 2017, 132, 10.1140/epjp/i2017-11452-6
Simonelli A, First deep underground observation of rotational signals from an earthquake at teleseismic distance using a large ring laser gyroscope, Annals of Geophysics 59, Fast Track, 2016, 4, 1, 6, 10.4401/ag-6970
Simonelli A, Rotational motions from the 2016, central italy seismic sequence, as observed by an underground ring laser gyroscope, Geophys. J. Int., 2018, 214, 705, 7152018GeoJI.214..705S, 10.1093/gji/ggy186
Kano Y, Yanagidani T, Broadband hydroseismograms observed by closed borehole wells in the Kamioka mine, central Japan: Response of pore pressure to seismic waves from 0.05 to 2 Hz, J. Geophys. Res., 2006, 111, B034102006JGRB..111.3410K, 10.1029/2005JB003656
Jónsson S, Segall P, Pedersen R, Björnsson G, Post-earthquake ground movements correlated to pore-pressure transients, Nature, 2003, 424, 179, 1832003Natur.424..179J, 10.1038/nature01776
Chiodini G, Carbon dioxide Earth degassing and seismogenesis in central and southern Italy, Geophys. Res. Lett., 2004, 31, L076152004GeoRL..31.7615C, 10.1029/2004GL019480
Doglioni C, Carminati E, Petricca P, Riguzzi F, Normal fault earthquakes or graviquakes, Scientific Reports, 2015, 52015NatSR...512110D, 10.1038/srep121104501010
Miller AM, Aftershocks driven by a high-pressure CO, 2 source at depth, Nature, 2004, 427, 724, 7272004Natur.427..724M, 10.1038/nature02251
Crews JB, Cooper CA, Experimental evidence for seismically initiated gas bubble nucleation and growth in groundwater as a mechanism for coseismic borehole water level rise and remotely triggered seismicity, J. Geophys. Res. Solid Earth, 2014, 119, 7079, 70912014JGRB..119.7079C, 10.1002/2014JB011398
Toutain J-P, Baubron J-C, Gas geochemistry and seismotectonics: a review, Tectonophysics, 1999, 304, 1, 271999Tectp.304....1T, 10.1016/S0040-1951(98)00295-9
King C-Y, Zhang W, Zhang Z, Earthquake-induced groundwater and gas changes, Pure Appl. Geophys., 2006, 163, 633, 6452006PApGe.163..633K, 10.1007/s00024-006-0049-7
Di Luccio F, Seismic signature of active intrusions in mountain chains, Sci. Adv., 2018, 4, 1, 10.1126/sciadv.1701825
Girault F, Persistent CO, 2 emissions and hydrothermal unrest following the 2015 earthquake in Nepal, Nat. Commun., 2018, 92018NatCo...9.2956G, 10.1038/s41467-018-05138-z6063904
Hotovec-Ellis AJ, Deep fluid pathways beneath Mammoth Mountain, California, illuminated by migrating earthquake swarms. Sci, Adv., 2018, 4, eaat5258, 10.1126/sciadv.aat5258
]
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
By Gaetano De Luca; Giuseppe Di Carlo and Marco Tallini
