Seismic potential of Southern Italy

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.     

MODES-2D—Method of detecting strike-shift shear zones

One of the goals of using the Global Positioning System (GPS) and other geodetic survey techniques in tectonics has been to detect boundaries such as faults or shear zones between rigid or mildly deforming crustal masses. The calculation of infinitesimal strains and rotations with GPS data has been widely used to detect shear zones but it has been largely unsuccessful because infinitesimal strain and rotation, although useful in many other ways, is non-diagnostic of shear zones or faults. Our approach is to work with components of deformation, not strain, and to design specifically a diagnostic method of detecting shear zones. This paper introduces the first part of our method, the detection of two-dimensional, strike-shift shear zones (MODES-2D). The MODES-2D method has three elements: (1) determination of the orientation of a suspected strike-shift shear zone by analyzing components of a deformation tensor derived from a data set of displacements in an arbitrary coordinate system; (2) resolution of the deformation tensor into the coordinate system parallel and normal to the detected shear zone; and (3) exploration of the resolved data set for evidence for a belt of inhomogeneous deformation, which is an essential characteristic of a shear zone. The operation of MODES-2D is illustrated herein with a theoretical survey network across an ideal shear zone developed with a buried dislocation-fault and with a survey network afforded by the crossing of the Kaynaşlı viaduct by the 1999 Düzce–Bolu earthquake rupture in Turkey.     

Aseismic deformation in the Alps: GPS vs. seismic strain quantification

Neotectonics of the Western and Central Alps is characterized by ongoing widespread extension in the highest zones of the chain and transcurrent/compressive tectonics at the external limits of the belt. The overall geodetically measured deformations also indicate extension across the Western Alps. There is a good qualitative coherency between seismotectonic and geodetic approaches. Here we attempt to quantify the seismic part of the deformation. The seismic strain is compared to the deformation derived from geodesy. In sub‐areas of homogeneous seismic stress/strain, we computed the total seismic moment tensor and related strain tensor. This study provides new quantitative elements about the ongoing geodynamic processes in the alpine belt. The important discrepancies obtained between seismic strains and geodetically‐measured deformations raise the issue of aseismic deformation in the Alps, which could be related to elastic loading, creeping and/or a slower ductile‐style deformation.     

On the Bayesian analysis of the earthquake hazard in the North-East Indian peninsula

The Bayesian extreme-value distribution of earthquake occurrences has been used to estimate the seismic hazard in 12 seismogenic zones of the North-East Indian peninsula. The Bayesian approach has been used very efficiently to combine the prior information on seismicity obtained from geological data with historical observations in many seismogenic zones of the world. The basic parameters to obtain the prior estimate of seismicity are the seismic moment, slip rate, earthquake recurrence rate and magnitude. These estimates are then updated in terms of Bayes’ theorem and historical evaluations of seismicity associated with each zone. From the Bayesian analysis of extreme earthquake occurrences for North-East Indian peninsula, it is found that for T = 5 years, the probability of occurrences of magnitude (M _w = 5.0–5.5) is greater than 0.9 for all zones. For M _w = 6.0, four zones namely Z1 (Central Himalayas), Z5 (Indo-Burma border), Z7 (Burmese arc) and Z8 (Burma region) exhibit high probabilities. Lower probability is shown by some zones namely␣Z4, Z12, and rest of the zones Z2, Z3, Z6, Z9, Z10 and Z11 show moderate probabilities.     

Spatial-temporal variability of seismic hazard in peninsular India

This paper examines the variability of seismic activity observed in the case of different geological zones of peninsular India (10°N–26°N; 68°E–90°E) based on earthquake catalog between the period 1842 and 2002 and estimates earthquake hazard for the region. With compilation of earthquake catalog in terms of moment magnitude and establishing broad completeness criteria, we derive the seismicity parameters for each geologic zone of peninsular India using maximum likelihood procedure. The estimated parameters provide the basis for understanding the historical seismicity associated with different geological zones of peninsular India and also provide important inputs for future seismic hazard estimation studies in the region. Based on present investigation, it is clear that earthquake recurrence activity in various geologic zones of peninsular India is distinct and varies considerably between its cratonic and rifting zones. The study identifies the likely hazards due to the possibility of moderate to large earthquakes in peninsular India and also presents the influence of spatial rate variation in the seismic activity of this region. This paper presents the influence of source zone characterization and recurrence rate variation pattern on the maximum earthquake magnitude estimation. The results presented in the paper provide a useful basis for probabilistic seismic hazard studies and microzonation studies in peninsular India.     

Seismic activity changes progressing simultaneously with slow-slip in the Tokai area

Anomalous movements were detected simultaneously in both the seismic and the GPS observations in the Tokai area, the central part of the Japanese islands from the late 1990s to 2000. The anomalies are of great concern since the pending risk of a large megathrust earthquake in this area has been predicted for more than 20 years. The GPS data revealed that a slow-slip on the plate interface had commenced beneath Lake Hamana, the center of which is positioned around the edge of the assumed focal zone. On the other hand, the seismic data indicated that a delicate but clear quiescence appeared over a wide area that spreads into the main focal zone. Analyses of the seismicity changes in space and time confirmed that the contrast in the seismicity rate is distinct inside the focal zone. While the integrated seismicity indicated quiescence, some locations were distinguished as activated zones, possibly indicating the appearance of asperities. The rise of the seismicity rate in a quasi-stationary manner suggests an increase in the stress rate at that location. The following hypothesis is proposed based on the simultaneously detected evidences. The slow-slip progressing beneath Lake Hamana will induce a stress shift that invades the interior of the main locked zone, which will increase the contrast of the seismicity rate, possibly reflecting inhomogeneity in the locking strength. Even in this stage, the activated zones still maintain a locked state to prevent overall breakage. Investigations of the b-value changes and of tidal dependence in seismicity that reveal the stress-concentrated state also support the hypothesis. If this is the case, the observed change in seismicity would indicate the process of stress redistribution in the locking state, which represents the preparatory process toward final breakage. Tracking such seismicity changes would yield valid information for predictions of the next Tokai earthquake.     

Co-seismic displacements associated to the Molise (Southern Italy) earthquake sequence of October–November 2002 inferred from GPS measurements

The 2002 earthquake sequence of October 31 and November 1 (main shocks Mw = 5.7) struck an area of the Molise region in Southern Italy. In this paper we analyzed the co-seismic deformation related to the Molise seismic sequence, inferred from GPS data collected before and after the earthquake, that ruptured a rather deep portion of crust releasing a moderate amount of seismic energy with no surface rupture. The GPS data have been reduced using two different processing strategies and softwares (Bernese and GIPSY) to have an increased control over the result accuracy, since the expected surface displacements induced by the Molise earthquake are in the order of the GPS reliability. The surface deformations obtained from the two approaches are statistically equivalent and show a displacement field consistent with the expected deformation mechanism and with no rupture at the surface. In order to relate this observation with the seismic source, an elastic modeling of fault dislocation rupture has been performed using seismological parameters as constraints to the model input and comparing calculated surface displacements with the observed ones. The sum of the seismic moments (8.9 × 10¹⁷ Nm) of the two main events have been used as a constraint for the size and amount of slip on the model fault while its geometry has been constrained using the focal mechanisms and aftershocks locations. Since the two main shocks exhibit the same fault parameters (strike of the plane, dip and co-seismic slip), we modelled a single square fault, size of 15 km × 15 km, assumed to accommodate the whole rupture of both events of the seismic sequence. A vertical E–W trending fault (strike = 266°) has been modeled, with a horizontal slip of 120 mm. Sensitivity tests have been performed to infer the slip distribution at depth. The comparison between GPS observations and displacement vectors predicted by the dislocation model is consistent with a source fault placed between 5 and 20 km of depth with a constant pure right-lateral strike-slip in agreement with fault slip distribution analyses using seismological information. The GPS strain field obtained doesn't require a geodetic moment release larger than the one inferred from the seismological information ruling out significant post-seismic deformation or geodetic deformation released at frequencies not detectable by seismic instruments. The Molise sequence has a critical seismotectonic significance because it occurred in an area where no historical seismicity or seismogenic faults are reported. The focal location of the sequence and the strike-slip kinematics of main shocks allow to distinguish it from the shallower and extensional seismicity of the southern Apennines being more likely related to the decoupling of the southern Adriatic block from the northern one.     

Global strain rates in western to central Himalayas and their implications in seismic hazard assessment

The Himalayas has experienced varying rates of earthquake occurrence in the past in its seismo-tectonically distinguished segments which may be attributed to different physical processes of accumulation of stress and its release, and due diligence is required for its inclusion for working out the seismic hazard. The present paper intends to revisit the various earthquake occurrence models applied to Himalayas and examines it in the light of recent damaging earthquakes in Himalayan belt. Due to discordant seismicity of Himalayas, three types of regions have been considered to estimate larger return period events. The regions selected are (1) the North-West Himalayan Fold and Thrust Belt which is seismically very active, (2) the Garhwal Himalaya which has never experienced large earthquake although sufficient stress exists and (3) the Nepal region which is very seismically active region due to unlocked rupture and frequently experienced large earthquake events. The seismicity parameters have been revisited using two earthquake recurrence models namely constant seismicity and constant moment release. For constant moment release model, the strain rates have been derived from global strain rate model and are converted into seismic moment of earthquake events considering the geometry of the finite source and the rates being consumed fully by the contemporary seismicity. Probability of earthquake occurrence with time has been estimated for each region using both models and compared assuming Poissonian distribution. The results show that seismicity for North-West region is observed to be relatively less when estimated using constant seismicity model which implies that either the occupied accumulated stress is not being unconfined in the form of earthquakes or the compiled earthquake catalogue is insufficient. Similar trend has been observed for seismic gap area but with lesser difference reported from both methods. However, for the Nepal region, the estimated seismicity by the two methods has been found to be relatively less when estimated using constant moment release model which implies that in the Nepal region, accumulated strain is releasing in the form of large earthquake occurrence event. The partial release in second event of May 2015 of similar size shows that the physical process is trying to release the energy with large earthquake event. If it would have been in other regions like that of seismic gap region, the fault may not have released the energy and may be inviting even bigger event in future. It is, therefore, necessary to look into the seismicity from strain rates also for its due interpretation in terms of predicting the seismic hazard in various segments of Himalayas.     

2015年尼泊尔8.1级大地震前后喜马拉雅带地震活动图像演变

针对2015年4月25日尼泊尔M8.1地震后喜马拉雅造山带的未来地震危险性问题,通过对喜马拉雅带历史大地震应变能释放和在尼泊尔地震发震前后的区域地震活动图像进行了分析研究。结果发现喜马拉雅带很可能已进入新-轮的地震活跃期。此次尼泊尔大地震不足以将喜马拉雅带中段的地壳应变能全部释放,喜马拉雅带中段的地震活动和藏南裂谷带地震活动具有密切的关联,在喜马拉雅带中段和藏南裂谷带还将有大地震活动。同时研究结果还显示现今在喜马拉雅带的东段存在阿萨姆围空区和不丹围空区,在喜马拉雅的西段出现噶尔围空区,喜马拉雅西段新德里和西藏接壤地区以及喀喇昆仑断裂上噶尔县地区地震危险性很高,喜马拉雅东段林芝山南地区以南的阿萨姆和不丹地区危险性很高,应引起重视。     

Recent crustal deformation and seismicity studies at Abu-Dabbab area, Red Sea, Egypt

The Abu-Dabbab area is characterized by high seismicity and complex tectonic setting; for these facts, a local geodetic network consisting of 11 geodetic benchmarks has been established. The crustal deformation data in this area are collected using the GPS techniques. Five campaigns of GPS measurements have been collected, processed, and adjusted to get the more accurate positions of the GPS stations. The horizontal velocity vectors, the dilatational, the maximum shear strains, and the principal strain rates were estimated. The horizontal velocity varies in average between 3 and 6 mm per year across the network. The results of the deformation analysis indicate a significant contraction and extension across the southern central part of the study area which is characterized by high seismic activity represented by the clustering shape of the microearthquakes that trending ENE. The north and northeastern parts are characterized by small strain rates. This study is an attempt to provide valuable information about the present state of the crustal deformation and its relationship to seismic activity and tectonic setting at the Abu-Dabbab area. The present study is the first work demonstrating crustal deformation monitoring at the Abu-Dabbab area. The time interval is relatively short. Actually, these results are preliminary results. So, the continuity of GPS measurements is needed for providing more information about the recent crustal deformation in that area.     

Seismic moment release data in earthquake catalogue: Application of Hurst statistics in delineating temporal clustering and seismic vulnerability

Sequential cumulative moment release data of macroearthquakes (Mw≥4.3) of seventeen seismic zones (A to Q) belonging to NE-Himalaya, Burmese-Andaman arc and West- Sunda arc are analysed by Hurst analysis, a non-parametric statistical procedure to identify clustering of low and high values in a time series. The moment release in a zone occurs in alternate positive, negative and positive sloping segments forming a wave like pattern with intervening small horizontal segment. The negative sloping segments indicate decelerated moment release pattern or temporal slackening of elastic strain release with high b–value (>0.95). The horizontal segment indicates temporal clustering of moderate magnitude events/seismic moments with moderate b-values (0.8–0.95). The positive segment is characterised by accelerated moment release within a short span of time indicating temporal clustering of larger magnitude earthquakes/seismic moments and exhibit lowest b–value (<0.7). All zones attest moderate to high Hurst K values, range 0.7-0.86. The pattern in Hurst plots, specially a reversal of trend after prolong negative slope is used for earthquake prognostication in the seismic zones. Our analysis shows that most of the zones register a notable reversal of Hurst clustering trend after a prolonged negative slope which is accompanied by a major earthquake near its end. However, South Burma region (Zone-I) and Tripura fold belt and Bangladesh Plain (Zone-K) do not show any moderate or large shock around the end of the negative sloping trend in Hurst plot. Hence, these two zones can be considered more prone to produce moderate to larger earthquakes in future.     

Coexisting tectonic settings: the example of the southern Tyrrhenian Sea

We performed geodetic strain rate analyses in southern Italy, using new GPS velocities. Two-dimensional strain and rotation rate fields were estimated and results show that most of the shortening is distributed in the northern Sicily offshore. Extension becomes more evident and comparable with shortening on the eastern side of the same margin, and greater in the eastern Sicily offshore. Principal shortening and extension rate axes are consistent with long-term geological features: seismic reflection profiles show both active compressive and extensional faults affecting Pleistocene strata. We show evidence for contemporaneous extension and transtension in the Cefalù Basin. Combining geodetic data and geological features point to the coexistence of independent geodynamic processes, i.e., the active E–W backarc spreading in the hangingwall of the Apennines subduction zone and shortening along the southern margin of the Tyrrhenian backarc basin operated by the NNW-motion of Africa relative to Eurasia.     

2008年汶川大地震孕震、同震及震后变形和应力演化全过程的数值模拟

2008年MS 8.0级汶川大地震发生在具有复杂的地质构造背景、 强烈的地表起伏、 不均匀的弹性和黏性结构的龙门山断裂带上.由于震前地震活动性不够强烈且地表构造变形较小,龙门山断裂带的地震危险性在汶川地震之前被低估.从数值模拟的角度,建立黏弹性有限元模型,考虑了初始地形、 重力、构造加载、 黏弹性松弛等因素对2008年...     

郯庐断裂带中南段闭锁特征与地震危险性分析

利用华北地区1999～2007、2013～2017两期GPS水平运动速度场数据,采用块体负位错模型,分别反演了郯庐断裂带中南段不同段的断层闭锁程度和滑动亏损速率分布;结合地表应变结果,综合分析了郯庐断裂带前后两期的变形差异特征,并探讨了其与日本3·11地震间的可能关系。研究结果表明:日本地震后,郯庐断裂带中南段郯城以北的段落闭锁程度有所减弱,中南段东部地区主张应变率增强,处于拉张状态;日本大地震的发生对郯庐断裂带中南段的应变积累起到一定的缓解作用。2013～2017最新一期反演结果显示莒县以北断层闭锁程度仍较高,闭锁深度较深,为右旋挤压亏损,是1668年郯城地震的未破裂段;莒县以南到泗洪附近断层闭锁程度较低,无滑动亏损积累;泗洪以南到嘉山段断层闭锁程度较高,是历史地震的未破裂段,同时该地区小震不活跃,易于应力积累,地震危险性值得关注。     

祁连山北缘—河西走廊地区大地形变与地震的关系

河西走廊—祁连山北缘地区地处青藏高原北缘,受印度板块与欧亚板块中生代末—新生代早期的碰撞及持续至今的向北推挤作用的远程效应的影响,该地区是现今的地壳活动地区,其中地壳形变是最主要的表现形式。地壳形变监测显示,隆起区垂直位移速率最大可达15mm/a,沉降区最大位移速率为-15mm/a。祁连山和河西走廊的相对隆升变化与该区地震具有密切的关系,河西走廊相对下降、祁连山相对隆升的后期是地震多发时期,河西走廊相对隆升、祁连山相对下降的后期是地震少发时期,这与该区处于挤压体制下的区域构造背景密切相关。GPS水平位移监测显示,河西走廊—祁连山北缘地区全区都一致向东位移,且位移速率非常大,大者大于10mm/a;位移速率具有南部大于北部、东部大于西部的特点,水平位移速率变化与现代活动断裂具有非常密切的关系,并以主要断裂构造为区带的边界;水平位移速率矢量与2002年玉门地震的震源机制解所显示的沿地震破裂面发生的滑动方向非常一致。     

Transtensional deformation in the Lake Tahoe region, California and Nevada, USA

Dextral transtensional deformation is occurring along the Sierra Nevada–Great Basin boundary zone (SNGBBZ) at the eastern edge of the Sierra Nevada microplate. In the Lake Tahoe region of the SNGBBZ, transtension is partitioned spatially and temporally into domains of north–south striking normal faults and transitional domains with conjugate strike-slip faults. The normal fault domains, which have had large Holocene earthquakes but account only for background seismicity in the historic period, primarily accommodate east–west extension, while the transitional domains, which have had moderate Holocene and historic earthquakes and are currently seismically active, primarily record north–south shortening. Through partitioned slip, the upper crust in this region undergoes overall constrictional strain.Major fault zones within the Lake Tahoe basin include two normal fault zones: the northwest-trending Tahoe–Sierra frontal fault zone (TSFFZ) and the north-trending West Tahoe–Dollar Point fault zone. Most faults in these zones show eastside down displacements. Both of these fault zones show evidence of Holocene earthquakes but are relatively quiet seismically through the historic record. The northeast-trending North Tahoe–Incline Village fault zone is a major normal to sinistral-oblique fault zone. This fault zone shows evidence for large Holocene earthquakes and based on the historic record is seismically active at the microearthquake level. The zone forms the boundary between the Lake Tahoe normal fault domain to the south and the Truckee transition zone to the north.Several lines of evidence, including both geology and historic seismicity, indicate that the seismically active Truckee and Gardnerville transition zones, north and southeast of Lake Tahoe basin, respectively, are undergoing north–south shortening. In addition, the central Carson Range, a major north-trending range block between two large normal fault zones, shows internal fault patterns that suggest the range is undergoing north–south shortening in addition to east–west extension.A model capable of explaining the spatial and temporal partitioning of slip suggests that seismic behavior in the region alternates between two modes, one mode characterized by an east–west minimum principal stress and a north–south maximum principal stress as at present. In this mode, seismicity and small-scale faulting reflecting north–south shortening concentrate in mechanically weak transition zones with primarily strike-slip faulting in relatively small-magnitude events, and domains with major normal faults are relatively quiet. A second mode occurs after sufficient north–south shortening reduces the north–south S_hmax in magnitude until it is less than S_v, at which point S_v becomes the maximum principal stress. This second mode is then characterized by large earthquakes on major normal faults in the large normal fault domains, which dominate the overall moment release in the region, producing significant east–west extension.     

Spatio-temporal stress states estimated from seismicity rate changes in the Tokai region, central Japan

Since unprecedented large-scale silent slip was detected by GPS in 2001 in the Tokai region, evaluating whether such movement is uniquely connected to the expected Tokai earthquake or repeatedly occurs in this area becomes vitally important. Because of short history of GPS observations and the limited areal coverage surrounding the Suruga trough, we take advantage of continuously recorded seismicity that is presumed to be sensitive to the deformation at seismogenic depth. Together with the well-maintained NIED earthquake data, we employ the seismicity-to-stress inversion approach of rate/state friction to infer the spatio-temporal stress changes in and around the presumed hypocentral zone of the future Tokai earthquake. Mapping stress changes inverted from microseismicity year by year, we find that the stress under Lake Hamana, the western expected future Tokai source, has been decreasing since 1999, during which the GPS data showed a normal trend of plate coupling. In contrast, stresses in the surrounding regions are calculated to have increased by transfer from Lake Hamana region. We interpret that this continuous process is associated with the 2000–2004 Tokai slow slip event. The characteristic patterns related to aseismic stress-release are also identified in the early 1980s and during 1987–1989, when slow events are inferred to have occurred on the basis of conventional geodetic measurements. Revisiting the seismotectonics and taking into account the mechanical implications of the inversion results, we argue that the transition zone situated between a deep stable creeping zone and a locked zone undergoes episodic creep and plays an important role in the transfer of stress to the locked zone. Consequently, even though we speculate that the current (2000 to present-day) silent slip event might be one of the repeating events, the inferred enlargement of the stress releasing area is significant and possibly raises the likelihood of the next Tokai earthquake.     

Crustal deformation in northern Central America

Evaluation of the seismic moment tensor for earthquakes on plate boundary is a standard procedure to determine the relative velocity of plates, which controls the seismic deformation rate predicted from the slip on a single fault. The moment tensor is also decomposed into an isotropic and a deviatoric part to discover the relationship between the average strain rate and the relative velocity between two plates. We utilize this procedure to estimate the rates of deformation in northern Central America where plate boundaries are seismically well defined. Four different tectonic environments are considered for modelling of the plate motions. The deformation rates obtained here compare well with those predicted from the plate motions models and are in good agreement with actual observations. Deformation rates on faults are increasingly being used to estimate earthquake recurrence from information on fault slip rate and more on how we can incorporate our current understanding into seismic hazard analyses.     

Double-sided subduction with contrasting polarities beneath the Pamir-Hindu Kush: Evidence from focal mechanism solutions and stress field inversion

The Pamir-Hindu Kush region at the western end of the Himalayan-Tibet orogen is one of the most active regions on the globe with strong seismicity and deformation and provides a window to evaluate continental collision linked to two intra-continental subduction zones with different polarities. The seismicity and seismic tomography data show a steep northward subducting slab beneath the Hindu Kush and southward subducting slab under the Pamir. Here, we collect seismic catalogue with 3988 earthquake events to compute seismicity images and waveform data from 926 earthquake events to invert focal mechanism solutions and stress field with a view to characterize the subducting slabs under the Pamir-Hindu Kush region. Our results define two distinct seismic zones: a steep one beneath the Hindu Kush and a broad one beneath the Pamir. Deep and intermediate-depth earthquakes are mainly distributed in the Hindu Kush region which is controlled by thrust faulting, whereas the Pamir is dominated by strike-slip stress regime with shallow and intermediate-depth earthquakes. The area where the maximum principal stress axis is vertical in the southern Pamir corresponds to the location of a high-conductivity low-velocity region that contributes to the seismogenic processes in this region. We interpret the two distinct seismic zones to represent a double-sided subduction system where the Hindu Kush zone represents the northward subduction of the Indian plate, and the Pamir zone shows southward subduction of the Eurasian plate. A transition fault is inferred in the region between the Hindu Kush and the Pamir which regulates the opposing directions of motion of the Indian and Eurasian plates.     

Late Cretaceous-Paleogene deepwater basin of North Afghanistan and the Central Pamirs: Issue of Hindu Kush earthquakes

The within-Iranian backarc basins, including the largest Sebzawar Basin, opened in the Mid-Cretaceous. Spreading in this basin was completed by the end of the Cretaceous. The basin closed in the Eocene with the formation of subduction zones and volcanic-plutonic belts. Data on North Afghanistan and the Central Pamirs have allowed us to reconstruct the eastern continuation of the Sebzawar Basin up to the west of the Central Pamirs. No fragments of oceanic crust are retained in Afghanistan and the Pamirs, but by analogy with the Sebzawar Basin, thick Paleogene flysch sequences and volcanic-plutonic complexes indicate setting of the active margin and subduction. It is suggested that the belt of mantle seismicity that extends for 550 km to the south of the Central Pamirs is related to the plunging and deformation of the lithosphere once underlying the Cretaceous-Paleogene basin. The extremely vigorous seismicity of the Hindu Kush megasource at the western termination of the seismic belt is caused by a number of specific tectonic features that predetermined the early onset of plunging of the subducted sheet (slab). In the megasource, the slab sank to a depth of 300 km and became vertical; its active deformation has proceeded up to the present. In the eastern part of the seismic belt, the slab started to plunge much later and therefore has retained a gentle slope, so that the depth of the hypocenters is shallower (down to 200 km), and earthquakes are less strong.