大陆动力学数值模拟:问题、进展与展望

大陆动力学是地球动力学的基本组成部分,是板块构造理论的重要拓展,是固体地球科学的核心命题之一。在系统的地质、地球物理和地球化学观测的基础上,数值模拟是探讨大陆动力学过程和机制的有效手段。本文主要基于大陆动力学数值模拟,围绕大陆形成和演化过程的四个关键科学问题进行总结和探讨。(1)大陆起源与早期地球动力学演化。现今观测到的最早的大陆地壳岩石来自冥古宙,说明地球早期就已经开始大陆地壳的形成;关于当时的构造体制,存在多种不同的模式及其过渡和转换(岩浆洋、热管、滞盖等),该问题的约束甚少,是一个地球动力学的前沿科学问题。(2)大陆岩石圈的稳定性与破坏。大陆岩石圈形成之后经历几十亿年尺度的长时间演化,有些克拉通可以保持大致的稳定性直至现今,而有的克拉通却在显生宙期间经历显著的改造和破坏,无论其稳定存在还是改造破坏,都是值得深入探讨的科学问题。(3)大陆深俯冲与极限折返。超高压岩石折返代表了地球表层物质由浅入深而又由深及浅的物质循环,系统的观测和模拟对该过程和机制已有比较清晰的理解;而近年来观测到的200～350 km的超深折返岩石又提出新的挑战,其模式和机制有待进一步探索。(4)大陆碰撞造山差异...     

Ocean floor and subduction record in the Zermatt‐Saas rodingites,Valtournanche, Western Alps

Multiscale structural analysis and petrological modelling were used to establish the pressure‐peak mineral assemblages and pressure–temperature (P–T) conditions recorded in the rodingites of the upper Valtournanche portion of the oceanic Zermatt‐Saas Zone (ZSZ; Western Alps, northwestern Italy) during Alpine subduction. Rodingites occur in the form of deformed dykes and boudins within the hosting serpentinites. A field structural analysis showed that rodingites and serpentinites record four ductile deformation stages (D1–D4) during the Alpine cycle, with the first three stages associated with new foliations. The most pervasive fabric is S2 that is marked by mineral assemblages in serpentinite indicating pressure‐peak conditions, involving mostly serpentine, clinopyroxene, olivine, Ti‐clinohumite and chlorite. Three rodingite types can be defined: epidote‐bearing, garnet–chlorite–clinopyroxene‐bearing and vesuvianite‐bearing rodingite. In these, the pressure‐peak assemblages coeval with S2 development involve: (i) epidote_II + clinopyroxene_II + Mg‐chlorite_II + garnet_II ± rutile ± tremolite_I in the epidote‐bearing rodingite; (ii) Mg‐chlorite_II + garnet_II clinopyroxene_II ± vesuvianite_II ± ilmenite in the garnet–chlorite–clinopyroxene‐bearing rodingite; (iii) vesuvianite_II + Mg‐chlorite_II + clinopyroxene_II + garnet_II ± rutile ± epidote in vesuvianite‐bearing rodingite. Despite the pervasive structural reworking of the rodingites during Alpine subduction, the mineral relicts of the pre‐Alpine ocean floor history have been preserved and consist of clinopyroxene porphyroclasts (probable igneous relicts from gabbro dykes) and Cr‐rich garnet and vesuvianite (relicts of ocean floor metasomatism). Petrological modelling using thermocalc in the NCFMASHTO system was used to constrain the P–T conditions of the S2 mineral assemblages. The inferred values of 2.3–2.8 GPa and 580–660 °C are consistent with those obtained for syn‐S2 assemblages in the surrounding serpentinites. Multiscale structural analysis indicates that some ocean floor minerals remained stable under eclogite facies conditions suggesting that minerals such as vesuvianite, which is generally regarded as a low‐P phase, could also be stable in favourable chemical systems under high‐P/ultra‐high‐pressure (HP/UHP) conditions. Finally, the reconstructed P–T–d–t path indicates that the P/T ratio characterizing the D2 stage is consistent with cold subduction as estimated in this part of the Alps. The estimated pressure‐peak values are higher than those previously reported in this part of ZSZ, suggesting that the UHP units are larger and/or more abundant than those previously suggested.     

Polygenetic titanite records the composition of metamorphic fluids during the exhumation of ultrahigh‐pressure metagranite in the Sulu orogen

An integrated study of U–Pb ages and trace elements was carried out for titanite and zircon from ultrahigh‐pressure (UHP) metagranites in the Sulu orogen, east‐central China. The results provide constraints on the composition of metamorphic fluids during the exhumation of deeply subducted continental crust. Titanite has two domain types based on REE patterns and trace element variations, Ttn‐I and Ttn‐II respectively. These two domains show indistinguishable U–Pb ages of 232 ± 14 to 220 ± 8 Ma, in general agreement with anatectic zircon U–Pb ages of 223 ± 4 to 219 ± 2 Ma for the partial melting event during early exhumation. The Ttn‐I domains have significantly higher REE, Th, Ta and Sr, and higher Th/U ratios than the Ttn‐II domains, indicating that the two domains have grown from metamorphic fluids with different compositions. For the Ttn‐I domains, Zr‐in‐titanite thermometry yields high temperatures of 773–851 °C at 2.5 GPa, and petrographic observations reveal the presence of melt pseudomorphs. Thus, they are interpreted to have grown from hydrous melts in the early exhumation stage. In contrast, the Ttn‐II domains were texturally equilibrated with amphibolite facies minerals such as biotite and plagioclase and contain inclusions of plagioclase and quartz. The Zr‐in‐titanite thermometry yields lower temperatures of 627–685 °C at 1.0 GPa. In combination with their REE patterns, they are interpreted to have grown from aqueous solutions at amphibolite facies metamorphic conditions during further exhumation. The differences in Th and Sr contents are prominent between the Ttn‐I and Ttn‐II domains, signifying the compositional difference between the hydrous melts and aqueous solutions. Therefore, the polygenetic titanite in the UHP metamorphic rocks provides insights into the geochemical property of metamorphic fluids during the continental subduction‐zone processes.     

Ultrahigh‐pressure and high‐P lawsonite eclogites in Muzhaerte,Chinese western Tianshan

Lawsonite eclogites are crucial to decipher material recycling along a cold geotherm into the deep Earth and orogenic geodynamics at convergent margins. However, their tectono‐metamorphic role and record especially at ultrahigh‐pressure (UHP) conditions are poorly known due to rare exposure in orogenic belts. In a ~4 km long cross‐section in Muzhaerte, China, at the western termination of the HP‐UHP metamorphic belt of western Tianshan, metabasite blocks contain omphacite and lawsonite inclusions in porphyroblastic garnet, although matrix assemblages have been significantly affected by overprinting at shallower structural levels. Two types of lawsonite eclogites occur in different parts of the section and are distinguished based on inclusion assemblages in garnet: Type 1 (UHP) with the peak equilibrium assemblage garnet+omphacite±jadeite+lawsonite+rutile+coesite±chlorite±glaucophane and Type 2 (HP) with the assemblage garnet+omphacite±diopside+lawsonite+titanite+quartz±actinolite±chlorite+glaucophane. Pristine coesite and lawsonite and their pseudomorphs in Type 1 are present in the mantle domains of zoned garnet, indicative of a coesite‐lawsonite eclogite facies. Regardless of grain size and zoning profiles, garnet with Type 1 inclusions systematically shows higher Mg and lower Ca contents than Type 2 (prp_4–25grs_13–24 and prp_1–8grs_20–45 respectively). Phase equilibria modelling indicates that the low‐Ca garnet core and mantle of Type 1 formed at UHP conditions and that there was a major difference in peak pressures (i.e., maximum return depth) between the two types (2.8–3.2 GPa at 480–590°C and 1.3–1.85 GPa at 390–500°C respectively). Scattered exposures of Type 1 lawsonite eclogite is scatteredly exposed in the north of the Muzhaerte section with a structural thickness of ~1 km, whereas Type 2 occurs throughout the rest of the section. We conclude from this regular distribution that they were derived from two contrasting units that formed along two different geothermal systems (150–200°C/GPa for the northern UHP unit and 200–300°C/GPa for the southern HP unit), with subsequent stacking of UHP and HP slices at a kilometre scale.     

Petrogenesis and tectonic implication of lavas from the Yap Trench,western Pacific

We present major and trace element data of lava recovered from the northern Yap Trench in the western Pacific and discuss their petrogenesis and tectonic implications within the framework of interactions between the Caroline Ridge and Yap Trench. Rocks were collected from both landward and seaward trench slopes and exhibited geochemical characteristics similar to backarc basin basalt (BABB) and mid-ocean ridge basalt (MORB), including high Fe content, tholeiitic affinity, high TiO2 value at a given FeOT/MgO ratio, Ti/V ratio between 20 and 50, low Ba/Nb ratio and Th/Nb ratio, and trace element patterns commonly displayed by BABB and MORB, which are distinct from arc lava. These rocks seem to have been generated during mantle upwelling and decompression melting at a spreading center. However, compared with typical forearc lava produced by seafloor spreading in the Mariana forearc region, such as the early Eocene forearc basalts and late Neogene forearc lava in the southernmost Mariana Trench, the Yap Trench lava is derived from a more fertile mantle and feature a more minor subduction component; thus, they cannot be the products of forearc mantle decompression melting. We suggest that the landward slope lava represents backarc basin crust that was overthrust onto the forearc lithosphere during the collision of the Caroline Ridge with the Yap Trench (20–25 Ma), which played a key role in the evolution of the Yap subduction system. Moreover, the seaward slope lava represents the subduction plate crust that accreted onto the deep trench during the collision. This collision event resulted in the cessation of Yap Arc magmatism; thus, the Yap Trench volcanic rocks (<25 Ma) previously suggested to be arc magma products may actually represent the nascent island arc lava with a lower subduction component than in the mature Mariana Arc lava.     

新疆西准噶尔巴尔鲁克蛇绿混杂岩带发现蓝闪片岩

通过1∶25万区调修测工作,确认了巴尔鲁克蛇绿构造混杂岩(带)的存在.在野外调查与研究工作中,发现该蛇绿岩(带)组分较为齐全,在其中还发现了蓝闪片岩,这对研究巴尔鲁克蛇绿构造混杂岩(带)及本地区的地质构造演化具有重要意义.     

Mineral inclusions in diamonds track the evolution of a Mesozoic subducted slab beneath West Gondwanaland

Three major suites of silicate inclusions in sublithospheric diamonds show evidence of formation at depths > 250 km, and for each suite there is evidence of their formation from subducted material. Two of these are the well known basic (majoritic garnet) and ultrabasic (MgSi-perovskite + ferropericlase) suites. The third, the recently recognised Ca-rich suite, is characterised by carbonate, Ca-Si-Ti minerals and some aluminous material. Carbon isotope ratios in the host diamonds and geochemical-petrological features of the inclusions themselves provide evidence for their derivation from subducted lithosphere materials. The diamonds hosting the basic and ultrabasic suites are suggested to form in fluids/melts resulting from the release of water caused by dehydration reactions affecting both the crustal and mantle portions of a subducting slab of ocean lithosphere. Conversely, the diamonds containing the Ca-rich suite are linked with the formation of carbonatitic melts. In the Juina kimberlite province of Brazil, all three suites have been found in close proximity. A model is presented whereby the formation of the suites occurs progressively during the subduction and stagnation of a single lithospheric slab, with all three suites being transported to the lithosphere by a plume with which the carbonatitic melts of the Ca-rich suite are associated. Nd-Sr isotopic data are presented for the Juina majoritic-garnet inclusions, which supports their formation from oceanic crust of Mesozoic age. In conjunction with published age data for a Ca-Si-Ti inclusion, the Juina (Brazil) sublithospheric inclusions document a series of events involving diamond formation during and following the emplacement of a subducted slab between ca 190 and 90 Ma beneath west Gondwanaland. This slab and related subducted slabs dating from the Palaeozoic at the Gondwanan margin may be the source of the widespread DUPAL geochemical anomaly in the South Atlantic and Indian Oceans. The kimberlites bringing the diamonds to the Earth's surface may have arisen from a superplume, developed from a graveyard of former Gondwanan stagnant slabs, at the Core-Mantle-Boundary.     

Application of Organic Geochemistry to Coastal Tar Residues from Central California

Tar residues are common on the coastline of the Monterey Bay National Marine Sanctuary. These coastal tar residues have been washed ashore and usually occur on headlands near the high-tide line. In this study, 18 coastal tar residues were collected and analyzed to determine their carbon isotopic compositions and the values of selected biomarker ratios. All of the residues have very heavy (¹³C-enriched) carbon isotopic compositions spanning a narrow range (δ¹³C = ?22.2 to ?23.4‰), and 28,30-bisnorhopane is present in all samples. These same geochemical characteristics are found in Monterey Formation oils from which the coastal tar residues were likely derived. These coastal residues could result from natural seeps or from accidental spills. Statistically the coastal tar residues can be organized into three groups, each of which may represent different spill or seep events. Seven samples of potential local representative sources for the tar residues were examined, but none could account for the coastal tars.     

Asia: a frontier for a future supercontinent Amasia

Asia is the world’s largest but youngest continent, in which Pacific-type (P-type) and collision-type (C-type) orogenic belts coexist with numerous amalgamated continental blocks. P-type orogens represent major sites of continental growth through tonalite-trondhjemite-granodiorite type (TTG-type) juvenile granitoid magmatism and accretion of oceanic crust and intra-oceanic arcs. The Asian continent includes several P-type orogenic belts, of which the largest are the Central Asian and Western Pacific. The Central Asian Orogenic Belt is dominated by P-type fossil orogens arranged with a regular northward subduction polarity. The Western Pacific is characterized by ongoing P-type orogeny related to the westward subduction of the Pacific plate. Asia has a multi-cratonic structure and its post-Palaeozoic history has witnessed amalgamation of the Laurasia composite continent and Pangaea supercontinent. Nowadays, Asia is surrounded by double-sided subduction zones, which generate new TTG-type crust and supply oceanic crust and microcontinents to its active margins. The TTG-crust can be tectonically eroded and subducted down to the mantle transition zone to form a ‘second’ continent, which may generate mantle upwelling, plumes, and extensive intra-plate volcanism. Moreover, recent plate movements around Asia are dominated by northward directions, which resulted in the India–Eurasia and Arabia–Eurasia collisions beginning at 50–45 and 23–20 Ma, respectively, and will result in Africa–Eurasia collision in the near future. Therefore, Asia is the best candidate to serve as the nucleus for a future supercontinent ‘Amasia’, likely to form 200–250 Ma in the future. In this paper we unravel a puzzle of continental growth in Asia through P-type orogeny by discussing its tectonic history and geological structure, subduction polarity in P-type orogens, tectonic erosion of TTG-type crust and arc subduction at convergent margins, generation of mantle plumes, and prospects of Asia growth and overgrowth.     

Early Cenozoic Tectonics of the Tibetan Plateau

Geological mapping at a scale of 1:250000 coupled with related researches in recent years reveal well Early Cenozoic paleo-tectonic evolution of the Tibetan Plateau. Marine deposits and foraminifera assemblages indicate that the Tethys-Himalaya Ocean and the Southwest Tarim Sea existed in the south and north of the Tibetan Plateau, respectively, in Paleocene-Eocene. The paleooceanic plate between the Indian continental plate and the Lhasa block had been as wide as 900km at beginning of the Cenozoic Era. Late Paleocene transgressions of the paleo-sea led to the formation of paleo-bays in the southern Lhasa block. Northward subduction of the Tethys-Himalaya Oceanic Plate caused magma emplacement and volcanic eruptions of the Linzizong Group in 64.5-44.3 Ma, which formed the Paleocene-Eocene Gangdise Magmatic Arc in the north of Yalung-Zangbu Suture (YZS), accompanied by intensive thrust in the Lhasa, Qiangtang, Hoh Xil and Kunlun blocks. The Paleocene-Eocene depression of basins reached to a depth of 3500-4800 m along major thrust faults and 680-850 m along the boundary normal faults in central Tibetan Plateau, and the Paleocene-Eocene depression of the Tarim and Qaidam basins without evident contractions were only as deep as 300-580 m and 600-830 m, respectively, far away from central Tibetan Plateau. Low elevation plains formed in the southern continental margin of the Tethy-Himalaya Ocean, the central Tibet and the Tarim basin in Paleocene-Early Eocene. The Tibetan Plateau and Himalaya Mts. mainly uplifted after the Indian-Eurasian continental collision in Early-Middle Eocene.