Magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> in metropolitan France

The recent seismicity catalogue of metropolitan France Sismicité Instrumentale de l’Hexagone (SI-Hex) covers the period 1962–2009. It is the outcome of a multipartner project conducted between 2010 and 2013. In this catalogue, moment magnitudes (M _w) are mainly determined from short-period velocimetric records, the same records as those used by the Laboratoire de Détection Géophysique (LDG) for issuing local magnitudes (M _L) since 1962. Two distinct procedures are used, whether M _L-LDG is larger or smaller than 4. For M _L-LDG >4, M _w is computed by fitting the coda-wave amplitude on the raw records. Station corrections and regional properties of coda-wave attenuation are taken into account in the computations. For M _L-LDG ≤4, M _w is converted from M _L-LDG through linear regression rules. In the smallest magnitude range M _L-LDG <3.1, special attention is paid to the non-unity slope of the relation between the local magnitudes and M _w. All M _w determined during the SI-Hex project is calibrated according to reference M _w of recent events. As for some small events, no M _L-LDG has been determined; local magnitudes issued by other French networks or LDG duration magnitude (M _D) are first converted into M _L-LDG before applying the conversion rules. This paper shows how the different sources of information and the different magnitude ranges are combined in order to determine an unbiased set of M _w for the whole 38,027 events of the catalogue.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

The 20-second regional magnitude <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript>(20R) for the Russian Far East

A new modified magnitude scale M _S (20R) is elaborated. It permits us to extend the teleseismic magnitude scale M _S (20) to the regional epicenter distances. The data set used in this study contains digital records at 12 seismic stations of 392 earthquakes that occured in the northwest Pacific Ocean in the period of 1993–2008. The new scale is based on amplitudes of surface waves of a narrow range of the periods (16–25 s) close to the period of 20 s, for distances of 80–3000 km. The digital Butterworth filter is used for processing. On the basis of the found regional features concerning distance dependence for seismic wave attenuation, all the stations of the region have been subdivided into two groups, namely, “continental” and “island-arc.” For each group of stations, its own calibration function is proposed. Individual station corrections are used to compensate for the local features.     

Rupture behavior and deformation localization of the Kunlunshan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>7.8) and their tectonic implications

Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M _w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M _w =6.8, M _w =6.2 and M _w ⩽=7.8 events, respectively. The M _w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)     

The earthquake of July 22, 2011 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 4.5) in a low-seismicity area of the Argun region

The paper considers the Argun earthquake of July 22, 2011 (M _w = 4.5), which occurred in the Argun River valley in a low-seismicity territory in China. The focal parameters of the earthquake (depth of the hypocenter, moment magnitude, scalar seismic moment, and focal mechanism) were determined by calculating the seismic moment tensor from the amplitude spectra of surface waves and the data on the signs of the first arrivals of body waves at regional stations. The solution of the focal mechanism makes it possible to assume a relationship between the earthquake focus and a fault with a northeastern strike bordering the southeastern side of the Argun Basin (in Chinese territory). The Argun earthquake was felt in Russia with an intensity of II–III to V at the epicentral distances up to 255 km. The intensity of shaking did not exceed values suggested by new GSZ-2012 and GSZ-2014 seismic zoning maps of Russian territory. Nevertheless, the question on the possible occurrence of stronger earthquakes in the studied region remains open.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> and <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript> comparison for earthquakes in Greece for the period 1969 to 2007

We use 576 earthquakes of magnitude, M _w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M _w, and local magnitude, M _L, for earthquakes in Greece and the surrounding regions. We compare M _w to M _L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M _L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M _w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M _L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M _w and M _L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M _L’s. We propose scaling relations to convert the M _L values reported from the two centers to M _w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.     

Stochastic finite-fault simulation of ground motion from the August 11, 2012, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.4 Ahar earthquake,northwestern Iran

In this study, the 11 August 2012 M _w 6.4 Ahar earthquake is investigated using the ground motion simulation based on the stochastic finite-fault model. The earthquake occurred in northwestern Iran and causing extensive damage in the city of Ahar and surrounding areas. A network consisting of 58 acceleration stations recorded the earthquake within 8–217 km of the epicenter. Strong ground motion records from six significant well-recorded stations close to the epicenter have been simulated. These stations are installed in areas which experienced significant structural damage and humanity loss during the earthquake. The simulation is carried out using the dynamic corner frequency model of rupture propagation by extended fault simulation program (EXSIM). For this purpose, the propagation features of shear-wave including \( {Q}_s \) value, kappa value \( {k}_0 \), and soil amplification coefficients at each site are required. The kappa values are obtained from the slope of smoothed amplitude of Fourier spectra of acceleration at higher frequencies. The determined kappa values for vertical and horizontal components are 0.02 and 0.05 s, respectively. Furthermore, an anelastic attenuation parameter is derived from energy decay of a seismic wave by using continuous wavelet transform (CWT) for each station. The average frequency-dependent relation estimated for the region is \( Q=\left(122\pm 38\right){f}^{\left(1.40\pm 0.16\right)}. \) Moreover, the horizontal to vertical spectral ratio \( H/V \) is applied to estimate the site effects at stations. Spectral analysis of the data indicates that the best match between the observed and simulated spectra occurs for an average stress drop of 70 bars. Finally, the simulated and observed results are compared with pseudo acceleration spectra and peak ground motions. The comparison of time series spectra shows good agreement between the observed and the simulated waveforms at frequencies of engineering interest.     

Water-level changes in the wells of Kamchatka at the time of the <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.6, April 20, 2006 Olyutorskii earthquake

The seismic waves excited by the M _w 7.6 Olyutorskii earthquake that occurred on April 20, 2006 in the Koryak Upland gave rise to water-level changes in five wells situated in continental areas of Kamchatka at hypocentral distances of 750–1150 km. We describe the effects due to seismic waves, as well as the water-level anomalies for February–April 2006 before the earthquake. We used an original technique for the processing of water-level records based on the study of barometric and tidal water-level responses in order to estimate the volume strain in water-saturated rocks during synchronous level variations at two wells. We discuss possible mechanisms for producing anomalous water-level changes due to elastic deformation of monitored groundwater reservoirs and to crack dilatancy in the water-saturated rocks.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

On the Ultra Long Propagation of Felt Ground Motion Due to the <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 8.3 Deep-Focus Sea-of-Okhotsk Earthquake of May 24, 2013

This study uses macroseismic data and wave equations to solve the problem of ultra long propagation of felt ground motion (over 9000 km from the epicenter) due to the Sea-of-Okhotsk earthquake. We show that the principal mechanism of this phenomenon could be excitation of a previously unknown standing radial wave as a mode of the Earth’s free oscillations, ₀S₀, due to the superposition of an incident and a reflected spherical P wave in the epicentral area of the Sea-of-Okhotsk earthquake. The standing wave generates slowly attenuating P waves that travel over the earth’s surface that act as carrying waves; when superposed on these, direct body waves acquire the ability to travel over great distances. We show previously unknown parameters of the radial mode ₀S₀ for the initial phase of earth deformation due to the large deep-focus earthquake. We used data on the Sea-of-Okhotsk and Bolivian earthquakes to show that large deep-focus earthquakes can excite free oscillations of the Earth that are not only recorded by instrumental means, but are also felt by people, with the amplification of the macroseismic effect being directly related to the phenomenon of resonance for multistory buildings.     

Application of Fisher method to discriminating earthquakes and explosions using criterion <Emphasis Type="Italic">m</Emphasis><Subscript>b</Subscript>/<Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>

We try to give a quantitative and global discrimination function by studying m _b/M _S data using Fisher method that is a kind of pattern recognition methods. The reliability of the function is also analyzed. The results show that this criterion works well and has a global feature, which can be used as first-level filtering criterions in event identification. The quantitative and linear discrimination function makes it possible to identify events automatically and achieve the goal to react the events quickly and effectively.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

Dispersion Measurements of <Emphasis Type="Italic">P</Emphasis> Waves and their Implications for Mantle <Emphasis Type="Italic">Q</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript>

We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q _p ), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q _p with an interface at 660 km depth and take Q _p constant with period. Our data set is too small to solve for a possible frequency dependence of Q _p . The upper mantle Q ₁ is 476 [299–1176] and the lower mantle Q ₂ is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q _p ^–1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135.     

Rupture Directivity During the September 7, 1999 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 5.9) Athens (Greece) Earthquake Inferred from Forward Modeling of Strong Ground Motion

The empirical Green's functions technique is applied to simulate strong ground motion records from the September 7, 1999, Athens earthquake. Information on the fault parameters from previous independent studies has been used and several scenarios were examined, in regard to the location of the starting point of the rupture, by comparing the synthetic records with the corresponding observed ones, through a residual function and a correlation function. The results show that the rupture started at the deepest, 4–5 km, part of the fault from its western edge. This hypocenter was then used, in combination with the initial fault model, to stochastically simulate the strong ground motion during the Athens main shock, in terms of peak-ground acceleration at hard rock. The results show that directivity might have significantly contributed to the destructiveness of this earthquake at specific parts of the Athens Metropolitan area.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

Slope topography-induced spatial variation correlation with observed building damages in Corso during the May 21, 2003, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.8, Boumerdes earthquake (Algeria)

During the May 21, 2003 M _w 6.8 Boumerdes earthquake, in the “Cité des 102 Logements” built on a hilltop, in Corso, heavy damages were observed: near the crest, a four-story RC building collapsed while others experienced severe structural damage and far from the crest, slight damage was observed. In the present paper, we perform a 2D slope topography seismic analysis and investigate its effects on the response at the plateau as well as the correlation with the observed damage distribution. A site-specific seismic scenario is used involving seismological, geological, and geotechnical data. 2D finite element numerical seismic study of the idealized Corso site subjected to vertical SV wave propagation is carried out by the universal code FLUSH. The results highlighted the main factors that explain the causes of block collapse, located 8-26 m far from the crest. These are as follows: (i) a significant spatial variation of ground response along the plateau due to the topographic effect, (ii) this spatial variation presents high loss of coherence, (iii) the seismic ground responses (PGA and response spectra) reach their maxima, and (iv) the fundamental frequency of the collapsed blocks coincides with the frequency content of the topographic component. For distances far from the crest where slight damages were observed, the topographic contribution is found negligible. On the basis of these results, it is important to take into account the topographic effect and the induced spatial variability in the seismic design of structures sited near the crest of slope.     

The <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 4.3 January 17, 2014, earthquake: very rare seismic event on the Siberian platform

We described the 2014 January 17 earthquake (M _w?=?4.3) occurred on the Siberian platform in the area of sublongitudinal part of the Angara river in the zone of possible influence of two large reservoirs—Ust-Ilimsk and Boguchan. This is the first event of such magnitude recorded in this previously aseismic area during the whole period of instrumental observations. A seismic moment, a moment magnitude, a hypocentral depth, and a focal mechanism of the event were calculated on the basis of surface wave amplitude spectra. Analysis of the geological and geophysical data showed that the earthquake origin is connected with high-velocity gradient zone located at the border of the Late Proterozoic cover and the Precambrian basement of the Siberian platform. Some evidences for a natural character of the earthquake were considered.     

The 2014 Kefalonia Doublet (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>6.1 and <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>6.0), Central Ionian Islands,Greece: Seismotectonic Implications along the Kefalonia Transform Fault Zone

The 2014 Kefalonia earthquake sequence started on 26 January with the first main shock (M_W6.1) and aftershock activity extending over 35 km, much longer than expected from the causative fault segment. The second main shock (M_W6.0) occurred on 3 February on an adjacent fault segment, where the aftershock distribution was remarkably sparse, evidently encouraged by stress transfer of the first main shock. The aftershocks from the regional catalog were relocated using a 7-layer velocity model and station residuals, and their distribution evidenced two adjacent fault segments striking almost N-S and dipping to the east, in full agreement with the centroid moment tensor solutions, constituting segments of the Kefalonia Transform Fault (KTF). The KTF is bounded to the north by oblique parallel smaller fault segments, linking KTF with its northward continuation, the Lefkada Fault.     

Structural features and seismotectonic implications of coseismic surface ruptures produced by the 2016 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.1 Kumamoto earthquake

Field investigations and analyses of satellite images and aerial photographs reveal that the 2016 M _w 7.1 (Mj 7.3) Kumamoto earthquake produced a ～40-km surface rupture zone striking NE-SW on central Kyushu Island, Japan. Coseismic surface ruptures were characterized by shear faults, extensional cracks, and mole tracks, which mostly occurred along the pre-existing NE-SW-striking Hinagu–Futagawa fault zone in the southwest and central segments, and newly identified faults in the northeast segment. This study shows that (i) the Hinagu–Futagawa fault zone triggered the 2016 Kumamoto earthquake and controlled the spatial distribution of coseismic surface ruptures; (ii) the southwest and central segments were dominated by right-lateral strike-slip movement with a maximum in-site measured displacement of up to 2.5 m, accompanied by a minor vertical component. In contrast, the northeast segment was dominated by normal faulting with a maximum vertical offset of up to 1.75 m with a minor horizontal component that formed graben structures inside Aso caldera; (iii) coseismic rupturing initiated at the jog area between the Hinagu and Futagawa faults, then propagated northeastward into Aso caldera, where it terminated. The 2016 M _w 7.1 Kumamoto earthquake therefore offers a rare opportunity to study the relationships between coseismic rupture processes and pre-existing active faults, as well as the seismotectonics of Aso volcano.     

Stress changes on major faults caused by <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>7.9 Wenchuan earthquake,May 12, 2008

On May 12, 2008, a magnitude 7.9 earthquake ruptured the Longmenshan fault system in Sichuan Province, China, collapsing buildings and killing tens of thousands people. As predicted, aftershocks may last for at least one year, and moreover, large aftershocks are likely to occur. Therefore, it is critical to outline the areas with potential aftershocks before reconstruction and re-settling people as to avoid future disasters. It is demonstrated that the redistribution of stress induced by an earthquake should trigger successive seismic activity. Based on static stress triggering theory, we calculated the coseismic stress changes on major faults induced by the Wenchuan earthquake, with elastic dislocation theory and the multilayered crustal model. We also discuss the stress distribution and its significance for future seismic activity under the impact of the Wenchuan earthquake. It is shown that coulomb failure stress (CFS) increases obviously on the Daofu-Kangding segment of the Xianshuihe Fault, the Maqu and Nanping segment of the Eastern Kunlun Fault, the Qingchuan Fault, southern segment of the Minjiang Fault, Pengxian-Guanxian Fault, Jiangyou-Guangyuan Fault, and Jiangyou-Guanxian Fault. The increased stress raises the probability of earthquake occurrence on these faults. Since these areas are highly populated, earthquake monitoring and early disaster alarm system are needed. CFS increases with a magnitude of 0.03–0.06 MPa on the Qingchuan Fault, which is close to the northern end of the rapture of Wenchuan earthquake. The occurrence of some strong aftershocks, including three events with magnitude higher than 5.0, indicates that the seismic activities have been triggered by the main shock. Aftershocks seem to migrate northwards. Since the CFS change on the Lueyang-Mianxian Fault located on the NEE of the Qingchuan Fault is rather small (±0.01 MPa), the migration of aftershocks might be terminated in the area near Hanzhong City. The CFS change on the western Qinling Fault is around 10 Pa, and the impact of static triggering can be neglected. The increment of CFS on the Pengxian-Guanxian Fault and Beichuan-Yingxiu Fault southwest to the main rupture is 0.005–0.015 MPa, which would facilitate earthquake triggering in these areas. Very few aftershocks in these areas indicate that the accumulated stress has not been released sufficiently. High seismic risk is predicated in these areas due to co-seismic CFS loading. The Wenchuan earthquake released the accumulated CFS on the Fubianhe Fault, the Huya Fault, the Ha’nan-Qingshanwan Fault, and the Diebu-Bailongjiang Fault. The decrement of CFS changes on the Longquanshan Fault east to Chengdu City is about 0.002 MPa. The seismic activity will be depressed by decrement of CFS on these faults. Supported by Knowledge Innovation Program of Chinese Academy of Sciences (Grant No. KZCX-SW-153), National Natural Science Foundation of China (Grant Nos. 40574011 and 40474028)