The Scientific World Journal

Volume 2014 (2014), Article ID 531212, 8 pages

http://dx.doi.org/10.1155/2014/531212

Research Article

Seismicity Pattern Changes before the M = 4.8 Aeolian Archipelago (Italy) Earthquake of August 16, 2010

Salvatore Gambino,¹ Antonino Laudani,² and Salvatore Mangiagli¹

¹Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Catania Osservatorio Etneo, Piazza Roma 2, 95123 Catania, Italy

²Dipartimento di Ingegneria, Università di Roma Tre, Via V. Volterra 62, 00146 Roma, Italy

Received 12 November 2013; Accepted 28 November 2013; Published 8 January 2014

Academic Editors: H. Steffen and Y.-M. Wu

Copyright © 2014 Salvatore Gambino et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigated the seismicity patterns associated with an earthquake recorded in the Aeolian Archipelago on 16, August, 2010, by means of the region-time-length (RTL) algorithm. This earthquake triggered landslides at Lipari; a rock fall on the flanks of the Vulcano, Lipari, and Salina islands, and some damages to the village of Lipari. The RTL algorithm is widely used for investigating precursory seismicity changes before large and moderate earthquakes. We examined both the spatial and temporal characteristics of seismicity changes in the Aeolian Archipelago region before the earthquake. The results obtained reveal 6-7 months of seismic quiescence which started about 15 months before the earthquake. The spatial distribution shows an extensive area characterized by seismic quiescence that suggests a relationship between quiescence and the Aeolian Archipelago regional tectonics.

1. Introduction

The quiescence of seismic activity has been defined as the notable decrease in the seismic activity against the average background. Temporal seismic observations have shown trends of seismic quiescence preceding large and moderate events [1, 2]. Successively Sobolev and Tyupkin [3, 4] proposed the region-time-length algorithm (RTL algorithm), a statistical method for the investigation of the seismic activity level preceding large earthquakes.

This method may evidence a decrease (quiescence) or an increase (activation) in the seismic activity against the average background [5]. The RTL method has previously been applied to earthquakes in Kamchatka and Caucasus [4, 6, 7], Japan (e.g., [7–9]), China [10, 11], Greece [12], Turkey [13], Taiwan [14], and India [15]. Some moderate Italian earthquakes have been studied by Di Giovanbattista and Tyupkin [16–18], Gentili and Bressan [19], and Gentili [20] by using this technique.

The Aeolian Archipelago (Figure 1) is located in the Southern Tyrrhenian Sea (Italy) and represents the manifestation of a submarine volcanic arc originating in the central sectors of the Tyrrhenian Sea during the Pliocene and successively migrating towards the southeast.

Figure 1: Map of the investigated area with main structural features and seismic network. The largest earthquakes, which occurred in Aeolian Archipelago in the last 50 years, are reported.

It can be subdivided into three sectors with a different structural and tectonic evolution [21]. In the western sector, comprising the Alicudi and Filicudi islands, the volcanic activity started at about 1.3 Ma [22] and ended at about 30–40 kyr. At present, the seismicity occurs in the crust along the WNW-ESE Sisifo fault system (Figure 1). The eastern sector, which comprises Panarea and Stromboli islands and where volcanism developed from 0.8 Ma ago and is still active and is affected by a prevailing NE-SW striking fault system. The central sector includes the islands of Salina, Lipari, and Vulcano. Here, the volcanism began at 0.4 Myr [23] and is still active (last eruption 1888–1890) at Lipari and Vulcano (e.g., [24–27]).

These volcanoes are aligned along a lithospheric NNW-SSE fault system, the Aeolian-Tindari-Letojanni fault system (Figure 1) with right-lateral to oblique kinematics along which the seismicity is roughly aligned (e.g., [28–30]). Earthquakes occur mostly in the upper 20 km of the crust [31]; in particular, the seismicity west of Tindari-Letoianni fault system is distributed in a 7–18 km interval of depth, whereas earthquakes of Lipari-Vulcano eastern area are not deeper than 12-13 km [32].

The southern Tyrrhenian area is characterized by seismicity with maximum magnitude in the range of 5-6; in the last 50 years two strong events have been recorded: (15/04/1978) and (28/05/1980) (Figure 1, [33]). These two moderate earthquakes have marked an increase of the regional dynamics that, according to Chiodini et al., 1992 [34], and Montalto, 1996 [35], caused the reactivation of the volcanic system on Vulcano. Moreover the occurrence of an earthquake of regional significance shortly before the last eruption of Vulcano [36] confirms that a moderate seismic event could initiate a rapid magma ascent.

The seismicity recorded from 1999 to 2011 comprises events with and the 16, August, 2010 (), one represents the event with the highest magnitude recorded. In this study, the region-time-length (RTL) algorithm has been implemented to the catalogue of earthquakes which occurred in the period from 2000 to 2010 and we discuss the phases of seismic activation and quiescence preceding the event in 2010.

2. Data

Since the late ‘70s, continuous seismic monitoring activity in the Aeolian Archipelago has been performed by a permanent seismic network made up of a few analogical 3C stations. Starting from the ‘80s, the network was augmented with other stations deployed over the entire Aeolian Archipelago and equipped with short-period seismometers, having a natural frequency of 1 Hz. During 2005 and 2007, almost all the stations were replaced by new digital 24-bit ones, equipped with broadband (40 s) three-component sensors, with a dynamic range of 144 dB. To date, the Aeolian permanent seismic network, managed by INGV-CT (Istituto Nazionale di Geofisica e Vulcanologia-Sezione di Catania), consists of 12 three-component digital seismic stations (Figure 1). In order to reduce the azimuthal gap, events location is obtained also using the stations deployed in the Calabro-Peloritan area and on the northern flanks of Mt. Etna (Figure 1). Furthermore, where possible, we added data from the INGV national permanent seismic network.

We considered an area of  km with a latitude between 38.00 and 38.85 and a longitude of 14.00 and 15.30. The dataset used in this study comprises 1680 crustal earthquakes recorded from August 1999 to 2011 with magnitude , whose location, performed for surveillance purposes, is obtained by using the Hypoellipse code [37] (Figure 2). The mean errors of the analytical locations are, respectively, 0.95 km for the epicentral coordinates and 1.15 km for focal depth; the mean root mean square (RMS) is 0.16 s.

Figure 2: Seismic activity recorded in the Aeolian Archipelago area during the August 1999–December 2011 period. The dataset comprises 1680 earthquakes recorded from August 1999 to 2011; the red circle shows location of the 16, August, 2010, earthquake.

The main event of the catalogue occurred at 12.54 GMT of the 16, August, 2010, when Aeolian Archipelago was shaken by an earthquake of an estimated 4.8 magnitude. The hypocenter of the earthquake was situated 8 km west-south-west of the island of Vulcano at a depth of 13.0 km b.s.l. (Figure 1). The earthquake was felt on the northern coast of Sicily and in the cities of Palermo, Catania, and Messina. The earthquake triggered some landslides at Lipari, a rock fall on the flanks of Vulcano Lipari and Salina. In the village of Lipari minor damages to buildings and roads were reported and some beaches were closed for safety reasons.

3. Method

The analysis of the earthquake dataset has been performed by using the well-established method known as RTL algorithm [4, 8] which uses three parameters, namely, (region around the earthquake epicenter), (time), and (rupture length). The fundamental idea of RTL algorithm is to assign a weighting RTL value to a given spatiotemporal value (), which comes from events occurring in a prescribed space-time window within the characteristic distance and time. An RTL parameter is defined as the product of , , and describing the influence weights of location, occurrence time, and magnitude as

where is the rupture dimension (a function of magnitude ); is the occurrence time of the th earthquake; is the distance from the position ) to the epicenter of the th event; and are the characteristic distance and time associated with the spatiotemporal criteria; is the cut-off depth; and is the number of events satisfying the following criterion:

( is the magnitude of th event and is the cut-off magnitude ensuring the completeness of the earthquake catalogue); ; and .

For the rupture dimension the following expression is used [16]:

, , and are the background values of , , and , respectively, obtained as the expected values in the time interval considered for the analyzed position. RTL parameter describes the deviation from the background level of seismicity and is expressed in units of the standard deviation. A negative RTL value indicates a lower seismicity and a positive RTL value indicates a higher seismicity compared to the background. Clearly, both a temporal and a spatial analysis of RTL can be performed and some authors often use the spatial average value for the RTL parameters calling it parameter [11]. The algorithm for the computation of , , , and parameters has been implemented in the MATLAB environment (Figure 3), allowing simple management and plotting of the results.

Figure 3: Flowchart that documents the algorithm used on Matlab for RTL analyses. The input variable “data” contains all the information about the events to be analyzed. test and , are respectively, the position and the time to be studied, whereas IT1 and IT2 individuate the interval of time to be considered for the RTL analysis. Clearly the constants “ref_,” “ref_,” “ref,” and “ref” correspond, respectively, to , , , and .

4. RTL’s Calculation and Results

The RTL analysis needs to be applied to declustered catalogues, where aftershocks are removed [9]. In order to decluster the INGV catalogue, we applied the Reasemberg [38] algorithm implemented in Zmap software [39].

The Reasenberg algorithm defines a seismic sequence as a chain of events linked to each other by spatial and temporal windows. The variables are , the factor for the interaction radius of dependent event, , the look-ahead time for un-clustered events in days, , the maximum look-ahead time for clustered events in days, and , a measure of the confidence that the next event in the sequence is being observed.

For declustering we used the default parameters (, , , and ) obtaining a catalogue of 1324 events.

Moreover the magnitude of completeness has been evaluated for the catalogue by using the Gutenberg-Richter relation of earthquake frequency and magnitude.

The completeness of the data entries depends on the characteristics of the seismic network. The geometry, sensitivity, and resolution of the seismic network quantify different parts of the region in order to judge the behavior of the seismic regimes based on the representation of minimum magnitude.

We applied the Gutenberg-Richter relation to search the , representative magnitude for the present earthquake time series of the Aeolian Archipelago region.

The power law of Gutenberg-Richter fits the earthquake energy distribution as a linear plot of recurrence. The bending of the linear plot for the smaller magnitude earthquake gives an indication of incompleteness of the catalogue below a specified magnitude. This specified magnitude is the minimum or the threshold magnitude for the studied area. Figure 4 shows the earthquake frequency magnitude plot. It may be noted from here that the data are complete for earthquakes of . After removing events with the remaining data comprised 838 events which were used for the present study to estimate the RTL variation.

Figure 4: -value estimation of the complete earthquake dataset (August 1999 to 2011) after removing the aftershocks in the Aeolian Archipelago region.

We calculated the RTL and parameters [9, 13], that is, the possible time and spatial variation of the seismic quiescence, and to this end we made some choices about the input parameters; if we consider, for example, the earthquake, which occurred in the western region of Tottori prefecture, Japan, on 6, October, 2000, Huang [9] adopted a distance  km, year, and a focal depth () of 30 km.

Shashidhar et al. [15], for moderate earthquakes () in a small area (20 km × 30 km), tried different values for parameters adopting the following values:  km, days, and  km.

In order to obtain the RTL, we have considered location of the 16, August, earthquake () and the 1, January, 2008–15, August, 2010 (958 days), period.

We set a focal depth of 30 km considering that almost all earthquakes (98.7%) are not more than 30 km deep (Figure 5) and we tried different values of and (Figure 6). RTL algorithm does not show large differences between the different curves (Figure 6); we adopted the following model parameters: a characteristic distance = 25 km and = 50 days. Finally, we also ran (Figure 7) the RTL algorithm for the entire 2000–August 2010 period.

Figure 5: Cross-section of the selected seismicity catalogue to show depth distribution of all seismicity (grey diamonds) and after removing aftershocks and (purple diamonds).

Figure 6: Temporal variation of the RTL at the epicenter of the 16, August, 2010,  earthquake. The variation in RTL anomaly is obtained at this location with different values of parameters reported on the left of the graphs. The three values are referred to as , , and . Starting point is 1, January, 2008 and duration is 958 days. Some examples of seismic activation are indicated by arrows.

Figure 7: 2000–August 2010 time variation of the RTL at the epicenter of the 16, August, 2010,  earthquake by using , , and as input parameters. Red diamonds identify the occurrence times of the earthquakes located nearby (10–15 km) the 16, August, 2010, epicenter.

5. Discussion

In many parts of the world, the RTL method has been used for larger regions and longer seismic catalogues, obtaining valid observations on the quiescence phenomenon prior to large earthquakes [11, 40]. In this study, the RTL algorithm has been implemented in the MATLAB environment and tested to a triggered earthquake time series occurring in the Aeolian Archipelago (Italy), a relatively small area (100 km × 80 km) characterized by a moderate seismicity. A phase of seismic quiescence (between the 700th and the 850th day in Figure 6) was detected by the parameter around the epicenter of the 16, August, 2010, earthquake. The seismic quiescence spans the June–December 2009 period, ending 8-9 months before the earthquake.

Figures 6 and 7 show the presence of short-time positive changes (activations) during the 2000–2010 period. These short changes are linked to the occurrence of earthquakes located nearby (10–15 km) the epicenter. Moreover some modest and rapid negative changes are also visible.

All these short-time variations may be related to an imperfect removing of the aftershocks in the catalogue and/or to the adopted parameters.

Six-seven-month quiescence ending 8-9 months before the earthquake is in accordance with duration of the quiescence (0.6 to 3 years) and time shift from the end of the quiescence to the earthquake (0 to 2.9 years) was found by Gentili [20] for several Italian sectors.

Finally, in order to investigate its possible spatial variation, we calculated the parameter in the quiescence period (June–December 2009) for the entire area (Figure 8). To this end, the territory has been divided into 2,500 cells, each one with an area of ca. 4 km².

Figure 8: Spatial variation of RTL in the Aeolian Archipelago area during the observed quiescence period (June–December 2009). The scale on the right corresponds to the RTL value in the units of the standard deviation. The white circle shows location of the 16 August 2010, earthquake.

The main area covered by the quiescence (ca. 200 km²) comprises a sector around Vulcano island. The earthquake occurred in this area, which is about 4 km west of the pixel with lower (−3.3) value. However, an area between Salina and Filicudi of ca. 80 km² and a small sector (ca. 40 km²) in the north of Sicily also show negative values. Areas covered by the quiescence agree with sectors affected by Sisifo and Tindari-Letojanni fault system.

Considering the location (Figures 1 and 8) we would expect an quiescence focussed only on Tindari-Letojanni fault system. As different tectonic zones have different background seismicity [13] our results suggest a continuity between the two structures and a relationship between the quiescence recorded before the and overall Aeolian Archipelago regional tectonics.

6. Conclusions

The results obtained reveal a seismic quiescence phase before an earthquake with an extensive (ca. 320 km²) spatial distribution which comprises the triggered zone. However the area covered by the quiescence seems large in order to affirm that the method presented here is an effective tool of improving the significance and reliability of earthquake precursors.

These features are consistent with the results obtained by different authors by using RTL worldwide and encourage us to improve RTL analyses on other earthquakes by testing different parameters in order to evaluate the future possibility of moderate earthquake occurrence in this region that could also have an impact on the volcanic system of Vulcano.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors thank the “INGV-CT Gruppo Analisi Dati Sismici” for their help in elaborating the earthquakes.

References

M. Wyss and R. E. Habermann, “Precursory seismic quiescence,” Pure and Applied Geophysics, vol. 126, no. 2–4, pp. 319–332, 1988.

S. Wiemer and M. Wyss, “Seismic quiescence before the Landers (M = 7.5) and Big Bear (M = 6.5) 1992 earthquakes,” Bulletin of the Seismological Society of America, vol. 84, no. 3, pp. 900–916, 1994.

G. A. Sobolev and Y. S. Tyupkin, “New method of intermediate-term earthquake prediction,” in Seismology in Europe, pp. 229–234, ESC, Reykjavik, Iceland, 1996.

G. A. Sobolev and Y. S. Tyupkin, “Low-seismicity precursors of large earthquakes in Kamchatka,” Volcanology and Seismology, vol. 18, no. 4, pp. 433–446, 1997.

Y. M. Wu, C. C. Chen, L. Zhao, and C. H. Chang, “Seismicity characteristics before the 2003 Chengkung, Taiwan, earthquake,” Tectonophysics, vol. 457, no. 3-4, pp. 177–182, 2008.

G. A. Sobolev and Y. S. Tyupkin, “Precursory phases, seismicity precursors, and earthquake prediction in Kamchatka,” Volcanology and Seismology, vol. 20, no. 6, pp. 615–627, 1999.

Q. Huang, “Seismicity pattern changes prior to large earthquakes—an approach of the RTL algorithm,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 15, no. 3, pp. 469–491, 2004.

Q. Huang, G. A. Sobolev, and T. Nagao, “Characteristics of the seismic quiescence and activation patterns before the M = 7.2 Kobe earthquake, January 17, 1995,” Tectonophysics, vol. 337, no. 1-2, pp. 99–116, 2001.

Q. Huang, “Search for reliable precursors: a case study of the seismic quiescence of the 2000 Western Tottori prefecture earthquake,” Journal of Geophysical Research B, vol. 111, no. 4, 2006.

D. Rong and Y. Li, “Estimation of characteristic parameters in region-time-length algorithm and its application,” Acta Seismologica Sinica, vol. 20, no. 3, pp. 265–272, 2007.

Q. Huang, “Seismicity changes prior to the Ms 8.0 Wenchuan earthquake in Sichuan, China,” Geophysical Research Letters, vol. 35, no. 23, Article ID L23308, 2008.

G. A. Sobolev, Y. S. Tyupkin, and A. Zavialov, “Map of expectation earthquakes algorithm and RTL prognostic parameter: joint application,” in Proceedings of the 29th General Assembly of IASPEI, Abstracts 77, Thessaloniki, Greece, 1997.

Q. Huang, A. O. Öncel, and G. A. Sobolev, “Precursory seismicity changes associated with the Mw = 7.4 1999 August 17 Izmit (Turkey) earthquake,” Geophysical Journal International, vol. 151, no. 1, pp. 235–242, 2002.

C. C. Chen and Y. X. Wu, “An improved region-time-length algorithm applied to the 1999 Chi-Chi, Taiwan earthquake,” Geophysical Journal International, vol. 166, no. 3, pp. 1144–1148, 2006.

D. Shashidhar, N. Kumar, K. Mallika, and H. Gupta, “Characteristics of seismicity patterns prior to the M~5 earthquakes in the Koyna Region, Western India—application of the RTL algorithm,” Episodes, vol. 33, no. 2, pp. 83–90, 2010.

R. Di Giovambattista and Y. Tyupkin, “The fine structure of the dynamics of seismicity before M ≥ 4.5 earthquakes in the area of Reggio Emilia (Northern Italy),” Annali di Geofisica, vol. 42, no. 5, pp. 897–909, 1999.

R. Di Giovambattista and Y. S. Tyupkin, “Saptial and temporal distribution of seismicity before the Umbria-Marche September 26,1997 eartquakes,” Journal of Seismology, vol. 4, no. 4, pp. 589–598, 2000.

R. Di Giovambattista and Y. S. Tyupkin, “Seismicity patterns before the M = 5.8 2002, Palermo (Italy) earthquake: seismic quiescence and accelerating seismicity,” Tectonophysics, vol. 384, no. 1–4, pp. 243–255, 2004.

S. Gentili and G. Bressan, “Seismicity patterns before MD C 4.1 earthquakes in the Friuli-Venezia Giulia (Northeastern Italy) and Western Slovenia areas,” Bollettino di Geofisica Teorica e Applicata, vol. 48, pp. 33–51, 2007.

S. Gentili, “Distribution of Seismicity Before the Larger Earthquakes in Italy in the Time Interval 1994–2004,” Pure and Applied Geophysics, vol. 167, no. 8-9, pp. 933–958, 2010.

G. de Astis, G. Ventura, and G. Vilardo, “Geodynamic significance of the Aeolian volcanism (Southern Tyrrhenian Sea, Italy) in light of structural, seismological and geochemical data,” Tectonics, vol. 22, no. 4, 2003.

P. Y. Gillot, “Histoire volcanique des Iles Eoliennes: arc insulaire or complexe orogenique anulaire?” Documents et Travaux Institute Géologique Albert-de-Lapparent, vol. 11, pp. 35–42, 1987.

L. Beccaluva, G. Gabbianelli, F. Lucchini, P. L. Rossi, and C. Savelli, “Petrology and K Ar ages of volcanics dredged from the Eolian seamounts: implications for geodynamic evolution of the Southern Tyrrhenian basin,” Earth and Planetary Science Letters, vol. 74, no. 2-3, pp. 187–208, 1985.

A. Peccerillo, M. L. Frezzotti, G. De Astis, and G. Ventura, “Modeling the magma plumbing system of Vulcano (Aeolian Islands, Italy) by integrated fluid-inclusion geobarometry, petrology, and geophysics,” Geology, vol. 34, no. 1, pp. 17–20, 2006.

S. Gambino, O. Campisi, G. Falzone et al., “Tilt measurements at Vulcano Island,” Annals of Geophysics, vol. 50, no. 2, pp. 233–247, 2007.

S. Gambino and F. Guglielmino, “Ground deformation induced by geothermal processes: a model for La Fossa Crater (Vulcano Island, Italy),” Journal of Geophysical Research B, vol. 113, no. 7, 2008.

S. Alparone, A. Cannata, S. Gambino, S. Gresta, V. Milluzzo, and P. Montalto, “Time-space variation of volcano-seismic events at La Fossa (Vulcano, Aeolian Islands, Italy): new insights into seismic sources in a hydrothermal system,” Bulletin of Volcanology, vol. 72, no. 7, pp. 803–816, 2010.

G. Neri, G. Barberi, B. Orecchio, and A. Mostaccio, “Seismic strain and seismogenic stress regimes in the crust of the Southern Tyrrhenian region,” Earth and Planetary Science Letters, vol. 213, no. 1-2, pp. 97–112, 2003.

G. Neri, G. Barberi, G. Oliva, and B. Orecchio, “Spatial variations of seismogenic stress orientations in Sicily, South Italy,” Physics of the Earth and Planetary Interiors, vol. 148, no. 2–4, pp. 175–191, 2005.

A. Billi, G. Barberi, C. Faccenna, G. Neri, F. Pepe, and A. Sulli, “Tectonics and seismicity of the Tindari fault system, Southern Italy: crustal deformations at the transition between ongoing contractional and extensional domains located above the edge of a subducting slab,” Tectonics, vol. 25, no. 2, 2006.

A. Cannata, I. S. Diliberto, S. Alparone et al., “Multiparametric approach in investigating hydrothermal systems: the case of study of Vulcano (Aeolian Islands, Italy),” Pure and Applied Geophysics, vol. 169, no. 1-2, pp. 167–182, 2012.

S. Gambino, V. Milluzzo, A. Scaltrito, and L. Scarfì, “Relocation and focal mechanisms of earthquakes in the south-central sector of the Aeolian Archipelago: new structural and volcanological insights,” Tectonophysics, vol. 524-525, pp. 108–115, 2012.

G. Neri, D. Caccamo, O. Cocina, and A. Montalto, “Geodynamic implications of earthquake data in the Southern Tyrrhenian sea,” Tectonophysics, vol. 258, no. 1–4, pp. 233–249, 1996.

G. Chiodini, R. Cioni, S. Falsaperla, A. Montalto, M. Guidi, and L. Marini, “Geochemical and seismological investigations at Vulcano (Aeolian Islands) during 1978–1989,” Journal of Geophysical Research, vol. 97, no. 7, pp. 11–32, 1992.

A. Montalto, “Signs of potential renewal of eruptive activity at La Fossa (Vulcano, Aeolian Islands),” Bulletin of Volcanology, vol. 57, no. 7, pp. 483–492, 1996.

G. Mercalli and O. Silvestri, “L'eruzione dell’Isola di Vulcano incominciata il 3 agosto 1888 e terminata il 22 marzo 1890,” Annali dell’ufficio Centrale di Meteorologia e Geodinamica, vol. 10, pp. 71–281, 1891.

Gruppo Analisi Dati Sismici, “Catalogo dei terremoti della Sicilia Orientale—Calabria Meridionale (1999–2012),” INGV, CT, 2013, http://www.ct.ingv.it/ufs/analisti/catalogolist.php

P. A. Reasemberg, “Second-order moment of Central California seismicity, 1969–1982,” Journal of Geophysical Research, vol. 90, no. 7, pp. 5479–5495, 1985.

S. Wiemer, “A software package to analyze seismicity: ZMAP,” Seismological Research Letters, vol. 72, no. 3, pp. 373–382, 2001.

G. A. Sobolev, “Microseismic variations prior to a strong earthquake,” Izvestiya, Physics of the Solid Earth, vol. 40, no. 6, pp. 455–464, 2004.