Static stress changes as a triggering mechanism of a shallow earthquake: case study of the 1999 Chi-Chi (Taiwan) earthquake

The role of static stress changes in triggering an earthquake has long been debated in the fields of geophysics and fault mechanics. Valuable data sets for the study of static triggering were provided within the 1-year period following the devastating 1999 Chi-Chi (Taiwan) earthquake (M_W=7.6), during which more than 20,000 aftershocks occurred. In this study, stress waves generated by the Chi-Chi earthquake were calculated using a source rupture model in conjunction with a layered elastic model. Static (permanent) stress changes were extracted from the long-period offsets in the stressgrams. Correlations between the calculated stress changes and seismicity were analyzed at different depths and over varying time intervals to ascertain the impact effects of stress changes on triggering aftershocks. Correlations between prior seismicity rates and static stress changes imposed by the Chi-Chi event were low, while correlations between late seismicity rates and static stress changes were much higher. This indicates that static stress changes did affect the occurrence of the Chi-Chi aftershock sequence. The percentage of early aftershocks at shallow depths (0-10 km) in static stress-enhanced areas within 2 weeks of the main shock was high but decreased considerably at greater depths (>10 km) and over longer time periods. It is concluded that static stress changes at depths of 0-10 km played a major role in triggering crustal aftershocks, especially those that occurred within 2 weeks of the main shock. In the deeper crust, static stress changes may have been modified by viscous flow, and at later times, perturbed by earlier, larger aftershocks. Although the correlations between seismicity rate changes and static stress changes are imperfect, a region that was anti-triggered is detected when these two results are compared. Static stress changes are presumably not the only aftershock triggering mechanism, but they definitively play a major role in triggering shallow aftershocks.     

Anomalous Seismicity and Accelerating Moment Release Preceding the 2001 and 2002 Earthquakes in Northern Baja California, Mexico

— An algorithm recently developed by RUNDLE et al. (2002) to find regions of anomalous seismic activity associated with large earthquakes identified the location of an M _w = 5.6 earthquake near Calexico, Mexico. In this paper we analyze the regional seismicity before this event, and a nearby M _w = 5.7 event, using time-to-failure algorithms developed by BOWMAN et al. (1998) and BOWMAN and KING (2001a,b). The former finds the radius of a circular region surrounding the epicenter that optimizes the time-to-failure acceleration of seismic release. The latter optimizes acceleration based on the expected stress accumulation pattern for a dislocation source. Both methods found a period of accelerating seismicity in an optimal region, the size of which agrees with previously proposed scaling relations. This positive result suggests that the Rundle algorithm may provide a useful technique to identify regions of accelerating seismicity, which can then be analyzed using signal optimization time-to-failure techniques.     

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.     

Seismotectonic of Southern Apennines from recent passive seismic experiments

We used data of local earthquakes collected during two recent passive seismic experiments carried out in southern Italy in order to study the seismotectonic setting of the Lucanian Apennine and the surrounding areas. Based on continuous recordings of the temporary stations we extracted over 15,600 waveforms, which were hand-picked along with those recorded by the permanent stations of the Italian national seismic network obtaining a dense, high-quality dataset of P- and S-arrival times. We examined the seismicity occurring in the period 2001–2008 by relocating 566 out of 1047 recorded events with magnitudes M_L ≥ 1.5 and computing 162 fault-plane solutions. Earthquakes were relocated using a minimum one-dimensional velocity model previously obtained for the region and a V_p/V_s ratio of 1.83. Background seismicity is concentrated within the upper crust (between 5 and 20 km of depth) and it is mostly clustered along the Lucanian Apennine chain axis. A significant feature extracted from this study relates to the two E–W trending clusters located in the Potentino and in the Abriola–Pietrapertosa sector (central Lucania region). Hypocentral depths in both clusters are slightly deeper than those observed beneath the Lucanian Apennine. We suggest that these two seismic features are representative of the transition from the inner portion of the chain to the external margin characterized by dextral strike-slip kinematics. In the easternmost part of the study area, below the Bradano foredeep and the Apulia foreland, seismicity is generally deeper and more scattered. The sparse seismicity localized in the Sibari Plain, in the offshore area along the northeastern Calabrian coast and in the Taranto Gulf is also investigated thanks to the new recordings. This seismicity shows hypocenters between 12 and 20 km of depth below the Sibari Plain and is deeper (foci between 10 and 35 km of depth) in the offshore area of the Taranto Gulf. 102 well-constrained fault-plane solutions, showing predominantly normal and strike-slip character with tensional axes (T-axes) generally NE oriented, were selected for the stress tensor analysis. We investigated stress field orientation inverting focal mechanism belonging to the Lucanian Apennine and the Pollino Range, both areas characterized by a more concentrated background seismicity.     

Seismicity and stress field in Okinawa Trough and Ryukyu regions

（黄培华）（苏维加）（陈金波）ＳｅｉｓｍｉｃｉｔｙａｎｄｓｔｒｅｓｓｆｉｅｌｄｉｎＯｋｉｎａｗａＴｒｏｕｇｈａｎｄＲｙｕｋｙｕｒｅｇｉｏｎｓ￥Ｐｅｉ－ＨｕａＨＵＡＮＧ；Ｗｅｉ－ｊｉａＳＵａｎｄＪｉｎ－ＢｏＣＨＥＮ（ＤｅｐａｒｔｍｅｎｔｏｆＥａｒｔｈａ...     

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.     

Microseismicity Induced by Heavy Rainfall Around Flooded Vertical Ore Veins

—?Microseismicity (M?M?M?M??6 were also monitored at a crustal movement monitoring station located several hundred meters from the veins. It was found that the opening of the vertical ore veins primary led to significant strain and tilt, but not to seismicity, because the delay and the longer duration of the seismicity were significant. Most seismic events involve thrusting mechanisms that are consistent with the present stress state of E-W-oriented tectonic compression, but are not consistent with the opening of the deepest ore vein. Interstingly, all the events within a few months of the heavy rainfall occurred near the faults that offset the deepest ore veins, wheareas all those events located away from the deepest ore veins occurred many months after the heavy rainfall. Consequently, the delayed diffusion of water appears to have played a dominant role in reducing rock strength, which led to seismicity in the Ikuno mine.     

Induced stresses due to fluid extraction from axisymmetric reservoirs

Earthquakes can be induced by fluid extraction, as well as by fluid injection.Segall (1989) proposed that poroelastic stresses are responsible for inducing earthquakes associated with fluid extraction. Here, I present methods for computing poroelastic stress changes due to fluid extraction for general axisymmetric reservoir geometries. The results ofGeertsma (1973) for a thin disk reservoir with uniform pressure drop are recovered as a special case. Predicted surface subsidence agrees very well with measured leveling changes over the deep Lacq gas field in southwestern France. The induced stresses are finite if the reservoir pressure changes are continuous. Computed stress changes are on the order of several bars, suggesting that the preexisting stress states in regions of extraction induced seismicity are very close to frictional instability prior to production.     

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.     

Properties of the Aftershock Sequences of the 2003 Bingöl, M D = 6.4, (Turkey) Earthquake

Aftershock sequences of the magnitude M _{W =}6.4 Bingöl earthquake of 1 May, 2003 (Turkey) are studied to analyze the spatial and temporal variability of seismicity parameters of the b value of the frequency-magnitude distribution and the p value describing the temporal decay rate of aftershocks. The catalog taken from the KOERI contains 516 events and one month’s time interval. The b value is found as 1.49 ± 0.07 with Mc =3.2. Considering the error limits, b value is very close to the maximum b value stated in the literature. This larger value may be caused by the paucity of the larger aftershocks with magnitude M _D ≥ 5.0. Also, the aftershock area is divided into four parts in order to detect the differences in b value and the changes illustrate the heterogeneity of the aftershock region. The p value is calculated as 0.86 ± 0.11, relatively small. This small p value may be a result of the slow decay rate of the aftershock activity and the small number of aftershocks. For the fitting of a suitable model and estimation of correct values of decay parameters, the sequence is also modeled as a background seismicty rate model. Constant background activity does not appear to be important during the first month of the Bingöl aftershock sequences and this result is coherent with an average estimation of pre-existing seismicity. The results show that usage of simple modified Omori law is reasonable for the analysis. The spatial variability in b value is between 1.2 and 1.8 and p value varies from 0.6 to 1.2. Although the physical interpretation of the spatial variability of these seismicity parameters is not straightforward, the variation of b and p values can be related to the stress and slip distribution after the mainshock, respectively. The lower b values are observed in the high stress regions and to a certain extent, the largest b values are related to Holocene alluvium. The larger p values are found in some part of the aftershock area although no slip occurred after the main shock and it is interpreted that this situation may be caused by the alluvium structure of the region. These results indicate that the spatial distribution in b and p values are generally related to the rupture mechanism and material properties of an aftershock area.     

Source Parameters and Scaling Relations for Reservoir Induced Seismicity in the Longtan Reservoir Area

We investigate the relationship between the impoundment and seismicity in the Longtan reservoir, southwestern China and find evidence that the seismicity was reservoir induced. After the reservoir impoundment, a pronounced increase in seismicity was observed in five clusters mainly concentrated in the areas where few earthquakes had occurred before the first filling. The observed induced seismicity shows a strong correlation with the filling cycles. The activity levels in the five clusters are different due to differences in the structures and permeabilities of the faults. Source parameters for 1,616 earthquakes with M _L 0.1–4.2 recorded by 24 fixed and temporary stations deployed around the reservoir were calculated after applying corrections for geometrical spreading, frequency-dependent Q, and site effects. The static stress drop and apparent stress in this area both appear to increase with increasing seismic moment over the entire magnitude range. Our results show that reservoir induced earthquakes have ten times lower average stress drop than natural tectonic earthquakes. These results may indicate that the reservoir induced seismicity can occur with a lower tectonic stress due to the high pore pressures of the underground medium, and that the effect of the water decreases the coefficient of friction.     

Precursory seismic quiescence along the Sumatra-Andaman subduction zone: past and present

In this study, the seismic quiescence prior to hazardous earthquakes was analyzed along the Sumatra-Andaman subduction zone (SASZ). The seismicity data were screened statistically with mainshock earthquakes of M _w?≥?4.4 reported during 1980–2015 being defined as the completeness database. In order to examine the possibility of using the seismic quiescence stage as a marker of subsequent earthquakes, the seismicity data reported prior to the eight major earthquakes along the SASZ were analyzed for changes in their seismicity rate using the statistical Z test. Iterative tests revealed that Z factors of N?=?50 events and T?=?2?years were optimal for detecting sudden rate changes such as quiescence and to map these spatially. The observed quiescence periods conformed to the subsequent major earthquake occurrences both spatially and temporally. Using suitable conditions obtained from successive retrospective tests, the seismicity rate changes were then mapped from the most up-to-date seismicity data available. This revealed three areas along the SASZ that might generate a major earthquake in the future: (i) Nicobar Islands (Z?=?6.7), (ii) the western offshore side of Sumatra Island (Z?=?7.1), and (iii) western Myanmar (Z?=?6.7). The performance of a stochastic test using a number of synthetic randomized catalogues indicated these levels of anomalous Z value showed the above anomaly is unlikely due to chance or random fluctuations of the earthquake. Thus, these three areas have a high possibility of generating a strong-to-major earthquake in the future.     

A Probabilistic Assessment of Earthquake Hazard Parameters in NW Himalaya and the Adjoining Regions

The maximum likelihood estimation method is applied to study the geographical distribution of earthquake hazard parameters and seismicity in 28 seismogenic source zones of NW Himalaya and the adjoining regions. For this purpose, we have prepared a reliable, homogeneous and complete earthquake catalogue during the period 1500–2010. The technique used here allows the data to contain either historical or instrumental era or even a combination of the both. In this study, the earthquake hazard parameters, which include maximum regional magnitude (M _max), mean seismic activity rate (λ), the parameter b (or β?=?b/log e) of Gutenberg–Richter (G–R) frequency-magnitude relationship, the return periods of earthquakes with a certain threshold magnitude along with their probabilities of occurrences have been calculated using only instrumental earthquake data during the period 1900–2010. The uncertainties in magnitude have been also taken into consideration during the calculation of hazard parameters. The earthquake hazard in the whole NW Himalaya region has been calculated in 28 seismogenic source zones delineated on the basis of seismicity level, tectonics and focal mechanism. The annual probability of exceedance of earthquake (activity rate) of certain magnitude is also calculated for all seismogenic source zones. The obtained earthquake hazard parameters were geographically distributed in all 28 seismogenic source zones to analyze the spatial variation of localized seismicity parameters. It is observed that seismic hazard level is high in Quetta-Kirthar-Sulaiman region in Pakistan, Hindukush-Pamir Himalaya region and Uttarkashi-Chamoli region in Himalayan Frontal Thrust belt. The source zones that are expected to have maximum regional magnitude (M _max) of more than 8.0 are Quetta, southern Pamir, Caucasus and Kashmir-Himanchal Pradesh which have experienced such magnitude of earthquakes in the past. It is observed that seismic hazard level varies spatially from one zone to another which suggests that the examined regions have high crustal heterogeneity and seismotectonic complexity.     

Accelerated Seismic Release and Related Aspects of Seismicity Patterns on Earthquake Faults

—Observational studies indicate that large earthquakes are sometimes preceded by phases of accelerated seismic release (ASR) characterized by cumulative Benioff strain following a power law time-to-failure relation with a term (t _f?t) ^m , where t _f is the failure time of the large event and observed values of m are close to 0.3. We discuss properties of ASR and related aspects of seismicity patterns associated with several theoretical frameworks. The subcritical crack growth approach developed to describe deformation on a crack prior to the occurrence of dynamic rupture predicts great variability and low asymptotic values of the exponent m that are not compatible with observed ASR phases. Statistical physics studies assuming that system-size failures in a deforming region correspond to critical phase transitions predict establishment of long-range correlations of dynamic variables and power-law statistics before large events. Using stress and earthquake histories simulated by the model of Ben-Zion (1996) for a discrete fault with quenched heterogeneities in a 3-D elastic half space, we show that large model earthquakes are associated with nonrepeating cyclical establishment and destruction of long-range stress correlations, accompanied by nonstationary cumulative Benioff strain release. We then analyze results associated with a regional lithospheric model consisting of a seismogenic upper crust governed by the damage rheology of Lyakhovsky et al. (1997) over a viscoelastic substrate. We demonstrate analytically for a simplified 1-D case that the employed damage rheology leads to a singular power-law equation for strain proportional to (t _f?t)^?1/3, and a nonsingular power-law relation for cumulative Benioff strain proportional to (t _f?t)^1/3. A simple approximate generalization of the latter for regional cumulative Benioff strain is obtained by adding to the result a linear function of time representing a stationary background release. To go beyond the analytical expectations, we examine results generated by various realizations of the regional lithospheric model producing seismicity following the characteristic frequency-size statistics, Gutenberg-Richter power-law distribution, and mode switching activity. We find that phases of ASR exist only when the seismicity preceding a given large event has broad frequency-size statistics. In such cases the simulated ASR phases can be fitted well by the singular analytical relation with m = ?1/3, the nonsingular equation with m = 0.2, and the generalized version of the latter including a linear term with m = 1/3. The obtained good fits with all three relations highlight the difficulty of deriving reliable information on functional forms and parameter values from such data sets. The activation process in the simulated ASR phases is found to be accommodated both by increasing rates of moderate events and increasing average event size, with the former starting a few years earlier than the latter. The lack of ASR in portions of the seismicity not having broad frequency-size statistics may explain why some large earthquakes are preceded by ASR and other are not. The results suggest that observations of moderate and large events contain two complementary end-member predictive signals on the time of future large earthquakes. In portions of seismicity following the characteristic earthquake distribution, such information exists directly in the associated quasi-periodic temporal distribution of large events. In portions of seismicity having broad frequency-size statistics with random or clustered temporal distribution of large events, the ASR phases have predictive information. The extent to which natural seismicity may be understood in terms of these end-member cases remains to be clarified. Continuing studies of evolving stress and other dynamic variables in model calculations combined with advanced analyses of simulated and observed seismicity patterns may lead to improvements in existing forecasting strategies.     

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.     

A parametric representation of Kamchatka seismicity over time

This paper presents a set of seismicity parameters that are estimated at the Kamchatka Branch of the Geophysical Service, Russian Academy of Sciences based on the regional catalog data with the purpose of routine monitoring of the current seismic situation in the region. The focus is on the identification of changes in the seismic regime (seismic quiescences and seismicity increases) in earth volumes adjacent to the maturing rupture zone of a large earthquake. The techniques we use include estimation of the seismicity level for the region using the SOUS’09 scale; calculation of the variations in the slope of the recurrence relation, identification of statistically significant anomalies in the slope using the Z test, and calculation of the seismic activity A ₁₀; monitoring the RTL parameter and variations in the area of seismogenic ruptures; using the Z test to detect areas of statistically significant decreases in the rate of seismicity; and identification of earthquake clusters. We furnish examples of such anomalies in these seismicity parameters prior to large earthquakes in Kamchatka.     

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.     

The distribution of earthquakes with depth and stress in subducting slabs

We explain the global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250–300 km depth and increases thereafter. Stress magnitude in our models shows the same pattern, in accord with experiments showing N proportional to e^kσ, with k a constant and σ stress magnitude. The models predict downdip compression in the slabs at depths below 300–400 km, as observed for earthquake stress axes.     

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.     

Tidal effects on intra-continental earthquake initiation in the Tibetan Plateau and its adjacent areas*

To investigate the tidal effects on intra-continental earthquake initiation in the Tibetan Plateau and its surrounding areas, we selected over 1,500 focal mechanism solutions of inland earthquakes (epicenter locates at least 100 km to the coastlines) from Global Centroid Moment Tensor (GCMT) project and analyzed the p-values of tidal normal and shear stress as well as tidal Coulomb failure stress. For Coulomb failure stress calculation, we used Coulomb 3.40 software. We find that: (1) p-values of tidal stress change suggests a high tidal correlation of earthquake imitations with tidal normal stress change; (2) when tidal normal stress reached the local maximum values of compression and when tidal shear stress were closed to the positive peaks, earthquakes generated more frequently; (3) particular seismogenic environments such as strong continental plate interactions and the existence of fluids or rheologic substance possibly raise the tidal correlations and (4) higher sensitivity of earthquake initiation to earth tide presents along with higher seismicity, suggesting the rate of rain energy accumulation somehow has a dominating effect on the tidal correlation of earthquake initiation.