唐山地震序列应力触发的粘弹性力学模型研究

根据前人对唐山地震破裂分布、地壳波速和粘性结构的研究,考虑局部应力场、孔隙流体压力和断层附近软介质的影响,计算了唐山地震产生的,投影到后续大余震断层面和滑动方向上的库仑破裂应力变化。结果表明,随后发生的滦县地震和宁河地震均发生在唐山地震产生的库仑破裂应力变化为正的区域。为研究唐山地震、滦县地震和宁河地震对后续小震的触发作用,根据前人对该地区构造应力场和地震破裂分布的研究,假定构造应力量值为10 MPa,求得了震源附近各处可能的小震震源机制。将上述3次地震产生的应力变化投影到可能的小震破裂面和滑动方向上,发现唐山地震、滦县地震和宁河地震产生的正库仑破裂应力变化的ldquo;蝴蝶rdquo;形分布与后续小震发生的空间分布具有较好的一致性,95%的余震发生在库仑破裂应力变化增加的区域,说明唐山地震序列中前面的大震对后续小震的发生起到了调制作用。该研究结果对大震后余震的危险性快速评估具有一定意义。如果大地震发生后能够快速确定详细的破裂分布和震源区域详细断层及滑动特性资料,本文方法可用来预测未来大余震的发震位置。      

日本1995年1月17日兵库县南部7.2级地震触发的区域地震活动

研究了日本1995年1月17日兵库县南部7.2级大地震触发的区域地震活动时空特征．研究区是主震周围几百公里的区域(33～37N, 133～138E)．研究区被划分为16个子研究区．对16个子研究区自1976年1月～1996年6月地震活动的研究结果表明:① 1995年1月17日兵库县南部7.2级主震触发了10个子区显著的地震变化,速率变化通过置信水平大于0.95的Z-统计检验,这些子区中最远的距主震震中约280 km;② 近震源区子区中的3个,地震活动在主震发生后1～5天内被触发,其它远离震源区的在几十天之后被触发;③ 被触发的区域地震活动最大震级为5.4,与主震诱发的最大余震震级相当;④ 区域地震活动触发的可能机制是主震引起的区域应力变化．      

几次复杂地震中不同破裂事件之间的应力触发问题

很多地震都具有复杂的震源破裂过程,作为一个近似,可以用两次或两次以上子破裂事件来描述,子破裂事件之间的时间间隔从几秒到数十天．为探讨不同的子破裂事件之间的关系,计算了邢台地震、唐山地震、青海共和地震和云南丽江地震的第一次子破裂事件所产生的静态库仑破裂应力变化．与断层滑动有关的静态应力的计算采用Okada（1992）给出的公式．结果表明:上述地震中第一次子破裂事件所产生的静态库仑应力变化均对后续的子破裂事件有触发作用,即后续的子破裂事件落在第一次子破裂事件引起的库仑破裂应力增加的区域中．      

2010年4月14日玉树MS7.1地震对余震的触发研究

研究了2010年4月14日青海玉树7.1级地震产生的静态库仑破裂应力变化对余震及区域地震活动空间分布的影响.在考虑震源区附近区域构造应力场的基础上,由震源区构造应力和主震破裂产生的应力叠加,计算得到最容易破裂的余震断层面,进而计算玉树主震在上述余震断层面上产生的库仑破裂静态应力变化.结果表明,库仑应力变化图像与余震分布较为吻合,说明玉树主震对大部分余震有触发作用.利用相同方法计算了玉树主震对周边更大范围内地震的应力触发影响,发现大部分区域地震活动空间分布与库仑应力变化图像基本一致,中小地震大多发生在正向触发区域,但量值较小.对余震断层面与主震完全一致的情形进行了同样的计算,与前述方法计算结果的对比显示,考虑区域应力场的计算方式所得库仑破裂应力变化图像,能够更好地与余震及区域中小地震空间分布图像相吻合.     

1976年唐山地震成因解释的一个新模式

指出已有对唐山地震成因解释的不足，综合分析了唐山地震前后的地震现象．认为唐山及周围地区1973～1975年ML4地震密集，而同一时空范围内３级和２级地震并不活跃是最为显著的地震活动特征，并认为这一特征应与区域高应力及地壳介质相对较为完整这一物理实质相对应．假定中国大陆周边的强震活动大体上反映着相应板块对中国大陆作用的强弱，推论20世纪70年代华北地区存在多个动力过程，为唐山地震提供了最基本的力源．进而提出多动力过程和地壳局部弱化模式，以解释唐山地震的成因．这一模式能够较好地解释与唐山地震有关的诸多地震现象．      

2014年于田MS7.3地震对后续余震和远场小震活动的动态应力触发

本文采用离散波数法,计算了2014年于田M_S7.3地震的断层破裂在近场和远场产生的库仑破裂应力变化,并结合地震活动特征,讨论了M_S7.3地震对后续余震活动和远场区域小震活动的动态应力触发作用.结果表明, ① M_S7.3地震产生的库仑破裂应力变化对其西南侧主体余震区的地震活动起到了抑制作用,这可能是本次M_S7.3地震序列余震活动水平不高的主要原因;距主震约30 km的北东方向余震区后续地震活动受到了主震产生的动态和静态应力变化的共同触发作用,动态应力变化峰值为2.78 MPa,静态应力变化为0.80 MPa,这与该区余震较为活跃相一致;距主震约45 km的北部余震区受到动态应力触发作用,应力变化峰值为0.72 MPa. ② M_S7.3地震产生的动态库仑应力变化空间分布呈非对称性,其中北东方向、北部余震分布与动态应力变化正值区存在相关性,从应力变化的角度解释了M_S7.3地震的后续余震空间活动特征. ③ M_S7.3地震在沙雅、伽师地区的远场接收点产生的动态应力变化峰值分别为0.09 MPa、0.1 MPa,对两个区域的小震活动具有动态触发作用.     

2001年1月2 6日印度古杰拉特(Gujarat)Ｍem>Ssub>7.8地震时空破裂过程

从全球数字地震台网的长周期记录中，选择了震中距小于90的27个台站的54个P波震相和44个S波震相资料．首先，用波形反演方法确定了2001年1月26日印度古杰拉特(Gujarat)ＭS7.8地震的地震矩张量、震源机制、震源时间函数和时空破裂过程等震源参数．通过矩张量反演，并根据Kutch Mainland断层的走向、地震烈度的空间分布、余震震源的空间分布和震害的空间分布，确认2001年1月26日印度古杰拉特ＭS7.8地震的发震断层的走向为92、倾角为58、滑动角为62，即一走向近东-西向、断层面向南倾斜、以逆冲为主的左旋-逆断层．这次地震所释放的地震矩为3.51020 Nm，矩震级ＭW＝7.6．然后，借助合成地震图，采用频率域求谱商的方法，得到了依赖于台站方位的27个P波震源时间函数、22个S波震源时间函数以及平均的P波震源时间函数和S波震源时间函数．对震源时间函数的分析表明，这次地震是一次连续的破裂事件，开始比较急遽，但结束比较迟缓，总持续时间约19 s．最后，以所提取的P波和S波震源时间函数为资料，采用时间域的反演技术得到了断层面上滑动的时空分布．滑动量在断层面上的静态分布表明，断层面上的最大滑动量约为7 m．断层面上的最大应力降约为30 MPa，平均应力降约为7 MPa．滑动量大于0.5 m的区域在走向方向长85 km，在断层面倾斜方向宽约60 km(相应地，在深度方向约51 km)．破裂向东扩展约50 km，向西扩展约35 km．滑动量大于0.5 m的区域的主要部分呈椭圆形，其长轴取向与断层滑动方向一致．表明此区域破裂扩展的方向即是断层错动的方向．这种现象对于走滑断层情形是多见的，但对逆冲断层情形却少见．断层面上初始破裂点以东、以上部分面积大于初始破裂点以西、以下部分的面积，这是破裂非对称性的表现，表明破裂具有自西向东、自下向上单侧破裂的特征．从滑动率随时空变化的快照可以看出，滑动率在第4 s达到最大值，此时滑动率约为0.2 m/s，滑动基本上发生在破裂起始点及其周围．从第6 s开始，起始点的破裂基本结束，破裂开始向外围扩展．破裂向西的扩展速度明显小于向东的扩展速度．在第15 s，这种环形的扩展基本结束．自16 s以后，主要是一些零星的破裂点分布在破裂区的外围．从滑动量随时空变化的快照看，破裂自起始点开始后，逐渐向四周扩展．主要的破裂(滑动量大于5 m的区域)在6～10 s，具有明显的自西向东、向上的单侧破裂特征．在第11～13 s，破裂的西端向西、向下有所扩展．整个破裂过程持续约19 s．在整个破裂过程中的平均破裂速度约为3.3 km/s．      

2007年云南宁洱Ms6．4地震震源过程

通过反演全球范围内20个地震台的宽频带波形资料,获得了2007年6月3日在云南宁洱发生的Ms6．4地震的矩张量解、震源时间函数和断层面上滑动随时间和空间的变化过程．根据反演结果,这次地震的标量地震矩为5．51×10^18 Nm,相当于矩震级Mw6．4．震源机制解中,最佳双力偶对应的节面Ⅰ的走向、倾角和滑动角分别为152°,54°和166°,节面Ⅱ的走向、倾角和滑动角分别为250°,79°和37°．结合震后考察得到的烈度等震线分布特征以及当地的地质构造,可以判定这次地震的发震断层的走向为152°,倾角为54°,滑动角为166°,是一次以右旋走滑为主的地震．从震源时间函数的形态来看,震源破裂持续时间为14s,地震矩的释放主要集中在前11s,在11-14s之间释放的地震矩很少．震源的时空破裂过程图像表明,破裂过程分为3个阶段,在前4S的时间段内,破裂主要沿着走向方向和朝深处发展;在4～7s间,破裂呈扇形向着深处扩展;在7s之后的时间段,破裂点比较零散．地震破裂总体上表现为双侧破裂方式,但在走向方向和深度方向上的滑动略占优势．破裂较强的区域呈菱形,长约为19km．地震断层面上最大滑动量为1．2m,平均滑动量为0．1m,最大滑动速率为0.4m／s,平均滑动速率为0．1m/s．由反演得到的静态位错模型计算的震中区地表位移场的特征与地震的烈度分布特征具有很好的一致性．     

大震地震波对云南地震活动的远场动态应力触发作用

基于云南数字地震台网记录,计算了印尼大震、昆仑山口西大震和古吉拉特邦大震的地震波在云南地区产生的远场动态库仑破裂应力变化,研究了大震产生的动态库仑破裂应力变化对云南地震活动的远场动态应力触发作用及其特征。结果显示,印尼大震产生的库仑破裂应力对云南地震活动有明显的触发作用,昆仑山口西大震产生的库仑破裂应力对云南地震活动有触发作用,但不显著,古吉拉特邦大震在云南地区产生的库仑破裂应力没有触发作用。最后分析和探讨了这3次大震触发作用差异的原因。     

阿拉斯加迪纳利M=7．9级地震瞬时变形引起的地震成核

推断导致地震触发的永久和动态(瞬时)应力变化常比地震本身释放的应力数量级小，这意味着触发发生在应力处于临界状态的断层上(Hill et al，1993；Gomberg et al，2000；Huc and Main，2003；Belardinelli et al，2003)。所以，触发的地震活动速率增加可能出现在加载速率最高以及可能由高温流体引起孔隙压力升高，使孽擦力减小促使断层破裂的区域(Coccoand Rice，2002；Streitand Cox，2001；Sturtevant et，al，1996)。本文说明2002年阿拉斯加州迪纳利M=7．9级地震在整个不列颠哥伦比亚省和美国西部触发了广泛的地震活动速率增加。地震波引起的动态触发应该沿着辐射能集中的破裂方向增加，我们利用迪纳利主震的地震波和新的高精度GPS记录证实了这个结论。这些观测结果和1992年加利福尼亚州兰德斯M=7．4级地震(Hill et al，1993)引起的触发仅在震级上具有可比性，并说明兰德斯地震的触发没有反映该地区的特性和地震活动性。可是，由迪纳利地震触发的地震活动速率增加的地区构造活动并不明显，意味着即使在周围应力加载速率低的地区，断层可能仍然处于应力临界状态，动态的触发很可能是普遍存在和无法预测的。     

On the synchronism in the events within the core and on the surface of the earth: the changes in the organic world and in the polarity of the geomagnetic field in the phanerozoic

The data on geomagnetic reversals are compared with the changes in the organic world and with the lower-mantle plumes. The times of the formation of plumes and the times of their appearance on the Earth’s surface relate to the intervals characterized by the different frequencies of geomagnetic reversals, i.e., there is no interrelation between the formation of plumes and the frequency of the changes in the geomagnetic field polarity. At the same time, a certain synchronism is observed between the frequency of the geomagnetic reversals and the boundaries of the biostratigraphic ages, i.e., the changes in the organic world in the long-period range. A hypothesis is proposed, which explains the change in the sign of the geomagnetic field by the combined effect of the irregular rotation of the internal core relative to the mantle and the changes in the slope angle of the axis of the Earth’s rotation, which, in turn, results in synchronous events on the Earth’s surface: the rates of changes in the organic world.     

On the maintenance of the axisymmetric part of the flow in the atmosphere

Summary The maintenance of the axisymmetric component of the flow in the atmosphere is investigated by means of a steady-state, quasi-geostrophic formulation of the meteorological equations. It is shown that the meridional variations in the time-averaged axisymmetric variables can be expressed as the sum of three contributions, one being due to the eddy heat transport, another to the eddy momentum transport, and a third to the convective-radiative equilibrium temperature which enters the problem through the specification of a Newtonian form of diabatic heating. The contributions by the large scale eddies are evaluated through the use of observed values for the eddy heat and momentum transports.The contributions from each of the three forcing mechanisms to the temperature and zonal wind fields are invstigated individually and found to be of about equal importance. The sum of the three contributions are also presented for the temperature, the zonal wind, the stream function associated with the mean meridional circulation and the corresponding vertical motion. Although the results fail to reproduce the main observed features of the lower stratosphere, they are found to be in good agreement with observations in the middle latitude troposphere. At any pressure level, for example, the computed mean zonal wind has a jet-like profile and the axis of the jet is found to slope to the south with height, as observed in the atmosphere.Based in part on a thesis submitted by the first author as partial fulfillment of the requirements for the Ph.D. degree at the University of Michigan. — Publication No. 194 from the Department of Meteorology and Oceanography, The University of Michigan.     

On the gravimetric undulations of the geoid and the deflections of the vertical

Summary The author mentions the aims of the World-wide gravity project he established in the Ohio State University in Columbus, in 1950. He outlines the practical procedure of the gravimetric computations of the undulationsN and the vertical deflection components and and emphasizes that only by the global international cooperation and additional gravity observations at sea carried out during the last decade it has been possible to gather to Columbus the needed gravity material. Since there exist still large gravimetrically unsurveyed areas it is of vital significance to study what gravity anomalies are best to be used for these regions. The given figures concerning the accuracy of theN, and , estimated theoretically and obtained in practice, indicate that in the gravimetrically well surveyed parts of the world like in Europe and the United States we can get gravimetrically on basis of existing gravity material theN-values with accuracy of about 5–10 meters, and and with the accuracy of about 1. The geoid undulationsN are already computed in Columbus for more than 6000 points of the northern hemisphere. The sample maps show the interesting geoid of Europe and vicinity between the latitudes 60° and 30° and longitude 5° W and 30° E, drafted on basis of more than 1000N-values computed at the corners of 1°×1° squares. It is interesting to realize that the geoid undulations in all this area are positive, the extreme values being between 40 and 50 meters. The geophysical significance of the geoid maps of this kind is pointed out.     

太平洋中部地磁台11年周期变化的研究

综合APIA、HONOLULU、PAMATAI 3个地磁台三分量的11年周期变化的形态后,可定性地提出该变化源于内场。如内源为核幔边界的电涡流,则其中心点应位于3台之间。     

Magnetic stratigraphy of the Ordovician in the lower reach of the Kotuy River: the age of the Bysy-Yuryakh stratum and the rate of geomagnetic reversals on the eve of the superchron

Until recently, the existing data prevented the geophysicists from accurately dating the Bysy-Yuryakh stratum, which outcrops in the middle reach of the Kotuy River, constraining the time of its formation to a wide interval from the end of the Late Cambrian to the beginning of the Silurian. The obtained paleomagnetic data unambiguously correlate the Bysy-Yuryakh stratum to the Nyaian regional stage and constrain its formation, at least a considerable part of it, by the Tremadocian. This result perfectly agrees with the data on the Bysy-Yuryakh conodonts studied in this work and yields a spectacular example of the successful application of paleomagnetic studies in solving important tasks of stratigraphy and, correspondingly, petroleum geology. Within the Bysy-Yuryakh stratum, we revealed a large normal-polarity interval corresponding to the long (>1 Ma) period when the geomagnetic reversals were absent. This result, in combination with the data for the Tremadocian and Middle–Upper Cambrian sequences of the other regions, indicates that (1) the rate of occurrence of the geomagnetic reversals on the eve of the Ordovician Moyero superchron of reversed polarity was at most one reversal per Ma; (2) the superchron does not switch on instantaneously but is preceded by a certain gradual change in the operation conditions of the dynamo mechanism which, inter alia, manifests itself by the reduction of the frequency of geomagnetic reversals with the approach of the superchron. This finding supports the views according to which a process preparing the establishment of the superchrons takes place at the core–mantle boundary.     

On the variation in the geomagnetic dipole over the geological history of the Earth

The global database on the paleointensity, containing determinations of the virtual dipole moment (VDM) for a stable (normal) regime of the geomagnetic field in a time interval of up to 3.5 Ga, is supplemented by new VDM determinations and analyzed. The field generation process started no later than 3–3.5 Ga (earlier data are absent) at the stage of the Earth’s core formation. Since that time, the dipole value has differed from its present value by no more than an order of magnitude, and the deviations that have already been detected tend toward smaller values. The distribution of VDM values in the time interval 0–400 Ma is bimodal, which apparently reflects the presence of two different generation levels of the geomagnetic field distinguished by a relatively large value (close to the present field value) and a relatively small value (approximately half as large as the present value). The total duration of decreased VDM values appreciably exceeds that of increased VDM values (179.1 and 28.6 Myr, respectively). On the whole, data on the paleointensity do not contradict the hypothesis about the dipole nature of the field over the last 400 Myr; however, the number of determinations at high paleolatitudes is too small to draw decisive conclusions on the validity (or invalidity) of the dipole field approximation based solely on paleointensity data.     

Comparison Between the Nutations of the Planet Mars and the Nutations of the Earth

The nutations of the planets Mars andEarth are investigated and compared. Alarge number of interior structureparameters are involved in the nutationcomputations. The comparison between the observations and the computationsprovides several constraints on these parmeters andtherefore allows a better understanding of the physics of the interior of theplanet. For the Earth, the high precision of the observations of the nutationshas led to a very good determination of interior properties of the planet. ForMars, observations of nutations are not yet available, and we review how theamplitude of the Martian nutations depends on the hypotheses consideredfor its interior. Although Mars is very similar to the Earth, its interior is not well known;for example, we don't knowif its core is liquid or solid. Only if the core is liquid,the Free Core Nutation (FCN) normal mode exists and can alter the nutationswhich are close to the resonance. From the observed geoids, it is known thatboth planets are not in hydrostatic equilibrium. The departure is larger forMars than for the Earth, and consequently, the implication of considering a convective mantle instead of a mantle in hydrostatic equilibrium described byClairaut's equation for the initial equilibrium state of the planet is largeron the Martian nutations than on the Earth nutations. The consequences of theuncertainty in the core dimensions are also examined and shown to be of a veryhigh influence for Mars if the core is liquid, due to the potential changes inthe FCN resonance. The influence of the presence of an inner core, which isknown to exist for the Earth, could be more important for Mars than for theEarth if the inner core is large. Due to the presence of Tharsis on Mars, thetriaxiality of this planet has, additionally, larger effects than on Earth.     

Evolutionary Changes in the Configuration of the Magnetic Field of the Solar Corona in the Epoch of the Minimum

Geomagnetism and Aeronomy - The evolution of the magnetic configuration in the solar corona is studied. The curvature of the K-corona helmets obtained from an analytical approximation of the...