Seismogenesis of clustered seismicity beneath the Kangra–Chamba sector of northwest Himalaya: Constraints from 3D local earthquake tomography

To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km² centered NE to the epicenter of Kangra earthquake.     

Strength and stability of frictional sliding of gabbro gouge at elevated temperatures

To investigate the strength of frictional sliding and stability of mafic lower crust, we conducted experiments on oven-dried gabbro gouge of 1 mm thick sandwiched between country rock pieces (with gouge inclined 35° to the sample axis) at slip rates of 1.22 × 10^− 3 mm/s and 1.22 × 10^− 4 mm/s and elevated temperatures up to 615 °C. Special attention has been paid to whether transition from velocity weakening to velocity strengthening occurs due to the elevation of temperature.Two series of experiments were conducted with normal stresses of 200 MPa and 300 MPa, respectively. For both normal stresses, the friction strengths are comparable at least up to 510 °C, with no significant weakening effect of increasing temperature. Comparison of our results with Byerlee's rule on a strike slip fault with a specific temperature profile in the Zhangbei region of North China shows that the strength given by experiments are around that given by Byerlee's rule and a little greater in the high temperature range.At 200 MPa normal stress, the steady-state rate dependence a − b shows only positive values, probably still in the “run-in” process where velocity strengthening is a common feature. With a normal stress of 300 MPa, the values of steady-state rate dependence decreases systematically with increasing temperature, and stick-slip occurred at 615 °C. Considering the limited displacement, limited normal stress applied and the effect of normal stress for the temperatures above 420 °C, it is inferred here that velocity weakening may be the typical behaviour at higher normal stress for temperature above 420 °C and at least up to 615 °C, which covers most of the temperature range in the lower crust of geologically stable continental interior. For a dry mafic lower crust in cool continental interiors where frictional sliding prevails over plastic flow, unstable slip nucleation may occur to generate earthquakes.     

3D modeling of earthquake cycles of the Xianshuihe fault,southwestern China

We perform 3D modeling of earthquake generation of the Xianshuihe fault, southwestern China, which is a highly active strike-slip fault with a length of about 350 km, in order to understand earthquake cycles and segmentations for a long-term forecasting and earthquake nucleation process for a short-term forecasting. Historical earthquake data over the last 300 years indicates repeated periods of seismic activity, and migration of large earthquake along the fault during active seismic periods. To develop the 3D model of earthquake cycles along the Xianshuihe fault, we use a rate- and state-dependent friction law. After analyzing the result, we find that the earthquakes occur in the reoccurrence intervals of 400–500 years. Simulation result of slip velocity distribution along the fault at the depth of 10 km during 2694 years along the Xianshuihe fault indicates that since the third earthquake cycle, the fault has been divided into 3 parts. Some earthquake ruptures terminate at the bending part of the fault line, which may means the shape of the fault line controls how earthquake ruptures. The change of slip velocity and displacement at 10 km depth is more tremendous than the change of the shallow and deep part of the fault and the largest slip velocity occurs at the depth of 10 km which is the exact depth of the seismic zone where fast rupture occurs.     

Dynamic process analysis for the initiation and movement of the Donghekou landslide-debris flow triggered by the Wenchuan earthquake

The Donghekou landslide-debris flow was a remarkable geological disaster triggered by the Wenchuan earthquake in 2008. The dynamic process of a rapid landslide-debris flow is very complicated and can be divided into two aspects: the slope dynamic response of the earthquake and the mass movement and accumulation process. A numerical method combined with a finite difference method (FDM) and discrete element method (DEM) for simulation of landslide-debris flow under seismic loading is presented. The FDM and DEM are coupled through the critical sliding surface, initiation time and velocity. The dynamic response of the slope is simulated by the finite difference method, and critical sliding surface is determined using the earthquake response spectrum method. The landslide initiation time and the velocity are determined by time–history analysis. The mass movement and accumulation process is simulated using the discrete element method. Simulation results demonstrate that the maximum amplification coefficient of dynamic acceleration for the Donghekou slope is approximately 3.909, the initiation time of landslide is approximately 6.0 s, and the average initial velocity of the sliding mass is approximately 0.85 m/s. The failure of the slope is the result of elevation-orientated amplification effect and the sliding mass triggered with a small initial velocity. The numerical simulated result of the maximum sliding velocity is approximately 66.35 m/s, and the mass is disintegrated rapidly because of collision and free fall. The landslide velocity decreases when the flowing mass reaches a lower slope angle and gradually comes to a stop, and the total travel distance is approximately 2400 m.     

Frictional properties of saponite-rich gouge from a serpentinite-bearing fault zone along the Gokasho-Arashima Tectonic Line,central Japan

We studied a serpentinite-bearing fault zone in Gokasho-Arashima Tectonic Line, Mie Prefecture, central Japan, characterizing its internal structures, mineral assemblage, permeability, and frictional properties. The fault core situated between the serpentinite breccia and the adjacent sedimentary rocks is characterized by a zone locally altered to saponite. The clayey gouge layer separates fault rocks of serpentinite origin containing talc and tremolite from fault rocks of sedimentary origin containing chlorite but no quartz. The minerals that formed within the fault are the products of metasomatic reaction between the serpentinite and the siliceous rocks. Permeability measurements show that serpentinite breccia and fault gouge have permeability of 10⁻¹⁴–10⁻¹⁷ m² and 10⁻¹⁵–10⁻¹⁸ m², respectively, at 5–120 MPa confining pressure. Frictional coefficient of the saponite-rich clayey fault gouge ranged between 0.20 and 0.35 under room-dry condition, but was reduced to 0.06–0.12 when saturated with water. The velocity dependence of friction was strongly positive, mostly ranging between 0.005 and 0.006 in terms of a–b values. The governing friction law is not constrained yet, but we find that the saponite-rich gouge possesses an evolutional behavior in the opposite direction to that suggested by the rate and state friction law, in addition to its direct velocity dependence.     

Heterogeneous bubble nucleation in highly viscous silicate melts during instantaneous decompression from high pressure

The nucleation of H₂O bubbles in magmas has been proposed as a trigger for volcanic eruptions. To determine how bubbles nucleate heterogeneously in silicate melts, experiments were carried out in which high-silica rhyolitic melts were hydrated at 800 °C and either 50 or 125 MPa, and then decompressed by 20–91 MPa at temperatures that ranged between 550 and 700 °C, and held at the lower pressures for 10–720 s before being quenched. Bubbles nucleated in number densities (N_B) that vary between 3 × 10⁷ and 2 × 10⁸ cm^− 3. Blocky shaped magnetite or the ends of needle-shaped hematite acted as sites for nucleation, but only if a minimum super-saturation was exceeded, which increases with increasing melt viscosity. Bubbles did not nucleate along the lengths of hematite needles nor on plagioclase. Both the beginning and ending times of the nucleation event increases with increasing melt viscosity. Using nucleation theory predictions, neither the slower nucleation rates nor the changing activation of nucleation sites can be adequately explained by the differences in temperature, water diffusivity, or viscosity. Instead, the variations in nucleation kinetics are best explained by changes in surface tension between melt and vapor, resulting from the increasing polymerization of the melt at lower temperatures and water contents. Because only ∼ 10⁸ bubbles cm^− 3 nucleate on magnetite in the rhyolite melt used, almost regardless of experimental conditions, results from this study may not be directly comparable to vesicle numbers in volcanic pumice of different compositions.     

Seismic imaging,crustal stress and GPS data analyses: Implications for the generation of the 2008 Wenchuan Earthquake (M7.9), China

Seismic imaging together with global positioning system (GPS) and crustal stress data analyses show that the Mw7.9 2008 Wenchuan earthquake occurred within a distinct area of high crustal stress (～ 17.5 MPa) and high Poisson's ratio (7–10%) anomalies centered on the Longmen-Shan (Shan means Mountain in Chinese) tectonic fault belt. Low P-wave and S-wave velocities in the southwest (SW) segment contrast with high-velocity anomalies in the central portion (CP) and northeast (NE) segment within the uppermost ～ 15 km depths along the tectonic fault belt, though a presumably ductile zone with low-velocity anomalies separates the CP and NE segment. The rupture initiated near the southwestern end of the CP at a zone of high Poisson's ratio (σ) which extends down into the lower crust. These low-velocity and high-σ anomalies immediately below the source hypocenter, together with the high crustal stress, indicate the presence of high-pressure fluids from the lower crust, which might have reduced the mechanical strength of the fractured rock matrix and triggered the earthquake. Our study suggests that the structural heterogeneity and high crustal stress played an important role in the nucleation of the Wenchuan earthquake and its rupture process.     

Tectonics and thermal structure of western Satpura,India

Increased seismicity and occurrences of hot springs having surface temperature of 36–58 °C are observed in the central part of India (74–81° E, 20–25° N), where the NE trending Middle Proterozoic Aravalli Mobile Belt meets the ENE trending Satpura Mobile Belt. Earlier Deep Seismic Sounding (DSS) studies along Thuadara-Sendhwa-Sindad profile in the area has showed Mesozoic Sediments up to around 4 km depth covered by Deccan Trap and the Moho depth with a boundary velocity (P_n) of 8.2 km/s. In the present study, surface heat flow of 48 ± 4 mW m^?2 has been estimated based on P_n velocity, which agrees with the value of heat flow of 52 ± 4 mW m^?2 based on Curie point isotherms estimates. The calculated temperature-depth profile shows temperature of 80–120 °C at the basement, which is equivalent to oil window temperature in Mesozoic sediments and around 570–635 °C at Moho depth of 38–43 km and the thermal lithosphere is about 110 km thick, which is comparatively higher than those of adjoining regions. The present study reveals the brittle–ductile transition zone at 14–41 km depth (temperature around 250–600 °C) where earthquake nucleation takes place.     

ESR dating of the Wangsan fault,South Korea

The Wangsan fault is exposed at Kyeongju, Korea. The andesite is unconformably covered by Quaternary alternating conglomerate and sandstone deposits. The unconformity is cut by the thrust fault which displaces a hanging wall block of about 30 m. Exposed at the surface of this fault is a light gray and brown fault gouge, about 40–120 cm thick. In order to test the consistency of ESR ages from a single fault gouge zone, we collected six gouge samples systematically along and across the Wangsan fault. We found that six samples collected from the same gouge zone show consistent ESR date estimates (average 550 ka). Because fault rock is rare along faults in unconsolidated sediments due to low confining stress near the surface, we consider that this fault gouge has been moved up with the hanging wall block along the fault. The estimated average uplift rate of the hanging wall block is about 0.04 cm/year based on the age of the displaced Quaternary deposits (vertical separation; about 20 m) dated by OSL dating methods. The depth of the fault gouge at the time of reactivation, which was estimated from uplift rate and the ESR ages, is about 220 m. Therefore, we conclude that the results of ESR age estimates represent the time of reactivation of the fault gouge at a depth of 220 m in the past during fault movement, because later movements which occurred during uplift near the surface may not have zeroed ESR signals significantly for ESR dating of fault movements.     

Fracture termination and step-over at bedding interfaces due to frictional slip and interface opening

Three types of fracture intersection with bedding contacts have been investigated within numerical experiments: fracture transection through bed contacts, termination (abutment) at contacts and step-over of fractures at bedding contacts. To evaluate the mechanisms responsible for different fracture intersections with bed contacts, the numerical experiments explored deformation associated with end-member conditions of sliding-only interfaces and opening-only interfaces. A third suite of models explored the combined influence of both sliding and opening, as a fracture approached the interface. In contrast to our initial supposition that interface sliding promotes fracture termination, the sliding-only interfaces encouraged propagation of fractures straight through the modeled interface. In contrast, the opening-only interfaces yielded either fracture termination or initiation of a new fracture near the ends of the open interface segment (several centimeters from parent fracture in these models). These results suggest that local interface opening near the tip of approaching fractures, rather than sliding, is responsible for fracture termination and step-over at bedding contacts. Combined sliding and opening yielded fracture termination in models with weak interfaces (μ=0; c=0 MPa; T=0 MPa) and either fracture step-over or termination at moderate-strength interfaces (μ=0.65; c=3.25 MPa; T=5 MPa). Fracture termination occurs at moderate-strength interfaces when the stresses along the interface are not great enough to initiate a new step-over fracture. Fracture termination is more likely under conditions of shallower burial depth, lower layer-parallel effective tension and fluid-driven fracture propagation rather than remote layer-parallel tension. Furthermore, thicker beds and greater layer-parallel effective tension may produce greater distances of fracture step-over than thinner beds and more compressive layers. These results may assist in the prediction of subsurface fracture networks and associated fluid flow paths.     

Mechanisms of continental subduction and exhumation of HP and UHP rocks

We discuss possible scenarios of continental collision, and their relation to mechanisms of exhumation of HP and UHP rocks, inferred from thermo-mechanical numerical models accounting for thermo-rheological complexity of the continental lithosphere. Due to this complexity, mechanisms of continental convergence are versatile and different, in many aspects from those that control oceanic subduction. Elucidating these mechanisms from conventional observations is difficult, and requires additional constraints such as those derived from petrological data. Indeed, exhumation of HP/UHP rocks is an integral part of convergent processes, and burial/exhumation dynamics inferred from metamorphic P–T–t paths provides strong constraints on the collision scenarios. Metamorphic rocks also play an active role due to their contrasting physical properties (rheology, density, fluid transport capacity). Numerical thermo-mechanical experiments suggest that HP/UHP exhumation can only be produced in subduction contexts, as well as that long-lasting (> 10 Myr) continental subduction can only occur in case of cold strong lithospheres (T_Moho < 550 °C, the equivalent elastic thickness T_e > 50 km) and of relatively high convergence rates (> 3–5 cm yr^− 1 ). In this case, high density UHP material in the crustal part of subduction interface provides additional pull on the slab and is not always exhumed to the surface. In case of slower convergence and/or weaker lithosphere (T_e < 40 km), continental subduction is a transient process that takes a limited time span in the evolution of collision zone. Under these conditions, hot mechanically weak UHP rocks enhance decoupling between the upper and lower plate while their exhumation may be rapid (faster than convergence rate) and abundant. Therefore, the UHP exhumation paths can be regarded as sensitive indicators of subduction. Rheological changes and fluid exchanges associated with low-to-middle pressure phase transitions along the subduction interface, such as serpentinization during the oceanic phase and schisting, play a major role producing necessarily mechanical softening of the subduction interface and of the hydrated mantle wedge. The oceanic UHP rocks are exhumed thanks to mixing with low-density continental crustal units during transition from oceanic to continental subduction. At the continental phase, the UHP exhumation occurs as a result of a multi-stage process: at the deep stage (< 40 km depth) the exhumation is rapid and is driven by buoyancy of partly metamorphosed (or partly molten) UHP material often mixed with non-metamorphosed crustal volumes. At final stages, exhumation takes common slow path through the accretion prism mechanism and the erosional denudation. The experiments suggest that formation of UHP rocks requires that continental subduction starts at higher oceanic subduction rate. It then may progressively slow down until the lockup of the subduction interface and/or slab-break-off. A rate of ~ 1–2 cm yr^− 1 is generally sufficient to drive continental subduction during the first several Myr of convergence, but pertinent subduction requires faster convergence rates (> 3–5 cm yr^− 1). We suggest that most continental orogenic belts could have started their formation from continental subduction but this process has been generally limited in time.     

Regimes of subduction and lithospheric dynamics in the Precambrian: 3D thermomechanical modelling

Comparing the early Earth to the present day, geological–geochemical evidence points towards higher mantle potential temperature and a different type of tectonics. In order to investigate possible changes in Precambrian tectonic styles, we conduct 3D high-resolution petrological–thermomechanical numerical modelling experiments for oceanic plate subduction under an active continental margin at a wide range of mantle potential temperature T_P (∆ T_P = 0 − 250 K, compared to present day conditions). At present day mantle temperatures (∆ T_P = 0 K), results of numerical experiments correspond to modern-style subduction, whereas at higher temperature conditions important systematic changes in the styles of both lithospheric deformation and mantle convection occur. For ∆ T_P = 50 − 100 K a regime of dripping subduction emerges which is still very similar to present day subduction but is characterised by frequent dripping from the slab tip and a loss of coherence of the slab, which suggests a close relationship between dripping subduction and episodic subduction. At further increasing ∆ T_P = 150 − 200 K dripping subduction is observed together with unstable dripping lithosphere, which corresponds to a transitional regime. For ∆ T_P = 250 K, presumably equivalent to early Archean, the dominating tectonic style is characterised by small-scale mantle convection, unstable dripping lithosphere, thick basaltic crust and small plates. Even though the initial setup is still defined by present day subduction, this final regime shows many characteristics of plume-lid tectonics. Transition between the two end-members, plume-lid tectonics and plate tectonics, happens gradually and at intermediate temperatures elements of both tectonic regimes are present. We conclude, therefore, that most likely no abrupt geodynamic regime transition point can be specified in the Earth's history and its global geodynamic regime gradually evolved over time from plume-lid tectonics into modern style plate tectonics.     

U-series dating of co-seismic gypsum and submarine paleoseismology of active faults in Northern Chile (23°S)

The convergence of the Nazca and South American plates along the subduction margin of the central Andes results in large subduction earthquakes and tectonic activity along major fault systems. Despite its relevance, the paleoseismic record of this region is scarce, hampering our understanding about the relationship between the Andes building and earthquake occurrence. In this study, we used the U-series disequilibrium method to obtain absolute ages of paleoearthquake events associated with normal displacements along the active Mejillones and Salar del Carmen faults in the Coastal Range of the Atacama Desert of northern Chile. The ²³⁰Th–²³⁴U disequilibrium ages in co-seismic gypsum salts sampled along the fault traces together with marine evidences indicate that earthquakes occurred at ca. 29.7 ± 1.7 ka, 11 ± 4 ka and 2.4 ± 0.8 ka. When coupled with paleoseismic marine and radiocarbon (¹⁴C) records in the nearby Mejillones Bay evidencing large dislocations along the Mejillones Fault, the geochronological dataset presented here is consistent with the notion that gypsum salts formed during large earthquakes as a result of co-seismic dilatancy pumping of saline waters along the major faults. Based on maximum observed cumulative vertical offsets in the studied faults, this phenomena could have occurred episodically at a rate in the order of 1:40 to 1:50 with respect to the very large subduction earthquakes during the latest Pleistocene–Holocene period. The results presented here reveal that the U-series disequilibrium method can be successfully applied to date the gypsum salts deposited along faults during seismic events, and therefore directly constrain the age of large paleoearthquakes in hyperarid and seismically active zones.     

The large aftershocks triggered by the 2011 Mw 9.0 Tohoku-Oki earthquake,Japan

The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography.     

Shear localization,velocity weakening behavior,and development of cataclastic foliation in experimental granite gouge

Microstructural aspects of room-temperature deformation in experimental Westerly granite gouge were studied by a set of velocity stepping rotary-shear experiments at 25 MPa normal stress. The experiments were terminated at: (a) 44 mm, (b) 79 mm, and (c) 387 mm of sliding, all involving variable-amplitude fluctuations in friction. Microstructural attributes of the gouge were studied using scanning (SEM) and scanning transmission electron microscopy (STEM), image processing, and energy dispersive X-ray (EDX) analyses. The gouge was velocity weakening at sliding distances >10 mm as a core of cataclasites along a through-going shear zone developed within a mantle of less deformed gouge in all experiments. Unlike in experiment (a), the cataclasites in experiments (b) and (c) progressively developed a foliation defined by stacks of shear bands. The individual bands showed an asymmetric particle-size grading normal to shearing direction. These microstructures were subsequently disrupted and reworked by high-angle Riedel shears. While the microstructural evolution affected the effective thickness and frictional strength of the gouge, it did not affect its overall velocity dependence behavior. We suggest that the foliation resulted from competing shear localization and frictional slip hardening and that the velocity dependence of natural fault gouge depends upon compositional as well as microstructural evolution of the gouge.     

Observed source parameters for dynamic rupture with non-uniform initial stress and relatively high fracture energy

We have conducted dynamic rupture propagation experiments to establish the relations between in-source stress drop, fracture energy and the resulting particle velocity during slip of an unconfined 2 m long laboratory fault at normal stresses between 4 and 8 MPa. To produce high fracture energy in the source we use a rough fault that has a large slip weakening distance. An artifact of the high fracture energy is that the nucleation zone is large such that precursory slip reduces fault strength over a large fraction of the total fault length prior to dynamic rupture, making the initial stress non-uniform. Shear stress, particle velocity, fault slip and acceleration were recorded coseismically at multiple locations along strike and at small fault-normal distances. Stress drop increases weakly with normal stress. Average slip rate depends linearly on the fault strength loss and on static stress drop, both with a nonzero intercept. A minimum fracture energy of 1.8 J/m² and a linear slip weakening distance of 33 μm are inferred from the intercept. The large slip weakening distance also affects the average slip rate which is reduced by in-source energy dissipation from on-fault fracture energy.Because of the low normal stress and small per event slip (∼86 μm), no thermal weakening such as melting or pore fluid pressurization occurs in these experiments. Despite the relatively high fracture energy, and the very low heat production, energy partitioning during these laboratory earthquakes is very similar to typical earthquake source properties. The product of fracture energy and fault area is larger than the radiated energy. Seismic efficiency is low at ∼2%. The ratio of apparent stress to static stress drop is ∼27%, consistent with measured overshoot. The fracture efficiency is ∼33%. The static and dynamic stress drops when extrapolated to crustal stresses are 2–7.3 MPa and in the range of typical earthquake stress drops. As the relatively high fracture energy reduces the slip velocities in these experiments, the extrapolated average particle velocities for crustal stresses are 0.18–0.6 m/s. That these experiments are consistent with typical earthquake source properties suggests, albeit indirectly, that thermal weakening mechanisms such as thermal pressurization and melting which lead to near complete stress drops, dominate earthquake source properties only for exceptional events unless crustal stresses are low.     

Early Mesozoic southward subduction history of the Mongol–Okhotsk oceanic plate: Evidence from geochronology and geochemistry of Early Mesozoic intrusive rocks in the Erguna Massif,NE China

In this paper we present new zircon U–Pb ages, Hf isotope data, and whole-rock major and trace element data for Early Mesozoic intrusive rocks in the Erguna Massif of NE China, and we use these data to constrain the history of southward subduction of the Mongol–Okhotsk oceanic plate, and its influence on NE China as a whole. The zircon U–Pb dating indicates that Early Mesozoic magmatic activity in the Erguna Massif can be subdivided into four stages at ~ 246 Ma, ~ 225 Ma, ~ 205 Ma, and ~ 185 Ma. The ~ 246 Ma intrusive rocks comprise a suite of high-K calc-alkaline diorites, quartz diorites, granodiorites, monzogranites, and syenogranites, with I-type affinities. The ~ 225 Ma intrusive rocks consist of gabbro–diorites and granitoids, and they constitute a bimodal igneous association. The ~ 205 Ma intrusive rocks are dominated by calc-alkaline I-type granitoids that are accompanied by subordinate intermediate–mafic rocks. The ~ 185 Ma intrusive rocks are dominated by I-type granitoids, accompanied by minor amounts of A-types. These Early Mesozoic granitoids mainly originated by partial melting of a depleted and heterogeneous lower crust, whereas the coeval mafic rocks were probably derived from partial melting of a depleted mantle modified by subduction-related fluids. The rock associations and their geochemical features indicate that the ~ 246 Ma, ~ 205 Ma, and ~ 185 Ma intrusive rocks formed in an active continental margin setting related to the southward subduction of the Mongol–Okhotsk oceanic plate. The ~ 225 Ma bimodal igneous rock association formed within an extensional environment in a pause during the subduction process of the Mongol–Okhotsk oceanic plate. Every magmatic stage has its own corresponding set of porphyry deposits in the southeast of the Mongol–Okhotsk suture belt. Taking all this into account, we conclude the following: (1) during the Early Mesozoic, the Mongol–Okhotsk oceanic plate was subducted towards the south beneath the Erguna Massif, but with a pause in subduction at ~ 225 Ma; and (2) the southward subduction of the Mongol–Okhotsk oceanic plate not only caused the intense magmatic activity, but was also favorable to the formation of porphyry deposits.     

Exhumation of oceanic blueschists and eclogites in subduction zones: Timing and mechanisms

High-pressure low-temperature (HP–LT) metamorphic rocks provide invaluable constraints on the evolution of convergent zones. Based on a worldwide compilation of key information pertaining to fossil subduction zones (shape of exhumation P–T–t paths, exhumation velocities, timing of exhumation with respect to the convergence process, convergence velocities, volume of exhumed rocks,…), this contribution reappraises the burial and exhumation of oceanic blueschists and eclogites, which have received much less attention than continental ones during the last two decades.Whereas the buoyancy-driven exhumation of continental rocks proceeds at relatively fast rates at mantle depths (≥ cm/yr), oceanic exhumation velocities for HP–LT oceanic rocks, whether sedimentary or crustal, are usually on the order of the mm/yr. For the sediments, characterized by the continuity of the P–T conditions and the importance of accretionary processes, the driving exhumation mechanisms are underthrusting, detachment faulting and erosion. In contrast, blueschist and eclogite mafic bodies are systematically associated with serpentinites and/or a mechanically weak matrix and crop out in an internal position in the orogen.Oceanic crust rarely records P conditions > 2.0–2.3 GPa, which suggests the existence of maximum depths for the sampling of slab-derived oceanic crust. On the basis of natural observations and calculations of the net buoyancy of the oceanic crust, we conclude that beyond depths around 70 km there are either not enough serpentinites and/or they are not light enough to compensate the negative buoyancy of the crust.Most importantly, this survey demonstrates that short-lived (< ～ 15 My), discontinuous exhumation is the rule for the oceanic crust and associated mantle rocks: exhumation takes place either early (group 1: Franciscan, Chile), late (group 2: New Caledonia, W. Alps) or incidentally (group 3: SE Zagros, Himalayas, Andes, N. Cuba) during the subduction history. This discontinuous exhumation is likely permitted by the specific thermal regime following the onset of a young, warm subduction (group 1), by continental subduction (group 2) or by a major, geodynamic modification of convergence across the subduction zone (group 3; change of kinematics, subduction of asperities, etc).Understanding what controls this short-lived exhumation and the detachment and migration of oceanic crustal slices along the subduction channel will provide useful insights into the interplate mechanical coupling in subduction zones.     

A Late Cretaceous and Cenozoic reconstruction of the Southwest Pacific region: Tectonics controlled by subduction and slab rollback processes

A Cenozoic tectonic reconstruction is presented for the Southwest Pacific region located east of Australia. The reconstruction is constrained by large geological and geophysical datasets and recalculated rotation parameters for Pacific–Australia and Lord Howe Rise–Pacific relative plate motion. The reconstruction is based on a conceptual tectonic model in which the large-scale structures of the region are manifestations of slab rollback and backarc extension processes. The current paradigm proclaims that the southwestern Pacific plate boundary was a west-dipping subduction boundary only since the Middle Eocene. The new reconstruction provides kinematic evidence that this configuration was already established in the Late Cretaceous and Early Paleogene. From ∼ 82 to ∼ 52 Ma, subduction was primarily accomplished by east and northeast-directed rollback of the Pacific slab, accommodating opening of the New Caledonia, South Loyalty, Coral Sea and Pocklington backarc basins and partly accommodating spreading in the Tasman Sea. The total amount of east-directed rollback of the Pacific slab that took place from ∼ 82 Ma to ∼ 52 Ma is estimated to be at least 1200 km. A large percentage of this rollback accommodated opening of the South Loyalty Basin, a north–south trending backarc basin. It is estimated from kinematic and geological constraints that the east–west width of the basin was at least ∼ 750 km. The South Loyalty and Pocklington backarc basins were subducted in the Eocene to earliest Miocene along the newly formed New Caledonia and Pocklington subduction zones. This culminated in southwestward and southward obduction of ophiolites in New Caledonia, Northland and New Guinea in the latest Eocene to earliest Miocene. It is suggested that the formation of these new subduction zones was triggered by a change in Pacific–Australia relative motion at ∼ 50 Ma. Two additional phases of eastward rollback of the Pacific slab followed, one during opening of the South Fiji Basin and Norfolk Basin in the Oligocene to Early Miocene (up to ∼ 650 km of rollback), and one during opening of the Lau Basin in the latest Miocene to Present (up to ∼ 400 km of rollback). Two new subduction zones formed in the Miocene, the south-dipping Trobriand subduction zone along which the Solomon Sea backarc Basin subducted and the north-dipping New Britain–San Cristobal–New Hebrides subduction zone, along which the Solomon Sea backarc Basin subducted in the west and the North Loyalty–South Fiji backarc Basin and remnants of the South Loyalty–Santa Cruz backarc Basin subducted in the east. Clockwise rollback of the New Hebrides section resulted in formation of the North Fiji Basin. The reconstruction provides explanations for the formation of new subduction zones and for the initiation and termination of opening of the marginal basins by either initiation of subduction of buoyant lithosphere, a change in plate kinematics or slab–mantle interaction.     

Low b-value prior to the Indo-Myanmar subduction zone earthquakes and precursory swarm before the May 1995 M 6.3 earthquake

Some 455 events (m_b ⩾ 4.5) in the Indo-Myanmar subduction zone are compiled using the ISC/EHB/NEIC catalogues (1964–2011) for a systematic study of seismic precursors, b-value and swarm activity. Temporal variation of b-value is studied using the maximum likelihood method beside CUSUM algorithm. The b-values vary from 0.95 to 1.4 for the deeper (depth ⩾60 km) earthquakes, and from 0.85 to 1.3 for the shallower (depth <60 km) earthquakes. A sudden drop in the b-value, from 1.4 to 0.9, prior to the occurrence of larger earthquake(s) at the deeper depth is observed. It is also noted that the CUSUM gradient reversed before the occurrence of larger earthquakes. We further examined the seismicity pattern for the period 1988–1995 within a radius of 150 km around the epicentre (latitude: 24.96°N; longitude: 95.30°E) of a deeper event M 6.3 of May 6, 1995 in this subduction zone. A precursory swarm during January 1989 to July 1992 and quiescence during August 1992 to April 1995 are identified before this large earthquake. These observations are encouraging to monitor seismic precursors for the deeper events in this subduction zone.