Coseismic responses and the mechanism behind MW 5.1 earthquake of March 14, 2005 in the Koyna–Warna region, India

An earthquake of M_w 5.1 occurred on March 14, 2005, in the seismically active Koyna–Warna region in western India, the site known for the largest reservoir triggered seismicity (RTS) in the world. For more than four decades, earthquakes with M  4.0 have occurred in this region at regular intervals. Impoundment of reservoirs and changes in lake levels can trigger earthquakes by two processes of stress modifications, namely direct loading effect of the reservoir and diffusion through various faults and fractures. In this paper we analysed the reservoir water level data at Koyna and Warna reservoirs prior to the occurrence of the March 14, 2005 earthquake, to explain the dominant mechanism behind its occurrence and its correlation with the observed coseismic changes. We conclude that the diffusion process, not the reservoir load effect, is the dominating mechanism triggering earthquakes in the region. The coseismic changes in deep well water levels sensitive to earth tides are found to be to the order of 1–12 cm.     

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.     

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.

     

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.     

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.     

Estimating tsunami potential of earthquakes in the Sumatra–Andaman region based on broadband seismograms in India

In this paper, we report that the ratio of broadband energy (0.01?C2?Hz) to high-frequency energy (0.3?C2?Hz), E _r, estimated from regional seismograms of India, might be a useful parameter in estimating tsunami potential of earthquakes in the Sumatra?CAndaman region. E _r is expected to be sensitive to the depth as well as to the source characteristics of an earthquake. Since a shallow and slow earthquake has a greater tsunamigenic potential, E _r may be a useful diagnostic parameter. We base our analysis on broadband seismograms of the great earthquakes of Sumatra?CAndaman (2004, M _w?~?9.2) and Nias (2005, M _w 8.6), 41 of their aftershocks, and the earthquakes of north Sumatra (2010, M _w 7.8) and Nicobar (2010, M _w 7.4) recorded at VISK, a station located on the east coast of India. In the analysis, we also included the two recent, great strike-slip earthquakes of north Sumatra (2012, M _w 8.6, 8.2) recorded at VISK and three south Sumatra earthquakes (2007, M _w 8.5; 2007, M _w 7.9; 2010, M _w 7.8) recorded at PALK, a station in Sri Lanka. We find that E _r is a function of depth; shallower earthquakes have higher E _r values than the deeper ones. Thus, E _r may be indicative of tsunamigenic potential of an earthquake. As M _w and E _r increase so does the tsunami potential. In addition to the parameter E _r, the radiated seismic energy, E _s, may be estimated from the regional seismograms in India using empirical Green??s function technique. The technique yields reliable E _s for the great Sumatra and Nias earthquakes. E _r and E _s computed from VISK data, along with M _w and focal mechanism, may be useful in estimating tsunami potential along the east coast of India from earthquakes in the Sumatra?CAndaman region in less than ~20?min.     

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.     

Multifractality in seismic sequences of NW Himalaya

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.

     

Quantifying the role of sandy–silty sediments in generating slope failures during earthquakes: example from the Algerian margin

The Algerian margin is a seismically active region, where during the last century, several large magnitude earthquakes took place. This study combines geotechnical and sedimentological data with numerical modelling to quantitatively assess the present-day slope stability of the Algerian margin. Geotechnical laboratory tests, such as cyclic triaxial tests, oedometric tests and vane shear tests were carried out on sediment cores collected on the study area. The liquefaction potential of a sediment column located about 30 km from the Boumerdès earthquake epicentre of 21st May 2003 was evaluated theoretically for an earthquake of M _w  = 6.8. We show that thin sand and silt beds such as those described on recovered sediment cores are the main cause of sediment deformation and liquefaction during earthquakes. Numerical calculations showed that the slope failure may occur during an earthquake characterised by a PGA in excess of 0.1g, and also that, under a PGA of 0.2g liquefaction could be triggered in shallow silty–sandy deposits. Moreover, comparison of the predicted slope failure with failure geometries inferred from seafloor morphology showed that earthquakes and subsequent mass movements could explain the present-day morphology of the study area.     

Spatial characteristics of landslides triggered by the 2015 M<Subscript>w</Subscript> 7.8 (Gorkha) and M<Subscript>w</Subscript> 7.3 (Dolakha) earthquakes in Nepal

High magnitude earthquakes trigger numerous landslides and their occurrences are mainly controlled by terrain parameters. We created an inventory of 15,551 landslides with a total area of 90.2 km² triggered by the 2015 M_w 7.8 (Gorkha) and M_w 7.3 (Dolakha) earthquakes in Nepal, through interpretation of very high resolution satellite images (e.g. WorldView, Pleiades, Cartosat-1 and 2, Resourcesat-2). Our spatial analysis of landslide occurrences with ground acceleration, slope, lithology and surface defomation indicated ubiquitous control of steep slope on landslides with ground acceleration as the trigger. Spatial distribution of landslides shows increasing frequency away from the Gorkha earthquake epicentre up to 130 km towards east, dropping sharply thereafter, which is an abnormal phenomenon of coseismic landslides. Landslides are laterally concentrated in three zones which matches well with the seismic rupture evolution of Gorkha earthquake, as reported through teleseismic measurements.     

Paleoearthquakes in the Uimon basin (Gorny Altai)

Paleoseismological studies confirm that the Uimon basin is thrust by its northern mountain border along the active South Terekta fault. The latest motion along the fault in the 7-8th centuries AD induced an earthquake with a magnitude of M_w= 7.4-7.7 and a shaking intensity of I = 9-11 on the MSK-64 scale. The same fault generated another event (M > 7, I = 9-10), possibly, about 16 kyr ago, which triggered gravity sliding. The rockslide dammed the Uimon valley and produced a lake, where lacustrine deposition began about 14 ± 1 kyr ago, and a later M > 7 (I = 9-10) earthquake at ~ 6 ka caused the dam collapse and the lake drainage. Traces of much older earthquakes that occurred within the Uimon basin are detectable from secondary deformation structures (seismites) in soft sediments deposited during the drainage of a Late Pleistocene ice-dammed lake between 100 and 90 ka and in ~ 77 ka alluvium. The magnitude and intensity of these paleoearthquakes were at least M > 5.0-5.5 and I > 6-7.     

Features of seismic activation of the Middle Baikal region, 2008–2011

Scenarios for developing focal zones of strong (M _w = 5.3) earthquakes that occurred in the Middle Baikal region in 2008 and 2011 are considered. The new (submeridional and sublatitudinal) lines of destruction of the Earth’s crust in the water area of the lake are recorded. The facts of seismoactive structures forming in the surrounding mountains (to the southeast) under typical rift conditions of movements are established, which indicates that the basin of Lake Baikal is expanding and growing due to active capturing and processing of its mountain surroundings.     

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 earthquakes and related tsunamis of October 6, 1944 and March 7, 1867; NE Aegean Sea

During two distinct earthquakes occurred on March 7, 1867 and October 6, 1944, tsunami waves were also observed at some localities around the Gulf of Edremit, NE Aegean Sea. The first event (M _w = 6.8) mostly affected the city of Mitilini of Lesvos Island while the Gulf of Edremit-Ayvacık earthquake (M _S = 6.8) largely affected the northern and eastern coastal areas of the Gulf of Edremit. In 1944 earthquake, numerous surface cracks and water gushes were reported. The coastal neighborhoods of the town of Ayvalık in the east were flooded by tsunami waves. At the WSW extend of the main fault observed on land, which is parallel to the present-day slip vectors, some normal-oblique faults were observed close and subparallel to the northern coast. On the basis of historical documents, reports, interviews, geological setting, field observations and marine seismic reflection data, the 1944 earthquake was not triggered by one of the main fault segments but by a secondary fault or fault group which was described in this study. Depending on the distribution of tensional and compressional forces in the region, which rotates clockwise under the control of the middle strand of the North Anatolian fault, secondary fault groups become important. The moment tensor parameters of such small-size events have been determined and have obtained consistent results with the faults proposed in this study.     

Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Evaluation of earthquake-induced strain in promoting mud eruptions: the case of Shamakhi–Gobustan–Absheron areas,Azerbaijan

Although a relationship between the occurrence of large earthquakes and the eruptions of close mud volcanoes is well known, several uncertainties remain on understanding the triggering mechanisms. In the present study, we evaluate both the static and dynamic strains induced by earthquakes in the substratum of mud volcanoes. We studied the effects of two earthquakes of M _w 6.18 and 6.08 occurred in the Caspian Sea on 25 November 2000 close to Baku city, Azerbaijan. A total of 33 eruptions occurred at 24 mud volcanoes within a maximum distance of 108 km from the epicentres in the 5 years following the earthquakes. The overall eruption rate in the studied area of the 50 years before the 2000 earthquakes was 1.24 that is much smaller than the eruption rate of 6.6 of the 5 years following these earthquakes. The largest number of eruptions occurred within 2 years from the earthquakes with the highest frequency within 6 months. Our calculated earthquake-induced static effects show that crustal dilatation might have triggered only seven eruptions at a maximum distance of about 60 km from the epicentres and within 3 years. Based on our data, dynamic rather than static strain is likely to have been the dominating “promoting” factor because it affected all the studied unrest volcanoes and its magnitude was much larger.     

A comparison of the sedimentary records of the 1960 and 2010 great Chilean earthquakes in 17 lakes: Implications for quantitative lacustrine palaeoseismology

Seismically‐induced event deposits embedded in the sedimentary infill of lacustrine basins are highly useful for palaeoseismic reconstructions. Recent, well‐documented, great megathrust earthquakes provide an ideal opportunity to calibrate seismically‐induced event deposits for lakes with different characteristics and located in different settings. This study used 107 short sediment cores to investigate the sedimentary impact of the 1960 M_w 9·5 Valdivia and the 2010 M_w 8·8 Maule earthquakes in 17 lakes in South‐Central Chile (i.e. lakes Negra, Lo Encañado, Aculeo, Vichuquén, Laja, Villarrica, Calafquén, Pullinque, Pellaifa, Panguipulli, Neltume, Riñihue, Ranco, Maihue, Puyehue, Rupanco and Llanquihue). A combination of image analysis, magnetic susceptibility and grain‐size analysis allows identification of five types of seismically‐induced event deposits: (i) mass‐transport deposits; (ii) in situ deformations; (iii) lacustrine turbidites with a composition similar to the hemipelagic background sediments (lacustrine turbidites type 1); (iv) lacustrine turbidites with a composition different from the background sediments (lacustrine turbidites type 2) and (v) megaturbidites. These seismically‐induced event deposits were compared to local seismic intensities of the causative earthquakes, eyewitness reports, post‐earthquake observations, and vegetation and geomorphology of the catchment and the lake. Megaturbidites occur where lake seiches took place. Lacustrine turbidites type 2 can be the result of: (i) local near‐shore mass wasting; (ii) delta collapse; (iii) onshore landslides; (iv) debris flows or mudflows; or (v) fluvial reworking of landslide debris. On the contrary, lacustrine turbidites type 1 are the result of shallow mass wasting on sublacustrine slopes covered by hemipelagic sediments. Due to their more constrained origin, lacustrine turbidites type 1 are the most reliable type of seismically‐induced event deposits in quantitative palaeoseismology, because they are almost exclusively triggered by earthquake shaking. Moreover, they most sensitively record varying seismic shaking intensities. The number of lacustrine turbidites type 1 linearly increases with increasing seismic intensity, starting with no lacustrine turbidites type 1 at intensities between V½ and VI and reaching 100% when intensities are higher than VII½. Combining different types of seismically‐induced event deposits allows the reconstruction of the complete impact of an earthquake.     

On generation conditions of strong earthquakes in southern Baikal

In this paper the features of seismic process in the southern depression of Lake Baikal are considered. By the data on focal mechanisms of the earthquakes of February 25, 1999 (M _w = 6.0), and August 27, 2008 (M _w = 6.3), as well as based on configuration of their aftershock fields, it is determined that foci of strong seismic events in southern Baikal are controlled by the greatest structural elements of sublatitudinal and submeridional strikes. It has been shown that a substantial role in the formation of focal zones is played by low-scale destruction of the Earth’s crust, revealed by geological-geophysical data and proved by clustering of seismic shocks. New data on the August 27, 2008, earthquake have proved the high level of seismic danger of this part of the Baikal Rift Zone and allowed us to determine generation conditions of strong earthquakes more precisely.     

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.