Regional variation of Rayleigh wave attenuation coefficients in the Eastern Pacific

Seismograms recorded for five earthquakes on the east Pacific rise have been analyzed to obtain the attenuation coefficients of the fundamental Rayleigh mode for the eastern Pacific in the 15–110 second period range. The attenuation coefficients have been obtained using two new methods, a reference-station method, and an iterative method by which the seismic moment and regionalized attenuation coefficient values are obtained simultaneously after considering the effect of the source directivity and time-function. The reference-station method was applied to the entire eastern Pacific, excluding paths along the east Pacific rise. When using the iterative method we divided the eastern Pacific into three sub-regions, designated as the north-eastern Pacific, the Nazca plate and the east Pacific rise. Although much scatter is present, the data suggest that attenuation coefficients for the Nazca plate are higher than those for the northeastern Pacific, and both are substantially higher than average values obtained for the entire Pacific Ocean. Two paths that lie along or near the east Pacific rise are characterized by especially high attenuation coefficients. These values suggest that a low-Q zone exists beneath that narrow feature.     

S波影区地震精确定位的研究

长期的地震观测发现,北京地震台记录到的S波影区地震存在定位误差大的现象.由于来自震中距9°～36°的地震经过了低速层,造成了S波不发育.所以S波影区地震多用LG波来确定震中距.造成影区地震定位偏差大的现象有两个原因:①长期以来多是用短周期记录来区分LG1和LG2,定位偏差较大;②有的影区地震,S波也少许发育,但用S波来定位也存在定位偏差大的现象.为此搜集了北京地震台记录到的影区地震,通过NCDSN仿真成SK记录,通过理论震中距的比较,找到了准确定位影区地震的LG2识别位置.     

Using Long-Period Body Wave to Inverse Moment Tensors of the M_S6.8 Jiashi,Xinjiang Earthquake in 2003 and the Moderate and Small Earthquakes before and after It

An Ms6.8 strong earthquake took place in Jiashi, Xinjiang on February 24 of 2003. The digital wave form data recorded in Kashi and Wushi stations are selected to inverse the moment tensor solutions for the strong earthquake and the moderate and small earthquakes before and after it （ 108 earthquakes in 2001 - 2004）. 67 focal mechanism solutions have been calculated, and the results agree with those from Harvard University and USGS. The analysis reveals that before the strong earthquake, the moderate and small earthquake distribution was dispersed, and after the event the distribution was mainly concentrated around the strong earthquake. Before the strong earthquake, the seismic faults of the mid and small events had the character of strike-slip and normal faulting, and after the event, they exhibit strike-slip and thrust faulting. The region is dominated by near-NS horizontal compression from the southern block after the strong earthquake.     

Probability of radial anisotropy in the deep mantle

It is well established that the Earth's uppermost mantle is anisotropic, but observations of anisotropy in the deeper mantle have been more ambiguous. Radial anisotropy, the discrepancy between Love and Rayleigh waves, was included in the top 220 km of PREM, but there is no consensus whether anisotropy is present below that depth. Fundamental mode surface waves, for commonly used periods up to 200 s, are sensitive to structure in the first few hundred kilometers and therefore do not provide information on anisotropy below. Higher mode surface waves, however, have sensitivities that extend to and below the transition zone and should thus give insight about anisotropy at greater depths, but they are very difficult to measure. We previously developed a new technique to measure higher mode surface wave phase velocities with consistent uncertainties. These data are used here to construct probability density functions of a radially anisotropic Earth model down to approximately 1500 km. In the uppermost mantle, we obtain a high probability of faster horizontally polarized shear wave speed, likely to be related to plate motion. In the asthenosphere and transition zone, however, we find a high probability of faster vertically polarized shear wave speed. To a depth of 1500 km in the lower mantle, we see no significant shear wave anisotropy. This is consistent with results from laboratory measurements which show that lower mantle minerals are anisotropic but LPO is unlikely to develop in the pressure–temperature conditions present in the mid-mantle.     

The scattering of shear-waves in the crust

The two major sources of scattering for shear-waves in the crust, interactions with the topography at the surface and the effective anisotropy of aligned cracks throughout the rockmass, introduce first-order changes to the shear-wave particle-motion. At the surface, shear-waves are scattered by the topography within a wavelength or two of the recording site so that, unless the effective incidence angle is less than the critical angle sin^–1 V _S/V _P, the recorded waveforms may bear little relationship to the waveforms of the incident wave. Within the rockmass, shear-waves are scattered by extensive-dilatancy anisotropy (EDA), the distribution of stress-aligned fluid-filled cracks, microcracks, and preferentially oriented pore-space pervading most rocks in the crust. Analysis of this shear-wave splitting yields new information about the internal structure of thein situ rockmass which is not otherwise available.     

近30 a云贵高原湖泊表面水体面积变化遥感监测与时空分析

以云贵高原湖泊近30 a来的TM、ETM~+和OLI遥感影像为数据源,采用归一化水体指数(NDWI)、改进归一化水体指数(MNDWI)、新型水体指数(NWI)、增强型水体指数(EWI)和自动水体提取指数5种水体指数提取了1985—2015年云贵高原10个湖泊表面水体面积,并对各种算法进行精度对比分析.针对湖泊各自特点采用不同的水体指数提取其表面水体面积,并进行水体面积变化时空分析.结果表明:云贵高原湖泊表面水体面积总体呈现先增加后缩减趋势,1985—1995年湖泊表面水体面积增加了30.86 km~2,1995—2015年湖泊水体表面面积减少了48.12 km~2,其中,面积变化最大的是杞麓湖与异龙湖.对云贵高原湖泊表面水体面积变化与该区域的年降水量、蒸发量、平均气温、流域植被覆盖面积和人类活动时空进行相关分析,结果表明:1)高原湖泊对区域气候变化的响应具有明显的空间差异性,云贵高原湖泊的表面水体面积与气候相关性较显著,气温升高引起蒸发加速,降水量下降,湖面不断缩小,与逐年上升的气温呈负相关,与逐年波动上升的蒸发量呈负相关,与逐年减少的降水量呈正相关;2)云贵高原湖泊各流域的植被覆盖面积与湖泊面积变化相关性较弱;3)人类活动是影响湖泊面积变化的重要因素,大肆围湖造田、围湖养殖以及旅游开发等人类活动直接导致云贵高原湖泊面积的锐减,并对湖泊生态环境产生重要影响.     

Scattered Surface Wave Energy in the Seismic Coda

One of the many important contributions that Aki has made to seismology pertains to the origin of coda waves (Aki, 1969; Aki and Chouet, 1975). In this paper, I revisit Aki's original idea of the role of scattered surface waves in the seismic coda. Based on the radiative transfer theory, I developed a new set of scattered wave energy equations by including scattered surface waves and body wave to surface wave scattering conversions. The work is an extended study of Zeng et al. (1991), Zeng (1993) and Sato (1994a) on multiple isotropic-scattering, and may shed new insight into the seismic coda wave interpretation. The scattering equations are solved numerically by first discretizing the model at regular grids and then solving the linear integral equations iteratively. The results show that scattered wave energy can be well approximated by body-wave to body wave scattering at earlier arrival times and short distances. At long distances from the source, scattered surface waves dominate scattered body waves at surface stations. Since surface waves are 2-D propagating waves, their scattered energies should in theory follow a common decay curve. The observed common decay trends on seismic coda of local earthquake recordings particular at long lapse times suggest that perhaps later seismic codas are dominated by scattered surface waves. When efficient body wave to surface wave conversion mechanisms are present in the shallow crustal layers, such as soft sediment layers, the scattered surface waves dominate the seismic coda at even early arrival times for shallow sources and at later arrival times for deeper events.     

Attenuation of P and pP waves in Vrancea area – Romania

The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The amplitude spectral ratio method is used, assuming no frequency dependence of the Q _P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor is observed. Lateral variations of the attenuation are obtained, with a very low Q _P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different azimuth angles in relation to the epicentral area, Q _P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located immediately beneath the crust, but limited in depth to at least 100 km.     

Properties of the lithosphere-asthenosphere system in Europe with a view toward earth conductivity

Over the past 20 years the study ofP- andS-wave velocities in the upper mantle of the Mediterranean area and continental Europe has been the subject of intensive research work. We present a summary of results based on the inversion of available surface-wave dispersion data andP-wave trave time observations. For areas characterized by different tectonic settings and very large lateral variations, a discussion is made about structural models based on seismological, geothermal and electrical conductivity data.     

Comparison of techniques that use the single scattering model to compute the quality factorQ from coda waves

Using simulated coda waves, the resolution of the single-scattering model to extract codaQ (Q _c ) and its power law frequency dependence was tested. The back-scattering model ofAki andChouet (1975) and the single isotropic-scattering model ofSato (1977) were examined. The results indicate that: (1) The inputQ _c models are reasonably well approximated by the two methods; (2) almost equalQ _c values are recovered when the techniques sample the same coda windows; (3) lowQ _c models are well estimated in the frequency domain from the early and late part of the coda; and (4) models with highQ _c values are more accurately extracted from late code measurements.     

The lithosphere in the central-eastern Mediterranean area

The lithosphere beneath the central-eastern Mediterranean area has been investigated by the inversion of the regional dispersion relations derived from analysis of surface waves. It is possible to distinguish several types of crust with average S-wave velocities in the range 3.0–3.8 km/sec, and thicknesses varying from a minimum of about 30 km, which corresponds to the Apennines, Crete and Otranto Channel regions, to a maximum of about 51 km beneath the Ionian Sea, which can be considered as a submerged continent. Associated with these crustal features, large lateral variations have been detected in the lithosphere thickness, which varies from a minimum of about 30 km corresponding to the Tyrrhenian Sea and south of Crete to a maximum of about 130 km corresponding to south-eastern Alps and north-central Greece, while the sub-Moho S-wave velocity varies in the range 4.2–4.8 km/sec. The constraint furnished by our results to the geological-tectonic setting of the investigated area, characterized by the continent continent collision between Africa and Europe, is pointed out.Publication No. 405, P. F. Geodinamica, CNR, Roma, Italy.     

The mechanism of small mining tremors from amplitude inversion

A method of inverting seismic wave amplitudes from mining tremors is proposed to reveal their mechanism. A simple shear-implosive model of a mechanism is proposed. The shapes of theoretical and observed radiation patterns are compared in the least squares norm and the parameters describing the shear-slip together with the ratio of implosive and shear dislocations are optimized. The efficiency of the method, even in the case of a small number of stations, is demonstrated by applying it to weak mining tremors from two mine regions in Czechoslovakia.The weak mining tremors are found to be of a predominantly shear mechanism with a clear coincidence of nodal planes to local tectonics, which could imply their tectonic nature.     

利用S波高频衰减对天山中东段地区地壳Q值成像

利用天山中东段地区25个数字地震台站记录到的2009~2014年底5076个近震波形资料以及新疆地震台网提供的观测报告,采用遗传算法对S波位移谱的高频衰减进行拟合,得到19140条衰减算子t~*数据。根据这些t*数据,反演得到该区地壳Q_s值分布。结果表明,天山中东段地区平均Q_0值为520,其Q_s值分布及其所揭示的衰减变化特征与研究区的地表构造明显相关。天山南北两侧的山区盆地交汇部位Q_s值较低,而高Q_s值较为集中地分布于天山造山带内部。1900年以来天山中东段地区6级以上强震大多位于低Q_s值区域,该区24个高热流点也大多位于上述天山南北两侧的低Q_s值区域,热流值与衰减值呈负相关;此外,研究区的速度结构与衰减结构也呈一定的正相关,反映了二维衰减结构特征与速度结构、二维密度结构的一致性。     

Basin and crustal velocity structure models for the simulation of strong ground motions in the Kinki area,Japan

We constructed a prototype of the basin and crustal structure model for the Kinki area, southwest of Japan, for the simulation of strong ground motions of hypothetical crustal and subduction earthquakes. We collected results of the deep seismic velocity profiles obtained by the reflection experiments and seismic imaging results, which were conducted in the Kinki area. The obtained profiles give underground velocity structures of the crust, from the surface to the subducting slab. We also gather the basin velocity structure information of the Osaka, Kyoto, Nara, and Ohmi basins. To examine the applicability of the constructed velocity structure model to the ground motion simulation, we simulated waveforms of an intermediate size event occurred near the source area of the hypothetical subduction earthquakes. Simulated ground motions using the basin and crustal velocity structure model are fairly well reproducing the observations at most of stations, and the constructed basin and crustal velocity structure model is applicable for the long-period ground motion simulations.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

Reflections on measurement of residual stress in rock

Consideration of the behaviour of elastic bodies shows that it is not possible for a set of overcoring measurements that are made within isolated blocks to show residual strains or stresses that have a non-zero average, unless the size of the equilibrium volume over which the residual stresses balance is both considerably larger than the volume of the overcore and smaller than the size of the isolated block. Since some results have been reported that do not match these constraints, non-elastic behaviour must have occurred during overcoring. A possible explanation, is that oriented microcracks are somehow opened by the overcoring. In some cases, stress fields induced by microcracking near the overcoring cut may explain measured strain changes far from the cut. Consideration of various reported measurements in terms of relative sizes of overcoring and equilibrium volumes shows that care is necessary when interpreting residual stress measurements.     

Use of mode-converted waves in marine seismic data to investigate the lithology of the sub-bottom sediments in Isfjorden,Svalbard

The compression wavefield is efficiently converted to shear-wave energy at post-critical angles in areas of high impedance contrast at the sea floor. We have analysed mode-converted shear waves in a data set acquired with a hybrid marine/land geometry in Isfjorden, Svalbard. Through a kinematic 2D ray-tracing modellingV _p/V_s ratios for part of the uppermost 5km of the crust are obtained. Low values (V _p /V _s =1.65) are tentatively associated with the section of Devonian sandstones which appears to attain a minimum thickness of 1.5km below 3 km depth about 10km west of Kapp Thorden.     

Coherence of dispersed synthetic strong earthquake ground motion at small separation distances

It is shown that the observed loss of coherency of synthetic strong motion on ground surface for separation distances less than about 100 m can be described in terms of the dispersion of strong motion waves. Additional contributions to the loss of coherency from variations of material properties in the soil and from geometrical departures from perfectly flat ground surface and irregular layer geometries are not considered in this paper. It is also shown that the synthetic surface displacements over a large rectangular area (100×100 m) on ground surface can be described well by a flat surface undergoing translations and rotations only, and with only minor departures from the plane flat surface.     

内蒙古河套灌区3种沉水植物种子休眠及萌发特性

沉水植物在湖泊生态系统中具有重要作用,它不仅是湖泊生态系统食物链中重要的生产者,同时还对湖泊的营养物沉积和循环有着重要的意义,并由此影响着湖泊的富营养化进程.本文以黄河内蒙古段河套灌区湖泊中常见的3种沉水植物篦齿眼子菜(Potamogeton pectinatus)、穿叶眼子菜(Potamogeton perfoliatus)和狐尾藻(Myriophyllum verticillatum)为研究对象,对其种子休眠及萌发特性进行了初步研究.结果表明:狐尾藻种子休眠率较高但休眠易破除,通过切破种皮、低温层积、硝酸钾和赤霉素处理均可显著提高萌发率;篦齿眼子菜种子休眠性较强,低温层积及赤霉素浸泡处理能够显著提高其萌发率;穿叶眼子菜有较深的休眠特性,破除硬实、硝酸钾、赤霉素及低温层积处理均不能有效解除其休眠.     

Reliability of rayleigh wave dispersion curve obtained from f–k spectral analysis of microtremor array measurement

The paper utilizes previously developed microtremor simulation technique to evaluate the reliability of Rayleigh wave dispersion curve estimated by f–k spectral analysis of microtremor array measurement. The simulated microtremors are obtained for a fictitious (TEST) site. Attempt is also made to obtain the dispersion curve for two real sites (OHDATE and SKC) by inverse analysis of the microtremor array measurement using f–k spectra method. The estimated dispersion curve from simulated microtremors (TEST site) compares well with the theoretical dispersion curve, demonstrating the reliability of f–k spectra method and indicating that the estimated dispersion curve from microtremor measurements could be adequately used as the target for inverse analysis purposes. It is also demonstrated that the dispersion curve from microtremor measurements can be utilized to estimate the soil profile at OHDATE and SKC sites by inverse analysis. Results show that the theoretical dispersion curve of the fundamental mode of Rayleigh wave after the end of inverse analysis are in good agreement with the dispersion curve obtained by f–k spectral analysis of microtremor array measurement.