中国东部岩石圈热状态与流变学强度特征

根据均衡原理制约的地热计算得到中国东部岩石圈的温度分布状态,以40、70、100km和莫霍面深度等温线图以及600°C、1100°C等温面深度的形式表示.同时计算了以1350°C等温面深度表示的中国东部的热岩石圈厚度.结果显示:在扬子克拉通西部四川盆地之下存在160～200km厚的岩石圈根,但在整个华北克拉通之下缺失岩...     

西藏及其周围地区地壳、地幔地震层析成像——印度板块大规模俯冲于西藏高原之下的证据

文中介绍有关西藏—喜马拉雅碰撞带的一项地震层析成像研究。根据一个用天然地震数据产生的全球波速模型 ,印度板块有可能以近水平状俯冲于整个西藏高原之下至 16 5～ 2 6 0km深度。西藏岩石圈具有低波速地壳和高波速下岩石圈 (75～ 12 0km深 )。在 12 0～ 16 5km深度范围 ,西藏岩石圈与俯冲的印度板块之间有一层低速软流圈物质。高原中部从地表到 310km深处有一低速体 ,说明地幔物质有可能穿过俯冲板块的脆弱部位上隆。这些结果以及野外实测的地壳缩短值说明高原的抬升得助于印度板块的近水平俯冲。我们推论俯冲印度板块的升温上浮以及上覆软流层的存在是造成西藏高原高海拔抬升以及内部地表仍相对平坦的主要原因。2 0 0 1年 1月 2 6日在印度西部发生的毁灭性大地震有可能是俯冲应力在印度板块后缘薄弱处引发的岩石圈大断裂。     

长江中下游成矿带及邻区地壳结构——MASH成矿过程的P波接收函数成像证据？

利用长江中下游成矿带多学科深部探测剖面于2009年11月至2011年3月间采集的天然地震数据,通过天然地震接收函数成像等分析研究,得到了研究区地壳和上地幔结构的清晰图像。接收函数成像结果显示研究区内Moho面深度存在着明显的起伏变化,在长江中下游成矿带(指剖面穿过的长江中下游成矿带宁芜矿集区,下同)下方存在着"幔隆构造"。在剖面东南端(即扬子克拉通北缘),Moho面相对稳定,深度约为30km;在茅山和江南断裂附近,Moho面存在上下起伏现象;在剖面中部或宁芜矿集区下方,Moho面存在明显隆起,深度只有28km;在郯庐断裂带下方,Moho面明显加深,深度达到36km;进一步向北到华北地台南缘,Moho面深度逐渐恢复到了32km左右的平均深度水平。其次,我们在接收函数成像结果中发现,长江中下游成矿带与其周边下地壳结构存在着明显的差异,成矿带的下地壳具有显著的地震波方位各向异性。扬子克拉通北缘的下地壳呈高速的近水平状结构,地震波各向异性特征不明显;与此相比,长江中下游成矿带的下地壳虽然也呈近水平状结构特征,但是,对于沿成矿带走向方向传播的地震波,其下地壳具有高速特征,而对于垂直于成矿带走向方向上传播的地震波,其下地壳却又表现为低速特征,这意味着成矿带的下地壳存在着平行于成矿带走向(即近北东—南西)方向的地震波各向异性,我们解释其是下地壳熔融并沿成矿带走向水平流动导致矿物晶体定向排列的结果。最后,在郯庐断裂以西的华北地台南缘观测到一条从上地壳延伸到中下地壳的南南东向倾斜的转换震相,我们推测它可能是合肥盆地内地壳伸展构造的反映。此外,我们发现接收函数成像结果中观测到的"幔隆构造"与远震P波层析成像结果在成矿带下方150km深度上显示的上地幔低速异常(江国明等,另文发表)存在着良好的对应关系,我们解释它们是软流圈物质上涌的遗迹。综合天然地震接收函数成像、远震P波层析成像和前人关于岩浆岩等方面的研究成果,我们认为长江中下游成矿带现今的下地壳可能是中生代发生成矿作用的多级岩浆房系统的一部分,成矿带的形成可能是类似MASH过程的产物。首先,软流圈物质上涌导致了长江中下游成矿带及其周边拉张环境的形成,在其上部地壳中形成了一系列伸展构造;然后,软流圈物质通过底侵进入长江中下游成矿带的原下地壳并与原下地壳物质发生同化作用,形成类埃达克质岩浆;接着,类埃达克质岩浆沿着伸展、拆离构造上升到地壳浅部形成不同层次的岩浆房和侵入岩体,并与围岩作用形成矿床。     

华北克拉通中生代破坏前的岩石圈地幔与下地壳

华北克拉通是世界上最古老的克拉通之一,有 38亿年的古老陆壳存在,它经历了复杂的地质变迁,在太古宙末(约2500Ma)基本完成克拉通化,在古元古代(约1900~1850Ma)整体受到了高级变质作用,最终完成了克拉通化。它的东部在中生代发生了重大的构造机制的转变,克拉通基底发生了破坏、置换和再造。在太行山重力梯度带以西的华北克拉通受中生代构造转折的改造程度较低,它们的下地壳和岩石圈地幔结构,大致保持了华北克拉通破坏前的状态。前寒武纪麻粒岩地体代表了掀翻抬升到地表的古元古代下地壳,出露地表的时间大致在1850~1800Ma。中、新生代火山岩中的地幔和麻粒岩捕虏体代表了现代的岩石圈地幔和下地壳的岩石。岩石学、地球化学和地球物理的研究,推测华北克拉通西部的岩石圈厚约200km,地壳厚度约45km~50km,是在古元古代(约1.9Ga)时期终极克拉通化作用形成的,其厚度和结构与全球典型的元古宙克拉通岩石圈相同。而太行山重力梯度带以东的克拉通岩石圈地幔受到程度不等的交代、改造、置换和减薄,下地壳大规模重熔,地壳厚度也发生减薄,指示了强烈的壳幔解耦、物质交换和重新耦合的过程。     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

Crust-mantle interaction during the tectono-thermal reactivation of the North China Craton: constraints from SHRIMP zircon U–Pb chronology and geochemistry of Mesozoic plutons from western Shandong

Chronological, geochemical and Sr–Nd–Pb isotopic analyses have been carried out on the Mesozoic plutons in western Shandong with the aim of characterizing crustal–mantle evolution during the tectono-thermal reactivation of the craton. Detailed SHRIMP zircon U–Pb dating reveals two main periods of Mesozoic activity with contrasting compositions. The older magmatic pulse is manifested by monzonites and monzodiorites from Tongshi for which zircon rims yield a concordant age of 177±4 Ma and the cores have a discordant age of ca. 2.5 Ga. Low MgO and Cr, high Na₂O contents and especially their isotopic compositions (⁸⁷Sr/⁸⁶Sr < 0.7042, ²⁰⁶Pb/²⁰⁴Pb < 16.8 and _Nd ~ –12) are consistent with derivation from late Archean–Paleoproterozoic lower crust. Relatively high HREE contents in these Jurassic plutons indicate a garnet-free source (<32 km), in contrast to the garnet-bearing source (>40 km) of the late Mesozoic high Sr and low Y granitoids from the same region. Distinctively different depths of crustal melting suggest dynamic thickening of the crust by magmatic underplating during the Jurassic and Cretaceous. The younger dioritic plutons from Laiwu and Yinan were emplaced at 132–126 Ma and show relatively high MgO and Cr contents and large isotopic variability. They were likely derived from enriched lithospheric mantle source and were subjected to crustal contamination during magma evolution. Early Cretaceous mantle melting is coeval with the widespread late Yanshanian granitic magmatism in North China. Early Cretaceous time may correspond to a critical period when a temperature increase due to lithospheric thinning allowed the intersection of the local geotherm and the wet peridotite solidus. While some mantle-derived magmas were erupted, most were trapped at variable crustal depths, triggering large-scale concomitant melting of the crust. Lithospheric thinning must have continued until the late Cretaceous because of the change in the source of mafic magmas from lithospheric to asthenospheric at that time. It is proposed that removal of the lithospheric keel beneath the North China craton may have been initiated as early as the Jurassic, but with the most intense period in the Cretaceous between 130–75 Ma. Such a relatively long timescale (~100 Ma) emphasizes the role of thermomechanical erosion by convective mantle in lithospheric thinning beneath this region.     

中国大陆克拉通地体地壳密度结构特征

本文根据重力场小波变换的尺度—源深度转换律,进行地面重力异常场多尺度分解,取得了反映中国地壳不同埋藏深度的小波细节,揭示了克拉通地体地壳的密度结构。中国大陆克拉通地体地壳密度结构的总体特点是:地壳密度总体偏高,尤其是在经历了岩浆底垫作用后,下地壳形成了粗大的和高密度的壳根,而造山带的下地壳通常有低密度的山根。现今的克拉通地体有几个演化的阶段,由陆地核演化为一定规模的克拉通地体可能经历漫长的时间。不同演化的阶段地体的物理性质有明显的变化,要赋予大地构造单元以演化阶段的内涵,才能准确理解它们的行为和属性。同时,根据地壳密度结构成像的结果,可以了解克拉通地体的成熟度。中国大陆克拉通地体或者内部地块的地壳密度结构大致可以分为三类。上中下地壳密度都偏高的属于第一类,称为典型的克拉通,塔里木、阿拉善、华北东部、上中扬子等属于第一类。第二类克拉通地体上地壳密度偏低,而下地壳密度偏高,称为弱结晶基底的克拉通,鄂尔多斯和佳木斯属于第二类。第三类上中地壳密度偏高,而下地壳密度偏低,称为无壳根的克拉通残片,柴达木和准噶尔盆地属于此类,它们可能是古克拉通地体的残片。     

Geophysical constraints on mesozoic disruption of North China Craton by underplating‐triggered lower‐crust flow of the Archaean lithosphere

The mechanism of the disruption, both lithospheric thinning and oceanization of the commonly accepted long‐term‐stable Archaean craton, is still an open question. The available models, all imply a bottom to top process. With the construction of a 1660‐km‐long transect across the eastern North China Craton (NCC), we demonstrate that both the P‐wave velocity and density in the lowermost crust beneath the central section are significantly higher than in the corresponding parts of the south and north sections on the transect. These features are interpreted as geophysical signature of lower crustal underplating, which supplies sufficiently high gravitational potential energy to trigger lateral flow of the lower crust. This magma underplating‐triggered bilateral lower crust flow may facilitate the lithospheric thinning by means of asthenosphere upwelling and decompression melting, which infill the gap produced by the lower crust flow. The underplating‐triggered lower crustal flow can provide an alternative mechanism to explain the NCC lithosphere disruption, which highlights the crustal feedback to Archaean lithosphere disruption, from top to bottom.     

The three-dimensional shear velocity structure of lithosphere in the southern Baikal rift system and its surroundings

The three-dimensional shear velocity lithospheric structure at depths from 0 to 70 km beneath the southern Baikal rift system and its surroundings has been imaged by inversion of P-to-SV receiver functions from 46 digital stations operated in two teleseismic international projects in southern Siberia and Mongolia. The receiver functions were determined from teleseismic P waveforms and inverted to obtain depth dependences of S velocities at each station which were related to tectonic structures. The computed vertical and horizontal sections of the 3D shear velocity model imaged a transition from relatively thin crust of the southern Siberian craton to thicker crust in the folded area south and southeast of Lake Baikal, with a local zone of thin crust right underneath the South Baikal basin. The velocity structure beneath the Baikal rift, the mountains of Transbaikalia, Mongolia, and the southern craton margin includes several low-velocity zones at different depths in the crust. Some of these zones may record seismic anisotropy associated with mylonite alignment along large thrusts.     

Moho depths and three-dimensional velocity structure of the crust and upper mantle beneath the Baikal region,from local tomography

We studied the 3D velocity structure of the crust and uppermost mantle beneath the Baikal region using tomographic inversion of ∼25,000 P and S arrivals from more than 1200 events recorded by 86 stations of three local seismological networks. Simultaneous iterative inversion with a new source location algorithm yielded 3D images of P and S velocity anomalies in the crust and upper mantle, a 2D model of Moho depths, and corrections to source coordinates and origin times. The resolving power of the algorithm, its stability against variations in the starting model, and the reliability of the final results were checked in several tests. The 3D velocity structure shows a well-pronounced low-velocity zone in the crust and uppermost mantle beneath the southwestern flank of the Baikal rift which matches the area of Cenozoic volcanism and a high velocity zone beneath the Siberian craton. The Moho depth pattern fits the surface tectonic elements with thinner crust along Lake Baikal and under the Busiyngol and Tunka basins and thicker crust beneath the East Sayan and Transbaikalian mountains and under the Primorsky ridge on the southern craton border.     

Fossil flat-slab subduction beneath the Illinois basin, USA

The Illinois basin is one of several well-studied intracratonic sedimentary basins within the North American craton whose formational mechanisms and subcrustal structure are not well understood. We study the S-velocity structure of the upper mantle beneath the Illinois basin and its surrounding area through seismic tomography. We utilize continental scale waveform data of seismic S and surface waves, enhanced by regional earthquakes located near the Illinois basin. Our 3D tomographic model, IL05, confirms the existence of a slow S-velocity structure in the uppermost mantle beneath the Illinois basin region. This anomalously slow region exists from the base of the crust to depths of  90 km, and is slower than the North American cratonic average by about 200 m/s. This anomalous uppermost mantle beneath the Illinois basin is underlain by a faster lithosphere, typical of the surrounding craton, to depths of  200 km. Excluding the formation of the Reelfoot Rift, this area of North American has been stable for over 1.0 Gy. Thus, we do not expect thermal anomalies from before that time to persist into present day S-velocity anomalies and we consider a delamination origin as an explanation of Illinois basin subsidence unlikely. We cannot rule out that the slow mid-lithosphere beneath the Illinois basin is caused by an uppermost mantle enriched by a deep, but weak plume. We attribute the slow mid-lithosphere to the presence of either oceanic, hydrous crust, or, a relatively cool mantle wedge with preserved hydrous minerals in the Illinois basin's uppermost mantle, related to a fossilized flat subduction zone.     

Petrological-geophysical models of the internal structure of the lithospheric mantle of the Siberian Craton

Based on the simultaneous inversion of unique ultralong-range seismic profiles Craton, Kimberlite, Meteorite, and Rift, sourced by peaceful nuclear and chemical explosions, and petrological and geochemical data on the composition of xenoliths of garnet peridotite and fertile primitive mantle material, the first reconstruction was obtained for the thermal state and density of the lithospheric mantle of the Siberian craton at depths of 100–300 km accounting for the effects of phase transformation, anharmonicity, and anelasticity. The upper mantle beneath Siberia is characterized by significant variations in seismic velocities, relief of seismic boundaries, degree of layering, and distribution of temperature and density. The mapping of the present-day lateral and vertical variations in the thermal state of the mantle showed that temperatures in the central part of the craton at depths of 100–200 km are somewhat lower than those at the periphery and 300–400°C lower than the mean temperature of tectonically younger mantle surrounding the craton. The temperature profiles derived from the seismic models lie between the 32.5 and 35 mW/m² conductive geotherms, and the mantle heat flow was estimated as 11–17 mW/m². The depth of the base of the cratonic thermal lithosphere (thermal boundary layer) is close to the 1450 ± 100°C isotherm at 300 ± 30 km, which is consistent with published heat flow, thermobarometry, and seismic tomography data. It was shown that the density distribution in the Siberian cratonic mantle cannot be described by a single homogeneous composition, either depleted or enriched. In addition to thermal anomalies, the mantle density heterogeneities must be related to variations in chemical composition with depth. This implies significant fertilization at depths greater than 180–200 km and is compatible with the existence of chemical stratification in the lithospheric mantle of the craton. In the asthenosphere-lithosphere transition zone, the craton root material is not very different in chemical composition, thermal regime, and density from the underlying asthenosphere. It was shown that minor variations in the chemical composition of the cratonic mantle and position of chemical (petrological) boundaries and the lithosphere-asthenosphere boundary cannot be reliably determined from the interpretation of seismic velocity models only.     

地震层析地质解释原理及其在滇西深部构造研究中的应用

本文首先简要地介绍了地震层析方法的原理,着重探讨了三维速度异常地质解释的几个基本问题,包括划分构造层的物理基础、确定分辨尺度的依据以及估计地质时限依据的岩石热力学过程。然后介绍了应用地震层析及地质解释方法研究滇西造山带的地壳、下岩石圈及软流圈构造所得到的主要结果,并与人工地震测深及大地电磁测深结果进行了比较。这一工作说明地震层析可以作为岩石圈与软流圈构造研究的一种有效手段,在造山带深部构造的研究中将发挥独特的作用。     

2D model of the crust and uppermost mantle along rift profile, Siberian craton

A two-dimensional model of the crust and uppermost mantle for the western Siberian craton and the adjoining areas of the Pur-Gedan basin to the north and Baikal Rift zone to the south is determined from travel time data from recordings of 30 chemical explosions and three nuclear explosions along the RIFT deep seismic sounding profile. This velocity model shows strong lateral variations in the crust and sub-Moho structure both within the craton and between the craton and the surrounding region. The Pur-Gedan basin has a 15-km thick, low-velocity sediment layer overlying a 25-km thick, high-velocity crystalline crustal layer. A paleo-rift zone with a graben-like structure in the basement and a high-velocity crustal intrusion or mantle upward exists beneath the southern part of the Pur-Gedan basin. The sedimentary layer is thin or non-existent and there is a velocity reversal in the upper crust beneath the Yenisey Zone. The Siberian craton has nearly uniform crustal thickness of 40–43 km but the average velocity in the lower crust in the north is higher (6.8–6.9 km/s) than in the south (6.6 km/s). The crust beneath the Baikal Rift zone is 35 km thick and has an average crustal velocity similar to that observed beneath the southern part of craton. The uppermost mantle velocity varies from 8.0 to 8.1 km/s beneath the young West Siberian platform and Baikal Rift zone to 8.1–8.5 km/s beneath the Siberian craton. Anomalous high Pn velocities (8.4–8.5 km/s) are observed beneath the western Tunguss basin in the northern part of the craton and beneath the southern part of the Siberian craton, but lower Pn velocities (8.1 km/s) are observed beneath the Low Angara basin in the central part of the craton. At about 100 km depth beneath the craton, there is a velocity inversion with a strong reflecting interface at its base. Some reflectors are also distinguished within the upper mantle at depth between 230 and 350 km.     

Origin of high-velocity anomalies beneath the Siberian craton: A fingerprint of multistage magma underplating since the Neoarchean

Despite the violent eruption of the Siberian Traps at ~ 250 Ma, the Siberian craton has an extremely low heat flow (18–25 mW/m²) and a very thick lithosphere (300–350 km), which makes it an ideal place to study the influence of mantle plumes on the long-term stability of cratons. Compared with seismic velocities of rocks, the lower crust of the Siberian craton is composed mainly of mafic granulites and could be rather heterogeneous in composition. The very high V_p (> 7.2 km/s) in the lowermost crust can be fit by a mixture of garnet granulites, two-pyroxene granulites, and garnet gabbro due to magma underplating. The high-velocity anomaly in the upper mantle (V_p = 8.3-8.6 km/s) can be interpreted by a mixture of eclogites and garnet peridotites. Combined with the study of lower crustal and mantle xenoliths, we recognized multistage magma underplating at the crust-mantle boundary beneath the Siberian craton, including the Neoarchean growth and Paleoproterozoic assembly of the Siberian craton beneath the Markha terrane, the Proterozoic collision along the Sayan-Taimyr suture zone, and the Triassic Siberian Trap event beneath the central Tunguska basin. The Moho becomes a metamorphism boundary of mafic rocks between granulite facies and eclogite facies rather than a chemical boundary that separates the mafic lower crust from the ultramafic upper mantle. Therefore, multistage magma underplating since the Neoarchean will result in a seismic Moho shallower than the petrologic Moho. Such magmatism-induced compositional change and dehydration will increase viscosity of the lithospheric mantle, and finally trigger lithospheric thickening after mantle plume activity. Hence, mantle plumes are not the key factor for craton destruction.     

华北东部中生代晚期—新生代软流圈上涌与岩石圈减薄

现今的地幔是由软流圈地幔(热的,主元素饱满、微量元素亏损的,塑性流变性质的)、古老岩石圈地幔(地幔1,冷的,主元素贫瘠、微量元素富集的,刚性的,以方辉橄榄岩为代表)以及现今的岩石圈地幔(地幔2,主元素饱满、微量元素亏损,以二辉橄榄岩为代表,可能包括多时期形成的)组成。古老岩石圈地幔与地幔2样品的共存、100~4·3Ma在地幔内部持续发生的古老岩石圈与软流圈的相互作用以及上述作用的时空不均一性,都表明了岩石圈减薄是软流圈呈“蘑菇云状”大规模上涌的结果。上述事件发生于100Ma以后。软流圈来源的玄武岩大范围喷发并伴随了岩石圈的强烈拉伸是事件发生的主要标志,岩石圈减薄是一个深部地质过程,不像是突发事件。     

Integrated 3D density modelling and segmentation of the Dead Sea Transform

A 3D interpretation of the newly compiled Bouguer anomaly in the area of the “Dead Sea Rift” is presented. A high-resolution 3D model constrained with the seismic results reveals the crustal thickness and density distribution beneath the Arava/Araba Valley (AV), the region between the Dead Sea and the Gulf of Aqaba/Elat. The Bouguer anomalies along the axial portion of the AV, as deduced from the modelling results, are mainly caused by deep-seated sedimentary basins (D > 10 km). An inferred zone of intrusion coincides with the maximum gravity anomaly on the eastern flank of the AV. The intrusion is displaced at different sectors along the NNW–SSE direction. The zone of maximum crustal thinning (depth 30 km) is attained in the western sector at the Mediterranean. The southeastern plateau, on the other hand, shows by far the largest crustal thickness of the region (38–42 km). Linked to the left lateral movement of approx. 105 km at the boundary between the African and Arabian plate, and constrained with recent seismic data, a small asymmetric topography of the Moho beneath the Dead Sea Transform (DST) was modelled. The thickness and density of the crust suggest that the AV is underlain by continental crust. The deep basins, the relatively large intrusion and the asymmetric topography of the Moho lead to the conclusion that a small-scale asthenospheric upwelling could be responsible for the thinning of the crust and subsequent creation of the Dead Sea basin during the left lateral movement. A clear segmentation along the strike of the DST was obtained by curvature analysis: the northern part in the neighbourhood of the Dead Sea is characterised by high curvature of the residual gravity field. Flexural rigidity calculations result in very low values of effective elastic lithospheric thickness (t _e < 5 km). This points to decoupling of crust in the Dead Sea area. In the central, AV the curvature is less pronounced and t _e increases to approximately 10 km. Curvature is high again in the southernmost part near the Aqaba region. Solutions of Euler deconvolution were visualised together with modelled density bodies and fit very well into the density model structures. An erratum to this article can be found at     

Geophysical profiling across the Sulu ultra-high-pressure metamorphic belt, eastern China

The largest ultra-high pressure metamorphic (UHPM) belt in the world is located along the Dabie–Sulu region, which tectonically belongs to the east part of the central orogenic belt of China. Integrated geophysical investigations of using deep seismic reflection, MT, and geothermal observations have been carried out in the Sulu area since 1997. The results of integrated interpretation suggest the existence of three features: (1) a rift beneath the Lianshui basin by the Jiashan–Xionshui fault; (2) a special crustal pattern, called the magmatic multi-arch structure occurs beneath the northern Sulu UHPM zone; and (3) a northwest-dipping regional thrust crosses the Sulu crust, representing the intracontinental subduction of the Yangtze craton beneath the Sulu metamorphic belts after collision between the Yangtze and Sino-Korean cratons. A magmatic multi-arch structure consists of some arched reflectors that occur in both the lower and the upper crust where arched reflectors coincide with granitoid plutons. The multi-arch structures are common in eastern China where many Mesozoic granitoid plutons of different scales occur. The crustal structures in the Sulu metamorphic belts resulted from intensive dynamic processes following the Triassic collision between the Yangtze and Sino-Korean cratons. The formation and exhumation of UHPM rocks followed the collision, and then intracontinental subduction of the Yangtze craton beneath the Dabie–Sulu terranes took place in the early and middle Jurassic. In the late Jurassic, the Sulu lithosphere turned to an extensional regime, large-scale granitic intrusions occurred in eastern China; these likely resulted from lithospheric thinning and asthenospheric uplifting. The granitic intrusions came to a climax during the Cretaceous and were followed by rifting along existing faults in the early Eogene, resulting in many petroleum basins. The granitoid emplacement that generated the magmatic multi-arch structure and the rift were consequences of the lithospheric thinning process, and deep intracontinental subduction of the Yangtze craton beneath the Sulu metamorphic belt might partially contribute to the lithospheric thinning.     

华北地区新生代岩石圈伸展减薄机制的数值模拟

新生代时期华北东部裂谷的伸展减薄机制及其周边的构造应力场,西部鄂尔多斯克拉通的抬升和周边断陷盆地的形成机制是目前研究的热点问题,但是较少有人从数值模拟的角度进行探讨。笔者采用有限元程序FEVPLIB对该地区5个剖面进行了模拟,初步取得如下认识:①在太平洋俯冲带的附近岩石圈伸展减薄较强,这与剖面经过的冲绳海槽正在拉开是吻合的,而太平洋的俯冲对较远的华北盆地的伸展减薄的影响较弱;②火山喷发时期,华北盆地有大的软流圈物质上涌造成华北裂谷的伸展减薄,符合纯剪切的机制,现今华北地区已趋于均衡,动力正趋于稳态;③六盘山逆冲在鄂尔多斯块体之上,代表着青藏高原东北缘的挤压,对华北是一个大的推挤力,可诱发鄂尔多斯块体的隆升,而鄂尔多斯向东北方向移动时提供了周边盆地的拉张的背景;④华北地区岩石圈的伸展减薄是六盘山处的挤压和东部太平洋板块俯冲两者联合的影响。模拟的结果与研究区GPS、重力异常以及岩石圈三维结构是吻合的。     

藏北岩石圈厚度与减薄机制分析

天然地震S波和大地电磁测深给出了两种不同的藏北岩石圈厚度模型,两种测量结果的地质含义至今还不十分清楚。通过对地表高程与地壳厚度回归关系的研究,以回归直线的斜率和截距作为地壳和岩石圈地幔平均密度取值的约束,并考虑相变因素对软流圈密度的影响,采用均衡理论对藏北岩石圈厚度进行了计算。计算结果表明,在可能的软流圈温度取值范围内藏北岩石圈的平均厚度约为106～120km,地壳增厚前的岩石圈平均厚度约80km。藏北新生代火山作用和岩浆起源-分凝深度分析表明,藏北现今岩石圈厚度主要受金云母脱水深度所控制。增厚前岩石圈地幔底部温度高于橄榄岩湿固相线温度,并受闪石和金云母高压脱水作用的影响。加厚岩石圈地幔因其底部不断发生脱水低程度熔融而进入软流圈小尺度对流体系,使岩石圈加厚过程中伴随有底部的脉动减薄作用。