Site response analysis of Vartholomio W-Greece from singular spectrum analysis of microtremor and weak motion data

Twenty six sites were instrumented in the city of Vartholomio following the December 2, 2002 Ms 6.0 earthquake. Thirty weak events from the aftershock sequence as well as microtremors were used to identify amplifications due to geological site effects. Horizontal-to-vertical spectral ratios (HVSR—Nakamura estimates) and weak events ratios were calculated and the singular spectrum analysis (SSA) method was used. The results showed that the effects of SSA on the stability of the frequency peak and amplitude distribution of HVSR for both weak motion and microtremors. The data analysis confirms the role of near surface geology in causing locally significant variations of the predominant frequencies and amplitudes of ground shaking as already inferred from the distribution of damages. The site response spectra exhibited significant peaks within the range of 1.5–2.6 Hz and the amplification factor did not exceed 6.5. Finally the parts of the HVSR ratios from ～0.2 up to 10 Hz were used, in order to create an automatic optimal zonation of the study area using a genetic algorithm. This procedure resulted in the division of the city into 2 main zones.     

Geo-center movement caused by huge earthquakes

In this study we investigate co-seismic geo-center change based on a dislocation theory for a spherically symmetric, non-rotating, perfectly elastic and isotropic model. We first introduce the basic theory with emphasis on the dislocation Love numbers of degree 1, and then we present methods for computing co-seismic geo-center movement. It is found that the geo-center change reaches maximum value when δ = 45° and λ = 90°, i.e., a 45° declined dip fault causes the maximum geo-center movement. As an application, we apply the methods to compute the geo-center movement caused by the 2004 Sumatra earthquake (Mw9.3) and the 2011 Tohoku-Oki earthquake (Mw9.0). Results show that the maximum co-seismic geo-center movements for the two events are 0.87 mm and 0.43 mm, respectively.     

Seismic response of deep Quaternary sediments in historical center of L’Aquila City (central Italy)

On 6 April 2009 a M_w=6.1 earthquake produced severe destruction and damage over the historic center of L’Aquila City (central Italy), in which the accelerometer stations AQK and AQU recorded a large amount of near-fault ground motion data. This paper analyzes the recorded ground motions and compares the observed peak accelerations and the horizontal to vertical response spectral ratios with those revealed from numerical simulations. The finite element method is considered herein to perform dynamic modeling on the soil profile underlying the seismic station AQU. The subsurface model, which is based on the reviewed surveys that were carried out in previous studies, consists of 200–400 m of Quaternary sediments overlying a Meso-Cenozoic carbonate bedrock. The Martin-Finn-Seed's pore-water pressure model is used in the simulations. The horizontal to vertical response spectral ratio that is observed during the weak seismic events shows three predominant frequencies at about 14 Hz, 3 Hz and 0.6 Hz, which may be related to the computed seismic motion amplification occurring at the shallow colluvium, at the top and base of the fluvial-lacustrine sequence, respectively. During the 2009 L’Aquila main shock the predominant frequency of 14 Hz shifts to lower values probably due to a peculiar wave-field incidence angle. The predominant frequency of 3 Hz shifts to lower values when the earthquake magnitude increases, which may be associated to the progressive softening of soil due to the excess pore-water pressure generation that reaches a maximum value of about 350 kPa in the top of fluvial-lacustrine sequence. The computed vertical peak acceleration underestimates the experimental value and the horizontal to vertical peak acceleration ratio that is observed at station AQU decreases when the earthquake magnitude increases, which reveals amplification of the vertical component of ground motion probably due to near-source effects.     

Multiresolution wavelets and natural time analysis before the January–February 2014 Cephalonia (Mw6.1 & 6.0) sequence of strong earthquake events

On January 26 and February 3, 2014, Cephalonia Island (Ionian Sea, Greece) was struck by two strong, shallow earthquakes (moment magnitudes Mw6.1 and Mw6.0, respectively) that ruptured two sub-parallel, strike-slip faults, with right-lateral kinematics. The scope of the present work is to investigate the complex correlations of the earthquake activity that preceded the Mw6.1 event in the broader area of the Cephalonia Island and identify possible indications of critical stages in the evolution of the earthquake generation process. We apply the recently introduced methods of Multiresolution Wavelet Analysis (MRWA) and Natural Time (NT) analysis and for the first time we combine their results in a joint approach that may lead to universal principles in describing the evolution of the earthquake activity as it approaches a major event. In particular, the initial application of MRWA on the inter-event time series indicates a time marker 12 days prior to the major event. By using this time as the initiation point of the NT analysis, the critical stage of seismicity, where the κ₁ parameter reaches the critical value of κ₁ = 0.070, is approached few days before the occurrence of the Mw6.1 earthquake.     

Refining the OSL age of the last earthquake on the Dheshir fault,Central Iran

In Central Iran there are several cities along the Dehshir fault, which have similar geological conditions to that of the city of Bam prior to the 2003 earthquake (Mw 6.5), during which more than 30,000 lives were lost. Optical stimulated luminescence (OSL) samples were collected from the Dehshir fault in order to place constraints on its seismic history. The single aliquot regenerative (SAR) dose measurement protocol on coarse grained quartz extracts was used for this study. This SAR protocol had to be optimized for the low OSL sensitivity by varying both the preheat temperatures and test doses used. Dose recovery tests showed that given laboratory dose could be successfully recovered. However, replicate palaeodose (D_e) data were scattered and consequently ages based on mean D_e's had large uncertainties. As this is thought to largely reflect poor bleaching conditions prior to sediment burial at the site, various statistical procedures were employed in conjunction with the stratigraphic knowledge of the site to try and extract more refined burial ages from the samples. From this the timing of the last earthquake was estimated around 2.0 ± 0.2 kyr. This refined age suggests that the earthquake catalogue of Iran is incomplete and more paleoseismological investigation is required to recognize and date the previous events of Dheshir fault.     

TEC variations over the Mediterranean before and during the strong earthquake (M = 6.5) of 12th October 2013 in Crete,Greece

In this paper, the total electron content (TEC) data from eight global positioning system (GPS) stations of the EUREF network, provided by IONOLAB (Turkey), were analyzed using discrete Fourier analysis to investigate the TEC variations over the Mediterranean before and during the strong earthquake of 12th October 2013, which occurred west of Crete, Greece. In accordance with the results of similar analyses in the area, the main conclusions of this study are the following: (a) TEC oscillations in a broad range of frequencies occur randomly over an area of several hundred km from the earthquake and (b) high frequency oscillations (f ⩾ 0.0003 Hz, periods T ⩽ 60 m) may point to the location of the earthquake with questionable accuracy. The fractal characteristics of the frequency distribution may point to the locus of the earthquake with higher accuracy. We conclude that the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism through acoustic or gravity waves could explain this phenomenology.     

Relocation of aftershocks of the 2001 Bhuj earthquake: A new insight into seismotectonics of the Kachchh seismic zone,Gujarat, India

In view of an anomalous crust–mantle structure beneath the 2001 Bhuj earthquake region, double-difference relocations of 1402 aftershocks of the 2001 Bhuj earthquake were determined, using an improved 1D velocity model constructed from 3D velocity tomograms based on data from 10 to 58 three-component seismograph stations. This clearly delineated four major tectonic features: (i) south-dipping north Wagad fault (NWF), (ii and iii) south-dipping south Wagad faults 1 and 2 (SWF₁, SWF₂), and (iv) a northeast dipping transverse fault (ITF), which is a new find. The relocated aftershocks correlate satisfactorily with the geologically mapped and inferred faults in the epicentral region. The relocated focal depths delineate a marked variation to the tune of 12 km in the brittle–ductile transition depths beneath the central aftershock zone that could be attributed to a lateral variation in crustal composition (more or less mafic) or in the level of fracturing across the fault zone. A fault intersection between the NWF and ITF has been clearly mapped in the 10–20 km depth range beneath the central aftershock zone. It is inferred that large intraplate stresses associated with the fault intersection, deepening of the brittle–ductile transition to a depth of 34 km due to the presence of mafic/ultramafic material in the crust–mantle transition zone, and the presence of aqueous fluids (released during the metamorphic process of eclogitisation of lower crustal olivine-rich rocks) and volatile CO₂ at the hypocentral depths, might have resulted in generating the 2001 Bhuj earthquake sequence covering the entire lower crust.     

Active segment of the 12 November 2003 Mw 5.6 earthquake at Salsipuedes oceanic basin,Gulf of California. Mexico

We analyzed records of eight seismic stations of the autonomous broadband seismograph network of a joint project between Utrecht University (the Netherlands), California Institute of Technology, and Centro de Investigación Científica y de Estudios Superiores de Ensenada (CICESE). These stations recorded the Mw 5.6 earthquake that occurred on 12 November 2003 at Salsipuedes basin in the middle of the Gulf of California 2 km west of the island Angel de la Guarda. This event was located at 29.16º N and 113.37º W, 30 km northeast of Bahia de los Angeles. A foreshock and hundreds of aftershocks were recorded in the 48 hours after its origin time. With the location of 29 earthquakes we identified the active segment, perpendicular to the main transform fault NW–SE of Canal de Ballenas, representing the transtensional boundary between the Pacific and North American plates. The direction of the active fault described is consistent with the normal fault mechanism reported by the National Earthquake Information Center (strike=39º, dip=34º, slip=–44º).From the duration magnitude of 456 aftershocks, we calculated a b-value of 1.14±0.28; furthermore, we calculated a seismic moment of (3.5 ±3.3) X10¹⁷Nm, a source radius of 3.7 ± 2.63 km, and a static stress drop of 3.94 ± 1.15 MPa (39.4 ± 11.5 bar.).     

Estimation of Shallow S-Wave Velocity Structure in Damascus City,Syria, Using Microtremor Exploration

Array measurements of microtremors were carried out at thirty sites in Damascus city, Syria to estimate S-wave velocity structures of shallow soil formations for site effect analysis. The microtremor data were recorded by 6 vertical-component seismometers distributed along the circumferences of two circles as well as a 3-component seismometer deployed in the center. The phase velocities were estimated at each site from the vertical components of recorded microtremor data by using the Spatial Autocorrelation method. Then, Genetic Simulated Annealing Algorithm technique was applied for inversion of the phase velocities to estimate 1-D S-wave velocity structures beneath the sites. The inverted Vs profiles are not uniform in Damascus city and the results show that a shallow soft layer (∼200 m/s) appears in the eastern part of the city as well as the central part along Barada River. This layer controls the amplification distribution in the city with a high amplification mainly observed at the locations having this layer. The inversion results also show that the depth to the engineering bedrock (∼750 m/s) is very shallow along the foothills of Mt. Qasyoun in the north-west. Then the depth increases towards the east and the south. The maximum depth to the engineering bedrock (∼80 m) was observed in the southern part of Damascus. To validate the results of the inversions, the spectral ratios between the horizontal and vertical components (H/V) of the recorded microtremor data at the central seismometer were compared with the computed ellipticities of the fundamental-mode Rayleigh-waves based on the respective Vs structure. The results show a good agreement in a period range of 0.05 s to 0.5 s. In this period range, the dominant peaks of the H/V ratios are due to the overall effect of the velocity contrasts between the shallow layers representing the subsurface S-wave velocity structure. Moreover, the average S-wave velocity for the top 10 m of soils (V_S10) shows a better correlation with the averaged site amplification in a period range of 0.05 s to 0.5 s than V_S30 which indicates that V_S10 can be a better proxy for high-frequency site amplification in the case of Damascus city.     

Analysis of observations backing up the existence of VLF and ionospheric TEC anomalies before the Mw6.1 earthquake in Greece,January 26, 2014

The present work integrates ground-based ionosphere measurements using very-low-frequency radio transmissions with satellite measurements of the total electron content to draw common conclusions about the possible impact that the Mw6.1 earthquake that took place in Greece on January 26, 2014, had on the ionosphere.Very-low-frequency radio signals reveal the existence of an ∼4-day anomaly in the wavelet spectra of the signals received inside the earthquake preparation zone and a significant increase in the normalized variance of the signals prior to the earthquake (approximately 1 day before).Through total electron content analysis, it was possible to identify a clear anomaly from 15:00 until 20:00 UT on the day before the earthquake that appears again on the day of the earthquake between 07:00 UT and 08:00 UT. The anomalous values reach TEC1Sigma ∼4.36 and 3.11, respectively. Their spatial and temporal distributions give grounds to assume a possible link with the earthquake preparation. The geomagnetic, solar and weather conditions during the considered period are presented and taken into account.This work is an initial and original step towards a multi-parameter approach to the problem of the possible earthquake-related effects on the ionosphere joining observations made from both ground stations and satellites. A well-founded knowledge of these phenomena is clearly necessary before dealing with their application to earthquake prediction purposes.     

Post-seismic relaxation: An example of earthquake triggering in the Apennine belt (1915–1920)

The strong earthquake (M = 7) that occurred in the Fucino basin (central Italy) on January 13, 1915 was followed by six earthquakes of M > 5.5 and several other shocks of M > 5 in the major seismic zones of the northern Apennines from 1916 to 1920. This seismicity pattern is consistent with the implications of the present tectonic setting in the study area, which suggests that strong decoupling earthquakes in the central Apennines cause a significant increase of tectonic load, and possibly of seismicity, in the northern Apennines. A numerical simulation, carried out by an elastic-viscous model, of the stress diffusion induced by the Fucino and successive largest earthquakes, shows that each of the above shocks occurred when the respective zone was reached by the highest values of the strain and strain rate perturbation triggered by the previous events. Furthermore, the computed strain regime at each earthquake site is consistent with the known faulting pattern. The results provide important insights into the physical mechanism that controls the interaction of seismic sources in the central and northern Apennines.     

Inferred fault geometry and slip distribution of the 2010 Jiashian,Taiwan, earthquake is consistent with a thick-skinned deformation model

We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 10¹⁸ N-m is equivalent to a M_w = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.     

Spatio-temporal variability of size fractionated phytoplankton on the shelf adjacent to the Ebro river (NW Mediterranean)

The mesoscale distribution and seasonal variation of the size structure of phytoplankton biomass, as measured by chlorophyll a (chl a), was studied in the Ebro shelf area (NW Mediterranean) during three different seasons: autumn, winter and summer. In autumn and summer, when the water column was, respectively, slightly or strongly stratified and nutrient concentrations were low at surface, average total chl a values were 0.31 and 0.29 mg m⁻³, respectively. In winter, the intrusion of nutrients into the photic zone by intense vertical mixing and strong riverine inputs, produced an increase of the total autotrophic biomass (0.76 mg m⁻³). In the three seasons, the main contributor to total chl a was the picoplanktonic (<2 μm) size fraction (42% in winter and around 60% in autumn and summer). The nanophytoplankton (2–20 μm) contribution to total chl a showed the lowest variability amongst seasons (between 29% and 39%). The microplanktonic (>20 μm) chl a size fraction was higher in winter (27%) than in the other seasons (less than 13%). The maximum total chl a concentrations were found at surface in winter, at depths of 40 m in autumn and between 50 and 80 m in summer. The relative contribution of the <2 μm size fraction at these levels of the water column tended to be higher than at other depths in autumn and winter and lower in summer. In autumn and winter, nutrient inputs from Ebro river discharge and mixing processes resulted in an increase on the >2 μm contribution to total chl a in the coastal zone near the Ebro Delta area. In summer, the contribution of the <2 and >2 μm chl a size fractions was homogeneously distributed through the sampling area. In autumn and summer, when deep chl a maxima were observed, the total amount of the autotrophic biomass in the superficial waters (down to 10 m) of most offshore stations was less than 10% of the whole integrated chl a (down to 100 m or to the bottom). In winter, this percentage increased until 20% or 40%. The >2 μm chl a increased linearly with total chl a values. However, the <2 μm chl a showed a similar linear relationship only at total chl a values lower than 1 mg m⁻³ (in autumn and summer) or 2 mg m⁻³ (winter). At higher values of total chl a, the contribution of the <2 μm size fraction remained below an upper limit of roughly 0.5 mg m⁻³. Our results indicate that the picoplankton fraction of phytoplankton may show higher seasonal and mesoscale variability than is usually acknowledged.     

Seismic microzonation studies for the city of Ragusa (Italy)

The seismic history of the city of Ragusa (Italy), the geotechnical characterisation of the subsoil and the site response analysis should be correctly evaluated for the definition of the Seismic Geotechnical Hazard of the city of Ragusa, through geo-settled seismic microzoning maps. Basing on the seismic history of the city of Ragusa, the following earthquake scenarios have been considered: the “Val di Noto” earthquake of January 11, 1693 (with intensity X–XI on MCS scale, magnitude M_W=7.41 and epicentral distance of about 53 km); the “Etna” earthquake of February 20, 1818 (with intensity IX on MCS scale, magnitude M_W=6.23 and epicentral distance of about 64 km); the Vizzini earthquake of April 13, 1895 (with intensity I=VII–VIII on MCS scale, magnitude M_W=5.86 and epicentral distance of about 26 km); the “Modica” earthquake of January 23, 1980 (with intensity I=V–VI on MCS scale, magnitude M_W=4.58 and epicentral distance of about 10 km); the “Sicilian” earthquake of December 13, 1990 (with intensity I=VII on MCS scale, magnitude M_W=5.64 and epicentral distance of about 50 km). Geotechnical characterisation has been performed by in situ and laboratory tests, with the definition of shear wave velocity profiles in the upper 30 m of soil. Soil response analyses have been evaluated for about 120 borings location by some non-linear 1-D models. Finally the seismic microzonation of the city of Ragusa has been obtained in terms of maps with different peak ground acceleration at the surface; shaking maps for the central area of the city of Ragusa were generated via GIS for the earthquake scenarios.     

Bayesian analysis on earthquake magnitude related to an active fault in Taiwan

It is understood that sample size could be an issue in earthquake statistical studies, causing the best estimate being too deterministic or less representative derived from limited statistics from observation. Like many Bayesian analyses and estimates, this study shows another novel application of the Bayesian approach to earthquake engineering, using prior data to help compensate the limited observation for the target problem to estimate the magnitude of the recurring Meishan earthquake in central Taiwan. With the Bayesian algorithms developed, the Bayesian analysis suggests that the next major event induced by the Meishan fault in central Taiwan should be in M_w 6.44±0.33, based on one magnitude observation of M_w 6.4 from the last event, along with the prior data including fault length of 14 km, rupture width of 15 km, rupture area of 216 km², average displacement of 0.7 m, slip rate of 6 mm/yr, and five earthquake empirical models.     

Patterns of significant seismic quiescence on the Mexican Pacific coast

Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, M_w = 7.8), Petatlán (March 1979, M_w = 7.6), Michoacán (September 1985, M_w = 8.0, and M_w = 7.6) and Colima (October 1995, M_w = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

Dating Holocene lavas on Stromboli,Italy using cosmogenic He

In-situ cosmogenic ³He exposure ages of pyroxene phenocrysts from basalts from the Upper Neostromboli formation in southwest Stromboli date its eruption at 7.0 ± 0.3 ka (1σ, n = 3, Ginostra site) and 6.8 ± 0.2 ka (1σ, n = 10, Timpone del Fuoco site) respectively. Correlation of our new data to previous K/Ar and palaeomagnetic ages from the northwestern Neostromboli phase suggests that it erupted within a confined period between roughly 6 and 14 ka. The low uncertainty on the ³He_cos ages as well as on individual exposure ages (4.4–8.7%) demonstrates that ³He_cos exposure dating is a viable tool for dating Holocene basalt lavas. The ages compare favourably to uncertainties obtained for radiocarbon dating of similar rocks.     

Temporal and spatial variations of seismicity scaling behavior in Southern México

R/S analysis is used in this work to investigate the fractal correlations in terms of the Hurst exponent for the 1998–2011 seismicity data in Southern Mexico. This region is the most seismically active area in Mexico, where epicenters for severe earthquakes (e.g., September 19, 1985, M_w = 8.1) causing extensive damage in highly populated areas have been located. By only considering the seismic events that meet the Gutenberg–Ritcher law completeness requirement (b = 0.97, M_GR = 3.6), we found time clustering for scales of about 100 and 135 events. In both cases, a cyclic behavior with dominant spectral components at about one cycle per year is revealed. It is argued that such a one-year cycle could be related to tidal effects in the Pacific coast. Interestingly, it is also found that high-magnitude events (M_w ≥ 6.0) are more likely to occur under increased interevent correlations with Hurst exponent values H > 0.65. This suggests that major earthquakes can occur when the tectonic stress accumulates in preferential directions. In contrast, the high-magnitude seismic risk is reduced when stresses are uniformly distributed in the tectonic shell. Such cointegration between correlations (i.e., Hurst exponent) and macroseismicity is confirmed for spatial variations of the Hurst exponent. In this way, we found that, using the Hurst exponent standpoint, the former presumed Michoacan and the Guerrero seismic gaps are the riskiest seismic zones. To test this empirical finding, two Southern Mexico local regions with large earthquakes were considered. These are the Atoyac de Alvarez, Guerrero (M_w = 6.3), and Union Hidalgo, Oaxaca (M_w = 6.6), events. In addition, we used the Loma Prieta, California, earthquake (October 17, 1989, M_w = 6.9) to show that the high-magnitude earthquakes in the San Andreas Fault region can also be linked to the increments of determinism (quantified in terms of the Hurst exponent) displayed by the stochastic dynamics of the interevent period time series. The results revealed that the analysis of seismic activity by means of R/S analysis could provide further insights in the advent of major earthquakes.