华北克拉通及东邻西太平洋活动大陆边缘地区的P波速度结构：对岩石圈减薄动力学过程的探讨

本文通过地震层析成像研究获得了华北克拉通及其东邻地区（30°N-50°N,95°E -145°E）1°×1°的P波速度扰动图像.结果显示,在西太平洋俯冲带地区,上地幔中西倾的板片状高速异常体与其上方的低速异常区构成俯冲带与上覆地幔楔的典型速度结构式样.俯冲板片高速体在约300~400 km深度范围内被低速物质充填,暗示俯冲板片可能发生了断离.在华北克拉通地区的上地幔中发现三个东倾排列的高速异常带.在此基础上,本文构建了华北克拉通及其东邻西太平洋活动大陆边缘地区的上地幔速度结构模式图,并据此探讨克拉通岩石圈减薄与西太平洋活动大陆边缘的深部动力学联系.本文认为,太平洋板片的俯冲（断离）,触发热地幔物质上涌并在上覆地幔楔中形成对流,使克拉通岩石圈受到改造（底侵与弱化）.随着俯冲板片后撤,地幔楔中的对流场以及对岩石圈改造的影响范围均随之东移,最终导致华北克拉通岩石圈自下而上、从西向东分三个阶段依次拆沉减薄.这一模式能很好地解释现今克拉通岩石圈自西向东呈台阶状减薄的深部现象.     

Cretaceous episodic growth of the Japanese Islands

Abstract The Japanese Islands formed rapidly in situ along the eastern Asian continental margin in the Cretaceous due to both tectonic and magmatic processes. In the Early Cretaceous, huge oceanic plateaus created by the mid-Panthalassa super plume accreted with the continental margin. This tectonic interaction of oceanic plateau with continental crust is one of the significant tectonic processes responsible for continental growth in subduction zones. In the Japanese Islands, Late Cretaceous-Early Paleogene continental growth is much more episodic and drastic. At this time the continental margin uplifted regionally, and intra-continent collision tectonics took place in the northern part of the Asian continent. The uplifting event appears to have been caused by the subduction of very young oceanic crust (i.e. the Izanagi-Kula Plate) along the continental margin. Magmatism was also very active, and melting of the young oceanic slab appears to have resulted in ubiquitous plutons in the continental margin. Regional uplift of the continental margin and intra-continent collision tectonics promoted erosion of the uplifted area, and a large amount of terrigenous sediment was abruptly supplied to the trench. As a result of the rapid supply of terrigenous detritus, the accretionary complexes (the Hidaka Belt in Hokkaido and the Shimanto Belt in Southwest Japan) grew rapidly in the subduction zone. The rapid growth of the accretionary complexes and the subduction of very young, buoyant oceanic crust caused the extrusion of a high-P/T metamorphic wedge from the deep levels of the subduction zone. Episodic growth of the Late Cretaceous Japanese Islands suggests that subduction of very young oceanic crust and/or ridge subduction are very significant for the formation of new continental crust in subduction zones.     

Subduction of continental margins and the uplift of high-pressure metamorphic rocks

The mechanism by which high-pressure metamorphosed continental material is emplaced at high structural levels is a major unsolved problem of collisional orogenesis. We suggest that the emplacement results from partial subduction of the continental margin which, because of its high flexural rigidity, produces a rapid change in the trajectory of the descending slab. We assume a two-fold increase in effective elastic thickness of the lithosphere as the continental margin approaches the subduction zone, and calculate the flexural profile of a thin plate for progressive downward migration of the zone of increased rigidity. We assess the effect of changes in the flexural profile on the overlying accretionary prism and mantle wedge as the continent approaches by estimating the extra stresses that are imposed on the wedge due to the bending moment exerted by the continental part of the plate. The wedges overlying the subduction zones, and the subducting slab itself, experience substantial extra compressional stress at depths of around 100 km, and extensional stress at shallower depths, as the continental margin passes through the zone of maximum curvature. The magnitudes of such extra stresses are probably adequate to effect significant deformation of the wedge and/or the descending plate, and are experienced in a time interval of less than 5 m.y. for typical subduction rates. The spatial variation of yield stresses in the region of the wedge and descending slab indicates that much of this deformation may be taken up in the crustal part of the descending slab, which is the weakest region in the deeper parts of the subduction zone. This may result in rapid upward migration of the crust of the partially subducted continental margin, against the flow of subduction. High-pressure metamorphosed terranes emplaced by the mechanism envisaged in this paper would be bounded by thrust faults below and normal faults above. Movement on the faults would have been coeval, and would have resulted in rapid unroofing of the high-pressure terranes, synchronous with arrival of the continental margin at the subduction zone and, therefore, relatively early in the history of a collisional orogen.     

The seismogenic zone of subduction thrust faults

Abstract Subduction thrust faults generate earthquakes over a limited depth range. They are aseismic in their seaward updip portions and landward downdip of a critical point. The seaward shallow aseismic zone, commonly beneath accreted sediments, may be a consequence of unconsolidated sediments, especially stable-sliding smectite clays. Such clays are dehydrated and the fault may become seismogenic where the temperature reaches 100--150°C, that is, at a 5--15 km depth. Two factors may determine the downdip seismogenic limit. For subduction of young hot oceanic lithosphere beneath large accretionary sedimentary prisms and beneath continental crust, the transition to aseismic stable sliding is temperature controlled. The maximum temperature for seismic behavior in crustal rocks is ～ 350°C, regardless of the presence of water. In addition, great earthquake ruptures initiated at less than this temperature may propagate with decreasing slip to where the temperature is ～ 450°C. For subduction beneath thin island arc crust and beneath continental crust in some areas, the forearc mantle is reached by the thrust shallower than the 350°C temperature. The forearc upper mantle probably is aseismic because of stable-sliding serpentinite hydrated by water from the underthrusting oceanic crust and sediments. For many subduction zones the downdip seismogenic width defined by these limits is much less than previously assumed. Within the narrowly defined seismic zone, most of the convergence may occur in earthquakes. Numerical thermal models have been employed to estimate temperatures on the subduction thrust planes of four continental subduction zones. For Cascadia and Southwest Japan where very young and hot plates are subducting, the downdip seismogenic limit on the subduction thrust is thermally controlled and is shallow. For Alaska and most of Chile, the forearc mantle is reached before the critical temperature, and mantle serpentinite provides the limit. In all four regions, the seismogenic zones so defined agree with estimates of the extent of great earthquake rupture, and with the downdip extent of the interseismic locked zone.     

Comparison of gravimetric and seismic constraints on the structure of the Aegean lithosphere in the forearc of the Hellenic subduction zone in the area of Crete

The island of Crete is located in the forearc of the Hellenic subduction zone, where the African lithospheric plate is subducting beneath the Eurasian one. The depth of the plate contact as well as the internal structure of the Aegean plate in the area of Crete have been a matter of debate. In this study, seismic constrains obtained by wide-angle seismic, receiver function and surface wave studies are discussed and compared to a 3D density model of the region.The interface between the Aegean continental lithosphere and the African one is located at a depth of about 50 km below Crete. According to seismic studies, the Aegean lithosphere in the area of Crete is characterised by strong lateral, arc–parallel heterogeneity. An about 30 km thick Aegean crust is found in central Crete with a density of about 2850 kg/m³ for the lower Aegean continental crust and a density of about 3300 kg/m³ for the mantle wedge between the Aegean crust and the African lithosphere. For the deeper crust in the area of western Crete two alternative models have been proposed by seismic studies. One with an about 35 km thick crust and another one with crustal velocities down to the plate contact. A grid search is performed to test the consistency of these models with gravimetric constraints. For western Crete a model with a thick lower Aegean crust and a density of about 2950 kg/m³ is favoured. The inferred density of the lower Aegean crust in the area of Crete correlates well with S-wave velocities obtained by surface wave studies.Based on the 3D density model, the weight of the Aegean lithosphere is estimated along an E–W oriented profile in the area of Crete. Low weights are found for the region of western Crete.     

Time dependence of negative buoyancy and the subduction of continental lithosphere

The time evolution of negative buoyancy of a subducting slab is modelled from the beginning of subduction under various kinematic conditions (dip angle and subduction velocity). The calculations take into account the thermal and density effects of the variations of the thermophysical parameters with temperature and pressure, and of phase transitions. The magnitude of the negative buoyancy increases during subduction of oceanic lithosphere, up to values in the (2–4) × 10¹³ N m⁻¹ range when the tip of the slab reaches a depth of 600–700 km. If continental material arrives at the trench and is subducted, the downward buoyancy decreases by an amount proportional to the volume of the subducted continental crust. Assuming that subduction stops when the buoyancy becomes zero, and that delamination of the continental crust or slab breakoff do not occur, the maximum downdip length of the subductable continental crust is estimated as a function of the dip angle, subduction velocity and geometry of the margin. In most cases, subduction of continental material down to depths of 100–250 km is possible, and continental subduction can continue for times up to 10–15 Ma if the velocity is low. These estimates are not significantly affected by the hypothetical occurrence of a metastable olivine wedge within the slab, and could be lower bounds if the lower continental crust is mafic and transforms to eclogite.     

The great Sumatra–Andaman earthquakes — Imaging the boundary between the ruptures of the great 2004 and 2005 earthquakes

Segmentation along convergent margins controls earthquake magnitude and location, but the physical causes of segment boundaries, and their impact on earthquake rupture dynamics, are still poorly understood. One aspect of the 2004 and 2005 great Sumatra–Andaman earthquakes is their abrupt termination along a common boundary. This has led to speculation on the nature of the boundary, its origin and why it was not breached.

For the first time the boundary has been imaged and, with newly acquired marine geophysical data, we demonstrate that a ridge on the subducting Indo-Australian oceanic crust may exert a control on margin segmentation. This suggests a lower plate influence on margin structure, particularly its segmentation. The ridge is masked by the sedimentary cover in the trench. Its most likely trend is NNE–SSW. It is interpreted as a fracture zone on the subducting oceanic plate. A ramp or tear along the eastern flank of the subducting fracture zone beneath Simeulue Island may be considered as an intensification factor in terms of rupture propagation barrier.     

中国大陆东南缘地震接收函数与地壳和上地幔结构

从2008-2011年,分别在中国大陆东南缘沿海和内陆两条NE向剖面上进行了宽频地震观测,利用记录到的远震波形资料提取得到1446个远震P波接收函数,用H-κ叠加扫描和CCP偏移叠加方法研究了中国大陆东南缘地壳及上地幔过渡带的结构及其变化特征.结合固定台网25个台站的H-κ结果,获得中国大陆东南缘(福建地区)地壳厚度从内陆到沿海逐渐减薄的图像:地壳从闽西北山区的33 km减薄到厦门沿海一带的29 km以下,平均地壳厚度为31.3 km,具有陆地向洋壳过渡的特征;地壳泊松比从内陆到沿海显示出分带特征,闽中西部内陆地区小于0.26,沿海地带高于0.26,且在断裂带的交汇区域表现为相对异常高值.地壳上地幔顶部(0~200 km)的CCP偏移叠加成像结果显示闽江断裂等NW向断裂深切Moho界面,在断裂两侧Moho面急剧抬升或下沉,产状改变,这些特征向内陆地区逐渐变得不明显.闽江等NW向断裂对研究区地壳厚度、地震等有明显控制作用.上地幔尺度(300~700 km)的CCP偏移叠加成像,未见410 km和660 km速度间断面突变和起伏异常,其绝对深度略大于IASP91模型的,上地幔转换带厚度正常(250±5 km),表明中国大陆东南缘上地幔转换带未受欧亚与菲律宾板块碰撞的明显影响,推断中国大陆东南缘及台湾海峡下方不存在俯冲板块,或俯冲前缘未扰动到410 km的深度.     

马尼拉海沟俯冲带热结构的模拟研究

俯冲带热结构的数值模拟研究是对地表观测研究的重要补充,也是验证地球动力学模型的重要方法.本文沿马尼拉海沟俯冲带东火山链(EVC)和西火山链(WVC)各取一条剖面,依据地质、地球物理条件,进行了有限元热模拟计算.计算过程中,分析了摩擦和剪切热对俯冲带热结构的影响,模拟了EVC和WVC两条测线下俯冲带的热结构,并结合岩石学实验结果预测了俯冲板块发生脱水和部分熔融的位置.模拟结果表明,在100 km深度处,考虑摩擦和剪切热时,俯冲板块表面的温度约为865 ℃;而不考虑摩擦和剪切时,俯冲板块表面的温度仅为770 ℃,二者温差可达95 ℃.在相同深度处,考虑摩擦和剪切热时,在EVC和WVC测线下俯冲板块表面的温度分别为865 ℃和895 ℃,俯冲洋壳底部温度分别为560 ℃和605 ℃.俯冲板块表面少量矿物开始脱水的深度小于50 km,但大量脱水和部分熔融主要发生在深度100 km左右,这与地表观测的火山活动位置一致.     

海山的倾斜俯冲对上覆板块变形的影响

伴随洋壳的俯冲,驼伏其上的海山会导致上覆板块的强烈变形.为解释该构造变形特征,本文运用物理模拟实验的方法,着重分析海山的斜向俯冲对上覆板块变形的影响,并将模拟结果与正向俯冲过程进行对比.实验结果显示:海山开始进入俯冲,前缘楔体的增生会被阻止,同时楔体被抬升并出现脱顶构造,未被海山破坏的楔体会出现后冲断层的激活,后冲断层轴平行于海山的俯冲方向.海山进一步俯冲,突起项部发育一系列张扭性质的微断裂和走滑性质的共轭断裂,尾随突起之后的楔体由于重力会产生正断层系统.比起正向俯冲,斜向俯冲过程中所产生的后逆冲体、海山两侧的叠瓦状逆冲推覆构造都出现不对称分布,断裂和微断裂束的走向不规则散开,后冲断层的轴向及海山俯冲过后在楔体上产生的凹槽的轨迹都不断斜向迁移,且凹槽两侧的地势不一致等.最后利用文中的物理模拟结果,很好的解释了马尼拉海沟中段俯冲构造的构造特征,同时对其他俯冲大陆边缘的构造解释具有指导意义.     

Origin of the early Cenozoic belt boundary thrust and Izanagi–Pacific ridge subduction in the western Pacific margin

The belt boundary thrust within the Cretaceous–Neogene accretionary complex of the Shimanto Belt, southwestern Japan, extends for more than ~ 1 000 km along the Japanese islands. A common understanding of the origin of the thrust is that it is an out of sequence thrust as a result of continuous accretion since the late Cretaceous and there is a kinematic reason for its maintaining a critically tapered wedge. The timing of the accretion gap and thrusting, however, coincides with the collision of the Paleocene–early Eocene Izanagi–Pacific spreading ridges with the trench along the western Pacific margin, which has been recently re‐hypothesized as younger than the previous assumption with respect to the Kula‐Pacific ridge subduction during the late Cretaceous. The ridge subduction hypothesis provides a consistent explanation for the cessation of magmatic activity along the continental margin and the presence of an unconformity in the forearc basin. This is not only the case in southwestern Japan, but also along the more northern Asian margin in Hokkaido, Sakhalin, and Sikhote‐Alin. This Paleocene–early Eocene ridge subduction hypothesis is also consistent with recently acquired tomographic images beneath the Asian continent. The timing of the Izanagi–Pacific ridge subduction along the western Pacific margin allows for a revision of the classic hypothesis of a great reorganization of the Pacific Plate motion between ~ 47 Ma and 42 Ma, illustrated by the bend in the Hawaii–Emperor chain, because of the change in subduction torque balance and the Oligocene–Miocene back arc spreading after the ridge subduction in the western Pacific margin.     

Implications from the seismic crustal structure of the northern Izu–Bonin arc

The results of a controlled source seismic reflection–refraction experiment carried out in 1992 reveal the following characteristics of the northern Izu–Bonin (Ogasawara) oceanic island arc–trench system. (1) The crust rapidly thickens from the Shikoku back-arc basin to the arc, is thickest beneath the active rifts, and then gradually thins to the forearc. The thickness of the crust beneath the arc rift zone and the back-arc basin are ∼ 20 km and 8 km, respectively. (2) The Moho vanishes beneath the forearc. Velocities rapidly decrease eastwards beneath the inner trench wall. (3) The velocity of the lower crust of the arc and the back-arc basin is 7.1–7.3 km/s. This velocity is higher than the typical oceanic lower crust whose velocity is ∼ 6.7 km/s. (4) The velocity of the middle crust of the arc is ∼ 6 km/s. This layer does not exist beneath the back-arc basin. (5) A slight difference in the velocity gradient of the middle crust exists between the arc rift zone and the forearc. Based on these findings and previous studies, it is inferred that: (i) the middle crust is probably granitic rock and formed in more than two episodes; (ii) the lower crust formed by igneous underplating which may also have affected part of the back-arc basin; and (iii) the root of the serpentinite diapir on the inner trench wall is a low-velocity mantle wedge that was probably caused by large amounts of water released from the subducting Pacific plate at depths shallower than 30 km.     

Terrane analysis and accretion in North-East Asia

Abstract A terrane map of North-East Asia at 1:5 000 000 scale has been compiled. The map shows terranes of different types and ages accreted to the North-Asian craton in the Mesozoic–Cenozoic, sub-and superterranes, together with post-amalgamation and post-accretion assemblages. The great Kolyma-Omolon superterrane adjoins the north-east craton margin. It is composed of large angular terranes of continental affinity: craton fragments and fragments of the passive continental margin of Siberia, and island arc, oceanic and turbidite terranes that are unconformably overlain by shallow marine Middle-Upper Jurassic deposits. The superterrane resulted from a long subduction of the Paleo-Pacific oceanic crust beneath the Alazeya arc. Its south-west boundary is defined by the Late Jurassic Uyandina-Yasachnaya marginal volcanic arc which was brought about by subduction of the oceanic crust that separated the superterrane from Siberia. According to paleomagnetic evidence the width of the basin is estimated to be 1500–2000 km. Accretion of the superterrane to Siberia is dated to the late Late Jurassic-Neocomian. The north-east superterrane boundary is defined by the Lyakhov-South Anyui suture which extends across southern Chukotka up to Alaska. Collision of the superterrane with the Chukotka shelf terrane is dated to the middle of the Cretaceous. The Okhotsk-Chukotka belt, composed of Albian-Late Cretaceous undeformed continental volcan-ites, defines the Cretaceous margin of North Asia. Terranes eastward of the belt are mainly of oceanic affinity: island arc upon oceanic crust, accretion wedge and turbidite terranes, as well as cratonic terranes and fragments of magmatic arcs on the continental crust and metamorphic terranes of unclear origin and age. The time of their accretion is constrained by post-accretionary volcanic belts that extend parallel to the Okhotsk-Chukotka belt but are displaced to the east: the Maastrichtian-Miocene Kamchatka-Koryak belt and the Eocene-Quaternary Central Kamchatka belt which mark active margins of the continent of corresponding ages.     

马尼拉俯冲带北段增生楔前缘构造变形和精细结构

马尼拉俯冲带是南海的东部边界,记录了南海形成演化的关键信息,同时也是地震和海啸多发区域.本文利用过马尼拉俯冲带北段的高分辨率多道地震剖面,分析了研究区内海盆和海沟的沉积特征,精细刻画了区内增生楔前缘的构造变形、结构以及岩浆活动特征.研究区内增生楔下陆坡部分由盲冲断层、构造楔和叠瓦逆冲断层构成,逆冲断层归并于一条位于下中新统的滑脱面上,滑脱面向海方向的展布明显受到增生楔之下埋藏海山和基底隆起的影响;上陆坡的反射特征则因变形强烈和岩浆作用而难以识别;岩浆活动开始于晚中新世末期并持续至第四纪.马尼拉俯冲带北段增生楔的形成时间早于16.5 Ma,并通过前展式逆冲向南海方向扩展;马尼拉俯冲带的初始形成时间可能在晚渐新世,而此时南海海盆扩张仍在持续.南海东北缘19°N-21°N区域为南海北部陆坡向海盆的延伸,高度减薄的陆壳的俯冲造成马尼拉海沟北段几何形态明显地向东凹进.     

Pacific-type orogeny revisited: Miyashiro-type orogeny proposed

Abstract The concept of Pacific-type orogeny is revised, based on an assessment of geologic data collected from the Japanese Islands during the past 25 years. The formation of a passive continental margin after the birth of the Pacific Ocean at 600 Ma was followed by the initiation of oceanic plate subduction at 450 Ma. Since then, four episodes of Pacific-type orogeny have occurred to create an orogenic belt 400 km wide that gradually grew both oceanward and downward. The orogenic belt consists mainly of an accretionary complex tectonically interlayered with thin (<2 km thick), subhorizontal, high-P/T regional metamorphic belts. Both the accretionary complex and the high-P/T rocks were intruded by granitoids ～100 million years after the formation of the accretionary complex. The intrusion of calc-alkaline (CA) plutons was synchronous with the exhumation of high-P/T schist belts. Ages from microfossils and K-Ar analysis suggest that the orogenic climax happened at a time of mid-oceanic ridge subduction. The orogenic climax was characterized by the formation of major subhorizontal orogenic structures, the exhumation of high-P/T schist belts by wedge extrusion and subsequent domed uplift, and the intrusion-extrusion of CA magma dominantly produced by slab melting. The orogenic climax ended soon after ridge subduction, and thereafter a new Pacific-type orogeny began. A single Pacific-type orogenic cycle may correspond to the interaction of the Asian continental margin with one major Pacific oceanic plate. Ophiolites in Japan occur as accreted material and are not of island-arc but of plume origin. They presumably formed after the birth of the southern Pacific superplume at 600 Ma, and did not modify the cordilleran-type orogeny in a major way. Microplates, fore-arc slivers, intra-oceanic arc collisions and the opening of back-arc basins clearly contributed to cordilleran orogenesis. However, they were of secondary importance and served only to modify pre-existing major orogenic components. The most important cause of cordilleran-type orogeny is the subduction of a mid-oceanic ridge, by which the volume of continental crust increases through the transfer of granitic melt from the subducting oceanic crust to an orogenic welt. Accretionary complexes are composed mainly of recycled granitic sediments with minor amounts of oceanic material, which indicate that the accretion of oceanic material, including huge oceanic plateaus, was not significant for orogenic growth. Instead, the formation and intrusion of granitoids are the keys to continental growth, which is the most important process in Pacific-type orogeny. Collision-type orogeny does not increase the volume of continental crust. The name ‘Miyashiro-type orogeny’ is proposed for this revised concept of Pacific-type or cordilleran-type orogeny, in order to commemorate Professor A. Miyashiro's many contributions to a better understanding of orogenesis.     

三维板块几何形态对大陆深俯冲动力学的制约

大陆深俯冲及超高压变质作用是大陆动力学的重要研究内容,前人进行了系统的地质、地球物理观测以及数值模拟研究.然而,自然界中大陆板块的俯冲、碰撞及造山过程大部分具有明显的沿走向的差异性,这种典型的三维特征可能很大程度上依赖于会聚大陆板块的初始几何学和运动学特征.本文采用三维高分辨率的动力学数值模拟方法,建立了方形大陆板块和楔形大陆板块两种不同的俯冲-碰撞模型,并且俯冲大陆板块侧面与大洋俯冲带相邻.数值模拟结果揭示大洋板块可以持续地俯冲到地幔之中,而大陆板块俯冲到一定深度处,其前端的俯冲板块将发生断离,并进而造成残余的大陆板块俯冲角度的减小.方形大陆俯冲板块的断离深度约为150km,而楔形大陆俯冲板块的断离深度较大,约250~300km,这很大程度上取决于俯冲带中大洋板块的牵引力和大陆板块的负浮力之间的竞争关系.同时,无论方形还是楔形大陆板块俯冲模型中,板块断离后,侧向的大洋俯冲板块仍可以拖曳约60~70km宽的大陆边缘岩石圈持续向下俯冲,揭示了新西兰东部的洋-陆空间转换俯冲带的动力学机制.并且,数值模型与喜马拉雅造山带和秦岭—大别—苏鲁造山带进行了对比,进而对其高压-超高压岩石空间展布沿走向的差异性特征和机制提供了一定的启示.     

On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective

The origin of El Chichón volcano is poorly understood, and we attempt in this study to demonstrate that the Tehuantepec Ridge (TR), a major tectonic discontinuity on the Cocos plate, plays a key role in determining the location of the volcano by enhancing the slab dehydration budget beneath it. Using marine magnetic anomalies we show that the upper mantle beneath TR undergoes strong serpentinization, carrying significant amounts of water into subduction. Another key aspect of the magnetic anomaly over southern Mexico is a long-wavelength (∼ 150 km) high amplitude (∼ 500 nT) magnetic anomaly located between the trench and the coast. Using a 2D joint magnetic-gravity forward model, constrained by the subduction P–T structure, slab geometry and seismicity, we find a highly magnetic and low-density source located at 40–80 km depth that we interpret as a partially serpentinized mantle wedge formed by fluids expelled from the subducting Cocos plate. Using phase diagrams for sediments, basalt and peridotite, and the thermal structure of the subduction zone beneath El Chichón we find that ∼ 40% of sediments and basalt dehydrate at depths corresponding with the location of the serpentinized mantle wedge, whereas the serpentinized root beneath TR strongly dehydrates (∼90%) at depths of 180-200 km comparable with the slab depths beneath El Chichón (200-220 km). We conclude that this strong deserpentinization pulse of mantle lithosphere beneath TR at great depths is responsible for the unusual location, singularity and, probably, the geochemically distinct signature (adakitic-like) of El Chichón volcano.     

Mesozoic mafic magmatism in North China: Implications for thinning and destruction of cratonic lithosphere

The North China Craton (NCC) has been thinned from >200 km to <100 km in its eastern part. The ancient subcontinental lithospheric mantle (SCLM) has been replaced by the juvenile SCLM in the Meoszoic. During this period, the NCC was destructed as indicated by extensive magmatism in the Early Cretaceous. While there is a consensus on the thinning and destruction of cratonic lithosphere in North China, it has been hotly debated about the mechanism of cartonic destruction. This study attempts to provide a resolution to current debates in the view of Mesozoic mafic magmatism in North China. We made a compilation of geochemical data available for Mesozoic mafic igneous rocks in the NCC. The results indicate that these mafic igneous rocks can be categorized into two series, manifesting a dramatic change in the nature of mantle sources at ~121 Ma. Mafic igneous rocks emplaced at this age start to show both oceanic island basalts (OIB)-like trace element distribution patterns and depleted to weakly enriched Sr-Nd isotope compositions. In contrast, mafic igneous rocks emplaced before and after this age exhibit both island arc basalts (IAB)-like trace element distribution patterns and enriched Sr-Nd isotope compositions. This difference indicates a geochemical mutation in the SCLM of North China at ~121 Ma. Although mafic magmatism also took place in the Late Triassic, it was related to exhumation of the deeply subducted South China continental crust because the subduction of Paleo-Pacific slab was not operated at that time. Paleo-Pacific slab started to subduct beneath the eastern margin of Eruasian continent since the Jurrasic. The subducting slab and its overlying SCLM wedge were coupled in the Jurassic, and slab dehydration resulted in hydration and weakening of the cratonic mantle. The mantle sources of ancient IAB-like mafic igneous rocks are a kind of ultramafic metasomatites that were generated by reaction of the cratonic mantle wedge peridotite not only with aqueous solutions derived from dehydration of the subducting Paleo-Pacific oceanic crust in the Jurassic but also with hydrous melts derived from partial melting of the subducting South China continental crust in the Triassic. On the other hand, the mantle sources of juvenile OIB-like mafic igneous rocks are also a kind of ultramafic metasomatites that were generated by reaction of the asthenospheric mantle underneath the North China lithosphere with hydrous felsic melts derived from partial melting of the subducting Paleo-Pacific oceanic crust. The subducting Paleo-Pacific slab became rollback at ~144 Ma. Afterwards the SCLM base was heated by laterally filled asthenospheric mantle, leading to thinning of the hydrated and weakened cratonic mantle. There was extensive bimodal magmatism at 130 to 120 Ma, marking intensive destruction of the cratonic lithosphere. Not only the ultramafic metasomatites in the lower part of the cratonic mantle wedge underwent partial melting to produce mafic igneous rocks showing negative ε_Nd(t) values, depletion in Nb and Ta but enrichment in Pb, but also the lower continent crust overlying the cratonic mantle wedge was heated for extensive felsic magmatism. At the same time, the rollback slab surface was heated by the laterally filled asthenospheric mantle, resulting in partial melting of the previously dehydrated rocks beyond rutile stability on the slab surface. This produce still hydrous felsic melts, which metasomatized the overlying asthenospheric mantle peridotite to generate the ultramafic metasomatites that show positive ε_Nd(t) values, no depletion or even enrichment in Nb and Ta but depletion in Pb. Partial melting of such metasomatites started at ~121 Ma, giving rise to the mafic igneous rocks with juvenile OIB-like geochemical signatures. In this context, the age of ~121 Ma may terminate replacement of the ancient SCLM by the juvenile SCLM in North China. Paleo-Pacific slab was not subducted to the mantle transition zone in the Mesozoic as revealed by modern seismic tomography, and it was subducted at a low angle since the Jurassic, like the subduction of Nazca Plate beneath American continent. This flat subduction would not only chemically metasomatize the cratonic mantle but also physically erode the cratonic mantle. Therefore, the interaction between Paleo-Pacific slab and the cratonic mantle is the first-order geodynamic mechanism for the thinning and destruction of cratonic lithosphere in North China.     

洋脊和年轻洋壳俯冲与埃达克岩生成的热模拟

新生代以来 ,环太平洋周边分布的埃达克岩 (Adakite)主要与年轻洋壳俯冲时在 70～ 90km深处的部分熔融有关。利用数值方法 ,模拟了洋壳俯冲的热演化过程并讨论了脱水、熔融对埃达克岩浆活动的影响。结果表明 :仅在活动海岭俯冲前后约 10Ma内 ,年轻的、热的俯冲海洋板片在 75～85km深度范围内 ,温度升高至 82 5～ 10 0 0℃脱水 ,导致年轻洋壳中角闪岩部分熔融 ,形成埃达克岩(Adakite)。而一般洋壳俯冲在 10 0km以下深度才脱水 ,由于脱水区压力较高洋壳自身不能熔融 ,水进入上覆地幔楔状体导致部分熔融 ,形成安山岩 (Andesite     

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

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