西北太平洋岩石圈有效弹性厚度及其构造意义

本文引入滑动窗口导纳技术(MWAT),计算西北太平洋岩石圈有效弹性厚度(Te).首先,基于SIO V15.1海底地形模型,模拟研究了MWAT法计算Te的精度,表明当TeTe≥5 km时,相对误差在10%以内.分别采用GEBCO、SIO V15.1和BAT_VGG海底地形模型,构建了西北太平洋Te,通过对获得的洋壳密度参数和实测导纳与模型导纳之差的均方根进行分析,结果表明,BAT_VGG模型更适用于Te计算.西北太平洋Te均值为13.2 km,标准差为6.9 km,以板块冷却模型为参考,主要分布在150℃~450℃等温线深度范围内.白垩纪和侏罗纪时期岩石圈Te分布在150℃~300℃等温线深度范围内,且未随海山加载时岩石圈年龄增大而增大,说明海山加载时岩石圈年龄不是影响其强度的唯一因素.南太平洋超级海隆活动,以及研究区域广泛存在的断裂带构造,都曾对本区域岩石圈演化产生过重要影响,可能是本地区岩石圈Te较小的构造原因.     

海洋岩石圈板块有效弹性厚度研究

本文在前人研究大陆岩石圈板块有效弹性厚度的基础上,建立研究海洋岩石圈板块有效弹性厚度的理论模型,推导出与大陆岩石圈不同的海洋岩石圈板块响应函数 Z(k,Te) 理论计算公式.并分析海洋岩石圈板块响应函数 Z(k,Te) 的特点.文中对实际的海洋测量数据的响应函数 Z(k,Te) 进行计算和分析,估算我国南海南沙海域和南海中央海盆岩石圈板块有效弹性厚度分别约为10 km和6~7 km.     

Effective elastic thickness of island arc lithosphere under Japan

Abstract Using topography and observed gravity anomalies, we have estimated the effective elastic thickness as a measure of strength of Japanese island arc lithosphere. The thickness is found to range from about 3 km to >20 km. The thickness seems to be controlled primarily by the thermal state of the lithosphere. The higher the heat flow, the thinner is the elastic plate. However, several areas show significant deviations. The smaller effective elastic thickness in the northern Ryukyu arc than that inferred from heat flow may be attributed to the stress regime. In Japan, extensional tectonics are going on only in the Ryukyu arc region. Shallow subducting slab under the south-western Japan frontal arc probably increases the effective thickness by several kilometers. The determined effective elastic thickness suggests that when we consider vertical movements in the volcanic arc, we should take account of topographic and subsurface loading over a few hundred kilometers. However, if the dip of the slab is shallow, the flexural responses of the underlying slab, not only that of the island arc lithosphere, should be taken into account for the compensation, as is the case of the south-western Japan frontal arc.     

The effective elastic thickness of the lithosphere in the collision zone between Arabia and Eurasia in Iran

The effective elastic thickness, T_e, has been calculated in the collision zone between Arabia and Eurasia in Iran from the wavelet coherence. The wavelet coherence is calculated from Bouguer anomalies and topography data using the isotropic fan wavelet method, and gives T_e values between 14.2 and 62.2 km. The lower value is found in the Central Iranian Blocks and the East Iranian Belt which are bounded by several large strike-slip faults with lithospheric origin. The higher value occurs in the east of the South Caspian Sea Basin. The resulting T_e map shows positive and negative correlation with shear wave velocity and surface heat flow, respectively. A comparison between the seismogenic thickness (T_s) and T_e in Iran suggests that T_e > T_s. Results of the load ratio in Iran indicate that in most of the study area surface loads are much more prevalent than subsurface loads, except in the Central Iranian Blocks and NW of Iran. Intermediate to low T_e values in Iran were inherited from multiple rifting and orogenic activities from Late Precambrian (∼650 Ma) to present day which are not only reflected in thin and warm lithosphere but also an increasing seismicity rate.     

Temperature control of continental lithosphere elastic thickness,Te vs Vs

The elastic thickness of the continental lithosphere is closely related to its total strength and therefore to its susceptibility to tectonic deformation and earthquakes. Recently it has been questioned whether the lithosphere thickness and strength are dependent on crustal and upper mantle temperatures and compositions as predicted by laboratory data. We test this dependence regionally by comparison in northwestern North America of the effective elastic thickness Te, from topography–gravity coherence, with upper mantle temperatures mapped by shear wave tomography velocities Vs and other temperature indicators. The Te values are strongly bimodal as found globally, less than 20 km for the hot Cordillera backarc and over 60 km for the cold stable Craton. These Te correspond to low Vs beneath the Cordillera and high Vs beneath the Craton. Model temperature-depth profiles are used to estimate model Te for comparison with those observed. Only limited areas of intermediate thermal regimes, i.e., thermotectonic ages of ~ 300 Ma, have intermediate Te that suggest a weak lower crust over a stronger upper mantle. There are large uncertainties in model Te associated with composition, water content, strain rate, and decoupling stress threshold. However, with reasonable parameters, model yield stress envelopes correspond to observed Te for thermal regimes with 800–900 °C at the Cordillera Moho and 400–500 °C for the Shield, in agreement with temperatures from Vs and other estimators. Our study supports the conclusion that lithosphere elastic thickness and strength are controlled primarily by temperature, and that laboratory-based rheology generally provides a good estimate of the deformation behaviour of the crust and upper mantle on geological time scales.     

利用Moho地形导纳法(MDDF)反演中国大陆岩石圈有效弹性厚度

本文用数值实验方法验证了,在Moho地形导纳法(MDDF)中使用先由重力数据反演的Moho面相对起伏数据,能较传统的重力地形导纳法获得更高的反演精度.之后作者应用Moho地形导纳法(MDDF)反演获得了精度较高的中国大陆区域岩石圈有效弹性厚度.结果显示:(1)中国地区岩石圈有效弹性厚度Te从东向西大体上呈阶梯状上升;(2)Te与岩石圈地震-热厚度、地表热流、上地幔顶部的地震波速度等数据密切相关;(3)在中国中东部, Te较低的区域以及Te的高低值转换带对应较强的地震活动性.但在西部,Te与地震活动性的相关性并不明显.     

Effective elastic thickness of the lithosphere from joint inversion in western China and its implications*

The western China lies in the convergence zone between Eurasian and Indian plates. It is an ideal place to study the lithosphere dynamics and tectonic evolutions on the continental Earth. The lithospheric strength is a key factor in controlling the lithosphere dynamics and deformations. The effective elastic thickness (T_e) of the lithosphere can be used to address the lithospheric strength. Previous researchers only used one of the admittance or coherence methods to investigate the T_e in the western China. Moreover, most of them ignored the internal loads of the lithosphere during the T_e calculation, which can produce large biases in the T_e estimations. To provide more reliable T_e estimations, we used a new joint inversion method that integrated both admittance and coherence techniques to compute the T_e in this study, with the WGM2012 gravity data, the ETOPO1 topographic data, and the Moho depths from the CRUST1.0 model. The internal loads are considered and investigated using the load ratio (F). Our results show that the joint inversion method can yield reliable T_e and F values. Based on the analysis of T_e and F distributions, we suggest (1) the northern Tibetan Plateau could be the front edge of the plate collision of Eurasian and Indian plates; (2) the southern and part of central Tibetan Plateau have a strong lithospheric mantle related to the rigid underthrusting Indian plate; (3) the southeastern Tibetan Plateau may be experiencing the delamination of lithosphere and upwelling of asthenosphere.     

利用Moho面起伏及地表地形数据反演岩石圈有效弹性厚度的莫霍地形导纳法(MDDF)

区域岩石圈有效弹性厚度Te的确定对认识其力学性质及演化等动力学问题具有重要意义.长期以来,大陆地区有效弹性厚度的确定一直存在着争议.尽管Pérez-Gussinyé等2004年的研究消除了不同方法(自由空气导纳法和布格相关法)之间的分歧,但其确定有效弹性厚度的误差仍然非常大.本文提出了利用Moho面起伏及地表地形反演岩石圈有效弹性厚度的Moho地形导纳法(MDDF),并利用合成的地表地形和Moho面起伏模型验证了该方法的可行性.结果表明,与传统的重力地形导纳法相比,使用Moho地形导纳法(MDDF)反演岩石圈有效弹性厚度,能较好地提高反演精度.     

青藏高原东北缘岩石圈有效弹性厚度与地震分布关系

岩石圈有效弹性厚度（T_e）表征岩石圈综合力学强度,对理解区域深部构造孕震环境具有重要意义.青藏高原东北缘地质构造复杂,强烈地震活动频发.为进一步了解该区域深部岩石圈力学性质特征及其与地震活动的关系,本文利用基于Fan小波的谱方法,使用WGM2012重力异常数据、ETOPO1地形数据和CRUST1.0模型,通过导纳和相关联合反演计算了青藏高原东北缘的岩石圈有效弹性厚度（T_e）与荷载比（F）.结果显示研究区T_e整体呈明显的东高西低分布,青藏块体T_e变化剧烈,西部高值（>40 km）和东部低值（< 20 km）共存;鄂尔多斯地块T_e较高（>30 km）,变化相对平缓;荷载比F存在局部西南高,巴颜喀拉和羌塘地块荷载比F较高（>0.5）,说明以地下荷载为主,其他地块荷载比F较低（< 0.2）,以地表荷载为主.鄂尔多斯地块结构稳定,岩石圈强度较大,T_e值较高;内部构造活动性微弱,深部物质密度横向变化较小,岩石圈所受荷载以地表为主,荷载比F较低.柴达木地块东部与巴颜喀拉地块东部可能发生有下地壳流,岩石圈较为软弱,T_e值较低;巴颜喀拉地块内部与鲜水河断裂的荷载比F较高,岩石圈所受荷载以地下为主,可能是地壳物质向东流动导致岩石圈深部物质密度横向变化较大引起的.通过对比分析研究区T_e和T_e梯度、F和F梯度与地震之间的对应关系,发现地震一般发生在区域相对低T_e值、T_e梯度值和荷载比F梯度值较大的区域,震级较小的地震多发生在荷载比F较小的地区,6级以上强震则较多发生在荷载比F约为0.1和0.8的区域.因此,岩石圈有效弹性厚度T_e值较低、T_e梯度和荷载比F梯度较大的地区,更易出现地震.

     

A comparison of the estimated effective elastic thickness of the lithosphere using terrestrial and satellite-derived data in Iran

The effective elastic thickness of the lithosphere has an important role in constraining compositional structure, geothermal gradient and tectonic forces within the lithosphere and the thickness of this layer can be used to evaluate the earthquakes’ focal depth. Hence, assessment of the elastic thickness of the lithosphere by gravitational admittance method in Iran is the main objective of this paper. Although the global geopotential models estimated from the satellite missions and surface data can portray the Earth’s gravity field in high precision and resolution, there are some debates about using them for lithosphere investigations. We used both the terrestrial data which have been provided by NCC (National Cartographic Center of Iran) and BGI (Bureau Gravimetrique International), and the satellite-derived gravity and topography which are generated by EIGEN-GL04C and ETOPO5, respectively. Finally, it is concluded that signal content of the satellite-derived data is as rich as the terrestrial one and it can be used for the determination of the lithosphere bending.     

The effective elastic lithosphere under the Cook-Austral and Society islands

We have determined the elastic thicknessT_e of the oceanic lithosphere along two volcanic chains of the South Central Pacific: Cook-Austral and Society islands. We used a three-dimensional spatial method to model the lithospheric flexure assuming a continuous elastic plate. The model was constrained by geoid height data from the SEASAT satellite.Along the Cook-Austral chain the elastic thickness increases westward, from 2–4 km at McDonald hot spot to 14 km at Rarotonga. At McDonald seamount, however, the data are better explained by a local compensation model. The observed trend shows an increase ofT_e with age of plate at loading time. However, the elastic layer under the Cook-Austral appears systematically thinner by several kilometers than expected for “normal” seafloor, suggesting that substantial thermal thinning has taken place in this region. Considering the apparent thermal age of the plate instead of crustal age improves noticeably the results. Along the Society chainT_e varies from 20 km under Tahiti to 13 km under Maupiti which is located 500 km westward. When plotting together the Society and Cook-AustralT_e results versus age of load, we notice that within the first five million years after loading,T_e decreases significantly while tending rapidly to an equilibrium value. This may be interpreted as the effect of initial stress relaxation which occurs just after loading inside the lower lithosphere and suggests that the presently measured elastic thickness under the very young Tahiti load ( 0.8 Ma) is not yet the equilibrium thickness.     

西太平洋岩石圈有效弹性厚度的空间分布与控制因素

西太平洋地区板块间相互作用强烈,热演化和构造演化过程复杂.为了揭示构造相互作用对岩石圈强度的影响,本文使用自由空气重力异常模型WGM2012和地形模型ETOPO1,基于小波变换的导纳法计算得到了该地区的岩石圈有效弹性厚度(Te).西太平洋区域的Te主要分布在5～85 km之间,南海等张裂环境地区Te普遍小于20 km,俯冲带附近Te一般大于80 km,与俯冲板片年龄呈正相关.参照平板冷却模型,弹性岩石圈底界面主要分布在200～500℃等温面之间,随洋壳年龄增大逐渐趋于平稳,热点及年轻洋壳部分地区弹性岩石圈底界面处于200℃等温面之上.西太平洋海山与年轻海盆等区域Te与居里点深度一般呈正相关,与地表热流一般呈负相关,但由于强烈的构造运动、热液循环、岩浆活动、地幔流变性等因素的影响,整体Te与居里点深度和地表热流所反映的岩石圈热结构相关性不高.     

Heat flow and thickness of the oceanic lithosphere

Existing thermal models of the oceanic lithosphere predict too sharp an increase of heat flow towards the ridge axis. A new mathematical model of a thickening lithosphere is presented. The temperature distribution is computed by the use of observed surface heat flow as a boundary condition. If observed heat flow values represent flow of heat from the mantle, the model predicts a rather rapid growth of the lithosphere within the first 30 m.y. and a nearly steady state after 100 m.y. Heat flow from the asthenosphere to the lithosphere shows a minimum near the ridge axis, suggesting a down-going convective flow in the asthenosphere at both sides of a spreading center.     

中国陆地热岩石圈厚度及其地球动力学意义

依靠最新的中国大地热流数据、精细的地壳分层结构, 通过求解一维稳态热传导方程获得各个热流测量点对应的热岩石圈厚度, 通过克里金插值法绘制中国陆地热岩石圈厚度分布等值线图.计算结果表明, 中国陆地各构造区的热岩石圈厚度差异较大, 稳定的克拉通地区最厚, 可达200 km以上, 造山系次之, 多在100~200 km之间, 破坏的克拉通地区岩石圈最薄, 可以低于100 km. 通过对比三大克拉通地区的热岩石圈厚度和地震岩石圈厚度, 得出了四点认识: (1) 塔里木克拉通西部、中上扬子克拉通、华北克拉通西部以及南华北基本保留了稳定的克拉通巨厚岩石圈特征, 而华北克拉通东部的渤海湾盆地、下扬子克拉通以及塔里木克拉通东南部则发生了大规模的减薄; (2) 华北克拉通西部从鄂尔多斯向东北的银川—河套凹陷及向东南的汾渭凹陷的岩石圈厚度和流变边界层厚度逐渐变薄, 主要受控于地幔对流强度的增强; (3) 华北克拉通东部的南华北依然保持稳定, 而渤海湾的岩石圈厚度减薄显著, 体现了华北克拉通破坏在空间上的不均匀性; (4) 扬子克拉通自西向东岩石圈厚度和流变边界层厚度逐渐变薄, 可能受控于太平洋板块的俯冲, 和华北克拉通东部经历了相似的地球动力学过程.

     

Lithospheric thickness in the western Pacific

Rayleigh wave phase velocities were determined by the single-station method for ten paths in the western Pacific. The data show that even 100 m.y. after formation, the phase velocity and upper-mantle structure are still dependent upon age. Inversion of the data gave a model with a lithospheric thickness of 76 km at 100 m.y., increasing to 104 km at 150 m.y., measured from the base of the crust.     

A model of the lithosphere thickness in the region of the Carpathians

Summary Directionally independent average P residuals computed for waves of teleseismic events arriving under various azimuths and incidence angles provided the basis for estimating the lithosphere thickness beneath the Carpathians and their surroundings. A thin lithosphere (60–80 km) was determined for the Pannonian Basin and the Transylvanian Basin, the thickest lithosphere(about 180 km) beneath the South Carpathians at the contact with the Moesian Platform. In other parts of the Carpathian belt the lithosphere thickens beneath the outer parts towards the SW margin of the East-European and the Moldavian Platforms. The lithosphere thicknesses derived from P residuals correlate well with the magnetotelluric determinations of a layer of increased electrical conductivity in the upper mantle.

---

Резюме Осре?rt;ненные временные невязкu Р волн, незaвuсuмые оm нanрaвленuя u рaссчumaнные ?rt;ля у?rt;aленных землеmрясенuŭ uз волн, nрuхо?rt;ящuх nо?rt; рaзнымu aзuмуmaмu u у лaмu na?rt;енuя, ?rt;aюm основaнuе ?rt;ля оnре?rt;еленuя мощносmu лumосферы Кaрnam u uх окресmносmеŭ. Тонкaя лumосферa(60–80 км) оnре?rt;еленa ?rt;ля Пaннонско о u Трaнсuльвaнско о бaссеŭнов, a сaмaя мощнaя лumосферa(около 180 км) нaхо?rt;umся nо?rt; Южнымu Кaрnamaмu. В ?rt;ру uх чaсmях Кaрnamско о nоясa мощносmь лumосферы нaрaсmaеm во внешнuх зонaх в нanрaвленuu к ю о-зana?rt;ноŭ окрauне Восmочно-Евроnеŭскоŭ u Мол?rt;aвскоŭ nлamформ. Нaблю?rt;aеmся хорошaя корреляцuя мощносmu лumосферы, оnре?rt;еленноŭ нa основе невязок временu Р волн, с мa нumоmеллурuческuмu ?rt;aннымu nо лубuне слоя nовышенноŭ элекmроnрово?rt;uмосmu в верхнеŭ мaнmuu.

     

Distribution ofBe andBe in the Pacific Ocean

The vertical distributions of¹⁰Be and⁹Be at three locations in the Pacific (25°N, 170°E; 17°N, 118°W; 3°S, 117°W) are presented. The results show that both isotopes exhibit nutrient-like profiles. From the surface to the bottom, the increase for¹⁰Be is two- to threefold and that for⁹Be is about fivefold. While the inter-station variations in surface water concentrations may reach a factor of two, deep-water values tend to be much more uniform averaging about 2000 atoms/g for¹⁰Be and 30 pM for⁹Be. A similar situation applies to the¹⁰Be/⁹Be ratio; it varies approximately from 1 to 3 × 10⁻⁷ (atom/atom) at shallow depths but tends toward a value close to 1.1 × 10⁻⁷ in the deep ocean. The variation of¹⁰Be/⁹Be can be viewed as resulting from the fact that¹⁰Be in a given parcel of water consists of two components: recycled and primary. The recycled component is that part of¹⁰Be which has reached tracer equilibrium with⁹Be, as opposed to the primary component which, upon entering the sea from the atmosphere, has yet to equilibrate with⁹Be through particle cycling and mixing processes. It is estimated that 70% to nearly 100% of¹⁰Be at the three stations are being recycled, and the recycled beryllium bears an atomic ratio of¹⁰Be/⁹Be close to 1 × 10⁻⁷. The oceanic residence time of Be is of the order of 1000–4000 years, comparable to or slightly longer than the ocean mixing time.     

Geoid anomalies across fracture zones and the thickness of the lithosphere

Recent advances in the measurement and interpretation of geoid height anomalies provide a new way to estimate the thickness of the oceanic lithosphere as a function of crustal age. GEOS-III satellite altimetry measurements show abrupt changes in sea level across fracture zones which separate areas of lithosphere with different ages. These changes have the correct location, amplitude, and wavelength to be caused by the combined gravitational attraction of the relief across the fracture zone and the isostatic support of this relief. Eight profiles of geoid height and bathymetry across the Mendocino fracture zone are inverted to determine the depth of the isostatic compensation, assuming that the compensation occurs in a single layer. These depths are then interpreted with a thermal boundary layer model of lithospheric growth. To explain satisfactorily the geoid measurements, the thermal diffusivity of the upper mantle must be 3.3 × 10^?3 cm² s^?1 and the thickness of the lithosphere, defined as the depth at which the geotherm reaches 95% of its maximum value, must be9.1km m.y.^?1/2 × t^1/2, where t is lithospheric age.     

Volcano spacing,fractures, and thickness of the lithosphere

Both “hot-spot” type and possibly island-arc volcanoes may form at the intersections of fractures whose spacing is near the thickness of the lithosphere and increases with increasing thickness. An approximate equality between layer thickness and spacing of major fractures observed in some sedimentary rocks and clay cake models may thus extend to the “mega-joints” that have fractured the lithosphere and controlled volcano spacing on the earth, and possibly on Mars. If the hot-spot fractures are interpreted as due to shear, many hot-spot fracture systems suggest roughly north-south least principal stress, or, alternatively in some instances, a 90° rotation of this pattern.     

The role of elastic stored energy in controlling the long term rheological behaviour of the lithosphere

The traditional definition of lithospheric strength is derived from the differential stresses required to form brittle and ductile structures at a constant strain rate. This definition is based on dissipative brittle and ductile deformation and does not take into account the ability of the lithosphere to store elastic strain. Here we show the important role of elasticity in controlling the long-term behaviour of the lithosphere. This is particularly evident when describing deformation in a thermodynamic framework, which differentiates between stored (Helmholtz free energy) and dissipative (entropy) energy potentials. In our model calculations we stretch a continental lithosphere with a wide range of crustal thickness (30–60 km) and heat flow (50–80 mW/m²) at a constant velocity. We show that the Helmholtz free energy, which in our simple calculation describes the energy stored elastically, converges for all models within a 25% range, while the dissipated energy varies over an order of magnitude. This variation stems from complex patterns in the local strain distributions of the different models, which together operate to minimize the Helmholtz free energy. This energy minimization is a fundamental material behaviour of the lithosphere, which in our simple case is defined by its elastic properties. We conclude from this result that elasticity (more generally Helmholtz free energy) is an important regulator of the long-term geological strength of the lithosphere.