New constraints on the current stress field and seismic velocity structure of the eastern Yilgarn Craton from mechanisms of local earthquakes

The Yilgarn Craton has hosted some of the largest earthquakes within the Australian continent in the last 100 years. Earthquakes have mainly been studied in the western part of the craton, and are thought to result from the reactivation of Precambrian structures in an E–W compressive regional stress field imposed by plate-scale processes. Here we present moment tensor solutions for three recent moderate-sized earthquakes around the town of Kalgoorlie that are inconsistent with E–W compression, but instead suggest E–W extension in the eastern Yilgarn Craton. Waveforms of earthquakes at Boulder (M_W = 4.0, 20 April 2010), Kalgoorlie (M_W = 4.3, 26 February 2014) and Coolgardie (M_W = 3.9, 31 October 2014) were inverted for moment tensors. All three earthquakes were shallow (centroid depth ≤4 km) normal-faulting events that occurred along roughly N–S-striking planes, either with a steep westward or a relatively shallow eastward dip. The robustness of the retrieved mechanisms has been thoroughly tested, employing different earth models, assuming different locations for the earthquakes and using different period bands for the inversion. The fit of synthetic long-period waveforms to the observations was in all cases substantially improved by assuming a two-layered crust with high S wavespeeds (about 3.9–4 km/s) overlying substantially slower material. Since there is independent evidence from active source profiles for a P velocity increase between the upper and lower crust, a large difference in v_p/v_s ratio between upper and lower crust is the only way to explain both lines of evidence. This vertical contrast could represent a dominance of felsic material in the upper crust, and substantially more mafic material in the lower crust. Taken together, our results also appear to imply that the regional stress field is E–W extensive in the Kalgoorlie area, and possibly for the entire Kalgoorlie Terrane. This is contrary to current assumptions from continent-scale stress modelling. That the orientations of rupture planes roughly align with the regional structural grain could indicate that Archean structures are reactivated in response to the current stress field.     

3-D GRACE gravity model for the 2011 Japan earthquake

The GRACE mission has contributed to the seismic characterization of major earthquakes in offshore regions of the world. Here, we isolate satellite gravity signal (μGal range) for the Japan Earthquake of 2011 using a difference method. Contrary to the existing gravity models, we propose a unit vertical pyramid based five-layer 3-D thrust fault model, which extends to the hypocenter and honors the ocean water layer and sea floor upheaval also. Our model partly uses existing seismological information (hypocenter depth of 32 km, rupture length of 300 km and vertical slip of 4 m), provides a snapshot of episodic subduction of the Pacific Plate below the Atlantic Plate and its gravity response closely matches the observed gravity (RMS error of 3.4012×10^?13μGal), fully accounting for co-seismic mass redistribution including sea surface deformation. Our inferred rupture length, rupture velocity, average seismic moment magnitude and momentum, respectively, are 300 km, 4.49 km/s, 1.152×10²¹?1.8816×10²¹ N m and 2.319×10⁶ GNs, which fairly agree with the literature. Further, our model inferred momentum at the sea floor corresponds to an area pulse that led to Tsunami generation.     

Teleseismic finite-fault inversion of two <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> = 6.4 earthquakes along the East Anatolian Fault Zone in Turkey: the 1998 Adana and 2003 Bingöl earthquakes

The East Anatolian Fault Zone is a continental transform fault accommodating westward motion of the Anatolian fault. This study aims to investigate the source properties of two moderately large and damaging earthquakes which occurred along the transform fault in the last two decades using the teleseismic broadband P and SH body waveforms. The first earthquake, the 27 June 1998 Adana earthquake, occurred beneath the Adana basin, located close to the eastern extreme of Turkey’s Mediterranean coast. The faulting associated with the 1998 Adana earthquake is unilateral to the NE and confined to depths below 15 km with a length of 30 km along the strike (53°) and a dipping of 81° SE. The fixed-rake models fit the data less well than the variable-rake model. The main slip area centered at depth of about 27 km and to the NE of the hypocenter, covering a circular area of 10 km in diameter with a peak slip of about 60 cm. The slip model yields a seismic moment of 3.5?×?10¹⁸ N-m (M_w???6.4). The second earthquake, the 1 May 2003 Bingöl earthquake, occurred along a dextral conjugate fault of the East Anatolian Fault Zone. The preferred slip model with a seismic moment of 4.1?×?10¹⁸ N-m (M_w???6.4) suggests that the rupture was unilateral toward SE and was controlled by a failure of large asperity roughly circular in shape and centered at a depth of 5 km with peak displacement of about 55 cm. Our results suggest that the 1998 Adana earthquake did not occur on the mapped Göksun Yakap?nar Fault Zone but rather on a SE dipping unmapped fault that may be a split fault of it and buried under the thick (about 6 km) deposits of the Adana basin. For the 2003 Bingöl earthquake, the final slip model requires a rupture plane having 15° different strike than the most possible mapped fault.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

Source parameters of some recent earthquakes in the Gulf of Aqaba, Egypt

Gulf of Aqaba is recognized as an active seismic zone where many destructive earthquakes have occurred. The estimation of source parameters and coda Q attenuation are the main target of this work. Fifty digital seismic events in eight short-period seismic stations with magnitude 2.5–5.2 are used. Most of these events occurred at hypocentral depths in the range of 7–20 km, indicating that the activity was restricted in the upper crust. Seismic moment, M _o, source radius, r, and stress drop, Δσ, are estimated from P- and S-wave spectra using the Brune’s seismic source model. The average seismic moment generated by the whole sequence of events was estimated to be 4.6E?+?22 dyne/cm. The earthquakes with higher stress drop occur at 10-km depth. The scaling relation between the seismic moment and the stress drop indicates a tendency of increasing seismic moment with stress drop. The seismic moment increases with increasing the source radius. Coda waves are sensitive to changes in the subsurface due to the wide scattering effects generating these waves. Single scattering model of local earthquakes is used to the coda Q calculation. The coda with lapse times 10, 20, and 30 s at six central frequencies 1.5, 3, 6, 12, 18, 24 Hz are calculated. The Q _c values are frequency dependent in the range 1–25 Hz, and are approximated by a least squares fit to the power law [ $ {Q_c}(f) = {Q_o}{(f/{f_o})^\eta } $ ]. The average of Q _c values increases from 53?±?10 at 1.5 Hz to 700?±?120 at 24 Hz. The average of Q _o values ranges from 13?±?1 at 1.5 Hz to 39?±?4 at 24 Hz. The frequency exponent parameter η ranges between 1.3?±?0.008 and 0.9?±?0.001.     

Extension and direction of ruptures from teleseismic P-wave patterns for the earthquakes of Gazli 1976, April 8 and Vrancea 1977, March 4

For the earthquakes of Gazli (April 8, 1976) and Vrancea (March 4, 1977) seismogram interpretation was carried out by means of a new technique taking into account the finite dimensions of the source. This interpretation is based on theoretical and experimental considerations of wave patterns generated by a moving rupture. The amplitude of P-waves up to the maximum phases, P_max, were taken into account. The arrival time of P_max on the seismogram was considered as corresponding to the time of the maximum energy release. Source extension and direction of rupture were estimated by means of the azimuthal travel-time curve of the P_max-wave.For the Vrancea earthquake the process at the focus occurred in the form of bi-directional rupturing at azimuths of 180° and 290° from the initial hypocenter. For the Gazli earthquake rupturing developed in several branching directions with azimuths 180°, 240° and 300°. A northwestern direction was prevailing in this event. Records along the rupture and in the opposite direction are shown for both earthquakes.     

Transients of Giggenbach’s ‘Na‐K‐Mg‐Ca Geoindicators’ preceding the 27 October 2004, Mw = 6.0 earthquake in Vrancea area (Romania)

The Na, K, Mg and Ca contents of certain deep‐origin groundwater discharges have been used by Giggenbach (1988) to define a series of ‘geoindicators’, which may provide hints on the up‐flow depth of origin, on the duration of the fluid ascent to the ground surface and on the associated CO₂ flux. On occurrence of a M_w = 6.0 Vrancea earthquake, significant fluctuations of Giggenbach’s geoindicators have been recorded in a saline spring, some 50 km away from the epicentre. A pre‐seismic overall anomaly was monitored for 1 year and a half, the sharpest variations occurring about 3 months before the earthquake. Processes controlling the geoindicator fluctuations assumedly took place at 7–8 km depth, while the earthquake hypocenter depth was about 100 km. This could be an evidence for a mechanical coupling still existing between the seismogenic body in the lithosphere and the overlying crust.     

Earthquake source parameters at the sumatran fault zone: Identification of the activated fault plane

Fifteen earthquakes (M_w 4.1–6.4) occurring at ten major segments of the Sumatran Fault Zone (SFZ) were analyzed to identify their respective fault planes. The events were relocated in order to assess hypocenter uncertainty. Earthquake source parameters were determined from three-component local waveforms recorded by IRIS-DMC and GEOFON broadband lA networks. Epicentral distances of all stations were less than 10°. Moment tensor solutions of the events were calculated, along with simultaneous determination of centroid position. Joint analysis of hypocenter position, centroid position, and nodal planes produced clear outlines of the Sumatran fault planes. The preferable seismotectonic interpretation is that the events activated the SFZ at a depth of approximately 14–210 km, corresponding to the interplate Sumatran fault boundary. The identification of this seismic fault zone is significant to the investigation of seismic hazards in the region.     

Crustal heterogeneities beneath the 2011 Talala,Saurashtra earthquake,Gujarat, India source zone: Seismological evidence for neo-tectonics

During the 1st decade of the 21st century, the study area of Talala, Saurashtra of western India witnessed three damaging earthquakes of moderate magnitude, year 2007 [Mw 5.0; Mw 4.8] and in the year 2011 [Mw 5.1] that generated public panic in the region. The last damaging moderate earthquake of the 20th October 2011 in Talala region (21.09°N;70.45°E), located at about 200 km south to the devastating 2001 Bhuj (23.412°N, 70.232°E) mainshock (Mw 7.6), jolted the entire Saurashtra region of Gujarat. A long series of aftershocks followed hereafter, recorded at nine seismograph/accelerograph stations. Hypocenters of aftershocks were relocated accurately using absolute and relative travel time (double-difference) method. In this study, we, for the first time, determined 3-D tomographic images of the upper crust beneath the 2011 Talala earthquake source zone by inverting about 1135 P and 1125 S wave arrival time data. Estimates of seismic velocities (Vp, Vs) and Poisson’s ratio (σ) structures offer a reliable interpretation of crustal heterogeneities and their bearing on geneses of moderate earthquakes and their aftershock sequences beneath the source zone. It is found that the 2011 Talala mainshock hypocenter depth (6 km) is located near the boundary of the low and high velocity (Vp, Vs) and the source zone is associated with low-σ anomalies guarded by the prominent high-σ anomalies along the active fault zone having strike-slip motion beneath the earthquake source zone. The pattern of distribution of (Vp, Vs, σ) and its association with occurrences of aftershocks provide seismological evidence for the neo-tectonics in the region having left lateral strike-slip motion of the fault.     

The 1980, 1997 and 1998 Azores earthquakes and some seismo-tectonic implications

We have studied the focal mechanisms of the 1980, 1997 and 1998 earthquakes in the Azores region from body-wave inversion of digital GDSN (Global Digital Seismograph Network) and broadband data. For the 1980 and 1998 shocks, we have obtained strike–slip faulting, with the rupture process made up of two sub-events in both shocks, with total scalar seismic moments of 1.9 × 10¹⁹ Nm (M_w = 6.8) and 1.4 × 10¹⁸ Nm (M_w = 6.0), respectively. For the 1997 shock, we have obtained a normal faulting mechanism, with the rupture process made up of three sub-events, with a total scalar seismic moment of 7.7 × 10¹⁷ Nm (M_w = 5.9). A common characteristic of these three earthquakes was the shallow focal depth, less than 10 km, in agreement with the oceanic-type crust. From the directivity function of Rayleigh (LR) waves, we have identified the NW–SE plane as the rupture plane for the 1980 and 1998 earthquakes with the rupture propagating to the SE. Slow rupture velocity, about of 1.5 km/s, has been estimated from directivity function for the 1980 and 1998 earthquakes. From spectral analysis and body-wave inversion, fault dimensions, stress drop and average slip have been estimated. Focal mechanisms of the three earthquakes we have studied, together with focal mechanisms obtained by other authors, have been used in order to obtain a seismotectonic model for the Azores region. We have found different types of behaviour present along the region. It can be divided into two zones: Zone I, from 30°W to 27°W; Zone II, from 27°W to 23°W, with a change in the seismicity and stress direction from Zone I. In Zone I, the total seismic moment tensor obtained corresponded to left-lateral strike–slip faulting with horizontal pressure and tension axes in the E–W and N–S directions, respectively. In Zone II, the total seismic moment tensor corresponded to normal faulting, with a horizontal tension axis trending NE–SW, normal to the Terceira Ridge. The stress pattern for the whole region corresponds to horizontal extension with an average seismic slip rate of 4.4 mm/yr.     

Seismic moment tensor for intermediate depth earthquakes at regional distances in Southern Spain

A method has been developed to obtain the seismic moment tensor components by linear inversion of P waves recorded at regional distance for intermediate depth earthquakes. The seismic moment tensor is separated into double couple (DC) and compensated linear vector dipole (CLVD) parts. The method has been applied to four earthquakes (64<h<95 km) which occurred in the Malaga region (southern Spain). Solutions for the 1987 event show a percentage of CLVD of 20% with a short source time function and DC part corresponding to vertical motion. For the 1989, 1990 and 1992 earthquakes, percentages of CLVD between 0% and 6% have been found. Comparison with the results obtained in a previous study [Buforn et al., J. Seismol. 1 (1997) 113] by modelling of P waves using a DC model, shows that the use of a more general representation of the source (seismic moment tensor) gives a fit of data for the 1987 event.     

LP/HT metamorphism as a temporal marker of change of deformation style within the Late Palaeozoic accretionary wedge of central Chile

A Late Palaeozoic accretionary prism, formed at the southwestern margin of Gondwana from Early Carboniferous to Late Triassic, comprises the Coastal Accretionary Complex of central Chile (34–41°S). This fossil accretionary system is made up of two parallel contemporaneous metamorphic belts: a high‐pressure/low temperature belt (HP/LT – Western Series) and a low pressure/high temperature belt (LP/HT – Eastern Series). However, the timing of deformation events associated with the growth of the accretionary prism (successive frontal accretion and basal underplating) and the development of the LP/HT metamorphism in the shallower levels of the wedge are not continuously observed along this paired metamorphic belt, suggesting the former existence of local perturbations in the subduction regime. In the Pichilemu region, a well‐preserved segment of the paired metamorphic belt allows a first order correlation between the metamorphic and deformational evolution of the deep accreted slices of oceanic crust (blueschists and HP greenschists from the Western Series) and deformation at the shallower levels of the wedge (the Eastern Series). LP/HT mineral assemblages grew in response to arc‐related granitic intrusions, and porphyroblasts constitute time markers recording the evolution of deformation within shallow wedge material. Integrated P–T–t–d analysis reveals that the LP/HT belt is formed between the stages of frontal accretion (D₁) and basal underplating of basic rocks (D₂) forming blueschists at c. 300 Ma. A timeline evolution relating the formation of blueschists and the formation and deformation of LP/HT mineral assemblages at shallower levels, combined with published geochronological/thermobarometric/geochemistry data suggests a cause–effect relation between the basal accretion of basic rocks and the deformation of the shallower LP/HT belt. The S₂ foliation that formed during basal accretion initiated near the base of the accretionary wedge at ~30 km depth at c. 308 Ma. Later, the S₂ foliation developed at c. 300 Ma and ~15 km depth shortly after the emplacement of the granitoids and formation of the (LP/HT) peak metamorphic mineral assemblages. This shallow deformation may reflect a perturbation in the long‐term subduction dynamics (e.g. entrance of a seamount), which would in turn have contributed to the coeval exhumation of the nearby blueschists at c. 300 Ma. Finally, ⁴⁰Ar–³⁹Ar cooling ages reveal that foliated LP/HT rocks were already at ~350 °C at c. 292 Ma, indicating a rapid cooling for this metamorphic system.     

Average seismic source-parameter relations for plate-margin earthquakes

Empirical m_b, M_s and M_o data are used to develop an average spectral scaling relation for plate-margin earthquakes. Using equations based upon a rectangular, bilateral dislocation model with uniform rupture velocity, the spectra give values of fault rupture length and width, static stress drop and average fault displacement as a function of m_b, M_s and M_o. Compared to mid-plate earthquakes of the same seismic moment, the large average plate-margin earthquake has a bigger rupture length, rupture area and average fault displacement and a smaller rupture width and static stress drop.     

Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

Seismic moment release along selected major transforms on the Central Indian Ridge,Western Indian Ocean,and its tectonic implications

A detailed seismicity map of the Central Indian Ridge for the period 1912–1993 is presented, and the earthquakes pertaining to four major transforms offsetting the ridge are utilized to study the moment release pattern. The scalar moment release for the period 1912–1993, and the summed moment rate tensors for both short period (1977–1993) and long period (1912–1993) bring out a unified picture of moment release pattern. The fraction of seismic slip calculated based on depths of 100°C and 400°C limiting temperatures suggests that the Marie-Celeste transform requires a slip almost to a depth of 400°C isotherm to account for the observed moment, and the Argo transform requires depth of faulting much above the 400°C isotherm. A very small fraction of slip is accounted seismically for Vema (53%) and 12° 12′S (23%) even to depths of 100°C isotherm, suggesting a very low order of moment release along these transforms. The horizontal plate velocities and the corresponding strain rates obtained from moment tensor summation of long period data (82 years) give rise to (V _y ^y ; V _y ^x mm. yr⁻¹) of 6.0 and 6.1 along Marie-Celeste, 1.3 and 0.50 along Argo, 0.06 and 0.06 along 12° 12′S, 1.6 and 0.25 along Vema transforms. The corresponding strain rates (ε _y ^y ;ε _y ^x × 10⁻¹⁵ S⁻¹) are 12.7 and 6.8 along MarieCeleste, 6.9 and 1.4 along Argo, 0.27 and 0.14 along 12° 12′S, 7.3 and 0.58 along Vema transforms. These results suggest that the strain rates were highest and almost all predicted motion is taken up seismically along the Marie-Celeste transform. The strain rates are lower along Argo transform and the observed moment release require shallower depth of faulting in order to slip to be accounted seismically. The Vema and 12° 12′S transforms are characterized by low strain rates and less than 15 per cent of motion is accommodated seismically within the seismogenic layer. It is proposed that the deficiency of moment release along the Vema and 12° 12′S multiple transform system may be due to most of the plate motion occurring aseismically.     

Near-surface earthquakes in the Baltic shield, notably Sweden

A recent series of Swedish earthquakes at a focal depth not exceeding 2–3 km, the largest with I_O = V + (MSK scale) andM_L = 3.2 shows that relatively strong seismic activity can occur in the uppermost part of the Baltic Shield. During the last 15 years several near-surface earthquakes have occurred in this region, as indicated by recorded Rg-waves and/or macroseismic data. Many events are located along the coast of central Sweden, suggesting a seismic belt of minor, near-surface activity, which should be considered in connection with the radioactive waste storage in the Swedish bedrock. The appearance of Rg, common in seismic records of explosions and rockbursts, is not a sufficient discriminator between artificial events and earthquakes.     

Determination of 3-D Velocity Anomalies of the Nanbei Tectonic Zone of China Based on Local Earthquakes

In this paper, 3-D velocity images of the crust and upper mantle beneath the Nanbei tectonic zone of China are constructed using P-wave travel time residuals of earthquakes, with the data supplied by China's seismic networks and the International Seismic Centre.During the model parameterization in the tomographic inversion the interpolation function of grid node velocities is used as the space function of velocity, and the velocity function is allowed to be discontinuous. The difficulties caused by large memory demand and high computing cost in solving the system of equations are avoided by utilizing the least squares QR decomposition algorithm. Thus, the stability of the algorithm is ensured. Though resolution images are not derived in the inversion process based on the resolution theory of Backus and Gilbert, a covariance resolution method is put forward by the authors. The resolution analyses have proved the reliability of the image results. 3-D ray tracing is conducted to obtain the ray paths in a s     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Correction to: Impact of climate change on landslides frequency: the Esino river basin case study (Central Italy)

We refine the 1-D velocity model of the Central India Tectonic Zone (CITZ) using well-selected arrival times of P- and S-phases of 354 local earthquakes of magnitude (M_w) between 2.0 and 5.8, recorded by national seismic network from May 1997 to March 2016. Further, we have determined the source mechanisms of 26 selected local events using moment tensor inversion to characterize the dynamics beneath the CITZ. The best-fit simulation between observed and synthetic waveforms obtained the nodal and auxiliary planes of the each faults associated with the earthquake moment magnitude (M_w) for each event. Depth of the fault plane along the CITZ varies from 5 to 38 km. From this study, we found that the western part along the CITZ shows minimum focal depth and reaches maximum 38 kms at Jabalpur in the eastern part. This complex nature of earthquake dynamics occurrence along the CITZ. We propose that the curviplanar the CITZ dominated with sinistral curvature is subjected to compression along the longer ~E–W segments and transtension along shorter segments with ~NE–SW orientations. The occurrences of normal faulting, intrusion of mafic plutons and CLVD mechanisms for earthquakes are interpreted to be linked to the transtension zones and reverse mechanisms associated with the compressions along ~E–W segments.