Temporal and spatial distributions of precursory seismicity rate changes in the Thailand-Laos-Myanmar border region: implication for upcoming hazardous earthquakes

To study the prospective areas of upcoming strong-to-major earthquakes, i.e., M _w  ≥ 6.0, a catalog of seismicity in the vicinity of the Thailand-Laos-Myanmar border region was generated and then investigated statistically. Based on the successful investigations of previous works, the seismicity rate change (Z value) technique was applied in this study. According to the completeness earthquake dataset, eight available case studies of strong-to-major earthquakes were investigated retrospectively. After iterative tests of the characteristic parameters concerning the number of earthquakes (N) and time window (T _w ), the values of 50 and 1.2 years, respectively, were found to reveal an anomalous high Z-value peak (seismic quiescence) prior to the occurrence of six out of the eight major earthquake events studied. In addition, the location of the Z-value anomalies conformed fairly well to the epicenters of those earthquakes. Based on the investigation of correlation coefficient and the stochastic test of the Z values, the parameters used here (N = 50 events and T _w  = 1.2 years) were suitable to determine the precursory Z value and not random phenomena. The Z values of this study and the frequency-magnitude distribution b values of a previous work both highlighted the same prospective areas that might generate an upcoming major earthquake: (i) some areas in the northern part of Laos and (ii) the eastern part of Myanmar.     

Coal Mining Induced Seismicity in the Ruhr Area, Germany

Over the last 25 years mining-induced seismicity in the Ruhr area has continuously been monitored by the Ruhr-University Bochum. About 1,000 seismic events with local magnitudes between 0.7 ≤ M _L ≤ 3.3 are located every year. For example, 1,336 events were located in 2006. General characteristics of induced seismicity in the entire Ruhr area are spatial and temporal correlation with mining activity and a nearly constant energy release per unit time. This suggests that induced stresses are released rapidly by many small events. The magnitude–frequency distribution follows a Gutenberg–Richter relation which is a result from combining distributions of single longwalls that themselves show large variability. A high b-value of about 2 was found indicating a lack of large magnitude events. Local analyses of single longwalls indicate that various factors such as local geology and mine layout lead to significant differences in seismicity. Stress redistribution acts very locally since differences on a small scale of some hundreds of meters are observed. A regional relation between seismic moment M ₀ and local magnitude M _L was derived. The magnitude–frequency distribution of a single longwall in Hamm was studied in detail and shows a maximum at M _L = 1.4 corresponding to an estimated characteristic source area of about 2,200 m². Sandstone layers in the hanging or foot wall of the active longwall might fail in these characteristic events. Source mechanisms can mostly be explained by shear failure of two different types above and below the longwall. Fault plane solutions of typical events are consistent with steeply dipping fracture planes parallel to the longwall face and nearly vertical dislocation in direction towards the goaf. We also derive an empirical relation for the decay of ground velocity with epicenter distance and compare maximum observed ground velocity to local magnitude. This is of considerable public interest because about 30 events larger than M _L ≥ 1.2 are felt each month by people living in the mining regions. Our relations, for example, indicate that an event in Hamm with a peak ground velocity of 6 mm/s which corresponds to a local magnitude M _L between 1.7 and 2.3 is likely to be felt within about 2.3 km radius from the event.     

Forecasting the Locations of Future Large Earthquakes: An Analysis and Verification

The objective of this paper is to quantify the use of past seismicity to forecast the locations of future large earthquakes and introduce optimization methods for the model parameters. To achieve this the binary forecast approach is used where the surface of the Earth is divided into l° × l° cells. The cumulative Benioff strain of m ≥ m _c earthquakes that occurred during the training period, ΔT _tr, is used to retrospectively forecast the locations of large target earthquakes with magnitudes ≥m _T during the forecast period, ΔT _for. The success of a forecast is measured in terms of hit rates (fraction of earthquakes forecast) and false alarm rates (fraction of alarms that do not forecast earthquakes). This binary forecast approach is quantified using a receiver operating characteristic diagram and an error diagram. An optimal forecast can be obtained by taking the maximum value of Pierce’s skill score.     

Did the M S7.0 Lushan earthquake dynamically trigger earthquakes in the Datong volcanic region (Shanxi Province)?

Immediately following the M _S7.0 Lushan earthquake on April 20, 2013, using high-pass and low-pass filtering on the digital seismic stations in the Shanxi Province, located about 870–1,452 km from the earthquake epicenter, we detected some earthquakes at a time corresponding to the first arrival of surface waves in high-pass filtering waveform. The earthquakes were especially noticed at stations in Youyu (YUY), Shanzizao (SZZ), Shanghuangzhuang (SHZ), and Zhenchuan (ZCH), which are located in a volcanic region in the Shanxi Province,but they were not listed in the Shanxi seismic observation report. These earthquakes occurred 4–50 min after the passage of the maximum amplitude Rayleigh wave, and the periods of the surface waves were mainly between 15 and 20 s following. The Coulomb stresses caused by the Rayleigh waves that acted on the four stations was about 0.001 MPa, which is a little lower than the threshold value of dynamic triggering, therefore, we may conclude that the Datong volcanic region is more sensitive to the Coulomb stress change. To verify, if the similar phenomena are widespread, we used the same filtering to observe contrastively continuous waveform data before, and 5 h after, the M _S7.0 Lushan earthquake and M _S9.0 Tohoku earthquake in 2011. The results show that the similar phenomena occur before the earthquakes, but the seismicity rates after the earthquakes are remarkably increased. Since these weak earthquakes are quite small, it is hard to get clear phase arrival time from three or more stations to locate them. In addition, the travel time differences between P waves and S waves (S–P) are all less than 4 s, that means the events should occur in 34 km around the stations in the volcanic region. The stress of initial dynamic triggering of the M _S9.0 Tohoku earthquake was about 0.09 MPa, which is much higher than the threshold value of dynamic triggering stress. The earthquakes after the M _S9.0 Tohoku earthquake are related to dynamic triggering stress, but the events before the earthquake cannot be linked to seismic events, but may be related to the background seismicity or from other kinds of local sources, such as anthropogenic sources (i.e., explosions). Using two teleseismic filtering, the small background earthquakes in the Datong volcanic region occur frequently, thus we postulate that previous catalog does not apply bandpass filter to pick out the weak earthquakes, and some of the observed weak events were not triggered by changes in the dynamic stress field.     

Seismic Hazard Evaluation in Western Turkey as Revealed by Stress Transfer and Time-dependent Probability Calculations

Western Turkey has a long history of destructive earthquakes that are responsible for the death of thousands of people and which caused devastating damage to the existing infrastructures, and cultural and historical monuments. The recent earthquakes of Izmit (Kocaeli) on 17 August, 1999 (M _w  = 7.4) and Düzce (M _w  = 7.2) on 12 November, 1999, which occurred in the neighboring fault segments along the North Anatolian Fault (NAF), were catastrophic ones for the Marmara region and surroundings in NW Turkey. Stress transfer between the two adjacent fault segments successfully explained the temporal proximity of these events. Similar evidence is also provided from recent studies dealing with successive strong events occurrence along the NAF and parts of the Aegean Sea; in that changes in the stress field due to the coseismic displacement of the stronger events influence the occurrence of the next events of comparable size by advancing their occurrence time and delimiting their occurrence place. In the present study the evolution of the stress field since the beginning of the twentieth century in the territory of the eastern Aegean Sea and western Turkey is examined, in an attempt to test whether the history of cumulative changes in stress can explain the spatial and temporal occurrence patterns of large earthquakes in this area. Coulomb stress changes are calculated assuming that earthquakes can be modeled as static dislocations in elastic half space, taking into account both the coseismic slip in large (M ≥ 6.5) earthquakes and the slow tectonic stress buildup along the major fault segments. The stress change calculations were performed for strike-slip and normal faults. In each stage of the evolutionary model the stress field is calculated according to the strike, dip, and rake angles of the next large event, whose triggering is inspected, and the possible sites for future strong earthquakes can be assessed. A new insight on the evaluation of future seismic hazards is given by translating the calculated stress changes into earthquake probability using an earthquake nucleation constitutive relation, which includes permanent and transient effects of the sudden stress changes.     

Seismicity anomalies before the great earthquake of <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 in the Kunlun Pass and its significance to earthquake prediction

A great earthquake of M _S=8.1 took place in the west of Kunlun Pass on November 14, 2001. The epicenter is located at 36.2°N and 90.9°E. The analysis shows that some main precursory seismic patterns appear before the great earthquake, e.g., seismic gap, seismic band, increased activity, seismicity quiet and swarm activity. The evolution of the seismic patterns before the earthquake of M _S=8.1 exhibits a course very similar to that found for earthquake cases with M _S≥7. The difference is that anomalous seismicity before the earthquake of M _S=8.1 involves in the larger area coverage and higher seismic magnitude. This provides an evidence for recognizing precursor and forecasting of very large earthquake. Finally, we review the rough prediction of the great earthquake and discuss some problems related to the prediction of great earthquakes.     

Long aftershock sequences in North China and Central US: implications for hazard assessment in mid-continents

Because seismic activity within mid-continents is usually much lower than that along plate boundary zones, even small earthquakes can cause widespread concerns, especially when these events occur in the source regions of previous large earthquakes. However, these small earthquakes may be just aftershocks that continue for decades or even longer. The recent seismicity in the Tangshan region in North China is likely aftershocks of the 1976 Great Tangshan earthquake. The current earthquake sequence in the New Madrid seismic zone in central United States, which includes a cluster of M ~ 7.0 events in 1811–1812 and a number of similar events in the past millennium, is believed to result from recent fault reactivation that releases pre-stored strain energy in the crust. If so, this earthquake sequence is similar to aftershocks in that the rates of energy release should decay with time and the sequence of earthquakes will eventually end. We use simple physical analysis and numerical simulations to show that the current sequence of large earthquakes in the New Madrid fault zone is likely ending or has ended. Recognizing that mid-continental earthquakes have long aftershock sequences and complex spatiotemporal occurrences are critical to improve hazard assessments.     

Stress Release and Source Scaling of the 2010–2011 Canterbury,New Zealand Earthquake Sequence from Spectral Inversion of Ground Motion Data

The Canterbury earthquake sequence beginning with the 2010 M _W 7.2 Darfield earthquake is one of the most notable and well-recorded crustal earthquake sequences in a low-strain-rate region worldwide and as such provides a unique opportunity to better understand earthquake source physics and ground motion generation in such a tectonic setting. Ground motions during this sequence ranged up to extreme values of 2.2 g, recorded during the February 2011 M _W 6.2 event beneath the city of Christchurch. A better understanding of the seismic source signature of this sequence, in particular the stress release and its scaling with earthquake size, is crucial for future ground motion prediction and hazard assessment in Canterbury, but also of high interest for other low-to-moderate seismicity regions where high-quality records of large earthquakes are lacking. Here we present a source parameter study of more than 200 events of the Canterbury sequence, covering the magnitude range M _W 3–7.2. Source spectra were derived using a generalized spectral inversion technique and found to be well characterized by the ω ^?2 source model. We find that stress drops range between 1 and 20 MPa with a median value of 5 MPa, which is a factor of 5 larger than the median stress drop previously estimated with the same method for crustal earthquakes in much more seismically active Japan. Stress drop scaling with earthquake size is nearly self-similar, and we identify lateral variations throughout Canterbury, in particular high stress drops at the fault edges of the two major events, the M _W 7.2 Darfield and M _W 6.2 Christchurch earthquakes.     

Stress Drops of the 1997–1998 Colfiorito,Central Italy Earthquakes: Hints for a Common Behaviour of Normal Faults in the Apennines

Stress drop estimates of moderate-magnitude earthquakes in the Umbria–Marche region, in the northern Apennines, exhibit a large scatter. For the two M _w 5.7 and 6.0 main shocks of 26 September 1997 near Colfiorito, several papers resulted in stress drop estimates of 20 MPa, but values as low as 2–3 MPa were proposed as well. Also for the largest aftershocks (M _w > 4), estimates spread from < 1 MPa up to values ten times larger. We have critically revisited methods and data used in the literature. We have specifically faced the trade-off between source and propagation effects, as we believe that it is responsible for a part of the large scatter. To keep this trade-off under control, we have applied a methodology that combines the best fit of both source spectra after Empirical Green’s Function (EGF) deconvolution and observed ground motion spectra, finding that the results of the two different data sets converge independently at the same solution. We have used ground motions observed in the Colfiorito basin, where an accelerograph and a co-located seismological broad-band station recorded three clusters of earthquakes in a broad magnitude interval (1.7 ≤ M _w ≤ 6.0). We have found that the mainshock–aftershock sequences result in stress drops of 2–5 MPa at M _w ≥ 5.6, with an average tendency to decrease at smaller magnitudes where stress drop variability increases. These findings confirm the source scaling recently assessed through Empirical Green’s Function deconvolution for another well-monitored seismic sequence of normal-faulting earthquakes, which struck the city of L’Aquila in the central Apennines in April 2009. The similar scaling law of the two areas suggests common mechanisms of stress release for the shallow normal faults in the Apennines. The propensity of smaller earthquakes to increase in variability, with a tendency toward smaller stress drops, may reflect an effect of fault strength heterogeneities for smaller size ruptures.     

Semi-Periodic Sequences and Extraneous Events in Earthquake Forecasting: I. Theory and Method,Parkfield Application

We present a new method to identify semi-periodic sequences in the occurrence times of large earthquakes, which allows for the presence of multiple semi-periodic sequences and/or events not belonging to any identifiable sequence in the time series. The method, based on the analytic Fourier transform, yields estimates of the departure from periodicity of an observed sequence, and of the probability that the sequence is not due to chance. These estimates are used to make and to evaluate forecasts of future events belonging to each sequence. Numerous tests with synthetic catalogs show that the method is surprisingly capable of correctly identifying sequences, unidentifiable by eye, in complicated time series. Correct identification of a given sequence depends on the number of events it contains, on the sequence’s departure from periodicity, and, in some cases, on the choice of starting and ending times of the analyzed time window; as well as on the total number of events in the time series. Some particular data combinations may result in spectra where significant periods are obscured by large amplitudes artifacts of the transform, but artifacts can be usually recognized because they lack harmonics; thus, in most of these cases, true semi-periodic sequences may not be identified, but no false identifications will be made. A first example of an application of the method to real seismicity data is the analysis of the Parkfield event series. The analysis correctly aftcasts the September 2004 earthquake. Further applications to real data from Japan and Venezuela are shown in a companion paper.     

Large Earthquake Hazard of the San Jacinto Fault Zone,CA, from Long Record of Simulated Seismicity Assimilating the Available Instrumental and Paleoseismic Data

We investigate spatio-temporal properties of earthquake patterns in the San Jacinto fault zone (SJFZ), California, between Cajon Pass and the Superstition Hill Fault, using a long record of simulated seismicity constrained by available seismological and geological data. The model provides an effective realization of a large segmented strike-slip fault zone in a 3D elastic half-space, with heterogeneous distribution of static friction chosen to represent several clear step-overs at the surface. The simulated synthetic catalog reproduces well the basic statistical features of the instrumental seismicity recorded at the SJFZ area since 1981. The model also produces events larger than those included in the short instrumental record, consistent with paleo-earthquakes documented at sites along the SJFZ for the last 1,400 years. The general agreement between the synthetic and observed data allows us to address with the long-simulated seismicity questions related to large earthquakes and expected seismic hazard. The interaction between m ≥ 7 events on different sections of the SJFZ is found to be close to random. The hazard associated with m ≥ 7 events on the SJFZ increases significantly if the long record of simulated seismicity is taken into account. The model simulations indicate that the recent increased number of observed intermediate SJFZ earthquakes is a robust statistical feature heralding the occurrence of m ≥ 7 earthquakes. The hypocenters of the m ≥ 5 events in the simulation results move progressively towards the hypocenter of the upcoming m ≥ 7 earthquake.     

Earthquakes Along the Passive Margin of Greenland: Evidence for Postglacial Rebound Control

—?An intriguing observation in Greenland is a clear spatial correlation between seismicity and deglaciated areas along passive continental margins, a piece of evidence for earthquake triggering due to postglacial rebound. Another piece of evidence for induced seismicity due to deglaciation derives from earthquake source mechanisms. Sparse, low magnitude seismicity has made it difficult to determine focal mechanisms from Greenland earthquakes. On the basis of two normal faulting events along deglaciated margins and from the spatial distribution of epicenters, earlier investigators suggested that the earthquakes of Greenland are due to postglacial rebound. This interpretation is tested here by using more recent data. Broadband waveforms of teleseismic P waves from the August 10, 1993 (m _b = 5.4) and October 14, 1998 (m _b = 5.1) earthquakes have been inverted for moment tensors and source parameters. Both mechanisms indicate normal faulting with small strike-slip components: the 1993 event, strike = 348.9°, dip = 41.0°, rake =?56.3°, focal depth = 11?km, seismic moment = 1.03?×?10²⁴ dyne-cm, and M _w = 5.3; the 1998 event, strike = 61.6°, dip = 58.0°, rake =?95.5°, focal depth = 5?km, seismic moment = 5.72?×?10²³ dyne-cm, and M _w = 5.1. These and the two prior events support the theory that the shallow part of the lithosphere beneath the deglaciated margins is under horizontal extension. The observed stress field can be explained as flexural stresses due to removal of ice loads and surface loads by glacial erosion. These local extensional stresses are further enhanced by the spreading stress of continental crust and reactivate preexisting faults. Earthquake characteristics observed from Greenland suggest that the dominant seismogenic stresses are from postglacial rebound and spreading of the continental lithosphere.     

Attenuation of High-Frequency Seismic Waves in Eastern Iran

We investigated the frequency-dependent attenuation of the crust in Eastern Iran by analysis data from 132 local earthquakes having focal depths in the range of 5–25 km. We estimated the quality factor of coda waves (Q _c) and body waves (Q _p and Q _s) in the frequency band of 1.5–24 Hz by applying the single backscattering theory of S-coda envelopes and the extended coda-normalization method, respectively. Considering records from recent earthquakes (Rigan M _w 6.5, 2010/12/20, Goharan M _w 6.2, 2013/5/11 and Sirch M _w 5.5, 2013/1/21), the estimated values of Q _c, Q _p and Q _s vary from 151 ± 49, 63 ± 6, and 93 ± 14 at 1.5 Hz to 1,994 ± 124, 945 ± 84 and 1,520 ± 123 at 24 Hz, respectively. The average frequency-dependent relationships (Q = Q _o f ⁿ ) estimated for the region are Q _c = (108 ± 10)f ^(0.96±0.01), Q _p = (50 ± 5)f ^(1.01±0.04), and Q _s = (75 ± 6)f ^(1.03±0.06). These results evidenced a frequency dependence of the quality factors Q _c, Q _p, and Q _s, as commonly observed in tectonically active zones characterized by a high degree of heterogeneity, and the low value of Q indicated an attenuative crust beneath the entire region.     

Multi-Methods Combined Analysis of Future Earthquake Potential

Prior to an earthquake, natural seismicity is correlated across multiple spatial and temporal scales. Many studies have indicated that an earthquake is hard to accurately predict by a single time-dependent precursory method. In this study, we attempt to combine four earthquake prediction methods, i.e. the Pattern Informatics (PI), Load/Unload Response Ratio (LURR), State Vector (SV), and Accelerating Moment Release (AMR) to estimate future earthquake potential. The PI technique is founded on the premise that the change in the seismicity rate is a proxy for the change in the underlying stress. We first use the PI method to quantify localized changes surrounding the epicenters of large earthquakes to objectively quantify the anomalous areas (hot spots) of the upcoming events. Next, we delineate the seismic hazard regions by integrating with regional active fault zones and small earthquake activities. Then, we further evaluate the earthquake potential in the seismic hazard regions using the LURR, SV and AMR methods. Retrospective tests of this new approach on the large earthquakes (M > 6.5) which have occurred in western China over the last 3 years show that the LURR and SV time series usually climb to an anomalously high peak months to years prior to occurrence of a large earthquake. And, the asymptote time, t _c, “predicted” by the AMR method correspond to the time of the actual events. The results may suggest that the multi-methods combined approach can be a useful tool to provide stronger constraints on forecasts of the time and location of future large events.     

Seismicity study in Pournari reservoir area (W. Greece) 1981–2010

The spatial-temporal evolution of seismicity is examined, during the initial impoundment of Pournari reservoir located on Arachthos River (Western Greece), as well as for the next 30 years. The results show that, despite the relatively moderate-to-high seismicity from west to east, there is no remarkable earthquake in the vicinity before the first reservoir impoundment. Immediately after the impoundment (January 1981), and during the first 4 months, a considerable number of low-magnitude seismic events were recorded in the broader area of the dam. Moreover, two independent major events occurred on March 10, 1981 (M _L ?=?5.6) and April 10, 1981 (M _L ?=?4.7) with focal depths 13 and 10 km, respectively. The detailed analysis of the two corresponding aftershock sequences shows that they present different behaviors (e.g., larger b-value and lower magnitude of the main aftershock) than that of other aftershock sequences in Greece. This seismicity is probably due to triggering, via the water loading mechanism and the undrained response due to a flysch appearance on the reservoir basement. The activation of the thrust fault may be attributed to the bulging of evaporites that characterize the disordered structure of W. Greece, via possible water intake. The detailed processing of the recorded seismicity during the period 1982–2010, in comparison with the variations of Pournari Dam water level, shows an increase of shallow seismicity (h?≤?5 km) in the vicinity of the reservoir up to a 10-km distance—in contrast to the initial period, characterized by a number of deeper events due to the background response change from undrained to drained status.     

Location Performance and Detection Threshold of the Spanish National Seismic Network

Spain is a low-to-moderate seismicity area with relatively low seismic hazard. However, several strong shallow earthquakes have shaken the country causing casualties and extensive damage. Regional seismicity is monitored and surveyed by means of the Spanish National Seismic Network, maintenance and control of which are entrusted to the Instituto Geográfico Nacional. This array currently comprises 120 seismic stations distributed throughout Spanish territory (mainland and islands). Basically, we are interested in checking the noise conditions, reliability, and seismic detection capability of the Spanish network by analyzing the background noise level affecting the array stations, errors in hypocentral location, and detection threshold, which provides knowledge about network performance. It also enables testing of the suitability of the velocity model used in the routine process of earthquake location. To perform this study we use a method that relies on P and S wave travel times, which are computed by simulation of seismic rays from virtual seismic sources placed at the nodes of a regular grid covering the study area. Given the characteristics of the seismicity of Spain, we drew maps for M _L magnitudes 2.0, 2.5, and 3.0, at a focal depth of 10 km and a confidence level 95 %. The results relate to the number of stations involved in the hypocentral location process, how these stations are distributed spatially, and the uncertainties of focal data (errors in origin time, longitude, latitude, and depth). To assess the extent to which principal seismogenic areas are well monitored by the network, we estimated the average error in the location of a seismic source from the semiaxes of the ellipsoid of confidence by calculating the radius of the equivalent sphere. Finally, the detection threshold was determined as the magnitude of the smallest seismic event detected at least by four stations. The northwest of the peninsula, the Pyrenees, especially the westernmost segment, the Betic Cordillera, and Tenerife Island are the best-monitored zones. Origin time and focal depth are data that are far from being constrained by regional events. The two Iberian areas with moderate seismicity and the highest seismic hazard, the Pyrenees and Betic Cordillera, and the northwestern quadrant of the peninsula, are the areas wherein the focus of an earthquake is determined with an approximate error of 3 km. For M _L 2.5 and M _L 3.0 this error is common for almost the whole peninsula and the Canary Islands. In general, errors in epicenter latitude and longitude are small for near-surface earthquakes, increasing gradually as the depth increases, but remaining close to 5 km even at a depth of 60 km. The hypocentral depth seems to be well constrained to a depth of 40 km beneath the zones with the highest density of stations, with an error of less than 5 km. The M _L magnitude detection threshold of the network is approximately 2.0 for most of Spain and still less, almost 1.0, for the western sector of the Pyrenean region and the Canary Islands.     

An improved earthquake catalogue (<Emphasis Type="Italic">M</Emphasis> ≥ 4.0) for Turkey and its near vicinity (1900–2012)

Many catalogues, agency reports and research articles have been published on seismicity of Turkey and its surrounding since 1950s. Given existing magnitude heterogeneity, erroneous information on epicentral location, event date and time, this past published data however is far from fulfilling the required standards. Paucity of a standardized format in the available catalogues have reinforced the need for a refined and updated catalogue for earthquake related hazard and risk studies. During this study, ~37,000 earthquakes and related parametric data were evaluated by utilizing more than 41 published studies and, an integrated database was prepared in order to analyse all parameters acquired from the catalogues and references for each event. Within the scope of this study, the epicentral locations of M ≥ 5.0 events were firstly reappraised based on the updated Active Fault Map of Turkey. An improved catalogue of 12.674 events for the period 1900–2012 was as a result recompiled for the region between 32–45N° and 23–48E° by analyzing in detail accuracy of all seismological parameters available for each event. The events consist of M ≥ 4.0 are reported in several magnitude scales (e.g. moment magnitude, M_w; surface wave magnitude, M_S; body-wave magnitude m_b; local magnitude M_L and duration magnitude M_d) whereas the maximum focal depth reaches up to 225-km. In order to provide homogenous data, the improved catalogue is unified in terms of M_w. Fore-and aftershocks were also removed from the catalogue and completeness analyses were performed both separately for various tectonic sources and as a whole for the study region of interest. Thus, the prepared homogenous and declustered catalogue consisting of 6573 events provides the basis for a reliable input to the seismic hazard assessment studies for Turkey and its surrounding areas.     

An Evaluation of Coulomb Stress Changes from Earthquake Productivity Variations in the Western Gulf of Corinth,Greece

Spatial and temporal evolution of the stress field in the seismically active and well-monitored area of the western Gulf of Corinth, Greece, is investigated. The highly accurate and vast regional catalogues were used for inverting seismicity rate changes into stress variation using a rate/state-dependent friction model. After explicitly determining the physical quantities incorporated in the model (characteristic relaxation time, fault constitutive parameters, and reference seismicity rates), we looked for stress changes across space and over time and their possible association with earthquake clustering and fault interactions. We focused our attention on the Efpalio doublet of January 2010 (M = 5.5 and M = 5.4), with a high aftershock productivity, and attempted to reproduce and interpret stress changes prior to and after the initiation of this seismicity burst. The spatial distribution of stress changes was evaluated after smoothing the seismological data by means of a probability density function (PDF). The inverted stress calculations were compared with the calculations derived from an independent approach (elastic dislocation model) and this comparison was quantified. The results of the two methods are in good agreement (up to 80 %) in the far field, with the inversion technique providing more robust results in the near field, where they are more sensitive to the uncertainties of coseismic slip distribution. It is worth mentioning that the stress inversion model proved to be a very sensitive stress meter, able to detect even small stress changes correlated with spatio–temporal earthquake clustering. Data analysis was attempted from 1975 onwards to simulate the stress changes associated with stronger earthquakes over a longer time span. This approach revealed that only M > 5.5 events induce considerable stress variations, although in some cases there was no evidence for such stress changes even after an M > 5.5 earthquake.     

Gorkha earthquake (MW7.8) and aftershock sequence: A fractal approach

On April 25, 2015, Nepal was struck by the M_W7.8 Gorkha earthquake followed by an intense aftershock sequence. It was one of the most destructive earthquakes in the Himalayan arc, causing more than 8900 fatalities. In this study, we analyzed the dataset (429 events, magnitude of completeness (M_c) ≥ 4.2 local magnitude) of the first 45 days after the Gorkha earthquake to estimate the seismicity parameters b-value, D-value, and p-value. We used the maximum likelihood method to estimate the b-value and Omori-Utsu parameters, whereas the correlation integral method was applied to estimate the fractal dimension (D-value). The analysis was carried out using running and sliding window techniques. The lowest b-value (0.57 ± 0.04) and the highest D-value (1.65 ± 0.02) were computed at the time of the Gorkha earthquake, after which the b-value significantly increased to a maximum of 1.57. It again dropped to 0.93 at the time of the major aftershock on May 12, 2015. The D-value showed an initial quick drop and then decreased in a wavy pattern until the end of the study period, indicating the clustering and scattering of earthquakes in a fault region. The b-value contour map identified the eastern part of the study area as a high stress region (b = ～0.8), implying that the stress shifted to that region. The D-value contour map reveals that the seismogenic structure shifted from linear to planar in the region. The rate of aftershock decay (p = 0.86 ± 0.04) for a short period reflects that the level of stress decreased rapidly. This study helps to understand the level of stress and seismicity pattern of a region, which could be useful for aftershock studies.