Abstract

Many Italian rock slopes are characterized by unstable rock masses that cause constant rock falls and rockslides. To effectively mitigate their catastrophic consequence thorough studies are required. Four velocimeters have been placed in the Torgiovannetto quarry area for an extensive seismic noise investigation. The study area (with an approximate surface of 200×100 m) is located near the town of Assisi (Italy) and is threatened by a rockslide. In this work, we present the results of the preliminary horizontal to vertical spectral ratio analysis of the acquired passive seismic data aimed at understanding the pattern of the seismic noise variation in case of stress state and/or weathering conditions (fluid content and microfracturing). The Torgiovannetto unstable slope has been monitored since 2003 by Alta Scuola of Perugia and the Department of Earth Sciences of the University of Firenze, after the observation of a first movement by the State Forestry Corps. The available data allowed an extensive comparison between seismic signals, displacement, and meteorological information. The measured displacements are well correlated with the precipitation trend, but unfortunately no resemblance with the seismic data was observed. However, a significant correlation between temperature data and the horizontal to vertical spectral ratio trend of the seismic noise could be identified. This can be related to the indirect effect of temperature on rock mass conditions and further extensive studies (also in the time frequency domain) are required to better comprehend this dependency. Finally, the continuous on-line data reveal interesting applications to provide near-real time warning systems for emerging potentially disastrous rockslides.

1. Introduction

For many years, researchers have turned their attention to the massive problem of landslides in Italy. The topic is high on the agenda because roughly 70% of all the landslides in the European continent are concentrated in Italy [1]. As a consequence of steep slopes, high seismic activity, and soil and bedrock properties, many hillsides of the Italian valleys are characterized by unstable rock masses causing constant rock falls and rockslides of various sizes and types [2]. A thorough understanding of failure types, mechanisms, and possible causes of landslides is required to effectively mitigate their catastrophic consequences. Moreover, currently early warning systems (EWS) can be implemented in order to prevent loss of life and to reduce the economic and material impact of landslide events [3, 4]. Nevertheless, frequently enough, it is not easy to find a technique able to provide an immediate alert [5]. Therefore, slope failure of rock masses represents an interesting case study for verifying the feasibility of using passive seismic monitoring in EWS. By means of the observation of the changes which occurred in the acquired signal, in fact, it could be possible to detect variations in the elastic parameters of the rock body related to changes in pore-fluid pressure, consolidation, and microfracturing that could forecast failure [6].

In the last years, besides the traditional geotechnical and structural monitoring (e.g., topographic total stations, extensometers, and inclinometers, [7]), new techniques have been used to characterize and monitor landslides: aerial photos and LiDAR [8, 9], GPS monitoring [1012], InSAR and GB-InSAR technique [3, 1315], laser scanner [16, 17], infrared thermography [1820], and optic fiber strain sensors [21]. Shallow geophysical methods represent a valid complement to the aforementioned techniques [6, 2227].

To verify the performance of a small-scale seismic network as part of an EWS, a pilot scale experiment was arranged to monitor an unstable rock mass. The test site is the Torgiovannetto quarry located in Umbria Region, one of the Italian Regions that is more prone to landslide. In general, quarries can be characterized as remarkably vulnerable areas, since their natural geomorphology is altered by excavating activities [28]. The data were collected during a 7-month-period monitoring. In this paper we present the results of the preliminary analysis carried out on the acquired data by means of the horizontal to vertical spectral ratio analysis (HVSR or H/V). These analyses were aimed at understanding the pattern of seismic noise variation in case of stress state and/or weathering conditions (fluid content and microfracturing) as the first step to set up a reliable EWS. The studied quarry rockslide was also extensively monitored since 2003 with traditional methods. Therefore, the multiparameter analysis was useful to understand the mechanisms that control the rockslide dynamics and to evaluate possible connection between rainfall/temperature/displacement and rockslide seismic activity. Thus, a comparison between the seismic data and both temperature and precipitation data is discussed, in order to highlight a correlation between them.

2. The Study Area

The Torgiovannetto test site is located in a micritic limestone former quarry (dismissed since the late '90s), 2 km NE from Assisi (Umbria Region in Central Italy) in the northward facing slope of Mount Subasio (red square in Figure 1). Landslides in Umbria occupy about 14% of the entire land cover (8456 km2) and affect many urban areas.

Mount Subasio (1109 m a.s.l.) is part of the Umbria-Marche Apennines, a complex fold and arcuate thrust belt that occupies the outer zones of the Northern Apennines of Italy. The belt developed during the Neogene as a result of the Ligurian Ocean closure, followed by the continental collision between the European Corsica-Sardinia Margin and the African Adria Promontory [29]. A northeast-directed compressional tectonic phase started during the middle Miocene and is still active near the Adriatic coast [30]. During the upper Pliocene an extensional phase started with a principal stress oriented about NE-SW that resulted in the dissection of the Umbria-Marche Apennines and the opening of a NW-SE-trending set of continental basins. Mount Subasio area consists in a SSE-NNW trending anticline [31, 32] with layers dipping almost vertically in the NE side of the mountain with several NW-SE striking normal faults on the eastern and western flanks. The local geological formations, belonging to the Umbro-Marchigiana Sequence (from Calcare Massiccio to Marnoso Arenacea), represent the progressive sinking of a marine environment.

The study area consists mainly of micritic limestone belonging to the Maiolica Formation (Upper Jurassic-Lower Cretaceous) that widely outcrops in the area. The thickness of the Formation is about 100 m and is composed by white or light grey well-stratified micritic limestone layers, whose thickness ranges between 10 cm and 1 m, and thin clay interlayers may sporadically occur. The site is also partially covered by very heterometric debris (from pebble- to cobble-sized angular clasts, with scattered boulders, in a silty or coarse grained sandy matrix), some of which are of anthropogenic nature. The dip direction varies between 350° and 5°, while the dip of layers from 25° to 35°, which means that, in general, the layers’ dip is in the same direction of the slope but with a gentler angle.

First deformations within the quarry site were observed in May 2003 by the State Forestry Corps, in the form of tension cracks in the vegetated area above and within the quarry front. From then, several monitoring campaigns were carried out by means of different techniques (topographic total station, inclinometers, extensometers, ground-based interferometric radar, laser scanner, and infrared thermal camera [7]). It is assured that the main predisposing factor of instability was the quarrying activity that heavily altered the original front. Actually the quarry is structured in four main terraces, that the dense vegetation prevents distinguish well them (refer to [4] for the quarry view from north to south), with an overall height of about 140 m (Figure 2). Nevertheless, earthquakes-induced landslides cannot be neglected among the instability factors. In fact, the link between earthquake and landslide is well documented in the literature, especially in the cases of high-magnitude seismic event [3336]. For example, the seismic sequence that affected the area southeast of the quarry (Colfiorito basin) in the 1997-98 reached the Assisi area in a macroseismic intensity (MCS) Io = 8-9 [37]. Therefore, the seismicity of the area surrounding the quarry is another important instability factor.

The main rockslide [38] in the Torgiovannetto quarry has a rough trapezoidal shape and covers about 200 m x 100 m in surface and 550 m a.s.l. and 680 m a.s.l. in altitude. The geometry and other soil parameters (such as densities and body wave velocities) are well known thanks to the geotechnical and geophysical investigations carried out on the site by Alta Scuola of Perugia and by the University of Firenze [4, 7]. Among these investigations, a passive seismic network in continuous recording was installed on this rockslide from December 2012 to July 2013. The “traditional” monitoring network was composed by 13 wire extensometers, 1 accelerometer, 1 meteorological station (composed of 1 thermometer and 1 rain gauge), and 3 inclinometers (Figure 2). The monitoring network, progressively enhanced and improved throughout the years, was completed with hydrological data [39], modelling computation analysis [7, 40], and the seismological stations. Nowadays, the active volume of Torgiovannetto rockslide is estimated to be about 182,000 m3. The upper boundary is defined by a big open subvertical fracture (Figure 2), a tension crack with an EW strike, which in some places displays a width up to 2 m and depth of about 20 meters [40].

3. Methods

The HVSR technique was introduced for the first time by [41, 42]. It is based on the ratio between the horizontal and vertical components of ground motion and it requires a 3-component sensor to acquire data. According to [43] microtremor energy consists mainly of SH waves, while, according to other authors, as discussed in [44], H/V peaks are related to Rayleigh waves. One of the striking features of the HVSR ratio is its stability in time, documented in many papers [27, 45, 46]. The HVSR curve allows gaining additional information about the underlying velocity profile at the site, especially when a strong different shear wave velocity exists between the shallow layer and the bedrock [47, 48]. The site effect amplification, in fact, could be caused by several geological conditions and one of them is the presence of a soft soil layer overlying a rigid half space. Nowadays, the HVSR is widely used both for environmental [49, 50] and for structural [5153] problems. For a more detailed discussion about the seismic noise method please refer to the wide literature [27, 43, 5456]. The main application of HVSR technique on landslide concerns the possibility to reconstruct the geometry of the sliding mass and to detect the depth of the shear surface [27, 5759] with a good approximation. This point is beyond the scope of this paper that, instead, aimed at the evaluation of the dynamic behaviour of the rock mass affected by the presence of fractures linked to the sliding wedge, searching for changes in its internal characteristics detectable by the HVSR shape [6, 23] that could be used in early-warning procedure.

At the Torgiovannetto quarry, seismic measurements were performed using a small-scale network composed by four seismic stations (TOR1, TOR2, TOR3, and TOR4; locations are shown in Figure 2). Due to the geomorphological characteristics of the site and the lack of access to the eastern part of the slope, the installation was really challenging. Station TOR4 was located on the sliding mass while the other three stations (TOR1, TOR2, and TOR3) were located at the edge of the quarry arranged in pairs with diametrically opposite position with respect to the centre of the landslide. This configuration (a reverse Y with respect to the sliding orientation) allowed us to retrieve punctual information both inside and outside the landslide. Each station with a SARA 24bit A/D converter (SL06) coupled with a SS45 tri-axial velocimeter sensor with a natural frequency of 4.5 Hz and transduction factor of 78 V/m/s. Instruments response is flat down to 2 Hz, with an upper-corner frequency of 100 Hz. All of them were equipped with Global Positioning System (GPS) receivers for time synchronization. The sensors were placed on a concrete base with supporting plinth, isolated from the exterior in order to attain protection from severe weather conditions. Battery supply and digitizer, connected to the sensors through a connector cable, were housed in a separate case. Data were recorded in continuous mode at 200 Hz sampling frequency, as the best compromise between signal resolution and data storage. Data acquisition was continuous for 210 days from December 7, 2012, to July 3, 2013, except for some short intervals due to the batteries change. Data format of the seismic records retrieved from the converters SL06 is miniSEED (‘Data-only’ volume; http://ds.iris.edu/ds/nodes/dmc/data/formats/miniseed/). Nevertheless, this format was mainly designed for the exchange of geophysical data and not for analysis. Therefore, first of all, recorded data were converted into a more suitable format for elaborations like SAC (Seismological Analysis Code; https://ds.iris.edu/files/sac-manual/manual/file_format.html). For each station, every 6 hours, three separate files were generated (Figure 3), which correspond to the east-west (SHE), north-south (SHN), and vertical or up-down (SHZ) components of ground velocity. The amplitude (y-axis) was expressed in counts, while the x-axis in time (hours).

Data analysis was performed by means of Geopsy software (http://www.geopsy.org; cf. [53] as an example of application). For all the 3 components of ground motion the acquired data were detrended, mean-removed, and filtered. Then, each trace was divided into windows of 120 s length, and each window was tapered with a Tukey window and padded with zeros. The amplitude spectrum was evaluated via the Fast Fourier Transform (FFT); individual spectra were finally smoothed using a boxcar of 0.1 Hz width. The H/V ratio was calculated for each window, and the final HVSR function was given by the average of the HVSRs over 6 h intervals. In this work the horizontal (H) spectra have been computed by averaging E-W and N-S components using a quadratic mean, which shows a lower bias with respect to the simple arithmetic mean [52]. Finally, a special filtering process was not applied since it did not significantly affect nonstationary noises as happened in other studies [6062].

4. Results

The resonance frequency peaks, determined using the H/V method described above, were analysed for stations TOR1, TOR3, and TOR4 throughout the whole monitored period. The HVSR analysis of station TOR2 is not presented here because of the typical flat shape of the outcropping seismic bedrock [55]. The TOR1 HVSR (Figure 4) exhibits the highest amplitudes over the frequency band of 4.5 Hz to 13 Hz with a stable peak around 10.5 Hz whose amplitude is generally slightly above 2. Throughout the observation period, both, peak amplitude and frequency, did not exhibit any particular trend, with the exception of a slight increase in amplitude within the period of January to mid-March, 2013. The TOR3 HVSR (Figure 5) is characterized by more closely spaced peaks of amplitude higher than 2, coalescent in the spectral band spanning from 2.5 to 6 Hz. Two main peaks are present: one, more frequent at 2.5 Hz and the other at 5 Hz. This behaviour suggesting that the medium properties are likely subjected to slight, periodic variations potentially related to temporary fluctuation in water content that influences the propagation velocity.

At TOR4 (Figure 6) the HVSR exhibits two main peaks at frequencies of roughly 2.7 Hz and 5.5 Hz. The amplitude of these peaks varies according to a characteristic and systematic daily, and therefore weekly, behaviour, in which the largest amplitudes of noise are higher. This could be associated with (a) an artefact related to the internal electronic noise of the instrument, whose effects become relevant when the ground vibrations have very low amplitude, such as night time or during the weekend, or (b) the variations of the noise wavefield, as a consequence of the activation of different sources related to anthropogenic activities. Beginning in April, 2013, the amplitudes of these peaks start increasing from the values of 3, and by the end of the monitoring period they attain values around 5, that is, about 65% greater than those observed during the early phases of the experiment. Such amplitude increase is likely to reflect a corresponding increase of the impedance contrast between the unstable mass and the underlying solid bedrock. However, the peak frequency remains stable in time, indicating that both thickness and velocity of the shallowest layer remain substantially unchanged. Thus, an increment in the velocity and/or density of the underlying layer must be invoked in order to explain the inferred impedance variations. Potential phenomena provoking this possible velocity increase will be discussed in the following.

Also, the HVSR directivity throughout the 7 months of recording was analysed. As an example, contour maps in Figures 7(a) and 7(b) compare the medium directionality of data acquired in December 2012 and July 2013, respectively. For the two different intervals, the directivity at stations TOR1, TOR2, and TOR3 stayed substantially unchanged. On the other hand, the TOR4 polarization direction between the two periods change slightly (Figure 8), even though for the later interval directivity is clearer as a consequence of the amplitude increase of the horizontal components, as also manifested by the growing number of the HVSR peaks (Figure 6). This suggests that the observed temporal variations in the HVSR plots are not due to changes in the distribution of active sources; if this would be the case, consequently the polarization direction should most likely have changed.

5. Discussion

Assuming that the HVSR is strictly related to the dynamic properties of the medium and that it is supposed to be stable if no change occurred in the velocity and/or density of the ground [63], results from HVSR analyses can be summarized as follows: (i) there are clear configurations of quasiconstant or slowly varying contiguous frequencies whose H/V peak values depend on the considered station; (ii) the stations located on the sliding mass (TOR4) and at its head (TOR3), on potentially loose section, show an amplitude peak which is sharper and larger than those observed at the stations settled downstream. At TOR4, the amplitude variations of the HVSR cannot be unequivocally interpreted. However, the overall stationarity of the polarization properties suggests that those changes most likely reflect a variation in the acoustic properties of the medium rather than a change in the distribution of noise sources.

As mentioned in the Introduction section the quarry rockslide was extensively monitored since 2003 with traditional methods. Among these, as shown in Figure 2, there were 13 extensometers. All the extensometer data (E1-E15 in Figure 2) were individually normalized and compared with the measured cumulative rainfall in order to highlight a possible linear correlation between two different time series. The corr function in MATLAB was employed to both evaluate the linear (or rank) correlation (Rho) and perform a hypothesis test. The hypothesis was of no correlation against the alternative that there is a nonzero correlation (Pval) assuming by the authors that the correlation between two data is significant if Pval is sufficiently small (< 0.05). Table 1 shows the values obtained for each comparison.

The results of the correlation analysis clearly show that the deformational fields in the upper section of the quarry (E7, E8, E9, E10, E13, E14, located on the main cracks whose widths enlarge up to 2 m from East to West) and in the western part of the quarry (E2, located on the lateral crack) are strictly related to the seasonal rainfall, since the Pval values are very small (exponent lower than -100). This behaviour could be explained taking into account that, at sites where opening of the fractures is significant, pore water pressures in the fractures/cracks can critically influence the stability of rock. Unfortunately, because of a problem in the instrumentation, no data are available on the water level in the cracks. Moreover, from a qualitatively point of view, looking at Figure 9, it is possible to assess that periods characterized by the main soil movements (highlighted by the vertical sections of the extensometer curves) follow periods with higher rainfall (highlighted by the vertical sections of the cumulated rain curve). In particular, this behaviour is clear at the half of January, and at the end of February and May. The rainfall also seems to have a weaker but still significant influence on the deformations measured by E11, E3, E4, and E15 (Figure 9) while an inverse correlation exists with the data recorded by E1 (located in correspondence of the basal plane). Finally, there is no evidence of correlation for E12 data neither with the rain trend nor with the temperature variation. Unfortunately, there are no superficial evidences that could justify this behaviour, apart from the fact that E12 is located in correspondence of a junction between two main fractures (Figure 2). Perhaps its behaviour is caused by this junction (i.e., the highest movements are recorded by the extensometer around), or otherwise it could be possible that some errors in the data registration occurred.

A good match (Rho: 0.7147; Pval 9.3923 e−30) is obtained by comparing the TOR4 HVSR amplification value and temperature variation (Figure 10). This could be caused by the water content variation in the medium (water content alternatively empties and fills the rock pores) and consequently the relative variation, related to changes in the HVSR amplification value. The saturation of pores with water, in fact, tends to increase the velocity of P-waves (which propagate more efficiently through water than air), also increasing the Poisson ratio. This has a strong influence on the Rayleigh waves and, in particular, on the ellipticity of the particle motion with a consequent increase of the ratio between the horizontal component H and the vertical one V of ground motion [6466].

To justify the strong resemblance that emerges by comparing the temperature trend and that of the H/V recorded at TOR4, a direct dependency of this latter parameter on the meteorological conditions could be supposed. This behaviour is suggested in [67] that ascribes a fundamental role to barometric conditions variation concerning the composition of the noise wavefield. In [68] this behaviour is related to microseismic frequencies (lower than 1 Hz) and is related to oceanic storm waves. Reference [69] observes similar phenomena at very high latitudes: in that case the variations could be explained by cycles of freezing and thawing that crumble the rock surface and change their acoustic properties. Moreover, [70] points out how the HVSR amplitude could be affected by the local meteorological conditions (e.g., the wind). If this would be the case, an extended dataset (> 1 yr) would be necessary in order to clarify whether the observed variations at TOR4 are part of cyclical phenomenon occurring over longer periods as a consequence of seasonal changes. Unfortunately, at the Torgiovannetto quarry it was not possible to extend the experiment over longer time intervals because of the hard acquisition conditions. Nevertheless, the hypothesis that the HVSR amplitude value is directly related to meteorological factors can be excluded in Torgiovannetto area. All the stations, in fact, given their small spacing, should have shown the same amplitude increase. At the contrary, TOR1 shows a minimum in January 2013 and relatively constant values in the other months; TOR3 shows pronounced maxima on December 2012 and February 2013; TOR4 shows an increasing trend from January to July 2013. These H/V frequency variations could be associated with a different depth of fracturing (i.e., at site TOR1 fractures are shallower/near to the surface while at sites TOR3 and TOR4 they could be observed at depth) since the penetration of the surface waves is related to the frequency, but there are no experimental data on the depth of the fractures. Moreover, the surface geology at site TOR1 is characterized by stiffer or thicker geological unit, as indicated by the H/V frequency at 4.5 Hz. The surface geology at site TOR3 is characterized by softer/or thinner geological unit, as indicated by the H/V frequency at 2.5 Hz.

Probably the temperature variation does not directly affect the H/V amplitude but is responsible of other mechanisms like: (i) increasing of the fracturing degree of the medium acting directly on the dilatancy of the rocks (an increase in the medium fracturing may result, directly or, more often, indirectly, in density or velocity of propagation variations); (ii) influencing the water content of the superficial layer leading to an increase of the wave velocity of this portion of the slope. This hypothesis is supported by the variation in the impedance contrast which occurred with the approach of the hot season (early April, 2013) and therefore higher temperature (i.e., the water in the superficial layers is more prone to evaporate with the higher temperature).

Unfortunately, at the present state the lack of evidence of surface displacements corresponding to the observed variations in the HVSR amplitude trend foreclosed any possibility of threshold identification that could be used as an EWS. It could be interesting to evaluate this technique as a surveillance method when it can be calibrated on monitoring intervals characterized by a high rate of surface activity or over longer period in order to explain the cyclical variation of that parameter.

6. Conclusion

Implementing an EWS is a challenging issue in landslide monitoring. To verify the usefulness of seismic noise analysis as part of an EWS, a pilot scale experiment was arranged to monitor an unstable rock mass. A 7-month period of passive seismic data was analysed by means of the H/V method. Possible connection between rainfall/temperature/displacement and rockslide seismic activity was evaluated, and the hypothesis that the HV amplitude value is directly related to meteorological factors can be excluded. On the contrary, the H/V observed variations with time are interesting, in that they potentially reveal changes of subsoil site conditions and have also implications for the assessment of site response to seismic shaking. The presented analysis was just the first step to employ H/V variations in an EWS. Many efforts, in fact, have to be employed both to understand how the observed variations are correlated with slope stability conditions and to set up a reliable EWS. For the first point (a) a longer time acquisition period and (b) a comparison with many other parameters to model and interpret in a quantitative way are needed. There are many factors (like cracks, joints, rock diagenesis, and saturation), in fact, that could cause velocity or density variations and therefore influence the ellipticity and/or polarization of the surface waves. For the second point, there are some open questions like, (a) what is the main information that the EWS will receive, (b) how this information will be processed, (c) what are the preferred time responses, and (d) how the potential variations and/or errors from (a) and (b) will affect the false alarm/no alarm ratios of the EWS. Nonetheless, the rapid technological advances increasing the speed in acquisition, transmission, and processing of data suggest that it is clearly worthy to proceed in the field of seismic monitoring of unstable slopes.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The Department of Earth Sciences (Unifi-DST) supported this research as part of its program to improve rockslide early warning system (PRIN 2009 – Advanced monitoring techniques for the development of early warning procedures on large rockslides - prot. 20084FAHR7_001). We gratefully acknowledge Sara Electronic Instrument for providing the devices installed at the four seismic stations. We thank Massimiliano Nocentini and Luca Lombardi (Unifi-DST) for the huge efforts deployed to install, maintain, and make available the microseismic data and Francesco Ponziani (Centro Funzionale Regione Umbria) for providing meteorological data and very helpful comments and suggestions for the manuscript. Thanks are also due to Prof. Nicola Casagli (Unifi-DST) whose initial review was very helpful. The authors are also grateful to the anonymous reviewers for providing very helpful comments to improve the manuscript.