龙门山构造带壳幔电性结构特征及其与汶川、芦山强震关系

青藏高原东缘龙门山构造带是研究青藏高原地壳物质向东侧向挤出的焦点地区.为探索龙门山构造带活动构造特征及其与发震构造的关系,本文通过布置垂直龙门山构造带南段芦山地震震源区的大地电磁测深剖面,运用多种数据处理手段,得到研究区可靠的电性结构,并通过与已有龙门山中段和北段剖面进行对比分析.研究表明：（1）青藏高原东缘岩石圈存在明显的低阻异常带--松潘岩石圈低阻带,该低阻异常带沿龙日坝断裂-岷山断裂-龙门山后山断裂分布,形成松潘-甘孜地块向扬子地块俯冲的深部动力学模式,通过统计研究区的历史强震,发现震源主要沿低阻异常带东侧分布,同时,低阻异常带也是低速度、低密度异常带,松潘岩石圈低阻带可能是扬子地块的西缘边界;（2）青藏高原物质东移过程中,受到克拉通型四川盆地的强烈阻挡,龙门山构造带表层岩块和物质发生仰冲推覆,表现为逆冲推覆特征的薄皮构造,中下地壳和上地幔顶部物质向龙门山构造带岩石圈深部俯冲,印支运动晚期,扬子古板块持续向华北板块俯冲,在上述构造运动作用下,呈现出刚性的上扬子地块西缘高阻楔形体向西插入柔性青藏块体的楔状构造;（3）根据电性结构推断,芦山地震受到深部上里隐伏壳幔韧性剪切带向上扩展的影响,构成芦山地震的深部主要动力来源;汶川地震的发生,在龙门山南段形成应力加载区,是触发或加快芦山地震孕育发生的另一个动力来源.     

利用虚震源反射成像法研究江南造山带东段上地壳结构

江南造山带是华南地区扬子地块与华夏地块碰撞的产物,其地壳构造记录了两地块的碰撞过程,研究江南造山带的地壳构造有助于重建扬子地块与华夏地块的碰撞过程。本研究在江南造山带上布设了两条流动地震台阵,利用虚震源反射法提取其所记录的远震事件初至P波在地表的反射波（PPdp）波形,重构了研究区内两条测线下方的上地壳结构。结果显示:江绍断裂两侧上地壳沉积层的厚度变化明显,推断该断裂是扬子地块与华夏地块的东边界;相较华夏地块,江南造山带与扬子地块的层位连续性更强,符合江南造山带先与扬子地块合为整体后再与华夏地块碰撞的多期构造过程及其对应产生的亲扬子地块属性;江绍断裂西北侧的地层不整合以及赣东北断裂区域的断陷构造,可为了解古华夏洋向扬子陆块俯冲及碰撞和随后的构造运动过程提供参考依据。      

青藏高原东缘及四川盆地的壳幔导电性结构研究

自从2008年M_S8.0级汶川大地震发生以来,青藏高原东缘便成为地质与地球物理研究的热点区域.该区域的龙门山断裂带标志着青藏高原东缘与四川盆地的边界.汶川地震即发生于龙门山断裂带内的映秀-北川断裂上.该地区现有的研究工作多集中于青藏高原东缘及四川盆地的西部,对四川盆地东部构造情况的研究目前较少.在SinoProbe项目的资助下,完成了一条跨越青藏高原东缘及整个四川盆地的大地电磁测深剖面.该剖面自西北始于青藏高原内部的松潘-甘孜地块,向东南延伸穿过龙门山断裂带、四川盆地内部及四川盆地东部的华蓥山断裂,最终止于重庆东南的川东滑脱褶皱带附近.维性分析表明剖面数据整体二维性较好,通过二维反演得到了最终的电性结构模型.该模型表明,从电性结构上看,沿剖面可分为三个主要的电性结构单元,分别为:浅部高阻、中下地壳低阻的松潘-甘孜地块,浅部低阻、中下地壳相对高阻的四川盆地,以及华蓥山以东整体为高阻特征的扬子克拉通地块.龙门山断裂带在电性结构上表现为倾角较缓、北西倾向的逆冲低阻体,反映了青藏高原东缘相对四川盆地的推覆作用.其在地下向青藏高原内部延伸,深度约为20 km左右.在标志逆冲推覆滑脱面的低阻层下存在一电性梯度带,表征着低阻的青藏高原中下地壳与高阻的扬子地壳之间的电性转换.位于四川盆地东边界的华蓥山断裂在电性结构上表现为一倾向为南东向的低阻体插入高阻的扬子克拉通结晶基底,切割深度约为30 km左右.这一结构反映出华蓥山向西的推覆作用.在电性结构模型的基础上,进一步讨论了青藏高原东缘的壳内物质流、青藏块体与扬子块体的深部关系以及青藏高原东部的隆升机制等构造问题.     

青藏高原东北缘祁连山造山带至阿拉善地块壳幔电性结构研究

在青藏高原东北缘祁连山造山带至阿拉善地块之间完成了一条372km的大地电磁剖面,通过二维反演计算,获得了沿剖面180km深的壳幔电性结构模型,结合研究区地质和地球物理资料开展综合分析,研究结果表明:(1)剖面自南向北所经过的祁连山造山带、走廊过渡带和阿拉善地块对应3种壳幔电性结构模型:东祁连壳幔高-低-高阻似层状电性结构、河西走廊壳幔低阻带状电性结构和阿拉善南缘壳幔高-低-高阻层状电性结构.(2)剖面所经过的主要断裂带在电性结构上表现为低阻异常带或电性梯度带,并且止于中上地壳或消失于下地壳低阻层中.除这些分布于中上地壳的断裂系统以外,在下地壳至上地幔顶部还存在两条切割莫霍面的壳幔韧性剪切带:西华山北缘壳幔韧性剪切带和阿拉善南缘壳幔韧性剪切带.其中,西华山北缘壳幔韧性剪切带可能是1920年海原8.6级地震发生的深部背景之一;而阿拉善南缘壳幔剪切带可能是卫宁北山燕山晚期和喜山期幔源岩浆上升到地壳浅部或喷出到地表的通道,为在该区域寻找晚中生代至新生代含矿隐伏岩体提供了深部电性结构依据.(3)由若干形状不规则、彼此不相连的"碎块状"极高阻块体组成的中上地壳与"似层状"的中下地壳低阻层共同构成的地壳电性结构,是引起青藏高原东北缘强烈破坏性地震最佳的地壳电性结构组合之一.印度板块向欧亚板块俯冲碰撞楔入引起青藏高原块体向北东方向运移与阿拉善地块向南的俯冲碰撞楔入,是青藏高原东北缘强震活动带产生的动力学背景.     

华北坳陷中南部深部结构大地电磁探测

为研究华北坳陷中南部河南地区的深部电性结构和该区3次破坏性地震的深部孕震背景，采用宽频带大地电磁测深法对该区深部的电性结构进行探测。使用NLCG二维反演方法对TE、TM模式进行联合反演，最终获得研究区的深部电性结构。结果显示，以新乡-商丘断裂为界，其北部的北华北坳陷中的内黄隆起、东濮凹陷的电性结构相对简单，为低阻、高阻2层结构，与该区地壳由沉积盖层和基底整体性较好的坚硬地块结晶基底组成相对应。新乡-商丘断裂以南的南华北坳陷中的构造单元深部的电性结构相对复杂，纵向上表现为低—高—低3层结构，横向上高阻、低阻交互出现，分析其可能与受到扬子板块与华北板块南缘相互俯冲碰撞影响有关。而北华北坳陷因边界断裂———新乡-商丘断裂的调控，受到扬子板块和秦岭-大别造山带造山作用的影响较小。1342年和1918年河南通许2次破坏性地震的震源位于太康隆起下方的高、低阻边界交接区域，1737年河南封丘■级地震的震源位于壳内高、低阻梯度带附近。     

青藏高原东缘龙门山逆冲构造深部电性结构特征

通过对汶川地震前观测的碌曲—若尔盖—北川—中江大地电磁剖面的数据处理和反演解释,揭示了沿剖面的松潘—甘孜地块、川西前陆盆地、龙门山构造带及秦岭构造带50 km深度的电性结构特征及相互关系,表明青藏高原东缘向东挤压,迫使向东流动的地壳物质沿高原东缘堆积,并向扬子陆块逆冲推覆.龙门山恰好位于松潘—甘孜地块与扬子陆块对挤部位,主要受松潘—甘孜地块壳内高导层滑脱和四川盆地基底高阻体阻挡的约束,地壳深部存在着西倾且连续展布的壳内低阻层,表明龙门山深部确实存在着逆冲推覆构造,其逆冲断裂系中的三条断裂不仅以不同的倾角向西北倾斜,并且向深部逐渐汇集,但茂县—汶川断裂可能在深部与北川—映秀断裂是分离的.龙门山两翼的四川盆地和松潘甘孜褶皱带的电性结构既具有明显差异性,又具有一定的相关性.四川盆地显示巨厚的低阻沉积盖层和连续稳定的高阻基底的二元电性结构,而松潘—甘孜地块则表现为反向二元结构,即上部大套高阻褶皱带,下部整体为低阻的变化带,龙门山逆冲构造带本身又表现为松潘地块逆冲上覆在四川盆地之上,构成上部高阻褶皱带、中部低阻逆冲断裂带和底部盆地高阻基底的三层电性结构.对比龙门山逆冲构造断裂带的西倾延伸上下盘两侧的两个反对称的二元电性结构,松潘区块深部推断的结晶基底与龙门山断裂带下盘推断的下伏盆地结晶基底又存在某种内在对应关系,推断可能存在一个西延至若尔盖地块的泛扬子陆块.因此,龙门山构造带地壳电性结构研究对于揭示青藏高原东缘陆内造山动力过程,探索汶川大地震的深部生成机理都具有重要意义.     

中朝与扬子地块结合带东部的卫星重力异常特征研究

中朝板块与扬子板块碰撞结合带过黄海向东延伸方式一直是国内外学者关注的问题.本文采用最新的地球重力场模型(GGM03C)计算了中国东部海域及朝鲜半岛的360阶卫星布格重力异常,通过异常分离、水平梯度及其垂直导数等位场资料处理和边界识别技术的应用,获得与研究区地块边界及断裂相关的各种信息.结合地质资料,对中国东部海域及朝鲜半岛进行地质构造单元划分与断裂系统分析,探讨了中朝与扬子地块结合带在海区位置与东部延伸模式.对京畿地块的深部结构及地质构造的认识是探讨中朝-扬子地块结合带在朝鲜半岛所在位置的关键.分析认为,五莲-青岛-荣城断裂和嘉山-响水断裂构成中朝-扬子地块结合带的海区段,南黄海南部断裂带与绍兴-勿南沙断裂带构成中朝-华南地块结合带;在京畿地块中部可能存在近南北向的断裂带,与临津江构造带、沃川构造带和上述两个结合带共同构成中朝、扬子与华南地块结合带;南黄海属于扬子地块,而北黄海与朝鲜半岛大部分仍属于中朝地块,其南界位于济州岛南缘断裂带附近,济州岛南缘断裂带南部海域为华南地块.     

切穿兴蒙造山带的深断面——岩石圈电性结构及改造模式

兴蒙造山带自早古生代形成稳定的佳蒙地块以来,受到古亚洲洋闭合、蒙古—鄂霍茨克洋闭合,尤其是太平洋板块俯冲构造事件的影响.为了揭示兴蒙造山带岩石圈电性结构及改造模式,本文基于横切额尔古纳地块、兴安地块、松辽地块及佳木斯地块约1300 km的大地电磁测深剖面,对81个长周期测点进行非线性共轭梯度二维反演,获得了兴蒙造山带东西向完整的电性结构模型.结果显示：对于岩石圈尺度的结构,不同构造单元具有明显的横向差异性,额尔古纳地块岩石圈整体呈低阻特征;兴安地块岩石圈呈高阻特征;松辽地块的松辽盆地和小兴安岭岩石圈分别呈低阻-高阻特征;佳木斯地块岩石圈呈高、低阻交替的特征.在软流圈尺度,兴安地块和松辽地块存在大规模高导异常C2、C3,呈“U”型状与上方岩石圈相连,连接处对应东北地区盆山结合部位或主要的缝合带.研究表明兴蒙造山带岩石圈遭受了强烈的改造作用,结合区域构造背景,本文提出“自下而上”的岩石圈改造模式：在地幔转换带处水平滞留的太平洋板块持续扰动兴蒙造山带下方的软流圈,导致软流圈物质的上涌以及古缝合线的活化,形成的薄弱带为软流圈物质的侵入提供了有利通道,太平洋板块的俯冲深刻影响着兴蒙造山带岩石圈的...     

青藏高原东缘地壳上地幔电性结构研究进展

经过数十年的努力,中国学者针对青藏高原东缘地壳上地幔探测,累积完成超过20000 km的大地电磁测深剖面,取得了一系列重要科学数据和认识,为青藏高原东缘构造格局、地壳上地幔电性结构、地震机制和动力学研究奠定了基础.根据青藏高原东缘的主要构造和断裂分布特征,本文重点对龙门山构造带、川滇构造带和三江构造带三个构造带分区进行研究,主要依据大地电磁探测工作成果和壳幔电性结构特征,系统地对青藏高原东缘地壳上地幔电性结构、与扬子西缘接触关系、汶川地震和芦山地震的电性孕震环境及弱物质流通道等几个方面进行了梳理和分析.一是青藏高原东缘地壳表层岩块和物质沿壳内高导层向龙门山造山带仰冲推覆,表现为逆冲推覆特征的薄皮构造;二是高原东部地壳中下部及上地幔顶部向龙门山造山带和上扬子地块西缘岩石圈深部俯冲,呈现刚性的上扬子地块西缘高阻楔形体向西插入柔性青藏块体的楔形构造;三是将汶川地震和芦山地震的震源投影到大地电磁剖面上,发现震源位于剖面下方的高阻块体与低阻体之间靠近高阻体的一侧,龙门山构造带岩石圈表现出高阻、高密度和高速的"三高"特征,这种非均匀电性结构可能构成地震孕育发生条件;四是川滇和三江地区的多条大地电磁剖面探测结果表明,在青藏高原东缘中下地壳存在下地壳流和局部管道流,大地电磁结果对其空间分布形态、位置及大小进行了较好的刻画.根据研究区壳幔电性结构特征的构造解析和综合实例分析,总结了青藏高原东缘六类壳幔电性结构模型,提出了下一步重点研究领域和目标.总之,青藏高原东缘壳幔电性结构的研究对揭示研究区岩石圈结构和构造格局提供了重要依据,对油气及矿产资源远景评价提供了背景资料,对"Y"型多地震区的构造关系和发震机理研究具有重要指导意义.     

扇形边界条件下的龙门山壳幔电性结构特征

沿甘肃碌曲-四川龙门山-重庆合川布设了长周期大地电磁剖面,对龙门山及邻区进行了壳幔电性结构探测,采用更直观合理的扇形边界条件下的反演算法对长周期大地电磁资料进行二维反演.该剖面电性结果揭示了自北西向南东岩石圈深部的若尔盖壳幔高阻块体、松潘壳幔低阻带、龙门山壳幔高阻块体和川中壳幔高阻块体电性结构特征;龙门山逆冲推覆构造带下方的龙门山壳幔高阻体显示为向北西延伸的楔形构造,推断龙门山及松潘-甘孜地块由于受青藏高原东缘和上扬子地块双向挤压,松潘-甘孜地块地壳物质向龙门山逆冲推覆,中下地壳至上地幔向下向南东俯冲,呈现上扬子地块西缘壳幔高阻楔形体插入青藏高原东缘的态势;初步认为上扬子地块西缘深部以松潘壳幔韧性剪切带作为中新生代以来的边界.     

Contribution to the knowledge of the magmatic evolution by the study of the variation of the coefficient of viscosity

The determination of the coefficient of viscosity of eruptive products gives useful elements to the knowledge of possible variations of composition and physical conditions of the magmas.     

Pictorial Series of the Manifestations of the Dynamics of the Earth

Studia Geophysica et Geodaetica -     

Evolution of the Source of the Active Region of the Sun

Geomagnetism and Aeronomy - The magnetic evolution of the active region of the Sun at the level of the photosphere has been studied. Magnetic synoptic maps of the 2007 Kitt Peak Observatory were...     

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.     

The effect of the pollution of the ground layer of the atmosphere on the albedo

Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.     

On the solution of the reverse problem of geodesics with the help of the chord

Резюме Дается вывод формул для вычислений длины хорды, длины и азимутов нормаляных сечений и длины геодезической линии эллипсоида. Применяются тригонометрические функции данных величин и постоянные принятого ргфгргнц-эллипсоида, что является целесообразным при использовании вычислителQjных машин, особенно автоматов. Рабочие формулы рекомендуются в форме (8–11).

---

---

Dedicated to Professor František Fiala on His 85th Birthday

---

Address: Veveři 95, Brno.     

On the determination of the non-linearity of the gravimeter scale

Summary The results of two different methods of the determination of the non-linearity of the gravimeter scale are compared; the economy is discussed as well.     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

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.