汶川地震断裂带科学深钻WFSD-3附近上地壳S波分裂特征

为配合汶川地震断裂带科学钻探工程(WFSD),中国地震局地球物理研究所在四川省绵竹市天池乡和灌县—安县断裂附近分别架设了15套微震仪器和17个短周期地震台.基于WFSD-3附近的微震、短周期和区域台网的固定台站记录的近震数据,通过横波窗内S波分裂计算,得到其上地壳各向异性参数,即快波偏振方向和慢波的时间延迟,并分析了研究区的上地壳各向异性特征.结果显示,研究区大部分区域的快波偏振方向为NE向,与龙门山断裂带走向一致,但在研究区微震台阵布设小区域内,快波偏振方向表现出东西分区特征,东部为NE向,西部为NW向.上地壳各向异性主要是受到岩层中随应力分布排列的微裂隙和岩石或矿物结构的影响,研究区内快波偏振方向主要表现为NE方向,与断裂走向一致,反映了研究区上地壳各向异性主要受控于结构控制的各向异性,局部区域的快波偏振方向为NW向,与区域最大主压应力方向一致,说明区域应力场对研究区上地壳各向异性也有影响.通过分析微震台阵的归一化时间延迟随时间的变化情况可以反映区域应力场的变化情况.微震台阵的慢S波时间延迟在2012、2013年较为离散,在2014年有收敛的趋势,反映了强震后区域应力场逐渐稳定的趋势.     

龙门山断裂带北段南坝地区上地壳S波分裂特征

基于汶川科钻4号井孔（WFSD-4）附近的较小尺度的南坝微震台阵以及较大尺度的川西流动台站和区域台网的固定台站记录到的近震波形资料,通过横波窗内的S波分裂计算,分析了龙门山断裂带北段南坝地区的上地壳介质各向异性特征,并对区域应力场及构造特征展开讨论.S波分裂计算的结果显示研究区快波偏振方向主要表现为NE-NEE向,与北川断裂的走向一致,也与区域主压应力方向一致.小尺度密集分布的南坝微震台阵的计算结果进一步显示,靠近北川断裂的台站,其快波偏振方向与断裂走向一致,而距北川断裂较远的西北部台站的快波偏振方向与断裂走向不一致,反映了上地壳各向异性特征与地表的活动断裂结构密切相关.从南坝微震台阵的归一化时间延迟随时间的变化情况可以发现,在一个震中位于台阵布设范围内的地震事件发生前后,时间延迟有明显的变化,表明时间延迟随时间变化较为敏感,地震的能量影响了介质的性质.对比龙门山断裂带中段的科钻3号井孔周边地区的S波分裂计算结果,归一化时间延迟与本文结果一致,表明龙门山断裂带中段和北段的上地壳介质各向异性强度基本相同.     

青藏高原东北缘上地幔地震各向异性:来自SKS、PKS和SKKS震相分裂的证据

基于青藏高原东北缘甘肃区域台网41个宽频带地震台站的远震记录资料,通过PKS、SKS和SKKS震相的剪切波分裂分析,获取了台站下方介质的各向异性分裂参数,得到该地区上地幔各向异性分布图像,并结合GPS速度场和地壳剪切波各向异性分析青藏高原东北缘各向异性形成机制及壳幔各向异性特征.分析结果认为,在阿尔金断裂带西侧,各向异性快波偏振呈NWW-SEE方向,与断裂带走向有一定夹角,与塔里木盆地向柴达木盆地俯冲方向一致,说明该地区上地幔物质变形主要受古构造运动的影响,属于"化石"各向异性.在祁连山-河西走廊构造区,XKS快波偏振呈NW-SE方向,一致性较好,与区域断层走向方向相同;由区域小震的地壳剪切波分裂分析得到的地壳剪切波快波偏振在该区域呈NE-SW方向,与相对于稳定欧亚大陆GPS运动速率一致,地壳和地幔快波偏振方向的差异表明壳幔变形可能有不同的形变机制.在陇中盆地及其周缘,由于处于活跃青藏地块与稳定鄂尔多斯地块之间的过渡带,相对于其他区域具有更加复杂的构造背景,地壳快波偏振和地幔快波偏振总体上呈NWW-SEE方向,说明壳幔变形机制可能相同;但不同台站结果之间存在一定离散性,推测是由于受局部构造特征差异性造成.     

青藏高原东南缘三江地区上地壳各向异性初步研究

三江地区位于青藏高原东南缘,川滇菱形块体的西侧.受陆-陆板块碰撞作用的影响,构造活动强烈,地震活动频繁.为研究该区的应力环境、构造特征及二者间的关系,本文使用三江流动线性地震台阵SL-Array（2016-12—2017-05）和国家固定地震台网（2015-01—2017-12）记录到的地震波形资料,运用剪切波分裂分析技术研究三江地区上地壳各向异性研究.计算得到该区域快剪切波优势偏振方向为NNW向,与区域主压应力方向一致.结果显示各向异性特征有分区性,以维西—乔后断裂和小金河—丽江断裂为界,将线性台阵划分为A、B和C三个区域.A区快波优势偏振方向表现出与区域主压应力方向的一致性.B区局部构造复杂,快波优势偏振方向表现为近NS向.C区结果比较离散,无明显快波优势方向.自西向东,研究区域快波优势偏振方向表现为NNW至近NS向的变化趋势.计算得到研究区域慢剪切波时间延迟为2.8±1.7 ms·km^-1,其中B区最大,A区最小,反映了该区地壳介质各向异性强度的不均匀分布,也揭示了区域构造复杂程度与地震各向异性强度的关系.     

腾冲地区近震S波分裂研究

本文利用布设在云南腾冲地区的15个固定和流动地震台站记录的近震波形数据,采用剪切波分裂分析方法得到了593对高质量的各向异性分裂参数.结果显示,腾冲火山区地震台站下方的近震各向异性的慢波延迟时间为0.02～0.37 s,平均延迟时间0.2 s.结合已有接收函数地壳各向异性研究结果,推测研究区地壳各向异性的主要贡献源自中上地壳.研究区不同台站的快波偏振方向变化很大,似乎反映了构造和区域应力场的共同作用.其中腾冲火山断裂西侧多数台站的快波偏振方向呈近N-S向,而东部多数台站的快速偏振方向呈NE-SW向,与区域主压应力方向一致,暗示研究区中上地壳各向异性主要是受主压应力引起定向排列的裂隙所致.基于近震走时得到的研究区平均V_P/V_S为1.68,推测腾冲火山区地壳应力场的局部变化可能与上地壳中富含气体的中酸性岩浆膨胀活动有关.另一方面,在腾冲火山区外围个别台站(MIZ、MZT)观测到了快波偏振方向与主压应力、已知断层等构造走向不一致的现象,暗示其各向异性是构造或构造和区域应力场共同作用的结果.     

2020年济南ML4.4地震前后近场S波分裂变化特征

利用山东省地震局布设的流动地震台站和区域台网固定台站记录的2020年济南M_L4.4地震序列资料,通过横波窗内的近震S波分裂计算,获取了各台站上地壳各向异性空间分布和时间变化的结果,对上地壳变形特征以及S波分裂参数(快S波偏振方向和慢S波时间延迟)的时空变化现象进行了分析.结果表明,上地壳各向异性在空间分布上存在分区特征,并且表现出随时间的趋势性变化特征.发震断层东侧的台站LJIN和L3704只有一个快S波优势偏振方向NW向,其上地壳各向异性主要受控于断裂结构;发震断层西侧的L3705台,以及震源区外的FEIC和PYI台主要表现为NE向,反映出区域应力场控制的上地壳各向异性;而靠近长清断裂带的L3706和CHQ台,表现出两个近于垂直的快S波偏振方向,反映了不同射线路径上S波分裂的不同特征.震源区内台站的慢S波时间延迟平均值大于震源区外各台站,表明济南M_L4.4地震孕育过程中震源区的应力积累水平高于周边地区.震源区内CHQ台的慢S波时间延迟在主震前较低,在主震后迅速增大,一段时间后逐渐减小并趋于恢复到震前水平,反映了济南M_L...     

中国地震科学台阵两期观测资料近场记录揭示的南北地震带地壳剪切波分裂特征

利用中国地震科学台阵第一期(2011-01-2014-06)及部分中国地震科学台阵第二期(2013-02-2015-12)的流动地震台阵记录到的小震波形资料,运用剪切波分裂系统分析(SAM)方法,分析南北地震带的地壳各向异性,对剪切波分裂参数所反映的区域应力环境及构造特征,以及区域内主压应力方向与断裂分布的关系展开讨论.研究结果表明,南北地震带快剪切波偏振方向自北向南由NE向逐渐转变为NNW向,与南北地震带区域主压应力的方向变化具有一致性.区域内分布的大量NE及WNW或NW向断裂构造同样对快波偏振方向有比较大的影响,位于走滑断裂附近的台站,其快波方向与断裂走向大致平行,部分位于走滑断裂附近的台站其快波方向几乎垂直于断裂走向,而与构造应力场方向一致性较好.个别台站表现出复杂快波优势方向特征,反映出研究区内构造环境的复杂性.慢波时间延迟结果显示,南北地震带南段的平均时间延迟高于北段,反映了受印度板块和欧亚板块的碰撞挤压作用,南段地壳介质各向异性程度更大,构造变形更加剧烈.对比南北地震带上地幔各向异性特征,推测在川滇菱形块体内部可能存在复杂的壳幔耦合现象,地壳剪切波分裂除了反映区域应力特征,还可以揭示出区域构造信息.     

南迦巴瓦地区近震剪切波分裂研究

南迦巴瓦地区位于喜马拉雅山脉构造结东端,地壳各向异性具有其特殊性和复杂性.本研究收集了布设在南迦巴瓦地区的流动地震台站和固定台站于2017—2019年记录的近震波形数据,利用剪切波分裂偏振分析方法得到了研究区内10个地震台站的剪切波分裂参数快波偏振方向和慢波延迟时间,分析了剪切波分裂参数在空间上和时间上的分布,探讨了南迦巴瓦地区上地壳地震各向异性特征,以及2017年11月18日米林M_S6.9地震发生后地壳应力状态的变化.结果显示,南迦巴瓦地区快波偏振优势方向有两个,即北西向和北东向.位于米林断裂、嘉黎断裂、墨脱断裂和雅鲁藏布江断裂附近的台站,快波偏振方向基本平行于断裂走向.研究区域上地壳各向异性方向受区域应力场的控制,另外还受到断裂构造的影响.南迦巴瓦地区慢波延迟时间均值为4.13 ms·km^-1.距离米林地震主震震中位置较近的L0230台站慢波延迟时间较大.本研究没有观测到快波偏振方向在时间上的规律性变化.L0230台站的慢波延迟时间在米林M_S6.9地震之后表现出明显的降低,反映了大地震之后地壳应力的释放与调整.     

鄂尔多斯块体西缘地壳介质各向异性:从银川地堑到海原断裂带

通过收集鄂尔多斯块体西缘固定地震台网2010年6月至2017年8月的近场地震资料,选择符合剪切波分裂分析的14个台站记录的共137个有效事件波形,得到了剪切波分裂参数,即快剪切波（简称快波）偏振方向和慢剪切波（简称慢波）时间延迟.结果表明,研究区的快波偏振方向和慢波时间延迟具有明显的分区特征,快波偏振方向主要与构造应力场方向或者断层走向大体一致.鄂尔多斯西缘紧邻块体边界的台站,快波偏振方向自北向南呈现NS、NNE、NE向的变化,与青藏高原东北缘主压应力方向变化基本一致.银川地堑东西两侧的快波偏振方向有差异,东侧区域主要受青藏高原NNE向挤压和黄河-灵武断裂共同影响,而西侧区域可能受到阿拉善块体与鄂尔多斯块体之间的NW方向的主张应力和阿拉善块体内部应力分布的影响;鄂尔多斯块体、阿拉善块体与青藏高原的交汇区快波优势偏振方向为NE向,与青藏高原东北缘主压应力方向一致;海原断裂带及以南区域快剪切波优势偏振方向为WNW向,与断裂走向基本一致,较好的说明了海原断裂带为活跃的活动断裂.构造与断裂分布都是控制快波偏振方向的主要因素,走滑断裂上的台站快波偏振方向与断裂走向一致,表明这些台站主要受到断裂的强烈影响;走滑断裂附近的个别台站快波偏振方向呈现与构造应力场一致的方向,表明几乎没有受到断裂的影响.鄂尔多斯、阿拉善与青藏高原的交汇区平均时间延迟高于其他地区,反映了青藏高原在NE向运动过程中,受到稳定的鄂尔多斯块体阻挡作用,导致了交汇区地壳介质各向异性程度增加.以海原断裂带到六盘山断裂带为界,其两侧区域的各向异性差异性明显,揭示了应力与介质特性的差异,暗示其邻近区域,特别在海原断裂带东端到六盘山断裂带与鄂尔多斯块体西缘交汇区域,可能有较高的强震危险背景.本研究还对该区域的地壳和上地幔的耦合问题进行了初步讨论.     

华南块体中部上地壳剪切波分裂特征

本文采用剪切波分裂分析方法,使用华南流动地震台阵(HN Array)和国家地震台网记录的近场波形资料,分析华南块体上地壳介质的地震各向异性特征.结果显示,华南块体快剪切波优势偏振方向整体表现为近东西向,与区域最大主压应力方向基本一致.结合华北块体地震各向异性特征,推测太平洋板块西向俯冲和菲律宾板块北西向俯冲作用,是中国大陆东部中上地壳各向异性的重要构造背景,导致中国大陆东部表现出近东西向的构造应力特征.苏鲁—大别造山和江南造山带的快剪切波偏振方向发生明显区域性变化,揭示了块体之间相互作用产生局部复杂构造特征和应力分布特点.研究显示,各向异性观测结果会因地壳结构的多样性变得复杂,通过江汉盆地的地震各向异性特征揭示出盆地内存在不同尺度的应力场.研究认为,岩石组分对中上地壳地震各向异性有重要贡献,华夏块体东部广泛存在的变质岩导致东西华夏块体各向异性偏振方向存在明显差异.     

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

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.     

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

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

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.     

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.     

Structure of the poleward wall of the trough and the inclination of the geomagnetic field above the EISCAT radar

A special high-resolution routine of the EISCAT radar has been used to investigate the structure and development of the poleward wall of a deep trough in electron density. The feature was tracked by the radar during a 7-hour period under very quiet geomagnetic conditions. The field-aligned nature of the structure enabled an estimate to be made of the inclination of the geomagnetic field above EISCAT that was in good agreement with the current model. Observations of narrow field-aligned enhancements in electron temperature demonstrated that the wall of this trough is a dynamic feature, reforming regularly as the electron density responds on a time scale of tens of minutes to energy input from soft-particle precipitation.