三维数值模拟研究俯冲区火山的时空活动性

俯冲板块的深部脱水使得上覆地幔含水, 从而降低含水地幔的熔点, 导致上覆地幔部分熔融。 部分熔融的地幔柱一旦喷发到地表就是俯冲带火山, 也形成新的地壳。 相对于周围的地幔来讲, 具有较小密度和黏度的部分熔融地幔的时空活动性就控制着俯冲带火山的时空分布特征。 本文主要回顾近年来运用三维热力学岩石力学模型数值模拟研究与板片脱水相关的俯冲带火山活动的时空分布特性。 结果表明, 部分熔融地幔的有效黏度和密度是影响俯冲板片之上的三维地幔柱横向分布特征的主要因素。 高黏度的部分熔融地幔(10²⁰~10²¹ Pa·s )易于形成近平行于海沟的、 长波长(70~100 km)的、 薄的波状地幔柱; 低黏度(10¹⁸~10¹⁹ Pa·s )的熔融地幔易于形成平行于海沟的, 短波长(30~50 km)的蘑菇状地幔柱和垂直于海沟的山脊状地幔柱。 当部分熔融地幔和周围地幔的密度相差小于50 kg/m³时, 在俯冲板片之上只能形成长波长低幅度(宽50~100 km, 高10~15 km)的地幔山丘。 岩浆产率随着时间的变化反映了火山活动的生命周期性。 板块俯冲速度会影响地幔柱形成的深度和范围大小。 高效率熔融提取会增加新地壳增长总量。 低的板块俯冲速度和低的熔融提取效率会增加上地壳(花岗岩质)和中地壳(英安岩质)化学成分的比例。 数值模拟结果可以很好地解释如日本东北、 新西兰、 南阿拉斯加俯冲区火山的横向分布特征。     

西藏吉定蛇绿岩铂族元素地球化学及其对地幔过程的制约

吉定蛇绿岩地幔橄榄岩ＰＧＥ总量略高于原始地幔，但其玄武岩比大洋中脊玄 武岩中的ＰＧＥ含量高，且Ｐｄ／Ｉｒ比值明显偏低．堆晶岩、岩墙群和玄武岩具明显 Ｐｔ负异 常和 Ｒｈ正异常，形成特殊的“Ｎ”型 ＰＧＥ模式，地幔橄榄岩与壳层岩石具有相似的分布 模式．提出吉定蛇绿岩ＰＧＥ分布与本区地幔较高程度部分熔融有关，岩浆结晶分异作 用使ＰＧＥ分异，导致堆晶岩中ＰＧＥ含量向上变低．吉定蛇绿岩中Ａｕ高是由于后期交 代蚀变的结果； Ｐｔ－Ｐｄ分异说明两者分别受合金相和硫化物相控制；Ｒｈ正异常可能与熔 体中较高的ｆｏ＿２有关．     

活动海岭俯冲与岛弧火山活动的热模拟研究

为解释活动海岭的俯冲会造成岛弧火山活动的间断这一现象,本文采用有限单元法对活动海岭俯冲的热演化过程进行了模拟计算．一般情况下,摩擦剪切生热使岛弧下100km左右深度形成地温反转,俯冲板片海洋地壳内角门岩等含水矿物脱水,释放的水进入其上覆板块,降低了地幔岩石的熔点,使热的地幔楔状体内发生部分熔融,形成岛弧火山活动．高温的活动海岭俯冲时不再出现这种温度反转,俯冲板片在较浅深度达到较高温度而脱水,水进入上覆相对较冷的地幔楔状体不能造成熔融,因此岛弧火山活动会中断．     

年青的九州俯冲带的地震结构和岩浆作用

地震层析成像表明，在九州北部深达俯冲的菲律宾海板块的地幔楔存在低速异常，且扩展到弧前区。我们还用岩石学资料采用数值模拟方法估计地幔楔的流体分布。地震学的数值模拟结果表明，年青的板块俯冲带的熔融和岩浆作用不同于老板块。脱水和熔融发生在年青（且热）板块如九州北部（≤26Ma）之上的板块弧和弧前底下，而流体（水溶液和熔融）主要发生在老板块区的弧后。     

西藏高原地热活动,温泉分布与地球物理场特征

本文根据西藏高原地热活动与温泉分布和各种地球物理场的特征,进行了综合分析与研究,给出了地壳深部要素的统一解释。说明了西藏高原水热分布及其活动原因是由于地壳介质中温度很高,地壳中存在着物质的熔融或部分熔融所致。水热显示是受深部地壳构造与物质运移所制约,其热能是由印度板块与欧亚板块碰撞过程中所产生的热量和岩浆上涌所提供的。     

鄂尔多斯及其周缘地热分布的某些特征

利用鄂尔多斯及其周缘的大地热流、地温梯度和地温分布等资料研究了其展布特征和差异性变化，并讨论了地热分布与地壳深浅构造的关系。认为，鄂尔多斯块体内大地热流、地温梯度值整体偏低，而其周缘地区则呈现出环状高值异常带。该区地震活动主要集中分布在鄂尔多斯周缘地区，即鄂尔多斯周缘存在着高热流、高地温梯度和地震丛集分布相统一的现象。这一现象说明，地热和地震活动之间存在着某种必然的联系。     

洋脊和年轻洋壳俯冲与埃达克岩生成的热模拟

新生代以来 ,环太平洋周边分布的埃达克岩 (Adakite)主要与年轻洋壳俯冲时在 70～ 90km深处的部分熔融有关。利用数值方法 ,模拟了洋壳俯冲的热演化过程并讨论了脱水、熔融对埃达克岩浆活动的影响。结果表明 :仅在活动海岭俯冲前后约 10Ma内 ,年轻的、热的俯冲海洋板片在 75～85km深度范围内 ,温度升高至 82 5～ 10 0 0℃脱水 ,导致年轻洋壳中角闪岩部分熔融 ,形成埃达克岩(Adakite)。而一般洋壳俯冲在 10 0km以下深度才脱水 ,由于脱水区压力较高洋壳自身不能熔融 ,水进入上覆地幔楔状体导致部分熔融 ,形成安山岩 (Andesite     

张裂大陆边缘形成演化的数值模拟

本文在等黏态角落流模型的基础上,建立了上升离散地幔流场,将该流场作用于大陆岩石圈底部,能解释岩石圈减薄、裂解并最终形成海底扩张和张裂大陆边缘等一系列过程.数值模拟的结果表明,岩石圈在上升离散地幔流的作用下发生依赖于深度的伸展减薄,表现为不同深度的拉张因子不同,地表热流显著升高,热扰动引起的均衡调整造成地表沉降,同时热扰动造成岩石圈流变强度尤其是在变形中心处显著减小,脆性变形临界深度变浅,而韧性变形范围扩大.在上升离散地幔流的持续作用下最终导致大陆裂解,岩石圈地幔出露,形成海底扩张和张裂大陆边缘.     

山西断陷带地热分布的某些特征

利用山西断陷带的大地热流、地温梯度和温泉分布的资料，研究了地热分布与重力场、地壳深部构造和地震分布的关系。认为，区域性重力负异常多反映的是沉积盆地或凹陷；而在温泉附近，由于地壳深部高密度的熔融物质沿断裂上涌而形成重力正异常；在同一深度上，地壳和岩石圈薄的断陷盆地内部，其热流和温度都为高值；指出山西断陷带内的构造活动性、壳幔突变带、重磁力等值线密集带、大地热流和地温梯度高值区以及温泉密集分布带与地震活动之间都有着内在的联系。     

地震层析成像—地幔对流新模型的研究

利用地震层析成像数据计算地幔对流模型的方法。假设地幔地震层析成像数据对应的地幔横向不均匀结构是地幔热对流的结果，将地震层析成像数据转化为地幔温度不均匀分布。该方法考虑热流体动力学的3个基本方程，以及热输运方程中的非线性项，直接将上述地震层析成像转化的地幔温度不均匀分布作为内部荷载直接引入基本方程，进而反演计算地幔对流的基本格局。理论结果表明，对流格局不仅依赖地震层析成像数据，而且在很大程度上受地幔动力学框架和热动力参数（如Rayleigh数）以及边界条件和系统响应函数的影响。     

Arc-type and intraplate-type ridge basalts formed at the trench-trench-ridge tripIe junction: Implication for the extensive sub-ridge mantle heterogeneity

Abstract ' In situ basalts' represent the ridge magmatism at and close to the ancient trench-trench-ridge triple junction. Such basalts in the Amami, Mugi, and Setogawa accretionary complexes, Southwest Japan, were described and analysed. The geochemical data show that the ' in situ basalts' include all the types of basalts, ranging from alkali basalts and high-alumina basalts to tholeiites, and the compositions tend towards intermediate and silicic rocks. The data also reveal that the ridge basalts are indistinguishable both from the island arc and intraplate basalts, no affinities with mid-ocean-ridge basalts. The sub-ridge mantle adjacent to the triple junction had a component of sub-arc mantle, and this mantle heterogeneity can be generated by the formation of a slab window.     

Effects of plate and slab viscosities on the geoid

The effects of plate rheology (strong plate interiors and weak plate margins) and stiff subducted lithosphere (slabs) on the geoid and plate motions, considered jointly, are examined with three-dimensional spherical models of mantle flow. Buoyancy forces are based on the internal distribution of subducted lithosphere estimated from the last 160 Ma of subduction history. While the ratio of the lower mantle/upper mantle viscosity has a strong effect on the long-wavelength geoid, as has been shown before, we find that plate rheology is also significant and that its inclusion yields a better geoid model while simultaneously reproducing basic features of observed plate motion. Slab viscosity can strongly affect the geoid, depending on whether a slab is coupled to the surface. In particular, deep, high-viscosity slabs beneath the northern Pacific that are disconnected from the surface as a result of subduction history produce significant long-wavelength geoid highs that differ from the observation. This suggests that slabs in the lower mantle may be not as stiff as predicted from a simple thermally activated rheology, if the slab model is accurate.     

Rheological properties of deep subducted oceanic lithosphere and their geodynamic implications

According to the experimental studies on the rheology of two important mantle rocks (eclogite and harzburgite), the rheological properties of the deep subducted oceanic lithosphere are investigated by assuming a simplified harzburgite type slab model with moderate thickness of basaltic layer. When the mantle convergence rate is small or the subducting slab has been trapped in the mantle for an enough long time, the strength profile of the slab is characterized by a strong subducting crustal component lying on a weak subducting upper mantle. However, if the convergence rate is large enough, the subducting slab will be featured only by a rigid cold center. Our study suggests that the detachment of the subducting crust component from the underlying upper mantle is only likely to happen in hot slow subducting slabs, but not the cold fast subducting lithosphere. Rheological properties of the harzburgitic and the eclogitic upper mantle vary with depths. The eclogitic upper mantle is stronger than the peridotitic upper mantle across the upper mantle. Transition zone is the high strength and high viscosity layer in the upper mantle except the lithosphere.     

Mesozoic mafic magmatism in North China: Implications for thinning and destruction of cratonic lithosphere

The North China Craton (NCC) has been thinned from >200 km to <100 km in its eastern part. The ancient subcontinental lithospheric mantle (SCLM) has been replaced by the juvenile SCLM in the Meoszoic. During this period, the NCC was destructed as indicated by extensive magmatism in the Early Cretaceous. While there is a consensus on the thinning and destruction of cratonic lithosphere in North China, it has been hotly debated about the mechanism of cartonic destruction. This study attempts to provide a resolution to current debates in the view of Mesozoic mafic magmatism in North China. We made a compilation of geochemical data available for Mesozoic mafic igneous rocks in the NCC. The results indicate that these mafic igneous rocks can be categorized into two series, manifesting a dramatic change in the nature of mantle sources at ~121 Ma. Mafic igneous rocks emplaced at this age start to show both oceanic island basalts (OIB)-like trace element distribution patterns and depleted to weakly enriched Sr-Nd isotope compositions. In contrast, mafic igneous rocks emplaced before and after this age exhibit both island arc basalts (IAB)-like trace element distribution patterns and enriched Sr-Nd isotope compositions. This difference indicates a geochemical mutation in the SCLM of North China at ~121 Ma. Although mafic magmatism also took place in the Late Triassic, it was related to exhumation of the deeply subducted South China continental crust because the subduction of Paleo-Pacific slab was not operated at that time. Paleo-Pacific slab started to subduct beneath the eastern margin of Eruasian continent since the Jurrasic. The subducting slab and its overlying SCLM wedge were coupled in the Jurassic, and slab dehydration resulted in hydration and weakening of the cratonic mantle. The mantle sources of ancient IAB-like mafic igneous rocks are a kind of ultramafic metasomatites that were generated by reaction of the cratonic mantle wedge peridotite not only with aqueous solutions derived from dehydration of the subducting Paleo-Pacific oceanic crust in the Jurassic but also with hydrous melts derived from partial melting of the subducting South China continental crust in the Triassic. On the other hand, the mantle sources of juvenile OIB-like mafic igneous rocks are also a kind of ultramafic metasomatites that were generated by reaction of the asthenospheric mantle underneath the North China lithosphere with hydrous felsic melts derived from partial melting of the subducting Paleo-Pacific oceanic crust. The subducting Paleo-Pacific slab became rollback at ~144 Ma. Afterwards the SCLM base was heated by laterally filled asthenospheric mantle, leading to thinning of the hydrated and weakened cratonic mantle. There was extensive bimodal magmatism at 130 to 120 Ma, marking intensive destruction of the cratonic lithosphere. Not only the ultramafic metasomatites in the lower part of the cratonic mantle wedge underwent partial melting to produce mafic igneous rocks showing negative ε_Nd(t) values, depletion in Nb and Ta but enrichment in Pb, but also the lower continent crust overlying the cratonic mantle wedge was heated for extensive felsic magmatism. At the same time, the rollback slab surface was heated by the laterally filled asthenospheric mantle, resulting in partial melting of the previously dehydrated rocks beyond rutile stability on the slab surface. This produce still hydrous felsic melts, which metasomatized the overlying asthenospheric mantle peridotite to generate the ultramafic metasomatites that show positive ε_Nd(t) values, no depletion or even enrichment in Nb and Ta but depletion in Pb. Partial melting of such metasomatites started at ~121 Ma, giving rise to the mafic igneous rocks with juvenile OIB-like geochemical signatures. In this context, the age of ~121 Ma may terminate replacement of the ancient SCLM by the juvenile SCLM in North China. Paleo-Pacific slab was not subducted to the mantle transition zone in the Mesozoic as revealed by modern seismic tomography, and it was subducted at a low angle since the Jurassic, like the subduction of Nazca Plate beneath American continent. This flat subduction would not only chemically metasomatize the cratonic mantle but also physically erode the cratonic mantle. Therefore, the interaction between Paleo-Pacific slab and the cratonic mantle is the first-order geodynamic mechanism for the thinning and destruction of cratonic lithosphere in North China.     

智利三联点南部扩张洋脊俯冲区域岩石层热结构的数值模拟

智利三联点(CTJ)位于纳兹卡板块、南极洲板块与南美板块的交界处,由南极洲—纳兹卡板块之间的智利洋脊俯冲到智利海沟而形成.巴塔哥尼亚板片窗的发展是智利洋脊长期扩张俯冲的结果之一.随着纳兹卡板块的不断东向俯冲,纳兹卡板块范围逐渐变小,CTJ同时向北移动.本文采用数值模拟方法,建立了关于洋脊海沟碰撞的简单二维模型,来研究智利三联点南部扩张洋脊俯冲区域岩石层的热结构.模拟结果表明,洋脊的位置、板块相对汇聚速度及上覆大陆板块的存在均对俯冲区域海洋板块的温度结构有着很大影响,并且大陆板块下方海洋板块温度变化最大的位置距洋脊的水平距离与洋脊到板片窗范围的水平距离两者之间具有较好的一致性.同时,当存在两两板块间的相对汇聚时,洋脊右侧大陆板块下表面的温度升高,俯冲带内海洋板块温度接近于地幔温度.纳兹卡板块以7.8 cm·a~(-1)的速度急速俯冲于南美板块之下的过程中,同时伴随着智利洋脊的持续扩张俯冲,在智利三联点南部,南美板块之下纳兹卡板块的温度因而可以更快地达到地幔软流层的约1300℃温度,并最终消亡于地幔之中.     

Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Mariana island arcs

An increasing number of seismological studies indicate that slabs of subducted lithosphere penetrate the Earth's lower mantle below some island arcs but are deflected, or, rather, laid down, in the transition zone below others. Recent numerical simulations of mantle flow also advocate a hybrid form of mantle convection, with intermittent layering. We present a multi-disciplinary analysis of slab morphology and mantle dynamics in which we account explicitly for the history of subduction below specific island arcs in an attempt to understand what controls lateral variations in slab morphology and penetration depth. Central in our discussion are the Izu-Bonin and Mariana subduction zones. We argue that the differences in the tectonic evolution of these subduction zones—in particular the amount and rate of trench migration—can explain why the slab of subducted oceanic lithosphere seems to be (at least temporarily) stagnant in the Earth's transition zone below the Izu-Bonin arc but penetrates into the lower mantle below the Mariana arc. We briefly speculate on the applicability of our model of the temporal and spatial evolution of slab morphology to other subduction zones. Although further investigation is necessary, our tentative model shows the potential for interpreting seismic images of slab structure by accounting for the plate-tectonic history of the subduction zones involved. We therefore hope that the ideas outlined here will stimulate and direct new research initiatives.     

深部俯冲板片三维构造重建及其几何学、运动学研究——以汤加—克马德克地区俯冲板片为例

西南太平洋板块与澳大利亚板块之间的汤加—克马德克俯冲带,是研究地球动力学最重要的区域之一.本文研究根据MIT-P08地震数据,结合板块构造边界、地震活动分布、海岸地形数据等,基于GOCAD软件平台建立三维地震层析成像,对西南太平洋板块的汤加—克马德克俯冲板片进行三维解释.地震层析成像显示汤加—斐济地区地幔至少存在三个"高速"异常体.早期汤加—克马德克俯冲板片穿过地幔转换带,并进入下地幔,最大深度达1600 km.三维构造模型揭示了汤加—克马德克板片在深度600~800 km处存在断折形变,该俯冲板片去褶皱恢复后,测量其俯冲的最大位移达2600 km.汤加—克马德克板片开始快速俯冲的时间至少在30 Ma之前,平均移动速率约为68~104 mm /a.俯冲板片三维构造重建和恢复,可以有效揭示俯冲板片几何学、运动学,为研究深源地震成因、地球深部变化过程和动力学机制提供约束.     

Thermal Structure of the Ionian Slab

We propose a thermal model of the subducting Ionian microplate. The slab sinks in an isothermal mantle, and for the boundary conditions we take into account the relation between the maximum depth of seismicity and the thermal parameter L_th of the slab, which is a product of the age of the subducted lithosphere and the vertical component of the convergence rate. The surface heat-flux dataset of the Ionian Sea is reviewed, and a convective geotherm is calculated in its undeformed part for a surface heat flux of 42 mW m^–2, an adiabatic gradient of 0.6 mK m^–1, a mantle kinematic viscosity of 10¹⁷ m² s^–1 and an asthenosphere potential temperature of 1300°C. The calculated temperature-depth distribution compared to the mantle melting temperature indicates the decoupling limit between lithosphere and asthenosphere occurs at a depth of 105 km and a temperature of 1260°C. A 70–km thick mechanical boundary layer is found. By considering that the maximum depth of the seismic events within the slab is 600 km, a L_th of 4725 km is inferred. For a subduction rate equal to the spreading rate, the corresponding assimilation and cooling times of the microplate are about 7 and 90 Myr, respectively. The thermal model assumes that the mantle flow above the slab is parallel and equal to the subducting plate velocity of 6 cm yr^–1, and ignores the heat conduction down the slab dip. The critical temperature, above which the subduced lithosphere cannot sustain the stress necessary to produce seismicity, is determined from the thermal conditions governing the rheology of the plate. The minimum potential temperature at the depth of the deepest earthquake in the slab is 730°C.     

Deep melting of old subducted oceanic crust recorded by superchondritic Nb/Ta in modern island arc lavas

Niobium–tantalum systematics of slab-derived melts are powerful tracers that discriminate residual high-pressure rutile-bearing eclogite from low-pressure garnet-bearing amphibolite in subducting plates. Previously reported low Nb–Ta ratios in modern slab melts suggested a predominance of shallow melting in the presence of residual amphibole and that deep melting of rutile-bearing eclogitic slabs, devoid of residual amphibole, is volumetrically insignificant. This study evaluates Nb/Ta in combination with other trace element systematics of modern intra-oceanic and slab melt-related arc lavas from the south-western volcanic chain of the Solomon Islands that cover over 1000 km of the SW Pacific plate border. After a change of subduction polarity, an old subducted Pacific slab and a recently subducting Indian–Australian slab are both present beneath the arc. Solomon arc lavas show sub- to superchondritic Nb–Ta ratios (ca. 10 to 27) which is the largest range ever reported in modern island arc lavas. The large range of Nb/Ta likely results from enrichment of the depleted sub-arc mantle by two distinct slab-derived melts in addition to fluids. One minor slab melt component is derived from the shallow and recent subducting Indian–Australian plate where amphibole is still a significant residual phase. The second slab melt component is predominant in Solomon arc lavas and can be attributed to deep rutile–eclogite-controlled melting of old subducted Jurassic Pacific oceanic crust where residual amphibole is entirely absent or insignificant. The deep Pacific slab melt component is the most likely origin of the extremely high and superchondritic Nb/Ta signatures that produce the upper half of the observed range of Nb/Ta in Solomon arc lavas. The slab melt component that enriched the sub-arc mantle with an unusually high Nb/Ta signature is derived from an initially intact Pacific plate that was probably subject to a slab break-off event and subsequent melting at depths exceeding 100 km. The geochemical evidence presented here shows that old and cold subducted oceanic crust, which is initially not torn, may resist shallow melting but can melt at greater depths instead. The resulting slab melts are generated in the presence of residual rutile-bearing eclogite and significantly fractionate Nb–Ta ratios which may be of relevance at a global scale.