Climate extremes: an observation and projection of its impacts on food production in ASEAN

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.     

Precursory seismicity rate changes prior to the large and major earthquakes along the Sagaing fault zone,Central Myanmar

In this study, the seismicity rate changes that can represent an earthquake precursor were investigated along the Sagaing Fault Zone (SFZ), Central Myanmar, using the Z value technique. After statistical improvement of the existing seismicity data (the instrumental earthquake records) by removal of the foreshocks and aftershocks and man-made seismicity changes and standardization of the reported magnitude scales, 3574 earthquake events with a M _w ≥ 4.2 reported during 1977–2015 were found to directly represent the seismotectonic activities of the SFZ. To find the characteristic parameters specifically suitable for the SFZ, seven known events of M _w ≥ 6.0 earthquakes were recognized and used for retrospective tests. As a result, utilizing the conditions of 25 fixed earthquake events considered (N) and a 2-year time window (T _w), a significantly high Z value was found to precede most of the M _w ≥ 6.0 earthquakes. Therefore, to evaluate the prospective areas of upcoming earthquakes, these conditions (N = 25 and T _w = 2) were applied with the most up-to-date seismicity data of 2010–2015. The results illustrate that the vicinity of Myitkyina and Naypyidaw (Z = 4.2–5.1) cities might be subject to strong or major earthquakes in the future.     

Magnitude conversion problem for the Turkish earthquake data

Earthquake catalogues which form the main input in seismic hazard analysis generally report earthquake magnitudes in different scales. Magnitudes reported in different scales have to be converted to a common scale while compiling a seismic data base to be utilized in seismic hazard analysis. This study aims at developing empirical relationships to convert earthquake magnitudes reported in different scales, namely, surface wave magnitude, M _S, local magnitude, M _L, body wave magnitude, m _b and duration magnitude, M _d, to the moment magnitude (M _w). For this purpose, an earthquake data catalogue is compiled from domestic and international data bases for the earthquakes occurred in Turkey. The earthquake reporting differences of various data sources are assessed. Conversion relationships are established between the same earthquake magnitude scale of different data sources and different earthquake magnitude scales. Appropriate statistical methods are employed iteratively, considering the random errors both in the independent and dependent variables. The results are found to be sensitive to the choice of the analysis methods.     

Seismicity of Gujarat

Paper describes tectonics, earthquake monitoring, past and present seismicity, catalogue of earthquakes and estimated return periods of large earthquakes in Gujarat state, western India. The Gujarat region has three failed Mesozoic rifts of Kachchh, Cambay, and Narmada, with several active faults. Kachchh district of Gujarat is the only region outside Himalaya-Andaman belt that has high seismic hazard of magnitude 8 corresponding to zone V in the seismic zoning map of India. The other parts of Gujarat have seismic hazard of magnitude 6 or less. Kachchh region is considered seismically one of the most active intraplate regions of the World. It is known to have low seismicity but high hazard in view of occurrence of fewer smaller earthquakes of M????6 in a region having three devastating earthquakes that occurred during 1819 (M _w7.8), 1956 (M _w6.0) and 2001 (M _w7.7). The second in order of seismic status is Narmada rift zone that experienced a severely damaging 1970 Bharuch earthquake of M5.4 at its western end and M????6 earthquakes further east in 1927 (Son earthquake), 1938 (Satpura earthquake) and 1997 (Jabalpur earthquake). The Saurashtra Peninsula south of Kachchh has experienced seismicity of magnitude less than 6.     

The critical behaviour of a power law between preseismic electric signals and earthquakes of different mechanism in Greece and Japan

The consistency of the critical exponent in the power law relation between the stress drop of the earthquake and the lead time of the precursory seismic electric signal is checked using new data from the recent M _w 4.9 earthquake of strike-slip mechanism that occurred on 12 November 2013 in northern Evia island, Greece and the megathrust M _w 9.0 Tohoku earthquake on 11 March 2011, in Japan. For the first case, the derived exponent is in excellent agreement with previous ones obtained from all non thrust events analysed by the author and matches the value of critical exponent for fracture. On the other hand, the megathrust Tohoku earthquake follows the behaviour of all thrust events studied by the author, and thus, the calculated exponent significantly deviates from this critical value. The different behaviour between non thrust and thrust-type events could be attributed to the fact that thrust mechanism earthquakes usually occur in collision or subduction zones which are characterised by high accumulation of strain. However, a larger number of thrust events are required in order to obtain reliable results and shed light in the above experimental findings.     

A complete and homogeneous magnitude earthquake catalogue of Iraq

A complete and homogeneous magnitude earthquake catalogue spanning the period 1900 to 2010 was created. The catalogue covers the area 29° to 37.5° N and 39° to 48° E. Entries in the new earthquake catalogue were cross checked and additions made from various sources of earthquake records to ensure that repetitions are not included in this analysis. Events were considered duplicates if they had a time difference of 10 s or less and space origin difference of 0.5° or less. In a given set of duplicate events, an event, which had a magnitude and International Seismological Center source, was retained as the record of the event. The unified magnitude scale, the moment magnitude (M _w), was applied throughout the catalogue. The M _w for 18 events was reported. The M _w for other events was estimated using empirical relations between m _b, M _s, M _L, and M _w. Magnitude of completeness, M _c, was estimated using the maximum curvature. It was 4.3 M _w. Finally, a list of 213 events from 1900 to 2010 with M _w?≥?4.3 is presented. The list is considered complete for the period from 1962 to 2010.     

Andean earthquakes triggered by the 2010 Maule,Chile (Mw 8.8) earthquake: Comparisons of geodetic,seismic and geologic constraints

The Maule, Chile, (M_w 8.8) earthquake on 27 February 2010 triggered deformation events over a broad area, allowing investigation of stress redistribution within the upper crust following a mega-thrust subduction event. We explore the role that the Maule earthquake may have played in triggering shallow earthquakes in northwestern Argentina and Chile. We investigate observed ground deformation associated with the M_w 6.2 (GCMT) Salta (1450 km from the Maule hypocenter, 9 h after the Maule earthquake), M_w 5.8 Catamarca (1400 km; nine days), M_w 5.1 Mendoza (350 km; between one to five days) earthquakes, as well as eight additional earthquakes without an observed geodetic signal. We use seismic and Interferometric Synthetic Aperture Radar (InSAR) observations to characterize earthquake location, magnitude and focal mechanism, and characterize how the non-stationary, spatially correlated noise present in the geodetic imagery affects the accuracy of our parameter estimates. The focal mechanisms for the far-field Salta and Catamarca earthquakes are broadly consistent with regional late Cenozoic fault kinematics. We infer that dynamic stresses due to the passage of seismic waves associated with the Maule earthquake likely brought the Salta and Catamarca regions closer to failure but that the involved faults may have already been at a relatively advanced stage of their seismic cycle. The near-field Mendoza earthquake geometry is consistent with triggering related to positive static Coulomb stress changes due to the Maule earthquake but is also aligned with the South America-Nazca shortening direction. None of the earthquakes considered in this study require that the Maule earthquake reactivated faults in a sense that is inconsistent with their long-term behavior.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

Tsunami hazard evaluation for Kuwait and Arabian Gulf due to Makran Subduction Zone and Subaerial landslides

Given the recent historical disastrous tsunamis and the knowledge that the Arabian Gulf (AG) is tectonically active, this study aimed to evaluate tsunami hazards in Kuwait from both submarine earthquakes and subaerial landslides. Despite the low or unknown tsunami risks that impose potential threats to the coastal area’s infrastructures and population of Kuwait, such an investigation is important to sustain the economy and safety of life. This study focused on tsunamis generated by submarine earthquakes with earthquake magnitudes (M _w ) of 8.3–9.0 along the Makran Subduction Zone (MSZ) and subaerial landslides with volumes of 0.75–2.0 km³ from six sources along the Iranian coast inside the AG and one source at the Gulf entrance in Oman. The level of tsunami hazards associated with these tsunamigenic sources was evaluated using numerical modeling. Tsunami model was applied to conduct a numerical tsunami simulation and predict tsunami propagation. For landslide sources, a two-layer model was proposed to solve nonlinear longwave equations within two interfacing layers with appropriate kinematic and dynamic boundary conditions. Threat level maps along the coasts of the AG and Kuwait were developed to illustrate the impacts of potential tsunamis triggered by submarine earthquakes of different scales and subaerial landslides at different sources. GEBCO 30 arc-second grid data and others were used as bathymetry and topography data for numerical modeling. Earthquakes of M _w 8.3 and M _w 8.6 along the MSZ had low and considerable impacts, respectively, at the Gulf entrance, but negligible impacts on Kuwait. An earthquake of M _w 9.0 had a remarkable impact for the entire Gulf region and generated a maximum tsunami amplitude of up to 0.5 m along the Kuwaiti coastline 12 h after the earthquake. In the case of landslides inside the AG, the majority impact occurred locally near the sources. The landslide source opposite to Kuwait Bay generated the maximum tsunami amplitudes reaching 0.3 m inside Kuwait Bay and 1.8 m along the southern coasts of Kuwait.

     

The record of historic earthquakes in lake sediments of Central Switzerland

Deformation structures in lake sediments in Central Switzerland can be attributed to strong historic earthquakes. The type and spatial distribution of the deformation structures reflect the historically documented macroseismic intensities thus providing a useful calibration tool for paleoseismic investigations in prehistoric lake sediments.The Swiss historical earthquake catalogue shows four moderate to strong earthquakes with moment magnitudes of M_w=5.7 to M_w=6.9 and epicentral intensities of I₀=VII to I₀=IX that affected the area of Central Switzerland during the last 1000 years. These are the 1964 Alpnach, 1774 Altdorf, 1601 Unterwalden, and 1356 Basel earthquakes. In order to understand the effect of these earthquakes on lacustrine sediments, four lakes in Central Switzerland (Sarner See, Lungerer See, Baldegger See, and Seelisberg Seeli) were investigated using high-resolution seismic data and sediment cores. The sediments consist of organic- and carbonate-rich clayey to sandy silts that display fine bedding on the centimeter to millimeter scale. The sediments are dated by historic climate and environmental records, ¹³⁷Cs activity, and radiocarbon ages. Deformation structures occur within distinct zones and include large-scale slumps and rockfalls, as well as small-scale features like disturbed and contorted lamination and liquefaction structures. These deformations are attributed to three of the abovementioned earthquakes. The spatial distribution of deformation structures in the different lakes clearly reflects the historical macroseismic dataset: Lake sediments are only affected if they are situated within an area that underwent groundshaking not smaller than intensity VI to VII. We estimate earthquake size by relating the epicentral distance of the farthest liquefaction structure to earthquake magnitude. This relationship is in agreement with earthquake size estimations based on the historical dataset.     

Time-dependent analysis of aftershock events and structural impacts on intraplate crustal seismicity of the Van earthquake (Mw 7.1, 23 October 2011), E-Anatolia

The Van earthquake (M _W 7.1, 23 October 2011) in E-Anatolia is typical representative of intraplate earthquakes. Its thrust focal character and aftershock seismicity pattern indicate the most prominent type of compound earthquakes due to its multifractal dynamic complexity and uneven compressional nature, ever seen all over Turkey. Seismicity pattern of aftershocks appears to be invariably complex in its overall characteristics of aligned clustering events. The population and distribution of the aftershock events clearly exhibit spatial variability, clustering-declustering and intermittency, consistent with multifractal scaling. The sequential growth of events during time scale shows multifractal behavior of seismicity in the focal zone. The results indicate that the extensive heterogeneity and time-dependent strength are considered to generate distinct aftershock events. These factors have structural impacts on intraplate seismicity, suggesting multifractal and unstable nature of the Van event. Multifractal seismicity is controlled by complex evolution of crustal-scale faulting, mechanical heterogeneity and seismic deformation anisotropy. Overall seismicity pattern of aftershocks provides the mechanism for strain softening process to explain the principal thrusting event in the Van earthquake. Strain localization with fault weakening controls the seismic characterization of Van earthquake and contributes to explain the anomalous occurrence of aftershocks and intraplate nature of the Van earthquake.     

Tectonic structures across the East African Rift based on the source parameters of the 20 May 1990 M7.2 Sudan earthquake

Earthquakes in Kenya are common along the Kenya Rift Valley because of the slow divergent movement of the rift and hydrothermal processes in the geothermal fields. This implies slow but continuous radiation of seismic energy, which relieves stress in the subsurface rocks. On the contrary, the NW-SE trending rift/fault zones such as the Aswa-Nyangia fault zone and the Muglad-Anza-Lamu rift zone are the likely sites of major earthquakes in Kenya and the East African region. These rift/fault zones have been the sites of a number of strong earthquakes in the past such as the M _w = 7.2 southern Sudan earthquake of 20 May 1990 and aftershocks of M _w = 6.5 and 7.1 on 24 May 1990, the 1937 M _s = 6.1 earthquake north of Lake Turkana close to the Kenya-Ethiopian border, and the 1913 M _s = 6.0 Turkana earthquake, among others. Source parameters of the 20 May 1990 southern Sudan earthquake show that this earthquake consists of only one event on a fault having strike, dip, and rake of 315°, 84°, and ?3°. The fault plane is characterized by a left-lateral strike slip fault mechanism. The focal depth for this earthquake is 12.1 km, seismic moment M _o = 7.65 × 10¹⁹ Nm, and moment magnitude, M _w = 7.19 (?7.2). The fault rupture started 15 s earlier and lasted for 17 s along a fault plane having dimensions of ?60 km × 40 km. The average fault dislocation is 1.1 m, and the stress drop, ?σ, is 1.63 MPa. The distribution of historical earthquakes (M _w ≥ 5) from southern Sudan through central Kenya generally shows a NW-SE alignment of epicenters. On a local scale in Kenya, the NW–SE alignment of epicenters is characterized by earthquakes of local magnitude M _l ≤ 4.0, except the 1928 Subukia earthquake (M _s = 6.9) in central Kenya. This NW–SE alignment of epicenters is consistent with the trend of the Aswa-Nyangia Fault Zone, from southern Sudan through central Kenya and further southwards into the Indian Ocean. We therefore conclude that the NW–SE trending rift/fault zones are sites of strong earthquakes likely to pose the greatest earthquake hazard in Kenya and the East African region in general.     

Probabilities for the occurrences of medium to large earthquakes in northeast India and adjoining region

The return periods and occurrence probabilities related to medium and large earthquakes (M _w 4.0–7.0) in four seismic zones in northeast India and adjoining region (20°–32°N and 87°–100°E) have been estimated with the help of well-known extreme value theory using three methods given by Gumbel (1958), Knopoff and Kagan (1977) and Bury (1999). In the present analysis, the return periods, the most probable maximum magnitude in a specified time period and probabilities of occurrences of earthquakes of magnitude M ≥ 4.0 have been computed using a homogeneous and complete earthquake catalogue prepared for the period between 1897 and 2007. The analysis indicates that the most probable largest annual earthquakes are close to 4.6, 5.1, 5.2, 5.5 and 5.8 in the four seismic zones, namely, the Shillong Plateau Zone, the Eastern Syntaxis Zone, the Himalayan Thrusts Zone, the Arakan-Yoma subduction zone and the whole region, respectively. The most probable largest earthquakes that may occur within different time periods have been also estimated and reported. The study reveals that the estimated mean return periods for the earthquake of magnitude M _w 6.5 are about 6–7 years, 9–10 years, 59–78 years, 72–115 years and 88–127 years in the whole region, the Arakan-Yoma subduction zone, the Himalayan Thrusts Zone, the Shillong Plateau Zone and the Eastern Syntaxis Zone, respectively. The study indicates that Arakan-Yoma subduction zone has the lowest mean return periods and high occurrence probability for the same earthquake magnitude in comparison to the other zones. The differences in the hazard parameters from zone to zone reveal the high crustal heterogeneity and seismotectonics complexity in northeast India and adjoining regions.     

Seismic hazard across Bulgaria and neighbouring areas: regional and site-specific maximum credible magnitudes and earthquake perceptibility

A probabilistic seismic hazard assessment is developed here using maximum credible earthquake magnitude statistics and earthquake perceptibility hazard. Earthquake perceptibility hazard is defined as the probability a site perceives ground shaking equal to or greater than a selected ground motion level X, resulting from an earthquake of magnitude M, and develops estimates for the most perceptible earthquake magnitude, M _P(max). Realistic and usable maximum magnitude statistics are obtained from both whole process and part process statistical recurrence models. These approaches are extended to develop relationships between perceptible earthquake magnitude hazard and maximum magnitude recurrence models that are governed by asymptotic and finite return period properties, respectively. Integrated perceptibility curves illustrating the probability of a specific level of perceptible ground motion due to all earthquakes over the magnitude range extending from ?∞ to a magnitude M _i are then developed from reviewing site-specific magnitude perceptibility. These lead on to achieving site-specific annual probability of exceedance hazard curves for the example cities of Sofia and Thessaloniki for both horizontal ground acceleration and ground velocity. Both the maximum credible earthquake magnitude M ₃ and the most perceptible earthquake magnitude M _P(max) are of importance to the earthquake engineer when approaching anti-seismic building design. Both forms of hazard are illustrated using contoured hazard maps for the region bounded by 39°–45°N, 19°–29°E. Patterns are observed for these magnitude hazard estimates—especially M _P(max) specific to horizontal ground acceleration and horizontal ground velocity—and compared to inferred patterns of crustal deformation across the region. The full geographic region considered is estimated to be subject to a maximum credible earthquake magnitude M ₃—estimated using cumulative seismic moment release statistics—of 7.53 M _w, calculated from the full content of the adopted earthquake catalogue, while Bulgaria’s capital, Sofia, is estimated a comparable value of 7.36 M _w. Sofia is also forecast most perceptible earthquake magnitudes for the lowest levels considered for horizontal ground acceleration of M _PA(50) = 7.20 M _w and horizontal ground velocity of M _PV(5) = 7.23 M _w for a specimen focal depth of 15 km.     

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.     

Correlation Fractal and Multifractal Characteristics of Seismic Activity in the Taiwan Area, China

Based on the analysis of newly collected data of plate tectonics, distribution of active faults and crustal deformation, the Taiwan area is divided into two seismic regions and six seismic belts. Then, correlation fractal dimensions of all the regions and belts are calculated, and the fractal characteristics of hypocenteral distribution can be quantitatively analyzed. Finally, multifractal dimensions Dq and f(α) are calculated by using the earthquake catalog of the past 11 years in the Taiwan area. This study indicates that (1) there exists a favorable corresponding relationship between spatial images of seismic activity described with correlation fractal dimension analysis and tectonic settings; (2) the temporal structure of earthquakes is not single but multifractal fractal, and the pattern of Dq variation with time is a good indicator for predicting strong earthquake events.     

Global regression relations for conversion of surface wave and body wave magnitudes to moment magnitude

A homogenous earthquake catalog is a basic input for seismic hazard estimation, and other seismicity studies. The preparation of a homogenous earthquake catalog for a seismic region needs regressed relations for conversion of different magnitudes types, e.g. m _b , M _s , to the unified moment magnitude M _w. In case of small data sets for any seismic region, it is not possible to have reliable region specific conversion relations and alternatively appropriate global regression relations for the required magnitude ranges and focal depths can be utilized. In this study, we collected global events magnitude data from ISC, NEIC and GCMT databases for the period 1976 to May, 2007. Data for m_b magnitudes for 3,48,423 events for ISC and 2,38,525 events for NEIC, M _s magnitudes for 81,974 events from ISC and 16,019 events for NEIC along with 27,229 M _w events data from GCMT has been considered. An epicentral plot for M _w events considered in this study is also shown. M _s determinations by ISC and NEIC, have been verified to be equivalent. Orthogonal Standard Regression (OSR) relations have been obtained between M _s and M _w for focal depths (h < 70 km) in the magnitude ranges 3.0 ≤ M _s  ≤ 6.1 and 6.2 ≤ M _s  ≤ 8.4, and for focal depths 70 km ≤ h ≤ 643 km in the magnitude range 3.3 ≤ M _s  ≤ 7.2. Standard and Inverted Standard Regression plots are also shown along with OSR to ascertain the validation of orthogonal regression for M _s magnitudes. The OSR relations have smaller uncertainty compared to SR and ISR relations for M _s conversions. ISR relations between m _b and M _w have been obtained for magnitude ranges 2.9 ≤ m _b  ≤ 6.5, for ISC events and 3.8 ≤ m _b  ≤ 6.5 for NEIC events. The regression relations derived in this study based on global data are useful empirical relations to develop homogenous earthquake catalogs in the absence of regional regression relations, as the events catalog for most seismic regions are heterogeneous in magnitude types.     

Thermal infrared anomalies as a precursor of strong earthquakes in the distant future

Satellite thermal infrared images contain valuable earthquake precursor information. Past studies concluded that such information appeared only a few days or dozens of days before an earthquake would occur. In our study, though, we observed that the time intervals between the thermal infrared precursor and an earthquake??s occurrence can be up to 10?years. An infrared image can also synchronously indicate the locations of additional future earthquakes with different epicenters within a region. The shape, area, intensity, and movement of thermal infrared anomaly areas are a combination of all the future strong earthquakes within a region. These distant future earthquakes are generally located near the edges, endpoints, or corners of the main structure, fine structures or periphery structures of a thermal infrared anomaly area and play a role in confining the anomaly area. There have not been any exceptions among the strong earthquakes we analyzed, which have included the 2011 Japan M _w 9 event, the 2010 Yushu M _S 7.1 event, the 2008 Wenchuan M _S 8 event, and many other strong events following the 2004 Sumatra M _S 9 event. Surprisingly, some of the earthquakes can outline an area of elevated temperature observed many months ago. If we can roughly locate these potential epicenters through the analysis of thermal infrared images and combining the analysis with other information, and then dynamically monitor them, it may be easier to observe the precursor of an earthquake and predict its occurrence.     

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.