Preliminary estimation of the tsunami hazards associated with the Makran subduction zone at the northwestern Indian Ocean

We present a preliminary estimation of tsunami hazard associated with the Makran subduction zone (MSZ) at the northwestern Indian Ocean. Makran is one of the two main tsunamigenic zones in the Indian Ocean, which has produced some tsunamis in the past. Northwestern Indian Ocean remains one of the least studied regions in the world in terms of tsunami hazard assessment. Hence, a scenario-based method is employed to provide an estimation of tsunami hazard in this region for the first time. The numerical modeling of tsunami is verified using historical observations of the 1945 Makran tsunami. Then, a number of tsunamis each resulting from a 1945-type earthquake (M _w 8.1) and spaced evenly along the MSZ are simulated. The results indicate that by moving a 1945-type earthquake along the MSZ, the southern coasts of Iran and Pakistan will experience the largest waves with heights of between 5 and 7 m, depending on the location of the source. The tsunami will reach a height of about 5 m and 2 m in northern coast of Oman and eastern coast of the United Arab Emirates, respectively.     

Aftershock sequences of two great Sumatran earthquakes of 2004 and 2005 and simulation of the minor tsunami generated on September 12, 2007 in the Indian Ocean and its effect

We present the seismic energy, strain energy, frequency–magnitude relation (b-value) and decay rate of aftershocks (p-value) for the aftershock sequences of the Andaman–Sumatra earthquakes of December 26, 2004 (M _w 9.3) and March 28, 2005 (M _w 8.7). The energy released in aftershocks of 2004 and 2005 earthquake was 0.135 and 0.365% of the energy of the respective mainshocks, while the strain release in aftershocks was 39 and 71% for the two earthquakes, respectively. The b-value and p-value indicate normal value of about 1. All these parameters are in normal range and indicate normal stress patterns and mechanical properties of the medium. Only the strain release in aftershocks was considerable. The fourth largest earthquake in this region since 2004 occurred in September 2007 off the southern coast of Island of Sumatra, generating a relatively minor tsunami as indicated by sea level gauges. The maximum wave amplitude as registered by the Padang, tide gauge, north of the earthquake epicenter was about 60 cm. TUNAMI-N2 model was used to investigate ability of the model to capture the minor tsunami and its effect on the eastern Indian Coast. A close comparison of the observed and simulated tsunami generation, propagation and wave height at tide gauge locations showed that the model was able to capture the minor tsunami phases. The directivity map shows that the maximum tsunami energy was in the southwest direction from the strike of the fault. Since the path of the tsunami for Indian coastlines is oblique, there were no impacts along the Indian coastlines except near the coast of epicentral 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.

     

Probabilistic tsunami hazard assessment along Oman coast from submarine earthquakes in the Makran subduction zone

The Sultanate of Oman is among the Indian Ocean countries that were subjected to at least two confirmed tsunamis during the twentieth and twenty-first centuries: the 1945 tsunami due to an earthquake in the Makran subduction zone in the Sea of Oman (near-regional field tsunami) and the Indian Ocean tsunami in 2004, caused by an earthquake from the Andaman Sumatra subduction zone (far - field tsunami). In this paper, we present a probabilistic tsunami hazard assessment for the entire coast of Oman from tectonic sources generated along the Makran subduction zone. The tsunami hazard is assessed taking into account the contribution of small- and large-event magnitudes. Results of the earthquake recurrence rate studies and the tsunami numerical modeling for different magnitudes were used through a logic-tree to estimate the tsunami hazard probabilities. We derive probability hazard exceedance maps for the Omani coast considering the exposure times of 100, 250, 500, and 1000 years. The hazard maps consist of computing the likelihood that tsunami waves exceed a specific amplitude. We find that the probability that a maximum wave amplitude exceeds 1 m somewhere along the coast of Oman reaches, respectively, 0.7 and 0.85 for 100 and 250 exposure times, and it is up to 1 for 500 and 1000 years of exposure times. These probability values decrease significantly toward the southern coast of Oman where the tsunami impact, from the earthquakes generated at Makran subduction zone, is low.     

Simulation of the Arabian Sea Tsunami propagation generated due to 1945 Makran Earthquake and its effect on western parts of Gujarat (India)

The 1945 Tsunami generated due to Makran Earthquake in the Arabian Sea was the most devastating tsunami in the history of the Arabian Sea and caused severe damage to property and loss of life. It occurred on 28th November 1945, 21:56 UTC (03:26 IST) with a magnitude of 8.0 (M _w), originating off the Makran Coast of Pakistan in the Arabian Sea. It has impacted as far as Mumbai in India and was noticed up to Karvar Coast, Karnataka. More than 4,000 people were killed as a result of the earthquake and the tsunami. In this paper an attempt is made for a numerical simulation of the tsunami generation from the source, its propagation into the Arabian Sea and its effect on the western coast of India through the use of a numerical model, referred to as Tunami-N2. The present simulation is carried out for a duration of 300 min. It is observed from the results that the simulated arrival time of tsunami waves at the western coast of India is in good agreement with the available data sources. The paper also presents run-up elevation maps prepared using Shuttle Radar Topographic Mission (SRTM) data, showing the possible area of inundation due to various wave heights along different parts of the Gujarat Coast. Thus, these results will be useful in planning the protection measures against inundation due to tsunami and in the implementation of a warning system.     

Insights on the small tsunami from January 28, 2020, Caribbean Sea MW7.7 earthquake by numerical simulation and spectral analysis

A major left-lateral strike-slip Mw7.7 earthquake occurred in the vicinity of the Caribbean Sea on January 28, 2020. As a result, a small-scale tsunami was generated. The properties of the seismogenic source were described using observational data gathered for the earthquake and tsunami, as well as information on the regional tectonic setting. The tsunami was simulated with the COMCOT model and Okada’s dislocation model from finite fault solutions for M_W7.7 Caribbean Sea earthquakes published by the United States Geological Survey. The simulation results were compared to tide gauge records to validate whether the seafloor’s vertical displacement generated by the strike-slip fault caused a small-scale tsunami. We conducted a spectral analysis of the tsunami to better understand the characteristics of tsunami records. The tsunami simulation results showed that the co-seismic vertical displacement caused by a strike-slip M_W7.7 earthquake could have contributed to the small-scale tsunami, but the anomalously large high-frequency tsunami waves recorded by the George Town tide gauge 11 min after the earthquake were unrelated to the earthquake-generated tsunami. According to the spectrum analysis, the predominant period of noticeable high-frequency tsunami waves recorded by the George Town tide gauge occurred only two minutes after the earthquake. This indicates that the source of the small-scale tsunami was close to the George Town station and the possible tsunami source was 150 km away from George Town station. These facts suggest that a submarine landslide was caused by the strike-slip earthquake. The comprehensive analysis showed that the small-scale tsunami was not caused solely by co-seismic seafloor deformation from the strike-slip event but that an earthquake-triggered submarine landslide was the primary cause. Therefore, the combined impact of two sources led to the small-scale tsunami.

     

Discovery of possible mega-thrust earthquake along the Seram Trough from records of 1629 tsunami in eastern Indonesian region

Arthur Wichmann’s “Earthquakes of the Indian Archipelago” documents several large earthquakes and tsunami throughout the Banda Arc region that can be interpreted as mega-thrust events. However, the source regions of these events are not known. One of the largest and well-documented events in the catalog is the great earthquake and tsunami affecting the Banda Islands on August 1, 1629. It caused severe damage from a 15-m tsunami that arrived at the Banda Islands about a half hour after violent shaking stopped. The earthquake was also recorded 230 km away in Ambon, but no tsunami is mentioned. This event was followed by at least 9 years of uncommonly frequent seismic activity in the region that tapered off with time, which can be interpreted as aftershocks. The combination of these observations indicates that the earthquake was most likely a mega-thrust event. We use an inverse modeling approach to numerically reconstruct the tsunami, which constrains the likely location and magnitude of the 1629 earthquake. Only, linear numerical models are applied due to the low resolution of bathymetry in the Banda Islands and Ambon. Therefore, we apply various wave amplification factors (1.5–4) derived from simulations of recent, well-constrained tsunami to bracket the upper and lower limits of earthquake moment magnitudes for the event. The closest major earthquake sources to the Banda Islands are the Tanimbar and Seram Troughs of the Banda subduction/collision zone. Other source regions are too far away for such a short arrival time of the tsunami after shaking. Moment magnitudes predicted by the models in order to produce a 15-m tsunami are M_w of 9.8–9.2 on the Tanimbar Trough and M_w 8.8–8.2 on the Seram Trough. The arrival times of these waves are 58 min for Tanimbar Trough and 30 min for Seram Trough. The model also predicts 5-m run-up for Ambon from a Tanimbar Trough source, which is inconsistent with the historical records. Ambon is mostly shielded from a wave generated by a Seram Trough source. We conclude that the most likely source of the 1629 mega-thrust earthquake is the Seram Trough. Only one earthquake >M_w 8.0 is recorded instrumentally from the eastern Indonesia region although high rates of strain (50–80 mm/a) are measured across the Seram section of the Banda subduction zone. Enough strain has already accumulated since the last major historical event to produce an earthquake of similar size to the 1629 event. Due to the rapid population growth in coastal areas in this region, it is imperative that the most vulnerable coastal areas prepare accordingly.     

Submarine landslides at the eastern Sunda margin: observations and tsunami impact assessment

Our analysis of new bathymetric data reveals six submarine landslides at the eastern Sunda margin between central Java and Sumba Island, Indonesia. Their volumes range between 1 km3 in the Java fore-arc basin up to 20 km3 at the trench off Sumba and Sumbawa. We estimate the potential hazard of each event by modeling the corresponding tsunami and its run-up on nearby coasts. Four slides are situated remarkably close to the epicenter of the 1977 tsunamigenic Sumba M _w  = 8.3 earthquake. However, comparison of documented tsunami run-up heights and arrival times with our modeling results neither allows us to confirm nor can we falsify the hypothesis that the earthquake triggered these submarine landslides.     

Reconstructing hydrodynamic flow parameters of the 1700 tsunami at Cannon Beach, Oregon, USA

Coastal communities in the western United States face risks of inundation by distant tsunamis that propagate across the Pacific Ocean as well as local tsunamis produced by great (M_w?>?8) earthquakes on the Cascadia subduction zone. In 1964, the M_w 9.2 Alaska earthquake launched a Pacific-wide tsunami that flooded Cannon Beach, a small community (population 1640) in northwestern Oregon, causing over $230,000 in damages. However, since the giant 2004 Indian Ocean tsunami, the 2010 Chile tsunami and the recent 2011 Tohoku-Oki tsunami, renewed concern over potential impacts of a Cascadia tsunami on the western US has motivated closer examination of the local hazard. This study applies a simple sediment transport model to reconstruct the flow speed of the most recent Cascadia tsunami that flooded the region in 1700 using the thickness and grain size of sand layers deposited by the waves. Sedimentary properties of sand from the 1700 tsunami deposit provide model inputs. The sediment transport model calculates tsunami flow speed from the shear velocity required to suspend the quantity and grain size distribution of the observed sand layers. The model assumes a steady, spatially uniform tsunami flow and that sand settles out of suspension forming a deposit when the flow velocity decreases to zero. Using flow depths constrained by numerical tsunami simulations for Cannon Beach, the sediment transport model calculated flow speeds of 6.5?C7.6?m/s for sites within 0.6?km of the beach and higher flow speeds (~8.8?m/s) for sites 0.8?C1.2?km inland. Flow speed calculated for sites within 0.6?km of the beach compare well with maximum velocities estimated for the largest tsunami simulation. The higher flow speeds calculated for the two sites furthest landward contrast with much lower maximum velocities (<3.8?m/s) predicted by numerical simulations. Grain size distributions of sand layers from the most distal sites are inconsistent with deposition from sediment falling out of suspension. We infer that rapid deceleration in tsunami flow and convergences in sediment transport formed unusually thick deposits. Consequently, higher flow speeds calculated by the sediment model probably overestimate the actual wave speed at sites furthest inland.     

Structural elements of the Makran region, Oman sea and their potential relevance to tsunamigenisis

The character of convergence along the Arabian–Iranian plate boundary changes radically eastward from the Zagros ranges to the Makran region. This appears to be due to collision of continental crust in the west, in contrast to subduction of oceanic crust in the east. The Makran subduction zone with a length of about 900 km display progressively older and highly deformed sedimentary units northward from the coast, together with an increase in elevation of the ranges. North of the Makran ranges are large subsiding basins, flanked to the north by active volcanoes. Based on 2D seismic reflection data obtained in this study, the main structural provinces and elements in the Gulf of Oman include: (i) the structural elements on the northeastern part of the Arabian Plate and, (ii) the Offshore Makran Accretionary Complex. Based on detailed analysis of these data on the northeastern part of the Arabian Plate five structural provinces and elements—the Musendam High, the Musendam Peneplain, the Musendam Slope, the Dibba Zone, and the Abyssal Plain have been identified. Further, the Offshore Makran Accretionary Complex shown is to consist Accretionary Prism and the For-Arc Basin, while the Accretionary Prism has been subdivided into the Accretionary Wedge and the Accreted/Colored Mélange. Lastly, it is important to note that the Makran subduction zone lacks the trench. The identification of these structural elements should help in better understanding the seismicity of the Makran region in general and the subduction zone in particular. The 1945 magnitude 8.1 tsunamigenic earthquake of the Makran and some other historical events are illustrative of the coastal region’s vulnerability to future tsunami in the area, and such data should be of value to the developing Indian Ocean Tsunami Warning System.     

Tsunami impact to Washington and northern Oregon from segment ruptures on the southern Cascadia subduction zone

This paper explores the size and arrival of tsunamis in Oregon and Washington from the most likely partial ruptures of the Cascadia subduction zone (CSZ) in order to determine (1) how quickly tsunami height declines away from sources, (2) evacuation time before significant inundation, and (3) extent of felt shaking that would trigger evacuation. According to interpretations of offshore turbidite deposits, the most frequent partial ruptures are of the southern CSZ. Combined recurrence of ruptures extending ~490 km from Cape Mendocino, California, to Waldport, Oregon (segment C) and ~320 km from Cape Mendocino to Cape Blanco, Oregon (segment D), is ~530 years. This recurrence is similar to frequency of full-margin ruptures on the CSZ inferred from paleoseismic data and to frequency of the largest distant tsunami sources threatening Washington and Oregon, ~M _w 9.2 earthquakes from the Gulf of Alaska. Simulated segment C and D ruptures produce relatively low-amplitude tsunamis north of source areas, even for extreme (20 m) peak slip on segment C. More than ~70 km north of segments C and D, the first tsunami arrival at the 10-m water depth has an amplitude of <1.9 m. The largest waves are trapped edge waves with amplitude ≤4.2 m that arrive ≥2 h after the earthquake. MM V–VI shaking could trigger evacuation of educated populaces as far north as Newport, Oregon for segment D events and Grays Harbor, Washington for segment C events. The NOAA and local warning systems will be the only warning at greater distances from sources.     

Rupture length and velocity for earthquakes in the Mid-Atlantic Ridge from directivity effect in body and surface waves

The dimensions and rupture velocities of four earthquakes, two in the Mid-Atlantic Ridge and two in Iceland with strike–slip mechanisms and magnitudes (M_w) between 6.2 and 6.8 were studied using the directivity effects of Rayleigh and body waves. For Rayleigh waves we used the directivity function for different pairs of stations and for body waves the waveforms of P and SH waves corresponding to a simple extended line source. We have found that three have very shallow depths about 3 km and one 8 km, fault lengths between 12 km and 21 km, and a low rupture velocity of about 1.5 km/s to 2.0 km/s which supports the idea of the presence of slow earthquakes in transform faults.     

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.     

大南海地区新生代板块构造活动

在新生代澳大利亚板块和欧亚板块之间的大洋中，存在一些地块(微板块)；同时，澳大利亚板块北部边缘的一些地块先后和澳大利亚板块分离，向北运动，与一些和欧亚板块分离出来的地块先后发生碰撞缝合。在此期间，由于地块分离而发生海底扩张，产生许多小洋盆，如南海、苏录海、苏拉威西海、安达曼海等，最后形成了东南亚地区今日的构造景观。笔者从大南海地区新生代的构造演化史之框架来研究南海地区新生代的构造演化历史，认为南海地区新生代的构造活动既与印度板块和欧亚板块的碰撞有关，也与太平洋板块向欧亚板块的俯冲活动有联系；同时，还受到澳大利亚板块向北运动之影响。南海地区在新生代发生过两次海底扩张，第一次海底扩张发生在42～35Ma前．是受印度板块和欧亚板块碰撞而引起欧亚大陆之下向东南方向之地幔流的影响而发生的，其海底扩张方向为NWSE，产生了南海西南海盆；第二次海底扩张发生于32～17Ma前。由于太平洋板块向欧亚板块俯冲，俯冲的大洋岩石圈已达700km深处，阻挡了欧亚大陆的上地幔向东南方向之流动，从而转向南流动。引起南海地区南北向海底扩张，即新生代第二次海底扩张，产生了南海中央海盆。南海新生代洋盆诞生之后，由于大南海地区继续有地块碰撞和边缘海海底扩张，对南海南部地区产生挤压，从而使这里的沉积发生变形，这就引起万安运动(南海南部)。     

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.     

Modeling the 2004 Sumatra tsunami at Seychelles Islands: site-effect analysis and comparison with observations

The M _w = 9.1 mega-thrust Sumatra–Andaman earthquake that occurred on December 26, 2004, was followed by a devastating tsunami that killed hundreds of thousands of people and caused catastrophic effects on human settlements and environments along many coasts of the Indian Ocean, where even countries very far from the source were affected. One of these cases is represented by the Republic of Seychelles, where the tsunami reached the region about 7 h after the earthquake and produced relevant damages, despite the country was more than 4,500 km far from the seismic source. In the present work, we present and discuss a study of the 2004 Sumatra tsunami by means of numerical simulations with the attention focused on the effects observed at the Seychelles Archipelago, a region never previously investigated with this approach. The case is interesting since these islands lay on a very shallow oceanic platform with steep slopes so as the ocean depth changes from thousands to few tens of meters over short distances, with significant effects on the tsunami propagation features: the waves are strongly refracted by the oceanic platform and the tsunami signal is modified by the introduction of additional frequencies. The study is used also to validate the UBO-TSUFD numerical code on a real tsunami event in the far field, and the results are compared with the available observations, i.e., the sea level time series recorded at the Pointe La Rue station, Mahé Island, and run-up measurements and inundation lines surveyed few weeks after the tsunami at Praslin Island, where the tsunami hit during low tide. Synthetic results are found in good agreement with observations, even though some of the observations remain not fully solved. Moreover, simulations have been run in high-tide condition since the 2004 Sumatra tsunami hitting at high tide can be taken as the worst-case scenario for the Seychelles islands and used for tsunami hazard and risk assessments.     

New discoveries in dynamics of an M8 earthquake-phenomena and their implications from the 2003 Tokachi-oki earthquake using a long term monitoring cabled observatory

At the 2003 Tokachi-oki earthquake of M8, seafloor phenomena such as a generation process of tsunami, seafloor uplifts, turbidity current, etc., were observed using a cabled observatory installed on the seafloor. The turbidity current was observed as a benthic storm caused presumably by the mainshock. The seafloor uplifts were observed at the mainshock and continuously after the mainshock. The uplifts were 0.35, 0.37, and 0.12 m for epicentral distances of 25.5, 31.4, and 81.7 km, respectively. After the mainshock, a continuous uplift of the seafloor is observed at all three pressure gauge locations indicating that there was a change in the state of friction on the plate boundary interface by the mainshock. In this paper, we first show what was observed using the cabled observatory installed right above the focal area of the earthquake, and then we discuss to summarize these phenomena associated with the earthquake, its possible causes, and future directions in long term monitoring of seismogenic processes.     

A probabilistic tsunami hazard assessment for the Makran subduction zone at the northwestern Indian Ocean

A probabilistic tsunami hazard assessment is performed for the Makran subduction zone (MSZ) at the northwestern Indian Ocean employing a combination of probability evaluation of offshore earthquake occurrence and numerical modeling of resulting tsunamis. In our method, we extend the Kijko and Sellevoll’s (1992) probabilistic analysis from earthquakes to tsunamis. The results suggest that the southern coasts of Iran and Pakistan, as well as Muscat, Oman are the most vulnerable areas among those studied. The probability of having tsunami waves exceeding 5 m over a 50-year period in these coasts is estimated as 17.5%. For moderate tsunamis, this probability is estimated as high as 45%. We recommend the application of this method as a fresh approach for doing probabilistic hazard assessment for tsunamis. Finally, we emphasize that given the lack of sufficient information on the mechanism of large earthquake generation in the MSZ, and inadequate data on Makran’s paleo and historical earthquakes, this study can be regarded as the first generation of PTHA for this region and more studies should be done in the future.     

Tsunami magnitudes in Taiwan, the Philippines and Indonesia

The regional characteristics of tsunami magnitudes in the SE Asia region are discussed in relation to earthquake magnitudes during the period from 1960 to 1994. Tsunami magnitudes on the Imamura-Iida scale are investigated by the author's method (Hatori 1979, 1986) using the data of inundation heights near the source area and tide-gauge records observed in Japan. The magnitude values of the Taiwan tsunamis showed relatively to be small. On the contrary, the magnitudes of tsunamis in the vicinities of the Philippines and Indonesia exceed more than 1–2 grade (tsunami heights: 2–5 times) compared to earthquakes with similar size on the circum-Pacific zone. The relation between tsunami magnitude, m, and earthquake magnitude, M _s, is expressed as m = 2.66 M _s– 17.5 for these regions. For example, the magnitudes for the 1976 Mindanao tsunami (M _s= 7.8, 3702 deaths) and the 1992 Flores tsunami (M _s= 7.5, 1713 deaths) were determined to be m = 3 and m = 2.5, respectively. The focal depth of tsunamigenic earthquakes is shallower thand< 36 km, and the detectively of tsunamis is small for deep earthquakes being d > 40 km. For future tsunamis, it is indispensable to take precautions against shallow earthquakes having the magnitudes M _s> 6.5.     

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.