首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The implications of subdividing the lithosphere into upper elastic and lower viscoelastic layers are investigated by finite-element analysis. Application of uniform horizontal boundary stresses at the ends of a lithospheric plate leads to amplification of the stress in the elastic layer by a factor about equal to the ratio of lithospheric to elastic layer thicknesses, and the visco-elastic layer becomes nearly stress-free except near its ends. The time constant for approach to equilibrium is proportional to viscosity, being for our model 0.21 My for 1023 N s m−2, and there is some accompanying flexure of the lithosphere. Local variation in the thickness of the elastic layer causes inverse variation in the stress, in part explaining the stability of shield regions and the tectonic activity of hot plateau uplift regions. It is shown that stress amplification also occurs where the stresses arise from body forces such as differential loading and isostatic compensation across continental margins.  相似文献   

2.
The Aegean–Anatolian region is characterised by an inhomogeneous deformation pattern with high strain rates and a high seismicity both at the boundaries and in the plate interior. This pattern is controlled by the geometry and rheology of the structural units involved and their tectonic setting. A numerical analysis with a finite-element model of the region is used to quantify the influence of different rheological parameters. Viscoelastic material behaviour is implemented for the mantle lithosphere, whereas the crust is modelled with an elastic–plastic rheology. The variation of the inelastic material properties (viscosity and plastic strength) quantifies the influence of these material parameters on the deformation, stress, and strain patterns. Comparison of the modelled results with geodetic and geophysical observations reveals that the viscosity of the mantle lithosphere is the key to explaining the inhomogeneous deformation pattern. The best-fit model yields a viscosity of 1020 Pa s beneath Anatolia, whereas adjacent regions have viscosities between 1021 and 1023 Pa s. The model also explains the intra-plate seismicity and the stress field as well as its partitioning into regions with strike-slip and normal faulting. The final model is in good agreement with seismological, geodetic, and geological observations. Local deviations can be tracked down to small-scale structures, which are not included in the model.  相似文献   

3.
The relatively low elevation and thick crust in the Altiplano, in comparison to the higher elevation, but thinner crust in the Puna plateau, together with geophysical data, suggests that isostatic equillibrium is achieved by cooler and denser lithospheric mantle in the Altiplano. Excess density in the Altiplano mantle could create differential horizontal stress in the order of 25 MPa between both lithospheric columns. Numerical models accounting for pressure and temperature-dependent rheology show that such stress can induce horizontal ductile flow in the lower crust, from the Puna towards the Altiplano. With a minimum viscosity of 1019 Pa s, this flow reaches 1 cm/year, displacing more than 50 km of material within 5 Ma. If the lower crust viscosity is smaller, the amount of orogeny-parallel lower crustal flow can be even greater. Such a mechanism of channel flow may explain that different amounts of crustal material have been accommodated by shortening in the Altiplano and in the Puna. Because of the strength of the elastic-brittle upper crust, this channel flow does not necessitate large amounts of surface deformation (except vertical uplift), making it difficult to detect from the geology.  相似文献   

4.
The estimates of rheological thickness and total lithospheric strength for the Indian continental lithosphere have been obtained based on the representative rheological properties of upper crust, lower crust and upper mantle, and some of the available heat flow and heat generation data. The rheological thickness, computed at different locations in the Indian shield, shows lateral variation ranging from 79km in the southern part to 65 km in the northern part for a strain rate of 10-14 s-1. The total strength of the continental lithosphere is of the order of 1013 Nm-1 for the same value of strain rate and decreases northward. The computations carried out for a range of strain rates show an increase in the rheological thickness and strength of the lithosphere with increasing strain rate. These results would be important in understanding the flexural response of the Indian continental lithosphere to surface and subsurface loading, and response to tectonic forces acting on it.  相似文献   

5.
Conditions for indentation and channelised flow are investigated with two-dimensional thermomechanical models of Alpine-type continental collision. The models mimic the development of an orogen at an initial central portion of weakened lithosphere 150 km wide, coherent with several geological reconstructions. We study in particular the role of lower crustal strength in developing peculiar geometries after 20 Ma of shortening at 1 cm/year. Crustal layers produce geometries of imbricate layers, which result from two contrasted mechanisms of either channelised ductile lateral flow or horizontal rigid-like indentation:
– Channelised lateral flow develops when the lateral lower crust has a viscosity less than 1021 Pa s, exhibiting velocities opposite to the direction of convergence. This mechanism of deformation produces subhorizontal shear zones at the boundaries between the lower crust and the more competent upper crust and lithospheric mantle. It is also associated with a topographic plateau that equilibrates with a wide (about 200 km) but quasi-constant crustal root about 50 km deep.
– In contrast, indentation occurs with lateral lower crust layers that have a viscosity greater than about 1023 Pa s, producing significant shortening and thickening of the central crust. In this case topography develops steep and narrow (around 100 km wide), associated with a thickened crust exceeding 60 km depth. A crustal-scale pop-up forms bounded by subvertical shear zones that root into the mantle lithosphere.
Keywords: Continental collision indentation; Channel flow; Lower crust; Elastic–viscous–plastic rheology  相似文献   

6.
Deformation models used to explain the triggering mechanism often assume pure elastic behaviour for the crust and upper mantle. In reality however, the mantle and possibly the lower crust behave viscoelastically, particularly over longer time scales. Consequently, the stress field of an earthquake is in general time-dependent. In addition, if the elastic stress increase were enough to trigger a later earthquake, this triggered event should occur instantaneously and not many years after the triggering event. Hence, it is appropriate to include inelastic behaviour when analysing stress transfer and earthquake interaction.In this work, we analyse a sequence of 10 magnitude Ms > 6.5 events along the North Anatolian Fault between 1939 and 1999 to study the evolution of the regional Coulomb stress field. We investigate the triggering of these events by stress transfer, taking viscoelastic relaxation into account. We evaluate the contribution of elastic stress changes, of post-seismic viscoelastic relaxation in the lower crust and mantle, and of steady tectonic loading to the total Coulomb stress field. We analyse the evolution of stress in the region under study, as well as on the rupture surfaces of the considered events and their epicentres. We study the state of the Coulomb stress field before the 1999 İzmit and Düzce earthquakes, as well as in the Marmara Sea region.In general, the Coulomb stress failure criterion offers a plausible explanation for the location of these events. However, we show that using a purely elastic model disregards an important part of the actual stress increase/decrease. In several cases, post-seismic relaxation effects are important and greater in magnitude than the stress changes due to steady tectonic loading. Consequently, viscoelastic relaxation should be considered in any study dealing with Coulomb stress changes.According to our study, and assuming that an important part of the rupture surface must be stressed for an earthquake to occur, the most likely value for the viscosity of the lower crust or mantle in this region is 5 · 1017–1018 Pa · s. Our results cannot rule out the possibility of other time-dependent processes involved in the triggering of the 1999 Düzce event. However, the stress increase due to viscoelastic relaxation brought 22% of the 1999 Düzce rupture area over the threshold value of Δσc ≥ 0.01 MPa (0.1 bar), and took the whole surface closer to failure by an average of 0.2 MPa. Finally, we argue that the Marmara Sea region is currently being loaded with positive Coulomb stresses at a much faster rate than would arise exclusively from steady tectonic loading on the North Anatolian Fault.  相似文献   

7.
粘弹性数值模拟龙门山断裂带应力积累及大震复发周期   总被引:4,自引:0,他引:4  
柳畅  朱伯靖  石耀霖 《地质学报》2012,86(1):157-169
2008年5月12日在低地形变速率的龙门山断裂带上突发汶川强震,引发人们对该地震孕震机制的思考。本文根据GPS观测资料确定边界条件,通过三维粘弹性数值模拟探讨了汶川地震的孕震机理,计算了该区域岩石圈的应力增加速率和积累过程,以及汶川地震同震应力变化与震后应力松弛,在此基础上估算了汶川8.0级大地震的复发周期。数值模拟结果表明:印度板块对欧亚板块的推挤造成青藏高原的物质东流,高原中、下地壳物质在龙门山断裂带处遭到相对坚硬的四川盆地的阻挡之后,部分中、下地壳物质在龙门山断裂带下堆积产生应力集中。两个重要因素为应力集中提供了重要控制作用:其一是青藏高原中、下地壳较低的粘滞系数与四川盆地中、下地壳较高的粘滞系数的差异,其二是从青藏高原到四川盆地的Moho面深度在龙门山断裂带的突变。低应变速率的龙门断裂带岩石圈在数千年时间尺度的应力积累过程中,脆性上地壳的应力随时间近乎线性增长,并且上地壳深部的应力增长率超过浅部,6000年内应力积累最大量达到-21.6MPa,应力增长速率为-0.0036MPa/a;而柔性的中、下地壳以及岩石圈上地幔的应力在增长一段时间之后趋于稳定。在空间上,龙门山断裂带受到的压应力从断层西南向北东方向逐渐减小,而剪应力从西南到北东方向逐渐增大,应力状态有利于地震发生时断层的破裂方式从西南的逆冲运动向北东的逆冲兼走滑运动的方式发展。通过应力积累与地震应力降的计算得到汶川8.0级大地震的复发周期约为5400年。  相似文献   

8.
Pinch and swell structures occur where a more competent layer in a weaker matrix is subjected to layer-parallel extension. In this contribution, we use numerical models to explore the use of pinch and swell structure shape symmetry and asymmetry as a determinant of relative viscosity between layers. Maximum asymmetry is attained when the matrix viscosity on one side is subtly weaker than the competent layer, while the other side is significantly weaker.Our numerical results are directly applied to asymmetrically developed pinch and swell structures in exposed lower continental crust. Here, shape geometries observed in a shear zone comprised of plagioclase-dominated, garnet-dominated and mixed amphibole-plagioclase-dominated bands, reveals that the plagioclase-dominated band is the most competent band and is marginally stronger (2×) and significantly stronger (10–40×) than the fine grained garnet-dominated and mixed amphibole-plagioclase-dominated band, respectively. Based on the experimentally determined viscosity of a plagioclase-dominated material and quantitative microstructural analysis, the viscosity range of the natural rock bands is 2.8 × 1015 to 1.1 × 1017 Pa s. Consequently, the assumption that the experimentally-derived plagioclase flow law is an appropriate proxy for the middle to lower continental crust may lead to a viscosity over-estimation by up to forty times.  相似文献   

9.
An empirical formula for the fracture strength of the principal rock type in the lithosphere is obtained based on the experimental data from previous studies, in which the effects of the confining pressure, size of the rock sample, temperature, strain rate and the pore pressure are taken into account, the empirical formulae for the effects of them are also presented. By comparing the frictional strength to the fracture strength calculated using the new empirical formula, it is shown that frictional sliding is dominant in the upper crust but brittle fracture is dominant in the lower part of the crust and the lithosphere beneath the crust. Therefore the fracture mechanism must be taken into account in the study of the rheological structure of the lithosphere. The empirical formula for the fracture strength is applied to study the rheological structure of the lithosphere in the Ordos block. Brittle regime in the rheological structure can be divided into two sub-regions, in one of which brittle fracture and in the other frictional sliding are dominant, respectively, unlike previous conventional studies in which frictional sliding is assumed to be the only factor; the magnitude of the rheological strength of the lithosphere calculated by the empirical formula is also lower than that obtained in previous conventional studies.  相似文献   

10.
略论引潮力的交变应力效应及其地球动力学意义   总被引:2,自引:1,他引:1  
地球在日、月引潮力作用下,产生周期性引潮力与西向引潮力。周期性引潮力所致构造应力达到104Pa左右的量级,属对称循环交变应力,对应于105Pa量级的静态构造应力,对地震等构造过程能起显著的触发作用。西向引潮力所致西向剪应力达到105Pa左右的量级,属脉动循环交变应力,对应于106Pa量级的静态构造应力,对岩石圈、软流圈水平运动有显著的驱动意义,是构造运动基本驱动力之一。   相似文献   

11.
Deformation within the Earth's lithosphere is largely controlled by the rheology of the rock. Fracture and faulting are characterized by elastic rheologies with brittle mechanisms, while folding and flow are characterized by plastic and/or viscous rheologies due to ductile mechanisms. However, it has been recognized that deformation that resembles ductile behavior can be produced within the confines of the brittle lithosphere. Specific examples are folds that form in the shallow crust, steep hinges at subduction zones that are accompanied by seismicity, and large-scale deformation at plate boundaries. In these cases, the brittle lithosphere behaves elastically with fracture and faulting yet produces ductile behavior. In this paper, we attempt to simulate such ductile behavior in elastic materials using continuum damage mechanics. Engineers utilize damage mechanics to model the continuum deformation of brittle materials. We utilize a modified form of damage mechanics that represents a reduction in frictional strength of preexisting fractures and faults. We use this empirical approach to simulate the bending of the lithosphere under the application of a constant moment.We use numerical simulations to obtain elastostatic solutions for plate bending and where the longitudinal stress at a particular node exceeds a yield stress, we apply damage to reduce Young's modulus at the node. Damage is calculated at each time step by a power-law relationship of the ratio of the yield stress to the longitudinal stress and the yield strain to the longitudinal strain. This results in the relaxation of the material due to increasing damage. To test our method, we apply our damage rheology to an infinite plate deforming under a constant bending moment. We simulate a wide range of behaviors from slow relaxation to instantaneous failure, over timescales that span six orders of magnitude. Using this method, stress relaxation produces elastic-perfectly plastic behavior in cases where failure does not occur. For cases of failure, we observe a rapid increase in damage leading to failure, analogous to the acceleration of microcrack formation and acoustic emissions prior to failure. The changes in the rate of damage accumulation in failure cases are similar to the changes in b-values of acoustic emissions observed in triaxial compression tests of fractured rock and b-value changes prior to some large earthquakes. Thus continuum damage mechanics can simulate the phenomenon of ductile behavior due to brittle mechanisms as well as observations of laboratory experiments and seismicity.  相似文献   

12.
亚洲中东部岩石圈塑性流动网络与构造应力场   总被引:2,自引:0,他引:2  
:板块边缘驱动作用主要通过岩石圈下层的网络状塑性流动实现其向板内的远程传递 ,并控制上层的构造变形。相应地造成构造应力场的的多层分布格局 ,主要包括下层 (含下地壳和岩石圈地幔 )、多震层 (属于上部地壳 )和浅层 (浅表地壳 )应力场。这些层中最大压应力方向可分别采用塑性流动网络共轭角平分线、震源机制解和井孔应力测量等方法加以确定。根据对亚洲中—东部下层与多震层应力场的对比表明 ,大部分地区二者的应力方向基本一致 ,表明了下层对上层的控制 ,这对了解板内应力场的一般趋势与格架提供了基础。靠近驱动边界的局部地段 ,由于上部地壳直接的传递作用 ,往往导致上、下层应力方向的显著偏差 ,有可能成为进一步了解深浅层联系和差异的一种方法  相似文献   

13.
J. Ebbing   《Tectonophysics》2004,380(1-2):89-104
The transition between European and Adriatic crust is an important feature related to the plate collision that formed the European Alps. The diversity of seismic and geological information allows the construction of two alternative 3D density models, which both match the observed gravity field. Different seismic experiments suggest a thickness for the Adriatic crust between 30 and 40 km. The thick crust model requires an unusually dense lower crust (>3050 kg/m3) to reproduce the observed Bouguer anomaly. To evaluate the two alternative models, the isostatic implications of the geometry and density distribution within both 3D models are investigated, using local (Airy) and regional (Vening Meinesz) isostasy.Airy isostatic investigations show that the Eastern Alps are not isostatically compensated and the residuals correlate strongly with exposed geological formations. Subsequently, subsurface loading is an important factor controlling isostatic processes. The different geometry and densities in the two 3D models imply different loading at the crust–mantle boundary. The subsurface loads calculated from the 3D density models were used to estimate regional isostasy by a convolution method. In general, small rigidity values (D<10×1021 Nm) are determined for the Eastern Alpine lithosphere. In the model with a 40-km-thick Adriatic crust, high flexural rigidities are inferred for the Adriatic plate (>100×1021 Nm), but these values are unusual for an active orogenic region. The results point to the interfingering of European and Adriatic crust that results in the squeezing of European crust between Adriatic crust and mantle with additional contamination by mantle material.  相似文献   

14.
大别造山带构造超压形成的碰撞力学机理   总被引:3,自引:0,他引:3  
提出了大别山构造超压形成的点碰撞模型,简要分析了大陆碰撞带构造运动引起的粘性介质中粘性应力和平均应力随岩石物性的变化规律。探讨了构造压力对超高压的贡献及对成岩深度的重要意义。研究表明:构造运动引起的岩石圈中的附加压力可能与静岩压力有相同的数量级,大陆造山带两陆块不规则边界的碰撞会引起局部应力集中,产生较大的构造压力,岩石介质的流变学分析表明,在相同外力作用下,岩石圈上部的高粘度性质决定了其在构造活动期间增温效果显著,但增压效果有限;而粘性较低的岩石圈下部则增压效果明显,为此,在下地壳与上地幔之间的低粘度带内有可能发生超高压变质作用。  相似文献   

15.
During continental collision in the middle Silurian, the thickness of the lithosphere under the Caledonides of S. Norway was doubled by subduction of the western margin of Baltica, including the Western Gneiss Region, under Laurentia. Crustal rocks of the Baltic plate reached sub-Moho depths of near 100 km or more as inferred from the presence of coesite in eclogites. Isostatic calculations indicate an average elevation of the mountain chain of about 3 km at this stage. The subducted lithosphere experienced vertical constrictional strains as a result of slab-pull by its heavy and cold root. Eduction of the deeply buried crustal material was initiated by decoupling of the Thermal Boundary Layer in the subducted lithosphere. Isostatic rebound resulted in very rapid uplift (1–2 mm yr-1), and the deep crust was exhumed, mainly by tectonic extensional stripping over a period of 30–40 Myr. The eduction was probably related to a rolling hinge, footwall uplift mechanism, and the early high-pressure coaxial fabrics were overprinted by extensional simple shear as the deep crust reached middle and upper crustal levels. The model explains the present-day normal crustal thickness under the exhumed deep rocks without necessarily invoking large-scale lateral flow of material in the lower crust or igneous underplating.  相似文献   

16.
拆沉作用(delamination)是地球科学中一个重要的科学问题。本文认为,大洋岩石圈拆沉和大陆下地壳拆沉是不一样的:(1)拆沉的物质不同。大洋岩石圈拆沉的物质包括大洋地壳、岩石圈地幔甚至一部分软流圈地幔,它们共同进入地幔深部;而大陆下地壳拆沉仅仅限制在下地壳,不包括岩石圈地幔。(2)拆沉的动力不同。大洋岩石圈拆沉是由板块俯冲引起的,是地幔对流的产物,因此是一种快速的主动的拆沉;而下地壳拆沉是由于下地壳加厚使下地壳密度增加引起的,还要求其下刚性的岩石圈地幔转变成塑性的软流圈地幔才有可能发生。因此下地壳拆沉要克服许多阻力才能实现,使拆沉成为一个漫长的过程,是慢速的和被动的拆沉。(3)拆沉的过程不同。大洋岩石圈拆沉是由板块俯冲触发的,俯冲导致碰撞,大洋岩石圈从根部断裂,拆沉进入地幔。大陆下地壳拆沉由地壳加厚开始,使下地壳转变为榴辉岩相;随后,岩石圈地幔减薄,直至全部转化为软流圈地幔;下地壳发生部分熔融,形成大规模的(埃达克质)岩浆,使下地壳榴辉岩的密度大于下伏的地幔,从而引发拆沉。大陆下地壳拆沉不大可能是整体进行的,可能是一块一块地被蚕食、被拆沉的。(4)拆沉后的效应不同。大洋岩石圈地幔拆沉,使热的软流圈地幔上涌,从而引发了一系列地质效应:如岩浆活动、地壳抬升、构造松弛以及随后的造山带垮塌等。而下地壳拆沉只引起地壳减薄,高原和山脉垮塌,并不伴有大规模的岩浆活动和地壳抬升等过程。(5)拆沉与岩浆活动的关系不同。主动拆沉导致大规模岩浆活动,而被动拆沉是在大规模岩浆活动的基础上开始的。此外,文中还对"下地壳 岩石圈地幔拆沉"模式提出了质疑,认为该模式有许多难以理解的问题和太多推测的成分,而且与现在保存的地质事实不符。  相似文献   

17.
华北克拉通晚中生代壳-幔拆离作用: 岩石流变学约束   总被引:6,自引:5,他引:1  
大陆岩石圈的流变学结构对于岩石圈深部过程(壳/幔过程)有着深刻的影响,直接表现在岩石圈壳-幔结构与浅部构造上.本文注意到华北克拉通晚中生代岩石圈减薄期间地壳的伸展、拆离与减薄在不同地区的宏观、微观构造及地壳岩石流变学等方面的差异表现与区域变化,以及现今和晚中生代时期岩石圈厚度的不均匀性.讨论了以水为主体的地质流体的存在对于岩石圈流变性的影响.综合克拉通东部与西部地壳/地幔厚度变化特点以及下地壳和上地幔含水性特点,阐述了晚中生代时期华北克拉通岩石圈内部壳幔耦合与解耦的规律,提出了华北岩石圈壳-幔拆离作用模型以解释华北克拉通晚中生代岩石圈减薄的基本现象与深部过程.提出区域性伸展作用是岩石圈减薄的主要动力学因素,东部地区在晚中生代伸展作用过程中壳-幔具有典型的解耦性,上部地壳、下部地壳和岩石圈地幔的变形具有显著差异性.而西部区壳幔总体具有耦合性,下地壳与岩石圈地幔共同构成流变学强度很高且难以变形的岩石圈根.  相似文献   

18.
用不可压缩固体模型对岩石圈内横向密度差异产生的构造应力和缓慢变形进行了三维有限元模拟。结果表明,对于达到均衡补偿状态的模型,5公里地形高差的重力势在地壳内产生数百巴的构造应力,高原下出现横向拉张应力,平原下出现挤压应力;地壳厚度变化大的过渡带,应力的增加更显著。在地壳局部增厚而地形高度未达到均衡补偿的情况下,深部侧向密度差异产生的重力调整作用使物质向地壳较厚的地方迁移,导致那里的地壳产生地表隆起。  相似文献   

19.
中国大陆岩石圈等效粘滞系数的计算和讨论   总被引:27,自引:0,他引:27  
大陆岩石圈的流变结构对岩石圈动力学过程有很大的影响,因此对岩石圈等效粘度的估计是大陆动力学研究中基础和重要的问题。文中对利用实验室流变实验结果估算岩石圈流变结构的计算方法中包含的多种不确定性进行了讨论,包括岩性、温度、应变速率、实验室速率数据外推到地质构造运动速率等因素对等效粘滞系数估算的影响,并以温度和应变速率的新研究成果为基础,对中国大陆地壳和上地幔等效粘滞系数做出了估计。中国中地壳等效粘滞系数一般在1021~1024Pa.s,下地壳等效粘滞系数在1021~1022Pa.s,其中青藏高原下地壳等效粘滞系数较低,约为1019~1020Pa.s;与前人研究认为青藏高原存在柔性下地壳流动的结论吻合。  相似文献   

20.
A two‐dimensional thermorheological model of the Central Alps along a north–south transect is presented. Thermophysical and rheological parameters of the various lithological units are chosen from seismic and gravity information. The inferred temperature distribution matches surface heat flow and results in Moho temperatures between 500 and 800 °C. Both European and Adriatic lithospheres have a ‘jelly‐sandwich’ structure, with a 15–20 km thick brittle upper crust overlying a ductile lower crust and a mantle lid whose uppermost part is brittle. The total strength of the lithosphere is of the order of 0.5–1.0 × 1013 N m−1 if the upper mantle is dry, or slightly less if the upper mantle is wet. In both cases, the higher values correspond to the Adriatic indenter.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号