Coseismic Coulomb stress changes imparted by the 1996 Minahasa Mw7.9 earthquake on the 2018 Palu Mw7.5 earthquake and expected seismicity rate changes

The Mw7.5 Palu earthquake that occurred on September 28, 2018, in Indonesia caused much damage to the city of Palu. Preliminary investigations indicated that the Palu‐Koro Fault (PKF) hosted this damaging event. We calculated the seismicity before and after the 1996 Minahasa Mw7.9 earthquake and found that the seismicity on the PKF was enhanced after this earthquake. The earlier earthquake added Coulomb stress changes (?CFS) to the seismogenic fault plane. We calculated the ?CFS produced by the Palu earthquake on a specified received nodal plane; the results suggest that many aftershocks occurred in the region of increased ?CFS. This region was consistent with the region of increased seismicity. The ?CFS on neighbouring faults increased, and up to 55.282 bar of stress was observed on the PKF. Furthermore, we calculated the expected seismicity rate and found that it will require ~50 years to recover to its original level.     

Possible evidence of remotely triggered and delayed seismicity due to the 2001 Bhuj earthquake (Mw?=?7.6) in western India

Remote triggering by large earthquakes at regional distances is a globally observed phenomenon. However, there are no reports of observations of dynamic triggering at regional distances of several source lengths associated with the large Mw?=?7.6 Bhuj earthquake of January 26, 2001, in western India. In the present study, a swarm of over 140 microearthquakes that occurred about 500?km southeast of Bhuj, in the geothermal province of the Western Ghats in the Deccan volcanic province (DVP) of India, immediately after the occurrence of the Bhuj earthquake in 2001 is investigated. The post-Bhuj seismicity (M?<?2.0) occurred in three bursts spread over 2?months with each burst of intense activity lasting for 2?C3?days. All the three bursts of seismicity occurred in the same volume along a 5-km-long NW?CSE trending fault. The temporal coincidence and the sudden rise in seismicity that interrupts the characteristically low background seismicity strongly suggest that the Bhuj earthquake may have remotely triggered this activity. The triggered seismicity began approximately 2.5?h after the onset of the Bhuj mainshock and continued well after the passage of the surface waves, suggesting that the dynamic stresses possibly gave rise to secondary time-dependent mechanisms leading to the triggering. It is proposed that the triggered and delayed seismicity is possibly a consequence of the redistribution in pore fluid pressure due to the Bhuj earthquake. This is the first documented observation of remotely triggered seismicity at regional distances due to the Bhuj earthquake.     

Lineament fabric from airborne LiDAR and its influence on triggered earthquakes in the Koyna-Warna region,western India

We present the results of the first airborne LiDAR survey flown in the Koyna-Warna region and examine the relationship between the lineament fabric and the ongoing seismicity in the region. Our studies reveal that earthquakes of M≥4.0 for the period 1968 to 2016 are strongly correlated with a 10 km wide N-S fracture zone, which not only represents the surface expression of seismically active basement faults, but also act as conduits for water percolation between the Koyna and Warna reservoirs. A decreasing trend in the annual distribution of earthquakes was observed from 1985. A new burst of seismic activity in 1993 followed the impoudment of the Warna reservoir. We report a change in annual seismicity pattern, where seismicity peaks during September and December in the pre-Warna period, with a new peak emerging during March-April subsequent to the impoundment of Warna reservoir. A model is proposed to explain the seismicity along dominant N-S lineaments and the impact of Warna reservoir impounding which altered the hydrogeologic regime in the region.     

Triggered/migrated seismicity due to the 2001 M w7.7 Bhuj earthquake, Western India

Paper describes triggered seismicity to 200?km distance and for a decade due to the 2001 M _w7.7 Bhuj earthquake. The Kachchh region is seismically one of the most active intraplate regions of the World due to the occurrence of two large earthquakes 1819 (M _w7.8) and 2001 (M _w7.7). Though, it has high hazard but was known to have low seismicity in view of the occurrence of fewer smaller shocks. However, the status seems to have changed after 2001. Besides the strong aftershock activity for over a decade, seismicity has spread to nearby faults in Kachchh peninsula and at several places southward for 200?km distance in Saurashtra peninsula. Beyond the rupture zone of the 2001 Bhuj earthquake, more than 40 mainshocks of M _w?~?3?C5 have occurred at 20 different locations, which is unusual. The increased seismicity is inferred to be caused by stress perturbation due to the 2001 Bhuj earthquake by viscoelastic process. In Saurashtra, over and above the viscoelastic stress increase, the transient stress increase by water table rise in monsoons seems to be affecting the timing of mainshocks and associated sequences of earthquakes.     

Seismicity of the Cyprus region

We have used a combination of macroseismic and instrumental techniques to re-evaluate the seismicity of Cyprus and the surrounding region for the period 1890–1900. We identified 56 events which we believe to have a surface-wave magnitude (Ms) of 4.5 or greater. Early events are located by felt effects alone, and macroseismic information is used to improve the location of many later events. We find the broad pattern of seismicity follows that given by standard listings: an area of weak seismicity to the north-west of Cyprus in the Gulf of Antalia includes intermediate-depth events as are found in the Hellenic Arc to the west; a gap of reduced activity separates this activity from that of the main region of Cyprus, where we could find no evidence for sub-crustal foci; there appears to be no present seismological connection along the presumed plate boundary to the north-east to the junction of the East Anatolian fault and the extension of the Dead Sea rift system. We find a much greater level of activity in the first part of this century than is given in the standard listings, such as that of the ISC. Up to 1963 we identified 47 events of Ms > 4.5, compared with only 10 given in the ISC files. Recent improved station coverage locates more events in the region, but these have not been of large magnitude, and the moment release in recent decades has been appreciably less than earlier in the century. For the events in our list the total moment release is 2.7x 1019 Nm with 75% of this originating in the main Cyprus region. We reaffirm the danger of uncritical use of standard earthquake listings for tectonic interpretation and in particular for hazard evaluation in areas such as this, where neither the historical nor the instrumental seismicity record is complete.     

An observational test of the stress accumulation model based on seismicity preceding the 1992 Landers, CA earthquake

We test the Bowman and King [Bowman, D.D., King, G.C.P., 2001a, Accelerating seismicity and stress accumulation before large earthquakes. Geophys. Res. Lett., 28 (21), 4039–4042, Bowman, D.D., King, G.C.P., 2001b. Stress transfer and seismicity changes before large earthquakes. C. R. Acad. Sci. Paris, 333, 591–599] Stress Accumulation model by examining the evolution of seismicity rates prior to the 1992 Landers, California earthquake. The Stress Accumulation (SA) model was developed to explain observations of accelerating seismicity preceding large earthquakes. The model proposes that accelerating seismicity sequences result from the tectonic loading of large fault structures through aseismic slip in the elasto-plastic lower crust. This loading progressively increases the stress on smaller faults within a critical region around the main structure, thereby causing the observed acceleration of precursory activity. A secondary prediction of the SA model is that the precursory seismicity rates should increase first at the edges of the critical region, with the rates gradually rising over time at closer distances to the main fault. We test this prediction by examining year-long seismicity rates between 1960 and 2004, as a function of distance from the Landers rupture. To quantify the significance of trends in the seismicity rates, we auto-correlate the data, using a range of spatial and temporal lags. We find weak evidence for increased seismicity rates propagating towards the Landers rupture, but cannot conclusively distinguish these results from those obtained for a random earthquake catalog. However, we find a strong indication of periodicity in the rate fluctuations, as well as high correlation between activity 130–170 km from Landers and seismicity rates within 50 km of the Landers rupture temporally offset 1.5–2 years. The implications of this spatio–temporal correlation will be addressed in future studies.     

Variations of Seismicity in the Avachinsky Gulf (Kamchatka, Russia)

Variations of seismic mode in the region of the Avachinsky Gulf (Kamchatka, Russia) are considered. Observed anomalies (seismic quiescence, the ring seismicity, reduction of the slope of the earthquake recurrence diagram) provide a basis to consider this region as a place of strong earthquake preparation. The Kamchatka regional catalogues of earthquakes between 1962–1995 were used in the analysis. A reduced seismicity rate is observed during 10 years in an area of 150 km × 60 km in size. During the last five years, in the vicinity of the area considered, earthquakes with M > 5 occurred three times more often than the average over thirty years. It is interpreted as ring seismicity. The block of 220 km × 220~km in size, including the quiescence zone, is characterized by a continuous decrease of the recurrence diagram slope, which has reached a minimum value for the last 33 years in this region.     

Seismotectonics and seismicity of the Silakhor region, Iran

This paper deals with seismotectonic and seismicity of the Silakhor region that shows high seismic activity in western Iran. Silakhor is a vast plain with several villages and cities of Dorud and Borujerd and a small town of Chalanchulan that were destroyed and/or damaged many times by large earthquakes. This paper addresses the historical and instrumental earthquakes and their causative faults, seismotectonic provinces and seismotectonic zones of the region. Available seismic data were normalized by means of time normalization technique that resulted in the magnitude-frequency relation for the Silakhor area and estimation of the return period of earthquakes with different magnitudes. Some active faults in this region include the Dorud fault, the main Zagros thrust, the Galehhatam fault, the Sahneh fault and others. Among them, the Dorud fault is an earthquake fault and is the cause for most of the large and intermediate earthquakes in the region. The return period of large earthquakes with magnitudes greater than 7.0 (Ms) is very low, however, the occurrence of destructive earthquakes is greater in the region than in the neighboring provinces. The study proves the high seismicity of this zone and it is required to develop an accurate national plan for future building and reinforcement of the existing buildings in this region.     

Seismicity and source parameters of moderate earthquakes in Sikkim Himalaya

In this study, we accurately relocate 360 earthquakes in the Sikkim Himalaya through the application of the double-difference algorithm to 4?years of data accrued from a eleven-station broadband seismic network. The analysis brings out two major clusters of seismicity??one located in between the main central thrust (MCT) and the main boundary thrust (MBT) and the other in the northwest region of Sikkim that is site to the devastating Mw6.9 earthquake of September 18, 2011. Keeping in view the limitations imposed by the Nyquist frequency of our data (10?Hz), we select 9 moderate size earthquakes (5.3????Ml????4) for the estimation of source parameters. Analysis of shear wave spectra of these earthquakes yields seismic moments in the range of 7.95?×?10²¹ dyne-cm to 6.31?×?10²³ dyne-cm and corner frequencies in the range of 1.8?C6.25?Hz. Smaller seismic moments obtained in Sikkim when compared with the rest of the Himalaya vindicates the lower seismicity levels in the region. Interestingly, it is observed that most of the events having larger seismic moment occur between MBT and MCT lending credence to our observation that this is the most active portion of Sikkim Himalaya. The estimates of stress drop and source radius range from 48 to 389?bar and 0.225 to 0.781?km, respectively. Stress drops do not seem to correlate with the scalar seismic moments affirming the view that stress drop is independent over a wide moment range. While the continental collision scenario can be invoked as a reason to explain a predominance of low stress drops in the Himalayan region, those with relatively higher stress drops in Sikkim Himalaya could be attributed to their affinity with strike-slip source mechanisms. Least square regression of the scalar seismic moment (M ₀) and local magnitude (Ml) results in a relation LogM ₀?=?(1.56?±?0.05)Ml?+?(8.55?±?0.12) while that between moment magnitude (M _w ) and local magnitude as M _w ?=?(0.92?±?0.04)Ml?+?(0.14?±?0.06). These relations could serve as useful inputs for the assessment of earthquake hazard in this seismically active region of Himalaya.     

A cost–benefit analysis of mitigation options for optimal management of risks posed by flow-like phenomena

We analyzed the fractal and multifractal properties of the earthquake time series occurred around the Enguri dam in West Georgia by applying the methods of detrended fluctuation analysis and multifractal detrended fluctuation analysis. We examined the interevent time series in two periods: (1) 1960–1980, in which the investigated area was characterized by the natural seismicity; and (2) 1981–2012, in which the quasi-periodic change of the reservoir water level affected the earthquake generation. Our findings show that the water level variation may influence the fractal properties of earthquake temporal distribution in the local area around the Enguri dam. In particular, it is observed that the time distribution features of seismicity occurred in the second period are more persistent than the natural seismicity occurred in the first period. Furthermore, the seismic process of the second period shows a lower multifractal degree than that of the first period, indicating that the influence of quasi-periodic fluctuation of water level features the seismicity as more regular compared to the natural seismicity.     

Premonitory seismicity changes prior to the Parkfield and Coalinga earthquakes in southern California

In this work, we apply the Pattern Informatics technique for evaluating one surface expression of the underlying stress field, the seismicity, in order to study the Parkfield–Coalinga interaction over the years preceding the 1983 Coalinga earthquake. We find that significant anomalous seismicity changes occur during the mid-1970s in this region prior to the Coalinga earthquake that illustrate a reduction in the probability of an event at Parkfield, while the probability of an event at Coalinga is seen to increase. This suggests that the one event did not trigger or hinder the other, rather that the dynamics of the earthquake system are a function of stress field changes on a larger spatial and temporal scale.     

A review of the seismicity and seismotectonics of Delhi and adjoining areas

Understanding of seismicity and seismotectonics of Delhi and adjoining areas is essential as these areas lie in the seismic zone IV and are geologically confined to the Delhi Fold Belt (DFB), juxtaposed to the Himalayan Frontal Thrust Fold Belt. Owing to the set-up, seismicity in this area is ascribed to the Himalayan Thrust System and activation of DFB Fault Systems. Considerably improved instrumental seismic monitoring in this area and data analysis had resolved three regions of pronounced seismicity that lie close to Sonepat, Rohtak and western part of the NCT Delhi, attributed to activation of various portions of the fault systems of the DFB. Based on seismic telemetry network data, the seismicity pattern analysis revealed that the Mahendragarh Dehradun Sub-Surface Fault (MDSSF) and Delhi Sargodha Ridge (DSR) are the two major zones of structural importance for the nucleation of seismicity in this region. These revelations were corroborated with the fault plane solution of the earthquakes. The dominant mechanism in nucleation of seismicity in DFB is the thrust with minor strike slip. The seismicity and seismotectonics of Delhi and adjoining areas endemic to activation of DFB is reviewed and presented in this paper.     

A model of characteristic earthquakes and its implications for regional seismicity

Regional seismicity (i.e. that averaged over large enough areas over long enough periods of time) has a size–frequency relationship, the Gutenberg–Richter law, which differs from that found for some seismic faults, the Characteristic Earthquake relationship. But all seismicity comes in the end from active faults, so the question arises of how one seismicity pattern could emerge from the other. The recently introduced Minimalist Model of Vázquez‐Prada et al. of characteristic earthquakes provides a simple representation of the seismicity originating from a single fault. Here, we show that a Characteristic Earthquake relationship together with a fractal distribution of fault lengths can accurately describe the total seismicity produced in a region. The resulting earthquake catalogue accounts for the addition of both all the characteristic and all the non‐characteristic events triggered in the faults. The global accumulated size–frequency relationship strongly depends on the fault length fractal exponent and, for fractal exponents close to 2, correctly describes a Gutenberg–Richter distribution with a b exponent compatible with real seismicity.     

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.     

Seismicity and seismotectonics of the Gulf of Aqaba region

The historical seismicity of the last ten centuries and the instrumental data that occurred in the Gulf of Aqaba region during the period 1982–2008 are evaluated. It is found that 12 historical earthquakes have occurred with average recurrence periods 70–90 and 333–500 years for M?≥?6.0 and 7.0, respectively. Those with M?≤?6.5 appear to be incomplete and require further investigation. More than 98 % of the instrumental data has occurred in the form of swarms and sequences. The first have released about 32 % of the total energy and are most likely related to subsurface volcanic activities. Their epicentral distribution indicates that all regional faults of the gulf area are active in the present, but with clear concentration within the area bound by latitudes 28.2°–29.8° and longitudes 34.4°–35.2°. Regional strike-slip faults of the northern two basins appear to be as twice active as the normal, or more. An appreciable level of seismic hazard is envisaged as the “a” value is 6.0–6.2 while the “b” value shows a temporal variation, mostly in the range 0.8–1.05. More than 95 % of the seismic energy was released from earthquakes shallower than 22 km. This indicates a brittle upper crust and a ductile lower crust and upper mantle. Tectonic movements at depths?>?22 km appear to be aseismic. The epicentral distribution of the five swarms indicates that the lengths of the causative faults varied in the range 45–70 km. The maximum expected magnitude is Mw?=?6.8–7.2. This implies a seismic slip rate of about 0.54–0.8 Cm/year and some 20–30 % of aseismic tectonic movements. This and the sequence nature of the seismicity of this region result in a noticeable hazard reduction. Combining the seismicity data of the Gulf of Aqaba region with other geophysical, geological, tectonic, and environmental data, clearly indicate that the seismicity of this region is as old as the initiation of the gulf itself. No apparent southward or northward migration of activity is observed.     

A seismotectonic study of the Southeastern Alaska Region

We compare relocations of recent (1973–2005) and historic (1919–1972) earthquakes to geologic and geophysical (gravity, aeromagnetic, and uplift) information to determine the relationship of seismicity to crustal deformation in southeastern Alaska. Our results suggest that along strike changes in the structure of the Pacific plate may control the location of the ends of rupture zones for large earthquakes along the offshore Queen Charlotte fault system in the southern portion of the study area. There is a marked increase in background seismicity in the northern portion of the study area where the Fairweather fault begins to bend toward the northwest and crustal uplift due to glacial unloading exceeds 20 mm/year. Focal mechanisms indicate that thrust and reverse mechanisms predominate in the region of maximum uplift, as might be expected by the decrease in ice sheet thickness. The diffuse nature of seismicity between the Fairweather and Denali faults in the northern study area suggests a complex interaction between plate/microplate interactions and glacial unloading, making it difficult to determine the optimal fault orientation for failure in moderate magnitude (5.5 to 6.5) earthquakes within this region.     

A new synthesis of Pyrenean seismicity and its tectonic implications

The Pyrenean range, which results from the convergence of the Iberian and Eurasian plates along the North Pyrenean fault, exhibits a permanent seismic activity with moderate magnitude events. From the end of the 1980s, seismic instrumentation has been deployed in the Pyrenees, making now possible the computation of improved seismicity maps. We have gathered all the arrival times published for the period 1989–1996 by the different Spanish and French institutions in charge of the seismic survey of the range, and reprocessed them in an homogeneous way, in order to obtain a coherent seismicity map over the whole range. Particular attention has been paid to the evaluation of the quality of the locations and to the focal depth determinations. The comparison with previous maps of the Pyrenean seismicity reveals significant improvements in both the quality of locations and the threshold of detection. The new seismicity map reveals that the North Pyrenean fault is active only in the western part of the range. In the central and eastern parts, the seismicity involves other tectonic units such as the Maladeta and Canigou granitic massifs, the North Pyrenean Frontal Thrust, the Tet fault and the volcanic units in Catalonia. Despite the short time interval considered, this new seismicity file may be a valuable tool for future tectonic studies.     

Multifractal analysis of earthquakes in Kumaun Himalaya and its surrounding region

Himalayan seismicity is related to continuing northward convergence of Indian plate against Eurasian plate. Earthquakes in this region are mainly caused due to release of elastic strain energy. The Himalayan region can be attributed to highly complex geodynamic process and therefore is best suited for multifractal seismicity analysis. Fractal analysis of earthquakes (mb ?? 3.5) occurred during 1973?C2008 led to the detection of a clustering pattern in the narrow time span. This clustering was identified in three windows of 50 events each having low spatial correlation fractal dimension (D _C ) value 0.836, 0.946 and 0.285 which were mainly during the span of 1998 to 2005. This clustering may be considered as an indication of a highly stressed region. The Guttenberg Richter b-value was determined for the same subsets considered for the D _C estimation. Based on the fractal clustering pattern of events, we conclude that the clustered events are indicative of a highly stressed region of weak zone from where the rupture propagation eventually may nucleate as a strong earthquake. Multifractal analysis gave some understanding of the heterogeneity of fractal structure of the seismicity and existence of complex interconnected structure of the Himalayan thrust systems. The present analysis indicates an impending strong earthquake, which might help in better hazard mitigation for the Kumaun Himalaya and its surrounding region.     

Co-seismic gravity changes in the Koyna-Warna region: Implications of mass redistribution

Koyna-Warna Region (KWR) is one of the known sites for reservoir triggered seismicity. The continued triggered seismicity over the five decades is restricted to a region of about 600–700 sq. km, which provides a unique opportunity to monitor geophysical anomalies likely to be associated with seismicity of the region. Present study confers temporal gravity changes recorded by gPhone and GRACE satellite and interprets observed changes in conjunction with seismological, geodetic (cGPS) observations and groundwater level measurements. GRACE data suggest that seasonal vertical deformation due to hydrological loading is ～ 2 cm, which corroborates with continuous GPS observations. Seasonal hydrological loading of the region, which is in a phase of reservoir loading, might be influencing the critically stressed KWR leading to the seasonal seismicity of the region. The gPhone gravity data distinctly show co-seismic gravity signals for eight earthquakes of M_w > 2 and gravity anomalies show positive correlation on a logarithmic scale with earthquake released energy. To investigate the cause of gravity changes, an estimate is made for 14^th April 2012 earthquake for M_w 4.8 using fault dislocation model. The recorded gravity changes of 189 μGal by gPhone located at a distance of 28 km from the hypocentre is much more than the estimate of ～0.1 μGal calculated for M_w 4.8 Koyna earthquake. Therefore, it is inferred that co-seismic gravity signals for eight earthquakes are primarily caused due to redistribution of mass at shallow depth.     

Annular seismicity structures and the March 11, 2011, earthquake (Mw = 9.0) in Northeast Japan

The characteristics of seismicity prior to the series of eight very strong earthquakes (M _w = 7.0–9.0) in Northeast Japan are discussed. Ring seismicity structures that appeared prior to all eight events in two depth ranges of 0–33 and 34–70 km are identified. The epicenters of the main shocks were located near areas of crossing or touching of shallow and deep rings. It was shown that the sizes of shallow rings and threshold magnitudes corresponding to seismicity rings grow with the energy of the main shocks. It was noted that the prognosis with respect to the place and magnitude of the catastrophic earthquake on March 11, 2011, had been made before it based on the data obtained prior to July 1, 2009. Use of the new data obtained prior to March 10, 2011, enabled us to specify this prognosis significantly. We obtained correlation dependences of threshold magnitudes on the energy of the main shocks (with a high correlation coefficients). It was shown that the duration of the period for seismicity rings to emerge in the considered region nearly did not depend on magnitude. The nature of annular structures and the possibility of application of their parameters for prognosis of strong earthquakes were discussed.