A Study of The Impact of Static Coulomb Stress Changes of Megathrust Earthquakes Along Subduction Zone On The Following Aftershocks
A Study of The Impact of Static Coulomb Stress Changes of Megathrust Earthquakes Along Subduction Zone On The Following Aftershocks
A Study of The Impact of Static Coulomb Stress Changes of Megathrust Earthquakes Along Subduction Zone On The Following Aftershocks
net/publication/264725152
CITATIONS READS
9 163
2 authors:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Miao Miao on 22 October 2018.
Abstract Research on static earthquake triggering has been carried out widely in the world, and achieved
remarkable results. But it is still unclear whether this model is effective to all large earthquakes. In this paper,
we investigated the Coulomb stress changes of 3 megathrust earthquakes along subduction zones (the 2011 Mw 9.1
Tohoku earthquake, the 2010 Mw 8.8 Chile earthquake, and the 2004 Mw 9.0 Sumatra-Andaman earthquake) to
test the triggering effects by examining the correlation between coulomb stress increases and spatial distribution
of the following aftershocks. The calculated results suggest that there is no obvious evidence that the Coulomb
stress changes caused by the 3 megathrust earthquakes promoted the occurrence of the aftershocks. There are
only 47% of the encouraged aftershocks following the Tohoku earthquake. And there are 47.6% and 49.8% for the
Sumatra-Andaman and the Chile earthquake, respectively. We also calculated the coulomb stress changes with
different focal models and parameters. It is still less than 60% of the promoted aftershocks in the optimal case.
However, the static triggering model is good for the Wenchuan earthquake and the Chi-Chi earthquake which have
enhanced more than 85% of the subsequent aftershocks. This model may be not reasonable for large subduction
earthquakes. Therefore, other model should be introduced in studying earthquake triggering in subduction zone
and further studies will be performed.
Key words Static earthquake triggering, Coulomb failure stress change, Subduction zone, Aftershock distribu-
tion
1 INTRODUCTION
In recent years, many researchers have focused on the role of static Coulomb stress change (∆CFS) to
study the impact of mainshock on the aftershocks and interaction between strong earthquakes[1−6] . King et
al.[7] calculated the ∆CFS caused by the 1992 Landers earthquake (Mw 7.3) on the optimum failure plane and
surrounding faults. They found that aftershocks were widely distributed in the stress increasing areas while less
activities in the reduction zones. Stress changes in Southern California were modeled by Deng et al.[8] for six
large (7.0 ≤ M ≤ 7.5) to great (M ≥ 7.5) earthquakes from 1812 to 1995. About 95% of the well-located M ≥ 6
earthquakes are consistent with the Coulomb stress evolutionary model; that is, they occurred in the areas of
positive ∆CFS. And above 85% of the M ≥ 5.0 aftershocks are enhanced by the large earthquakes. Stein et
al.[9] used the mapped surface slip and fault geometry to infer the transfer of stress throughout ten M ≥ 6.7
earthquakes in the North Anatolian fault (Turkey) during 1939–1992. Calculations of the change in Coulomb
failure stress reveal that 9 out of 10 ruptures were brought closer to failure by the preceding shocks, typically
by 0.1∼1 MPa, and rearrangement of stresses in the crust leads to subsequent damaging event, the 1999 Izmit
earthquake. Wan et al.[10] explored static stress triggering problem between several complex events since 1920
in the northeastern Qinghai-Tibetan Plateau based on the formula of Okada. The result suggests that the 2001
HOh Xil earthquake is promoted by the joint ∆CFS produced by the 1937 Huashixia, the 1963 Dulan, and the
1973, 1997 Mani events. Zhang et al.[11] constructed fault models for seven earthquakes (Ms ≥ 6.0) in Jiashi
region to calculate the ∆CFS produced by each of these events. Their calculations show that the ∆CFS caused
by the preceding events was around 0.05 MPa at the hypocenter of the subsequent events. It reveals a Coulomb
stress interactive cycle of earthquake triggering between two adjacent faults. Zhu et al.[12] studied the patterns
E-mail: Joe0810@163.com
*Corresponding author: zhushoubiao@gmail.com
Miao M et al.: A Study of the Impact of Static Coulomb Stress Changes of Megathrust Earthquakes · · · 1335
of ∆CFS induced by a sequence of 4 strong earthquakes occurred in Songpan County, Sichuan province from
1973 to 1976. Most aftershocks of the sequence occurred in the areas where Coulomb stress has been increased
along the near-field zone of the seismogenic fault-segments of the main shocks. In the 25 years following the
sequence and within about 200 km centered at the mid-segment of the seismogenic Huya fault of the sequence,
6 events with magnitudes 5.0 to 6.6 all occurred in those far field zone with very small amount of increase in
Coulomb stress induced by the sequence.
The above analysis shows that static triggering model has significant influence for many activities. However,
it remains unclear whether this model is effective to all large earthquakes. Previous researches pay less attention
on the megathrust earthquakes along the subduction zone. In this paper, we use the published earthquake source
models and aftershock catalogs to investigate the ∆CFS of 3 megathrust earthquakes (the 2011 Mw 9.1 Tohoku
earthquake, the 2010 Mw 8.8 Chile earthquake and the 2004 Mw 9.0 Sumatra-Andaman earthquake) to test the
triggering effects by examining the correlation between Coulomb stress increases and spatial distribution of the
following aftershocks.
σf = τ − µ(σn + P ), (1)
where τ is the shear stress on the failure plane, σn the normal stress, P the pore fluid pressure, µ and the
coefficient of friction. Since it is hard to directly measure the real crustal stress value, the Coulomb failure
stress change is used to describe the relative change of stress[1,7] . When the coefficient of friction is constant,
Eq.(1) can be rewritten as
∆σf = ∆τ − µ(∆σn + ∆P ). (2)
Pore fluid pressure P modifies the effective normal stress across the failure plane. When rock stress is
changed more rapidly than fluid pressure diffusion, P can be related to confining stress by Skemptons coefficient
B, where the effective coefficient of friction is defined by µ0 = µ(1−B), it can account for the pore fluid behavior
and property of failure plane, with the range 0∼1, then Eq.(2) can therefore be rewritten as
where the compressive stress is positive. The ∆CFS is defined on the specific failure plane which is so called
receive fault; it promotes the occurrence of earthquake when ∆σf > 0 or vice versa. Although ∆CFS caused by
earthquake is pretty smaller than the accumulated stress required for an earthquake, lots of seismic activities
show that increase of Coulomb stress of more than 0.01 MPa appears to be sufficient to trigger events[1,14] .
During the investigation of earthquake static triggering, we often define the optimally oriented fault plane on
which the value of stress change is greater than other oriented faults in the same place if the aftershocks’ focal
mechanisms are unknown. It could explain the aftershocks distribution, and forecast the spatial pattern of
subsequent earthquakes in the future[7,15] .
3 CALCULATION
3.1 2011 Tohoku Earthquake
On March 11, 2012, the Mw 9.1 Tohoku earthquake occurred in east coast offshore of Japan’s Honshu
Island. The earthquake is located near Miyagi, northeast Japan with a depth of 20 km which typically contains
1336 Chinese J. Geophys. Vol.55, No.5
foreshock-mainshock-aftershock events. The Pacific plate collides with the Eurasia Plate in the movement
towards northwest, and subducts below Japan at an average rate of about 8 to 10 cm/a. In the formation of
the Japan Trench, volcanoes and earthquake activities take place frequently due to the stress concentration
in the plate boundary[16] . The Tohoku earthquake took place in the subduction zone, and is the largest rare
megathrust earthquake of the instrumental era in Japan. It triggered devastating tsunami and caused huge
economic losses and casualties.
The inversion of the rupture dislocation distribution using GPS displacement and seismic waveform data
as constraints was made by Wei et al.[17] . Maximum slip is estimated up to 30 m which is located in the eastward
100 km of the epicenter. The source fault generated a 500-km-long and 200-km-wide surface rupture zone which
contains totally 19×10 slip patches with a 42 km maximum rupture depth. Table 1 shows the main rupture
plane of strike 201◦ , dip 9◦ and rake 93◦ . Focal mechanisms of the Chile earthquake and Sumatra-Andaman
earthquake are also listed in Table 1.
According to the method of Stein and King et al.[2,7] , we compute stress changes in elastic half-space
with a Poisson’s ratio of 0.25, effective coefficient of friction of 0.4, and a shear modulus of 3.3×104 MPa.
Given the great magnitude of this earthquake, the duration of its influence on surrounding regions would be
a longer period. In order to comprehensively analyze the correlation between stress changes and the following
earthquakes, we choose the subsequent events before July 11, 2011 to investigate the association. Data of
M ≥ 4.0 aftershocks given by USGS (http://neic.usgs.gov/) are used. Take the optimally oriented fault as the
defined failure plane.
The Tohoku earthquake is extremely destructive which results in great impact on the Japanese mainland.
The southern part of Honshu Island is loaded by more than 0.15 MPa while more than 0.01 MPa for the far
apart Hokkaido Island. The northern Honshu Island and the source fault plane are covered by stress shadow
whose Coulomb stress greatly reduces by –0.6 to –0.08 MPa.
Numerous aftershocks are generated after the mainshock. There are more than 3000 M ≥ 4.0 events
in the first three months which are mainly distributed on the eastern and southern parts of the fault plane.
Because depth has significant influence for the ∆CFS of thrust earthquake[18] , we investigate the relationship
between stress pattern and aftershocks distribution within the whole space. The spatial distribution of ∆CFS
induced by the Tohoku earthquake with aftershocks is shown in Fig. 1, and we can discover that a large amount
of aftershocks are not triggered by the mainshock. Through the statistics of the stress change pattern with
aftershocks, about 47% events are found in the regions where the stress change is positive. On the source failure
plane there are 1160 M ≥ 4.0 aftershocks and only 40.2% of them fall in the regions of increased stress. 50%
events lie in the areas of positive stress change in the southern part of the epicenter while more than 60% fall
into the negative zone in the eastern part.
Four M ≥ 7.0 great aftershocks occurred in the following 4 months, and we find that the March 11 Mw 7.9
and July 10 Mw 7.0 events are located in the stress shadow while the March 11 Mw 7.7 and April 7 Mw 7.1 events
in the positive zone. Only 26% of shallow (< 30 km) M ≥ 6.0 aftershocks are promoted while 67% in 30 to 50
km depth. Thus the Tohoku earthquake does not have a strong pushing effect on the following events.
Miao M et al.: A Study of the Impact of Static Coulomb Stress Changes of Megathrust Earthquakes · · · 1337
it is about 115 km from Concepcion and 320 km from the Santiago, the capital of Chile. This area has
undergone several Mw ≥ 8.0 strong earthquakes since the 20th century and the largest Mw 9.5 earthquake
occurred in 1960. Ruptures caused by repeated earthquakes, largely offshore, extend nearly 2000 km parallel
to the coast[20] . Millions people were killed, at least 2 million people were affected and more than 1.5 million
houses were damaged or destroyed by the earthquake and tsunami.
Observation suggests that the Chilean quake had an asymmetrical bilateral rupture with a dominant
direction from south to north. The rupture between Santiago and Concepcion extended about 650 km long
and 180 km wide. We choose the source model given by Hayes et al.[21] to do the calculation because of its
similarity to the models developed by Pollitz and Delouis et al.[22−23] through a joint inversion of GPS and
InSAR data. The northern asperity became predominant with maximum slippage reaching about 14.5 m. This
model consists of 20×11 slip patches with maximum depth of 60 km. Take the optimally oriented fault as the
defined failure plane, and other parameters are consistent with the Tohoku model.
As shown in Fig. 2, the Coulomb stress increase areas generated by the Chile earthquake are mainly
distributed in the northwest and southeast sides of the fault plane and the decrease areas in the northeast and
southwest. The fault zones in the north and south of Chile are loaded by more than 0.05 MPa. However, the
1338 Chinese J. Geophys. Vol.55, No.5
east and west sides of the fault plane are greatly relaxed in a range –18.4 MPa to –0.05 MPa.
Take the aftershocks of 4 months also from USGS to investigate the triggering effect. The subsequent
events are concentrated on the fault plane and the area of maximum slip. According to the distribution of
∆CFS, we particularly analyze the relationship between aftershock location and stress change pattern within
depth range of 0 to 65 km: only 47.6% of the aftershocks lie in the enhanced region. The northern part of the
model generated a large shadow zone with stress reducing 5 MPa where a large number of activities occurred.
And 41% of the aftershocks happened in the focal depth are distributed in the positive region. Most of the
events fall into the suppression zone.
3.3 2004 Sumatra-Andaman Earthquake
The devastating Mw 9.0 megathrust earthquake of December 26th, 2004 occurred in the Sumatra-Andaman
subduction zone off the west coast of northern Sumatra, Indonesia. The hypocenter lied in the interface of the
India-Australia and Burma plates and was caused by the fast release of strain energy that developed as the India
plate subducts beneath the overriding Burma plate and simultaneously triggered rapid uplift of the seafloor.
The fault-rupture extends over 1200 km and lasts for more than 1000 s. This is the largest earthquake in the
world during the last 40 years and triggered the most destructive tsunami in record that directly resulted nearly
0.3 million deaths.
The subduction zone manifests low-angle in shallow and high-angle in the deep. The depth of the Sumatra
earthquake is 28.6 km, consistent with low dip angle
of nodal plane I in Table 1. The dislocation model
published by Chlieh et al.[24] is adopted for the cal-
culation. This complex model is divided into three
sub-faults with maximum slippage of 17 m that are
located around the hypocenter and the confluence
of the Andaman Trench and Sumatra fault. Strikes
of the three sub-faults are 5◦ , 342◦ , 315◦ and dips
17◦ , 15◦ , 12◦ from north to south. There are totally
695 slip patches with spatially asymmetrical distri-
bution, and rupture extends to 52.6 km in depth.
∆CFS distribution pattern generated by the
Sumatra earthquake is similar to the Chile earth-
quake’s (Fig. 3). The loading areas are mainly lo-
cated in the northwest and southeast of the fault
zone as well as the fault plane. The southern seg-
ments of Sumatra fault are loaded by more than
0.03 MPa and these results have certain degree of
consistency with McCloskey’s study[25] . Although
there are 5 more strong earthquakes (Mw >7.0) oc-
curred in this region after the megathrust earth-
quake, they are separated by four years and three
of them lie in the negative area. Stress change re-
sults calculated by the full elastic half space model
cannot be a good description for the triggering ef-
fect.
Subsequent aftershocks are widely distributed
throughout the fault plane and plenty of them lie Fig. 3 Coulomb stress change of the Sumatra-
in the middle of the Andaman Trench (black box in Andaman earthquake and the spatial
Fig. 3). More than 70% aftershocks take place in distribution of the aftershocks
Miao M et al.: A Study of the Impact of Static Coulomb Stress Changes of Megathrust Earthquakes · · · 1339
30 km depth. In order to avoid the impact of the Mw 8.6 great earthquake (March 28, 2005) in this area, only
the events in subsequent 3 months (to March 26, 2005) would be taken into consideration. Although other
study pointed out that the Mw 8.6 strong earthquake is promoted by the megathrust earthquake[26] , we don’t
regard it as aftershock for the great magnitude.
Figure 3 introduces the relationship between spatial distribution of the aftershocks and ∆CFS pattern. It
is shown that about 49.8% of the aftershocks occur in the region where the stress is calculated to increase, the
percentage drops to 46% on the source failure plane. Numerous activities lie in the central section of Andaman
trench which is the boundary of stress raise and decline, and 60% of them are found where the stress change is
negative (–1 to –0.1 MPa). Majority of the aftershocks near the Sumatra fault fall in the failure plane and the
number of events in the region of stress increase is not significantly more than the number in the shadow zone.
No remarkable correlation is found between static stress changes produced by the mainshock and the following
quakes.
4 COMPARISON ANALYSIS
Static triggering effect is unsatisfactory for the three magathrust earthquake calculated above. Most of
the areas on the fault plane are the reduction zones which go against the occurrence of the subsequent seismic
activities while a great many of aftershocks distribute in these parts. Less than 50% of the seismic activities are
promoted by the three megathrust earthquakes. Mainshock stress change patterns have no obvious relationship
with aftershocks’ occurrence. To avoid the calculation error from the unique model and parameters, we select
different source models and parameters to analyze comparatively.
The unexpected Wenchuan earthquake occurred on the Longmenshan fault which is not active recently
in the eastern margin of the Tibetan Plateau, the middle of the North-south Seismic Belt. Previous studies
indicated that the Beichuan-Yingxiu fault along the Longmenshan thrust belt was the main source fault, which
generated a 240-km-long surface rupture zone, and its maximum vertical and horizontal slips were estimated up
to 6.2 and 4.9 m respectively[31] . An additional surface rupture zone developed along the Guanxian-Jiangyou
fault with a 72-km-long surface rupture and 3.5 m maximum horizontal slip. Great shock leads to landslides
and other geological disasters, besides, huge casualties and property losses.
The Wenchuan earthquake had a complicated
rupture process; Yingxiu and Beichuan were the most
damaged regions with high-slip areas. Ji et al. pub-
lished the slip distribution model seven hours after the
mainshock: This earthquake was mainly a thrust event
with right-lateral strike-slip component in the southern
section and mainly right-lateral strike-slip with thrust
component in the northern section. The model consists
of 21×8 slip patches with maximum rupture depth of
20 km.
We use the conditions above to calculate the
∆CFS, and the pattern is shown in Fig. 6. The spatial
distribution of stress change is intricate that alternates
positive and negative. Maximum increasing zone is lo-
cated in the north and south lobes of the fault plane
and maximum decreasing zone in the east and west
sides.
Over 1000 aftershocks happened within four Fig. 6 Coulomb stress change of the Wenchuan earth-
months after the Wenchuan earthquake, and are mostly quake and the spatial distribution of the aftershocks
concentrated in west side of Beichuan-Yingxiu fault, and a few aftershocks are distributed around vicinity of
Wenchuan-Yingxiu.
Figure 6 shows the spatial distribution of ∆CFS induced by the Wenchuan earthquake and its aftershocks.
In this model, more than 85% of aftershocks lie where the stress is calculated to increase, the percentage rises to
87% on the source failure plane. The stress variation is greater than triggering threshold 0.01 MPa. This result
indicates that the Wenchuan earthquake plays a good role in triggering the subsequent activities. Compared
with the aftershocks occurred in the first month (May 12,2008 to June 12, 2008), events in the second to
fourth month (June 13, 2008 to September 12, 2008) reduced sharply, and more than 90% of them fall into the
increasing zone. This result is consistent with other studies[32−33] and the triggering effect of the Wenchuan
earthquake would last more than six months in the study area.
According to previous studies, it is found that intraplate earthquake usually has a good triggering effect:
more than 20000 aftershocks occurred after the 1999 Chi-Chi earthquake (Mw 7.6), and the calculation result
shows that in shallow depths ∆CFS generated by the mainshock promote most of the subsequent events using
a half-space elastic model[34] . When it comes to multi-layered medium model there are about 72% of the
aftershocks in the region of Coulomb stress increasing over 0.01 MPa[35] . The study of Ma et al.[36] shows that
there is a widespread increase in seismicity rate in the surrounding areas. Besides, the earthquakes mentioned
above which gain good triggering effect are all intraplate earthquakes.
We can see from the comparison analysis above that regardless of the variation of the source models and
the parameters, the triggering effect of the three megathrust earthquakes is insignificant. In the optimal case
the promoted aftershocks are still less than 60%. However, the static triggering model is good for intraplate
earthquakes (the 2008 Wenchuan earthquake, the 1999 Chi-Chi earthquake) which have enhanced more than
1342 Chinese J. Geophys. Vol.55, No.5
can reduce the strength of the fault, resulting in fault rupture or reactivation of existing rupture with constant
shear stress, and finally earthquake; cracks and traps formed after the mainshock play strong control roles on
the pore pressure. Thermal expansion of the underground fluid would generate lots of gaps; these gaps become
good path for accumulation and migration of fluid to impact on the stability of the fault in return. For the West
Bohemia/Vogtland earthquake swarms, increased pore pressure of crustal fluids in the region takes a key part
in bringing the faults from the subcritical to critical state. The swarm activities are mainly driven by stress
changes due to co-seismic and post-seismic slips, which considerably depend on the frictional conditions at the
fault; crustal fluids keep the fault in a critical state[47] . Brodsky et al.[48] propose a new model for coseismic
pore pressure steps in which a temporary barrier deposited by ground water flow is entrained and removed by
the more rapid flow induced by the seismic waves. In hydrothermal areas, this mechanism could lead to 0.04
MPa pressure changes and trigger seismicity. Astiz et al.[49] also point out that the following-up activities of
the 1990 California Upland ML 5.5 earthquake can be related to the fluid-rich hypocenter.
(3) Different types of models. Earthquake occurrence is mainly a process of strain energy release, and
thus in small time scales we often simplify the earth as isotropic elastic half-space medium. This is a certain
approximation for stress transfer in short period after the mainshock. In order to get closer to the actual
situation, some scholars often adopted multi-layered medium model for the calculation of ∆CFS. Take Wang
and Shao’s researches[34−35] for example, both models have been done for the ∆CFS and the results reveal that
stress change patterns are consistent in general but some obvious differences in details cannot be ignored.
Afterslip generated by a great earthquake often gives rise to stress state perturbation of the hypocenter
and surrounding areas. It often acts in a short-term with rapid decay (a few months to one year). But it is
still difficult to figure out how much the afterslip would work (different afterslips may effect differently). So
far, the stress change evolution problem considering the afterslip has not been reported. The afterslip model
in subduction zone is hard to establish because of the lack of near-field surface deformation information. In
this paper we only considered the co-seismic ∆CFS. Stress perturbation from the afterslip would be taken into
consideration in further research.
For large time scales viscoelastic relaxation stress evolution is worth to explore due to the coupling of the
brittle upper-crust and viscoelastic lithosphere. Shen et al.[50] found that the ∆CFS generated by the East
Kunlun earthquake increase gradually with time because of the effect of viscoelastic relaxation in the lower
crust. And the ∆CFS caused by the viscoelastic relaxation can even greatly exceed the co-seismic ∆CFS.
Therefore, although energy accumulation is the key factor of the strong earthquake occurrence in the plate
boundary within a short period, crustal viscoelastic relaxation also plays a vital role for the subsequent events
close to the mainshock.
(4) The influence of the initial stress field. In the studies of static earthquake triggering, an empirical
threshold of ∆CFS (say, 0.01 MPa) is often taken to determine the risk of surrounding areas. In fact, ∆CFS> 0
only means the probability increase of the earthquake occurrence. The earthquake occurrence is related not only
to ∆CFS but also to the stresses before the mainshock[51] . On the one hand, the earthquake may be triggered
even if the ∆CFS is slightly greater than zero if the stresses before the mainshock are very close to Coulomb
failure. On the other hand, it may not be triggered even though ∆CFS is large if the stresses are far away from
the Coulomb failure. Hardbeck et al.[52] pointed out that the aftershocks in one month of the 1994 Northridge
earthquake can not be explained by the ∆CFS. The impact of the tectonic stress field and fault strength should
be considered. Hu et al.[51] used continuous evolution of stress field to investigate the triggering effect of the
1976 Tangshan earthquake. The results indicate that the effects of the stress field of the Tangshan earthquake
fault zone on the aftershock triggering are very important. The aftershocks fallen in the earthquake triggering
regions predicted by the new method are more than those fallen in the regions of ∆CFS≥ 0 predicted by seismic
dislocation theory.
(5) Specificity of the megathrust earthquake. Mechanism of the devastating earthquake in subduction
zone has some differences with intraplate earthquake. Due to the special geographical location, the seismogenic
1344 Chinese J. Geophys. Vol.55, No.5
environment, stress state and tectonic background are unknown, the nature specificity of the subduction zone
may cause the static triggering failure. In addition, Lin et al.[18] believed that when the fault aspect ratio L/W
is large, Coulomb stress changes would generate a large “earthquake shadow” on both sides of the fault. The
megathrust earthquake models are relatively long and narrow, such as the rupture of Sumatra earthquake is up
to 1000 km while the width is only 300 km. The large aspect ratio may be one reason of the reduction of the
triggered aftershocks.
Besides the three megathrust earthquakes there are several other examples of bad triggering effect on the
subduction zone. Such as the aftershock sequence of the 1989 Loma Prieta earthquake in Southern California has
no significant association with Coulomb stress change pattern[41] . The 1911 Morgan Hill earthquake’s hypocen-
ter falls in the 1906 San Francisco earthquake’s stress shadow zone[53] . Consequently, the static earthquake
triggering is a very complicated issue, and need further detailed study.
ACKNOWLEDGMENTS
Acknowledgement We thank two reviewers for their valuable suggestions. We are grateful to Jian Lin,
Ross Stein and Shinji Toda for their data and kind help. This study was supported by National Natural Science
Foundation of China (40974020).
REFERENCES
[1] Harris R A. Introduction to special section: Stress triggers, stress shadows, and implications for seismic hazard. J.
Geophys. Res., 1998, 103(B10): 24347-24358.
[2] Stein R S. The role of stress transfer in earthquake occurrence. Nature, 1999, 402(6762): 605-609.
[3] King G C P, Cocco M. Fault interaction by elastic stress changes: New clues from earthquake sequences. Advances
in Geophysics, 2001, 44: 1-38, I-VIII.
[4] Freed A M. Earthquake triggering by static, dynamic, and postseismic stress transfer. Annual Review of Earth and
Planetary Sciences, 2005, 33(1): 335-367.
[5] Toda S, Stein R S, Richards-Dinger K, et al. Forecasting the evolution of seismicity in southern California: Animations
built on earthquake stress transfer. J. Geophys. Res., 2005, 110(B5): B05S16.
[6] Shi Y L, Cao J L. Some aspects in static stress change calculation–case study on Wenchuan earthquake. Chinese J.
Geophys. (in Chinese), 2010, 53(1): 102-110.
[7] King G C P, Stein R S, Lin J. Static stress changes and the triggering of earthquakes. Bulletin of the Seismological
Society of America, 1994, 84(3): 935-953.
[8] Deng J S, Sykes L R. Evolution of the stress field in southern California and triggering of moderate-size earthquakes:
A 200-year perspective. J. Geophys. Res., 1997, 102(B5): 9859-9886.
[9] Stein R S, Barka A A, Dieterich J H. Progressive failure on the North Anatolian fault since 1939 by earthquake stress
triggering. Geophysical Journal International, 1997, 128(3): 594-604.
Miao M et al.: A Study of the Impact of Static Coulomb Stress Changes of Megathrust Earthquakes · · · 1345
[10] Wan Y G, Wu Z L, Zhou G W, et al. “Stress triggering” between different rupture events in several earthquakes.
Acta Seismologica Sinica (in Chinese), 2000, 22(6): 568-576.
[11] Zhang Z Q, Chen Y S, Lin J. Stress interaction between normal fault and strike-slip fault in Jiashi earthquake group
1997. Sci. China Earth Sci. (in Chinese), 2008, 38(3): 334-342.
[12] Zhu H, Wen X Z. Stress triggering process of the 1973 to 1976 Songpan, Sichuan, sequence of strong earthquakes.
Chinese J. Geophys. (in Chinese), 2009, 52(4): 994-1003.
[13] Okada Y. Internal deformation due to shear and tensile faults in a half-space. Bulletin of the Seismological Society
of America, 1992, 82(2): 1018-1040.
[14] Ziv A, Rubin A M. Static stress transfer and earthquake triggering: No lower threshold in sight? J. Geophys. Res.,
2000, 105(B6): 13631-13642.
[15] Stein R S, Lisowski M. The 1979 Homestead Valley earthquake sequence, California: control of aftershocks and
postseismic deformation. J. Geophys. Res., 1983, 88(B8): 6477-6490.
[16] DeMets C. Oblique convergence and deformation along the Kuril and Japan Trenches. J. Geophys. Res., 1992,
97(B12): 17615-17625.
[17] Wei S, Sladen A. Slip-History database: 2011 Tohoku-oki Earthquake. 2011: http://www.tectonics.caltech.edu/slip
history/2011 taiheiyo-oki/index.html.
[18] Lin J, Stein R S. Stress triggering in thrust and subduction earthquakes and stress interaction between the southern
San Andreas and nearby thrust and strike-slip faults. J. Geophys. Res., 2004, 109(B2): B02303.
[19] Barazangi M, Isacks B L. Spatial distribution of earthquakes and subduction of the Nazca plate beneath South
America. Geology, 1976, 4(11): 686-692.
[20] Kelleher J A. Rupture zones of large South American earthquakes and some predictions. J. Geophys. Res., 1972,
77(11): 2087-2103.
[21] Hayes G. Finite Fault Model, Updated Result of the Feb 27, 2010 Mw 8.8 Maule, Chile Earthquake 2010: http://earth-
quake.usgs.gov/earthquakes/eqinthenews/2010/us2010tfan/finite fault.php.
[22] Pollitz F F, Brooks B, Tong X P, et al. Coseismic slip distribution of the February 27, 2010 Mw 8.8 Maule, Chile
earthquake. Geophys. Res. Lett., 2010, 38(9): L09309.
[23] Delouis B, Nocquet J M, Vallée M. Slip distribution of the February 27, 2010 Mw =8.8 Maule Earthquake, central
Chile, from static and high-rate GPS, InSAR, and broadband teleseismic data. Geophys. Res. Lett., 2010, 37(17):
L17305.
[24] Chlieh M. Tectonics Observatory at Caltech Slip History database: Sumatra Earthquake, update 1. 2007: http://www.
tectonics.caltech.edu/slip history/2004 sumatra/sumatra.html.
[25] McCloskey J, Nalbant S S, Steacy S. Indonesian earthquake: Earthquake risk from co-seismic stress. Nature, 2005,
434(7031): 291-291.
[26] Pollitz F F, Banerjee P, Bürgmann R, et al. Stress changes along the Sunda trench following the 26 December 2004
Sumatra-Andaman and 28 March 2005 Nias earthquakes. Geophys. Res. Lett., 2006, 33(6): L06309.
[27] Shao G, Li X, Ji C, et al. Preliminary Result of the Mar 11, 2011 Mw 9.1 Honshu Earthquake. 2011: http://www.geol.
ucsb.edu/faculty/ji/ big earthquakes/2011/03/0311/Honshu main.html.
[28] Shao G, Li X, Liu Q, et al. Preliminary slip model of the Feb 27, 2010 Mw 8.9 Maule, Chile Earthquake. 2010:
http://www.geol.ucsb.edu/faculty/ji/big earthquakes/2010/02/27/chile 2 27.html.
[29] Toda S, Lin J, Meghraoui M, et al. 12 May 2008 M = 7.9 Wenchuan, China, earthquake calculated to increase failure
stress and seismicity rate on three major fault systems. Geophys. Res. Lett., 2008, 35(17): L17305.
[30] Ji C. Preliminary Result of the May 12, 2008 Mw 7.97 SiChuan Earthquake. 2008: http://www.geol.ucsb.edu/faculty/
ji/big earthquakes/2008/05/12/ShiChuan.html.
[31] Zhu S B, Zhang P Z. A study on the dynamical mechanisms of the Wenchuan MS 8.0 earthquake, 2008. Chinese J.
Geophys. (in Chinese), 2009, 52(2): 418-427.
[32] Parsons T, Ji C, Kirby E. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan
basin. Nature, 2008, 454(7203): 509-510.
[33] Xie C D, Zhu Y Q, Lei X L, et al. Pattern of stress change and its effect on seismicity rate caused by Ms 8.0 Wenchuan
earthquake. Sci. China Earth Sci. (in Chinese), 2010, 40(6): 688-698.
[34] Wang J C, Shieh C F, Chang T M. Static stress changes as a triggering mechanism of a shallow earthquake: case
study of the 1999 Chi-Chi (Taiwan) earthquake. Physics of the Earth and Planetary Interiors, 2003, 135(1): 17-25.
1346 Chinese J. Geophys. Vol.55, No.5
[35] Shao Z G, Fu R S, Xue T H, et al. Discussion about the relation of coulumb stress changes and distribution of
aftershocks for the case of the chi-chi earthquake. Progress in Geophys. (in Chinese), 2009, 24(2): 367-374.
[36] Ma K F, Chan C H, Stein R S. Response of seismicity to Coulomb stress triggers and shadows of the 1999 Mw = 7.6
Chi-Chi, Taiwan, earthquake. J. Geophys. Res., 2005, 110(B5): B05S19.
[37] Zhan Z, Jin B, Wei S, et al. Coulomb stress change sensitivity due to variability in mainshock source models and
receiving fault parameters: A case study of the 2010–2011 christchurch, New Zealand, Earthquakes. Seismological
Research Letters, 2011, 82(6): 800-814.
[38] D´lez M, La Femina P C, Connor C B, et al. Evidence for static stress changes triggering the 1999 eruption of Cerro
Negro Volcano, Nicaragua and regional aftershock sequences. Geophys. Res. Lett., 2005, 32(4): L04309.
[39] Mueller K, Hough S E, Bilham R. Analysing the 1811-1812 New Madrid earthquakes with recent instrumentally
recorded aftershocks. Nature, 2004, 429(6989): 284-288.
[40] Stramondo S, Kyriakopoulos C, Bignami C, et al. Did the September 2010 (Darfield) earthquake trigger the February
2011 (Christchurch) event? Sci. Rep., 2011, 1, doi:10.1038/srep00098.
[41] Kilb D, Gomberg J, Bodin P. Triggering of earthquake aftershocks by dynamic stresses. Nature, 2000, 408(6812):
570-574.
[42] Hori T, Kaneda Y. A simple explanation for the occurrence of the 1911 Morgan Hill Earthquake in the stress shadow
of the 1906 San Francisco Earthquake. Geophys. Res. Lett., 2001, 28(11): 2261-2264.
[43] Parsons T. A hypothesis for delayed dynamic earthquake triggering. Geophys. Res. Lett., 2005, 32(4): L04302.
[44] Tibi R, Wiens D A, Inoue H. Remote triggering of deep earthquakes in the 2002 Tonga sequences. Nature, 2003,
424(6951): 921-925.
[45] Wu X P, Hu X L, Bouchon M, et al. Complete coulomb stress change and aftershocks stress triggering of Ms 7.6
Lancang-Gengma Yunnan. Sci. China Earth Sci. (in Chinese), 2007, 37(6): 746-752.
[46] Miao M, Zhu S B. The impact of the underground fluid in the process of earthquake preparation and occurrence: A
review. Progress in Geophys. (in Chinese), 2012, 27(3): 950-959.
[47] Horálek J, Fischer T. Role of crustal fluids in triggering the West Bohemia/Vogtland earthquake swarms: Just what
we know (a review). Studia Geophysica et Geodaetica, 2008, 52(4): 455-478.
[48] Brodsky E E, Roeloffs E, Woodcock D, et al. A mechanism for sustained groundwater pressure changes induced by
distant earthquakes. J. Geophys. Res., 2003, 108(B8): 2390.
[49] Astiz L, Shearer P M, Agnew D C. Precise relocations and stress change calculations for the Upland earthquake
sequence in southern California. J. Geophys. Res., 2000, 105(B2): 2937-2953.
[50] Shen Z K, Wan Y G, Gan W J, et al. Viscoelastic triggering among large earthquakes along the east Kunlun fault
system. Chinese J. Geophys. (in Chinese), 2003, 46(6): 786-795.
[51] Hu C B, Zhou Y J, Cai Y E. A new finite element model in studying earthquake triggering and continuous evolution
of stress field. Sci. China Earth Sci. (in Chinese), 2009, 39(5): 546-555.
[52] Hardebeck J L, Nazareth J J, Hauksson E. The static stress change triggering model: Constraints from two southern
California aftershock sequences. J. Geophys. Res., 1998, 103(B10): 24427-24437.
[53] Gomberg J, Blanpied M L, Beeler N M. Transient triggering of near and distant earthquakes. Bulletin of the
Seismological Society of America, 1997, 87(2): 294-309.