昆仑山Ms 8.1级地震震后变形场数值模拟与成因机理探讨

2001年11月14日,在青海和新疆交界处发生了昆仑山M_s8.1级强烈地震,GPS后观测显示,此次地震震后形变不仅在断裂南北两侧存在很大的差异,而且在短时间调整后断裂南北两侧表现为同向运动.本文以观测的地震形变为约束,通过有限元数值模拟分析昆仑山地震震后形变的物理机制.建立有关的有限元虚功方程,通过有限元数值方法模拟震后形变,从理论上分析介质的非均匀性、黏滞性松弛、流体调整对震后形变的影响.采用网格搜索确定昆仑断裂南北两侧下地壳的黏滞系数分别为5.0×10¹⁷Pa·s, 9.0×10¹⁸Pa·s左右,正是这十余倍的差异引起了断裂两侧震后形变的非对称性和同向运动,这一差异既是长期地质作用的结果,又是现代地球动力学环境的决定因素之一.通过数值模拟定性讨论了断裂北侧地表形变在震后短期内的调整,对于靠近断裂附近的测点可能是黏弹性松弛和孔隙流体调整共同作用的结果,所以在分析短期震后形变时综合考虑黏弹松弛和孔隙流体调整是很有必要的.     

四川龙门山断裂带高精度地壳/岩石圈黏度结构及其动力学意义

岩石圈黏度是大陆动力学研究中一个重要参数,但是岩石圈黏度,尤其是横向小尺度(<100 km)黏度结构的确定是一个挑战.本文根据电阻率和黏度与它们控制因素的相似关系,直接把一条跨过青藏高原东缘和四川龙门山断裂带的大地电磁(MT)探测的电阻率剖面转换成黏度结构作为输入,在GPS速度和地表地形数据的约束下,利用地球动力学数值模拟获得了该剖面的二维地壳/岩石圈黏度结构.本文推断的黏度与前人获得的区域尺度的黏度值一致,但揭示出了更多的细节.本文的黏度结构揭示出研究区域内的地壳/岩石圈黏度存在较大的空间变化范围(约5量级),黏度值分布在1.48×10^17~8.44×10^22 Pa·s之间;龙门山断裂带下的黏度存在强烈的小尺度横向变化,其中、下地壳的黏度分别为1.99×10^18~8.21×10^20 Pa·s(平均1.17×10^20 Pa·s)和4.09×10^19~7.08×10^20 Pa·s(平均1.77×10^20 Pa·s).基于该黏度结构的地球动力学模型表明驱动青藏高原中-下地壳物质流动的可能是热-化学浮力,以及上地壳和中-下地壳可能处于解耦状态.本文获得的黏度结构可以为龙门山断裂带地震成因和机制、岩石圈小尺度变形和构造应力状态的深入研究提供重要的帮助.     

基于2001年MW7.8可可西里地震震后形变模拟研究藏北地区岩石圈流变学结构

青藏高原岩石圈的流变学结构和形变机制是地学界长期争论的重大科学问题.2001年发生在东昆仑断裂带的M_W7.8可可西里地震造成青藏高原北部地区岩石圈构造应力场的很大改变,引起下地壳与上地幔的快速弛豫形变,从而为研究这一问题提供了难得的机会.本研究采用该区域的GPS震后观测,反演这一地区岩石圈的流变学参数并探讨其形变机制.反演所采用的数据来自45个GPS观测点,其中包括一个中国地壳运动观测网络的基准站,数据最长时间跨度达6.4年.大地震震后形变场主要来源于地壳、上地幔的黏弹性松弛与断层面上的震后余滑,因此本研究同时反演介质的黏滞系数和断层的震后余滑.考虑到东昆仑断层南侧的巴颜喀拉-羌塘地区与北侧的柴达木盆地地区具有明显不同的地壳结构,断层南北两侧采用不同的Burgers体流变学结构,其下地壳-上地幔的短期和长期黏滞系数采用网格搜索法获得;断层震后余滑反演则同时施加近似正比于库仑应力的约束.最终结果显示：东昆仑断层北侧柴达木盆地地区下地壳-上地幔短期和长期黏滞系数分别为5×10¹⁸ Pa·s和1.5×10²⁰ Pa·s;东昆仑断层南侧巴颜喀拉-羌塘地区下地壳-上地幔短期和长期黏滞系数分别为1.5×10¹⁸ Pa·s和1.5×10¹⁹ Pa·s.这一结果表明：巴颜喀拉-羌塘地区下地壳-上地幔黏滞系数显著低于柴达木盆地,意味着巴颜喀拉-羌塘地区下地壳可能存在部分熔融,其地壳形变模式更趋近于连续形变,而柴达木盆地形变模式更趋近于块体运动.研究区下地壳长期黏滞系数比下地壳流模型所主张的黏滞系数高2~3个数量级,表明下地壳流在本地区可能不存在.     

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.     

Cenozoic sedimentary records and geochronological constraints of differential uplift of the Qinghai-Tibet Plateau

Geological mapping data (1:250000) in the Qinghai-Tibet Plateau and its adjacent regions reveal the sediment sequences, distribution and tectonic evolution of the 92 Tertiary remnant basins. Southern Tibet and the Yecheng area in Xinjiang, located at southern and northwestern margins of the Qinghai-Tibet Plateau, respectively, were parts of the Neo-Tethys remnant sea in the Paleogene. In southern Tibet, both the subabyssal and abyssal sequences occur at the Gyangze, Saga, Guoyala, and Sangmai areas. The deep-water facies successions outcrop in the west, whereas the shallow-water facies sequences in the east, indicating the east to the west retreat of the Neo-Tethys Ocean. The retreat of the Neo-Tethys Ocean in the east was contributed to the earlier tectonic uplift of the eastern Qinghai-Tibet Plateau. The uplift process of the Plateau from the Late Cretaceous to Pliocene is described as follows: During the Late Cretaceous, tectonic uplift of the Qinghai-Tibet Plateau occurred in the northeastern part and the configuration of the Qinghai-Tibet Plateau was characterized by rise in the northeast and depression in the west. In the Paleocene-Eocene interval, the Tengchong-Baingoin and Kuyake-Golmud areas experienced local tectonic uplifting, the West Kunlun uplift zone broadened easterly, the Qilian uplift zone broadened southerly, and the Songpan-Garzê uplift zone shrank easterly. The Oligocene configuration of the Qinghai-Tibet Plateau was characterized by mountain chains rising along its margins and sedimentary basins in the central part because of tectonic uplifts of the Gangdisê and the Himalaya blocks. Meanwhile, the Kunlun-Altyn-Qilian uplift zones have also broadened southerly and northerly. In contrast, the great uplift zones of the Gangdisê, the Himalaya, the Karakorum, and the Kunlun blocks characterize the paleogeographic contours of the Qinghai-Tibet Plateau during the Miocene-Pliocene. Additionally, the thermochronological data on tectonic uplift events in southern Tibet, West Kunlun Mountains, Altyn Tagh, eastern Tibet, and western Sichuan all suggest that the most intense deformation occurred at 13-8 Ma and since 5 Ma, respectively, corresponding to two great uplift periods in Neogene. As a result, turnover of paleogeographic configuration of the Qinghai-Tibet Plateau occurred during the Neogene, experiencing a change from high contours in the east in the pre-Oligocene to high contours in the west at the end-Pliocene. The uplift of the Qinghai-Tibet Plateau during the Cenozoic was episodic, and the uplifts of various blocks within the Plateau were spatially and chronologically different.     

震后GPS观测数据揭示的日本M_W9.0地震周边地区地幔黏滞性结构垂向变化

本文利用大范围的震后GPS数据和黏弹性球形地球位错理论,定量研究了日本M_W9.0地震周边地区地幔黏滞性结构的垂向变化.首先结合陆地和海底的GPS观测数据,以及基于球形地球位错理论格林函数和贝叶斯反演方法,反演了该地震的同震滑动分布,发现其最大错动量高达59m.然后在均一地幔黏滞性结构的假设前提下,确定了震源周边地区地幔黏滞因子的最优解,发现依据该地幔黏滞因子获得的理论远场震后位移和GPS观测结果之间的均方根误差高达0.81cm,不能解释远场观测结果.为解决上述问题,本文对震中周边地区地幔黏滞性结构沿垂向方向进行分层,建立了一个随深度变化的地幔黏滞性构造模型,然后综合利用远近场的GPS数据对该地区地幔黏滞因子进行反演研究,结果表明,震源周边地区岩石圈弹性层厚度最优解为40km,40～220km深度的地幔黏滞因子最优解为6×10~(18)Pa·s,220～670km深度之间的地幔黏滞因子最优解为1.5×10~(19)Pa·s.上述地幔黏滞性构造使远场的均方根误差降为0.12cm,仅为利用均一地幔黏滞性构造所得均方根误差值的15%,大大提高了远场模拟结果的准确性.最后,观测值和模拟值之间的均方根误差分析表明,近场震后形变数据主要约束浅层的地幔黏滞性结构,而远场震后形变数据主要约束深部的地幔黏滞性结构.     

InSAR约束下的2008年汶川地震同震和震后形变分析

2008年5月12日,青藏高原东缘的龙门山断裂带上发生了M_w7.9级汶川地震.本文通过分析覆盖汶川地震震中区域的ALOS/PALSAR像对的方位向偏移量来选择无明显电离层扰动影响的像对进行干涉处理,获取了高精度、连续的InSAR地表形变场.在此基础上,结合高精度GPS同震形变数据,采用同震、黏弹性松弛震后形变联合反演模型同时确定了汶川地震的同震滑动分布和龙门山地区的流变结构参数.研究结果表明,汶川地震是一个断层破裂非常复杂的地震事件,其中,北川段、岳家山段、虹口段和汉旺段的滑动以逆冲为主,而青川段以右旋走滑为主.滑动主要发生在10 km深度以上的区域,最大滑动量位于虹口段的东北端,达10.7 m.地震释放的总能量为9.28×10²⁰ N·m（M_w7.91）,与地震学的结果一致.联合反演模型确定的龙门山地区中下地壳的黏性系数为2×10¹⁸ Pa·s,为青藏高原东部地区的黏性系数提供了一个可靠的下限值.如果有更长时间的震后形变观测时间序列,将为该区域提供更为可靠的流变结构.     

ESTIMATING THE LOWER CRUSTAL VISCOSITY OF THE WESTERN QINLING-SONGPAN TECTONIC NODE AND ITS ADJACENT AREAS BY USING LANDFORM MORPHOLOGY

The western Qinling-Songpan tectonic node is located at the intersection of three major tectonic units of Tibetan plateau, the South China Block and the Ordos Block, and is at the forefront of the northeastern margin of Tibetan plateau. It has unique geological and dynamic characteristics from the surface to the deep underground. Based on the model for ductile flow in the lower crust, the geomorphological form is used to estimate the viscosity of the lower crust, and how the rheological process of the deep lithosphere acts on the upper crust deformation and structural geomorphology. And combined with GPS velocity field data, the current crustal deformation is analyzed to further study the regional dispersive deformation process. The results show that the viscosity of the north and northeast of the Zoige-Hongyuan Basin is smaller than that of the east and southeast. Therefore, the lower crust flow has a tendency of flowing to the northeastern low viscosity zone. We believe that when the lower crust flows from the central plain of the Qinghai-Tibet Plateau to the rigid Sichuan Basin with a higher viscosity of the lower crust, it cannot flow into the basin, and part of the lower crust flow accumulate here, causing the upper crust to rise, and the uplifting led to the formation of the Longmen Mountains and a series of NNE-striking faults as well. When the lower crust flows to the northeast direction with a low viscosity, the brittle upper crust is driven together. Because of the remote effects from the Ordos Basin and the Longxi Basin, the mountains in this region are built slowly and the stepped arc-shaped topography of the current 3 000-meter contour line and the 2 000-meter contour line are developed. At the same time, a series of NWW-trending left-lateral strike-slip faults are developed. This explains the seismogenic tectonic model of the western Qinling-Songpan tectonic node as from NWW-trending left-lateral strike-slip faulting to the NNE-trending right-lateral strike-slip faulting and both having a thrust component. The current crustal movement direction revealed by the GPS velocity field is consistent with the direction of historical crust evolution of the lower crust revealed by the viscosity, implying that there is a good coupling relationship between the lower crust and upper crust. The results provide a basis for studying the development of fault systems with different strikes and properties, the formation of orogenic belts, the macroscopic geomorphological evolution characteristics, and the rheological and uplift dynamics of the lithosphere in the northeastern margin of the Tibetan plateau. In addition, our research differs from the previous studies in the spatial and temporal scale. Previous studies included either the entire Qinghai-Tibet Plateau or only the eastern margin of the Qinghai-Tibet Plateau. However, our analysis on the contours and topographical differences in the topography of the western Qinling-Songpan tectonic knot reveals that the study area is controlled by the lower crust flow. Our results are confirmed by various observations such as seismology, magnetotellurics and geophysical exploration. Moreover, the previous studies did not point out enough that the elevation contours are elliptical, and the elliptical geomorphology further illustrates that the formation and evolution of the Qinghai-Tibet Plateau has rheological characteristics and also conforms to the continuous deformation mode. Meanwhile, in terms of time scale, the evolution time of the study area is divided into three types of simulation time according to geochronology. And the GPS velocity field is introduced to observe the present-day crustal deformation.     

联合GPS和GRACE观测研究日本MW9.0地震震后变形机制

利用GPS和GRACE观测数据研究了日本M_W9.0地震的震后变形特征.GPS观测显示,区域震后位移呈现随指数函数变化特征,变化速率符合大森公式的衰减特性;近五年的震后水平位移累积已达到东向60~165 cm,南向20~65 cm的量值,距震中较远站点已超过同震变化量,且震后变形仍然持续.GRACE观测到显著的震后重力变化,地震破裂两侧的重力变化总体均呈上升趋势,但海洋侧的变化速率较快.联合震后余滑和黏弹性位错理论对震后变形进行了模拟,探索了GPS和GRACE观测的综合应用方法.研究发现,综合考虑震后余滑和黏滞性松弛效应可以对日本地震的震后变形做出较合理的解释,震后初期余滑起主要作用,1至2年以后逐渐减弱,黏滞性松弛作用逐渐增强.在震后变形模拟和区域黏滞性结构反演中形成GPS和GRACE观测结合应用的方法,先基于震后GPS形变估算区域黏滞性结构,而后利用GRACE观测修正深部的黏滞系数,并综合利用这两种观测微调浅层黏滞系数,最终确定区域黏滞性结构.基于该方法反演了日本震源区的地幔黏滞性结构,地震断层破裂两侧的流变参数存在差异,大陆侧的地幔顶层黏滞系数在1.0×10¹⁹ Pa·s量级,而海洋侧的则略小于大陆的,在6.0×10¹⁸ Pa·s量级.     

Mechanisms of Transient Postseismic Deformation Following the 2001 Mw 7.8 Kunlun (China) Earthquake

Using global positioning system (GPS) technology, significant postseismic surface displacements were observed within the first 4 months after the 2001 Mw 7.8 Kunlun earthquake which occurred in China. In this study, we investigated the mechanisms that may have possibly contributed to the postseismic deformations that have been observed. Based on the modeling results, we find that an afterslip model can interpret postseismic displacements in the near field even when the fault plane is extended to the bottom of the crust (~70 km). Models based on the viscoelastic relaxation theory showed a large discrepancy in the spatial pattern of the deformation compared with what has been observed. Thus, we infer that both mechanisms cannot interpret the observed postseismic deformation independently. A combination of afterslip and viscoelastic relaxation can further improve the data fit, especially at sites far from the fault. With maximum afterslip of ~0.4 m occurring at a depth of 10 km in the central section, the combined model shows that the estimated afterslip occurred mostly on and below the coseismic rupture plane, as well as on its eastern extension. The estimated moment released by the afterslip in the first 4 months is almost 40% of that released by the coseismic slip. The best-fitting viscoelastic relaxation model shows a “weak” upper mantle with a viscosity of ~1.0 × 10¹⁸ Pa s. The combined model also suggests the existence of a lower crust with viscosity larger than 1.0 × 10¹⁸ Pa s, although it cannot be constrained accurately.     

Numerical simulation on the influences of Wenchuan earthquake on the surrounding faults

On 12 May 2008, the devastating Wenchuan earthquake struck the Longmenshan fault zone, which comprised the eastern margin of the Tibetan Plateau, and this fault zone was predominantly a convergent boundary with a right-lateral strike-slip component. After such a large-magnitude earthquake, it was crucial to analyze the influences of the earthquake on the surrounding faults and the potential seismic activity. In this paper, a complex viscoelastic model of western Sichuan and eastern Tibet regions was constructed including the topography. Based on the findings of co-seismic static slip distribution, we calculated the stress change caused by the Wenchuan earthquake with the post-seismic relaxation into consideration. Our preliminary results indicated that: (1) The tectonic stressing rate was relatively high in Kunlun mountain pass-Jiangcuo, Ganzi-Yushu, Xianshuihe and Zemuhe faults; while in the east Kunlun and Longriba was medium; also the value was less in the Minjiang, Longmenshan, Anninghe and Huya faults. As to the Longmenshan fault, the value was 0.28×10^-3 MPa/a to 0.35×10^-3 MPa/a, which is coincident with the previous long recurrence interval of Wenchuan earthquake; (2) The Wenchuan earthquake not only caused the Coulomb stress decrease in the source region, but also the stress increase in the two terminals, especially the northeastern segment, which is comparatively consistent with the aftershock distribution. Meanwhile, the high concentration areas of the static slip distribution were corresponding to the Coulomb stress reductions; (3) The Coulomb stress change caused by Wenchuan earthquake showed significant increase on five major faults, which were northwestern segment of Xianshuihe fault, eastern Kunlun fault, Longriba fault, Minjiang fault and Huya fault respectively; also the Coulomb stress on the fault plane of the Yushu earthquake was faintly increased; (4) We defined the recurrence interval as the time needed to accumulate the magnitude of the stress drop, and the recurrence interval of Wenchuan earthquake was estimated about 1 714 a to 2 143 a correspondingly.     

Dynamic mechanisms of the post-seismic deformation following large events: Case study of the 1999 Chi-Chi earthquake in Taiwan of China

The mechanism of postseismic deformation related to strong earthquakes is important in geodynamics, and presumably afterslip or viscoelastic relaxation is responsible for the postsesimic deformation. The 1999 Chi-Chi, Taiwan of China, earthquake occurred in the region where GPS observation station is most densely deployed in the world. The unprecedented GPS data provides a unique opportunity to study the physical processes of postseismic deformation. Here we assume that the interactions of viscoelastic relaxation, afterslip, fault zone collapse, poroelastic rebound, flow of underground fluids, and all these combined contribute to the surface displacements following the main shock. In order to know the essence of the postseismic deformation after the strong event, fault zone collapse, poroelastic rebound, flow of underground fluids, and so on, are represented equivalently by the variations of the focal medium properties. Therefore, the viscoelastic relaxation, afterslip, and the variations of the equivalent focal medium properties are inverted by applying the GPS temporal series measurement data with viscoelastic finite element method. Both the afterslip rate distribution along the fault and the afterslip evolution with time are obtained by means of inversion. Also, the preliminary result suggests that viscosities of the lower crust and the upper mantle in Taiwan region is 2.7×10¹⁸ and 4.2×10²⁰ Pa·s, respectively. Moreover, the inversion results indicate that the afterslip contributing to postseismic deformation of 44.6% in 450 days after the Chi-Chi earthquake, with 34.7% caused by the viscous relaxation and 20.7% by other factors such as fault zone collapse, poroelastic rebound, and the flow of liquids.     

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

深地震反射剖面揭示的海原断裂带深部几何形态与地壳形变

由于活动的青藏高原不断的隆升和推挤作用,在西南向东北的推挤作用和周缘块体的阻挡以及东北缘内部块体挤压形变的作用下,形成了多个走向不同的青藏高原东北缘构造体系.新生代构造变形和地震活动强烈,区内分布多条大型深断裂带.海原断裂是青藏高原东北缘发育的弧形活动断裂带中规模最大、活动最为强烈的一条左旋走滑型断裂带,是重要的大地构造区边界,也是控制现今强震活动的活断层.本文利用2009年完成的高分辨率深地震反射剖面的北段资料,对其进行初步构造解释,揭示出海原断裂带的深部几何形态和其两侧地壳上地幔细结构.结果显示海原断裂并不是简单的陡立或者较缓,其几何形态随着深度变化.在海原断裂之下的Moho并未错断的反射特征显示海原断裂并不是直接错断莫霍面的超壳断裂.海原断裂带及两侧岩石圈结构和构造样式的研究为探讨青藏高原东北缘岩石圈变形机制提供地震学依据.     

Rupture behavior and deformation localization of the Kunlunshan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>7.8) and their tectonic implications

Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M _w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M _w =6.8, M _w =6.2 and M _w ⩽=7.8 events, respectively. The M _w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)     

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

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

A preliminary model for 3-D rheological structure of the lithosphere in North China

3-D rheological structure is mainly the spatial distribution of lithospheric strength or viscos-ity, its strength and viscosity are indispensable parameters in quantitative study of the lithosphere deformation. Plate tectonics theory initially divided the…     

青藏高原东缘中下地壳流与地壳变形

地壳缩短导致青藏隆升造山是普遍的认识.然而,在青藏东部,越来越多的观测数据和研究支持了中下地壳流与隆升造山的关系.目前,地壳缩短造山机制和中下地壳流造山机制仍然处于争论之中.本文建立了二维黏弹塑性有限元模型,模拟了龙门山断层带的多个地震循环的应变与变形,探讨了无与有中下地壳流情况下,地壳地表的位移、速度与变形的分布和演化;以及有中下地壳流情况下,不同流动范围、速度与黏度对模型结果的影响;并结合地形变观测数据的约束,推测了青藏东缘中下地壳流的流动状态.模拟结果显示,通过对比有和无中下地壳流的模拟结果,发现青藏东部震间的地表垂向速度在变形样式及数值上存在较大差异,即存在地壳流的地表垂向抬升速率显著大于无地壳流;震间在龙门山断层西侧附近产生的垂向凸状隆起随中下地壳流的速度、黏度及通道长度的变化而变化.此外,本文研究结果对青藏其他地区可能存在的地壳流的研究也具有一定的参考意义.     

鲜水河断裂带跨断层变形分析和数值模拟

鲜水河断裂带积累了20多年的跨断层变形资料, 如何解释这种变形是一个重要的问题。 该文采用了弹性上地壳覆盖在粘弹性Maxwell体之上的模型。 炉霍1973年M_S7.6地震按该粘弹性模型计算的震后形变曲线能够在很大程度上解释实测观测结果, 显示与一些断层(如圣安德烈斯断层某些段落)存在缓慢蠕滑而无大地震发生不同, 鲜水河断裂带观测到的跨断层变形, 有相当一部分可以用大地震后的粘弹性变形来解释。 当然, 计算得到的震后粘弹性形变幅度同跨断层短基线和短水准的实测形变幅度比较还存在一些差异, 不能完全解释, 表明了大陆断层活动的复杂性, 并值得进一步深入研究。     

Application of rock creep experiment in calculating the viscoelastic parameters of earth medium

The viscosity of earth medium can explain many geophysical phenomena. The uplifting velocity of Ti- betan Plateau may be caused by extrusion of the low viscosity channel in middle and lower crust[1]. The post-deformation after earthquake also relates to the rheology of earth medium[2]. The aftershock sequence is considered a process of strain release after stress dropped[3]. There are two kinds of models in the study of rheological character of crustal medium, elastoplas- only when the forced…