Rupture behavior and deformation localization of the Kunlunshan earthquake (M_w7.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 earth- quake (Mw7.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 transten- sional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the Mw=6.8, Mw=6.2 and Mw≤7.8 events, respec- tively. The Mw≤7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earth- quake, 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 de- pends 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 rup- tures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Ba- yan 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.     

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)     

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 Boumerdes (Algeria) earthquake of May 21, 2003 (<Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> = 6.8): Ground deformation and intensity

A destructive earthquake of magnitude M_w = 6.8 hit the region of Boumerdes and Algiers (Algeria) on May 21, 2003. This is among the strongest seismic events of the mediterranean region and the most important event in the capital Algiers since 1716. It caused a widespread damage in the epicentral region, claimed 2271 human lives, injured 10000, about 20000 housing units affected and left about 160000 homeless. The main shock was felt about 250 km far from the epicenter and triggered sea waves of 1–3 m in amplitude in Balearic islands (Spain). Based on field observations and press report an intensity IX (MSK scale) is attributed to the epicentral area. The main shock was followed by many aftershocks among them several are of magnitude greater than 5.0, which added panic to inhabitants. The main shock triggered ground deformation, particularly liquefaction whose features are in different forms and sizes and caused damage and collapse of roads. The focal mechanism determined by worldwide institutions yield a pure reverse faulting with a compressional axis striking NE-SW. The epicenter is located offshore about 7 km from the Boumerdes-Dellys coast. Field observations show 0.7 m of coseismic uplift of shoreline between Boudouaou and Dellys. This uplift is about a half of the extracted coseismic slip from the seismic moment. On the other hand there is no clear surface break onshore, confirming hence, that the causative active fault is offshore. However, the rupture may propagate onshore to the SE near the Boudouaou region where ground cracks showing reverse faulting are observed a long a corridor of about 1 km wide. These fissures may correspond to a diffuse coseismic deformation.     

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)     

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)     

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.     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

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.     

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.     

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.     

Frequency-dependent rupture process,stress change,and seismogenic mechanism of the 25 April 2015 Nepal Gorkha <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.8 earthquake

On 25 April 2015, an M _w 7.8 earthquake occurred on the Main Himalaya Thrust fault with a dip angle of ~ 7° about 77 km northwest of Kathmandu, Nepal. This Nepal Gorkha event is the largest one on the Himalayan thrust belt since 1950. Here we use the compressive sensing method in the frequency domain to track the seismic radiation and rupture process of this event using teleseismic P waves recorded by array stations in North America. We also compute the distribution of static shear stress changes on the fault plane from a coseismic slip model. Our results indicate a dominant east-southeastward unilateral rupture process from the epicenter with an average rupture speed of ~3 km s^?1. Coseismic radiation of this earthquake shows clear frequency-dependent features. The lower frequency (0.05–0.3 Hz) radiation mainly originates from large coseismic slip regions with negative coseismic shear stress changes. In comparison, higher frequency (0.3–0.6 Hz) radiation appears to be from the down-dip part around the margin of large slip areas, which has been loaded and presents positive coseismic shear stress changes. We propose an asperity model to interpret this Nepal earthquake sequence and compare the frequency-dependent coseismic radiation with that in subduction zones. Such frequency-dependent radiation indicates the depth-varying frictional properties on the plate interface of the Nepal section in the main Himalaya thrust system, similar to previous findings in oceanic subduction zones. Our findings provide further evidence of the spatial correlation between changes of static stress status on the fault plane and the observed frequency-dependent coseismic radiation during large earthquakes. Our results show that the frequency-dependent coseismic radiation is not only found for megathrust earthquakes in the oceanic subduction environment, but also holds true for thrust events in the continental collision zone.     

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.     

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.     

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.     

Preliminary analysis on characteristics of coseismic deformation associated with <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 western Kunlunshan Pass earthquake in 2001

Based on the analysis of coseismic deformation in the macroscopic epicentral region extracted by Differential Interferometric Synthetic Aperture Radar (D-InSAR), and combined with the seismic activity, focal mechanism solutions of the earthquake and field investigation, the characteristic of coseismic deformation of M _S=8.1 western Kunlunshan Pass earthquake in 2001 was researched. The study shows that its epicenter lies in the northeast side of Hoh Sai Hu; and the seismogenic fault in the macroscopic epicentral region can be divided into two central deformation fields: the west and east segments with the lengths of 42 km and 48 km, respectively. The whole fault extends about 90 km. From the distribution of interferometry fringes, the characteristic of sinistral strike slip of seismogenic fault can be identified clearly. The deformations on both sides of the fault are different with an obviously higher value on the south side. In the vicinity of macroscopic epicenter, the maximum displacement in look direction is about 288.4 cm and the minimum is 224.0 cm; the maximum sinistral horizontal dislocation of seismogenic fault near the macroscopic epicenter is 738.1 cm and the minimum is 551.8 cm.     

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.     

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.     

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.