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.     

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.     

Study on potential tsunami by earthquake in subduction zone of Sulawesi Sea

Indonesia is one country in the world featuring a complex tectonic structure. This condition makes earthquakes often occur in many areas of this country and as an earthquake rages beneath the sea, it will potentially trigger tsunami. One of the areas in Indonesia with a high seismic activity is Sulawesi region particularly in the Sulawesi Sea subduction zone, making it important to carry out a study on the potential tsunami at this location. The purpose of this study was to analyze the existing huge potential energy in Sulawesi Sea subduction zone and to identify tsunami modeling likely to occur based on the potential energy of the region. The approach used in assessing the tsunami disaster was the calculation of the potential energy of an earthquake and tsunami modeling based on the potential energy. The method used in this research was the least squares method for the calculation of potential energy, and near-field tsunami modeling with the assistance of TUNAMI-N2 COD. The research finding has shown that the Sulawesi Sea subduction zone has potential energy of 1.35469?×?10²³ erg, equivalent to an earthquake with a magnitude of 7.6 Mw. The tsunami modeling made shown the average wave propagation reaching ashore within 12.3 min with a height varying between 0.1 and >?3 m. The tsunami modeling also indicated that there are seven sub-districts in Buol District, Central Sulawesi, which is affected by a significant tsunami.     

Deterministic seismic hazard assessment for Sultanate of Oman

The Sultanate of Oman forms the southeastern part of the Arabian plate, which is surrounded by relatively high active tectonic zones. Studies of seismic risk assessment in Oman have been an important on-going socioeconomic concern. Using the results of the seismic hazard assessment to improve building design and construction is an effective way to reduce the seismic risk. In the current study, seismic hazard assessment for the Sultanate of Oman is performed through the deterministic approach with particular attention on the uncertainty analysis applying a recently developed method. The input data set contains a defined seismotectonic model consisting of 26 seismic zones, maximum magnitudes, and 6 alternative ground motion prediction equations that were used in four different tectonic environments: obduction zone earthquake (Zagros fold thrust belt), subduction zone earthquakes (Makran subduction zones), normal and strike-slip transform earthquakes (Owen and Gulf of Aden zones), and stable craton seismicity (Arabian stable craton). This input data set yielded a total of 76 scenarios at each point of interest. A 10 % probability that any of the 76 scenarios may exceed the largest median ground acceleration is selected. The deterministic seismic hazards in terms of PGA, 5 % damped spectral acceleration at 0.1, 0.2, 1.0 and 2.0 s are performed at 254 selected points. The ground motion was calculated at the 50th and 84th percentile levels for selected probability of exceeding the median value. The largest ground motion in the Sultanate of Oman is observed in the northeastern part of the country.     

Tsunamis from the Arica-Tocopilla source region and their effects on ports of Central Chile

The last great earthquake in northern Chile took place in 1877, and the ensuing tsunami affected not only that region but also Central Chile. For example, the Bay of Concepción, which is located 1,500 km south of the tsunami source, experienced an inundation height of around 3 m. Ports are important in the Chilean economy, due to the fact that a large percentage of Chilean exports (excluding copper) use ports located in Central Chile. With this in mind, the authors investigated the potential effect of an 1877-like tsunami on the main ports of Central Chile. To do this, the dispersive wave model Non-hydrostatic Evolution of Ocean WAVEs was used. In addition, the first tsunami forecast model for Talcahuano, inside the Bay of Concepción, was developed by means of numerical simulation of several events of different moment magnitudes. The results showed that most of the important ports (Valparaiso, San Antonio, San Vicente and Coronel) had inundation heights on the order of just 1 m, while inundation levels in Talcahuano reached up to 3.5 m. The forecast model for Talcahuano uses only earthquake magnitude, focal depth and tide level to determine tsunami inundation heights. In addition, the tsunami arrival time was computed to be 3 h, and the maximum tsunami amplitude takes place at 4 h and 45 min after the earthquake.     

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.     

Confidence levels for tsunami-inundation limits in northern Oregon inferred from a 10,000-year history of great earthquakes at the Cascadia subduction zone

To explore the local tsunami hazard from the Cascadia subduction zone we (1) evaluate geologically reasonable variability of the earthquake rupture process, (2) specify 25 deterministic earthquake sources, and (3) use resulting vertical coseismic deformations for simulation of tsunami inundation at Cannon Beach, Oregon. Maximum runup was 9–30 m (NAVD88) from earthquakes with slip of ~8–38 m and M _w ~8.3–9.4. Minimum subduction zone slip consistent with three tsunami deposits was 14–15 m. By assigning variable weights to the source scenarios using a logic tree, we derived percentile inundation lines that express the confidence level (percentage) that a Cascadia tsunami will not exceed the line. Ninety-nine percent of Cascadia tsunami variation is covered by runup ≤30 m and 90% ≤16 m with a “preferred” (highest weight) value of ~10 m. A hypothetical maximum-considered distant tsunami had runup of ~11 m, while the historical maximum was ~6.5 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.     

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.     

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.     

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.     

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.     

Spatial and temporal variations in the incidence of dust events over Iran

Large near-field tsunamis pose a significant threat to the Canadian West Coast due to its proximity to the circum-Pacific belt where a significant tsunami-inducing earthquake event from the Cascadia subduction zone is expected. This study investigated the risks associated with such an event in terms of pedestrian evacuation needs and plans for the Town of Tofino, a small community located on the West Coast of Vancouver Island. The population-at-risk within the hazard zone and its ability to evacuate to safety is evaluated using anisotropic path-distance modelling. Mitigation measures, such as vertical evacuation buildings, are quantitatively evaluated. Site-specific inundation modelling was not performed as part of this study; tsunami hazard and safe zones were computed using a range of run-ups varying between 3 and 25 m. It was established that up to 80% of the population is within the maximum hazard zone considered. This evacuation modelling exercise indicates that a maximum of 13% of the population would have insufficient time to reach safety when using a mobility-impaired ambulatory speed. The use of three vertical evacuation buildings can reduce the risk of losing population in this category by 99%. Although some conservative assumptions were used (vertical datum at higher high water, reductions in safe zones by generalization process and mobility-impaired evacuation speeds), the evacuation potential is likely overestimated due to the coarseness of the topographic data used in the evacuation modelling and from an overestimated first wave arrival time. This is the first Canadian study which used anisotropic evacuation modelling to evaluate the vulnerability of a Canadian community to tsunami inundation.

     

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.     

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.     

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.     

Modeling a tsunami from the Nicoya,Costa Rica,seismic gap and its potential impact in Puntarenas

Although subduction zones around the world are known to be the source of earthquakes and/or tsunamis, not all segments of these plate boundaries generate destructive earthquakes and catastrophic tsunamis. Costa Rica, in Central America, has subduction zones on both the Pacific and the Caribbean coasts and, even though large earthquakes (Mw = 7.4–7.8) occur in these convergent margins, they do not produce destructive tsunamis. The reason for this is that the seismogenic zones of the segments of the subduction zones that produce large earthquakes in Costa Rica are located beneath land (Nicoya peninsula, Osa peninsula and south of Limón) and not off shore as in most subduction zones around the world. To illustrate this particularity of Costa Rican subduction zones, we show in this work the case for the largest rupture area in Costa Rica (under the Nicoya peninsula), capable of producing Mw ～ 7.8 earthquakes, but the tsunamis it triggers are small and present little potential for damage even to the largest port city in Costa Rica.The Nicoya seismic gap, in NW Costa Rica, has passed its ～50-year interseismic period and therefore a large earthquake will have to occur there in the near future. The last large earthquake, in 1950 generated a tsunami which slightly affected the southwest coast of the Nicoya Peninsula. We present here a simulation to study the possible consequences that a tsunami generated by the next Nicoya earthquake could have for the city of Puntarenas. Puntarenas has a population of approximately eleven thousand people and is located on a 7.5 km long sand bar with a maximum height of 2 m above the mean sea level. This condition makes Puntarenas vulnerable to tsunamis.     

Sedimentary evidence of prehistoric distant-source tsunamis in the Hawaiian Islands

Over the past 200 years of written records, the Hawaiian Islands have experienced tens of tsunamis generated by earthquakes in the subduction zones of the Pacific ‘Ring of Fire’ (for example, Alaska–Aleutian, Kuril–Kamchatka, Chile and Japan). Mapping and dating anomalous beds of sand and silt deposited by tsunamis in low-lying areas along Pacific coasts, even those distant from subduction zones, is critical for assessing tsunami hazard throughout the Pacific basin. This study searched for evidence of tsunami inundation using stratigraphic and sedimentological analyses of potential tsunami deposits beneath present and former Hawaiian wetlands, coastal lagoons, and river floodplains. Coastal wetland sites on the islands of Hawai΄i, Maui, O΄ahu and Kaua΄i were selected based on historical tsunami runup, numerical inundation modelling, proximity to sandy source sediments, degree of historical wetland disturbance, and breadth of prior geological and archaeological investigations. Sand beds containing marine calcareous sediment within peaty and/or muddy wetland deposits on the north and north-eastern shores of Kaua΄i, O΄ahu and Hawai΄i were interpreted as tsunami deposits. At some sites, deposits of the 1946 and 1957 Aleutian tsunamis are analogues for deeper, older probable tsunami deposits. Radiocarbon-based age models date sand beds from three sites to ca 700 to 500 cal yr bp , which overlaps ages for tsunami deposits in the eastern Aleutian Islands that record a local subduction zone earthquake. The overlapping modelled ages for tsunami deposits at the study sites support a plausible correlation with an eastern Aleutian earthquake source for a large prehistoric tsunami in the Hawaiian Islands.     

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.