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.     

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.     

An approach for estimating the largest probable tsunami from far-field subduction zone earthquakes

Following the recent unexpected earthquake events of 2004 and 2011, it can be cautiously extrapolated that all major subduction zones bearing the capacity to produce mega-earthquake events will eventually do so given enough time, irrespective of the lack of such in the relatively short historical record. This notion has led to an effort of assigning maximum earthquake magnitudes to all major subduction zones, either based on geological constraints or based on size–frequency relations, or a combination of both. In this study, we utilize the proposed maximum magnitudes to assess tsunami hazard in Central California in the very long return periods. We also assessed tsunami hazard following an alternative methodology to calculate maximum magnitudes, which uses scaling relations for subduction zone earthquakes and maximum fault rupture scenarios found in literature. A sensitivity analysis is performed for Central California that is applicable to any coastal site in the Pacific Rim and can readily provide a strong indication for which subduction zones beam the most energy toward a study area. The maximum earthquake scenarios are then narrowed down to a few candidates, for which the initial conditions are examined in more detail. The chosen worst-case scenarios for Central California stem from the Alaska–Aleutian subduction zone that beams more energy and generates the biggest amplitude waves toward the study area. The largest tsunami scenario produces maximum free surface elevations of 15 m and run-up heights greater than 20 m.     

Identification and possible recurrence of an oversized tsunami on the Pacific coast of northern Japan

Geological evidence of severe tsunami inundation has been discovered in northern Japan. In the dune fields of Shimokita, in northernmost Tohoku, we have found two distinctive sand layers that are tsunami deposits. The run-up height of >20 m and inland inundation of at least 1.4 km are notably larger than any known historical case in Japan. The tsunami-genic earthquake that resulted in these deposits is thought to have taken place in the Kuril Forearc-Trench system nearly 700 years ago. The recurrence interval of major tsunamis originating in the Kuril subduction zone is about 400 years. Given that the most recent unusually large earthquake took place in AD 1611 (corresponding to the Keicho earthquake tsunami), the findings presented here increase the potential and hazard for an outsized tsunami striking the Pacific coast of northern Japan.     

The Makassar Strait Tsunamigenic Region,Indonesia

The Makassar Strait region has had the highest frequency of historical tsunamievents for Indonesia. The strait has a seismic activity due to the convergenceof four tectonic plates that produces a complex mixture of structures. The maintsunamigenic features in the Makassar Strait are the Palu-Koro and Pasternostertransform fault zones, which form the boundaries of the Makassar trough.Analysis of the seismicity, tectonics and historic tsunami events indicatesthat the two fault zones have different tsunami generating characteristics.The Palu-Koro fault zone involves shallow thrust earthquakes that generatetsunami that have magnitudes that are consistent with the earthquakemagnitudes. The Pasternoster fault zone involves shallower strike-slipearthquakes that produce tsunami magnitudes larger than would normallybe expected for the earthquake magnitude. The most likely cause for theincreased tsunami energy is considered to be submarine landslidesassociated with the earthquakes. Earthquakes from both fault zonesappear to cause subsidence of the west coast of Sulawesi Island.The available data were used to construct a tsunami hazard map whichidentifies the highest risk along the west coast of Sulawesi Island.The opposite side of the Makassar Strait has a lower risk because it isfurther from the historic tsunami source regions along the Sulawesicoast, and because the continental shelf dissipates tsunami wave energy.The greatest tsunami risk for the Makassar Strait is attributed tolocally generated tsunami due to the very short travel times.     

The 1996 Sulawesi Tsunami

On 1 January, 1996 at 16:05 p.m. local time, an earthquake of magnitude M = 7.8 struck the central part of Sulawesi Island (Indonesia). It was accompanied by tsunami waves 2–4 m high. Nine people were killed and 63 were injured. A tsunami survey was conducted by Indonesian and Russian specialists. The measured tsunami runup heights and eyewitness accounts are reported and discussed. Historical data on the Sulawesi Island tsunamis are analysed and tsunami risk prediction in the central part of Sulawesi Island carried out for the first time.     

2018年帕卢MW7.6地震震源及苏拉威西地区构造应力场特征

为分析帕卢地震的发震构造及构造应力场特征,基于1976年1月至2018年4月全球CMT目录的震源机制解,反演了2018年9月28日印度尼西亚帕卢地震震源区及苏拉威西地区构造应力场,得到以下初步认识:赤道0°两侧南、北区域的应力场呈现出整体的一致性和局部的不均匀性,北部是以向北倾伏的低倾伏角挤压和向南倾伏的高倾伏角拉张作用为主;南部是近水平的NNW-EW向挤压及低倾伏角的近N-S向拉张的应力状态.摆堆易逆冲推挤及西里伯斯海阻挡导致北部地区的应力场显示俯冲带的特征,也导致北苏拉和马纳都块体的应力场出现差异性.在班达海自东向西的挤压及西部马卡萨盆地的阻挡作用下,南部地区呈现出近E-W向的挤压状态,并沿着NNE-SSW向呈现出地壳物质的构造逃逸现象.深部地幔物质上涌致使苏拉威西中部的R值偏小,形成活火山.受到桑吉双俯冲带和麦纳哈撒海沟的推挤及火山喷发作用,东部交汇区的应力场呈较大倾伏角的拉张为主及R值复杂突变.西部交汇地区的应力场是走滑为主兼有正断应力机制,该应力场有益于此次地震的触发及超剪切波破裂现象的出现.      

Tsunami hazard assessment of Canada

We present a preliminary probabilistic tsunami hazard assessment of Canadian coastlines from local and far-field, earthquake, and large submarine landslide sources. Analyses involve published historical, palaeotsunami and palaeoseismic data, modelling, and empirical relations between fault area, earthquake magnitude, and tsunami run-up. The cumulative estimated tsunami hazard for potentially damaging run-up (≥1.5 m) of the outer Pacific coastline is ~40–80 % in 50 years, respectively one and two orders of magnitude greater than the outer Atlantic (~1–15 %) and the Arctic (<1 %). For larger run-up with significant damage potential (≥3 m), Pacific hazard is ~10–30 % in 50 years, again much larger than both the Atlantic (~1–5 %) and Arctic (<1 %). For outer Pacific coastlines, the ≥1.5 m run-up hazard is dominated by far-field subduction zones, but the probability of run-up ≥3 m is highest for local megathrust sources, particularly the Cascadia subduction zone; thrust sources further north are also significant, as illustrated by the 2012 Haida Gwaii event. For Juan de Fuca and Georgia Straits, the Cascadia megathrust dominates the hazard at both levels. Tsunami hazard on the Atlantic coastline is dominated by poorly constrained far-field subduction sources; a lesser hazard is posed by near-field continental slope failures similar to the 1929 Grand Banks event. Tsunami hazard on the Arctic coastline is poorly constrained, but is likely dominated by continental slope failures; a hypothetical earthquake source beneath the Mackenzie delta requires further study. We highlight areas susceptible to locally damaging landslide-generated tsunamis, but do not quantify the hazard.     

Seismogenic sources in the Bay of Bengal vis-à-vis potential for tsunami generation and its impact in the northern Bay of Bengal coast

Geodynamic status, seismo-tectonic environment, and geophysical signatures of the Bay of Bengal do not support the occurrence of seismogenic tsunami. Since thrust fault and its intensity and magnitude of rupture are the key tectonic elements of tsunamigenic seismic sources, the study reveals that such characteristics of fault-rupture and seismic sources do not occur in most of the Bay of Bengal except a small segment in the Andaman–Nicobar subduction zone. The inferred segment of the Andaman–Nicobar subduction zone is considered for generating a model of the deformation field arising from fluid-driven source. The model suggests local tsunami with insignificant inundation potential along the coast of northern Bay of Bengal. The bathymetric profile and the sea floor configuration of the northern Bay of Bengal play an important role in flattening the waveform through defocusing process. The direction of motion of the Indian plate makes an angle of about 30° with the direction of the opening of Andaman Sea. The opening of Andaman Sea and the direction of plate motion of the Indian plate results in the formation of Andaman trench where the subducting plate dives more obliquely than that in the Sunda trench in the south. The oblique subduction reduces significantly the possibilities of dominant thrust faulting in the Andaman subduction zone. Further, north of Andaman subduction in the Bengal–Arakan coast, there is no active subduction. On the otherhand, much greater volume of sediments (in excess of 20 km) in the Bengal–Arakan segment reduces the possibilities of mega rupture of the ocean floor. The water depth (≈1,000 m) along most of the northern Bay of Bengal plate margin is not optimum for any significant tsunami generation. Hence, very weak possibility of any significant tsunami is suggested that based on the interpretation of geodynamic status, seismo-tectonic environment, and geophysical signatures of the Andaman subduction zone and the Bengal–Arakan coast.     

A data-driven approach to assess large fire size generation in Greece

The 2011 Tohoku earthquake and tsunami motivated an analysis of the potential for great tsunamis in Hawai‘i that significantly exceed the historical record. The largest potential tsunamis that may impact the state from distant, Mw 9 earthquakes—as forecast by two independent tsunami models—originate in the Eastern Aleutian Islands. This analysis is the basis for creating an extreme tsunami evacuation zone, updating prior zones based only on historical tsunami inundation. We first validate the methodology by corroborating that the largest historical tsunami in 1946 is consistent with the seismologically determined earthquake source and observed historical tsunami amplitudes in Hawai‘i. Using prior source characteristics of Mw 9 earthquakes (fault area, slip, and distribution), we analyze parametrically the range of Aleutian–Alaska earthquake sources that produce the most extreme tsunami events in Hawai‘i. Key findings include: (1) An Mw 8.6 ± 0.1 1946 Aleutian earthquake source fits Hawai‘i tsunami run-up/inundation observations, (2) for the 40 scenarios considered here, maximal tsunami inundations everywhere in the Hawaiian Islands cannot be generated by a single large earthquake, (3) depending on location, the largest inundations may occur for either earthquakes with the largest slip at the trench, or those with broad faulting over an extended area, (4) these extremes are shown to correlate with the frequency content (wavelength) of the tsunami, (5) highly variable slip along the fault strike has only a minor influence on inundation at these tele-tsunami distances, and (6) for a given maximum average fault slip, increasing the fault area does not generally produce greater run-up, as the additional wave energy enhances longer wavelengths, with a modest effect on inundation.     

Maritime tsunami evacuation guidelines for the Pacific Northwest coast of Oregon

Recent tsunamis affecting the West Coast of the USA have resulted in significant damage to ports and harbors, as well as to recreational and commercial vessels attempting to escape the tsunami. With the completion of tsunami inundation simulations for a distant tsunami originating from the Aleutian Islands and a locally generated tsunami on the Cascadia subduction zone (CSZ), the State of Oregon is now able to provide guidance on the magnitudes and directions of the simulated currents for the Oregon coast and shelf region. Our analyses indicate that first wave arrivals for an Aleutian Island event would take place on the north coast,?~?3 h 40 min after the start of the earthquake,?~?20 min later on the southern Oregon coast. The simulations demonstrated significant along-coast variability in both the tsunamis water levels and currents, caused by localized bathymetric effects (e.g., submarine banks and reefs). A locally generated CSZ event would reach the open coast within 7–13 min; maximum inundation occurs at?~?30–40 min. As the tsunami current velocities increase, the potential for damage in ports and harbors correspondingly increases, while also affecting a vessels ability to maintain control out on the ocean. Scientific consensus suggests that tsunami currents?<?1.54 m/s are unlikely to impact maritime safety in ports and harbors. No such guidance is available for boats operating on the ocean, though studies undertaken in Japan suggest that velocities in the region of 1–2 m/s may be damaging to boats. In addition to the effects of currents, there is the added potential for wave amplification of locally generated wind waves interacting with opposing tsunami currents in the offshore. Our analyses explore potential wave amplification effects for a range of generic sea states, ultimately producing a nomogram of wave amplification for a range of wave and opposing current conditions. These data will be useful for US Coast Guard and Port authorities as they evaluate maritime tsunami evacuation options for the Oregon coast. Finally, we identify three regions of hazard (high, moderate, and low) across the Oregon shelf, which can be used to help guide final designation of tsunami maritime evacuation zones for the coast.     

Tsunami inundation in Napier, New Zealand, due to local earthquake sources

Deterministic analysis of local tsunami generated by subduction zone earthquakes demonstrates the potential for extensive inundation and building damage in Napier, New Zealand. We present the first high-resolution assessments of tsunami inundation in Napier based on full simulation from tsunami generation to inundation and demonstrate the potential variability of onshore impacts due to local earthquakes. In the most extreme scenario, rupture of the whole Hikurangi subduction margin, maximum onshore flow depth exceeds 8.0 m within 200 m of the shore and exceeds 5.0 m in the city centre, with high potential for major damage to buildings. Inundation due to single-segment or splay fault rupture is relatively limited despite the magnitudes of M_W 7.8 and greater. There is approximately 30 min available for evacuation of the inundation zone following a local rupture, and inundation could reach a maximum extent of 4 km. The central city is inundated by up to three waves, and Napier Port could be inundated repeatedly for 12 h. These new data on potential flow depth, arrival time and flow kinematics provide valuable information for tsunami education, exposure analysis and evacuation planning.     

Effects of the October 25, 2010, Mentawai tsunami in La Réunion Island (France): observations and crisis management

On October 25, 2010, a Mw = 7.7 magnitude earthquake was recorded in the Kepulauan Mentawai archipelago, in Indonesia. Following the earthquake, the Pacific Tsunami Warning Center issued a tsunami bulletin to the countries of the Indian Ocean informing them that a local tsunami watch was in effect for Indonesia. La Réunion Island, a French territory located 5,000 km southwest of the seismic source, did not trigger its tsunami early warning protocol, but did suffer the impact of this tsunami 7h20 post-earthquake. This paper presents the post-tsunami survey conducted in La Réunion Island a few days after the event, providing run-up measurements that reach 1.72 m locally as well as information about the material losses and the crisis management. It provides a time sequence of the observations along the coast and the reactions of the population and authorities in dealing with this crisis. This paper also shows the local authorities’ lack of action in managing this event—thus increasing human and material vulnerabilities—the consequences for victims and the impact on the credibility of the local early warning system. These results provide geophysicists with information to better describe the source of this tsunami through far-field inversion and provide local decision makers with a comprehensive assessment of the weaknesses of the local early warning system and protocols.     

Assessing the threat to Western Australia from tsunami generated by earthquakes along the Sunda Arc

A suite of tsunami spaced evenly along the subduction zone to the south of Indonesia (the Sunda Arc) were numerically modelled in order to make a preliminary estimate of the level of threat faced by Western Australia from tsunami generated along the Arc. Offshore wave heights from these tsunami were predicted to be significantly higher along the northern part of the west Australian coast than for the rest of the coast south of the town of Exmouth. In particular, the area around Exmouth may face a higher tsunami hazard than other areas of the West Australian coast nearby. Large earthquakes offshore of Java and Sumbawa are likely to be a greater hazard to WA than those offshore of Sumatra. Our numerical models indicate that a magnitude 9 or above earthquake along the eastern part of the Sunda Arc has the potential to significantly impact a large part of the West Australian coastline. The Australian government reserves the right to retain a non-exclusive, royalty free license in and to any copyright.     

Tertiary evolution of the Eastern Indonesia Collision Complex

Eastern Indonesia is the zone of interaction between three converging megaplates: Eurasia, the Pacific and Indo-Australia. The geological basis for interpretations of the Tertiary tectonic evolution of Eastern Indonesia is reviewed, and a series of plate tectonic reconstructions for this region at 5 million year intervals covering the last 35 million years is presented.The oldest reconstruction predates the onset of regional collisional deformation. At this time a simple plate configuration is interpreted, consisting of the northward-moving Australian continent approaching an approximately E–W oriented, southward-facing subduction zone extending from the southern margin of the Eurasian continent eastwards into the Pacific oceanic domain. Beginning at about 30 Ma the Australian continental margin commenced collision with the subduction zone along its entire palinspastically-restored northern margin, from Sulawesi in the west to Papua New Guinea in the east. From this time until ca 24 Ma, the Australian continent indented the former arc trend, with the northward convergence of Australia absorbed at the palaeo-northern boundary of the Philippine Sea Plate (the present-day Palau-Kyushu Ridge).At ca 24 Ma the present-day pattern of oblique convergence between the northern margin of Australia and the Philippine Sea Plate began to develop. At about this time a large portion of the Palaeogene colliding volcanic arc (the future eastern Philippines) began to detach from the northern continental margin by left-lateral strike slip. From ca 18 Ma oblique southward-directed subduction commenced at the Maramuni Arc in northern New Guinea. At ca 12 Ma the Sorong Fault Zone strike-slip system developed, effectively separating the Philippines from the Indonesian tectonic domain. The Sorong Fault Zone became inactive at ca 6 Ma, since which time the tectonics of eastern Indonesia has been dominated by the anticlockwise rotation of the Bird’s Head structural block by some 30–40°.Contemporaneously with post-18 Ma tectonism, the Banda Arc subduction–collision system developed off the northwestern margin of the Australian continent. Convergence between Indo-Australia and Eurasia was accommodated initially by northward subduction of the Indian Ocean, and subsequently, since ca 8 Ma, by the development of a second phase of arc-continent collision around the former passive continental margin of NW Australia.     

A semi-probabilistic procedure for developing societal risk function

Among the coastal districts of mega city Istanbul, Bakirkoy is one of the most critical one with the importance of air and marine transportation and presence of many other coastal facilities and structures that are prone to suffer from marine hazards. In the history, the Sea of Marmara has experienced numerous earthquake and landslide events and associated tsunamis. Therefore, tsunami risk assessment is essential for all coastal districts of Istanbul, including Bakirkoy district. In this study, a further developed methodology for tsunami human vulnerability and risk assessment Metropolitan Tsunami Human Vulnerability Assessment (MeTHuVA) is applied for Bakirkoy district of Istanbul, considering earthquake generated tsunamis. High-resolution tsunami hazard analysis is performed with the integration of coastal inundation computation with tsunami numerical tool NAMI DANCE and tsunami human vulnerability assessment with GIS-based multi-criteria decision analysis methods (MCDA). Using analytical hierarchy process method of MCDA, a hierarchical structure is established, composed of two main elements of tsunami human vulnerability: Vulnerability at Location and Evacuation Resilience. Tsunami risk assessment for Bakirkoy district is calculated by integrating result of hazard and vulnerability assessments with a risk relation that includes a parameter (n), which represents the preparedness and awareness level of the community. Tsunami simulations revealed that the maximum inundation distance is over 350 m on land and water penetrates almost 1700 m along Ayamama stream. Inundation is observed in eleven neighborhoods of Bakirkoy district. In the inundation zone, maximum flow depth is found to be over 5.7 m. The inundated area forms 4.2% of whole Bakirkoy district, and 62 buildings are located in the inundation zone. Hazard, vulnerability and risk assessment results for different neighborhoods of Bakirkoy district are presented and discussed.     

Recurrence intervals of major paleotsunamis as calibrated by historic tsunami deposits in three localities: Port Alberni, Cannon Beach, and Crescent City, along the Cascadia margin, Canada and USA

We review geologic records of both historic and prehistoric tsunami inundations at three widely separated localities that experienced significant damage from the 1964 Alaskan tsunami along the Cascadia margin. The three localities are Port Alberni, Cannon Beach, and Crescent City, representing, respectively, the north, central, and south portions of the study area (1,000 km in length). The geologic records include anomalous sand sheets from marine surges that are hosted in supratidal peaty mud deposits. Paleotsunami sand sheets that exceed the thickness, continuity and/or extent of the 1964 historic tsunami are counted as major paleotsunami inundations. Major paleotsunamis (6–7 in number) at each locality during the last 3,000 years demonstrate mean recurrence intervals of 450–540 years, and within-cluster intervals (three events each) of 270–460 years. It has been 313 years since the last major paleotsunami from a great Cascadia earthquake in AD 1700. We compare the dated sequences of major paleotsunami inundations to the nearest regional records of coastal coseismic subsidence in Willapa Bay in the central margin, Waatch/Neah Bay in the northern margin, and Coquille in the southern margin. Similar numbers of events from both types of records suggest that the major paleotsunamis are locally derived (near-field) from ruptures of the Cascadia margin megathrust fault zone, rather than from transoceanic tsunamis (far-field) originating at other subduction zones around the Pacific Rim. Given the catastrophic hazard of the near-field Cascadia margin tsunamis, we propose a basic rule for reminding the general public of the need for self-initiated evacuation following a great earthquake at the Cascadia margin.     

Application of fragility curves to estimate building damage and economic loss at a community scale: a case study of Seaside,Oregon

Community-scale estimates of building damage and economic loss are modeled for Seaside, Oregon, for Cascadia subduction zone events ranging from 8.7 to 9.3 M_W with corresponding slip distances of 3–25 m considering only the effects of the tsunami. Numerical simulations are obtained from the National Oceanic and Atmospheric Administration’s method of splitting tsunami model which includes a source model, subsidence, and calculations of the propagation and inundation flow characteristics. The damage estimates are based on fragility curves from the literature which relate flow depth with probability of damage for two different structural materials of buildings. Calculations are performed at the parcel level for the inundation hazard without including damage caused by the earthquake itself. Calculations show that the severity of building damage in Seaside is sensitive to the magnitude of the event or degree of slip because the majority of the city is located on low-lying coastal land within the estimated inundation zone. For the events modeled, the percentage of building within the inundation zone ranges from 9 to 88 %, with average direct economic losses ranging from $2 million to $1.2 billion.     

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.     

Nine unusually large tsunami deposits from the past 4000 years at Kiritappu marsh along the southern Kuril Trench

Large earthquakes along the Kuril subduction zone in northern Japan are known to have caused damaging tsunami, although there is a little information on historical earthquakes and tsunami in this area because no documents exist before the 19th century that might refer to tsunami events. To determine the likely timing and size of future events we need information on their recurrence intervals and to do this for the prehistoric past we have investigated sediments located in the Kiritappu marsh in eastern Hokaido that we interpret as laid down by tsunami. Using reliable multiple lines of evidence from sedimentological, geomorphological, micropaleontological, and chronological results, we identify 13 tsunami sands. Two of these lie within a peat bed above a historical tephra, Ta-a (AD 1739); the upper one probably corresponds to the AD 1843 Tempo Tokachi-oki earthquake (M 8.2) tsunami, and the lower to either the AD 1952 Tokachi-oki earthquake (M 8.2) tsunami or the AD 1960 Chilean earthquake (M 9.5) tsunami. Underlying are 11 prehistoric tsunami sand beds (nine large sand beds and two smaller sand beds) deposited during the past 4000 years. Because of the wide spatial distribution of the large sand beds, and inundation distances inland of between 1200 to 3000 m, we suggest that they record unusually large tsunamis along the Kuril subduction zone. According to our analyses, these tsunami sands were derived from the coastal area and, although they do not show clear graded bedding, they commonly have gradational upper boundaries and erosional bases and include internal sedimentary structures such as plane beds, dunes, and current ripples, reflecting bedload transportation. Based on our results we calculate the recurrence interval of unusually large earthquakes (probably M 8.6) along the Kuril subduction zone as about 365–553 years and estimate the youngest large event to have occurred in the 17th century.