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.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

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.     

Dayside auroras during unusually large positive values of the IMF <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> component

The characteristics of dayside auroras during the large (16–24 nT) positive values of the IMF B _z component, observed on January 14, 1988, during the interaction between the Earth’s magnetosphere and the body of the interplanetary magnetic cloud, have been studied based on the optical observations on Heiss Island. A wide band of diffuse red luminosity with an intensity of 1–2 kilorayleigh (kR) was observed during 6 h in the interval 1030–1630 MLT at latitudes higher than 75° CGL. Rayed auroral arcs, the brightness of which in the 557.7 nm emission sharply increased to 3–7 kR in the postnoon sector immediately after the polarity reversal of the IMF B _y component from positive to negative, were continuously registered within the band. Bright auroral arcs were observed at the equatorward edge of red luminosity. It has been found out that the red auroral intensity increases and the band equatorward boundary shifts to lower latitudes with increasing solar wind dynamic pressure. However, a direct proportional dependence of the variations in the auroral features on the dynamic pressure variations has not been found. It has been concluded that the source of bright discrete auroras is located in the region of the low-latitude boundary layer (LLBL) on closed geomagnetic field lines. The estimated LLBL thickness is ∼3 R _e . It has been concluded that the intensity of the dayside red band depends on the solar wind plasma density, whereas the position of the position equatorward boundary depends on the dynamic pressure value and its variations.     

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.     

Array performances for ambient vibrations on a shallow structure and consequences over <Emphasis Type="Italic">V</Emphasis><Subscript><Emphasis Type="Italic">s</Emphasis></Subscript> inversion

Valuable information about one-dimensional soil structures can be obtained by recording ambient vibrations at the surface, in which the energy contribution of surface waves predominates over the one of other types of waves. The dispersion characteristics of surface waves allow the retrieval of the shear-wave velocity as a function of depth. Microtremor studies are usually divided in two stages: deriving the dispersion (or auto-correlation) curve from the recorded signals and inverting it to obtain the site velocity profile. The possibility to determine the dispersion curve over the adequate frequency range at one site depends on the array aperture and on the wavefield spectra amplitude that can be altered by filtering effects due to the ground structure. Microtremors are usually recorded with several arrays of various apertures to get the spectral curves over a wide frequency band, and different methods also exist for processing the raw signals. With the objective of defining a strategy to achieve reliable results for microtremor on a shallow structure, we analyse synthetic ambient vibrations (vertical component) simulated with 333 broadband sources for a 25-m deep soil layer overlying a bedrock. The first part of our study is focused on the determination of the reliable frequency range of the spectral curves (dispersion or auto-correlation) for a given array geometry. We find that the wavenumber limits deduced from the theoretical array re sponse are good estimates of the valid spectral curve range. In the second part, the spectral curves are calculated with the three most popular noise-processing techniques (frequency–wavenumber, high-resolution frequency–wavenumber and spa tial auto-correlation methods) and inverted indi vidually in each case. The inversions are performed with a tool based on the neighbour hood algorithm that offers a better estimation of the global uncertainties than classical linearised methods, especially if the solution is not unique. Several array apertures are necessary to construct the dispersion (auto-correlation) curves in the appropriate frequency range. Considering the final velocity profiles, the three tested methods are almost equivalent, and no significant advantage was found for one particular method. With the chosen model, all methods exhibit a penetration limited to the bedrock depth, as a consequence of the filtering effect of the ground structure on the vertical component, which was observed in numerous shallow sites.     

Super-resolution least-squares prestack Kirchhoff depth migration using the <Emphasis Type="Italic">L</Emphasis><Subscript>0</Subscript>-norm

Least-squares migration (LSM) is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L₂-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L₀-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.     

Dependence of geomagnetic disturbances on extreme values of the solar wind <Emphasis Type="Italic">E</Emphasis><Subscript><Emphasis Type="Italic">y</Emphasis></Subscript> component

The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E _y component, dynamic pressure P _d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E _y value at E _y < 12 mV m^?1 and are inversely proportional to this value at E _y > 12 mV m^?1 for the first-type events. On the contrary, the dependence of Dst on E _y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E _y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

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.     

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.     

<Emphasis Type="Italic">R/S</Emphasis> analysis of the <Emphasis Type="Italic">Dst</Emphasis> index

The paper addresses estimation of the Hurst exponent for time series of the hourly values of the Dst index for the period from 1957 to 2011. It is found that the Hurst exponent is 0.79–0.94 for yearly intervals and 0.8–1.0 for monthly intervals. Based on R/S graphs, the Dst cycles are identified; they range from 3–4 months to 2.2 years and from 8.5 to 22 years in length. It is shown that a Dst time series can be quite satisfactorily described by an α-stable Levy process.     

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 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.     

Behavior of <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2 after sunset

The time behavior of the foF2 and hmF2 values at the time moment T(ss + 2 h) 2 h after sunset is considered. It is assumed that at this moment, the horizontal winds in the thermosphere in the strongest way influence hmF2 and, therefore, foF2. It is found that a fairly well pronounced and statistically significant change (trend) is observed for the foF2(ss + 2)/foF2(14) ratio, the sign of the change being different for different stations and even different seasons at the same station. A similar picture is obtained for the value of hmF2(ss + 2). It is shown that a positive correlation between the trends of these two values is observed. This confirms the initial concept of the paper that the foF2 and hmF2 trends are caused by long-term trends in the thermospheric dynamics.     

Distribution of<Emphasis Type="Italic">L</Emphasis><Subscript>g</Subscript> coda<Emphasis Type="Italic">Q</Emphasis> in the Chinese continent and its adjacent regions

785 traces of vertical components from shallow earthquakes recorded by 10 CDSN (Chinese Digital Seismographic Network) stations and 5 GSN (Global Seismographic Network) stations were collected to study the attenuation characteristics ofL _g coda in the Chinese continent and its adjacent regions. The records were processed first using the stack spectral ratio method to obtain the average values ofQ ₀ (Q at 1Hz) and η, the frequency dependence, ofL _g coda in the ellipses corresponding to the paths. The back-projection technique was then employed to obtain the tomographic maps ofQ ₀ and η values, and the distribution of their errors. Results indicate that in the studied areaQ ₀ varies between 200 and 500. The lowest value ofQ ₀ exists in the Yun-nan-Tibetan region, while the highest value ofQ ₀ occurs in the southern edge of Siberian platform. η varies between 0.3 and 0.8. For most part of the studied area η varies inversely withQ ₀.     

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.     

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.     

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)