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.     

3D v_P and v_S models of southeastern margin of the Tibetan plateau from joint inversion of body-wave arrival times and surface-wave dispersion data

A new 3D velocity model of the crust and upper mantle in the southeastern(SE) margin of the Tibetan plateau was obtained by joint inversion of body-and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the twolocated to the northeast of 2013 M S7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes(M≥5) in the region occurring from2008 to 2015 are located around the high-velocity zones,indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.     

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.     

Multi-scale Tomography for Crustal <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Velocities in Southern California

Seismic tomography is a viable tool in building depth-velocity models in the presence of strong lateral velocity variations. In this study 3-D P- and S-velocity models for the crust of southern California are constrained using more than 1,000,000 P-wave first arrivals and 130,000 S-wave arrivals from local earthquakes. To cope with the uneven distribution of raypaths, a multi-scale tomography is applied with overlapping model cells of different sizes. Within the 300 × 480 × 39 km³ model volume, the smallest cell size is 10 × 10 × 3 km³. During the iterations of velocity updating, earthquake hypocenters are determined using both P and S arrivals, and full 3-D ray tracing is implemented. Except near the edges and in the lower crust, the resultant models are robust according to various tests on the effects of reference models, resolution and signal-to-noise ratio. The tomographic velocities at shallow depths correlate very well with the regional geology of southern California. In the upper crust the P-wave and S-wave models exhibit slow velocities in major sedimentary basins and fast velocities in areas of crystalline rocks. Mid-crustal low velocity zones are present under the Coso Range, San Gabriel Mountains, and a large portion of the Mojave Desert. P- and S-velocity patterns maintain their similarity in the lower crust though the models are less reliable there.     

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

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

On the <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Receiver Functions Used for Inverting the One-Dimensional Upper Mantle Shear-Wave Velocities

The teleseismic P receiver functions are customarily inverted to attain the seismic velocities beneath a seismic station. Surface wave dispersion data are often added to reduce the effect of the non-uniqueness. The combination of P receiver function and surface wave works well in resolving the structures in the crust and uppermost mantle, but is less effective in characterizing greater (lithosphere and asthenosphere) depths due to the interference from crustal multiples. A solution to this problem is jointly to model teleseismic S receiver functions with surface wave and P receiver functions. This study adopts a fast, one-dimensional (1-D) inversion scheme. To avoid the effect of multidimensional structures away from the seismic station, we eliminate multiples that reverberate between the surface and interfaces below a restriction depth (RD), as well as S-to-P conversions below an inversion depth (ID). P-to-S conversions off the interfaces above the half-space and S-to-P conversions above the ID and multiples above the RD are properly modelled. This approach favours ray paths travelling close to stations and is, therefore, more suitable for 1-D inversions. We perform numerical experiments with and without noise and highlight the advantages of a joint receiver function and surface wave analysis.     

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

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

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

Linear transformations of anomalous potential fields using <Emphasis Type="Italic">S</Emphasis> and <Emphasis Type="Italic">R</Emphasis> approximations

Methods for linear transformations of anomalous physical fields based on R and S approximations of the elements of the initial field are described. The results of the numerical experiment and analytical extension of the gravity field measured in a detailed gravity survey in two regions of Russia are presented.     

Attenuative Dispersion of <Emphasis Type="Italic">P</Emphasis> Waves and Crustal <Emphasis Type="Italic">Q</Emphasis> in Turkey and Germany

We have measured group delays of the spectral components of high-frequency P-waves along two portions of the North Anatolian Fault Zone (NAFZ) in Turkey and in a region of southern Germany. Assuming that the observed dispersion is associated with attenuation in the crust and that it can be described by a continuous relaxation model, we obtained Q and the high-frequency relaxation times for those waves for each of the three regions. Individual P-wave Q values exhibit large scatter, but mean values in the NAFZ increase from about 25 to 60 over the distance range 5–90 km. Mean Q values are somewhat higher in the eastern portion of the NAFZ than in the western portion for measurements made at distances between 10 and 30 km. P-wave Q values in Germany range between about 50 and 300 over the hypocentral distance range 20–130 km. In that region we separated the effects of Q for basement rock (2–10 km depth) from that of the overlying sediment (0–2 km depth) using a least-squares method. Q varies between 100 and 500 in the upper 8–10 km of basement, with mean values for most of the distance range being about 250. Q in the overlying sediments ranges between 6 and 10. Because of large scatter in the Q determinations we investigated possible effects that variations of the source-time function of the earthquakes and truncation of the waveform may have on Q determinations. All of our studies indicate that measurement errors are relatively large and suggest that useful application of the method requires many observations, and that the method will be most useful in regions where the number of oscillations following the initial P pulse is minimized. Even though there is large scatter in our Q determinations, the mean values that we obtained in Turkey are consistent with those found in earlier studies. Our conclusions that Q is significantly higher in the basement rock of Germany than in the basement rock of Turkey and that Q is lower in western Turkey than in eastern Turkey are also consistent with results of Q studies using Lg coda.     

Experimental examination on the heterogeneity parameter <Emphasis Type="Italic">C</Emphasis><Subscript>v</Subscript> of earthquake precursors

Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C _v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C _v value and to explore its physical meaning experimentally. The result shows that the change of C _v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C _v value increases firstly and then decreases before the instability, and the instability occurs when C _v value decreases to the level before increasing. This indicates that C _v value may be a useful parameter for earthquake prediction.     

Stable carbon isotope evidence for tracing the diet of the host <Emphasis Type="Italic">Hepialus</Emphasis> larva of <Emphasis Type="Italic">Cordyceps sinensis</Emphasis> in the Tibetan Plateau

Two types of Hepialus larvae with different diets were distinguished in the Sejila Mountain, Tibetan Plateau based on the stable carbon isotope data of the host Hepialus larva of Cordyceps sinensis and its closely adjacent tender plant roots and humus fractions. Type I is the larva chiefly fed by soil humus, and characterized by the δ ¹³C values of −22.6‰ to −23.4‰, and more than −23.4‰ in its heads. Type II is the larva chiefly fed by tender plant roots, and characterized by the δ ¹³C values of −24.6‰ to −27.6‰, and less than −24.6‰ in its heads. Our result has exceeded the traditional understanding that their food sources only come from the tender plant roots, and may provide evidence for choosing cheap and high-quality foods and further establishing artificial habitats in their large-scale reproduction. Supported by the National Key Technology R & D Program (Grant No. 2007BAI32B05)     

利用接收函数、背景噪声频散和Rayleigh波ZH振幅比联合反演青藏高原东北缘及周边地区地壳结构

本文对在青藏高原东北缘及邻近地区架设的571套宽频带流动地震台站记录到的连续波形数据进行处理,通过联合反演P波接收函数、背景噪声频散和Rayleigh波ZH振幅比,获得了研究区高分辨率三维地壳S波速度（V_S）结构.研究结果显示,松潘—甘孜块体和西秦岭造山带以及北祁连造山带下方15~40 km的深度范围内存在很显著的S波低速异常体.其中,松潘—甘孜块体和西秦岭造山带下方的低速异常体很可能与部分熔融有关,且造成部分熔融的原因除了剪切加热外,还有可能是软流圈物质上涌和地壳内部的放射生热.而北祁连造山带的低速异常体则可能由地表隆升和地壳增厚所造成.阿拉善块体内部分布有很多不太显著的低速异常体,这可能与阿拉善块体经历了复杂的构造变形有关.在银川—河套地堑下方20~35 km的深度范围内同样观测到了相对不太显著的低速异常体,这更可能与基性岩浆的底侵作用有关.

     

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.     

Combined Method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-Approximations of Increased Dimensionality in Solving the Problems of Geophysics and Geomorphology

The interrelation between different variants of the method of linear integral representations in the spaces of an arbitrary dimension is considered. The combined approximations of the topography and geopotential fields allows the selection of the optimal parameters of the method in solving a wide range of inverse problems in geophysics and geomorphology, as well as a most thorough use of the a priori information about the elevations and elements of anomalous fields. A method for numerically solving an inverse problem on finding the equivalent, in terms of the external field, mass distributions in the ordinary three-dimensional (3D) space and in the four-dimensional (4D) space is described.     

Combined method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-approximations in solving the problems of geophysics and geomorphology

The interrelation between different modifications of the method of linear integral representation is studied. Combined approximations of the topography and geopotential fields enable more refined tuning of the method in solving inverse problems of geophysics and geomorphology and provide a more complete allowance for the a priori information about the surface elevation data and elements of anomalous fields. A technique for finding the numerical solution for the inverse problem for determining the mass distributions equivalent in terms of the external field is presented. The results of the mathematical experiment are discussed.     

Gravity variation before Kunlun mountain pass western <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 earthquake

The relation between the gravity variation features and M _S=8.1 earthquake in Qinghai-Xizang monitoring area is analyzed preliminarily, by using spatial dynamic variation results of regional gravity field from absolute gravity and relative gravity observation in 1998 and 2000. The results show that: 1) M _S=8.1 earthquake in Kulun mountain pass western occurred in the gravity variation high gradient near gravity’s high negative variation; 2) The main tectonic deformation and energy accumulation before M _S=8.1 earthquake are distributed at south side of the epicenter; 3) The range of gravity’s high negative variation at east of the M _S=8.1 earthquake epicenter relatively coincides with that rupture region according to field geology investigation; 4) Gravity variation distribution in high negative value region is just consistent with the second shear strain’s high value region of strain field obtained from GPS observation.