Surface deformation related to the 2008 Wenchuan earthquake,and mountain building of the Longmen Shan,eastern Tibetan Plateau

The 12 May 2008 M_s 8.0 Wenchuan earthquake, China, was one of largest continental thrusting events worldwide. Based on interpretations of post-earthquake high-resolution remote sensing images and field surveys, we investigated the geometry, geomorphology, and kinematics of co-seismic surface ruptures, as well as seismic and geologic hazards along the Longmen Shan fold-and-thrust belt. Our results indicate that the Wenchuan earthquake occurred along the NE–SW-trending Yingxiu–Beichuan and Guanxian–Anxian faults in the Longmen Shan fold-and-thrust belt. The main surface rupture zones along the Yingxiu–Beichuan and Guanxian–Anxian fault zones are approximately 235 and 72 km in length, respectively. These sub-parallel ruptures may merge at depth. The Yingxiu–Donghekou surface rupture zone can be divided into four segments separated by discontinuities that appear as step-overs or bends in map view. Surface deformation is characterized by oblique reverse faulting with a maximum vertical displacement of approximately 10 m in areas around Beichuan County. Earthquake-related disasters (e.g., landslides) are linearly distributed along the surface rupture zones and associated river valleys.The Wenchuan earthquake provides new insights into the nature of mountain building within the Longmen Shan, eastern Tibetan Plateau. The total crustal shortening accommodated by this great earthquake was as much as 8.5 m, with a maximum vertical uplift of approximately 10 m. The present results suggest that ongoing mountain building of the Longmen Shan is driven mainly by crustal shortening and uplift related to repeated large seismic events such as the 2008 Wenchuan earthquake. Furthermore, rapid erosion within the Longmen Shan fold-and-thrust belt occurs along deep valleys and rupture zones following the occurrence of large-scale landslides triggered by earthquakes. Consequently, we suggest that crustal shortening related to repeated great seismic events, together with isostatic rebound induced by rapid erosion-related unloading, is a key component of the geodynamics that drive ongoing mountain building on the eastern Tibetan Plateau.     

The Zagros hinterland fold-and-thrust belt in-sequence thrusting,Iran

The collision of the Iranian microcontinent with the Afro-Arabian continent resulted in the deformation of the Zagros orogenic belt. The foreland of this belt in the Persian Gulf and Arabian platform has been investigated for its petroleum and gas resource potentials, but the Zagros hinterland is poorly investigated and our knowledge about its deformation is much less than other parts of this orogen. Therefore, this work presents a new geological map, stratigraphic column and two detailed geological cross sections. This study indicates the presence of a hinterland fold-and-thrust belt on northeastern side of the Zagros orogenic core that consists of in-sequence thrusting and basement involvement in this important part of the Zagros hinterland. The in-sequence thrusting resulted in first- and second-order duplex systems, Mode I fault-bend folding, fault-propagation folding and asymmetric detachment folding which indicate close relationships between folding and thrusting. Study of fault-bend folds shows that layer-parallel simple shear has the same role in the southeastern and northwestern parts of the study area (α_e = 23.4 ± 9.1°). A major lateral ramp in the basement beneath the Talaee plain with about one kilometer of vertical offset formed parallel to the SW movement direction and perpendicular to the major folding and thrusting.     

Plio-Quaternary stress states in NE Iran: Kopeh Dagh and Allah Dagh-Binalud mountain ranges

NE Iran, including the Kopeh Dagh and Allah Dagh-Binalud deformation domains, comprises the northeastern boundary of the Arabia–Eurasia collision zone. This study focuses on the evolution of the Plio-Quaternary tectonic regimes of northeast Iran. We present evidence for drastic temporal changes in the stress state by inversion of both geologically and seismically determined fault slip vectors. The inversions of fault kinematics data reveal distinct temporal changes in states of stress during the Plio-Quaternary (since ～ 5 Ma). The paleostress state is characterized by a regional transpressional tectonic regime with a mean N140 ± 10°E trending horizontal maximum stress axis (σ₁). The youngest (modern) state of stress shows two distinct strike-slip and compressional tectonic regimes with a regional mean of N030 ± 15°E trending horizontal σ₁. The change from the paleostress to modern stress states has occurred through an intermediate stress field characterized by a mean regional N trending σ₁. The inversion analysis of earthquake focal mechanisms reveals a homogeneous, transpressional tectonic regime with a regional N023 ± 5°E trending σ₁. The modern stress state, deduced from the youngest fault kinematics data, is in close agreement with the present-day stress state given by the inversions of earthquake focal mechanisms. According to our data and the deduced results, in northeast Iran, the Arabia–Eurasia convergence is taken up by strike-slip faulting along NE trending left-lateral and NNW trending right-lateral faults, as well as reverse to oblique-slip reverse faulting along NW trending faults. Such a structural assemblage is involved in a mechanically compatible and homogeneous modern stress field. This implies that no strain and/or stress partitioning or systematic block rotations have occurred in the Kopeh Dagh and Allah Dagh-Binalud deformation domains. The Plio-Quaternary stress changes documented in this paper call into question the extrapolation of the present-day seismic and GPS-derived deformation rates over geological time intervals encompassing tens of millions of years.     

Spatial mapping of earthquake hazard parameters in the Hindukush–Pamir Himalaya and adjacent regions: Implication for future seismic hazard

The study deals spatial mapping of earthquake hazard parameters like annual and 100-years mode along with their 90% probability of not being exceeded (NBE) in the Hindukush–Pamir Himalaya and adjoining regions. For this purpose, we applied a straightforward and most robust method known as Gumbel’s third asymptotic distribution of extreme values (GIII). A homogeneous and complete earthquake catalogue during the period 1900–2010 with magnitude M_W ⩾ 4.0 is utilized to estimate these earthquake hazard parameters. An equal grid point mesh, of 1° longitude X 1° latitude, is chosen to produce detailed earthquake hazard maps. This performance allows analysis of the localized seismicity parameters and representation of their regional variations as contour maps. The estimated result of annual mode with 90% probability of NBE is expected to exceed the values of M_W 6.0 in the Sulaiman–Kirthar ranges of Pakistan and northwestern part of the Nepal and surroundings in the examined region. The 100-years mode with 90% probability of NBE is expected to exceed the value of M_W 8.0 in the Hindukush–Pamir Himalaya with Caucasus mountain belt, the Sulaiman–Kirthar ranges of Pakistan, northwestern part of the Nepal and surroundings, the Kangra–Himanchal Pradesh and Kashmir of India. The estimated high values of earthquake hazard parameters are mostly correlated with the main tectonic regimes of the examined region. The spatial variations of earthquake hazard parameters reveal that the examined region exhibits more complexity and has high crustal heterogeneity. The spatial maps provide a brief atlas of the earthquake hazard in the region.     

The 10 June 2012 Fethiye Mw 6.0 aftershock sequence and its relation to the 24–25 April 1957 Ms 6.9–7.1 earthquakes in SW Anatolia,Turkey

The 10 June 2012 M_w 6.0 aftershock sequence in southwestern Anatolia is examined. Centroid moment tensors for 23 earthquakes with moment magnitudes (M_w) between 3.7 and 6.0 are determined by applying a waveform inversion method. The mainshock is a shallow focus strike-slip with reverse component event at a depth of 30 km. The seismic moment (M_o) of the mainshock is estimated as 1.28 × 10¹⁸ Nm and rupture duration of the Fethiye mainshock is 38 s. The focal mechanisms of the aftershocks are mainly strike-slip faulting with a reverse component. The geometry of the focal mechanisms reveals a strike-slip faulting regime with NE–SW trending direction of T-axis in the entire activated region. A stress tensor inversion of focal mechanism data is performed to obtain a more accurate picture of the Fethiye earthquake stress field. The stress tensor inversion results indicate a predominant strike-slip stress regime with a NW–SE oriented maximum horizontal compressive stress (S_H). According to variance of the stress tensor inversion, to first order, the Fethiye earthquake area is characterized by a homogeneous interplate stress field. The Coulomb stress change associated with the mainshock and the largest aftershock are also investigated to evaluate any significant enhancement of stresses along the Gulf of Fethiye and surrounding region. Positive lobes with stress more than 0.4 bars are obtained, indicating that these values are large enough to increase the Coulomb stress failure towards NNW–SSE and E–W directions.     

Geochemical signatures of uranium oxides in the Lufilian belt: From unconformity-related to syn-metamorphic uranium deposits during the Pan-African orogenic cycle

The Pan-African Lufilian belt (Zambia and Democratic Republic of Congo) is known for its world-class copper and cobalt deposits. In addition, the Lufilian Copperbelt hosts several uranium occurrences concentrated within deformed siliciclastic rocks of the basal Neoproterozoic Katanga Supergroup. We report LA-ICPMS and EMP analyses of the rare earth element (REE) and yttrium (Y) abundances (designated as the REY signatures) of uranium oxides from two uranium mineralizing events of the Lufilian belt previously dated at 652 ± 8 Ma and 530 ± 6 Ma by the U–Pb method on uraninite. Uranium oxides dated at ca. 650 Ma from the External fold-and-thrust belt are characterized by (i) bell shape REE patterns centered on middle REE (MREE), (ii) positive europium (Eu) anomalies and (iii) relatively low Y contents. In contrast, uranium oxides dated at ca. 530 Ma from the Domes region are characterized by (i) REE patterns but with a less pronounced light REE (LREE) fractionation, (ii) negative Eu anomalies and (iii) higher Y contents. Moreover, the External fold-and-thrust belt also contains uranium mineralization dated at ca. 530 Ma having the same characteristics as the ca. 530 Ma uranium oxides from the Domes region (a moderately fractionated REE pattern and a negative Eu anomaly).As REY signatures are known to reflect mineralizing processes, the distinct geochemical signatures of the two uranium oxide generations (ca. 650 Ma and ca. 530 Ma) provide meaningful information about the uranium cycle during the Pan-African orogeny. Compared to the REY signatures of the known worldwide uranium deposit types, the REY signature of uranium oxides dated at ca. 650 Ma of the External fold-and-thrust belt is similar to the REE patterns from unconformity-related U deposits (Athabasca in Canada and Kombolgie in Australia). Uranium oxides of the Domes region and some of the External fold-and-thrust belt display similar characteristics to syn-metamorphic U deposit (Mistamisk in Canada). Accordingly, we propose that the two stages of uranium oxide crystallizations within the Lufilian belt, at ca. 650 and ca. 530 Ma, occurred under distinct physico-chemical conditions. The first stage, at ca. 650 Ma, may be related to late diagenesis hydrothermal processes, at the basement/cover interface, with the circulation of highly saline basinal brines linked to evaporites of the Roan Group. This Pan-African unconformity-related uranium deposit is the youngest of this type described to date. The second stage may be connected to metamorphic fluid circulations, at about 530 Ma, during the Lufilian orogeny in the Domes region and also in the External fold-and-thrust belt.     

Miocene post-collisional shoshonites and their crustal xenoliths,Yarlung Zangbo Suture Zone southern Tibet: Geodynamic implications

The convergence between the Indian plate and the southern margin of the Eurasian continent created an active continental margin from Late Jurassic until about 40 Ma ago, which then evolved to form the Himalaya and the Tibetan Plateau during the continental collision stage. Post-collisional magmatism in southern Tibet, north of the Yarlung Zangbo Suture Zone (YZSZ) has been active since 45 Ma and is related to normal faulting and extensional tectonism. To date no such magmatism was reported within the YZSZ itself. This paper reports on the discovery of Miocene shoshonites within the YZSZ. They are significant because the magma traveled, at least in part, through oceanic crust, thus limiting interaction with the continental crust to the mid-crustal level and which affected the post-collisional magmatic rocks occurring in the northern part of the subduction system. In addition, xenoliths and xenocrysts of crustal origin in these rocks constrain the nature of metamorphic rocks underlying the YZSZ at mid-crustal level. The geochemical signatures of the shoshonitic rocks, including Nd and Sr isotope systematics, indicate derivation from a garnet-bearing middle continental crustal source. Crustal imprint complicates modeling of the petrogenetic processes which occurred prior to mid-crustal ponding of the magma which took place between 11 and 17 Ma at depths of 40 to 50 km. The significant role of crustal contamination raises serious concerns about models proposed for similar magmatic activity elsewhere in the Himalaya and the Tibetan Plateau.     

Mapping crustal structure beneath southern Tibet: Seismic evidence for continental crustal underthrusting

Receiver function imaging along a temporary seismic array (ANTILOPE-2) reveals detailed information of the underthrusting of the Indian crust in southern Tibet. The Moho dips northward from ~ 50 km to 80 km beneath the Himalaya terrane, and locally reaches ~ 85 km beneath the Indus–Yalung suture. It remains at ~ 80 km depth across the Lhasa terrane, and shallows to ~ 70 km depth under the Qiangtang terrane. An intra-crustal interface at ~ 60 km beneath the Lhasa terrane can be clearly followed southward through the Main Himalaya Thrust and connects the Main Boundary Thrust at the surface, which represents the border of the Indian crust that is underthrusting until south of the Bangong–Nujiang Suture. A mid-crustal low velocity zone is observed at depths of 14–30 km beneath the Lhasa and Himalaya terranes probably formed by partial melt and/or aqueous fluids.     

Isotopic characteristics of river sediments on the Tibetan Plateau

We systematically collected 40 modern clastic sediment samples from rivers in different tectonic units of the Tibetan Plateau and measured their Sr–Nd isotopic compositions. The isotopic characteristics provide insight into the controversial paleo-tectonic affinity of terranes of the Tibetan Plateau and the provenance of Songpan–Ganzi flysch complex. The Qilian Terrane and Himalaya Terrane have more negative ε_Nd(0) values (from ? 14.3 to ? 11.8 and from ? 20.64 to ? 13.26, respectively) and high ⁸⁷Sr/⁸⁶Sr values (from 0.719674 to 0.738818 and from 0.721020 to 0.824959, respectively), reflecting old and mature continental crust origin of these two terranes. The southern Lhasa Terrane is more radiogenic in ε_Nd(0) values (from ? 8.82 to ? 3.8) and low in ⁸⁷Sr/⁸⁶Sr values (from 0.711504 to 0.719489), implying the combined impact of the Neo-Tethys mantle and Himalaya old continental crust. Sr–Nd isotopic compositions of the Qilian Terrane are similar to those in the Yangtze Craton, indicating that the Qilian Terrane was probably separated from the Yangtze Craton. Sr–Nd isotopic characteristics of the Songpan–Ganzi Terrane are similar to the Yangtze Craton and are remarkably different to those in the North China Craton, eastern Kunlun–Qaidam and the central Qiangtang metamorphic belt, implying that the widely distributed flysch complex of the Songpan–Ganzi Terrane was sourced from the Yangtze Craton.     

New field and seismic data about the intraplate strike-slip deformation in Van region,East Anatolian plateau,E. Turkey

The study area is the Van earthquake region. It is located in the western section of the East Anatolian–Iranian plateau outside and to the east of the Karlıova triple junction. Based on the tectonic periods, the rock units exposed in the study area are classified into two common categories. These are the Pre-Late Pliocene paleotectonic units and the Plio-Quaternary neotectonic units. The Paleotectonic units are composed of the Yüksekova Complex of Campanian–Maastrichtian age and the Kırkgeçit Formation of Oligo-Miocene age. The paleotectonic units are intensely deformed (folded, thrust to reverse faulted and converted into an imbricate stack). The neotectonic units are composed of fluvio-lacustrine sedimentary facies with volcanic interclations. It is full of soft-sedimentary structures such as deltaic structure, slump fold, sand dikes to sills and normal to reverse types of growth faults which imply to a sedimentation accompanied by both a volcanic activity and active tectonics. Originally the Paleotectonic units are overlain with an angular unconformity by the nearly flat-lying neotectonic units. This angular unconformity and the big difference in the deformational patterns of both categories of rock units indicate an inversion in tectonic regime in Late Pliocene. The new tectonic regime is the strike-slip faulting-dominated neotectonic regime. It is governed by an approximately N–S-directed compression, and composed of NW- to NE-trending strike-slip faults, N–S trending oblique-slip normal faults to fissures and the E–W trending thrust to reverse faults. Most of thrust to reverse faults are inherited from the Pre-Late Pliocene paleotectonic regime. Some of them have reactivated and led to the occurrence of large and devastative earthquakes. The last devastative seismic event is the 23 October 2011 Tabanlı (Van) earthquake of Mw = 7.2 that caused 644 deaths and moderate to heavy damage of ¼ of structures (28,532) in Van earthquake region. The source of the Tabanlı earthquake is the Everek erosional reverse fault. In addition the Tabanlı earthquake is the largest seismic event occurred till now in Turkey. It was followed by a series (over 6000) of small-sized aftershocks and severeal moderate-sized indepentent earthquakes of reverse, normal and strike-slip faulting origin. Both the field and new seismic data strongly reveal that the prominent tectonic regime in the East Anatolian plateau is the strike-slip neotectonic regime, not the tensional tectonic regime as has been reported in some previous works. The strike-slip faulting and related deformation are confined into the upper shallowing part (up to 40 km) of the crust, whilst the extensional deformations are the subcrustal processes and being taking place in a squashy zone at the depths of approximately 40–60 km.     

Deglaciation and landscape history around Annapurna,Nepal, based on 10Be surface exposure dating

The High Himalaya is a key area for tectonic, geomorphological and climate studies, because of its extreme relief and location at the transition zone between areas with abundant monsoonal precipitation and the arid/semiarid Tibetan Plateau. We present ¹⁰Be surface exposure ages on 22 boulders from the Annapurna area in Nepal. The ages improve understanding of the Late Quaternary landscape history and the geomorphological processes operating in this part of the Himalaya.Although our study is reconnaissance in nature, it highlights the importance of catastrophic events, such as landslides and debris flows, for denudation of high mountains. Holocene exposure ages for the Dhampu–Chooya landslide (～4.1 ± 0.6 ka) and for 600 m of alluviation in Kali Gandaki Valley (～2.1 ± 0.6 ka), for example, indicate the frequent occurrence and extent of catastrophic events and their implications for natural hazards. We also offer an explanation for the differences in Late Quaternary glacial chronologies at closely spaced study sites in the Nepal Himalaya. Topographically controlled and spatially variable precipitation in the Himalaya determines the sensitivity of glaciers to changes in temperature and precipitation. Accordingly, some glaciers advanced in-phase with Northern Hemisphere ice sheets, whereas others reached their maximum extent at times of increased monsoonal precipitation during Marine Isotope Stage 3 and the early Holocene.     

Rapid cooling of the Yanshan Belt,northern China: Constraints from 40Ar/39Ar thermochronology and implications for cratonic lithospheric thinning

The Yanshan Orogenic Belt is located in the northern part of the North China Craton (NCC), which lost ∼120 km of lithospheric mantle during Phanerozoic tectonic reactivation. Mesozoic magmatism in the Yanshan fold-and-thrust belt began at 195–185 Ma (Early Jurassic), with most of the granitic plutons being Cretaceous in age (138–113 Ma). Along with this magmatism, multi-phase deformational structures, including multiple generations of folds, thrust and reverse faults, extensional faults, and strike-slip faults are present in this belt. Previous investigations have mostly focused on geochemical and isotopic studies of these magmatic rocks, but not on the thermal history of the Mesozoic plutons. We have applied ⁴⁰Ar/³⁹Ar thermochronology to biotites and K-feldspars from several Lower Cretaceous granitic plutons to decipher the cooling and uplift history of the Yanshan region. The biotite ⁴⁰Ar/³⁹Ar ages of these plutons range from 107 to 123 Ma, indicating that they cooled through about 350 °C at that time. All the K-feldspar step-heating results modeled using multiple diffusion domain theory yield similarly rapid cooling trends, although beginning at different times. Two rapid cooling phases have been identified at ca. 120–105 and 100–90 Ma. The first phase of rapid cooling occurred synchronously with widespread extensional deformation characterized by the formation of metamorphic core complexes, A-type magmatism, large-scale normal faults, and the development of half-graben basins. This suggests rapid exhumation took place in an extensional regime and was a shallow-crustal-level response to lithospheric thinning of the NCC. The second phase of rapid cooling was probably related to the regional uplift and unroofing of the Yanshan Belt, which is consistent with the lack of Upper Cretaceous sediments in most of the Yanshan region.     

The Isle of Wedmore relay ramp: how fault evolution created King Alfred's historic landmark

The Isle of Wedmore covers an area of ~ 19 km², rises up to ~ 65 m above the surrounding lowlands of the Somerset Levels, and was an island until the Middle Ages. The topography is interpreted as having been formed by a relay ramp between two right-stepping faults (the Weare Fault to the west and the Mudgley Fault to the east) which have tens of metres of downthrow to the south, and which are probably normal faults. The relay ramp has a dip of about 3° to the SW and is breached by the NW-striking Wedmore Fault, which has up to ~ 23 m downthrow to the NE. Several NE-trending faults occur in the relay ramp, which are interpreted as having formed when the relay ramp became a contractional step when the Weare and Mudgley faults underwent sinistral reactivation, or as N–S contraction occurred during the Cenozoic. Analogues for this behaviour are presented from the Liassic rocks on the coast between Lilstock and East Quantoxhead.     

Depth dependency of stress ratios during the sedimentation of NW Gyeongsang Basin (Cretaceous), southeast Korea: Estimate of stress parameters and timing of tectonic episodes

The natural stress state in the lithosphere consists of the vertical load and Poisson ratio, and then additional horizontal compression and extension (denoted by Δσ_H and Δσ_T, respectively) are assumed to be superimposed upon this gravitational stress field. The resulting stress state is composed of the maximum, medium and minimum stresses denoted by σ₁, σ₂, and σ₃, respectively. The stress ratio is given as Φ = (σ₂ − σ₃)/(σ₁ − σ₃). A linear relation is found between Φ or 1/Φ and the vertical load in wrench-faulting and extensional stress regimes, respectively. The slope and intercept of the linear relation result in the additional horizontal stresses and level of (average) paleo-surface, respectively. Stress ratio is also determinable by the stress tensor inversion of fault-slip data. The earliest tectonic event (T_1 Event) in the Cretaceous Gyeongsang Basin consists of coexisting E–W compression and N–S extensional faulting episodes. Plots of Φ or 1/Φ against the burial depth (or vertical load) display several linear trends: two clusters in extensional episodes, and two or three clusters in wrench-faulting episodes. Because Δσ_H is assumed to be null or negligible in the extensional regime, Δσ_T is determinable from the slopes of two linear clusters as being −2.5 to −4.0 MPa. In wrench-faulting episodes, the values of Δσ_H are given to be 61.6–101.4 MPa by applying determined additional horizontal extensions. Determined levels of average paleo-surfaces and those of syndepostional structures illustrate that more than five wrench-faulting or extensional episodes have occurred during the T_1 Event, whose active age, consequently, ranges from the Barremian to the Coniacian. This supports that the coexisting coaxial faulting episodes with the same extension may correspond to the alternation of wrench-faulting and extensional episodes.     

Paleoseismological evidence for non-characteristic behavior of surface rupture associated with the 2004 Mid-Niigata Prefecture earthquake, central Japan

The 2004 Mid-Niigata Prefecture earthquake sequence (mainshock magnitude, M_JMA 6.8), which occurred in an active fold-and-thrust belt in northern central Japan, generated a small thrust surface rupture (< 20 cm of vertical displacement) along a previously unmapped northern extension of the active Muikamachi–Bonchi–Seien fault zone, on the eastern margin of the epicentral region. To better understand past seismic behavior of the rupture, we conducted a paleoseismic trenching study across the 10-cm-high west-side-up surface rupture at the foot of a pre-existing 1.8-m-high east-facing scarp, which probably resulted from past earthquake(s). A well-defined west-dipping thrust fault zone accompanied by drag folding and displacing the upper Pliocene to lower Pleistocene strata and the unconformably overlying upper Pleistocene (?) to Holocene strata was exposed. The principal fault zone is connected directly to the 2004 surface rupture. From the deformational characteristics of the strata and radiocarbon dating, we inferred that two large paleoseismic events occurred during the past 9000 years prior to the 2004 event. These two pre-2004 events have a nearly identical fault slip (at minimum, 1.5 m), which is ≥ 15 times that of the 2004 event (∼ 10 cm). These paleoseismic data, coupled with the geological and geomorphological features, suggest that the 2004 event represented non-characteristic behavior of the fault, which can potentially generate a more destructive earthquake accompanied by meter-scale surface displacement. This study provides insight into the interpretation of past faulting events and increases our understanding of rupture behavior.     

Influence of the 2008 Wenchuan earthquake (Mw 7.9) on the occurrence probability of future earthquakes on neighboring faults

The Longquan–Shan fault and the Huya fault are two major neighboring faults of the Longmen–Shan fault zone where the 12 May 2008 Wenchuan earthquake (M_w 7.9) occurred. To study the influence of the Wenchuan event on these two active faults, we calculate changes of Coulomb stress on the Longquan–Shan fault and the Huya fault caused by the Wenchuan mainshock. Our results indicate that the Coulomb stress in the northern section (Zone A) of the Longquan–Shan fault is increased by 0.07–0.10 bars, that in the middle section (Zone B) by 0.04–0.11 bars, and that in the southern section (Zone C) shows almost no change. For the Huya fault, the Coulomb stress is decreased by 0.01–0.03 bars in the northern section (Zone A), 0.10–0.35 bars in the middle section (Zone B), and nearly 0.5 bars in the southern section (Zone C). The epicenter distribution of small earthquakes (M_L ⩾ 1.5) on the Longquan–Shan fault and the Huya fault after the Wenchuan earthquake is consistent with the distribution of the Coulomb stress change. This implies that the Wenchuan earthquake may have triggered small events on the Longquan–Shan fault, but inhibited those on the Huya fault. We then use the rate/state friction law to calculate the occurrence probability of future earthquakes in the study region for the next decade. They include the distribution of b-values, magnitude of completeness (M_c), the background seismicity rate, a value of Aσ_n and the duration for the transient effect (t_a) in the study region. We also estimate the earthquake occurrence probabilities on the neighboring faults after the Wenchuan earthquake. Our results show that, the occurrence probability of future earthquakes in the Longquan–Shan has a slight increase, being 7% for M ⩾ 5.0 shocks during the next decade, but the earthquake probability in the Huya region is reduced obviously, being 5–20%, 7–26% and 3–9% for M ⩾ 5.0 shocks during the next decade in sections A, B and C of the Huya fault, respectively.     

Differential timing of vertical-axis block rotations in the northern Ryukyu Arc: Paleomagnetic evidence from the Koshikijima Islands,Japan

Over 300 samples for paleomagnetic analysis and K–Ar dating were collected from 27 sites at NW–SE and NE–SW trending dike swarms (herein, NW dikes and NE dikes, respectively) in the Koshikijima Islands, northern Ryukyu Arc. The NW dikes are Middle Miocene in age and have directions (D = ? 37.7^°, I = 51.8^°, α₉₅ = 9.6^°, and κ = 40.8) that are deflected westward relative to the stable eastern Asian continent. Conversely, the NE dikes, of Late Miocene age, have directions (D = 16.1^°, I = 57.7^°, α₉₅ = 7.1^°, and κ = 41.9) that show no such deflection. These differences are interpreted as indicating that the Koshikijima Islands underwent approximately 40^° of counter-clockwise rotation during the Middle to Late Miocene. A synthesis of the paleomagnetic and structural data suggests a three-stage history of extensional deformation: (1) displacement upon normal faults (F₁ faults) without vertical-axis block rotation, (2) strike-slip reactivation of F₁ faults and oblique-normal displacement on NE–SW-trending faults (F₂ faults) with vertical-axis block rotation, and (3) oblique-normal displacement on F₂ faults without vertical-axis block rotation. Regional differences in the timing and amount of counter-clockwise vertical-axis block rotations indicate that the northern Ryukyu Arc rotated as several distinct rigid blocks.     

Quantitative analysis of the shallow crustal tectonic stress field in China mainland based on in situ stress data

The latest hydraulic fracturing and stress relief measurement data in the Chinese mainland were collected. The total of 3856 data entries are measured at 1474 locations. The measured area covers 75–130°E and 18–47°N, and the depth range varies from surface to 4000 meters depth, which generally includes each active tectonic block of China and each segment of North–South seismic belt. We investigated the tectonic stress field by removing the effect of gravity. For this, we assume lateral constraints and Heim’s rule. The gravity contribution is removed by using the assumption of lateral constraint and Heim’s rule. Our results show: (1) the maximum and the minimum horizontal principal stress σ_H, σ_h and the vertical stress σ_V in the shallow crust of China all increase linearly with depth: σ_H = 0.0229D + 4.738, σ_h = 0.0171D + 1.829, σ_V = 0.0272D. Maximum and minimum horizontal tectonic stress varies as a function of depth D linearly 4.738 < σ_T < 0.0139D + 4.738 and 1.829 < σ_t < 0.0162D + 1.829. The horizontal tectonic differential stress is σ_T − σ_t = 0.0058D + 2.912. (2) The intermediate value of σ_T1 (regression value of tectonic stress inferred from the assumption of lateral constraint at 2000 m depth) changes in different areas, the maximum value of which is 45.6 MPa, while the minimum value of which is 26.8 MPa. Horizontal tectonic differential stress σ_T − σ_t increases linearly with depth and the maximum and minimum of σ_T − σ_t is 25.3 MPa and 13.0 MPa, respectively. In general, the stress magnitude is much higher in western than in eastern China. This indicates that the strong Indo-Eurasian collision dominates the present tectonic stress field in Chinese mainland. (3) Compared with other study regions, the northward crustal compression to the Qinghai-Tibet block is relatively lower in magnitude in the shallow subsurface and higher at deeper depth. (4) The orientations of σ_T in China mainland generally form a radial scattering pattern centered in Tibetan Plateau. From western to eastern China, they rotate gradually clockwise from NS to NNE, NE, NEE, and SE, which is consistent with the result of focal mechanism solutions.     

Metallogeny,structural evolution,post-mineral cover distribution and exploration in concealed areas of the northern Chilean Andes

Mineral exploration of prospective areas concealed by extensive post-mineralization cover is growing, being very complex and expensive. The projection of rich and giant Paleocene to early Oligocene porphyry-Cu-Mo belts in northernmost Chilean Andes (17.5–19.5°S) has major exploration potential, but only a few minor deposits have been reported to date, due to the fact that the area is largely covered by post-mineral strata. We integrate the Cenozoic stratigraphic, structural and metallogenic evolution of this sector, in order to identify the most promising regions related to lesser post-mineral cover and the projection of different metallogenic belts. The Paleocene to early Eocene metallogenic belt extends along the Precordillera, with ca. 30 km wide, and includes porphyry-Cu prospects and small Cu (±Mo-Au-Ag) vein and breccia-pipe deposits. Geochronological data indicate an age of 55.5 Ma for an intrusion related to one deposit and ages from 69.5 to 54.5 Ma for hydrothermal alteration in one porphyry-Cu prospect and largest known Cu deposits. The middle Eocene to early Oligocene porphyry belt, in the Western Cordillera farther east, is associated with 46–44 Ma intrusions. It is estimated to be 40-km wide, but is largely concealed by thick post-mineral cover. The youngest Miocene to early Pliocene metallogenic belt, also in the Western Cordillera, is well-exposed and includes Au-Ag epithermal and polymetallic veins and manto-type deposits.The Oligocene-Holocene cover consists of a succession of continental sedimentary and volcanic rocks that overall increase in thickness from 0 to 5000 m, from west to east. These strata are subhorizontal in the west and folded-faulted towards the east. Miocene gentle anticlines and monocline flexures extend along strike for 30–60 km in the Precordillera and were generated by propagation of high-angle east-dipping blind reverse faults with at least 300–900 m of Oligocene bedrock offset. The thickness of cover exceeds 2000 m in the eastern Central Depression, whereas it is generally less than 1000 m in the Precordillera along the Paleocene to early Eocene porphyry-Cu belt and it can reach locally up to 5000 m in the Western Cordillera, above the middle Eocene to early Oligocene belt.In the studied Andean segment, the Miocene to early Pliocene metallogenic belt is superimposed on the Paleocene to Oligocene belts in a 40–50 km wide zone. This overlap may be explained by an accentuated migration of the magmatic front, from east to west, since ca. 25 Ma, as a consequence of subduction slab steepening after a period of magmatic lull and flat subduction from ca. 30–35 to 25 Ma. The identified areas of lesser cover thickness are prone to exploration for concealed deposits, especially along the projection of major porphyry-Cu-Mo belts.     

Detrital zircon U–Pb geochronology,trace-element and Hf isotope geochemistry of the metasedimentary rocks in the Eastern Himalayan syntaxis: Tectonic and paleogeographic implications

The origin of the Greater Himalayan Sequence in the Himalaya and the paleogeographic position of the Lhasa terrane within Gondwanaland remain controversial. In the Eastern Himalayan syntaxis, the basement complexes of the northeastern Indian plate (Namche Barwa Complex) and the South Lhasa terrane (Nyingchi Complex) can be studied to explore these issues. Detrital zircons from the metasedimentary rocks in the Namche Barwa Complex and Nyingchi Complex yield similar U–Pb age spectra, with major age populations of 1.00–1.20 Ga, 1.30–1.45 Ga, 1.50–1.65 Ga and 1.70–1.80 Ga. The maximum depositional ages for their sedimentary protoliths are ~ 1.0 Ga based on the mean ages of the youngest three detrital zircons. Their minimum depositional ages are ~ 477 Ma for the Namche Barwa Complex and ~ 499 Ma for the Nyingchi Complex. Detrital zircons from the Namche Barwa Complex and Nyingchi Complex also display similar trace-element signatures and Hf isotopic composition, indicating that they were derived from common provenance. The trace-element signatures of 1.30–1.45 Ga detrital zircons indicate that the 1.3–1.5 Ga alkalic and mafic rocks belt in the southeastern India is a potential provenance. Most 1.50–1.65 Ga zircons have positive ε_Hf(t) values (+ 1.2 to + 9.0), and most 1.70–1.80 Ga zircons have negative ε_Hf(t) values (− 7.1 to − 1.9), which are compatible with those of the Paleo- to Mesoproterozoic orthogneisses in the Namche Barwa Complex. Provenance analysis indicates that the southern Indian Shield, South Lhasa terrane and probably Eastern Antarctica were the potential detrital sources. Combined with previous studies, our results suggest that: (1) the Namche Barwa Complex is the northeastern extension of the Greater Himalaya Sequence; (2) the metasedimentary rocks in the Namche Barwa Complex represent distal deposits of the northern Indian margin relative to the Lesser Himalaya; (3) the South Lhasa terrane was tectonically linked to northern India before the Cambrian.