岩石圈中水的临界奇异性——大规模成矿的根源

水的临界奇异性是岩石圈中许多重大地质作用的关键和枢纽,有独特的地质意义、超乎常识。水在临界点发生二级相变,许多物理化学性质奇异性突变,特别是溶解行为和热压的突变,对成矿作用影响极大。     

水的临界奇异性及其对热液铀成矿作用的意义

水是地质流体的主体,随着温压条件的变化会发生相变,许多物理化学性质也随之改变.特别是水具有临界奇异性,这对热液成矿作用有巨大的影响.本文介绍了水的相变种类及特征,研究了岩石圈中水可能存在的相变,以等压热容为例研究了水的临界奇异性,指出:(1)当温压同时趋近于临界常数时,水的临界奇异性最明显;(2)在岩石圈中温压受地温线控制,只有在热液作用、中酸性岩浆活动及断裂作用同时出现时,才可能使上地壳中某一区域同时达到温压的临界值;(3)对于热液成矿作用而言,水的临界奇异性的意义至少包括溶解度发生急剧变化和产生瞬时高压,这些变化可导致成矿物质沉淀和断裂的产生.     

热液矿床水相变控矿理论初探

在综合水物理化学性质、热液矿床、溶解度实验研究等资料的基础上,曾提出热液矿床水相变控矿理论框架.在此进一步完善该理论:①当水发生相变时可引起热容、水中矿质溶解度、压缩系数、膨胀系数等物理化学性质突变,特别是在临界点处二级相变时还出现"临界奇异性"、部分参数趋于无穷大,与水在非相变区物理化学性质的小幅度渐变区剐明显,必将...     

岩石圈中的降压作用及其相变过程

岩石圈中许多地质过程,都曾经历过降压作用。岩石圈中的压力是一个极为重要的参数,压力的降低将改变地质作用过程的其他物理化学性质,包括导致相变。通过岩石圈中降压过程性质和机制的研究,重点探讨降压所引起的相变及其地质意义,试图从新的研究角度提出一些认识和思路。(1)降压是岩石圈中许多地质作用发生的必要条件,降压的性质不同对应不同的地质现象。(2)岩石圈物质的相变线(含两端点)及邻区的物理化学行为突变有重要研究价值。降压可导致岩石圈中的水发生一、二级相变,影响着热液成矿和油气成藏,是"构造降压-水相变耦合成矿(藏)"新模式的理论核心。降压导致岩石圈中固体熔融相变,是岩浆形成的必要条件,降压性质差异决定着岩浆岩的种类和结构。降压过程中,矿物组合随之调整,因各种矿物的相变线不同,矿物变形行为表现各异,地压梯度明显降低,影响超高压变质作用和变形、岩石圈减薄等的解释。(3)降压通常是岩石圈局部升温和聚集水的前提。岩石圈中若没有降压,就不会出现明显升温和聚集水的现象。     

长江中下游成矿带岩石圈结构与成矿动力学模型——深部探测(SinoProbe)综述

岩石圈结构和深部过程对理解成矿带和大型矿集区的形成十分重要。岩石圈尺度的地球动力学过程将在地壳中留下各种结构的或物质的"痕迹",这些"痕迹"可以通过地球物理的手段去探测。为深入理解长江中下游成矿带形成的深部动力学过程,作者在国家深部探测专项(SinoProbe)和国家自然科学基金重点项目支持下,在长江中下游成矿带开展了综合地球物理探测。方法包括宽频地震、深地震反射、广角反射/折射和大地电磁测深。数据处理和反演结果取得一系列新发现:(1)成矿带上地幔顶部存在低速体,在中心深度300km处有一向SW倾斜的高速体;(2)S波接收函数证实成矿带岩石圈较薄,只有50~70km;横波分裂结果显示,成矿带上地幔各向异性方向和强度与邻区有较大区别,显示平行成矿带(NE-SW向)的上地幔变形和流动;(3)深反射地震揭示成矿带上地壳曾发生强烈挤压变形,以紧闭褶皱、逆冲和推覆为特征;在宁芜火山岩盆地、长江断裂带和郯庐断裂之下出现"鳄鱼嘴"构造,指示上下地壳在挤压变形过程中解耦;深反射地震证实发生过陆内俯冲和叠瓦,并认为是岩石圈增厚和拆沉的主导机制;(4)广角反射和大地电磁反演给出了跨成矿带地壳剖面的速度和电性结构,速度和电阻率分布总体上与构造单元相吻合。本文分析和解释了这些发现的地质意义,并结合近年在长江中下游地区的地球化学研究进展,提出了成矿带地球动力学模型。该模型认为:中、晚侏罗世陆内俯冲、岩石圈拆沉、幔源岩浆底侵和MASH过程造就了长江中下游世界级成矿带的形成。     

大陆岩石圈有效弹性厚度与有关地球物理参数的关系--以泉州-黑水地学断面为例

大陆岩石圈有效弹性厚度(Te)是反映岩石圈综合强度的参数,它反映了岩石圈的整体特征。分析岩石圈有效样性厚度与反映深部地质特征的有关地球物理参数之间的关系,对研究控制Te的因素、各因素之间的关系以及探索大陆构造与大陆动力学等具有重要意义。泉州一黑水地学断面Te与地壳厚度、热岩石圈厚度、均衡重力异常、磁性构造层底面深度、上地幔低速层顶界面深度、上地幔低阻层顶面深度之间的关系研究表明：Te与大地热流关系密切的“热”地球物理参数磁性构造层底面深度、热岩石圈厚度相关性好;与地壳厚度有一定的相关性;上地幔低速层顶界面深度和上地幔低阻层顶面深度与大陆岩石圈Te相关性均较差。     

青藏高原莫霍面的研究进展

本文首先简要回顾了莫霍面的发现，介绍其基本性质，然后对青藏高原莫霍面研究的重要进展进行了评述。在区域尺度上，被动源地震(天然地震)方法研究结果勾勒出青藏高原地壳及岩石圈底部的深部构造轮廓。然而受分辨率的限制，天然地震结果给出的地壳及上地幔结构的细节不足。近年来已经用分辨率达到几千米甚至百米级的主动源地震(包括宽角反射与折射地震和深反射地震)方法，揭示出青藏高原地壳及上地幔的精细结构。本文对近30年来深地震探测获得的青藏高原各个地块的莫霍面深度、壳幔结构和上地幔盖层速度等基本数据进行了较系统的分析，并对青藏高原莫霍面研究存在的有关问题进行了讨论。     

长江中下游南京（宁）-芜湖（芜）段深部壳幔电性结构——宽频大地电磁测深研究

为了较全面、客观地认识南京(宁)-芜胡(芜)地区岩石圈深部结构,探讨岩石圈热结构和壳、幔物质状态等重要科学问题,为矿集区成矿作用和成矿规律的研究提供依据,我们完成了6条宽频大地电磁测深剖面。通过分析各剖面电性成像结果,讨论了研究区地壳-上地幔导电性的"拟三维"结构,发现测区普遍存在连续的下地壳-上地幔低阻层,仅在巢湖冲褶体断开。结合区域热流结果,我们认为下地壳-上地幔有较好的导电性,可能存在局部"熔融体"或"含水剪切断裂",其物质状态很可能是热的、软弱的。此外,我们通过物质状态和电性界面推断了上地幔隆起的位置和长江深断裂带的分布范围,认为长江断裂带不但存在,而且由多条北深南浅的断裂组成,较为复杂。     

长江中下游及其邻区深部地球物理背景与含金夕卡岩矿床的分布

长江中下游地区是我国重要的铁铜矿产基地。近十年来,我国所发现的含金夕卡岩矿床和铜伴生金夕卡岩矿床也主要产在长江中下游地区。作者借助12条大地电磁测深剖面、5条地震剖面、层析成像速度结构资料、重磁场等区域的和深部的地球物理资料进行综合对比研究,给出长江中下游及其邻区的三维深部构造格架及其与含金夕卡岩矿床和铜伴生金夕卡岩矿床的分布关系。作者认为,上地幔隆起带(岩石圈地幔减薄带)、上地幔异常区(相对低速区)、壳内高导层隆起带、深断裂(岩石圈剪切带)、地壳上地幔不均匀性块体的边缘、重力高反映的基底隆起区、跳跃磁场反映的岩浆岩带和构造交汇处等诸多因素的共同作用控制着含金夕卡岩矿床和铜伴生金夕卡岩矿床的分布。     

福建省上地幔结构与岩石矿物的基本特征

综合研究了福建省上地幔的矿物、岩石、地球化学、地球物理（重、震、电）等有关资料，全面论证上地幔B′与B″层的三维空间特征、可疑地幔岩性、古代与现代地幔岩岩石类型与矿物成分。继而从地幔岩矿物组合，化学成分、微观特征、包裹体、稳定同位素等特点，来探讨第三纪上地幔的成岩环境与动力学，并指出东南沿海岩石圈上地幔减薄斜坡带及其叠生的深切断裂是燕山期双系列岩浆岩形成的深部地质背景。     

青藏高原东南缘活动断层相互作用、应力触发与差别响应

岩石圈、地壳由众多板块、地块及层圈组合而成,是开放性复杂巨系统,活断层在地壳中犹如神经网络也是复杂的开放体系。因此活断层之间存在复杂的相互作用,例如地震断层破裂产生应力扰动,可能触发其他断层破裂,不仅在近处,也会在远处发生。不同的活断层构造产状、活动方式及应变阶段不同,对同一次触发响应不同。触发与差别响应二者的叠加效应在地震活动性上有重要表现。表现之一是区域大震后,余震区外地震活动显著增强处发生继发性大震的概率最大,也即率先出现"远余震、诱发前震、响应震"的地方地震危险性增大,对预测未来地震位置有效。对本区1950—2013年地震统计表明,预测成功率W=1-漏报率-虚报率=80%。同时,对当前地震危险区作了预测。     

青藏高原板内地震震源深度分布规律及其成因

青藏高原板内地震以浅源地震为主, 下地壳基本上没有地震, 地震震源多集中在15~40 km的深度范围, 主要在中地壳内, 呈似层状弥散分布.其中30~33 km深度是一个优势层, 与壳内分层有关.总体上青藏高原南、北部的震源面略呈相向倾斜特征.70~100 km深度区间出现了比较集中的震级较小的地震, 可能与壳幔过渡带的拆离作用有关.高原内部的正断层系与板内地震密切相关, 是板内浅源地震的主控构造.总之, 青藏高原地震震源沿着活动的上地壳脆性层与软弱层之间的脆-韧性过渡带分布.这些板内地震活动属于大陆动力学过程, 与板块碰撞和板块俯冲无关.初步认为青藏高原浅层到深层多震层的成因分别是韧性基底与脆性盖层、韧性下地壳与脆性上地壳、韧性下地壳与脆性上地幔的韧-脆性转换、拆离和解耦的产物.      

地震动力学:震源分布与岩石圈流变特性(英文)

地震震源分布强烈依赖于构造环境的温度和压力条件。震源机制可使大陆地壳脆性-韧性转换带(下部稳定性过渡边界)的温压条件复杂化。该过渡边界伴随速度弱化作用(有震活动)向速度强化作用(无震活动)的转化。在岩石圈流变和壳幔动力学的基础上研究了与板块边界有关的地震活动,包括板缘地震和俯冲板片地震。俯冲带板缘地震的深度分布受约于脆性摩擦动力机制,而摩擦剪切机制不能满意地解释深震活动,包括俯冲板片地震。这是因为深震震源机制可能与高压、高温条件下的固一固相变有关,而用脆性破裂或摩擦作用来解释就不近合理。以理论与实验研究为依据,本文对与震源物理和震源分布有关的岩石圈流变特性进行了较为深入的论述。     

A tectonophysical model of the Baikal seismic zone: testing and implications for medium-term earthquake prediction

The first tectonophysical model of the Baikal seismic zone represents a separate complex region of the lithosphere. It has a pinnate structure with a backbone belt of current deformation, which is a concentrator of largest earthquakes, and branching, repeatedly reactivated large and small faults. In its vertical section, the seismic zone is tree-like, the stem and the branches being faults of different size ranks which can generate earthquakes when reactivated. The real-time short-period fault motions and the respective seismicity occurring at a certain time and in certain places are triggered by strain waves, which disturb the metastable state of the faulted lithosphere subject to regional stress. The modeling work includes developing general requirements for tectonophysical models of continental rifts and special methods for identifying the faults that become active within short historic time spans, as well as techniques for locating potential events in space and time in specific active faults. The methods and model testing for medium-term earthquake prediction are described by the example of the well-documented Baikal seismic zone, which is the most active part of the Baikal rift system. The tectonophysical model for the Baikal zone is statistically supported by field data, and this allows estimating the velocities and periods of strain waves for different zone segments and faults, with implications for nearest-future earthquake prediction.     

Seismic coupling and uncoupling at subduction zones

Seismic coupling has been used as a qualitative measure of the “interaction” between the two plates at subduction zones. Kanamori (1971) introduced seismic coupling after noting that the characteristic size of earthquakes varies systematically for the northern Pacific subduction zones. A quantitative global comparison of many subduction zones reveals a strong correlation of earthquake size with two other variables: age of the subducting lithosphere and convergence rate. The largest earthquakes occur in zones with young lithosphere and fast convergence rates, while zones with old lithosphere and slow rates are relatively aseismic for large earthquakes. Results from a study of the rupture process of three great earthquakes indicate that maximum earthquake size is directly related to the asperity distribution on the fault plane (asperities are strong regions that resist the motion between the two plates). The zones with the largest earthquakes have very large asperities, while the zones with smaller earthquakes have small scattered asperities. This observation can be translated into a simple model of seismic coupling, where the horizontal compressive stress between the two plates is proportional to the ratio of the summed asperity area to the total area of the contact surface. While the variation in asperity size is used to establish a connection between earthquake size and tectonic stress, it also implies that plate age and rate affect the asperity distribution. Plate age and rate can control asperity distribution directly by use of the horizontal compressive stress associated with the “preferred trajectory” (i.e. the vertical and horizontal velocities of subducting slabs are determined by the plate age and convergence velocity). Indirect influences are many, including oceanic plate topography and the amount of subducted sediments.All subduction zones are apparently uncoupled below a depth of about 40 km, and we propose that the basalt to eclogite phase change in the down-going oceanic crust may be largely responsible. This phase change should start at a depth of 30–35 km, and could at least partially uncouple the plates by superplastic deformation throughout the oceanic crust during the phase change.     

Seismogenic structures along continental convergent zones: from oblique subduction to mature collision

We summarize seismogenic structures in four regions of active convergence, each at a different stage of the collision process, with particular emphases on unusual, deep-seated seismogenic zones that were recently discovered. Along the eastern Hellenic arc near Crete, an additional seismogenic zone seems to occur below the seismogenic portion of the interplate thrust zone—a configuration found in several other oblique subduction zones that terminate laterally against collision belts. The unusual earthquakes show lateral compression, probably reflecting convergence between the subducting lithosphere's flank and the collision zone nearby. Along oblique zones of recent collision, the equivalence between space and time reveals the transition from subduction to full collision. In particular, intense seismicity beneath western Taiwan indicates that along the incipient zone of arc–continent collision, major earthquakes occur along high-angle reverse faults that reach deep into the crust or even the uppermost mantle. The seismogenic structures are likely to be reactivated normal faults on the passive continental margin of southeastern China. Since high-angle faults are ineffective in accommodating horizontal motion, it is not surprising that in the developed portion of the central Taiwan orogen (<5 Ma), seismogenic faulting occurs mainly along moderate-dipping (20–30°) thrusts. This is probably the only well-documented case of concurrent earthquake faulting on two major thrust faults, with the second seismogenic zone reaching down to depths of 30 km. Furthermore, the dual thrusts are out-of-sequence, being active in the hinterland of the deformation front. Along the mature Himalayan collision zone, where collision initiated about 50 Ma ago, current data are insufficient to distinguish whether most earthquakes occurred along multiple, out-of-sequence thrusts or along a major ramp thrust. Intriguingly, a very active seismic zone, including a large (M_w=6.7) earthquake in 1988, occurs at depths near 50 km beneath the foreland. Such a configuration may indicate the onset of a crustal nappe, involving the entire cratonic crust. In all cases of collision discussed here, the basal decollement, a key feature in the critical taper model of mountain building, appears to be aseismic. It seems that right at the onset of collision, earthquakes reflect reactivation of high-angle faults. For mature collision belts, earthquake faulting on moderate-dipping thrust accommodates a significant portion of convergence—a process involving the bulk of crust and possibly the uppermost mantle.     

Magnetotelluric Constraints on the Occurrence of Lower Crustal Earthquakes in the Intra-plate Setting of Central Indian Tectonic Zone

Lower crustal earthquake occurrence in the Central Indian Tectonic Zone(CITZ) of the Indian sub-continent was investigated using magnetotelluric(MT) data. MT models across the CITZ, including the new resistivity model across the 1938 Satpura lower crustal earthquake epicenter, show low resistive(80 ?m) mid-lower crust and infer small volume(1 vol%) of aqueous fluids existing in most part of lower crust. This in conjunction with xenoliths and other geophysical data supports a predominant brittle/semi-brittle lower crustal rheology. However, the local deep crustal zones with higher fluid content of 2.2%–6.5% which have been mapped imply high pore pressure conditions. The observation above and the significant strain rate in the region provide favorable conditions(strong/moderate rock strength, moderate temperature, high pore pressure and high strain rate) for brittle failure in the lower crust. It can be inferred that the fluid-rich pockets in the mid-lower crust might have catalyzed earthquake generation by acting as the source of local stress(fluid pressure), which together with the regional stress produced critical seismogenic stress conditions. Alternatively, fluids reduce the shear strength of the rocks to favor tectonic stress concentration that can be transferred to seismogenic faults to trigger earthquakes.     

滇西南2014年景谷中-强震群的地质构造成因——茶房-普文断裂带贯通过程的构造响应

2014年10—12月期间,云南景谷接连发生了Ms6.6、Ms5.8、Ms5.9三次中-强地震。为确定地震的地质构造成因,在地表调查的基础上,综合该区的地质构造情况、烈度与余震分布、震源机制解等资料,确定此次震群活动的宏观震中位于永平盆地东南侧山地,发震断层为地质与地貌表现不显著的NW向右旋走滑断层。此次震群活动及余震迁移过程指示,由于断层斜接部位岩桥的临时阻碍,Ms6.6地震破裂在向南东扩展过程中发生短暂停滞,突破障碍后进一步引发了Ms5.8和Ms5.9地震,这符合震源破裂沿NW向发震断裂分段破裂的行为。区域活动断裂的遥感解译结果发现,发震断层位置恰好处于NW向右旋走滑的茶房断裂与普文断裂之间,区域上属于该断裂带的不连贯部位,指示此次中-强震群活动应该是茶房-普文断裂带贯通过程的构造活动表现。结合思茅地块的历史地震资料发现,思茅地块地震活动多以小于等于6.8级为主,发震构造多为NW向断裂。指示在现今构造应力场作用下,该区NW向断裂的活动性相对NE向断裂更加显著,属于该区主要控震构造,应在今后的地震地质工作中给予更多关注。     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.