Numerical simulation of the crustal stress field by using data of deep geophysical exploration in the northern part of North China

ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｃｒｕｓｔａｌｓｔｒｅｓｆｉｅｌｄｂｙｕｓｉｎｇｄａｔａｏｆｄｅｐｇｅｏｐｈｙｓｉｃａｌｅｘｐｌｏｒａｔｉｏｎｉｎｔｈｅｎｏｒｔｈｅｒｎｐａｒｔｏｆＮｏｒｔｈＣｈｉｎａＣＨＡＮＧＱＵＡＮＬＩＵ（刘...     

The current strain distribution in the North China Basin of eastern China by least-squares collocation

In this paper, the velocities of 154 stations obtained in 2001 and 2003 GPS survey campaigns are applied to formulate a continuous velocity field by the least-squares collocation method. The strain rate field obtained by the least-squares collocation method shows more clear deformation patterns than that of the conventional discrete triangle method. The significant deformation zones obtained are mainly located in three places, to the north of Tangshan, between Tianjing and Shijiazhuang, and to the north of Datong, which agree with the places of the Holocene active deformation zones obtained by geological investigations. The maximum shear strain rate is located at latitude 38.6°N and longitude 116.8°E, with a magnitude of 0.13 ppm/a. The strain rate field obtained can be used for earthquake prediction research in the North China Basin.     

华北渤海湾盆地大震的构造特征

渤海湾盆地是地震活动强烈的地区，记载有7级以上大地震6次。前人对几次大震的具体构造条件已做过大量研究。笔者基于已有的工作，通过对盆地构造发育早第三纪断陷和晚第三纪以来拗陷两个阶段的不同构造后性和区域地震构造格局的分析，以及大震区构造的剖析，综合研究了大震的构造特征。其中最根本的是，盆地发育的断陷阶段因北西—南东向水平拉张作用，于地壳上部形成了由北北东—北东向太行山山前断裂带等大型拆离断裂拉张滑脱控制的伸展断裂系统；拗陷阶段受北东东—南西西向水平挤压，正发育着一套地壳共轭剪切破裂系统。前者与地震关系不大，后者对地震孕育和发生有着重要的控制作用。     

Establishment of the GPS Crustal Deformation Monitoring Network in North China

The Global Positioning System (GPS) is a new method to monitor the crustal deformation for earthquake prediction.Now a four-dimensional crustal deformation monitoring network applying GPS techniques has been established in North China,which is an important seismic monitoring area.Results with high precision have been achieved in the first measurement since a series of advanced methods and measures were adopted,and thus a good foundation has been laid for future researches on earthquake prediction and geodynamics.     

鄂尔多斯盆地深部热结构特征及其对华北克拉通破坏的启示

基于二维稳态热传导方程,利用有限元数值模拟方法,选取东西向横穿鄂尔多斯盆地地质与地球物理解释大剖面进行了深部温度场数值模拟研究,得到了华北克拉通西部的鄂尔多斯盆地下伏岩石圈热结构特征.地幔热流变化范围:21.2~24.5 mW·m-2,体现为东高西低特征.壳幔热流比(Qc/Qm)介于1.51~1.84之间,为"热壳冷幔".与华北东部地幔热流对比表明,西部的鄂尔多斯盆地相对处于稳定的深部动力学环境.在岩石圈热结构研究基础上,对克拉通地震岩石圈与热岩石圈厚度差异进行了对比,研究表明:鄂尔多斯盆地西部地震岩石圈与热岩石圈厚度差异约达140 km,而东部的汾渭地堑,渤海湾盆地二者差异逐渐减小.华北克拉通自西向东,地震岩石圈厚度与热岩石圈厚度差异不断减小,意味着华北克拉通岩石圈下部的软流圈地幔黏性系数自西向东逐渐降低,本文从地热学角度可能印证了太平洋俯冲脱水作用对华北克拉通的影响.     

Crustal structure beneath the central and western North China from receiver function analysis

The North China Craton (NCC) is one of the oldest cratons on earth. Several important tectonic transformations of Mesozoic-Cenozoic tectonic regime led to the destruction of the North China craton. The knowledge of crustal structure can provide important constraints for the formation and evolution of cratons. New maps of sediment thickness, crustal thickness (H) and v_P/v_S (κ) in the central and western NCC were obtained using sequential H-κ stacking. P-wave receiver functions are calculated using teleseismic waveform data recorded by 405 stations from ChinArray project. Benefiting from the densely distribution of temporary seismic stations, our results reveal details of the crustal structure in the study area. The thickness of sedimentary layer in North China ranges from 0–6.4 km, and the thickest sedimentary layer is in Ordos block and its surroundings (about 2.8–6 km); The thickness of sedimentary layer in the Mongolia fold belt and Yinshan orogenic belt is relatively thin (less than 1 km). The crustal thickness of the study area varies between 27–48 km, of which the crust of the North China Plain is about 30–33 km, the central NCC is about 33–40 km, and the Ordos block is 40–48 km thick. The average v_P/v_S ratios in the study area is mostly between 1.66 and 1.90, and that in the Yanshan-Taihang mountain fold belt is between 1.70 and 1.85, and that in the Ordos block is between 1.65 and 1.90, with an average value of 1.77, indicating the absence of a thick basaltic lower crust. The obvious negative correlation between crustal thickness and average v_P/v_S ratio within Ordos and Central Asia orogenic belt may be related to magmatic underplating during the crustal formation. There is no significant correlation between the crustal thickness and the v_P/v_S ratio in the Lüliang-Taihang mountain fold belt, which may be related to the multiple geological processes such as underplating and crustal extension and thinning in this area. The lack of correlation between crust thickness and topography in the central orogenic belt and the North China Basin indicates the topography of these areas are controlled not only by crustal isostatic adjustment but also by the lithospheric mantle processes.     

华北盆地强震孕育的动力学机制研究

强震活跃区通常出现在地壳变形剧烈的地区.但是,华北盆地地壳变形速率小,而破坏性强震却十分频繁.文中利用有限元方法对华北盆地发震断层下方可能存在的速度异常体进行了力学分析,计算结果显示：无论这个异常区是低速异常还是高速异常,在水平构造力的作用下（无论是挤压还是引张力作用）,应力及应变能密度都在空间某些区域积聚.根据本文的工作及前人的研究结果,认为华北盆地强震孕育的物理机制如下：华北克拉通在减薄、破坏过程中,浸入的岩浆或上涌的地幔物质经过漫长的地质时期演化,成为低速异常体.当断裂带下面存在这种低速异常体时,水平向构造力的作用,会使异常体上方应力集中、应变能积聚.当异常体上方的应力水平达到某种限度,介质发生破坏.在破坏过程中,异常体释放的能量,使断裂带上的应变能集聚速度加快;另一方面,异常体的破坏使断层的闭锁面积减小.随着外部构造应力的持续加载,当发震断层上的应力及能量积累达到一定极限,断层产生错动,发生强震.同时,文中结果还显示,断层在从闭锁、积累能量到解锁、发生强震的过程中,其下方的低速异常体起着重要的推进作用.     

区域气候模式RegCM3对华北地区未来气候变化的数值模拟

采用高水平分辨率区域气候模式进行区域未来气候变化预估,对理解全球增暖对区域气候的潜在影响和科学评估区域气候变化有很好的参考价值.这里对国家气候中心使用25 km高水平分辨率区域气候模式RegCM3单向嵌套全球模式MIROC3.2_hires在观测温室气体(1951—2000)和IPCC A1B温室气体排放情景下(2001—2100)进行的共计150年长时间模拟结果,进行华北地区未来气温、降水和极端气候事件变化的分析.模式检验结果表明:模式对当代(1981—2000)气温以及和气温有关的极端气候事件(霜冻日数、生长季长度)的空间分布和数值模拟较好;对降水及和降水有关的极端气候事件(强降水日期、降水强度、五日最大降水量)能够模拟出它们各自的主要空间分布特征,但在模拟数值上存在偏大、偏强的误差.和全球模式驱动场相比,区域模式模拟的气温、降水和极端气候事件有明显的改进.2010—2100年华北地区随时间区域平均气温升高幅度逐渐增大,随之霜冻日数逐渐减少,生长季长度逐渐增多;同时随温室效应的不断加剧,未来降水呈增加的趋势,强降水日期和五日最大降水量逐渐增多、降水强度逐渐增大.从空间分布看,21世纪末期(2081—2100)气温、降水以及有关的极端气候事件变化比21世纪中期(2041—2060)更加明显.     

区域气候模式RegCM3对华北地区未来气候变化的数值模拟

采用高水平分辨率区域气候模式进行区域未来气候变化预估,对理解全球增暖对区域气候的潜在影响和科学评估区域气候变化有很好的参考价值.这里对国家气候中心使用25 km高水平分辨率区域气候模式RegCM3单向嵌套全球模式MIROC3.2_hires在观测温室气体(1951—2000)和IPCC A1B温室气体排放情景下(2001—2100)进行的共计150年长时间模拟结果,进行华北地区未来气温、降水和极端气候事件变化的分析.模式检验结果表明:模式对当代(1981—2000)气温以及和气温有关的极端气候事件(霜冻日数、生长季长度)的空间分布和数值模拟较好;对降水及和降水有关的极端气候事件(强降水日期、降水强度、五日最大降水量)能够模拟出它们各自的主要空间分布特征,但在模拟数值上存在偏大、偏强的误差.和全球模式驱动场相比,区域模式模拟的气温、降水和极端气候事件有明显的改进.2010—2100年华北地区随时间区域平均气温升高幅度逐渐增大,随之霜冻日数逐渐减少,生长季长度逐渐增多;同时随温室效应的不断加剧,未来降水呈增加的趋势,强降水日期和五日最大降水量逐渐增多、降水强度逐渐增大.从空间分布看,21世纪末期(2081—2100)气温、降水以及有关的极端气候事件变化比21世纪中期(2041—2060)更加明显.     

华北地区现今地壳运动动力学初步研究

本文基于GPS、断层形变等观测资料,实现华北地区构造运动有限元数值模拟,研究其现今地壳运动及形变动力学机理.结果表明,鄂尔多斯地块、华南地块、东北亚地块等周边构造块体的相对运动基本决定了华北地区现今表面运动及应力场格局.而另一方面,当考虑区域下部岩石层较快速的“拖动”作用时,表面速度场可以得到更好模拟,并同时形成共轭分布的剪应力梯度带.可见太平洋板块的俯冲作用、印-欧板块的碰撞挤压作用等可能造成岩石层深部、浅部运动差异,从而对研究区现今地壳运动产生深刻影响.此外,地形重力作用、断层分布及区域流变结构非均匀性也对现今地壳运动具有一定影响作用,但处于次要地位.     

Three dimensional tectonic stress field in North China

IntroductionThecrusttectonicstressfieldisoneofthekeyproblemsofgeodynamics,itsmaincontentsareaprobingintospatialdistributionanddynamicevaluationofcrusttectonicstressfield,thentracingtheevolutionprocessofcrusttectonicmovement,learningcontemporarycrusttectoniceventanditsmechanism,thuspredictingthefuturedevelopmentlawoftectonicmovement.Intheendof1970s,aseriesofpaperswerepublishedattheinternationalconferenceonMagnitudeofDeviatoricStressesintheEarth'SCrustandUppermostMantle.Theygavethebriefoutlin…     

Discussion on the dynamic mechanism of Great North China area based on the observed stress data

Introduction The Great North China, located at longitude 106E to 124E and latitude 31N to 42N, in-cludes three secondary active tectonic blocks, Ordos, Yanshan and North China plain (Figure 1). The geological tectonics of these three secondary blocks is much different from each other. As a stable block with high rigidity, the Ordos block is mostly surrounded by down-faulted basins with an inactive interior since Cenozoic, although the fault zones along its boundary are strongly active wi…     

三维黏弹性数值模拟华北盆地地震空间分布与构造应力积累关系

华北盆地为我国板内地震多发区域,历史以来相继发生多次破坏性大地震.前人地震勘探与震源定位结果揭示了华北地震的空间分布特征:横向上,华北地震基本发生在地壳的薄弱地带(Moho面上隆),或者地壳厚度的急剧变化带;纵向上,华北地震在地壳一定深度范围内呈现成层分布特征;主震一般在上地壳底部9~15 km深度范围,余震多发生在大约深5~25 km的上地壳与中地壳范围内,在中地壳下层与下地壳中仅有少量或者鲜见有余震发生.为研究解释华北盆地地震空间分布的以上特征,本文建立了华北盆地岩石圈三维黏弹性有限元模型.震源机制和GPS反映华北盆地处于NNE最大主压应力方向挤压,因此对模型边界施以恒定的位移速率边界条件;数值模拟华北岩石圈各层位在数百年以上长期匀速构造挤压作用下的应力积累特征,分析了华北地震空间分布与构造应力积累速率的关系,探讨了地壳结构与地壳分层流变性质对地壳应力积累的影响.计算结果表明,Moho面的隆起与地壳各层位岩石介质的黏滞系数是华北盆地地震孕育的重要因素.华北盆地在构造挤压的持续作用下,Moho面隆起处产生明显应力集中现象.该区域应力在长时期的积累过程中,在脆性的上地壳与中地壳上层,应力表现近于线性增长趋势,上地壳底部较其它深度有最大的应力增长率,其主震可以在应力积累至岩石破裂强度时发生;在脆、韧性转换的中地壳下层,应力增长速率次之,华北地震的大部分余震可能在该层位为主震所触发;而在柔性的下地壳应力增长近于指数形式,稳定状态之后其应力增长速率近于零,而鲜有地震发生.地壳各层位的应力增长率差异与地震成层分布的现象揭示了华北地壳的分层流变性质:脆性(上地壳)-较弱脆性(中地壳上层)-较弱韧性(中地壳下层)-较强韧性(下地壳)-韧性(岩石圈上地幔)的分层流变结构.     

大同地震与华北北部地壳异态垂直运动

本文根据笔者提取异态信息的新数据处理方法，分析了华北北部地区近十年来地壳异态垂直运动，结果显示出以大同地震为界，震前的整个华北北部地区均处于应变的积累阶段，而且震源区具有非常清晰的四象限分布图像，呈一种典型的水平应力作用下的点源理论垂直形变场图象，震后整个华北北部地区则出现了截然相反的异态垂直运动，这又清楚地表明了具有弹性回跳的理论特征。迄今为止，这是我们唯一获得中强地震前后地壳异态垂直运动在空间     

Evolution of 3D tectonic stress field and fault movement in North China

IntroductionThe research on stress field interior the Earth is an either relative old or young branch of geoscience. Early in the 20th century, in order to learn the state of inner stress field, a few scholars had used simple ideal model to analyze big scale tectonic conformations and activities according to tide, gravity, rotation and inner thermal materials. So far, lots of problems on the state of the Earth stress field have not been solved completely yet. The origin of stress field, the d…     

上海地区地壳应变能力数值模拟

上海市综合深井项目第2阶段将建设9—10个综合深井观测站。利用蠕变方程,通过ANSYS有限元软件进行数值分析,模拟上海市9个观测点在当前地壳运动水平下的形变信息,为综合深井第2阶段建设提供理论基础。     

Basin and crustal velocity structure models for the simulation of strong ground motions in the Kinki area,Japan

We constructed a prototype of the basin and crustal structure model for the Kinki area, southwest of Japan, for the simulation of strong ground motions of hypothetical crustal and subduction earthquakes. We collected results of the deep seismic velocity profiles obtained by the reflection experiments and seismic imaging results, which were conducted in the Kinki area. The obtained profiles give underground velocity structures of the crust, from the surface to the subducting slab. We also gather the basin velocity structure information of the Osaka, Kyoto, Nara, and Ohmi basins. To examine the applicability of the constructed velocity structure model to the ground motion simulation, we simulated waveforms of an intermediate size event occurred near the source area of the hypothetical subduction earthquakes. Simulated ground motions using the basin and crustal velocity structure model are fairly well reproducing the observations at most of stations, and the constructed basin and crustal velocity structure model is applicable for the long-period ground motion simulations.     

Rheology of the lower crust beneath the northern part of North China: Inferences from lower crustal xenoliths from Hannuoba basalts, Hebei Province, China

Lower crustal xenoliths brought up rapidly by basaltic magma onto the earth surface may provide di-rect information on the lower crust. The main purpose of this research is to gain an insight into the rheology of the lower crust through the detailed study of lower crustal xenoliths collected from the Hannuoba basalt, North China. The lower crustal xenoliths in this area consist mainly of two pyroxene granulite, garnet granulite, and light-colored granulite, with a few exception of felsic granulite. The equilibration temperature and pressure of these xenoliths are estimated by using geothermometers and geobarometers suitable for lower crustal xenoliths. The obtained results show that the equilibration temperature of these xenoliths is within the range of 785―900℃, and the equilibrium pressure is within the range of 0.8―1.2 GPa, corresponding to a depth range of 28―42 km. These results have been used to modify the previously constructed lower crust-upper mantle geotherm for the studied area. The dif-ferential stress during the deformation process of the lower crustal xenoliths is estimated by using recrystallized grain-size paleo-piezometer to be in the range of 14―20 MPa. Comparing the available steady state flow laws for lower crustal rocks, it is confirmed that the flow law proposed by Wilks et al. in 1990 is applicable to the lower crustal xenoliths studied in this paper. The strain rate of the lower crust estimated by using this flow law is within the range of 10-13―10-11 s-1, higher than the strain rate of the upper mantle estimated previously for the studied area (10-17―10-13 s-1); the equivalent viscosity is estimated to be within the range of 1017―1019Pa·s, lower than that of the upper mantle (1019―1021 Pa·s). The constructed rheological profiles of the lower crust indicate that the differential stress shows no significant linear relation with depth, while the strain rate increases with depth and equivalent vis-cosity decrease with depth. The results support the viewpoint of weak lower continental crust.     

华北克拉通岩石圈三维密度结构

显生宙以来华北克拉通岩石圈遭到破坏,这一现象的科学问题已受到世界地学家广泛关注.本文首先将地震层析成像反演得到的P波速度扰动转化为密度扰动,以此作为初始密度模型,然后利用布格重力异常反演得到了华北克拉通岩石圈高分辨三维密度结构.为了避开大型稀疏矩阵求逆计算,提高计算效率,我们将代数重构技术用于密度反演解算.反演结果表明:华北克拉通岩石圈密度在横向和纵向上均存在明显的不均匀性,密度分布形态与地表构造格局有很好的相关性;研究区地壳整体表现为低密度异常,地壳以下岩石圈部分则以高密度异常为主;鄂尔多斯块体地壳范围内以低密度异常为主,80~120 km深度上为呈南北两端集中分布的高密度异常,并分别与秦岭造山带和阴山造山带的高密度异常分布相连,这暗示了鄂尔多斯块体可能受到了来自其南北两端造山带深部动力学过程的影响;80~120 km深度上,华北克拉通东部地区呈现出显著的南北向非均匀的高密度异常,这表明遭到破坏后该地区上地幔物质分布具有强烈的南北向非均匀性.     

华北地区的前兆性地震系列

通过对1966年以来的华北地区地震活动期强震前的前性兆系列研究,得出华北地区在现今活动期内有统一的地壳应力场作用的结论。