研究现今地壳运动和强震机理的一种方法

本文阐述了地壳形变与密度变化耦合运动的基本理论，分析了地壳形变与密度变化耦合运动引起的重力时间变化.断层位错运动是地壳形变与密度变化耦合运动的特殊形式.利用重复重力观测数据，选用稳健算法反演研究了1985-2002年滇西地震实验场主要断裂运动的时间分布和强震响应特征，获得了不同类型强震孕育过程图像.     

地壳“轧展”效应对地震成因的解释

目前多数地球科学家认为,地球的演化（包括地壳运动和地震活动）是地球内部因素所决定的,地球外部因素是次要的,作认为,地球虽是庞大并具有确定边界的星体,从银河系看来,地球不能作为单独的星体存在,它必须依附太阳系才能生存下去,在太阳系的演化中,太阳是太阳系演化的核心和主导,太阳自身的演化控制着行星的演化是必然的。因此,研究太阳活动,太阳能量向地球的输送及太阳与地球间能量（包括动量和质量）相耦合的物理机     

On the physical model of earthquake precursor fields and the mechanism of precursors’ time-space distribution (III) — anomalies of seismicity and crustal deformation and their mechanisms when a strong earthquake is in preparation

By studying the seismicity pattern before 37 earthquakes withM?6.0 in North China and the pattern of crustal deformation in the Capital Area from 1954 to 1992, some abnormal characteristics of these patterns before strong earthquakes have been extracted. A comparison has been made between the anomalies of these two kinds of patterns. From the results we can know the following. (1) Before a strong earthquake, the seismicity will strengthen and the crustal deformation rate will increase. (2) Several years before a strong earthquake, there will be seismic gaps and deformation gaps around the epicenter of the quake. (3) The dynamic parameters of patterns all show a decrease in information dimension. This means that the crustal deformation has become more and more localized with time and it gives an important indication showing that a strong earthquake is in preparation. At the end of the paper, the physical mechanisms of the abnormal patterns of seismicity and crustal deformation have been explained in a unified way in terms of the earthquake-generating model of a inhomogeneous strongbody in inhomogeneous media.     

地壳应变速率与地震活动关系的研究——从地壳应变场探索强震活动的场源关系

通过研究应变积累速度（亦即地壳应变速率）与地震活动的关系,籍以探讨从地壳应变场中寻找近期强震危险区的方法。针对不同地区的地震活动水平特别是地震复发期的差异,丁国瑜（１９８４）指出,各地震区大震复发时间的差异主要取决于由大震孕育环境和构造物理条件所决定的构造活动速率（亦即应变累积速率）的差异。从近几年地震实况的初步（阶段性）检验看,应用大地测量中垂直形变速率梯度判定强震危险区,具有较好的效果和可喜的前景。作者从理论分析、实际资料对比及计算机模型研究等方面来讨论地壳应变速率与地震活动水平之间的关系,进而为场源关系研究中的“场中求源”的科学思路探索可能的技术途径,并讨论了从应变速率场中寻找强震危险区的方法。     

Thermal evolution of Venus

The theory of three-dimensional and finite-amplitude convection in a viscous spherical shell with temperature and pressure dependent physical parameters is developed on the basis of a modified Boussinesq fluid assumption. The lateral dependences of the variables are resolved through their spherical harmonic representations, whereas their radial and time dependences are determined by numerical procedures. The theory is then applied to produce thermal evolution models for Venus. The emphasis is on illustrating the effects of certain physical parameters on the thermal evolution rather than proposing a specific thermal history for the planet. The main conclusions achieved in this paper are (1) a significant portion of the present temperature in the mantle and heat flux at the surface of Venus is probably owing to the decay of a high temperature established in the planet at the completion of its core formation, (2) the effective Rayleigh number of the mantle is so high that even the lower order modes of convection cool the planet sufficiently and maintain an almost adiabatic temperature gradient in the convecting region and high temperature gradients in the thermal boundary layers, (3) the convection is oscillatory with avalanche type properties which induces oscillatory features to the surface heat flux and the thickness of the crustal layer, and (4) a planetary model with a recycling crust cools much faster than those with a permanently buoyant crust.The models presented in this paper suggest that Venus has been highly convective during its history until ～ 0.5 Ga ago. The vigorous convection was bringing hot and fresh material from the deep interior to the surface and dragging down the crustal slags, floating on the surface, in to the mantle. The rate of cooling of the planet was so high that its core has solidified. In the last 0.5 Ga the vigour of convection diminished considerably and the crustal slags developed into a global and permanently buoyant crustal layer. The tectonic style on Venus has, consequently, changed from the recycling of crustal plates to hot spot volcanics. At the present time the planet is completely solid, except in the upper part of its mantle where partial melting may occur.     

Cumulate nodules as evidence for convective fractionation in a phonolite magma chamber

The physical mechanism by which chemical zonation develops in magma chambers has been controversial partly because unambiguous geological constraints have been lacking. The 11,000 years B.P. eruption of Laacher See Volcano produced a zoned tephra deposit and also ejected crystal-rich nodules which provide a snapshot of the materials crystallising at the magma chamber margins. New data on petrography and chemical compositions of nodules, their cumulate minerals and interstitial glasses are used to deduce the chemical evolution of the phonolite melt due to fractional crystallisation of the mineral assemblages. These data, together with those on the vertical zonation of the melt in the bulk of the chamber, are shown to be consistent with a model of stratification of the chamber by convective fractionation, in which a thin boundary layer of residual melt from fractional crystallisation ascends at the chamber side and accumulates at the roof. Crystallisation could have provided buoyancy to drive convection by enriching incompatible volatile components (mainly water) in the residual melt. Available fluid dynamic studies of single- and double-diffusive boundary layers are used to assess convection in the Laacher See chamber. The boundary layer is likely to have been: (1) laminar, which implies that the density gradient in the chamber steepened upwards; (2) in the counterflow regime, in which compositional and thermal layers flow in opposite directions; and (3) thin ( 10 cm), estimated from theory for a flat wall, suggesting that wall morphology could be important in determining boundary layer characteristics. Estimates of mass transfer rates due to this mechanism suggest that the chamber could have become stratified in a time of the order of 10³ years.     

Physical mechanism and mathematical modeling of earthquake ionospheric precursors registered in total electron content

The physical mechanism by which the regions with increased or decreased total electron content, registered by measuring delays of GPS satellite signals before strong earthquakes, originate in the ionosphere has been proposed. Vertical plasma transfer in the ionospheric F ₂ region under the action of the zonal electric field is the main disturbance formation factor. This field should be eastward, generating the upward component of plasma electromagnetic drift, in the cases of increased total electron content at midlatitudes and deepened minimum of the F ₂ layer equatorial anomaly. Upward plasma drift increases electron density due to a decrease in the O⁺ ion loss rate at midlatitudes and decreases this density above the equator due to an enhancement of the fountain effect (plasma discharge into the equatorial anomaly crests). The pattern of the spatial distribution of the seismogenic electric field potential has been proposed. The eastward electric field can exist in the epicentral region only if positive and negative electric charges are located at the western and eastern boundaries of this region, respectively. The effectiveness of the proposed mechanism was studied by modeling the ionospheric response to the action of the electric field generated by such a charge configuration. The results of the numerical computations indicated that the total electron content before strong earthquakes at middle and low latitudes is in good agreement with the observations.     

时纬残差分析方法及残差变化与强震的相关性

本文改进了时纬线差的分析方法,改善了残差异常与强震的相关关系,减少了虚报率。综合分析两种方法的异常,可对预测地震的远近大小提供更多信息。给出了1983年以来的时纬线差异常与全部5级以上地震的对应结果。最后讨论了可能的物理机制。     

1995年7月河北沙城地震群的震源断层

利用在怀来附近布设的中－欧合作小孔径数字化地震台网的资料,以相对定位法测定了１９９５年７月河北沙城地震群５４个地震的震源机制。从三维空间上分析了沙城震群的震源位置,震源机制,认为震源断层是ＮＥ向高倾角右旋走滑断层,并将其与宏观调查结果,地质活断层等作了对比。研究后认为：沙城震群发生在地壳浅层ＮＥ向地质断层延伸和扩展的新破裂上,属于现代地壳应力场作用下浅层断层延伸和扩展的新破裂上,属于现代地壳应力场     

中国大陆6级强震前的TIP研究

使用地震频次及其变化，平均震级和地震加权和的变化作为地震流函数研究了发生在１９８０－１９８９年期间中国大陆的６级强震前的地震概率增长时间。结果表明，１０次强震中９次发生在ＴＩＰ之后，而且最长ＴＩＰ为３．５年，ＴＩＰ总时空占有率为３０％。用同样的原则对新区进行识别也获得了好的结果。     

A direct method of calculating interval velocities and layer thicknesses from wide-angle seismic reflection times

Wide-angle reflections are now routinely recorded in high resolution explosion seismics to study the crustal structure. Use of Dix's hyperbolic approximation to the nonhyperbolic wide-angle reflection travel times causes major errors in the determination of interval velocities and layer thicknesses of a stack of horizontal velocity layers. Here we propose a layer stripping method to directly calculate the interval velocities and layer thicknesses in a vertically heterogeneous earth from the strong and reliable wide-angle reflected events. Synthetic reflection travel times, at wide-angle range, for a given velocity model, contaminated by some random errors, have been used to demonstrate the reliability of the algorithms to determine the interval velocities and thicknesses of various layers. The method has also been tested on two field examples along two deep seismic sounding (DSS) profiles with well identified wide-angle reflection travel times, which illustrates the practical feasibility of the proposed method.     

THE THREE DIMENSIONAL DENSITY STRUCTURE OF CRUST AND UPPER MANTLE IN THE CENTRAL-SOUTHERN PART OF LONGMENSHAN

In recent years, strong earthquakes of M_S8.0 Wenchuan and M_S7.0 Lushan occurred in the central-southern part of Longmenshan fault zone. The distance between the two earthquakes is less than 80 kilometers. So if we can obtain the inner structure of the crust and upper mantle, it will benefit us to understand the mechanism of the two earthquakes. Based on the high resolution dataset of Bouguer gravity anomaly data and the initial model constrained by three-dimensional tomography results of P-wave velocity in Sichuan-Yunnan region, with the help of the preconditioned conjugate gradient(PCG)inversion method, we established the three dimensional density structure model of the crust and upper mantle of the central-southern segment of Longmenshan, the spatial interval of which is 10 kilometers along the horizontal direction and 5 kilometers along the depth which is limited to 0~65km, respectively. This model also provides a new geophysical model for studying the crustal structure of western Sichuan plateau and Sichuan Basin. The results show obvious differences in the crustal density structure on both sides(Songpan-Ganzê block and Sichuan Basin)of Longmenshan fault zone which is a boundary fault and controls the inner crustal structure. In Sichuan Basin, the sedimentary layer is represented as low density structure which is about 10km thick. In contrast, the upper crust of Songpan-Ganzê block shows a thinner sedimentary layer and higher density structure where bedrock is exposed. Furthermore, there is a wide scale low density layer in the middle crust of the Songpan-Ganzê block. Based on this, we inferred that the medium intensity of the Songpan-Ganzê block is significantly lower than that of Sichuan Basin. As a result, the eastward movement of material of the Qinghai-Tibet plateau, blocked by the Sichuan Basin, is inevitably impacted, resulting in compressional deformation and uplift, forming the Longmenshan thrust-nappe tectonic belt at the same time. The result also presents that the crustal structure has a distinct segmental feature along the Longmenshan fault zone, which is characterized by obviously discontinuous changes in crustal density. Moreover, a lot of high- and low-density structures appear around the epicenters of Wenchuan and Lushan earthquakes. Combining with the projection of the precise locating earthquake results, it is found that Longmenshan fault zone in the upper crust shows obvious segmentation, both Wenchuan and Lushan earthquake occurred in the high density side of the density gradient zone. Wenchuan earthquake and its aftershocks are mainly distributed in the west of central Longmenshan fault zone. In the south of Maoxian-Beichuan, its aftershocks occurred in high density area and the majority of them are thrust earthquake. In the north of Maoxian-Beichuan, its aftershocks occurred in the low density area and the majority of them are strike-slip earthquake. The Lushan earthquake and its aftershocks are concentrated near the gradient zone of crustal density and tend to the side of the high density zone. The aftershocks of Lushan earthquake ended at the edge of low-density zone which is in EW direction in the north Baoxing. The leading edge of Sichuan Basin, which has high density in the lower crust, expands toward the Qinghai-Tibet Plateau with the increase of depth, and is close to the west of the Longmenshan fault zone at the top of upper mantle. Our results show that there are a lot of low density bodies in the middle and lower crust of Songpan-Ganzê Block. With the increase of the depth, the low density bodies are moving to the south and its direction changed. This phenomenon shows that the depth and surface structure of Songpan-Ganzê Block are not consistent, suggesting that the crust and upper mantle are decoupled. Although a certain scale of low-density bodies are distributed in the middle and lower crust of Songpan-Ganzê, their connectivity is poor. There are some low-density anomalies in the floor plan. It is hard to give clear evidence to prove whether the lower crust flow exists.     

Intermediate and deep earthquakes in Spain

Recent improvements in the seismological networks on the Ibero-Maghrebian region have permitted estimation of hypocentral location and focal mechanisms for earthquakes which occurred at South Spain, Alboran Sea and northern Morocco of deep and intermediate depth, with magnitudes between 3.5 and 4.5. Intermediate depth shocks, range from 60 to 100 km, with greater concentration located between Granada and Málaga. Fault-plane solutions of 5 intermediate shocks have been determined; they present a vertical plane in NE-SW or E-W direction. Seismic moments of about 10¹⁵ Nm and dimensions of about 1 km have been determined from digital records of Spanish stations.P-wave forms are complex. This may be explained by the crustal structure near the station, discontinuities in the upper mantle and inhomogeneities near the source. Deep activity at about 650 km has only 3 shocks since 1954 (1954, 1973, 1990). Shocks are located at a very small region. Fault-plane solutions show a consistent direction of the pressure axis dipping 45° in E direction. For the 1990 shock seismic moment is 10¹⁶ Nm and dimensions 2.6 km. TheP-waves are of simpler form with a single pulse. The intermediate and deep activities are not connected and no activity has been detected between 100 and 650 km. The intermediate shocks may be explained in terms of a recent subduction from Africa under Iberia in SE direction. The very deep activity must be related to a sunk detached block of lithospheric material still sufficiently cold and rigid to generate earthquakes.     

青藏高原东缘的地壳流及动力过程

黏滞性地壳流对地壳及上地幔变形作用及动力机制,是大陆新生代造山带的一个重要研究内容.青藏高原中下地壳存在部分熔融或含水物质的黏滞性流体,已为一系列地球物理及岩石学研究所证实.为研究青藏高原东缘地壳流的动力作用,本文用密集的被动源宽频带地震台的观测数据,反演了地壳上地幔精细速度结构和泊松比.研究表明,川西及滇西北高原的中地壳内普遍存在低速层,而高泊松比的地壳只分布在川西北地区.位于中地壳的黏滞性地壳流从青藏高原腹地羌塘高原流出,自北西向南东流入青藏高原东缘.这些黏滞性地壳流带动了上地壳块体水平移动,当它们受到刚强的四川盆地及华南地块阻挡时将发生分层作用,地壳流将分为二或更多分支不同方向的分流,向上的一支地壳流将对上地壳产生挤压,引起地面隆升,向下的一支地壳流将使莫霍面下沉加厚下地壳·黏滞性地壳流的运动在地壳中产生应变破裂发生强烈地震活动,地震的空间分布与震源机制也受到地壳流动力作用控制.     

龙门山断裂带地壳密度结构

研究龙门山及邻区地壳密度结构对于认识该地区地震活动性具有重要意义.根据龙门山及邻区( 100°～105°E,28°～33°N)的布格重力异常资料,选取了跨越龙门山断裂带的6条重力测线,在深地震测深资料约束下,使用Geosoft软件分别反演出了龙门山地区地下的沉积层、康拉德界面和莫霍面的深度分布.研究结果表明:龙门山断裂带两侧的地壳结构明显不同,西面高原地区沉积层较薄,大部分为基岩出露;而东边盆地沉积层明显较厚,多在6km以上.莫霍面和康拉德面在两侧均相对平缓,康拉德面从东部的大约24km增加到青藏高原山区的35km左右;莫霍面深度从东部盆地的大约42km增加到西部青藏高原的67km左右.龙门山断裂带整体表现为一条近SN向的陡变重力梯度带,并在其地壳内各界面均发生错断,莫霍面和康拉德面错断距离分别达6 ～ 7km和3～ 5km.该区地壳的这种陡变和不均匀性是导致地震活动性强烈的主要原因之一.     

珍贝—黄岩海山链热-重力均衡动态调节机制

珍贝—黄岩海山链作为我国南海的残留扩张中心, 对其研究具有重要的科学意义. 本文运用均衡学方法, 通过重力异常数据反演了过珍贝—黄岩海山链剖面的地壳界面变化, 同时计算了岩石圈热结构状态, 在此基础上建立了珍贝—黄岩海山链的岩石圈地温结构模型. 通过均衡分析方法, 对剖面上测点的海底地形数据进行了热均衡和重力均衡分析, 得到了热均衡和重力均衡形变量. 结果表明, 在珍贝—黄岩海山链高热流区域, 热均衡作用可以产生最大约0.55 km的形变, 其重力均衡形变范围为0.77—1.89 km. 热均衡通过改变海底地形和地壳物质密度不断作用于重力均衡, 重力又反过来作用于热均衡, 形成了热均衡-重力均衡动态调节机制.      

基于地壳双层结构假设的大陆强震孕震机制探讨

通过对地壳低速高导层和地震震源深度等有关资料的分析研究,在对地震发生前兆、特征及物理机制研究的基础上,初步提出了一个简单的地壳双层孕震模式:地壳存在有两层孕震结构,上层为主孕震层,是地震应力集中积聚和发震后释放应变能的地方,主要是岩石发生脆性变形层,深度在中上地壳范围内;下层为低速高导层,是诱发和加强上层脆性层岩石变形、应变能积累和破裂的主要因素之一,处于岩石脆—韧性变形的过渡带,强震的发生与本层有着密切的关系,深度在中下地壳范围内。同时基于该双层孕震假设讨论了强震的孕育环境及作用。     

华北强震地形变及重力前兆特征与孕震机理

通过分析华北地区强震地形变及重务前兆特征,发现震前在较大范围内存在与构造密切相关的动力不过程;在震中区介质状态发生变化。在此基础上提出了改进的组合－硬化模式：调整单元的运动,牵动较大范围构造运动的加及断层两侧深层介质密度变化;应力集中单元由于能量的不断积累,介质出现硬化过程。上述机理形成了前兆的各个阶段。     

On the Relationship Between the Crustal Strain Rate and Seismicity:Detect the Source-Field Relation from the Crustal Strain Field

The relationship between the strain cumulative rate (i.e., the crustal strain rate, or CSR in short) and seismic activity is analyzed to develop a new method to determine risky regions for strong shocks within recent years by the recorded crustal strain field. Seismic activity, especially the recurrence period, is different in different areas. Ding Guoyu (1984) pointed out that, for different seismic regions, the difference in the recurrence period of strong earthquakes is mainly controlled by their difference in the rate of the tectonic movement, which is controlled by the seismogenic environment and the tectonic conditions. The method of determining the risky regions for strong shocks from the gradient of vertical strain rates observed in a geodetic survey is preliminarily tested with the earthquakes in recent years; the results show that this method is effective and useful for earthquake prediction. The relationship between CSR and seismicity in a specific region is studied with strain theories, obse