Heat loss from the earth: A constraint on Archaean tectonics from the relation between geothermal gradients and the rate of plate production

The models suggested for the oceanic lithosphere which best predict oceanic heat flow and depth profiles are the constant thickness model and a model in which the lithosphere thickens away from the ridge with a heat source at its base. The latter is considered to be more physically realistic. Such a model, constrained by the observed oceanic heat flow and depth profiles and a temperature at the ridge crest of between 1100°C and 1300°C, requires a heat source at the base of the lithosphere of between 0.5 and 0.9 h.f.u., thermal conductivities for the mantle between 0.005 and 0.0095 cal cm⁻¹ °C⁻¹ s⁻¹ and a coefficient of thermal expansion at 840°C between 4.1 × 10⁻⁵ and 5.1 × 10⁻⁵ °C⁻¹. Plate creation and subduction are calculated to dissipate about 45% of the total earth heat loss for this model. The efficiency of this mechanism of heat loss is shown to be strongly dependent on the magnitude of the basal heat source. A relation is derived for total earth heat loss as a function of the rate of plate creation and the amount of heat transported to the base of plates. The estimated heat transport to the base of the oceanic lithosphere is similar to estimates of mantle heat flow into the base of the continental lithosphere. If this relation existed in the past and if metamorphic conditions in late Archaean high-grade terrains can be used to provide a maximum constraint on equilibrium Archaean continental thermal gradients, heat flow into the base of the lithosphere in the late Archaean must have been less than about 1.2–1.5 h.f.u. The relation between earth heat loss, the rate of plate creation and the rate of heat transport to the base of the lithosphere suggests that a significant proportion of the heat loss in the Archaean must have taken place by the processes of plate creation and subduction. The Archaean plate processes may have involved much more rapid production of plates only slightly thinner than at present.     

高温高压下斜长角闪岩电导率研究及其地球物理启示

中下地壳和俯冲带区域的高电导率异常（0.01~1 S·m^-1）可能与地球内部的特定物质及其变化有关.斜长角闪岩是中下地壳以及俯冲带区域的重要组成之一,高温高压下斜长角闪岩的电导率研究对认识电导率异常具有重要意义.本研究采用交流阻抗谱法,在0.5,1.0,1.5 GPa和473~1073 K条件下测量了天然斜长角闪岩样品的复阻抗,实验结果表明压力对斜长角闪岩的电导率影响非常小,而温度对于电导率影响非常显著,其电导率在1073 K可以达到10^-0.5 S·m^-1;实验获得的活化能值为52.21 kJ·mol^-1,推断其导电机制可能为小极化子传导（Fe²⁺的氧化）主导.结合本实验获得的结果与大陆岩石圈和俯冲带的温度结构,我们计算得到相应的电性结构剖面,并与三种不同构造背景下的大陆岩石圈（克拉通、大陆裂谷和活动造山带）和俯冲带区域的电磁剖面结构进行了对比研究,结果发现斜长角闪岩可以解释大陆裂谷和活动造山带构造背景下的莫霍面附近的高电导率异常现象,同时可能是导致较热的俯冲带区域（例如卡斯卡迪地区）高电导率异常现象的原因.     

花岗岩非稳态流变实验

脆塑性转化带对于研究岩石圈变形、断层强度和变形机制以及强震的孕育和发生具有重要意义.文中采用汶川地震震源区彭灌杂岩中具有代表性的细粒花岗岩样品,在固体压力介质三轴实验系统上开展了高温高压非稳态流变实验研究.实验设计模拟了汶川地震区地壳10～30km深度的实际温度和压力,温度为190～490℃,压力为250～750MPa...     

Ar/Ar geochronology of the West Greenland Tertiary volcanic province

Paleocene volcanic rocks in West Greenland and Baffin Island were among the first products of the Iceland mantle plume, forming part of a larger igneous province that is now submerged beneath the northern Labrador Sea. A ⁴⁰Ar/³⁹Ar dating study shows that volcanism commenced in West Greenland between 60.9 and 61.3 Ma and that 80% of the Paleocene lava pile was erupted in 1 million years or less (weighted mean age of 60.5±0.4 Ma). Minimum estimates of magma production rates (1.3×10⁻⁴ km³ year⁻¹ km⁻¹) are similar to the present Iceland rift, except for the uppermost part of the Paleocene volcanic succession where the rate decreases to <0.7×10⁻⁴ km³ year⁻¹ km⁻¹ (rift). The timing of onset of volcanism in West Greenland coincides with the opening of the northern Labrador Sea and is also strikingly similar to the age of the oldest Tertiary volcanic rocks from offshore SE Greenland and the British–Irish province. This is interpreted as manifesting the impact and rapid (>1 m/year) lateral spreading of the Iceland plume head at the base of the Greenland lithosphere at 62 Ma. We suggest that the arrival, or at least a major increase in the flux, of the Iceland mantle plume beneath Greenland was a contributing factor in the initiation of seafloor spreading in the northern Labrador Sea. Our study has also revealed a previously unrecognised Early Eocene volcanic episode in West Greenland. This magmatism may be related to movement on the transform Ungava Fault System which transferred drifting from the Labrador Sea to Baffin Bay. A regional change in plate kinematics at 55 Ma, associated with the opening of the North Atlantic, would have caused net extension along parts of this fault. This would have resulted in decompression and partial melting of the underlying asthenosphere. The source of the melts for the Eocene magmatism may have been remnants of still anomalously hot Iceland plume mantle which were left stranded beneath the West Greenland lithosphere in the Early Paleocene.     

Link between SSZ ophiolite formation, emplacement and arc inception, Northland, New Zealand: U–Pb SHRIMP constraints; Cenozoic SW Pacific tectonic implications

New U–Pb age-data from zircons separated from a Northland ophiolite gabbro yield a mean ²⁰⁶Pb/²³⁸U age of 31.6 ± 0.2 Ma, providing support for a recently determined 28.3 ± 0.2 Ma SHRIMP age of an associated plagiogranite and  29–26 Ma ⁴⁰Ar/³⁹Ar ages (n = 9) of basalts of the ophiolite. Elsewhere, Miocene arc-related calc-alkaline andesite dikes which intrude the ophiolitic rocks contain zircons which yield mean ²⁰⁶Pb/²³⁸U ages of 20.1 ± 0.2 and 19.8 ± 0.2 Ma. The ophiolite gabbro and the andesites both contain rare inherited zircons ranging from 122–104 Ma. The Early Cretaceous zircons in the arc andesites are interpreted as xenocrysts from the Mt. Camel basement terrane through which magmas of the Northland Miocene arc lavas erupted. The inherited zircons in the ophiolite gabbros suggest that a small fraction of this basement was introduced into the suboceanic mantle by subduction and mixed with mantle melts during ophiolite formation.

We postulate that the tholeiitic suite of the ophiolite represents the crustal segment of SSZ lithosphere (SSZL) generated in the southern South Fiji Basin (SFB) at a northeast-dipping subduction zone that was initiated at about 35 Ma. The subduction zone nucleated along a pre-existing transform boundary separating circa 45–20 Ma oceanic lithosphere to the north and west of the Northland Peninsula from nascent back arc basin lithosphere of the SFB. Construction of the SSZL propagated southward along the transform boundary as the SFB continued to unzip to the southeast. After subduction of a large portion of oceanic lithosphere by about 26 Ma and collision of the SSZL with New Zealand, compression between the Australian Plate and the Pacific Plate was taken up along a new southwest-dipping subduction zone behind the SSZL. Renewed volcanism began in the oceanic forearc at 25 Ma producing boninitic-like, SSZ and within-plate alkalic and calc-alkaline rocks. Rocks of these types temporally overlap ophiolite emplacement and subsequent Miocene continental arc construction.     

Continental subduction and exhumation of high-pressure rocks: insights from thermo-mechanical laboratory modelling

Thermo-mechanical physical modelling of continental subduction is performed to investigate the exhumation of deeply subducted continental crust. The model consists of two lithospheric plates made of new temperature sensitive analogue materials. The lithosphere is underlain by liquid asthenosphere. The continental lithosphere contains three layers: the weak sedimentary layer, the crust made of a stronger material, and of a still stronger lithospheric mantle. The whole model is subjected to a constant vertical thermal gradient, causing the strength reduction with depth in each lithospheric layer. Subduction is driven by both push force and pull force. During subduction, the subducting lithosphere is heating and the strength of its layers reduces. The weakening continental crust reaches maximal depth of about 120 km and cannot subduct deeper because its frontal part starts to flow up. The subducted crust undergoes complex deformation, including indicated upward ductile flow of the most deeply subducted portions and localised failure of the subducted upper crust at about 50-km depth. This failure results in the formation of the first crustal slice which rises up between the plates under the buoyancy force. This process is accompanied by the delamination of the crustal and mantle layers of the subducting lithosphere. The delamination front propagates upwards into the interplate zone resulting in the formation of two other crustal slices that also rise up between the plates. Average equivalent exhumation rate of the crustal material during delamination is about 1 cm/year. The crust-asthenosphere boundary near the interplate zone is uplifted. The subducted mantle layer then breaks off, removing the pull force and thereby stopping the delamination and increasing horizontal compression of the lithosphere. The latter produces shortening of the formed orogen and the growth of relief. The modelling reveals an interesting burial/exhumation evolution of the sedimentary cover. During initial stages of continental subduction the sediments of the continental margin are dragged to the overriding plate base and are partially accreted at the deep part of the interplate zone (at 60-70 km-depth). These sediments remain there until the beginning of delamination during which the pressure between the subducted crust and the overriding plate increases. This results in squeezing the underplated sediments out. Part of them is extruded upwards along the interplate zone to about 30-km depth at an equivalent rate of 5-10 cm/year.     

华北克拉通东部地壳与地幔盖层结构——长观测距深地震测深剖面结果

利用文登—阿拉善左旗长观测距地震宽角反射/折射剖面东段资料,辩识出4组地壳震相和3组地幔盖层震相.采用二维射线追踪走时反演和正演拟合交替计算方法,得到了包括鲁东隆起和华北裂陷盆地在内的地壳和地幔盖层二维速度结构.研究结果表明:华北裂陷盆地基底深达6km以上,研究区壳内界面C1埋深约15km,C2界面深约25km,Moho面平均埋深约35km.上地壳速度6.0~6.1km·s-1,且横向变化较大;中地壳速度相对均匀约为6.2~6.4km·s-1;下地壳速度为6.5~7.0km·s-1,速度梯度较大.地壳平均速度与隆起和坳陷构造相关.研究区岩石圈底界面一般为75~80km,西端接近太行隆起构造时深至90km左右,向西呈明显加深趋势,地壳厚度呈现相同的增厚特征.地幔盖层上部速度8.0~8.2km·s-1,具明显正梯度特征.岩石圈平均速度在郯庐断裂带附近显著偏低.PmP和PLP震相存在不同程度的复杂性,意味着在本地区Moho界面和岩石圈界面有较为复杂的结构,可能具有一定厚度或过渡带性质.结合其他研究结果认为,地幔盖层和下地壳速度梯度、界面性质差异与华北克拉通破坏相关,意味着破坏是一个渐变、缓慢和不均匀的过程.郯庐断裂带附近的低速应是其为软弱带的证据.     

两相矿物蠕变:探讨俯冲带深源地震机制

本文利用数值模拟方法研究了橄榄石-单斜辉石、瓦兹利石-单斜辉石、林伍德石-石榴石三种两相矿物模型在俯冲带温度应力条件下的蠕变,分析了剪切应力在俯冲带两相矿物模型中的发展,结果显示矿物蠕变差异在俯冲带剪切应力发展过程中起到了重要作用.尤其是林伍德石的快速蠕变率及其与石榴石的蠕变差极大地加速了俯冲带剪切应力的发展过程,可能是600km深度附近地震数量急剧增多的重要原因.我们的结果可以从一个新的角度理解俯冲带深源地震机制.     

Strain localisation in bimineralic rocks: Experimental deformation of synthetic calcite–anhydrite aggregates

Deformation of synthetic calcite–anhydrite aggregates to large shear strains (up to γ = 12.4 at 600 °C, 300 MPa confining pressure and a constant angular displacement rate corresponding to a shear strain rate of 10^− 3 s^− 1) resulted in the first experimental observation of strain localisation from initially homogeneous rocks. In contrast to experiments on pure calcite and anhydrite, which deformed homogeneously to large strains (γ ≥ 5), all experiments on calcite–anhydrite mixtures resulted in heterogeneous deformation at γ > 1 and the formation of narrow localised bands in the microstructures at γ > 4. In these bands, the amount of strain is at least twice as large as in the rest of the sample and individual grains of the same phase cluster and align, thereby forming microstructural layering similar to planar fabrics in natural mylonites. A switch in deformation mechanism in anhydrite from dislocation creep to diffusion creep and/or grain boundary sliding occurs simultaneously with strain localisation. It is concluded that deformation-induced heterogeneous phase distributions cause local strength differences initiating strain localisation in the calcite–anhydrite mixtures. The study suggests that the presence of two phases in combination with a change in deformation mechanism may be responsible for strain localisation in natural poly-mineralic mylonites.     

A mechanical model for deformation and earthquakes on strike-slip faults

A two-dimensional model for stress accumulation and earthquake instability associated with strike-slip faults is considered. The model consists of an elastic lithosphere overlying a viscous asthenosphere, and a fault of finite width with an upper brittle zone having an elastoplastic response and a lower ductile zone having an elastoviscoplastic response. For the brittle, or seismic, zone the behavior of the fault material is assumed to be governed by a relation which involves strain hardening followed by a softening regime, with strength increasing with depth. For the fault material in the ductile, or aseismic, section, the viscous effect is included through use of a nonlinear creep law, and the strength is assumed to decrease with depth. Hence, because of the lesser strength and the viscous effect, continuous flow occurs at great depths, causing stress accumulation at the upper portion of the fault and leading to failure at the bottom of the brittle zone. The failure is initially due to localized strain softening but, with further flow, the material above the softened zone reaches its maximum strength and begins to soften. This process accelerates and may result in an unstable upward rupture propagation.Relations are developed for the history of deformation within the lithosphere, specifically for the velocity of particles within the fault and at the ground surface. The boundary-element method is used for a quantitative study, and numerical results are obtained and compared with the recorded surface deformation of the San Andreas fault. The effects of geometry and material properties on instability, on the history of the surface deformation, and on the earthquake recurrence time are studied. The results are presented in terms of variations of ground-surface shear strain and shear strain rate, and velocity of points within the fault at various times during the earthquake cycle.It is found that the location of rupture initiation, the possibility of a sudden rupture as opposed to stable creep, and also the ground deformation pattern and its history, all critically depend on the mechanical response of the material within the fault zone, especially that of the brittle section. Shorter earthquake recurrence times are obtained for shallower brittle zones and for a stiffer lithosphere. Lower viscosities of the aseismic zone and the absence of asthenospheric coupling tend to suppress instability and promote stable creep. The model results thus suggest that the overall viscosity of the ductile creeping zone must exceed a minimum value for a sudden upward propagating rupture to take place within the seismic section.     

利用横波分裂到时差确定北美夏洛特皇后群岛岩石圈各向异性分层参数

本文收集到北美夏洛特皇后群岛VIB和DIB台站2005年1月1日至2014年12月31日2700余条近震三分量横波观测记录.经严格筛选提取到55条35°横波窗内的快慢波到时差值.通过建立到时差与多界面横波分裂路径关系,实现了对横波不同深度分裂界面的可靠识别和深度确定.依据对横波到时差及归一化到时差的分析结果,确认地处俯冲带的夏洛特皇后群岛的地震各向异性分别存在于该地区的上层陆壳、中层洋壳俯冲带及俯冲带底部或下层岩石圈顶部附近.发现除地幔岩石圈外,上两个分层的各向异性强度由北向南减弱,同时向西运动的洋壳俯冲带向南加深.不仅如此,归一化到时差结果显示该地区2012年M7.8强震前后各向异性强度减弱,但在3个各向异性层存在较大差异.无论在分层界面较浅的北部格雷厄姆岛附近还是在分层界面较深的南部莫尔兹比岛附近,强震后最上层的陆壳和最下层地幔岩石圈的各向异性强度均没有变化.M7.8后各向异性减弱全部发生在中层洋壳俯冲带内.并且,距离强震震中越近各向异性强度减弱越大,在一定程度上揭示出强震前后应力变化的空间分布特征.表明该地区洋壳俯冲是引起M7.8强震和现今构造运动的主要动力源.     

汕头－吕宋岛岩石圈速度结构与震源构造──台湾浅滩7．3级地震构造之一

汕头－吕宋岛岩石圈速度结构剖面，划分出华南陆缘古生代陆壳、陆架区晚古生代－中生代陆壳、陆坡带中生代－早第三纪过渡壳、新生代南海海盆洋壳及吕宋岛中生代－新生代岛弧陆壳与东吕宋海槽洋壳等地壳构造组分，并确定了上述地壳构造之间的边界断裂构造及其性质。结合地震震源分布及机制，初步确定了华南陆架盆岭构造带北、南两侧地震构造的控震构造与发震构造性质及其震源力学特征；１）指出１９９４年９月１６日台湾浅滩７．３级地震属于板缘壳幔地震及造成一千公里有感范围的原因；２）马尼拉海沟的海底地堑构造与南海海盆岩石圈地幔上隆是马尼拉海沟俯冲带震源显示正断层性质的原因，且为被动的或转换俯冲带；３）东吕宋海槽仍属于菲律宾海俯冲带性质；吕宋岛东西两侧俯冲带岩石圈板片震源深度的准三层分布，可能表明俯冲带岩石圈板片存在相应的低速滑移层。     

The generation of plate tectonics from mantle convection

In the last decade, significant progress has been made toward understanding how plate tectonics is generated from mantle dynamics. A primary goal of plate-generation studies has been the development of models that allow the top cold thermal boundary layer of mantle convection, i.e. the lithosphere, to develop broad and strong plate-like segments separated by narrow, weak and rapidly deforming boundaries; ideally, such models also permit significant strike-slip (toroidal) motion, passive ridges (i.e. pulled rather than pried apart), and self-consistent initiation of subduction. A major outcome of work so far is that nearly all aspects of plate generation require lithospheric rheologies and shear-localizing feedback mechanisms that are considerably more exotic than rheologies typically used in simple fluid-dynamical models of mantle flow. The search for plate-generating behavior has taken us through investigations of the effects of shear weakening (‘stick-slip’) and viscoplastic rheologies, of melting at ridges and low-viscosity asthenospheres, and of grain-size dependent rheologies and damage mechanics. Many such mechanisms, either by themselves or in combination, have led to self-consistent fluid-mechanical models of mantle flow that are remarkably plate-like, which is in itself a major accomplishment. However, many other important problems remain unsolved, such as subduction intiation and asymmetry, temporal evolution of plate geometry, rapid changes in plate motion, and the Archaean initiation of the plate-tectonic mode of convection. This paper presents a brief review of progress made in the plate-generation problem over the last decade, and discusses unresolved issues and future directions of research in this important area.     

High concentration of solutes at the upper part of the saturated zone (water table) of a deep aquifer under sewage-irrigated land

The water-table region (upper 50 cm of the saturated zone) of a 25 m deep phreatic sandstone aquifer, lying under fields irrigated with sewage effluents for up to 22 yrs, was monitored in 1971 and 1984. Average concentrations of NO⁻₃, Cl⁻ and SO²⁻₄ of up to 225, 307 and 155 mg l⁻¹, respectively, were detected in the upper 50 cm of the saturated region in two research wells in 1984. These concentrations, which are related to effluent and fertilizer input to groundwater, were two to four times higher than those found deep (37–55 m) below the water table in nearby (1000 m distant) production wells. Nitrate data and the estimated transit time through the unsaturated zone (2 m yr⁻¹) support the model suggesting that the major source of nitrate pollution in the past should be related to the oxidation of soil organic matter. The SO²⁻₄/Cl⁻ ratio is found to be a useful indicator for the arrival of SO²⁻₄-fertilizers at the groundwater interface. The observations presented in this paper question the suitability of plans for using effluents as a water source for agriculture in regions which are the replenishment areas of phreatic aquifers.     

俯冲带耦合作用对苏门答腊地区应变场影响的三维数值模拟

采用有限元方法模拟了俯冲带耦合作用对巽他弧及其邻区的影响.根据模拟结果,对比GPS、地震和地质学观测数据,定量分析了苏门答腊及其周边地区的应变强度和主应变方向的分布特征,据此探讨了该区构造特征、地震发生模式与耦合面积之间的关系.模型由具有黏弹性性质的岩石圈和软流圈上地幔组成,其中岩石圈包括了大陆岩石圈和大洋岩石圈以及俯冲至上地幔中的俯冲板片.研究结果如下:(1) 通过对不同俯冲带耦合面积模拟,发现苏门答腊前弧伴随耦合面积的增加应变强度逐渐增大,而增大的应变强度又影响了其周边地区的应变分布特征,因此整个苏门答腊前弧呈现出明显的分段性,这与该区地震破裂模式有较好的对应.(2)苏门答腊北部地区主应变方向与南部相比存在一定的差异,该差异是俯冲带的俯冲方向、俯冲速度、俯冲形态以及不同区域间耦合面积共同作用的结果.(3)虽然苏门答腊2004年地震主震区处于弱耦合状态,但从本文模拟的结果中可以看到,在俯冲作用下该区依然存在垂直向下的位移,这为地震激发海啸提供了有利的构造环境.     

Collision, subduction and accretion events in the Philippines: A synthesis

Abstract The Philippines preserves evidence of the superimposition of tectonic processes in ancient and present‐day collision and subduction zone complexes. The Baguio District in northern Luzon, the Palawan–Central Philippine region and the Mati–Pujada area in southeastern Mindanao resulted from events related to subduction polarity reversal leading to trench initiation, continent‐arc collision and autochthonous oceanic lithosphere emplacement, respectively. Geological data on the Baguio District in Northern Luzon reveal an Early Miocene trench initiation for the east‐dipping Manila Trench. This followed the Late Oligocene cessation of subduction along the west‐dipping proto‐East Luzon Trough. The Manila Trench initiation, which is modeled as a consequence of the counter‐clockwise rotation of Luzon, is attributed to the collision of the Palawan microcontinental block with the Philippine Mobile Belt. In the course of rotation, Luzon onramped the South China Sea crust, effectively converting the shear zone that bounded them into a subduction zone. Several collision‐related accretionary complexes (e.g. Romblon, Mindoro) are present in the Palawan–Central Philippine region. The easternmost collision zone boundary is located east of the Romblon group of islands. The Early Miocene southwestward shift of the collision boundary from Romblon to Mindoro started to end by the Pliocene. Continuous interaction between the Palawan microcontinental block and the Philippine Mobile Belt is presently taken up again along the collisional boundary east of the Romblon group of islands. The Mati–Pujada Peninsula area, on the other hand, is underlain by the Upper Cretaceous Pujada Ophiolite. This supra‐subduction zone ophiolite is capped by chert and pelagic limestones which suggests its derivation from a relatively deep marginal basin. The Pujada Ophiolite could be a part of a proto‐Molucca Sea plate. The re‐interpretation of the geology and tectonic settings of the three areas reaffirm the complex geodynamic evolution of the Philippine archipelago and addresses some of its perceived geological enigmas.     

A comparative study of results obtained in magnetotelluric deep soundings in Villarrica active volcano zone (Chile) with gravity investigations, distribution of earthquake foci, heat flow empirical relationships, isotopic geochemistry Sr/Sr and SB systematics

The first magnetotelluric deep soundings in Chile were carried out during 1986 in the Villarrica active volcano zone (39°25′S, 71°57′W). In the TM mode of polarization, the curves show a distorted segment with dispersion. A static distortion at long periods is observed in curves in the TE mode of polarization; the segment was shifted vertically to fit the geomagnetic global model values at daily periods. This modified curve was used for 1D modelling to determine the electrical structure in the study area. The upper level of the intermediate conducting layer of resistivity 20–60 ω m is found to be at 35–50 km depth. A higher resistivity layer (600 ω m), starting at 100 km depth, may be resolved in the intermediate conducting layer. A sharp decrease in the resistivity is shown by the model at 500 km.

Large heterogeneities at the level of the conducting layer encountered in the 1D modelling, and increased resistivity of the ultimate layer, may account for distortion observed at long periods. Two-dimensional test models show that the conducting layer in the area of Villarrica volcano may be an anomalous heated layer surrounded by rocks of higher resistivity of about 2 × 10³ ω m. These features correspond to the interaction with a subsiding oceanic lithosphere resulting in a complex thermal structure and perturbed resistivity distribution in transition zones of the Pacific type such as Chile, and to the existence of a megafault and a system of fractures in the sounding area. These facts make it difficult to determine the conductance of the electrical asthenosphere.

The parameters of the model structure correlate well with geophysical and geochemical results obtained in the area by other workers. Gravity studies indicate a maximum crustal thickness of about 37 km, which implies a non-full compensation according to the Airy hypothesis. The morphology of the Wadati-Benioff zone clearly shows a sharp decrease of earthquake foci at 50 km depth, and a reinitiation of seismic activity from a depth of 100 km down to 160 km. Kaufman and Keller and Levi and Lysak empirical relationships between heat flow and thickness of the lithosphere are somewhat consistent when the parameters of the model structure are assumed. The strontium isotope ratios are indicative of minimal crustal contamination of mantle-derived magmas, thus allowing a relatively thin crust under the area. Furthermore, the SB index of partial melting of mantle peridotite may suggest the contribution of uprising material from the deep asthenosphere through fractures extending to depth.     

Dynamics of thinning and destruction of the continental cratonic lithosphere: Numerical modeling

Thinning of the cratonic lithosphere is common in nature, but its destruction is not. In either case, the mechanisms for both thinning and destruction are still widely under debate. In this study, we have made a review on the processes and mechanisms of thinning and destruction of cratonic lithosphere according to previous studies of geological/geophysical observations and numerical simulations, with specific application to the North China Craton (NCC). Two main models are suggested for the thinning and destruction of the NCC, both of which are related to subduction of the oceanic lithosphere. One is the “bottom-up” model, in which the deeply subducting slab perturbs and induces upwelling from the hydrous mantle transition zone (MTZ). The upwelling produces mantle convection and erodes the bottom of the overriding lithosphere by the fluid-melt-peridotite reaction. Mineral compositions and rheological properties of the overriding lithospheric mantle are changed, allowing downward dripping of lithospheric components into the asthenosphere. Consequently, lithospheric thinning or even destruction occurs. The other is the “top-down” model, characterized by the flat subduction of oceanic slab beneath the overriding cratonic lithosphere. Dehydration reactions from the subducting slab would significantly hydrate the lithospheric mantle and decrease its rheological strength. Then the subduction angle may be changed from shallow to steep, inducing lateral upwelling of the asthenosphere. This upwelling would heat and weaken the overriding lithospheric mantle, which led to the weakened lithospheric mantle dripping into the asthenosphere. These two models have some similarities, in that both take the subducting oceanic slab and relevant fluid migration as the major driving mechanism for thinning or destruction of the overriding cratonic lithosphere. The key difference between the two models is the effective depth of the subducting oceanic slab. One is stagnation and flattening in the MTZ, whereas the other is flat subduction at the bottom of the cratonic lithosphere. In the NCC, the eastern lithosphere was likely affected by subduction of the Izanagi slab during the Mesozoic, which would have perturbed the asthenosphere and the MTZ, and induced fluid migration beneath the NCC lithosphere. The upwelling fluid may largely have controlled the reworking of the NCC lithosphere. In order to discuss and analyze these two models further, it is crucial to understand the role of fluids in the subduction zone and the MTZ. Here, we systematically discuss phase transformations of hydrous minerals and the transport processes of water in the subduction system. Furthermore, we analyze possible modes of fluid activity and the problems to explore the applied feasibility of each model. In order to achieve a comprehensive understanding of the mechanisms for thinning and destruction of cratonic lithosphere, we also consider four additional possible dynamic models: extension-induced lithospheric thinning, compression-induced lithospheric thickening and delamination, large-scale mantle convection and thermal erosion, and mantle plume erosion. Compared to the subduction-related models presented here, these four models are primarily controlled by the relatively simple and single process and mechanism (extension, compression, convection, and mantle plume, respectively), which could be the secondary driving mechanisms for the thinning and destruction of lithosphere.     

流变边界层及其对华北克拉通热/地震岩石圈底界差异的意义

二维热传导/对流数值模型显示,纯传导的固体岩石圈与纯对流的流体软流圈之间存在一过渡层,即流变边界层,其间传导与对流共同作用来传递热量.流变边界层厚度主要由软流圈黏性系数（η）控制,而受固体岩石圈厚度及热状态影响很小.随着η从1×10²¹Pa·s降低至1×10¹⁹ Pa·s,流变边界层也随之减薄,流变边界层的厚度与lg（η）成正比. 流变边界层的存在是造成热岩石圈与地震岩石圈厚度差异的重要因素. 全球典型克拉通岩石圈的对比结果表明,地震岩石圈厚度普遍大于热岩石圈厚度,二者的差异多数在70~90 km,很好地验证了流变边界层的存在. 研究发现二者的差异在华北克拉通自西向东逐渐减小：由西部鄂尔多斯的约80 km减少至渤海湾盆地的约20 km. 反映出华北克拉通岩石圈下部流变边界层厚度自西向东减薄,意味着软流圈黏性系数自西向东逐渐降低.这可能与中生代太平洋俯冲脱水形成的低黏大地幔楔有关,从一侧面印证了太平洋俯冲对华北克拉通破坏的影响.     

中下地壳温度压力条件下岩石半脆性蠕变的实验研究

本文着重介绍了在中下地壳温度、压力条件下进行的几种岩石半脆性蠕变的实验研究结果。通过对这些岩石的蠕变曲线、蠕变类型、蠕变速度和时间、蠕变应变的关系、稳态蠕变速度-流动应力和蠕变破坏特征的观测、分析,阐述了地壳岩石半脆性蠕变的基本规律,并在温度、围压-流动应力座标系中给出了岩石半脆性蠕变的分区