Significant crustal structural variation across the Chaochou Fault, southern Taiwan: New tectonic implications for convergent plate boundary

The Chaochou Fault, a major geological boundary in southern Taiwan is considered to be a part of the convergent plate boundary between the Eurasia Plate and the Philippine Sea Plate. We applied the Common Conversion Point stacking technique to teleseismic radial receiver functions and obtained Moho variation and crustal structure across the Chaochou Fault. In the Eurasia Plate to its west, the Moho depth is about 37 km and the crust is subducting to the east beneath the Philippine Sea Plate with a dip angle of about 30° between the Backbone Belt and the Tananao Schist. In the Philippine Sea Plate, the Moho depth is about 17 km. The Longitudinal Valley marks the collision boundary between the Eurasia Plate and the Philippine Sea Plate. The results suggest that the depth extent of the Chaochou Fault is about 30–35 km and the fault becomes a “shallow-angle” thrust fault at depth. The Common Conversion Point image also shows several bending interfaces of velocity contrast in the crust. We proposed a simple model to explain the Philippine Sea Plate and Eurasia Plate collision process and the observed crustal deformations.     

Displacement and timing of left-lateral faulting in the Kunlun Fault Zone,northern Tibet,inferred from geologic and geomorphic features

The E-W to WNW-ESE striking Kunlun Fault Zone, extending about 1600 km, is one of the large strike-slip faults in the northern Tibet, China. As a major strike-slip fault, it plays an important role on the extrusion of Tibet Plateau in accommodating northeastward shortening caused by the India-Asia convergence. However, the time of initiation left-lateral faulting of the Kunlun Fault Zone is still largely debated, ranging from the Middle to Late Triassic (240–200 Ma) to early Quaternary (2 Ma). We document displaced basement rocks and geomorphic features along the Kunlun Fault Zone, based on tectono-geomorphic interpretation of satellite remote sensing images and field geologic and geomorphic observations. Our results show that the largest cumulative offset of basement rocks is likely to be 100 ± 20 km. Meanwhile, a series of pull-apart basins (Kusai, Xiugou and Tuosu lake basins) and pressure ridges (East Deshuiwai and Maji Snow Mountains), each 45–70 km long and ∼8–12 km wide, are developed along the Kunlun Fault Zone, which resulted from long-term tectono-geomorphic growth since the Late Miocene or Early Pliocene. Geologic evidence indicates that the Kunlun Fault Zone had a long-term slip rate of ca.10 mm/yr during the late Quaternary. This slip rate is similar to that shown by present-day GPS measurements. Thus, we estimate that the Kunlun Fault Zone probably began left-lateral faulting at 10 ± 2 Ma based on a total displacement of 100 ± 20 km, and assuming a constant long-term slip rate of ca.10 mm/yr for several millions of years. And this timing constraint on initiation of left-lateral faulting of the Kunlun Fault Zone is consistent with widespread tectonic deformation which occurred in the Tibetan Plateau.     

Variable behavior of the Dead Sea Fault along the southern Arava segment from GPS measurements

The Dead Sea Fault is a major strike-slip fault bounding the Arabia plate and the Sinai subplate. On the basis of three GPS campaign measurements, 12 years apart, at 19 sites distributed in Israel and Jordan, complemented by Israeli permanent stations, we compute the present-day deformation across the Wadi Arava fault, the southern segment of the Dead Sea Fault. Elastic locked-fault modelling of fault-parallel velocities provides a slip rate of 4.7 ± 0.7 mm/yr and a locking depth of 11.6 ± 5.3 km in its central part. Along its northern part, south of the Dead Sea, the simple model proposed for the central profile does not fit the velocity field well. To fit the data, two faults have to be taken into account, on both sides of the sedimentary basin of the Dead Sea, each fault accommodating ∼ 2 mm/yr. Locking depths are small (less than 2 km on the western branch, ∼ 6 km on the eastern branch). Along the southern profile, we are once again unable to fit the data using the simple model, similar to the central profile. It is very difficult to propose a velocity greater than 4 mm/yr, i.e. smaller than that along the central profile. This leads us to propose that a part of the relative movement from Sinai to Arabia is accommodated along faults located west of our profiles.     

Geodetic and seismological investigation of crustal deformation near Izmir (Western Anatolia)

The Aegean region including western Turkey, mainland Greece, and the Hellenic Arc is the most seismological and geodynamical active domain in the Alpine Himalayan Belt. In this study, we processed 3 years of survey-mode GPS data and present the analysis of a combination of geodetic and seismological data around Izmir, which is the third most populated city in Turkey. The velocities obtained from 15 sites vary between 25 mm/yr and 28 mm/yr relative to the Eurasian plate. The power law exponent of earthquake size distribution (b-value) ranges from 0.8 to 2.8 in the Izmir region between 26.2°E and 27.2°E. The lowest b-value zones are found along Karaburun Fault (b = 0.8) and, between Seferihisar and Tuzla Faults (b = 0.8). A localized stress concentration is expected from numerical models of seismicity along geometrical locked fault patches. Therefore, areas with lowest b-values are considered to be the most likely location for a strong earthquake, a prediction that is confirmed by the 2005 M_w = 5.9 Seferihisar earthquake sequences, with epicentres located to the south of the Karaburun Fault. The north–south extension of the Izmir area is corroborated by extension rates up to 140 nanostrain/yr as obtained from our GPS data. We combined the 3-year GPS velocity field with the published velocity field to determine the strain rate pattern in the area. The spatial distribution of b-value reflects the normal background due to the tectonic framework and is corroborated by the geodetic data. b-Values correlate with strain pattern. This relationship suggests that decrease of b-values signifies accumulating strain.     

Thermal evolution of the Hengshan extensional dome in central South China and its tectonic implications: New insights into low-angle detachment formation

The Hengshan massif is an exhumed, mid-crustal, plutonic–metamorphic dome formed during Cretaceous crustal extension in the Jiangnan orogenic belt, central South China. Multiple thermochronometers (mica ⁴⁰Ar/³⁹Ar, apatite fission track and zircon (U–Th)/He) are applied to its footwall along a slip-parallel transect to quantify its thermal history and cooling rate, and the slip magnitude, rate, initial geometry and kinematic evolution of the low-angle Hengshan detachment fault. Our thermochronological data, in conjunction with previous ages, indicate that (1) footwall rocks cooled from ~ 700 °C to ~ 60 °C in less than 60 Myr (136–80 Ma) at variable rates ranging from ~ 50 °C/Myr to ~ 13 °C/Myr, (2) the Hengshan detachment fault accommodated ~ 8–12 km of total slip at variable slip rates from 0.14 to 1 mm/yr during tectonic exhumation, (3) the footwall has been tilted ~ 26°–50° to the east since slip began, indicating that the low-angle Hengshan detachment fault initiated at a steep dip and was passively rotated to a more gentle orientation during subsequent normal slip. This study provides compelling evidence supporting that the low-angle detachment fault in the extensional dome can be generated by the reactivation and passive rotation of an initially steep reverse fault during normal slip. In addition, our thermochronological data constrain the time of extension in the Hengshan dome between 136 and 80 Ma, which implies that the back-arc extension within South China associated with the rollback of the Paleo-Pacific slab might have lasted until at least 80 Ma.     

Source complexity of the 4 March 2010 Jiashian,Taiwan, Earthquake determined by joint inversion of teleseismic and near field data

The Jiashian earthquake (M_L 6.4) occurred on 4 March 2010. It was the largest inland event in southern Taiwan of 2010. The mainshock location was unexpected since it occurred in an area with relatively low background seismicity. In addition, reports of earthquake focal mechanisms do not fit with any known active fault geometry. In order to understand the origin of this earthquake, especially its rupture process, we perform a joint source inversion by using teleseismic body wave, GPS coseismic displacements and near field ground motion data. In this study, we considered a northwest–southeast trending fault with a northeast dip retrieved from GPS coseismic data and aftershocks distribution. To analyze the detailed slip distribution in space and time, we used near field 3D Green’s functions provided by spectral-element method and a full time–space inversion technique. We find a complex rupture process with several slip patches distributed inside two main asperities. The slip map reveals a mean slip of 12.9 cm for a maximum slip of 27.3 cm leading to a M_w 6.47 for this event. The rupture initiates in the deepest portion of the fault at 20 km depth, and propagated upward up to 2 km depth to form the two asperities. The source time function of this event revealed two pulses corresponding to the two asperities, for a total duration time of about 16 s. Most aftershocks occurred near the upper boundary of the deepest asperity while no aftershocks are located close to the shallowest one. We infer that the locations of these slip patches are related to the surrounding fault systems that may have restricted the rupture propagation during the earthquake.     

Neotectonics of the Marikina Valley fault system (MVFS) and tectonic framework of structures in northern and central Luzon, Philippines

Recognition of neotectonic features along the Marikina Valley fault system (MVFS) in central Luzon, Philippines indicates a dominantly dextral strike-slip motion during its most recent activity believed to be Late Pleistocene to Holocene in age. Variations in the ratios of vertical to horizontal displacements for the segments imply a dominantly dextral motion of the West Marikina Valley fault (WMVF) and oblique dextral motion for the East Marikina Valley fault (EMVF). The displacement data further suggest that rupturing along the EMVF involved multiple segments and occurred separately from the events along the WMVF segments. Estimated earthquake magnitudes for the WMVF and EMVF based on single-event offsets fall within the range M 7.3–7.7. The vertical slip component in the northern part of the Marikina Valley is associated with the development of a basin between the EMVF and WMVF while the large vertical component in the southernmost segment of the EMVF (Talim) is attributed to volcanism-related extension. Lateral advection of the block bounded by the MVFS and the Philippine fault zone (PFZ), rather than pure shear resulting from an assumed east–west compression, best explains the observed kinematics of the MVFS. This is the result of compression during the westward drift of the Philippine Sea Plate and northern Luzon and occurs through slip along the WMVF and EMVF at rates of 5–7 mm/yr.     

Western Tibet relief evolution since the Oligo-Miocene

Western Tibet, between the Karakorum fault and the Gozha–Longmu Co fault system, is mostly internally drained and has a 1.5–2 km amplitude relief with km-large valleys. We investigate the origin of this peculiar morphology by combining a topography analysis and a study of the Cenozoic sedimentation in this area. Cenozoic continental strata correspond to a proximal, detrital fan deposition, and uncomformably rest on a palaeorelief similar to the modern one. Zircon U–Pb dating from trachytic flows interbedded within the Cenozoic continental sediments indicates that detrital sedimentation occurred at least between ca 24 and 20 Ma in the Shiquanhe basin, while K/Ar ages suggest it may have started since ~ 37 Ma in the Zapug basin. The distribution of continental deposits shows that present-day morphology features, including km-large, 1500 m-deep valleys, were already formed by Early Miocene times. We suggest that today's internally drained western Tibet was externally drained, at least during late Miocene, contemporaneously with early motion along the Karakorum Fault. Detailed study of the present day river network is compatible with a dextral offset on the Karakorum Fault of 250 km at a rate of ~ 10 ± 1 mm/yr. Displacement along the Karakorum fault possibly induced the shift from external to an internal drainage system, by damming of the Bangong Co ~ 4 Ma ago, leading to the isolation and preservation of the western Tibet relief.     

Field evidences from northern Dead Sea Fault Zone (South Turkey): New findings for the initiation age and slip rate

The left-lateral strike–slip Dead Sea Fault Zone (DSFZ) extends from the Red Sea in the south to the East Anatolian Fault Zone (EAFZ) in the north. This study examines the northern part of the DSFZ around Amik Basin and presents surface and subsurface geological evidence for the Quaternary activity and initiation age of the northernmost DSFZ. The DSFZ extends N–S in the south of the Amik Basin where clear geological and morphological evidence exists for faulting. Geological observations around Amik Basin, analyses of borehole data and electrical resistivity profiles within the Amik Basin indicate that the activity of the northern DSFZ started after Pliocene in the Amik Basin. Subsurface data in the basin suggest that the DSFZ offsets a pre-Quaternary basin sinistrally by about 7.9 km. The offset pre-Quaternary basin suggests at least 4.94 ± 0.13 mm/year slip rate for the northern part of the DSFZ. The Karasu Fault Zone (KFZ) extends in an en-echelon pattern along the western margin of the Karasu Valley and it transfers the significant amount of slip from DSFZ to the EAFZ.     

Paleomagnetic versus GPS determined tectonic rotation around eastern Himalayan syntaxis in East Asia

Northward indentation of the Indian Plate has brought about significant tectonic deformation into East Asia. A record of long-term tectonic deformation in this area for the past 50 M yr, particularly the vertical axis rotation, is available through paleomagnetic data. In order to depict rotational deformation in this area with respect to Eurasia, we compiled reliable paleomagnetic data sets from 79 localities distributed around eastern Himalayan syntaxis in East Asia. This record delineates that a zone affected by clockwise rotational deformation extends from the southern tip of the Chuan Dian Fragment to as far as the northwestern part of the Indochina Peninsula. A limited zone that experienced a significant amount of clockwise rotation after an initial India–Asia collision is now located at 23.5°N, 101°E, far away from an area (27.5°N, 95.5°E) where an intense rotational motion has been viewed by a snapshot of GPS measurements. This discrepancy in clockwise rotated positions is attributed to southeastward extrusion of the tectonic blocks within East Asia as a result of ongoing indentation of the Indian Plate. A quantitative comparison between the GPS and paleomagnetically determined clockwise rotation further suggests that following an initial India–Asia collision the crust at 30°N, 94°E paleoposition was subjected to southeastward displacement together with clockwise rotation, which eventually reached to present-day position of 23.5°N, 101°E, implying a crustal displacement of about 1000 km during the past 50 M yr.     

Mechanical relation between crustal rheology,effective fault friction,and strike-slip partitioning among the Xiaojiang fault system,southeastern Tibet

The north–south trending Xiaojiang fault system accommodates ～10–12 mm/yr sinistral motions between southeastern Tibet and south China. In the south segment, the fault system composes mainly of four parallel strike-slip faults, namely from west to east, the Luzhijiang fault, the Yimen fault, the Puduhe fault, and the Xiaojiang fault. Geological and Seismological observations have shown that these strike-slip faults are all of active, while the slip rates of the Luzhijiang, the Yimen, and the Puduhe faults are much less than that of the Xiaojiang fault. We use finite element modeling to explore the mechanical relation between crustal rheology, effective fault friction and long-term slip rate partitioning among the four parallel faults. The individual faults are simplified as vertical discontinuities embedded in the crust as geophysical explorations have predicted. A large number of models are tested, associating with variations of the crustal rheolohy and the effective fault friction of individual faults. Result shows that if crust bounding the faults trends to behave like rigid blocks and decoupled mechanically from underlying layer, the modeled result is hard to approximate slip rates of the individual faults. To better fit slip rates of the individual faults, viscous deformation of the lower crust should be included. With a heterogeneously viscous lower-crust model that is built upon thermal structure of the heat flow data, associating with relatively low effective friction of the Xiaojiang fault, the modeled results fit the geological slip rates well, with ～1–1.5 mm/yr for the Luzhijiang, the Yimen and the Puduhe faults, and ～6–6.5 mm/yr for the Xiaojiang fault. Thus, in the southward movement of the Tibetan plateau around the eastern Himalayan syntaxis, slip partitioning among the Xiaojiang fault system should be related to viscous deformation of the lower crust associated with different strength of the individual faults, highlighting that deformation of this fault system is coupled mechanically between the frictional upper crust and the viscous lower crust.     

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.     

Late Quaternary slip rates and slip partitioning on the Alpine Fault,New Zealand

Published geological data on Late Quaternary offsets on the Alpine Fault, New Zealand, have been assembled into a common format and analysed with respect to uncertainties. Uncertainties arise mainly from measurement of offset features, relating apparent offsets to actual fault slip, and dating the offset features. Despite the considerable uncertainties, the data form a coherent set consistent with a relatively constant rate of strike-slip of 27±5 mm/year between Milford Sound and Hokitika. This rate represents 70–75% of the fault-parallel interplate motion. North of the confluence with the Hope Fault, the rate drops substantially. Dip-slip rates, on the other hand, show considerable variation along strike, rising to a maximum of more than 10 mm/year in the central section and decreasing to zero at the southern end. Partitioning of c. 25% of the interplate slip on to structures east of the Alpine Fault occurs in the central section, consistent with predictions from critical wedge models. The partitioning of all the fault-normal component of displacement on to other structures in the south may be related, in part, to a doubling in width of the deforming wedge to the east. Most probably, however, the fault-normal displacement is mainly accommodated by underthrusting of the Australian plate offshore, due to a change in the nature of the crust from continental to oceanic.     

Late Jurassic intraplate faulting in eastern Australia: A link to subduction in eastern Gondwana and plate tectonic reorganisation

Eastern Gondwana was subjected to subduction processes during the Middle-Late Jurassic, but how these processes affected intraplate deformation in eastern Australia is poorly understood. Here we present ⁴⁰Ar/³⁹Ar, K-Ar, and Rb-Sr geochronological data from illitic clay-bearing fault gouges associated with the northern part of the 200 km long, N-striking, dextral strike-slip, Demon Fault in eastern Australia. We show a major range of geochronological ages at 162.99 ± 0.74–152.1 ± 4.8 Ma, indicating that the Demon Fault was active during the Late Jurassic. This period partially coincides with the Middle-Late Jurassic deposition of widespread ash-fall tuffs in the Clarence-Moreton, Surat, and Eromanga basins. We propose that Middle-Late Jurassic intraplate tectonism in eastern Australia was influenced by subduction processes farther east, which produced extensive calc-alkaline magmatism in New Zealand from ~170 Ma. A global plate reorganisation event, related to the development of Early-Middle Jurassic sea-floor spreading of the Pacific Plate, possibly acted as the driving mechanism responsible for the intensification of magmatism and intraplate faulting in eastern Gondwana.     

The 2007 Talala,Saurashtra, western India earthquake sequence: Tectonic implications and seismicity triggering

A damaging and widely felt moderate (M_w 5.0) earthquake occurred in the Talala region of Saurashtra, Gujarat (western India) on November 6, 2007. The highly productive sequence comprised about 1300 micro earthquakes (M > 0.5) out of which 325 of M ? 1.5 that occurred during November 6, 2007–January 10, 2008 were precisely located. The spatial aftershock distribution revealed a NE–SW striking fault in accordance with the centroid moment tensor solution, which in turn implies left-lateral motion. The orientation and sense of shear are consistent with similarly orientated geological fault identified in the area from satellite imagery and field investigation.The aftershocks temporal decay, b-value of frequency–magnitude distribution, spatial fractal dimension, D, and slip ratio (ratio of the slip occurred on the primary fault to the total slip) were examined with the purpose to identify the properties of the sequence. The high b-value (1.18 ± 0.01) may be attributed to the paucity of the larger (M ? 4.0) aftershocks and reveals crustal heterogeneity and low stress regime. The high p-value (1.10 ± 0.39), implying fast decay rate of aftershocks, evidences high surface heat flux. A value of the spatial fractal dimension (D) equal to 2.21 ± 0.02 indicates random spatial distribution and source in a two-dimensional plane that is being filled-up by fractures. A slip ratio of 0.42 reveals that more slip occurred on secondary fault systems.The static Coulomb stress changes due to the coseismic slip of the main shock, enhanced off fault aftershock occurrence. The occurrence of a moderate earthquake (M_w 4.3) on October 5, 2008 inside a region of positive Coulomb stress changes supports the postulation on aftershock triggering. When the stress changes were resolved on a cross section including the stronger (M4.8) foreshock plane that is positioned adjacent to the main fault, it became evident that the activity continued there due to stress transfer from the main rupture.     

Present-day crustal motion around the Pamir Plateau from GPS measurements

We have used geodetic techniques to improve constraints on the crustal motion of the Pamir Plateau. Three campaigns of Global Position System data acquisition between 2011 and 2015 demonstrate that, in association with the India–Asia collision, a complex pattern of crustal motion exists in the Pamir Plateau. In a north–south direction from the Indian Plate to the Hazak Block, the crust has absorbed ~ 35 mm/yr of shortening, of which ~ 35% is distributed around the Hindu Kush region (~ 12 mm/yr), and another ~ 35% is taken up around the Alai Valley (also ~ 12 mm/yr). Global Position System measurements also show ~ 5 mm/yr of shortening between the Pamir Plateau and the Tajik Basin, whereas between the Pamir and the Tarim Basin, an ~ 10 mm/yr extension rate is observed. With respect to the stable Eurasian Plate, the Pamir rotates counterclockwise at a rate of ~ 1.822°Myr^− 1, with an Euler pole positioned about the west end of the Tajik Basin (37.03 ± 0.74°N, 65.89 ± 0.12°E). The strain rate field calculated from Global Position System velocities reveals that the crustal motion is consistent with localized deformation around the Hindu Kush and the Alai Valley, the latter representing a zone with strong shallow seismic activity.     

5.9 Mw, 18th June 2010 earthquake and fault segment linkage at Andaman – A study based on macroseismic survey,GPS geodesy and Coulomb stress changes

18th June, 2010 5.9 M_w earthquake at North Andaman triggered along NW–SE pre-existing fault with reverse fault mechanism. Macroseismic survey and GPS geodesy reveal maximum damages following NE–SW trend due to normal fault mechanism. Coulomb stress modeling for post- and inter-seismic earthquakes after the 2004 mega-earthquake show different stages of fault segment linkage at North Andaman. The present earthquake has been explained as co-shock due to asiesmic soft linkage of fault propagation.     

Slip distribution of the 11 March 2011 Tohoku-Oki earthquake (Mw 9.0) inferred from the refined homomorphic deconvolution method

The kinematic rupture process of the 2011 Tohoku-Oki earthquake (Mw 9.0) was derived given apparent source time functions (ASTFs) retrieved from Rayleigh waves using a refined homomorphic deconvolution method. The total duration of the rupture process was about 165 s. Three slip-concentrated areas were identified based on images of static slip distribution. The largest asperity, located up-dip from the hypocenter with an area of 250 km × 110 km extending to the trench on the fault, had a maximum slip of about 54 m. The other asperities with smaller slip down-dip from the hypocenter were centered on the north and south of the hypocenter, respectively. The preferred average rupture expansion velocity was 1.2 km/s within 130 km from the hypocenter and up to 2.3 km/s over other areas on the fault. Thus, the region near the vicinity of the hypocenter with lower rupture velocity had higher slip amplitude, strongly suggesting brittle failure on a high friction fault. Based on ASTF results, our proposed method offers another way of directly detecting the kinematic source parameters of earthquakes.     

Indenter-tectonics in the Philippines: Example from the Palawan Microcontinental Block - Philippine Mobile Belt Collision

Abstract. The aseismic Palawan microcontinental block is an oceanic bathymetric high that has collided with the seismically-ac-tive Philippine Mobile Belt since the Early Miocene. Consequently, tectonic microblocks immediately north (Luzon) and south (Western Visayas Block) of the collision front rotated in opposite senses. The rotation led the microblocks to onramp adjacent strike-slip faults, and converted these to subduction zones, namely, the current Manila and Negros Trenches. In addition, the collision also initiated the southward propagation of a major left-lateral strike slip fault, the Philippine Fault Zone, and the Philippine Trench, which bounds the Philippine archipelago along its eastern boundary. Based on onshore and offshore data, the Philippine Fault Zone and the East Luzon Trough - Philippine Trench appears to also propagate northward. Furthermore, the opposite direction of propagation is also noted for the Manila and Negros Trenches from the locus of the collision in the Central Philippines to their northern and southern extensions, respectively. The ages of initiation of the Manila Trench (Early Miocene), Philippine Fault Zone (Middle Miocene) and Philippine Trench (Pliocene) as encountered along a west to east transect in the Central Philippines are consistent with the collision and subsequent indentation of Palawan with the rest of the Philippine Mobile Belt.     

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.