地幔热导率的选取对动力学数值模拟的影响——以岩石圈张裂过程为例

目前存在有多种地幔热导率模型,不同模型在数值和随温压变化的特征上有明显的差异.为探究不同热导率模型对动力学数值模拟结果的影响,本文对不同模型下的岩石圈张裂过程进行模拟研究,探讨地幔热导率对岩石圈热传输、变形和熔融过程的影响及其作用机理.结果显示,不同热导率模型下,岩石圈的变形和熔融特征表现出明显差异.高热导率模型下,岩石圈破裂较晚,形成陆缘较为宽阔,地壳熔融强烈而地幔熔融较弱;低热导率模型下,岩石圈破裂较早,形成陆缘较为狭窄,地幔熔融强烈而地壳熔融较弱.这种差异源于不同地幔热导率下岩石圈和地幔热状态的变化及相应力学性质的改变.高热导率下,热传导的增温效应显著,岩石圈呈现较热的状态,其强度整体较低,壳幔耦合减弱;而低热导率下,热对流的增温效应显著,岩石圈呈较冷的状态,其强度整体较高,壳幔耦合增强.基于模拟结果,本文认为地幔热导率的选取对动力学模拟的结果有着较为显著的影响,相对于随温压的变化,热导率数值的差异对动力学数值模拟的结果影响更大,尤其是对于地幔熔融过程的影响.     

利用剪切波分裂研究青藏高原东北缘及东南缘上地幔各向异性

利用从IRIS上下载的青藏高原东北缘238个台站以及中国地震科学探测台阵喜马拉雅一期350个台站记录到的远震波形数据,通过采用剪切波分裂方法,获取各个台站下方各向异性分裂参数——快波偏振方向（φ）和分裂时差（δt）,从而得到青藏高原东北缘、东南缘上地幔的各向异性特征。研究结果表明,祁连块体、阿拉善块体和鄂尔多斯块体北部,快波方向为NNW-SSE,明显不同于GPS测量得到的近NE-SW的地表位移场方向,延迟时间平均~0.85 s;羌塘块体以及松潘-甘孜褶皱带的西部,快波方向呈现沿顺时针方向旋转的趋势,并且与GPS测量得到的地表位移场方向一致,延迟时间平均为~1.24 s;松潘-甘孜褶皱带东部、秦岭造山带与鄂尔多斯块体南部的交界处,快波方向呈现无序分布,与GPS方向表现出不一致的分布,延迟时间平均~1.08 s;川滇块体北部,快波方向近似N-S方向,与GPS测量得到的地表位移场方向相同,平均延迟时间为~0.925 s;位于北纬26°以南的川滇块体南部,快波方向近E-W方向分布,明显不同于GPS地表位移场方向,平均延迟时间~1.065 s。综合分析推测,羌塘块体、松潘-甘孜褶皱带的西部以及川滇块体北部,地表形变和深部之间的变形是相互耦合的;祁连块体、阿拉善块体和鄂尔多斯块体北部,松潘-甘孜褶皱带东部、秦岭造山带与鄂尔多斯块体南部的交界处以及川滇块体南部,壳幔之间可能存在解耦现象。     

Evolution of lithospheric mantle in the north of Nain‐Baft oceanic crust (Neo‐Tethyan ophiolite of Ashin,Central Iran)

This study is focused on a plagioclase‐bearing spinel lherzolite from Chah Loqeh area in the Neo‐Tethyan Ashin ophiolite. It is exposed along the west of left‐lateral strike‐slip Dorouneh Fault in the northwest of Central‐East Iranian Microcontinent. Mineral chemistry (Mg#^olivine < ~ 90, Cr#^{clinopyroxene} < ~ 0.2, Cr#^spinel < ~ 0.5, Al₂O₃^{orthopyroxene} > ~ 2.5 wt%, Al₂O₃^{clinopyroxene} > ~ 4.5 wt%, Al₂O₃^spinel > ~ 41.5 wt%, Na₂O^{clinopyroxene} > ~ 0.11 wt%, and TiO₂^{clinopyroxene} > ~ 0.04 wt%) confirms Ashin lherzolite was originally a mid‐oceanic ridge peridotite with low degrees of partial melting at spinel‐peridotite facies in a lithospheric mantle level. However, some Ashin lherzolites record mantle upwelling and tectonic exhumation at plagioclase‐peridotite facies during oceanic extension and diapiric motion of mantle along Nain‐Baft suture zone. This mantle upwelling is evidenced by some modifications in the modal composition (i.e. subsolidus recrystallization of plagioclase and olivine between pyroxene and spinel) and mineral chemistry (e.g. increase in TiO₂ and Na₂O of clinopyroxene, and TiO₂ and Cr# of spinel and decrease in Mg# of olivine), as a consequence of decompression during a progressive upwelling of mantle. Previous geochronological and geochemical data and increasing the depth of subsolidus plagioclase formation at plagioclase‐peridotite facies from Nain ophiolite (~ 16 km) to Ashin ophiolite (~ 35 km) suggest a south to north closure for the Nain‐Baft oceanic crust in the northwest of Central‐East Iranian Microcontinent.     

日喀则弧前盆地的埋藏和剥蚀历史——来自低温热年代学的约束

日喀则弧前盆地紧邻印度板块与欧亚大陆碰撞带,研究其剥蚀历史对理解印度板块与欧亚大陆碰撞对造山带剥蚀的影响具有重要意义。文中利用磷灰石裂变径迹(AFT)及锆石和磷灰石的(U-Th)/He(ZHe和AHe)年龄数据,结合已发表的低温热年代数据探讨日喀则弧前盆地的热演化和剥露历史。日喀则弧前盆地磷灰石裂变径迹年龄存在明显的南北差异,南部磷灰石裂变径迹年龄为74～44Ma,对应的剥蚀速率为0. 03～0. 1km/Ma,剥蚀量≤2km;北部磷灰石裂变径迹年龄为27～15Ma,剥蚀速率为0. 09～0. 29km/Ma,但缺失早新生代的热演化历史。而磷灰石的(U-Th)/He年龄表明15Ma BP之后日喀则弧前盆地整体呈现一致的剥露历史。低温热年代数据表明日喀则弧前盆地南部自新生代以来尽管受到印度板块与欧亚大陆碰撞及后期断层活动的影响,海拔由海平面抬升至4. 2km,但一直保持缓慢的剥蚀,表明高原隆升并未直接促使该地区的岩石剥蚀速率加快,这与快速剥蚀即代表造山带开始隆升的假设不相符。此外,日喀则弧前盆地北部的低温热年代学研究表明晚渐新世—早中新世Kailas盆地仅发育于日喀则弧前盆地与冈底斯造山带之间的狭长地带,并在短期内经历了快速的埋藏和剥露。     

Nishinoshima volcano in the Ogasawara Arc: New continent from the ocean?

Nishinoshima, a submarine volcano in the Ogasawara Arc, approximately 1 000 km south of Tokyo, Japan, suddenly erupted in November 2013, after 40 years of dormancy. Olivine‐bearing phenocryst‐poor andesites found in older submarine lavas from the flanks of the volcano have been used to develop a model for the genesis of andesitic lavas from Nishinoshima. In this model, primary andesite magmas originate directly from the mantle as a result of shallow and hydrous melting of plagioclase peridotites. Thus, it only operates beneath Nishinoshima and submarine volcanoes in the Ogasawara Arc and other oceanic arcs, where the crust is thin. The primary magma compositions have changed from basalt, produced at considerable depth, to andesite, produced beneath the existing thinner crust at this location in the arc. This reflects the thermal and mechanical evolution of the mantle wedge and the overlying lithosphere. It is suggested that continental crust‐like andesitic magma builds up beneath submarine volcanoes on thin arc lithosphere today, and has built up beneath such volcanoes in the past. Andesites produced by this shallow and hydrous melting of the mantle could accumulate through collisions of plates to generate continental crust.     

利用远震研究四川盆地及其周缘瑞利面波相速度和方位各向异性

基于四川盆地及周边的245个宽频带台站2010年9月—2014年9月期间的远震记录,提取双台路径瑞利面波相速度频散资料,反演得到四川盆地20～120s的高分辨率瑞利面波相速度及各向异性空间分布.在丰富区域地球物理基础数据的同时,结合已有研究成果对地壳上地幔变形耦合进行探讨,结果表明短周期(20～30s)的相速度分布与四川盆地的地质构造特征相吻合,作为川滇地块、松潘—甘孜地块和四川盆地之间的边界——龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用;松潘—甘孜地块,特别是川滇地块中下地壳普遍表现为明显的低速异常,表明中下地壳相对软弱;而四川盆地的中下地壳整体呈现相对高速,表明四川盆地具有相对坚硬的中下地壳.研究区域东南角接近北扬子地块与南扬子地块的缝合部位,呈现高速异常.四川盆地南部和东南邻区不同周期均具有较强的各向异性,且快波方向较为一致,反映这些地区不同深度变形耦合较好.四川盆地西部、北部及东北部邻区,不同周期的各向异性快波方向变化较大,不同深度变形耦合较差.这些特征与绕喜马拉雅东构造结的物质流动被扬子地块的高速地壳阻挡的宏观认识基本一致.     

Constraining the characters of the upper mantle discontinuities beneath the NE margin of the Tibetan Plateau with a dense broadband seismic array

Based on observations from a dense broadband seismic array located along the northeastern (NE) margin of the Tibetan Plateau in southeastern Gansu Province,we use receiver functions (RFs) to pick the arrival times of P-to-S converted waves and bin the traces in different grids according to the piercing points of the 410 and 660 km discontinuities in the upper mantle.The depths of the two discontinuities are estimated by the ray tracing method with the IASP91 velocity model and a 3-D tomography model.The results indicate the following:(1) The arrival times of the P410s and P660s converted phases are delayed by approximately 1 s than those predicted by the IASP91 model.The mantle transition zone (MTZ) is thicker than that in the global model.(2) The synchronous lags in the P410s and P660s arrival times are consistent with low-velocity anomalies in the upper mantle,which are believed to result mainly from the eastward migration of materials beneath the NE margin of the Tibetan Plateau.(3) Combined with previous tomography results,the depression of the'660'discontinuity and the thickened MTZ are somewhat consistent with the big mantle wedge (BMW) model.However,due to data limitations,more studies are required to explore the BMW in the future.     

A helium stratified and ingassed lower mantle: resolving the helium paradoxes

It is generally believed a variation of 3He/4He isotopic ratios in the mantle is due to only the decay of U and Th,which produces4 He as well as heat.Here we show that not only3He/4He isotopic ratios but also helium contents can be fractionated by thermal diffusion in the lower mantle.The driving force for that fractionation is the adiabatic or convective temperature gradient,which always produces elemental and isotopic fractionation along temperature gradient by thermal diffusion with higher light/heavy isotopic ratio in the hot end.Our theoretical model and calculations indicate that the lower mantle is helium stratified,caused by thermal diffusion due to*400℃temperature contrast across the lower mantle.The highest3He/4He isotopic ratios and lowest He contents are in the lowermost mantle,which is a consequence of thermaldiffusion fractionation rather than the lower mantle is a primordial and undegassed reservoir.Therefore,oceanicisland basalts derived from the deepest lower mantle with high3He/4He isotopic ratios and less He contents—the long-standing helium paradox,is solved by our model.Because vigorous convection in the upper mantle had resulted in disordered or disorganized thermal-diffusion effects in He,Mid-ocean ridge basalts unaffected by mantle plume have a relatively homogenous and lower!3He/4He isotopic compositions.Our model also predicts that 3He/4He isotopic ratios in the deepest lower mantle of early Earth could be even higher than that of Jupiter,the initial He isotopic ratio in our solar system,because the temperature contrast across the lower mantle in the early Earth is the largest and less4 He had been produced by the decay of U and Th.Moreover,the early helium-stratified lower mantle owned the lowest He contents due to over-degassing caused by the largest temperature contrast.Consequently,succeeding evolution of the lower mantle is a He ingassed process due to secular cooling of the deepest mantle.This explains why significant amount of He produced by the decay of U and Th in the lower mantle were not released,another long-standing heat–helium paradox.     

Effect of drainage area on hypsometry from an analysis of small-scale drainage basins in the Siwalik Hills (Central Nepal)

Hypsometry of drainage basins (area–elevation analysis) has generally been used to infer the stage of geomorphic development and to study the influence of varying forcing factors (i.e. tectonics, climate, lithology) on topography. However, the scale dependence of hypsometry has generally been neglected. In order to assess the scale dependence of hypsometry, this study focuses on the sensitivity of hypsometry to different Digital Elevation Model (DEM) resolutions and on the influence of drainage area. Hypsometry inferred from different DEMs is shown to be robust against variations of their resolution. However, hypsometry appears to be dependent on drainage area. We propose that this scale dependence may reflect the varying importance of river and hillslope processes with basin area. Copyright © 1999 John Wiley & Sons, Ltd.     

Modelling aspects of structures isolated with the frictional pendulum system

Different modelling aspects of structures isolated using the frictional pendulum system and subjected to earthquake ground motions are studied herein. Although the vertical dynamics of these structures is given special emphasis, other effects such as large isolator deformations and bidirectional input motion are also considered. Different structural models of the FPS are developed and tested for single-storey structures and a real four-storey building frame; among them, an ‘exact’ formulation of the FPS force–deformation constitutive relationship is presented. Results show that global building responses can be computed within 20 per cent error in the mean using a simplified model that ignores the vertical motion of the building; however, structural member deformations and forces need to be computed using a model that considers such motion. This is of particular importance when there exist correlation between the horizontal and vertical components of ground motion. Further, a physical model of the FPS is introduced and used to determine the response of a real four-storey frame, including uplift and downward impact. Results from this analysis show that local column responses may vary substantially depending on the stiffness of the isolation storey and the presence of a mass at the isolation level. Such mass is capable of filtering the large increase in column shear that results from the impact of the structure after uplift. Uplift occurs at several instants of the response of the structure considered, leading to an increase in column base shear as large as 3 times the shear obtained by ignoring the vertical dynamics of the building. © 1998 John Wiley & Sons, Ltd.