(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

中强地震前后能量分布的自仿射特性

根据简化的地震波传播理论,“漫化”地震能量在空间的分配,减弱能量空间分布概率的极端情形及消除统计域为“空”的情形,由此构造类似于二维ｆＢｍ的地震能量空间分布概率曲面。与起伏不平的地形相似,由于基在水平方向及垂直方向的变化不是等比例的,因而这一分形曲面可能是统计自仿射的。     

北京及邻区微震活动的初步分析

通过对北京遥测地震台网近年来记录到的北京及邻区地震的震中分布,地震活动频度及能量释放强度的分析,得到本区地震活动在时间分布上具有“聚堆性”。在年发震频度,地震强度和能释放方面均具有双峰值特征,并且具有较好的一致性和同步性,在空间分布上具有条带特征,且形成北东～南西和北西～南东的两条相互交汇的条带。又通过统计分析得到本区发震概率最大的时间段是每年的１０月前后,而地震主要发生在北西～南东带上。     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

SOME NEW UNDERSTANDINGS OF LAND-SURFACE PROCESSES IN ARID AREA FROM THE HEIFE*

The survey for the HEIFE(Atmosphere-Land Surface Processes Experiment at Heihe River Basin,Western China) is given in the paper.The following basic subjects for land-surface processes in arid areas are studied:(1) the general characteristics of the energy budget on ground surface in arid areas;(2) the parameterization of the land surface processes;(3) the interaction between oasis and its desert circumstances,a special phenomenon in arid areas.The analysis shows that the sensible heat flux in the surface energy budget is in the majority,and the latent heat flux may be neglected.The influence of atmospheric stratification stability on the turbulent transfer of energy and substance must be considered in parameterization of land surface processes in arid areas.The "cold island effect" phenomenon in oasis and the "humidity inversion" phenomenon in desert near oasis are the result of the interaction between them.The results would improve the understanding of land surface processes in arid areas.     

Energy balance of simple elastodynamic sources

Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic     

Some comparisons between mining-induced and laboratory earthquakes

Although laboratory stick-slip friction experiments have long been regarded as analogs to natural crustal earthquakes, the potential use of laboratory results for understanding the earthquake source mechanism has not been fully exploited because of essential difficulties in relating seismographic data to measurements made in the controlled laboratory environment. Mining-induced earthquakes, however, provide a means of calibrating the seismic data in terms of laboratory results because, in contrast to natural earthquakes, the causative forces as well as the hypocentral conditions are known. A comparison of stick-slip friction events in a large granite sample with mining-induced earthquakes in South Africa and Canada indicates both similarities and differences between the two phenomena. The physics of unstable fault slip appears to be largely the same for both types of events. For example, both laboratory and mining-induced earthquakes have very low seismic efficiencies where _a is the apparent stress and is the average stress acting on the fault plane to cause slip; nearly all of the energy released by faulting is consumed in overcoming friction. In more detail, the mining-induced earthquakes differ from the laboratory events in the behavior of as a function of seismic momentM ₀. Whereas for the laboratory events 0.06 independent ofM ₀, depends quite strongly onM ₀ for each set of induced earthquakes, with 0.06 serving, apparently, as an upper bound. It seems most likely that this observed scaling difference is due to variations in slip distribution over the fault plane. In the laboratory, a stick-slip event entails homogeneous slip over a fault of fixed area. For each set of induced earthquakes, the fault area appears to be approximately fixed but the slip is inhomogeneous due presumably to barriers (zones of no slip) distributed over the fault plane; at constant , larger events correspond to larger _a as a consequence of fewer barriers to slip. If the inequality _a / 0.06 has general validity, then measurements of _a =µE _a /M ₀, where is the modulus of rigidity andE _a is the seismically-radiated energy, can be used to infer the absolute level of deviatoric stress at the hypocenter.     

The upper mantle structure of the yakutian kimberlite province,eastern Russia

Deep seismic sounding in the region of the Mirnyi kimberlite field indicates that the boundary velocity of the uppermost mantle is elevated (v _b=8.6–8.8 km/sec) and extremely variable near the Mir kimberlite pipe. These velocity heterogeneities are probably associated with the kimberlite magmatism and may be useful in the identification of other kimberlite fields.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

Textural segmentation of digital rock images into bedding units using texture energy and cluster labels

The texture of digital rock images, as recorded, for instance, with borehole imaging devices, is shown to reflect different bedding types. Textural segmentation of borehole images, therefore, subdivides the recorded sequence into bedding units. We show that a textural segmentation algorithm based on the concept of texture energy achieves good results when compared with synthetic as well as real data in which petroleum geologists have performed zonations on cores. Texture energy involves filtering of the original image with a set of texture sensitive masks. The filtering is done as a finite convolution over the size of the masks. On the resulting images the variance is computed over a relatively large sliding window, which, in its practical implementation, covers the full width of the image. The resulting nine one-dimensional curves are then clustered hierarchically into a user-determined number of image texture or lithological bedding classes. Principal component analysis previous to clustering can be used to reduce redundancy in the data. A recurring and relatively ill-defined problem in this field are macro-textures, i.e., the cyclic interbedding of two or more bedding types. We show that sliding Fourier transforms and variable mask scale can successfully address the zonation of macro-textures. In general, the method gives best results with mask sizes equivalent to 2–4 centimeters, reflecting the length scale at which the investigated geological bedding seems to have its highest variation.