Process and mechanism of mountain-root removal of the Dabie Orogen—Constraints from geochronology and geochemistry of post-collisional igneous rocks

Systematical studies of post-collisional igneous rocks in the Dabie orogen suggest that the thickened mafic lower crust of the oro- gen was partially melted to form low-Mg# adakitic rocks at 143-131 Ma. Delamination and foundering of the thickened mafic lower crust occurred at 130 Ma, which caused the mantle upwelling and following mafic and granitic magmatic intrusions. Mig- matite in the North Dabie zone, coeval with the formation of low-Mg# adakitic intrusions in the Dabie orogen, was formed by partial melting of exhumed ultrahigh-pressure metamorphic rocks at middle crustal level. This paper argues that the partial melting of thickened lower and middle crust before mountain-root collapse needs lithospheric thinning. Based on the geothermal gradient of 6.6~C/km for lithospheric mantle and initial partial melting temperature of ~1000~C for the lower mafic crust, it can be estimated that the thickness of lithospheric mantle beneath thickened lower crust has been thinned to 〈45 km when the thickened lower crust was melting. Thus, a two-stage model for mountain-root removal is proposed. First, the lithospheric mantle keel was partially removal by mantle convection at 145 Ma. Loss of the lower lithosphere would increase heat flow into the base of the crust and would cause middle-lower crustal melting. Second, partial melting of the thickened lower crust has weakened the lower crust and increased its gravity instability, thus triggering delamination and foundering of the thickened mafic lower crust or mountain-root collapse. Therefore, convective removal and delamination of the thickened lower crust as two mechanisms of lithospheric thin- ning are related to causality.     

Mantle compensation of active metamorphic core complexes at Woodlark rift in Papua New Guinea

In many highly extended rifts on the Earth, tectonic removal of the upper crust exhumes mid-crustal rocks, producing metamorphic core complexes. These structures allow the upper continental crust to accommodate tens of kilometres of extension, but it is not clear how the lower crust and underlying mantle respond. Also, despite removal of the upper crust, such core complexes remain both topographically high and in isostatic equilibrium. Because many core complexes in the western United States are underlain by a flat Moho discontinuity, it has been widely assumed that their elevation is supported by flow in the lower crust or by magmatic underplating. These processes should decouple upper-crust extension from that in the mantle. In contrast, here we present seismic observations of metamorphic core complexes of the western Woodlark rift that show the overall crust to be thinned beneath regions of greatest surface extension. These core complexes are actively being exhumed at a rate of 5-10 km Myr(-1), and the thinning of the underlying crust appears to be compensated by mantle rocks of anomalously low density, as indicated by low seismic velocities. We conclude that, at least in this case, the development of metamorphic core complexes and the accommodation of high extension is not purely a crustal phenomenon, but must involve mantle extension.     

Skew of mantle upwelling beneath the East Pacific Rise governs segmentation

Mantle upwelling is essential to the generation of new oceanic crust at mid-ocean ridges, and it is generally assumed that such upwelling is symmetric beneath active ridges. Here, however, we use seismic imaging to show that the isotropic and anisotropic structure of the mantle is rotated beneath the East Pacific Rise. The isotropic structure defines the pattern of magma delivery from the mantle to the crust. We find that the segmentation of the rise crest between transform faults correlates well with the distribution of mantle melt. The azimuth of seismic anisotropy constrains the direction of mantle flow, which is rotated nearly 10 degrees anticlockwise from the plate-spreading direction. The mismatch between the locus of mantle melt delivery and the morphologic ridge axis results in systematic differences between areas of on-axis and off-axis melt supply. We conclude that the skew of asthenospheric upwelling and transport governs segmentation of the East Pacific Rise and variations in the intensity of ridge crest processes.     

Imaging the Indian subcontinent beneath the Himalaya

The rocks of the Indian subcontinent are last seen south of the Ganges before they plunge beneath the Himalaya and the Tibetan plateau. They are next glimpsed in seismic reflection profiles deep beneath southern Tibet, yet the surface seen there has been modified by processes within the Himalaya that have consumed parts of the upper Indian crust and converted them into Himalayan rocks. The geometry of the partly dismantled Indian plate as it passes through the Himalayan process zone has hitherto eluded imaging. Here we report seismic images both of the decollement at the base of the Himalaya and of the Moho (the boundary between crust and mantle) at the base of the Indian crust. A significant finding is that strong seismic anisotropy develops above the decollement in response to shear processes that are taken up as slip in great earthquakes at shallower depths. North of the Himalaya, the lower Indian crust is characterized by a high-velocity region consistent with the formation of eclogite, a high-density material whose presence affects the dynamics of the Tibetan plateau.     

Destruction geodynamics of the North China craton and its Paleoproterozoic plate tectonics

Much attention has been paid in the last two decades to the physical and chemical processes as well as temporal-spatial variations of the lithospheric mantle beneath the North China Craton. In order to provide insights into the geodynamics of this variation, it is necessary to thoroughly study the state and structure of the lithospheric crust and mantle of the North China Craton and its adjacent regions as an integrated unit. Based on the velocity structure of the crust and upper mantle constrained from seismological studies, this paper presents various available geophysical results regarding the lithosphere thickness, the nature of crust-mantle boundary, the upper mantle structure and deformation characteristics as well as their tectonic features and evolution systematics. Combined with the obtained data from petrology and geochemistry, a mantle flow model is proposed for the tectonic evolution of the North China Craton during the Mesozoic-Cenozoic. We suggest that subduction of the Pacific plate made the mantle underneath the eastern Asian continent unstable and able to flow faster. Such a regional mantle flow system would cause an elevation of melt/fluid content in the upper mantle of the North China Craton and the lithospheric softening, which, subsequently resulted in destruction of the North China Craton in different ways of delamination and thermal erosion in Yanshan, Taihang Mountains and the Tan-Lu Fault zone. Multiple lines of evidence recorded in the crust of the North China Craton, such as the amalgamation of the Archean eastern and western blocks, the subduction of Paleo-oceanic crust and Paleo-continental residue, indicate that the Earth in the Paleoproterozoic had already evolved into the plate tectonic system similar to the present plate tectonics.     

Normal faulting in central Tibet since at least 13.5 Myr ago

Tectonic models for the evolution of the Tibetan plateau interpret observed east-west thinning of the upper crust to be the result of either increased potential energy of elevated crust or geodynamic processes that may be unrelated to plateau formation. A key piece of information needed to evaluate these models is the timing of deformation within the plateau. The onset of normal faulting has been estimated to have commenced in southern Tibet between about 14 Myr ago and about 8 Myr ago and, in central Tibet, about 4 Myr ago. Here, however, we report a minimum age of approximately 13.5 Myr for the onset of graben formation in central Tibet, based on mineralization ages determined with Rb-Sr and 40Ar-39Ar data that post-date a major graben-bounding normal fault. These data, along with evidence for prolonged activity of normal faulting in this and other Tibetan grabens, support models that relate normal faulting to processes occurring beneath the plateau. Thinning of the upper crust is most plausibly the result of potential-energy increases resulting from spatially and temporally heterogeneous changes in thermal structure and density distribution within the crust and upper mantle beneath Tibet. This is supported by recent geophysical and geological data, which indicate that spatial heterogeneity exists in both the Tibetan crust and lithospheric mantle.     

Evidence of power-law flow in the Mojave desert mantle

Studies of the Earth's response to large earthquakes can be viewed as large rock deformation experiments in which sudden stress changes induce viscous flow in the lower crust and upper mantle that lead to observable postseismic surface deformation. Laboratory experiments suggest that viscous flow of deforming hot lithospheric rocks is characterized by a power law in which strain rate is proportional to stress raised to a power, n (refs 2, 3). Most geodynamic models of flow in the lower crust and upper mantle, however, resort to newtonian (linear) stress-strain rate relations. Here we show that a power-law model of viscous flow in the mantle with n = 3.5 successfully explains the spatial and temporal evolution of transient surface deformation following the 1992 Landers and 1999 Hector Mine earthquakes in southern California. A power-law rheology implies that viscosity varies spatially with stress causing localization of strain, and varies temporally as stress evolves, rendering newtonian models untenable. Our findings are consistent with laboratory-derived flow law parameters for hot and wet olivine--the most abundant mineral in the upper mantle--and support the contention that, at least beneath the Mojave desert, the upper mantle is weaker than the lower crust.     

Melting in the Earth's deep upper mantle caused by carbon dioxide

The onset of partial melting beneath mid-ocean ridges governs the cycling of highly incompatible elements from the mantle to the crust, the flux of key volatiles (such as CO2, He and Ar) and the rheological properties of the upper mantle. Geophysical observations indicate that melting beneath ridges begins at depths approaching 300 km, but the cause of this melting has remained unclear. Here we determine the solidus of carbonated peridotite from 3 to 10 GPa and demonstrate that melting beneath ridges may occur at depths up to 330 km, producing 0.03-0.3% carbonatite liquid. We argue that these melts promote recrystallization and realignment of the mineral matrix, which may explain the geophysical observations. Extraction of incipient carbonatite melts from deep within the oceanic mantle produces an abundant source of metasomatic fluids and a vast mantle residue depleted in highly incompatible elements and fractionated in key parent-daughter elements. We infer that carbon, helium, argon and highly incompatible heat-producing elements (such as uranium, thorium and potassium) are efficiently scavenged from depths of approximately 200-330 km in the upper mantle.     

Evidence from three-dimensional seismic reflectivity images for enhanced melt supply beneath mid-ocean-ridge discontinuities

Quantifying the melt distribution and crustal structure across ridge-axis discontinuities is essential for understanding the relationship between magmatic, tectonic and petrologic segmentation of mid-ocean-ridge spreading centres. The geometry and continuity of magma bodies beneath features such as overlapping spreading centres can strongly influence the composition of erupted lavas and may give insight into the underlying pattern of mantle flow. Here we present three-dimensional images of seismic reflectivity beneath a mid-ocean ridge to investigate the nature of melt distribution across a ridge-axis discontinuity. Reflectivity slices through the 9 degrees 03' N overlapping spreading centre on East Pacific Rise suggest that it has a robust magma supply, with melt bodies underlying both limbs and ponding of melt beneath large areas of the overlap basin. The geometry of melt distribution beneath this offset is inconsistent with large-scale, crustal redistribution of melt away from centres of upwelling. The complex distribution of melt seems instead to be caused by a combination of vertical melt transport from the underlying mantle and subsequent focusing of melt beneath a magma freezing boundary in the mid-crust.     

Velocity structure of the crust and uppermost mantle in the boundary area of the Tianshan Mountains and the Tarim Basin

A portable 3-component broadband digital seismic array was deployed across the Tianshan orogenic belt (TOB) to investigate the lithospheric structure. Based on receiver function analysis of the teleseismic P-wave data, a 2-D S-wave velocity profile of the boundary area of the TOB and the Tarim Basin was obtained at the depths of 0--80 km.Our results reveal a vertical and lateral inhomogeneity in the crust and uppermost mantle. Four velocity interfaces divide the crystalline crust into the upper, middle and lower crust. A low velocity zone is widely observed in the upper-middle crust. The depth of Moho varies between 42 and 52 km. At the north end of the profile the Moho dips northward with a vertical offset of 4--6 km, which implies a subduction front of the Tarim Basin into the TOB. The Moho generally appears as a velocity transitional zone except beneath two stations in the northern Tarim Basin, where the Moho is characterized by a typical velocity discontinuity. The fine velocity structure and the deep contact deformation of the crust and upper most mantle delineate the north-south lithospheric shortening and thickening in the boundary area of the TOB and the Tarim Basin, which would be helpful to constructing the geodynamical model of the intracontinental mountain-basin-coupling system.     

Arc-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge or eroded fore-arc complexes but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (63-190 mm yr(-1)) and are comparable to the magnitude of subducting Cocos plate motion (approximately 85 mm yr(-1)). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.     

Uppermost mantle structure of the North China Craton: Constraints from interstation Pn travel time difference tomography

The uppermost mantle is the key area for exchange of heat flux and material convection between the crust and lithospheric mantle. Spatial variations of lithospheric thinning and dynamic processes in the North China Craton could inevitably induce the velocity heterogeneity in the uppermost mantle. In this study, we used Pn arrivals from permanent seismic stations in North China and surrounding regions to construct a tomographic image of the North China Craton. The tomographic method with Pn travel time difference data were used to study the velocity variations in the uppermost mantle. Pn velocities in the uppermost mantle varied significantly in the Eastern, Central and Western blocks of the North China Craton. This suggests that the lithosphere beneath different blocks of the North China Craton have experienced distinct tectonic evolutions and dynamic processes since the Paleozoic. The current uppermost mantle has been imprinted by these tectonic and dynamic processes. Fast Pn velocities are prominent beneath the Bohai Bay Basin in the Eastern Block of the North China Craton, suggesting residuals of the Archean lithospheric mantle. Beneath the Tanlu Fault Zone and Bohai Sea, slow Pn velocities are present in the uppermost mantle, which can be attributed to significant lithospheric thinning and asthenospheric upwelling. The newly formed lithospheric mantle beneath Yanshan Mountain may be the dominant reason for the existence of slow Pn velocities in this region. Conversely, the ancient lower crust and lithospheric mantle already have been delaminated. In the Central Block, significant slow Pn velocities are present in Taihangshan Mountain, which also extends northward to the Yinchuan-Hetao Rift on the northern margin of the Ordos Block and Yinshan Orogen. This characteristic probably is a result of hot asthenospheric upwelling along the active tectonic boundary on the margin of the Western Block. The protracted thermal erosion and underplating of hot asthenospheric upwelling may induce lithospheric thinning and significant slow velocities in the uppermost mantle. Fast velocities beneath the Western Block suggest that the thick, cold and refractory Archean lithospheric keel of craton still is retained without apparent destruction.     

A discontinuity in mantle composition beneath the southwest Indian ridge

The composition of mid-ocean-ridge basalt is known to correlate with attributes such as ridge topography and seismic velocity in the underlying mantle, and these correlations have been interpreted to reflect variations in the average extent and mean pressures of melting during mantle upwelling. In this respect, the eastern extremity of the southwest Indian ridge is of special interest, as its mean depth of 4.7 km (ref. 4), high upper-mantle seismic wave velocities and thin oceanic crust of 4-5 km (ref. 6) suggest the presence of unusually cold mantle beneath the region. Here we show that basaltic glasses dredged in this zone, when compared to other sections of the global mid-ocean-ridge system, have higher Na(8.0), Sr and Al2O3 compositions, very low CaO/Al2O3 ratios relative to TiO2 and depleted heavy rare-earth element distributions. This signature cannot simply be ascribed to low-degree melting of a typical mid-ocean-ridge source mantle, as different geochemical indicators of the extent of melting are mutually inconsistent. Instead, we propose that the mantle beneath approximately 1,000 km of the southwest Indian ridge axis has a complex history involving extensive earlier melting events and interaction with partial melts of a more fertile source.     

Typomorphism of the {211}-type zircon

A natural {211}-type zircon has been first discovered in the Pingtan gabbro in the Fujian Province, which is a new type of morphology reported in magmatic zircons. Investigation of morphology and trace-element chemistry shows two distinct growth stages, Stage 1 and Stage 2. The {211}-type zircon is the earliest crystallized phase of Stage 1 zircon. The low Hf contents and extreme depletion in U, Th and Y indicate that Stage 1 zircon was approximately in equilibration with melt during the crystallization, so that they are regarded as crystallized phase in deeper magmatic chamber. The Zr-saturated gabbroic magma is demonstrated to be derived from mantle magmas by differentiation. The magma differentiation and zircon crystallization at the lower crust level indicate the existence of underplating of mantle magma beneath this area, which results in strong magmatic activities in the late Mesozoic.     

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.     

Metastable garnet in oceanic crust at the top of the lower mantle

As oceanic tectonic plates descend into the Earth's lower mantle, garnet (in the basaltic crust) and silicate spinel (in the underlying peridotite layer) each decompose to form silicate perovskite-the 'post-garnet' and 'post-spinel' transformations, respectively. Recent phase equilibrium studies have shown that the post-garnet transformation occurs in the shallow lower mantle in a cold slab, rather than at approximately 800 km depth as earlier studies indicated, with the implication that the subducted basaltic crust is unlikely to become buoyant enough to delaminate as it enters the lower mantle. But here we report results of a kinetic study of the post-garnet transformation, obtained from in situ X-ray observations using sintered diamond anvils, which show that the kinetics of the post-garnet transformation are significantly slower than for the post-spinel transformation. Although metastable spinel quickly breaks down at a temperature of 1,000 K, we estimate that metastable garnet should survive of the order of 10 Myr even at 1,600 K. Accordingly, the expectation of where the subducted oceanic crust would be buoyant spans a much wider depth range at the top of the lower mantle, when transformation kinetics are taken into account.     

Distribution feature of terrestrial heat flow densities in the Bohai Basin,East China

Temperature logging curves at 8 boreholes and well-testing temperature data at 142 boreholes are used to determine geotemperature gradients in the Bohai Basin. The thermal conductivities of 86 rock samples are measured at laboratory and the effects of porosity and temperature are corrected to obtain conductivities in situ. Terrestrial heat flow densities at 76 wells are determined based on these data. The distribution of the heat flow indicates that the terrestrial heat flow in the Bohai Basin is relatively high with an average value of 65.8 mW/m². This characteristic is caused by the tectonic evolution of the basin. During Cenozoic, the lithosphere stretched intermittently and the crust thinned so that heat conducted from the mantle increased and formed thermal abnormity at depth beneath the basin.     

Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling

The Colorado plateau is a large, tectonically intact, physiographic province in the southwestern North American Cordillera that stands at ～1,800-2,000?m elevation and has long been thought to be in isostatic equilibrium. The origin of these high elevations is unclear because unlike the surrounding provinces, which have undergone significant Cretaceous-Palaeogene compressional deformation followed by Neogene extensional deformation, the Colorado plateau is largely internally undeformed. Here we combine new seismic tomography and receiver function images to resolve a vertical high-seismic-velocity anomaly beneath the west-central plateau that extends more than 200?km in depth. The upper surface of this anomaly is seismically defined by a dipping interface extending from the lower crust to depths of 70-90?km. The base of the continental crust above the anomaly has a similar shape, with an elevated Moho. We interpret these seismic structures as a continuing regional, delamination-style foundering of lower crust and continental lithosphere. This implies that Pliocene (2.6-5.3?Myr ago) uplift of the plateau and the magmatism on its margins are intimately tied to continuing deep lithospheric processes. Petrologic and geochemical observations indicate that late Cretaceous-Palaeogene (～90-40?Myr ago) low-angle subduction hydrated and probably weakened much of the Proterozoic tectospheric mantle beneath the Colorado plateau. We suggest that mid-Cenozoic (～35-25?Myr ago) to Recent magmatic infiltration subsequently imparted negative compositional buoyancy to the base and sides of the Colorado plateau upper mantle, triggering downwelling. The patterns of magmatic activity suggest that previous such events have progressively removed the Colorado plateau lithosphere inward from its margins, and have driven uplift. Using Grand Canyon incision rates and Pliocene basaltic volcanism patterns, we suggest that this particular event has been active over the past ～6?Myr.     

Recycling lower continental crust in the North China craton

Foundering of mafic lower continental crust into underlying convecting mantle has been proposed as one means to explain the unusually evolved chemical composition of Earth's continental crust, yet direct evidence of this process has been scarce. Here we report that Late Jurassic high-magnesium andesites, dacites and adakites (siliceous lavas with high strontium and low heavy-rare-earth element and yttrium contents) from the North China craton have chemical and petrographic features consistent with their origin as partial melts of eclogite that subsequently interacted with mantle peridotite. Similar features observed in adakites and some Archaean sodium-rich granitoids of the tonalite-trondhjemite-granodiorite series have been interpreted to result from interaction of slab melts with the mantle wedge. Unlike their arc-related counterparts, however, the Chinese magmas carry inherited Archaean zircons and have neodymium and strontium isotopic compositions overlapping those of eclogite xenoliths derived from the lower crust of the North China craton. Such features cannot be produced by crustal assimilation of slab melts, given the high Mg#, nickel and chromium contents of the lavas. We infer that the Chinese lavas derive from ancient mafic lower crust that foundered into the convecting mantle and subsequently melted and interacted with peridotite. We suggest that lower crustal foundering occurred within the North China craton during the Late Jurassic, and thus provides constraints on the timing of lithosphere removal beneath the North China craton.     

利用接收函数方法研究大盈江断裂两侧S波速度结构

 利用研究区(24.2°~25.2°N,97.5°~98.5°E)内大盈江断裂两侧5个流动数字地震台站记录到的宽频带远震P波波形数据进行接收函数反演,得到台站下方0~100km深度范围内地壳、上地幔S波速度细结构.结果表明:研究区内,以大盈江断裂为界,其西北侧Moho面深度约为38km;东南侧Moho面深度为40~42km.断裂两侧地壳、上地幔S波速度结构存在显著差异,东南侧台站下方地壳和上地幔均存在大范围低速区;西北侧台站下方地壳内存在低速层,而上地幔中无明显低速层.研究区内的S波速度结构存在明显的横向非均匀性.