Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

长江中下游成矿带三维S波速度结构及对深部过程的约束

利用73个固定台站记录的163个远震事件数据,采用多道互相关技术挑选了5524条S波到时数据,并对S波到时数据进行地壳校正,在此基础上采用天然地震层析成像方法和远震S波到时信息,获得了长江中下游成矿带上地幔的三维S波速度结构模型.研究结果表明:(1)研究区域上地幔存在着明显的低速异常,且走向与成矿带相同,可能为上涌的软流圈热物质;(2)研究区域地幔过渡带和上地幔底部存在着明显的高速异常,可能为俯冲的古太平洋板块和拆沉的岩石圈;(3)成矿带上地幔的低速异常呈现由南向北逐渐变浅的空间分布特征,该特征表明软流圈热物质由南向北上涌.综合分析认为,成矿带中生代大规模岩浆活动和成矿作用的深部过程主要与岩石圈的拆沉密切相关.     

Plutonism in the Variscan Odenwald (Germany): from subduction to collision

 Latest Devonian to early Carboniferous plutonic rocks from the Odenwald accretionary complex reflect the transition from a subduction to a collisional setting. For ∼362 Ma old gabbroic rocks from the northern tectonometamorphic unit I, initial isotopic compositions (ε_Nd=+3.4 to +3.8;⁸⁷Sr/⁸⁶Sr =0.7035–0.7053;δ¹⁸O=6.8–8.0‰) and chemical signatures (e.g., low Nb/Th, Nb/U, Ce/Pb, Th/U, Rb/Cs) indicate a subduction-related origin by partial melting of a shallow depleted mantle source metasomatized by water-rich, large ion lithophile element-loaded fluids. In the central (unit II) and southern (unit III) Odenwald, syncollisional mafic to felsic granitoids were emplaced in a transtensional setting at approximately 340–335 Ma B.P. Unit II comprises a mafic and a felsic suite that are genetically unrelated. Both suites are intermediate between the medium-K and high-K series and have similar initial Nd and Sr signatures (ε_Nd=0.0 to –2.5;⁸⁷Sr/⁸⁶Sr=0.7044–0.7056) but different oxygen isotopic compositions (δ¹⁸O=7.3–8.7‰ in mafic vs 9.3–9.5‰ in felsic rocks). These characteristics, in conjunction with the chemical signatures, suggest an enriched mantle source for the mafic magmas and a shallow metaluminous crustal source for the felsic magmas. Younger intrusives of unit II have higher Sr/Y, Zr/Y, and Tb/Yb ratios suggesting magma segregation at greater depths. Mafic high-K to shoshonitic intrusives of the southern unit III have initial isotopic compositions (ε_Nd=–1.1 to –1.8;⁸⁷Sr/⁸⁶Sr =0.7054–0.7062;δ¹⁸O=7.2–7.6‰) and chemical characteristics (e.g., high Sr/Y, Zr/Y, Tb/Yb) that are strongly indicative of a deep-seated enriched mantle source. Spatially associated felsic high-K to shoshonitic rocks of unit III may be derived by dehydration melting of garnet-rich metaluminous crustal source rocks or may represent hybrid magmas. Received: 7 December 1998 / Accepted: 27 April 1999     

Boron isotopic composition of subduction-zone metamorphic rocks

Many arc lavas contain material derived from subducted oceanic crust and sediments, but it remains unresolved whether this distinctive geochemical signature is transferred from the subducting slab by aqueous fluids, silicate melts, or both. Boron isotopic measurements have the potential to distinguish between slab transfer mechanisms because ¹¹B fractionates preferentially into aqueous fluids whereas little fractionation may occur during partial melting. Previous studies have shown that δ¹¹B values of island arc lavas (−6 to +7) overlap the range of δ¹¹B values for altered oceanic crust (−5 to +25) and pelagic sediments and turbidites (−7 to +11). Secondary ion mass spectrometry (SIMS) analyses of minerals in subduction-zone metamorphic rocks yield δ¹¹B=−11 to −3 suggesting that slab dehydration reactions significantly lower the δ¹¹B values of subducted oceanic crust and sediments. In order to explain the higher δ¹¹B values reported for arc lavas as compared to subduction-zone metamorphic rocks, the B-bearing component derived from the metamorphosed slab must be enriched in ¹¹B relative to the slab, favoring an aqueous fluid as the slab transfer mechanism.     

Constraints from high-pressure veins in eclogites on the composition of hydrous fluids in subduction zones

Hydrous high-pressure veins formed during dehydration of eclogites in two paleo-subduction zones (Trescolmen locality in the Adula nappe, central Alps and Münchberg Gneiss Massif, Variscan fold belt, Germany) constrain the major and trace element composition of solutes in fluids liberated during dehydration of eclogites. Similar initial isotopic compositions of veins and host eclogites at the time of metamorphism indicate that the fluids were derived predominantly from the host rocks. Quartz, kyanite, paragonite, phengite, zoisite and omphacite are the dominant minerals in the veins. The major element compositions of the veins are in agreement with experimental evidence indicating that the composition of solutes in such fluids is dominated by SiO₂ and Al₂O₃. Relative to N-MORB, the veins show enrichments of Cs, Rb, Ba, Pb, and K, comparable or slightly lower abundances of Sr, U, and Th, and very low abundances of Nd, Sm, Zr, Nb, Ti and Y. The differential fractionation of highly incompatible elements such as K, U and Th in the veins, as well as the presence of hydrous minerals in the eclogites rule out partial melting as a cause for vein formation. These results confirm previous suggestions that fluids derived from subducted basalt may have low abundances of high field strength elements, rare earth elements and Y. Variable vein-eclogite enrichment factors of incompatible alkalis and to a lesser extent Pb appear to reflect mineralogical controls (phengite, epidote-group minerals) on partitioning of these elements during dehydration of eclogite in subduction zones. However, abundance variations of incompatible elements in minerals from eclogites suggest that the composition of fluids released from eclogites at temperatures <700°C may not reflect true equilibrium partitioning during dehydration. Simple models for the trace elements U and Th indicate the relative importance of the basaltic and sedimentary portions of subducted oceanic crust in producing the characteristic chemical signatures of these elements in convergent plate margin volcanism.     

初议成矿地球物理场：献给刘光鼎院士七十华诞,祝愿刘光鼎院士健康长寿

中国大地构造格局和主要地球物理场分布十分吻合,一些主要地球物场的演化促使若干金属元素成矿,地幔热柱,幔隆和幔坳转换带,莫霍斜坡带,俯冲带等和成矿区关系密切。     

松嫩地块东缘和佳木斯地块西缘电性结构

横过松嫩地块东缘和佳木斯地块西缘的大地电磁测深剖面揭示了两块体结合带附近的深部电性结构.本文对剖面测点做了标准化数据处理,并对二维偏离度、构造走向进行了计算和分析,采用非线性共轭梯度(NLCG)二维反演方法对TM模式的数据进行了反演,获得了该剖面的地壳、上地幔电性结构模型,划分出三个典型构造单元:松嫩地块东缘、碰撞拼合带和佳木斯地块西缘.研究结果表明,研究区上地壳基本呈高阻特征,可能为岩浆岩,代表其经历了多期次岩浆作用,而松嫩地块东缘和佳木斯地块西缘的中下地壳的高导体可能与地幔物质的上涌有关;拼合带下方存在西倾的高导体和高阻体,可能是佳木斯地块向西俯冲到松嫩地块下方的构造遗迹;研究区可能发生了拆沉作用,与之伴随的地幔物质上涌可能是后期伸展作用的一个动力.     

Geochemical characteristics of a pre-Middle Jurassic oceanic crust fragment from the Central Pontides in northern Turkey: Geodynamic implications on intra-oceanic subduction initiation

The Central Pontides (northern Turkey) is one of the key localities to understand the geodynamic evolution of the Palaeo- and Neotethyan oceans. It consists of the pre-Jurassic basement units, the Early Jurassic and the Early Cretaceous accretionary complexes, the widespread Middle Jurassic continental arc magmatics and the Late Jurassic to Tertiary cover units. The Early Cretaceous accretionary complex is represented by the Central Pontide Structural Complex and includes the Middle Jurassic oceanic units, which were metamorphosed during the Early Cretaceous. Apart from these oceanic units, a few metaophiolite and serpentinite fragments have been recognized within the basement units, which may represent the remnants of an older ocean. The pre-Middle Jurassic Devrekani Metaophiolite is the largest oceanic fragment and tectonically intercalated within/between the Devrekani Metamorphics and the Çangaldağ Metamorphic Complex. It is mainly composed of harzburgites, dunites with chromite veins and metagabbros, and cut by metabasaltic andesites and metadacites. Petrographically, the gabbro consists mainly of plagioclase and clinopyroxene, and displays phaneritic/porphyritic texture. In contrast, the metabasaltic andesite includes plagioclase and mica phenocrysts within a fine-grained groundmass. Also, the metadacite is composed predominantly of quartz, plagioclase, and mica minerals. Two different magmatic groups belonging to completely different tectono-magmatic settings have been geochemically determined based on the immobile trace element systematics. The metadacites and metabasaltic andesites are akin to continental arc magmatics and characterized by negative Nb and Ta anomalies and depleted HFSE relative to Th and La contents. However, the metagabbro samples display the geochemical signatures of boninitic rocks and characterized by highly depletion in HFSEs and REEs relative to N-MORB. The Devrekani Metaophiolite in the Central Pontides may represent another remnant of pre-Middle Jurassic oceanic crust generation and can be north-eastward continuation of the Permian-aged Almacık complex and the Boğazköy Metaophiolite fragment in the western Sakarya Composite Terrane. It may have been cut by intrusions of the extensive Middle Jurassic continental arc magmatism after its imbrication within the basement unit. The presence of pre-Middle Jurassic oceanic units may indicate that the Paleozoic ocean may have survived as the Jurassic Intra-Pontide Ocean between the Scythian Platform and Sakarya Composite Terrane during the Mesozoic time. Thus, the Intra-Pontide Suture may normally include the Palaeozoic and Mesozoic remnants of the long-lived northward subducting Tethyan ocean.     

Redox evolution of western Tianshan subduction zone and its effect on deep carbon cycle

Knowing the phase relations of carbon-bearing phases at high-pressure(HP) and high-temperature(HT) condition is essential for understanding the deep carbon cycle in the subduction zones.In particular,the phase relation of carbon-bearing phases is also strongly influenced by redox condition of subduction zones,which is poorly explored.Here we summarized the phase relations of carbon-bearing phases(calcite,aragonite,dolomite,magnesite,graphite,hydrocarbon) in HP metamorphic rocks(marble,metapelite,eclogite) from the Western Tianshan subduction zone and high-pressure experiments.During prograde progress of subduction,carbonates in altered oceanic crust change from Ca-carbonate(calcite) to Ca,Mg-carbonate(dolomite),then finally to Mgcarbonate(magnesite) via Mg-Ca cation exchange reaction between silicate and carbonate,while calcite in sedimentary calcareous ooze on oceanic crust directly transfers to high-pressure aragonite in marble or amorphous CaCO3 in subduction zones.Redox evolution also plays a significant effect on the carbon speciation in the Western Tianshan subduction zone.The prograde oxygen fugacity of the Western Tianshan subduction zone was constrained by mineral assemblage of garnet-omphacite from FMQ-1.9 to FMQ-2.5 at its metamorphic peak(maximum P-T) conditions.In comparison with redox conditions of other subduction zones,Western Tianshan has the lowest oxygen fugacity.Graphite and light hydrocarbon inclusions were ubiqutously identified in Western Tianshan HP metamorphic rocks and speculated to be formed from reduction of Fe-carbonate at low redox condition,which is also confirmed by high-pressure experimental simulation.Based on petrological observation and high-pressure simulation,a polarized redox model of reducing slab but oxidizing mantle wedge in subduction zone is proposed,and its effect on deep carbon cycle in subduction zones is further discussed.     

Ancient deep roots for Mesozoic world-class gold deposits in the north China craton:An integrated genetic perspective

The North China Craton(NCC)hosts some of the world-class gold deposits that formed more than 2 billion years after the major orogenic cycles and cratonization.The diverse models for the genesis of these deposits remain equivocal,and mostly focused on the craton margin examples,although synchronous deposits formed in the interior domains.Here we adopt an integrated geological and geophysical perspective to evaluate the possible factors that contributed to the formation of the major gold deposits in the NCC.In the Archean tectonic framework of the NCC,the locations of the major gold deposits fall within or adjacent to greenstone belts or the margins of micro-continents.In the Paleoproterozoic framework,they are markedly aligned along two major collisional sutures-the Trans North China Orogen and the Jiao-Liao-Ji Belt.Since the Mesozoic intrusions hosting these deposits do not carry adequate signals for the source of gold,we explore the deep roots based on available geophysical data.We show that the gold deposits are preferentially distributed above zones of uplifted MOHO and shallow LAB corresponding to thinned crust and eroded sub-lithospheric mantle,and that the mineralization is located above regions of high heat flow representing mantle upwelling.The NCC was at the center of a multi-convergent regime during the Mesozoic which intensely churned the mantle and significantly en riched it.The geophysical data on Moho and LAB upwarp from the centre towards east of the craton is more consistent with paleo-Pacific slab subduction from the east exerting the dominant control on lithospheric thinning.Based on these results,and together with an evaluation of the geochemical and isotopic features of the Mesozoic magmatic intrusions hosting the gold mineralization,we propose a genetic model that invokes reworking of ancient Au archives preserved in the lower crust and metasomatised upper mantle and which were generated through multiple subduction,underplating and cumulation events associated with cratonization of the NCC as well as the subduction-collision of Yangtze Craton with the NCC.The heat and material input along zones of heterogeneously thinned lithosphere from a rising turbulent mantle triggered by Mesozoic convergent margins surrounding the craton aided in reworking the deep roots of the ancient Au reservoirs,leading to the major gold metallogeny along craton margins as well as in the interior of the NCC.