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.     

Receiver functions in northeast China - implications for slab penetration into the lower mantle in northwest Pacific subduction zone

Seismic studies of the subducting lithosphere and the upper mantle discontinuities in the northwest Pacific subduction zone beneath Japan and northeast China have suggested contrary subduction scenarios. There was little consensus on the issue whether the subducting slab penetrates the upper mantle discontinuities into the lower mantle or it is deflected atop of the 660-km discontinuity over several hundred kilometers. We calculate receiver functions from a recent seismic broadband station network located in northeast China and find topographic variations of the upper mantle discontinuities. A deeper-than-normal 660-km discontinuity is observed over an area of 400 km and it coincides with the stagnant slab imaged by seismic tomography. The 660-km discontinuity is locally depressed by more than 35 km and the transition zone is thickened by more than 20 km in the east of the region where it encounters the slab. These observations provide evidence of the slab accumulating in the mantle transition zone and locally penetrating into the lower mantle.     

Dynamics of unstable continental subduction: Insights from numerical modeling

Numerical experiments are used in this study to systematically investigate the effects of convergence rate, crustal rheological strength, and lithospheric thermal structure on the dynamics of continental collision. The study focuses on the types, conditions and processes of unstable continental subduction. Modelling results suggest that the development of unstable continental subduction can be promoted by conditions that tend to decrease rheological strength of the lithosphere, such as low crustal rheological strength, “hot” thermal structure of the lithosphere, or low convergence rate. Unstable subduction mode can be further categorized into three types: (1) multi-stage slab breakoff, (2) continuously “flowing” of fluid-like slab into the upper mantle, and (3) large-scale detachment of the thickened orogenic root. These three types of unstable continental subduction are respectively associated with (1) a low convergence rate, (2) “hot” thermal structure of the lithosphere with a high convergence rate, and (3) moderate-high crustal rheological strength with a low convergence rate. It is also revealed that the evolution of crustal melting is dominated by the deformation pattern of continental collision, which is mainly controlled by crustal rheological strength. The modelling results have important implications for understanding of continental subduction mode selection under specific geodynamic conditions.     

Phase transformations in a harzburgite composition to 26 GPa: implications for dynamical behaviour of the subducting slab

The mineralogy adopted by a depleted harzburgite composition has been studied over the pressure interval 5–26 GPa at temperatures of 1300–1400°C. The pyroxene-garnet component of the harzburgite composition (harzburgite minus 82 wt.% olivine) transforms to majorite garnet by 18–19 GPa, and further disproportionates to the assemblage of garnet + stishovite + Mg₂SiO₄ spinel above 20 GPa. At still higher pressures, first ilmenite (22–24 GPa) and then perovskite MgSiO₃ (24–26 GPa) are found to coexist with garnet. Garnet disappears at 26 GPa and almost complete transition to perovskite is achieved at this pressure. The mineral proportions and density profiles in the subducting oceanic lithosphere, modelled by a combination of 80% harzburgite + 20% primitive MORB compositions are calculated as a function of depth under conditions isothermal with surrounding pyrolite mantle, and also for a temperature distribution in which the slab is substantially cooler than surrounding mantle to below 700 km. Under isothermal conditions, the slab has a density similar to surrounding mantle to a depth of 600 km. However, between 600 and 700 km, the slab is up to 0.08 g/cm³ denser than surrounding mantle. This is caused primarily by the higher alumina content in pyrolite as compared to harzburgite, which causes the transition to perovskite in pyrolite to occur at substantially higher pressures than in harzburgite. The presence of alumina also smears out the garnet-perovskite transition in pyrolite over a depth interval of 50 km, whereas this transformation is much sharper in the harzburgite composition. Calculations based on the observed phase equilibria also show that a subducted cool slab remains much denser (by 0.1–0.3 g/cm³) than surrounding mantle to a depth of 700 km but possesses a density similar to surrounding mantle below this depth. These results have important implications for the dynamical behaviour of slabs possessing different thermal regimes when they encounter the 670 km discontinuity and also for the nature of this discontinuity.     

Dynamics of thinning and destruction of the continental cratonic lithosphere: Numerical modeling

Thinning of the cratonic lithosphere is common in nature, but its destruction is not. In either case, the mechanisms for both thinning and destruction are still widely under debate. In this study, we have made a review on the processes and mechanisms of thinning and destruction of cratonic lithosphere according to previous studies of geological/geophysical observations and numerical simulations, with specific application to the North China Craton (NCC). Two main models are suggested for the thinning and destruction of the NCC, both of which are related to subduction of the oceanic lithosphere. One is the “bottom-up” model, in which the deeply subducting slab perturbs and induces upwelling from the hydrous mantle transition zone (MTZ). The upwelling produces mantle convection and erodes the bottom of the overriding lithosphere by the fluid-melt-peridotite reaction. Mineral compositions and rheological properties of the overriding lithospheric mantle are changed, allowing downward dripping of lithospheric components into the asthenosphere. Consequently, lithospheric thinning or even destruction occurs. The other is the “top-down” model, characterized by the flat subduction of oceanic slab beneath the overriding cratonic lithosphere. Dehydration reactions from the subducting slab would significantly hydrate the lithospheric mantle and decrease its rheological strength. Then the subduction angle may be changed from shallow to steep, inducing lateral upwelling of the asthenosphere. This upwelling would heat and weaken the overriding lithospheric mantle, which led to the weakened lithospheric mantle dripping into the asthenosphere. These two models have some similarities, in that both take the subducting oceanic slab and relevant fluid migration as the major driving mechanism for thinning or destruction of the overriding cratonic lithosphere. The key difference between the two models is the effective depth of the subducting oceanic slab. One is stagnation and flattening in the MTZ, whereas the other is flat subduction at the bottom of the cratonic lithosphere. In the NCC, the eastern lithosphere was likely affected by subduction of the Izanagi slab during the Mesozoic, which would have perturbed the asthenosphere and the MTZ, and induced fluid migration beneath the NCC lithosphere. The upwelling fluid may largely have controlled the reworking of the NCC lithosphere. In order to discuss and analyze these two models further, it is crucial to understand the role of fluids in the subduction zone and the MTZ. Here, we systematically discuss phase transformations of hydrous minerals and the transport processes of water in the subduction system. Furthermore, we analyze possible modes of fluid activity and the problems to explore the applied feasibility of each model. In order to achieve a comprehensive understanding of the mechanisms for thinning and destruction of cratonic lithosphere, we also consider four additional possible dynamic models: extension-induced lithospheric thinning, compression-induced lithospheric thickening and delamination, large-scale mantle convection and thermal erosion, and mantle plume erosion. Compared to the subduction-related models presented here, these four models are primarily controlled by the relatively simple and single process and mechanism (extension, compression, convection, and mantle plume, respectively), which could be the secondary driving mechanisms for the thinning and destruction of lithosphere.     

Material transportation and fluid-melt activity in the subduction channel: Numerical modeling

The subduction channel is defined as a planar to wedge-like area of variable size,internal structure and composition,which forms between the upper and lower plates during slab subduction into the mantle.The materials in the channel may experience complex pressure,temperature,stress and strain evolution,as well as strong fluid and melt activity.A certain amount of these materials may subduct to and later exhume from100 km depth,forming high to ultra-high pressure rocks on the surface as widely discovered in nature.Rock deformation in the channel is strongly assisted by metamorphic fluids activities,which change composition and mechanical properties of rocks and thus affect their subduction and exhumation histories.In this study,we investigate the detailed structure and dynamics of both oceanic and continental subduction channels,by conducting highresolution petrological-thermomechanical numerical simulations taking into account fluid and melt activities.The numerical results demonstrate that subduction channels are composed of a tectonic rock melange formed by crustal rocks detached from the subducting slab and the hydrated mantle rocks scratched from the overriding plate.These rocks may either extrude sub-vertically upward through the mantle wedge to the crust of the upper plate,or exhume along the subduction channel to the surface near the suture zone.Based on our numerical results,we first analyze similarities and differences between oceanic and continental subduction channels.We further compare numerical models with and without fluid and melt activity and demonstrate that this activity results in strong weakening and deformation of overriding lithosphere.Finally,we show that fast convergence of orogens subjected to fluid and melt activity leads to strong deformation of the overriding lithosphere and the topography builds up mainly on the overriding plate.In contrast,slow convergence of such orogens leads to very limited deformation of the overriding lithosphere and the mountain building mainly occurs on the subducting plate.     

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.     

华北克拉通及东邻西太平洋活动大陆边缘地区的P波速度结构：对岩石圈减薄动力学过程的探讨

本文通过地震层析成像研究获得了华北克拉通及其东邻地区（30°N-50°N,95°E -145°E）1°×1°的P波速度扰动图像.结果显示,在西太平洋俯冲带地区,上地幔中西倾的板片状高速异常体与其上方的低速异常区构成俯冲带与上覆地幔楔的典型速度结构式样.俯冲板片高速体在约300~400 km深度范围内被低速物质充填,暗示俯冲板片可能发生了断离.在华北克拉通地区的上地幔中发现三个东倾排列的高速异常带.在此基础上,本文构建了华北克拉通及其东邻西太平洋活动大陆边缘地区的上地幔速度结构模式图,并据此探讨克拉通岩石圈减薄与西太平洋活动大陆边缘的深部动力学联系.本文认为,太平洋板片的俯冲（断离）,触发热地幔物质上涌并在上覆地幔楔中形成对流,使克拉通岩石圈受到改造（底侵与弱化）.随着俯冲板片后撤,地幔楔中的对流场以及对岩石圈改造的影响范围均随之东移,最终导致华北克拉通岩石圈自下而上、从西向东分三个阶段依次拆沉减薄.这一模式能很好地解释现今克拉通岩石圈自西向东呈台阶状减薄的深部现象.     

西太平洋俯冲板块对中国东北构造演化的影响及其动力学意义

中国东北地区在古生代期间以众多微陆块的拼合以及古亚洲洋的闭合为特征,其后又经历了中-新生代太平洋构造域及中生代蒙古—鄂霍茨克构造域的叠加与改造,以致东北地区的构造行迹显得极为复杂,而大兴安岭重力梯级带及其西部地区构造演化是否与西太平洋俯冲有关仍然存在争议.本研究利用分布于中国东北、华北地区以及韩国、日本等部分台网所接收的近震与远震走时数据获得了中国东北地区壳幔精细的三维P波速度结构.成像结果显示,太平洋板块持续西向俯冲,俯冲板片的前缘停滞在大兴安岭—太行山重力梯度带以东区域的地幔转换带之中;长白山火山区上地幔存在着显著的低速异常体,推测西太平洋板块的深俯冲脱水导致了上地幔底部岩石的熔点降低,从而形成了大范围的部分熔融物质上涌.通过分析上地幔的速度结构,我们认为由于太平洋板块的大规模西向深俯冲,在大地幔楔中发生板片脱水、低速热物质上涌等复杂的地球动力学过程;俯冲板片前缘带动上地幔中不均匀分布的地幔流强烈作用于上部的岩石圈,这对东北地区深部壳幔结构乃至大兴安岭重力梯级带的形成、演化有着重要的影响.     

The thermal structure of subduction zones constrained by seismic imaging: Implications for slab dehydration and wedge flow

Thermal models of subduction zones often base their slab–wedge geometry from seismicity at mantle depths and, consequently, cannot be used to evaluate the relationship between seismicity and structure. Here, high-resolution seismic observations from the recent Broadband Experiment Across the Alaska Range (BEAAR) constrain, in a rare instance, the subducting slab geometry and mantle wedge temperature independent of seismicity. Receiver functions reveal that the subducting crust descends less steeply than the Wadati-Benioff Zone. Attenuation tomography of the mantle wedge reveals a high Q and presumably cold region where the slab is less than 80 km deep. To understand these two observations, we generate thermal models that use the improved wedge geometry from receiver functions and that incorporate temperature- and strain-rate-dependent olivine rheology. These calculations show that seismicity within the subducting crust falls in a narrow belt of pressure–temperature conditions, illuminating an effective Clapeyron slope of 0.1 K/MPa at temperatures of 450–750 °C. These conditions typify the breakdown of high-pressure hydrous minerals such as lawsonite and suggest that a single set of dehydration reactions may trigger intermediate-depth seismicity. The models also require that the upper, cold nose of the mantle wedge be isolated from the main flow in the mantle wedge in order to sustain the cold temperatures inferred from the Q tomography. Possibly, sufficient mechanical decoupling occurs at the top of the downgoing slab along a localized shear zone to 80 km depth, considerably deeper than inferred from thrust zone seismicity.     

上地幔俯冲板块的动力学过程:数值模拟

大洋板块俯冲到地幔转换带,进而可形成不同的形态:板块可以停滞在660km不连续面,抑或穿过地幔转换带进入下地幔.这些不同的俯冲模式可进一步影响到海沟的运动.为更好地理解上地幔中俯冲板片的变形行为以及俯冲过程与海沟运动之间的关系,本文通过建立一系列高精度二维热-力学自由俯冲的数值模型,揭示了俯冲板块在上地幔中的变形方式及其与地幔转换带之间的相互作用过程.模拟结果显示,在俯冲板块与地幔转换带的相互作用过程中,其动力学过程可以分为以海沟后撤主导、海沟前进主导以及稳定型海沟等三种主要动力学类型.对于年龄较老,厚度较大的俯冲板块容易形成海沟后撤型俯冲,俯冲板块停滞在660km不连续面.相反,年龄较小,塑性强度较小的板块容易形成海沟前进型俯冲,俯冲板块穿越660km不连续面.     

Continental subduction and exhumation of high-pressure rocks: insights from thermo-mechanical laboratory modelling

Thermo-mechanical physical modelling of continental subduction is performed to investigate the exhumation of deeply subducted continental crust. The model consists of two lithospheric plates made of new temperature sensitive analogue materials. The lithosphere is underlain by liquid asthenosphere. The continental lithosphere contains three layers: the weak sedimentary layer, the crust made of a stronger material, and of a still stronger lithospheric mantle. The whole model is subjected to a constant vertical thermal gradient, causing the strength reduction with depth in each lithospheric layer. Subduction is driven by both push force and pull force. During subduction, the subducting lithosphere is heating and the strength of its layers reduces. The weakening continental crust reaches maximal depth of about 120 km and cannot subduct deeper because its frontal part starts to flow up. The subducted crust undergoes complex deformation, including indicated upward ductile flow of the most deeply subducted portions and localised failure of the subducted upper crust at about 50-km depth. This failure results in the formation of the first crustal slice which rises up between the plates under the buoyancy force. This process is accompanied by the delamination of the crustal and mantle layers of the subducting lithosphere. The delamination front propagates upwards into the interplate zone resulting in the formation of two other crustal slices that also rise up between the plates. Average equivalent exhumation rate of the crustal material during delamination is about 1 cm/year. The crust-asthenosphere boundary near the interplate zone is uplifted. The subducted mantle layer then breaks off, removing the pull force and thereby stopping the delamination and increasing horizontal compression of the lithosphere. The latter produces shortening of the formed orogen and the growth of relief. The modelling reveals an interesting burial/exhumation evolution of the sedimentary cover. During initial stages of continental subduction the sediments of the continental margin are dragged to the overriding plate base and are partially accreted at the deep part of the interplate zone (at 60-70 km-depth). These sediments remain there until the beginning of delamination during which the pressure between the subducted crust and the overriding plate increases. This results in squeezing the underplated sediments out. Part of them is extruded upwards along the interplate zone to about 30-km depth at an equivalent rate of 5-10 cm/year.     

Tomographic structure of East Asia:Ⅱ.Stagnant slab above 660 km discontinuity and its geodynamic implications

P-wave arrival times of both regional and teleseismic earthquakes were inverted to obtain mantle structures of East Asia.No fast(slab) velocity anomalies was not find beneath the 660-km discontinuity through tomography besides a stagnant slab within the transition zone.Slow P-wave velocity anomalies are present at depths of 100-250 km below the active volcanic arc and East Asia.The western end of the flat stagnant slab is about 1 500 km west to active trench and may also be correlated with prominent surface topographic break in eastern China.We suggested that active mantle convection might be operating within this horizontally expanded "mantle wedge" above both the active subducting slabs and the stagnant flat slabs beneath much of the North China plain.Both the widespread Cenozoic volcanism and associated extensional basins in East Asia could be the manifestation of this vigorous upper mantle convection.Cold or thermal anomalies associated with the stagnant slabs above the 660-km discontinuity have not only caused a broad depression of the boundary due to its negative Clapeyron slope but also effectively shielded the asthenosphere and continental lithosphere above from any possible influence of mantle plumes in the lower mantle.     

Numerical simulations of mantle flow around slab edges

A simple kinematic-dynamic model of mantle flow around the slab-edge is constructed in order to understand the flow complexity there. The flow velocity on the top and the small boundary region around the shallow plate boundary is kinematically imposed in order to achieve a subduction-like feature and the flow in other part is dynamically calculated. The geometry of the plate mimics the region around the junction of Aleutian Islands and Kamchatka, that are examples of the convergent-transform fault plate boundaries. In a simple model in which the overlying plate is almost stationary, the lateral flow from the mantle under the subducting slab to the mantle under the neighboring plate is of minor importance, once the slab penetrates into the high viscosity layer where the downward flow encounters the resistance. Similar situation was found when the trench is advancing, that is, the trench moves toward the overlying plate. For the case with retreating trench, that is, the trench moves toward the subducting plate, a lateral flow exists even after the slab penetrates into the high viscosity layer, although its magnitude is significantly smaller than that of the plate velocity. The presence of a low viscosity layer just beneath the subducting plate may promote the emergence of lateral flow. A significant lateral flow is observed when the high temperature anomaly, that is, buoyant and low viscosity block carried by the movement of subducting plate, approaches the slab. These results may have important implications for the possible existence of trench parallel flow in the sub-slab mantle.     

西北太平洋俯冲带大地幔楔及地幔过渡带各向异性证据和特征

地震层析成像研究清晰给出了地球深部俯冲板片的大尺度形态,但与俯冲过程相关的地幔流动特征仍不明确.在俯冲地幔楔系统中,前人观测到了与海沟平行和垂直的快波偏振方向.本文研究了西北太平洋俯冲板片在地幔过渡带中停滞形成的"大地幔楔"中的各向异性特征.对具有长期稳定观测数据的MDJ台站SKS震相和区域深源地震的直达S波震相进行了详细的剪切波分裂研究.采用多波形同时反演方法在保证数据后方位角覆盖范围的基础上得到了更加精确的SKS震相各向异性结果.剪切波分裂参数具有明显的频率依赖性和随初始极化方向变化的特征,表明台站下方存在多个可能的各向异性源区.对不同频段各向异性观测结果的双层各向异性模型拟合表明,上层各向异性参数为φ=32°±12°,δt=1.4±0.4 s,表现出受郯庐断裂带影响的岩石圈中"固化"各向异性特征;下层各向异性参数为φ=-71°±6°,δt=2.4±0.4 s,与海沟近乎垂直.下层较大的各向异性延时表明研究区域的地幔楔和地幔过渡带中存在明显各向异性,推测是由地幔中软流圈流动引起的橄榄石晶格优势取向和停滞的西北太平洋板片引起的地幔水平流动产生的瓦兹利石晶格优势取向共同作用所导致.     

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.     

中国东部海域岩石圈结构面波层析成像

本文通过面波层析成像得到了中国东部海域及邻近地区的地壳上地幔S波速度图像,给出了主要构造单元的区划及其结构特征,并讨论了速度结构与现今构造活动及构造演化历史的关系.研究区内中下地壳的平均速度与地震活动存在比较显著的关系,强震基本都发生在低速区内或高低速过渡区.太行山以东地壳内存在几条北西向低速带,其中张家口—渤海地震带下方的低速带最为显著.东部海域划分成北黄海、南黄海、东海、和冲绳海槽等4个构造块体.北黄海具有较薄较高速的岩石圈,与南华北盆地类似,推测是中生代特提斯洋向北俯冲造成岩石圈减薄的遗迹.北华北地区具有低速的地壳和较厚的岩石圈,岩石圈地幔速度偏低且上下比较均匀,可能反映中生代沿北方缝合带持续碰撞作用的特点.南黄海具有相对较厚的岩石圈,较多地保存了下扬子克拉通的特征.在下扬子与华北地块的拼合过程中,洋壳俯冲可能是北黄海和苏皖地区上地幔低速特征的成因.在125°E以东的朝鲜半岛地区未发现这一拼合过程的遗迹.有可能整个朝鲜半岛都是华北地块的一部分;但也有可能是太平洋俯冲和日本海张开的作用完全改造了朝鲜半岛的岩石圈上地幔,抹去了以往构造运动的痕迹.东海地区的地壳厚度,特别是岩石圈厚度向冲绳海槽方向减小,反映出菲律宾海板块俯冲在弧后广大地区都有影响.冲绳海槽地区可见俯冲的菲律宾海板片以及板片上方显著低速的地壳和上地幔,为冲绳海槽的弧后扩张机制提供了证据.     

Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography

We report a new model of the upper mantle structure beneath Italy obtained by means of P-wave teleseismic tomography. Besides the recent and remarkable development of the Italian Seismic Network, a high model resolution has been achieved improving the inversion method upon the ACH method used in previous investigations and picking high quality arrival times with the Multi-Channel Cross-Correlation technique. The finer details of our Vp model yield new insights into the heterogeneous structure of the Adria continental lithosphere involved in the collision between the Africa and Europe plates. A wide low Vp anomaly located in the northern Adria mantle, facing the Alpine high Vp slab, supports the idea that the Adria lithosphere has been hydrated and thinned during the Alpine subduction. We argue that this mantle softening may have played a key role in favoring the subsequent delamination of the Adria lithosphere in the northern Apennines. We hypothesize that delamination of continental lithosphere previously thinned in a back-arc setting may be considered a key process to favor subduction polarity reversal and recycling of continental material into the mantle circulation. Conversely, in the central-southern Apennines, the velocity structure is consistent with the existence of a deeper oceanic slab that flattens at the base of the upper mantle, in agreement with the widely accepted geodynamic evolution of the central Mediterranean by slab retreat and back-arc spreading. The oceanic slab is discontinuously detached from the surface plate, suggesting a different structure of the Adria lithosphere, which resists subduction instead of favoring delamination.     

中国西部及邻区岩石圈S波速度结构面波层析成像

本文利用瑞利波群速度频散资料和层析成像方法,研究了中国西部及邻近区域(20°N—55°N,65°E—110°E)的岩石圈S波速度结构.结果表明这一地区存在三个以低速地壳/上地幔为特征的构造活动区域:西蒙古高原—贝加尔地区,青藏高原,印支地区.西蒙古高原岩石圈厚度约为80 km,上地幔低速层向下延伸至300 km深度,说明存在源自地幔深部的热流活动.缅甸弧后的上地幔低速层下至200 km深度,显然与印度板块向东俯冲引起俯冲板片上方的热/化学活动有关.青藏高原地壳厚达70 km,边缘地区厚度也在50 km以上并且具有很大的水平变化梯度,与高原平顶陡边的地形特征一致.中下地壳的平均S波速度明显低于正常大陆地壳,在中地壳20~40 km深度范围广泛存在速度逆转的低速层,这一低速层的展布范围与高原的范围相符.这些特征说明青藏高原中下地壳的变形是在印度板块的北向挤压下发生塑性增厚和侧向流动.地幔的速度结构呈现与地壳显著不同的特点.在高原主体和川滇西部地区上地幔顶部存在较大范围的低速,低速区范围随深度迅速减小;100 km以下滇西低速消失,150 km以下基本完全消失.青藏高原上地幔速度结构沿东西方向表现出显著的分段变化.在大约84°E以西的喀喇昆仑—帕米尔—兴都库什地区,印度板块的北向和亚洲板块的南向俯冲造成上地幔显著高速;84°E—94°E之间上地幔顶部速度较低,在大约150~220 km深度范围存在高速板片,有可能是俯冲的印度岩石圈,其前缘到达昆仑—巴颜喀拉之下;在喜马拉雅东构造结以北区域,存在显著的上地幔高速区,可能阻碍上地幔物质的东向运动.川滇西部岩石圈底界深度与扬子克拉通相似,约为180 km,但上地幔顶部速度较低.这些现象表明青藏高原岩石圈地幔的变形/运动方式可能与地壳有本质的区别.