Coseismic Coulomb failure stress changes caused by the 2017M7.0 Jiuzhaigou earthquake,and its relationship with the 2008Wenchuan earthquake

On August 8, 2017, a M7.0 earthquake occurred in Jiuzhaigou County, Sichuan Province, China, resulting in significant casualties and property damage. Therefore, it is critical to identify the areas of potential aftershocks before reconstruction and re-settling people to avoid future disasters. Based on the elastic dislocation theory and a multi-layered lithospheric model, we calculate the Coulomb failure stress changes caused by the Wenchuan and Jiuzhaigou earthquakes, discuss the relationship between the Mw7.9 Wenchuan and M7.0 Jiuzhaigou earthquakes, and analyze the influence of the aftershock distribution and stress changes on the major faults in this region caused by the Jiuzhaigou earthquake. The co-and post-seismic stress changes caused by the Wenchuan earthquake significantly increased the stress accumulation at the hypocenter of the Jiuzhaigou earthquake. Therefore,the occurrence of the Jiuzhaigou earthquake was probably stimulated by the Wenchuan earthquake. The aftershock distribution is well explained by the co-seismic stress changes of the Jiuzhaigou earthquake. The stress accumulation and corresponding seismic hazard on the Maqu-Heye segment of the East Kunlun fault and the northern extremity of the Huya fault has been further increased by the Jiuzhaigou earthquake.     

Preliminary analysis on the source properties and seismogenic structure of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>7.0 Jiuzhaigou earthquake

At GMT time 13:19, August 8, 2017, an M_s7.0 earthquake struck the Jiuzhaigou region in Sichuan Province, China, causing severe damages and casualties. To investigate the source properties, seismogenic structures, and seismic hazards, we systematically analyzed the tectonic environment, crustal velocity structure in the source region, source parameters and rupture process, Coulomb failure stress changes, and 3-D features of the rupture plane of the Jiuzhaigou earthquake. Our results indicate the following: (1) The Jiuzhaigou earthquake occurred on an unmarked fault belonging to the transition zone of the east Kunlun fault system and is located northwest of the Huya fault. (2) Both the mainshock and aftershock rupture zones are located in a region where crustal seismic velocity changes dramatically. Southeast to the source region, shear wave velocity at the middle to lower crust is significantly low, but it rapidly increases northeastward and lies close to the background velocity across the rupture fault. (3) The aftershock zone is narrow and distributes along the northwest-southeast trend, and most aftershocks occur within a depth range of 5–20 km. (4) The focal mechanism of the Jiuzhaigou earthquake indicates a left-lateral strike-slip fault, with strike, dip, and rake angles of 152°, 74° and 8°, respectively. The hypocenter depth measures 20 km, whereas the centroid depth is about 6 km. The co-seismic rupture mainly concentrates at depths of 3–13 km, with a moment magnitude (M_w) of 6.5. (5) The co-seismic rupture also strengthens the Coulomb failure stress at the two ends of the rupture fault and the east segment of the Tazang fault. Aftershocks relocation results together with geological surveys indicate that the causative fault is a near vertical fault with notable spatial variations: dip angle varies within 66°–89° from northwest to southeast and the average dip angle measures ~84°. The results of this work are of fundamental importance for further studies on the source characteristics, tectonic environment, and seismic hazard evaluation of the Jiuzhaigou earthquake.     

Numerical simulation of dynamic Coulomb stress changes induced by <Emphasis Type="Italic">M</Emphasis>6.5 earthquake in Wuding,Yunnan and its relationship with aftershocks

Based on the discrete wavenumber method, we calculate the fields of dynamic Coulomb rupture stress changes and static stress changes caused by M6.5 earthquake in Wuding, and study their relationship with the subsequent aftershocks. The results show that the spatial distribution patterns of the positive region of dynamic stress peak value and static stress peak value are similarly asymmetric, which are basically identical with distribution features of aftershock. The dynamic stress peak value and the static stress in the positive region are more than 0.1 MPa and 0.01 MPa of the triggering threshold, respectively, which indicates that the dynamic and static stresses are helpful for the occurrence of aftershock. This suggests that both influences of dynamic and static stresses should be considered other than only either of them when studying aftershock triggering in near field.     

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)     

Deep tectonic setting of the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake in southwestern China

Teleseismic P-wave receiver functions at 20 broadband seismic stations in the Longmenshan fault zone (LMFZ) and its vicinity were extracted, and the crustal thickness and the P- and S-wave velocity ratio were calculated by use of the H-k stacking algorithm. With the results as constraints, the S-wave velocity structures beneath each station were determined by the inversion of receiver functions. The crustal structure of the Rear-range zone is similar to that of the Songpan-Garze Block, whereas the velocity structure of the Fore-range zone resembles that of Sichuan Basin, implying that the Central Principal Fault of LMFZ is the boundary between the eastern Tibetan Plateau and the Yangtze Block. Lower velocity zone exists in lower crust of the Songpan-Garze Block and the central-southern segment of the Rear-range zone, which facilitates the detachment of the material in upper and middle crust. Joint analysis of the receiver functions and the Bouguer gravity anomalies supports the thesis on the detachment-thrust mode of the LMFZ. A double-detachment pattern is suggested to the tectonic setting in the Songpan-Garze Block. The upper detachment occurs at the depth of 10-15 km, and represents a high-temperature ductile shear zone. There is a lower detachment at the depth of about 30 km, below which the lower crust flow exists in the eastern Tibetan Plateau. Interpretation of the Bouguer gravity anomalies indicates that the Sichuan Basin is of higher density in upper and middle crust in comparison with that of the Songpan-Garze Block. The LMFZ with higher density is the result from the thrusting of the Songpan-Garze Block over the Sichuan Basin. In the lower crust, higher P velocity and higher density in the Sichuan Basin are related to more rigid material, while lower S velocity and lower density in the Songpan-Garze Block are related to the softened and weakened material. The higher density block beneath the Sichuan Basin obstructs the eastward flow of lower crustal material from the Tibetan Plateau, which is driven by the compression of northward movement of Indian Plate. The eastward movement of upper and middle crustal material is also obstructed by the rigid Yangtze Block, resulting in the stress concentrated and accumulated along the LMFZ. When the stress releases sharply, the Wenchuan M _s8.0 earthquake occurs. Supported by the National Natural Science Foundation of China (Grant Nos. 40334041, 40774037) and Joint Foundation of Earthquake Science (Grant No. 1040062)     

2008年汶川地震同震-震后应力演化及其对2017年九寨沟M_S 7.0地震的影响

2017年8月8日四川九寨沟发生了MS7.0地震,距离汶川地震震中约为140km.地震的发生会调整区域应力分布和影响区域地震活动性变化.那么,汶川地震的发生对此次九寨沟地震有何影响?区域应力如何演化?为回答这些问题,本研究基于高性能并行有限元方法和网格自适应加密技术,采用含地形起伏的黏弹非均匀椭球型地球模型,分别计算了汶川地震引起的同震-震后应力时空动态演化及其对此次九寨沟地震的影响.同时,考虑同震静态应力和震后黏弹性应力调整,计算了汶川、芦山和九寨沟地震的发生对周边断层应力积累的影响.计算结果显示汶川地震对九寨沟发震断层的同震库仑应力加载约为0.008 MPa,震后9年的黏弹性应力加载约为0.012~0.016 MPa,可能会使九寨沟地震提前发生.此外,汶川、芦山和九寨沟地震的发生引起龙门山断裂映秀以南段、东昆仑断裂、龙日坝断裂东段和岷江断裂北段及鲜水河断裂康定段应力加载大于0.01 MPa,增加了这些断裂带发生地震的危险性.     

The seismotectonic and macroseismic features of the Wenchuan (Sichuan) earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript> = 8.0) of May 12, 2008

A disastrous earthquake with a magnitude M _S = 8.0 (M _W = 7.9), in China called “the 5.12 Wenchuan earthquake,” occurred on May 12, 2008, in Sichuan province on the border between the Sino-Tibetan Mountains and the Sichuan depression. The instrumental epicenter was registered in the southeastern part of Wenchuan county, and the hypocenter depth was 14 km. As the strongest and most destructive earthquake within mainland China, it caused numerous human losses and destruction of buildings and infrastructure. The seismic effect from the main shock and aftershocks was felt in many counties, towns, and villages, though Sichuan province suffered the most. The maximum intensity of the shocks was estimated at 11 degrees, according to the Chinese macroseismic scale. In the process of source opening, from the southern part of Wenchuan county to the vicinities of Quingchuan, a seismic fault system with a total length up to 240 km out-cropped on the earth’s surface, confined to the Longmenshan fault belt. The seismic fault system disturbed the original ground, resulting in the collapse or damage to various constructions, such as buildings, homes, bridges, roads, etc. Fault offsets had a dextral strike-slip and thrust kinematic combination. The earthquake generated several tens of thousands of landslides, rockfalls, and debris flows. Many dammed ponds appeared in the epicentral zone due to the activation of landslides. Thus, the geological effects turned out to be the most destructive factor in this case. At the same time, the seismic intensity of surface shaking was abnormally low even in direct proximity to the seismic fault system. Usually it was no more than 7–8 degrees. This macroseismic phenomenon may turn out to be rather typical for many major earthquakes.     

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.     

Numerical simulation of strong ground motion for the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground motion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the northwestern margin of the Sichuan Basin and caused by both the directivity of fault rupture and the amplification in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard assessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008. Supported by the U.S. National Science Foundation (Grant Nos. EAR 0738779 and OCE 0727919), the National Basic Research Program of China (Grant No. 2004CB418404), and partially by the National Nature Science Foundation of China (Grant No. 40521002)     

Study on a gravity precursor mode of the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>=7.0 Lijiang earthquake

To study the seismogeny process or the precursory behavior of the 1976 M _S=7.0 Lijiang earthquake, we analyze the repeat gravity data with high precision from the Western Yunnan Earthquake Prediction Experiment Area (WYEPEA) and the related results of geology and geophysics survey in this paper. Considering the gross errors caused by observation data and model difference, we have firstly inverted the slip distribution of the main active faults with time based on the robust Bayesian least squares estimation and multi-fault dislocation model. The results show that the slip changes of the faults with time from 1990 to 1997 obviously reflect the preparation process of the Lijiang earthquake. The images of main precursor mode have the characteristic of main shock-after shock type, which is agreement to the model of coupling movement between crust density and crust deformation (DD mode of coupling movement).     

A new technique for moment tensor inversion with applications to the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake sequence

A M_S8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (M_S5.0～6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was briefly introduced, and an aftershock of the Wenchuan earthquake sequence was employed to illustrate the real procedure for inverting the moment tensor; secondly, the moment tensor solutions and STFs of the eight events, including the main shock, were presented, and finally, the interpretation of the results was made. The agreement of our results with the GCMT results indicates the new approach is efficient and feasible. By using this approach, not only the moment tensor solution can be obtained but also the STF can be retrieved; the inverted STFs indicate that the source rupture process may be complicated even for the moderate earthquakes. The inverted focal mechanisms of the Wenchuan earthquake sequence show that the most of the aftershocks occurred in the main faults of the Longmenshan fault zone with predominantly thrustingwith minor right-lateral strike-slip component, but some of them may have occurred in the subfaults with strike-slip faulting in the vicinity of the main faults.     

Coseismic displacement estimate of the 2013 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>7.0 Lushan,China earthquake based on the simulation of near-fault displacement field

Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from InSAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the M _S7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were constructed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault, rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.     

Seismotectonic environment of occurring the February 3, 1996 Lijiang <Emphasis Type="Italic">M</Emphasis>=7.0 earthquake,Yunnan province

Lijiang-Daju fault, the seismogenic fault of the 1996 Lijiang M=7.0 earthquake, can be divided into Lijiang-Yuhu segment in the south and Yuhu-Daju segment in the north. The two segments show clear difference in geological tectonics, but have the similar dynamic features. Both normal dip-slip and sinistral strike-slip coexist on the fault plane. This kind of movement started at the beginning of the Quaternary (2.4–2.5 Ma B.P.). As to the tectonic types, the detachment fault with low angle was developed in the Early Pleistocene and the normal fault with high angle only after the Mid-Pleistocene (0.8 Ma B.P.). Based on the horizontal displacements of gullies and the vertical variance of planation surfaces cross the Lijiang-Daju fault at east piedmont of Yulong-Haba range, the average horizontal and vertical slip rates are calculated. They are 0.84 mm/a and 0.70 mm/a since the Quaternary and 1.56 mm/a and 1.69 mm/a since the Mid-Pleistocene. The movements of the nearly N-S-trending Lijiang-Daju fault are controlled not only by the regional stress field, but also by the variant movement between the Yulong-Haba range and Lijiang basin. The two kinds of dynamic processes form the characteristics of seismotectonic environment of occurring the 1996 Lijiang earthquake.     

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.     

Distribution and characteristics of gravelly soil liquefaction in the Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript> 8.0 earthquake

In this paper,a distribution map of gravelly soil liquefaction that was caused by the Wenchuan M_s 8.0 earthquake in China is proposed based on a detailed field investigation and an analysis of geological soil profiles. The geological background of the earthquake disaster region is summarized by compiling geological cross sections and borehole logs. Meanwhile,four typical liquefied sites were selected to conduct sample drillings,dynamic penetration tests (DPT),and shear wave velocity tests,to understand the features of liquefied gravelly soil. One hundred and eighteen (118) liquefied sites were investigated shortly after the earthquake. The field investigation showed:(1) sandboils and waterspouts occurred extensively,involving thousands of miles of farmland,120 villages,eight schools and five factories,which caused damage to some rural houses,schools,manufacturing facilities and wells,etc.; (2) the Chengdu plain is covered by a gravelly soil layer with a thickness of 0 m to 541 m according to the geological cross sections; (3) there were 80 gravelly soil liquefied sites in the Chengdu plain,shaped as five belt areas that varied from 20 km to 40 km in length,and about ten gravelly soil liquefied sites distributed within Mianyang area; and (4) the grain sizes of the sampled soil were relative larger than the ejected soil on the ground,thus the type of liquefied soil cannot be determined by the ejected soil. The gravelly soil liquefied sites are helpful in enriching the global database of gravelly soil liquefaction and developing a corresponding evaluation method in further research efforts.     

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>6.5 Jiuzhaigou,China, earthquake

The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 10¹⁸ N·m (M_w?≈?6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes >?6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 M_w7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.     

GPS time-series and its response to <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 Kunlunshan earthquake

In this paper, observation data in 25 GPS reference stations of China have been analyzed by calculating GPS position coordinate time-series with GIPSY. Result shows there is an obvious trend variation in such time-series. The trend variations of time series along the longitude and latitude coordinate reflect the motion of each position in the global-plate, in which the trend variation in the vertical direction reveals some large-scale construction information or reflects the local movement around the positions. The analysis also shows that such time-series have a variation cycle of nearly 1.02 a, but the reason still remains to be further studied. At the end of this paper, response of the time-series of M _S=8.1 Kunlunshan earthquake was analyzed, and the seismogenic process of M _S=8.1 Kunlunshan earthquake, according to the time proceeding and the feature of anomaly, was divided into 3 phases—changes in blocks with forces, strain accumulation, quick accumulation and slow release of energy. At the initial stage of seismogenic process of M _S=8.1 earthquake and at the imminent earthquake, coseismic process as well as during the post earthquake recovery, anomaly in vertical direction is always in a majority. The anomalous movement in vertical direction at the initial stage resulted in a blocking between faults, while at the middle stage of seismogenic process, the differential movement between blocks are in a majority, which is the major reason causing energy accumulating at the blocking stage of faults.     

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)     

Focal depths for moderate-sized aftershocks of the Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake and their implications

Sliding-window cross-correlation method is firstly adopted to identify sPn phase, and to constrain focal depth from regional seismograms, by measuring the time separation between sPn and Pn phases. We present the focal depths of the 17 moderate-sized aftershocks (M _S⩾5.0) of the Wenchuan M _S8.0 earthquake, using the data recorded by the regional seismic broadband networks of Shaanxi, Qinghai, Gansu, Yunnan and Sichuan. Our results show focal depths of aftershocks range from 8 to 20 km, and tend to cluster at two average depths, separate at 32.5°N, i.e., 11 km to the south and 17 km to the north, indicating that these aftershocks are origin of upper-to-middle crust. Combined with other results, we suggest that the Longmenshan fault is not a through-going crustal fault and the Pingwu-Qingchuan fault may be not the northward extension of the Longmenshan thrust fault. Supported by the National Natural Science Foundation of China (Grant Nos. 40604009 and 40574040) and Special Project for the Fundamental R & D of Institute of Geophysics, China Earthquake Administration (Grant No.DQJB08B20)     

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.