华南陆缘高热流区的壳幔温度结构与动力学背景

华南陆缘是我国重要的矿产、地热资源区.晚中生代以来,在太平洋板块西向俯冲,地幔热对流活动共同作用下,该区出现多期岩浆-热事件和大规模爆发式成矿作用.在前人研究基础上,本文利用地表热流观测资料、地震剪切波资料、重力位球谐系数,计算了壳-幔温度结构,分析了动力学背景.计算结果表明：华南陆缘东南沿海地带,地壳10 km以浅温度达200℃以上,居里点温度475℃,莫霍面平均温度550℃.地壳浅层较热,花岗岩中放射性元素衰变放热是地壳浅层地下水热活动的重要热源,但地壳总体温度不高,为"冷壳热幔"型热结构.地幔中,90 km深度,温度950~1250℃;120 km深度,温度1050~1400℃;150 km深度,温度1200~1450℃;220 km深度,温度1500~1700℃."热"岩石圈底界深度在110~150 km之间,西深东浅.岩石圈内,地幔应力场为挤压-伸展相间格局;岩石圈之下,地幔应力场为一个以南昌为中心、长轴NE-SW向的椭圆.分析认为,晚中生代以来,太平洋板块的西向俯冲,导致华南陆缘在区域性SE向地幔对流背景上叠加局域性不稳定热扰动,在175~85Ma期间,上地幔物质向上流动,形成不同的岩浆活动高峰期.同时,岩石圈地幔受俯冲洋壳流体的影响,含水量高,黏度小,在地幔流切向应力场作用下,岩石圈底界由西向东"波浪"状减薄.现今岩石圈之下仍具备地幔小尺度热对流温度条件,但除地表浅层外,地壳整体温度不高,岩石圈构造稳定.     

海洋板块俯冲作用下上覆大陆岩石层减薄机制的动力学模拟

几乎所有大陆岩石层的减薄现象,可能都与海洋板块的俯冲作用相关,但是两者之间的内在联系迄今仍不十分明确,为此,我们设计了一系列包含洋-陆俯冲系统的二维数值模型,来探讨海洋板块的俯冲作用对上覆大陆岩石层变形行为的影响,尤其对大陆岩石层减薄效应的制约.模型结果表明,海洋板块俯冲过程中的地幔楔熔体对大陆岩石层地幔的热侵蚀以及由熔体上升所诱发的地幔局部对流的强烈扰动会导致上覆大陆岩石层的减薄效应.这种效应不仅表现在横向上的向陆内蔓延,还表现在垂向上的向浅部发展.且多类动力学参数都能制约大陆岩石层的减薄效应.具体地,随着汇聚速率和洋壳厚度的增加,上覆大陆岩石层在横向上的减薄范围越大,在垂向上的减薄程度也越深;而随着俯冲海洋板块年龄的增加,上覆大陆岩石层在横向上的减薄范围增大,但在垂向上的减薄程度会减小;随着上覆大陆岩石层厚度的增加,其横向减薄范围会减小,但在垂向上的减薄程度会加深.本文研究成果能为揭示华北克拉通减薄/破坏的动力学过程提供一定的理论参考依据.     

Thermal-rheological structure of the lithosphere beneath Jiyang Depression: Its implications for geodynamics

Thermal regime of the lithosphere is the scenario of the lithospheric thermal evolution, and the thermo-mechanical state of lithosphere definitively controls its deformation style and mechanism. Better understanding of the lithospheric deep thermal-rheo- logical structure of sedimentary basin will shed light on the formation and evolution dynamic process of the basin. Surface tectonics is the response of the deep structure, and is controlled by the lithospheric ther-mal-rheological properties.…     

中国大陆及邻区岩石圈三维流变结构

依据地震波速得到的上地幔温度和气象台站记录的地表温度为约束,结合地表热流和热导率观测数据,利用有限元方法计算了中国大陆及邻区岩石圈三维热结构.基于此温度结果和GPS观测得到的应变率数据,以滑动摩擦、脆性破裂和蠕变三种强度机制为约束,计算得到了中国大陆及邻区岩石圈三维流变结构.结果显示:弱强度和低等效黏滞性系数的下地壳在中国大陆及邻区普遍存在,并且下地壳的流变强度和等效黏滞性系数比上地壳和岩石圈地幔一般要低1~2个数量级;中国大陆范围内青藏高原存在着厚度最大、强度最低的下地壳;青藏高原的岩石圈强度和等效黏滞性系数比华北、华南和印度板块的都要低;岩石圈流变结构的横向分布特征与重力梯度带和地形过渡带比较一致.     

增厚大陆岩石层热边界层对流剥离的数值模拟

用数值模拟方法模拟了增厚大陆岩石层热边界层被对流地幔剥离并为软流层物质替代的动力学过程．结果表明,在初始温度分层分布、侧向均匀但存在微小热扰动的流场中,80km厚的增厚岩石层热边界层约需60Ma才能被完全剥离,剥离的速率微弱地依赖扰动的强度;在已建立好的流场中,同样厚度的增厚热边界层只需约10Ma就可被剥离．模拟结果暗示青藏高原地壳及岩石层在岩石层增厚和剥离以前就很热,其下伏地馒中可能已存在建立好的上地幔小尺度对流系统,而该尺度的对流系统很可能是由特提斯海洋岩石圈俯冲和消减诱发的     

The characteristics,origins, and geodynamic settings of supergiant gold metallogenic provinces

There are six distinct classes of gold deposits, each represented by metallogenic provinces, having 100's to >1000 tonne gold production. The deposit classes are: (1) orogenic gold; (2) Carlin and Carlin-like gold deposits; (3) epithermal gold-silver deposits; (4) copper-gold porphyry deposits; (5) iron-oxide copper-gold deposits; and (6) gold-rich volcanic hosted massive sulfide (VMS) to sedimentary exhalative (SEDEX) deposits. This classification is based on ore and alteration mineral assemblages; ore and alteration metal budgets; ore fluid pressure(s) and compositions; crustal depth or depth ranges of formation; relationship to structures and/or magmatic intrusions at a variety of scales; and relationship to the P-T-t evolution of the host terrane. These classes reflect distinct geodynamic settings. Orogenic gold deposits are generated at mid-crustal (4–16 km) levels proximal to terrane boundaries, in transpressional subduction-accretion complexes of Cordilleran style orogenic belts; other orogenic gold provinces form inboard, by delamination of mantle lithosphere, or plume impingement. Carlin and Carlin-like gold deposits develop at shallow crustal levels (<4 km) in extensional convergent margin continental arcs or back arcs; some provinces may involve asthenosphere plume impingement on the base of the lithosphere. Epithermal gold and copper-gold porphyry deposits are sited at shallow crustal levels in continental margin or intraoceanic arcs. Iron oxide copper-gold deposits form at mid to shallow crustal levels; they are associated with extensional intracratonic anorogenic magmatism. Proterozoic examples are sited at the transition from thick refractory Archean mantle lithosphere to thinner Proterozoic mantle lithosphere. Gold-rich VMS deposits are hydrothermal accumulations on or near the seafloor in continental or intraoceanic back arcs.

The compressional tectonics of orogenic gold deposits is generated by terrane accretion; high heat flow stems from crustal thickening, delamination of overthickened mantle lithosphere inducing advection of hot asthenosphere, or asthenosphere plume impingement. Ore fluids advect at lithostatic pressures. The extensional settings of Carlin, epithermal, and copper-gold porphyry deposits result from slab rollback driven by negative buoyancy of the subducting plate, and associated induced convection in asthenosphere below the over-riding lithospheric plate. Extension thins the lithosphere, advecting asthenosphere heat, promotes advection of mantle lithosphere and crustal magmas to shallow crustal levels, and enhances hydraulic conductivity. Siting of some copper-gold porphyry deposits is controlled by arc parallel or orthogonal structures that in turn reflect deflections or windows in the slab. Ore fluids in Carlin and epithermal deposits were at near hydrostatic pressures, with unconstrained magmatic fluid input, whereas ore fluids generating porphyry copper-gold deposits were initially magmatic and lithostatic, evolving to hydrostatic pressures. Fertilization of previously depleted sub-arc mantle lithosphere by fluids or melts from the subducting plate, or incompatible element enriched asthenosphere plumes, is likely a factor in generation of these gold deposits. Iron oxide copper-gold deposits involve prior fertilization of Archean mantle lithosphere by incompatible element enriched asthenospheric plume liquids, and subsequent intracontinental anorogenic magmatism driven by decompressional extension from far-field plate forces. Halogen rich mantle lithosphere and crustal magmas likely are the causative intrusions for the deposits, with a deep crustal proximal to shallow crustal distal association. Gold-rich VMS deposits develop in extensional geodynamic settings, where thinned lithosphere extension drives high heat flow and enhanced hydraulic conductivity, as for epithermal deposits. Ore fluids induced hydrostatic convection of modified seawater, with unconstrained magmatic input. Some gold-rich VMS deposits with an epithermal metal budget may be submarine counterparts of terrestrial epithermal gold deposits. Real time analogs for all of these gold deposit classes are known in the geodynamic settings described, excepting iron oxide copper-gold deposits.

     

青藏高原上地幔速度结构及其动力学性质

利用地震层析成像结果分析了中国西部地区的上地幔速度结构,发现青藏高原北部至东南边缘上地幔顶部速度普遍偏低;随着深度的增加,低速区主要分布在羌塘、松潘—甘孜和云南西部地区,而印度大陆、塔里木、柴达木、鄂尔多斯和四川盆地均显示出较高的速度.上述速度分布与青藏高原及周边地区的岩石层结构和深部动力性质密切相关:其中羌塘地区的低速异常反映了青藏北部的地幔上涌和局部熔融,起因于印度大陆岩石层的向北俯冲;松潘—甘孜地区的低速异常与青藏东部的深层物质流动及四川盆地刚性岩石层的阻挡有关;而滇西地区的低速异常可能受到印缅块体向东俯冲作用的影响.以上三个区域构成青藏高原和周边地区的主要地幔异常区.相比之下,印度大陆、塔里木、柴达木、鄂尔多斯和四川盆地的高速异常反映了大陆构造稳定地区的岩石层地幔特点.根据速度变化推测,地幔上涌和韧性变形并非贯穿整个青藏高原,而是主要集中在羌塘、松潘—甘孜和滇西地区,上述构造效应不仅导致岩石层厚度减薄且引发了火山和岩浆活动.     

南海东北部及其邻近地区的Pn波速度结构与各向异性

利用中国地震台网和ISC台站1980～2004年的地震数据，反演了南海东北部及其邻近地区的Pn波速度结构和各向异性.上地幔顶部的速度变化揭示出区域地质构造的深部特征：华南地区速度较高并且变化平缓，具有构造稳定地区的岩石层地幔特征；华南沿海尤其是滨海断裂带附近出现低速异常，表明该断裂可能穿过壳幔边界深达上地幔顶部.南海北部至台湾海峡较高的速度与华南地区类似，反映出大陆边缘和陆架地区的岩石层地幔性质；西沙海槽附近较高的速度不仅反映了华南大陆向南的延伸，而且与海槽裂谷拉张引起的地幔上拱有关，整个南海北部没有发现大规模地幔热流的活动痕迹.相比之下，南海东部次海盆的上地幔顶部存在明显的低速异常，对应于海底扩张中心的地幔上涌区，表明岩石层地幔强烈减薄甚至缺失；台湾东部－吕宋－菲律宾北部的低速异常与地震、火山活动以及岩浆作用紧密相关，揭示了西太平洋岛弧俯冲带的活动特征；南海东北部的洋－陆边界清晰，南海东部和菲律宾海西部较高的速度代表了海洋岩石层地幔的性质.Pn波各向异性反映出区域性构造应力状态及岩石层地幔的变形痕迹：华南地区的各向异性较小，说明这一构造稳定地区的岩石层地幔变形程度较弱；南海北部的快波方向与地壳浅表层构造的伸展方向一致，主要反映了中、新生代以来的大陆边缘张裂和剪切作用对岩石层地幔结构的影响；琉球－台湾－吕宋岛弧两侧各向异性十分强烈，平行于海沟的快波方向表明菲律宾海板块和欧亚大陆的相互作用导致俯冲板块前缘的岩石层地幔强烈变形；台湾东南海域快波方向的变化可能与欧亚大陆和菲律宾海板块俯冲机制的转换以及岩石层被撕裂有关.     

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.     

Softening of the subcontinental lithospheric mantle by asthenosphere melts and the continental extension/oceanic spreading transition

The majority of ophiolitic peridotites in the Alpine–Apennine system show evidence of extensive interaction between subcontinental lithospheric mantle and fractional melts of asthenospheric origin. This interaction led to petrological, structural, and geochemical changes in the lithospheric mantle, and was accompanied by a temperature increase to near-asthenospheric values, resulting in the thermomechanical erosion of the lithosphere. We term the parts of mantle lithosphere thus affected the asthenospherized lithospheric mantle or ALM.The thermal and rheological consequences of thermomechanical erosion are explored by modelling the temperature and rheological properties of the thinned lithosphere as a function of thickness of ALM and time since asthenospherization (i.e., since the beginning of thermal relaxation). Results are given both in terms of rheological profiles (strength envelopes) and total lithospheric strength (TLS) for different lower crustal rheologies. The TLS decreases as a consequence of thermomechanical erosion. This decrease is a non-linear function of the thickness of the ALM. While practically negligible if less than 50% of lithospheric mantle is affected, it becomes significant (up to almost one order of magnitude) if thermomechanical erosion approaches the Moho. The maximum decrease in TLS is achieved within a short time span (1–2 Ma) after the end of the heating episode.As a working hypothesis, we propose that thermomechanical erosion of the lithospheric mantle, related to lithosphere/asthenospheric melts interaction, can be an important factor in a geologically rapid decrease in TLS. This softening could lead to whole lithospheric failure and consequently to a transition from continental extension to oceanic spreading.     

南海深部地球动力学特征及其演化机制

利用地热学、流变学和重力学方法，计算了南海岩石层温度结构、流变特征及地幔对流格局.南海莫霍面温度在600-1000℃之间.岩石层底界面温度在1150-1300℃之间，有效粘滞系数为1020-1021Pa·s，与冰期回弹资料确定的地幔粘度吻合，表明南海深部具备产生地幔热对流的物理条件.研究认为地幔物质由北西向南东方向的运移以及印澳-欧亚板块的碰撞，导致南海北部大陆边缘向洋扩张、离散和断裂解体.在向洋离散过程中，陆-洋岩石层底部地幔局部对流使中央海盆扩张和北部陆缘发生差异性块断运动.     

Seismic multi-arch structures in East China

The deep seismic reflection method has been play-ing an important role in revealing lithospheric struc-tures as well as geodynamic processes within the crustand uppermost mantle[1-3]. After careful analyzing the deep seismic reflection profiles along some of the Mesozoic plutons, the author found that the crustalstructures around these intrusions usually correlate with a special seismic fabric called the seismic multi-arch structure[1]. A seismic multi-arch structure consists of some arch refl…     

Implications of melting for stabilisation of the lithosphere and heat loss in the Archaean

The increased depth and volume of melting induced in a higher temperature Archaean mantle controls the stability of the lithosphere, heat loss rates and the thickness of the oceanic crust. The relationship between density distributions in oceanic lithosphere and the depth of melting at spreading centres is investigated by calculating the mineral proportions and densities of residual mantle depleted by extraction of melt fractions. The density changes related to compositional gradients are comparable to those produced by thermal effects for lithosphere formed from a mantle which is 200°C or more hotter than modern upper mantle. If Archaean continental crust formed initially above oceanic lithosphere, the compositional density gradients may be sufficient to preserve a thick Archaean continental lithosphere within which the Archaean age diamonds are preserved. The amount of heat advected by melts at mid-ocean ridges today is small but heat advected by melting becomes proportionally more important as higher mantle temperatures lead to a greater volume of melt and as the rate of production of oceanic plates increases. Archaean tectonics could have been dominated by spreading rates 2–3 times greater than now and with mantle temperatures between ca. 1600°C and 1800°C at the depth of the solidus. Mid-ocean ridge melting would produce a relatively thick but light refractory lithosphere on which continents could form, protected from copious volcanism and high mantle temperatures.     

Thermal consequences of the formation of a slab window beneath the Mid-Cretaceous southwest Japan arc: A 2-D numerical analysis

Abstract   The development of voluminous granitic magmatism and widespread high-grade metamorphism in Mid-Cretaceous southwest Japan have been explained by the subduction of a spreading ridge (Kula–Pacific or Farallon–Izanagi plate boundaries) beneath the Eurasian continent and the formation of a slab window. In the present study, the thermal consequences of the formation of a slab window beneath a continental margin are evaluated through a 2-D numerical simulation. The model results are evaluated by comparison with the Mid-Cretaceous geology of southwest Japan. Of particular interest are the absence of an amphibolite- to granulite-facies metamorphic belt near the Wadati–Benioff plane, and significant melting of the lower crustal-mafic rocks sufficient to form a large amount of granitic magma. Because none of the model results simultaneously satisfied these two geological interpretations, it is suggested that subduction of plate boundaries in Mid-Cretaceous southwest Japan was not associated with the opening of a slab window. According to previous studies, and the results of the present study, two different tectonic scenarios could reasonably explain the geological interpretations for Mid-Cretaceous southwest Japan: (i) The spreading ridge did not subduct beneath the Eurasian continent, but was located off the continental margin, implying the continuous subduction of very young oceanic lithosphere; (ii) ridge subduction beneath the continental margin occurred after active spreading had ceased. Consequently, in both tectonic scenarios, the subduction of plate boundaries at the Mid-Cretaceous southwest Japan was not associated with a slab window, but very young (hot) oceanic lithosphere.     

Enrichment of the continental lithosphere by OIB melts: Isotopic evidence from the volcanic province of northern Tanzania

Alkali basalts and nephelinites from the southern end of the East African Rift (EAR) in northern Tanzania have incompatible trace element compositions that are similar to those of ocean island basalts (OIB). They define a considerable range of Sr, Nd and Pb isotopic compositions (⁸⁷Sr/⁸⁶Sr= 0.7035−0.7058,ε_Nd = −5to+3, and²⁰⁶Pb/²⁰⁴Pb= 17.5−21.3), each of which partially overlaps the range found in OIB. However, they occupy a unique position in combined Nd, Sr and Pb isotopic compositional space. Nearly all of the lavas have radiogenic Pb, similar to HIMU with high time-integrated²³⁸U/²⁰⁴Pb coupled with unradiogenic Nd (+2 to −5) and radiogenic Sr (>0.704), similar to EMI. This combination has not been observed in OIB and provides evidence that these magmas predominantly acquired their Sr, Nd and Pb in the subcontinental lithospheric mantle rather than in the convecting asthenosphere. These data contrast with compositions for lavas from farther north in the EAR. The Pb isotopic compositions of basalts along the EAR are increasingly radiogenic from north to south, indicating a fundamental change to sources with higher time-integratedU/Pb, closer to the older cratons in the south. An ancient underplated OIB melt component, isolated for about 2 Ga as enriched lithospheric mantle and then remelted, could generate both the trace element and isotopic data measured in the Tanzanian samples. Whereas the radiogenic Pb in Tanzanian lavas requires a source with high time-integratedU/Pb, most continental basalts that are thought to have interacted with the continental lithospheric mantle have unradiogenic Pb, requiring a source with a history of lowU/Pb. Such lowU/Pb is readily accomplished with the addition of subduction-derived components, since the lower averageU/Pb of arc basalts (0.15) relative to OIB (0.36) probably reflects addition of Pb from subducted oceanic crust. If the subcontinental lithosphere is normally characterized by low time-integratedU/Pb it would appear that subduction magmatism is more important than OIB additions in supplying the Pb inventory of the lithospheric mantle. However,U/Pb ratios of xenoliths derived from the continental lithospheric mantle suggest that both processes may be important. This apparent discrepancy could be because xenoliths are not volumetrically representative of the subcontinental lithospheric mantle, or, more likely, that continental lithospheric mantle components in basalts are normally only identified as such when the isotopic ratios are dissimilar from MORB or OIB. Lithospheric enrichment from subaccreted OIB components appears to be more significant than generally recognized.     

上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响

本文以多孔介质中大尺度传热问题为基础,结合热平衡理论分析与数值计算,探讨了上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响.根据大陆岩石圈中孔隙波传热概念模型的初步理论分析结果,指出了采用理论分析和数值模拟相结合的方法在研究大陆岩石圈地幔-地壳热结构模式时的重要性.理论分析方法可用来确定岩石圈尺度范围内大陆岩石圈的厚度和大陆地壳相关的边界条件,从而为地壳范围内数值模型的建立提供一些重要信息.数值模拟方法可以用来模拟地壳尺度范围内地壳的详细结构和复杂几何形状.如果地壳内的热分布是所考虑的主要因素,采用具有地壳尺度的合理数值模型可以有效减少计算机工作量.利用理论分析方法求出的岩石圈尺度范围内大陆岩石圈厚度与地幔传导热流之间关系的理论解,不仅可以用来验证模拟大陆岩石圈内传热问题所采用的数值方法, 而且可以用来初步研究大陆岩石圈内热分布的基本规律,为研究岩石圈地幔热事件中大陆岩石圈热减薄过程提供相应的边界条件.本文从理论分析的观点初步探讨了中国大陆不同构造背景下大陆岩石圈的热结构模式,其结果与从地球物理和地质资料中获得的大陆岩石圈热结构模式十分吻合.研究结果表明由大陆岩石圈中孔隙波传播所导致的上通流是影响大陆岩石圈地幔-地壳热结构模式及大陆岩石圈地幔与地壳之间物质和能量交换的可能机制之一.     

Experimental models of extension of continental lithosphere weakened by percolation of asthenospheric melts

Analogue models are used to investigate extension of a continental lithosphere weakened by asthenospheric melts percolating through the upper mantle, a process that has been hypothesised to control the opening of the Ligurian Tethys. Models were performed in a centrifuge apparatus and reproduced, by using materials such as sand and viscous mixtures, extension of 60-km thick, three-layer continental lithosphere floating above the asthenosphere. The percolated lithospheric mantle was assumed to be characterised by a rheological behaviour similar to that of the asthenosphere. Two sets of experiments investigated the influence on deformation of (1) the thickness of the percolated mantle and the associated strength contrast between the normal and weakened lithosphere, and (2) the lateral width of the weakened zone. Model results suggest that mantle percolation by asthenospheric melts is able to promote strong localised thinning of the continental lithosphere, provided that a significant thickness of the lithospheric mantle is weakened by migrating melts within a narrow region. Strain localisation is maximised for percolation of the whole lithospheric mantle and strong strength contrast between the normal and weakened lithosphere. Under these conditions, the thickness of the lithosphere may be reduced to less than 12 km in 3 Ma of extension. Conversely, localised thinning is strongly reduced if the thickness of the percolated zone is ≤1/3 of the thickness of the whole lithospheric mantle and/or the lithosphere is weakened over wide regions. Overall, model results support the working hypothesis that mantle percolation by asthenospheric melts is a controlling factor in the transition from distributed continental deformation to localised oceanic spreading.     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

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.     

Plume-lithosphere interactions in the ocean basins: constraints from the source mineralogy

Trace element relationships of near-primary alkalic lavas from La Grille volcano, Grande Comore, in the Indian Ocean, as well as those of the Honolulu volcanic series, Oahu, Hawaii, show that their sources contain amphibole and/or phlogopite. Small amounts of each mineral (2% amphibole in the source of La Grille and 0.5% phlogopite plus some amphibole in the source of the Honolulu volcanics) and a range of absolute degrees of partial melting from 1 to 5% for both series are consistent with the observed trace element variation. Amphibole and phlogopite are not stable at the temperatures of convecting upper mantle or upwelling thermal plumes from the deep mantle; however, they are stable at pressure-temperature conditions of the oceanic lithospheric mantle. Therefore, the presence of amphibole and/or phlogopite in the magma source region of volcanics is strong evidence for lithospheric melting, and we conclude that the La Grille and the Honolulu series formed by melting of the oceanic lithospheric mantle.

The identification of amphibole ± phlogopite in the source region of both series implies that the metasomatism by fluids or volatile-rich melts occurred prior to melting. The presence of hydrous phases results in a lower solidus temperature of the lithospheric mantle, which can be reached by conductive heating by the thermal plumes. Isotope ratios of the La Grille and the Honolulu series display a restricted range in composition and represent compositional end-members for each island. Larger isotopic variations in shield lavas, represented by the contemporaneous Karthala volcano on Grande Comore and the older Koolau series on Oahu, reflect interaction of the upwelling thermal plumes with the lithospheric mantle rather than the heterogeneity of deep-seated mantle plume sources or entrainment of mantle material in the rising plume. Literature OsSr isotope ratio covariations constrain the process of plume-lithosphere interaction as occurring through mixing of plume melts and low-degree melts from the metasomatized oceanic lithospheric mantle.

The characterization of the lithospheric mantle signature allows the isotopic composition of the deep mantle plume components to be identified, and mixing relationships show that the Karthala and Koolau plume end-members have nearly uniform isotopic compositions. Based on independent arguments, isotopic variations on Heard and Easter islands have been shown to be a result of mixing between deep plume sources having distinct isotopic compositions with lithosphere or shallow asthenospheric mantle. To the extent that these case studies are representative of oceanic island volcanism, they indicate that interaction with oceanic lithospheric mantle plays an important role in the compositions of lavas erupted during the shield-building stage of plume magmatism, and that isotopic compositions of deep mantle plume sources are nearly uniform on the scale that they are sampled by melting.