The Storegga tsunami along the Norwegian coast, its age and run up

The statigraphy in 25 coastal lakes shows that most of the Norwegian coastline was impacted by a large tsunami about 7200 ¹⁴C BP. The methodology has been to core a staircase of lake basins above the contemporary sea level in several areas and to map the tsunami deposit to its maximum elevation. The tsunami was identified in the sedimentary record as an erosional unconformity overlain by graded or massive sand with shell fragments, followed by redeposited organic detritus. The greatest recorded runup along the coast (10–11 m above high tide) is found in areas most proximal to the Storegga slide scar on the Norwegian continental slope (Sunnmøre). To the north and south, runup is less, about 6–7 m at Bjugn (250 km north of Sunnmøre) and about 3–5 m in Austrheim (200 km to the south of Sunnmerre). This runup pattern supports the suggestion that the tsunami was generated by the Second Storegga Slide. The recorded runup heights are consistent within and between the investigated areas, and imply that the tsunami wave was not significantly influenced by the local topography, suggesting a very long wave length. The mapped runup estimates are in good agreement with a numerical model of the tsunami generated by the Second Storegga slide, and indicate that the slide was a single major event rather than a set of smaller slides.     

Population evacuation analysis: considering dynamic population vulnerability distribution and disaster information dissemination

Abnormal tsunami amplification and runup in narrow bays is studied with respect to the Samoa tsunami of 29 September 2009. The data of the tide gauge in Pago Pago harbour are used to calculate wave runup in the city of Pago Pago (Tutuila, American Samoa) for two approximations of the bottom topography: a plane beach and a narrow bay. Theoretical estimates of tsunami runup are compared with field survey data for the 2009 Samoa tsunami. It is shown that both formulations result in equally good estimates of runup, having approximately the same difference with the field measurements. However, the narrow bay model presents more wave amplification and, consequently, runup, which is the main observation of the field survey. The differences in estimated shoreline velocity, travel time and wave breaking regime, calculated in the framework of these two approximations, are also discussed. It is concluded that wave runup in narrow bays should be calculated by the corresponding formulas, which should be taken into account by tsunami early warning systems.     

Tsunami vulnerability evaluation in the Mentawai islands based on the field survey of the 2010 tsunami

On October 25, 2010, a large earthquake occurred off the coast of the Mentawai islands in Indonesia, generating a tsunami that caused damage to the coastal area of North Pagai, South Pagai, and Sipora islands. Field surveys were conducted soon after the event by several international survey teams, including the authors’. These surveys clarified the tsunami height distribution, the damage that took place, and residents’ awareness of tsunamis in the affected islands. Heights of over 5 m were recorded on the coastal area of the Indian Ocean side of North and South Pagai islands and the south part of Sipora island. In some villages, it was difficult to evacuate immediately after the earthquake because of the lack of routes to higher ground or the presence of rivers. Residents in some villages had taken part in tsunami drills or education; however, not all villages shared awareness of tsunami threats. In the present paper, based on the results of these field surveys, the vulnerability of these islands with regards to future tsunami threats was analyzed. Three important aspects of this tsunami disaster, namely the geographic disadvantage of the islands, the resilience of buildings and other infrastructure, and people’s awareness of tsunamis, are discussed in detail, and corresponding tsunami mitigation strategies are explained.     

印度洋地震海啸(2004-12-26)及其对中国的警示

2004年12月26日印度尼西亚苏门答腊岛西北近海发生ML9级强烈地震。地震的强度是100a来全球非常罕见的。地震引起了巨大海啸,浪高近10m,波及到东南亚、南亚和东非地区10多个国家,造成近30万人遇难。地震使印度尼西亚、泰国的部分岛屿发生了地形变化。海啸在受灾国留下了大片的盐碱地。苏门答腊板块边缘的一个长距离破裂带通过长时间积累,蓄积了巨大能量。这些能量在2004-12-26集中释放出来。导致了这次地震海啸的发生。地震海啸灾害本身规模巨大,发生异常突然,再加上受灾地区人员密集,缺乏海啸灾害逃生的知识和经验。印度洋沿岸国家没有海啸预警系统,是造成这次灾害巨大伤亡的原因。中国从台湾-海南岛一线的海区,存在地震海啸的可能性。因此应不断完善海啸预警系统,提高沿海地区建设工程的防灾抗灾标准,加强防波堤建设以及采取恢复红树林等生物工程措施,预防潜在的海啸灾害。     

Run-up and Inundation Pattern Developed During the Indian Ocean Tsunami of December 26, 2004 Along the Coast of Tamilnadu (India)

The tsunami run-up, inundation and damage pattern observed along the coast of Tamilnadu (India) during the deadliest Indian Ocean tsunami of December 26, 2004 is documented in this paper. The tsunami caused severe damage and claimed many victims in the coastal areas of eleven countries, bordering the Indian Ocean. Along the coast of Indian mainland, the damage was caused by the tsunami only. Largest tsunami run-up and inundation was observed along the coast of Nagapattinam district and was about 10–12 m and 3.0 km, respectively. The measured inundation data were strongly scattered in direct relationship to the morphology of the seashore and the tsunami run-up. Lowest tsunami run-up and inundation was measured along the coast of Thanjavur, Puddukkotai and Ramnathpuram districts of Tamilnadu in the Palk Strait. The presence of shadow of Sri Lanka, the interferences of direct/receded waves with the reflected waves from Sri Lanka and Maldive Islands and variation in the width of continental shelf were the main cause of large variation in tsunami run-up along the coast of Tamilnadu.     

Far-field impact and coastal sedimentation associated with the 2006 Java tsunami in West Australia

A detailed assessment of the impact of a far-field tsunami on the Australian coastline was carried out in the Steep Point region of Western Australia following the July 17 2006 Java tsunami. Tsunami inundation and run-up were mapped on the basis of eyewitness accounts, debris lines, vegetation damage and the occurrence of recently deposited fish, starfish, corals and sea urchins well above high-tide mark. A topographic survey using kinematic GPS with accuracies of 0.02 m in the horizontal and 0.04 m in the vertical recorded flow depths of between 1 and 2 m, inundation of up to 200 m inland, and a maximum recorded run-up of 7.9 m AHD (Australian Height Datum). The tsunami impacted the sparsely populated Steep Point coastline close to low tide. It caused widespread erosion in the littoral zone, extensive vegetation damage and destroyed several campsites. Eyewitnesses reported three waves in the tsunami wave train, the second being the largest. A sand sheet, up to 14 cm thick and tapering landwards over 200 m, was deposited over coastal dunes. The deposits are predominantly composed of moderately well-sorted, medium-grained carbonate sand with some gravel and organic debris. A basal unconformity defines the boundary between tsunami sediments and underlying aeolian dune sand. Evidence for up to three individual waves is preserved as normally graded sequences mantled by layers of dark grey, organic-rich fine silty sand. Given the strong wind regimes in the area and the similarity of the underlying dune deposits to the tsunami sediments, it is likely that seasonal erosion will remove all traces of these sediment sheets within years to decades.     

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.     

Tsunami bore runup

Nearshore behaviors of tsunamis, specifically those formed as a single uniform bore, are investigated experimentally in a laboratory environment. The transition process from tsunami bore to runup is described by the momentum exchange process between the bore and the small wedge-shaped water body along the shore: the bore front itself does not reach the shoreline directly, but the large bore mass pushes the small, initially quiescent water in front of it. The fluid motions near the runup water line appear to be complex. The complex flow pattern must be caused by irregularities involved in the driving bore and turbulence advected into the runup flow. Those experimental results suggest that the tsunami actions at the shoreline involve significant mean kinetic energy together with violent turbulence. Even though the behaviors of bore motion were found to be different from those predicted by the shallow-water wave theory, the maximum runup height appears to be predictable by the theory if the value of the initial runup velocity is modified (reduced). Besides the friction effect, this reduction of the initial runup velocity must be related to the transition process as well as the highly interacting three-dimensional runup motion.     

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.     

A new tsunami runup predictor

Natural Hazards - We introduce a new parameter, tsunami runup predictor (TRP), relating the accelerating phase of the wave to the length of the beach slope over which the wave is travelling. We...     

Distribution Functions of Tsunami Wave Heights

The problem of describing the distribution functions of tsunami wave heights is discussed. Data on runup heights obtained in field surveys of several tsunamis for the last decade are used to calculate the empirical distribution functions. It is shown that the log-normal distribution describes the observed data well. This means that the irregular topography and coastline are major factors which influence the height distribution. The power distribution related with the geometric decay of the propagated wave is a good approximation for one event (Sulawesi, January 1, 1996) only. Results of a numerical simulation of the tsunami event in the Japan (East) Sea on July 12, 1993 are presented. It is shown that the computed wave height distribution, obtained by using the runup correction in the framework of nonlinear shallow-water theory, is in good agreement with the observed height distribution. Simulations are used to study the transformation of the distribution function on different distances from the source.     

Sedimentology of the December 26, 2004, Sumatra tsunami deposits in eastern India (Tamil Nadu) and Kenya

The December 26, 2004 Sumatra tsunami caused severe damage at the coasts of the Indian ocean. We report results of a sedimentological study of tsunami run-up parameters and the sediments laid down by the tsunami at the coast of Tamil Nadu, India, and between Malindi and Lamu, Kenya. In India, evidence of three tsunami waves is preserved on the beaches in the form of characteristic debris accumulations. We measured the maximum run-up distance at 580 m and the maximum run-up height at 4.85 m. Flow depth over land was at least 3.5 m. The tsunami deposited an up to 30 cm thick blanket of moderately well to well-sorted coarse and medium sand that overlies older beach deposits or soil with an erosional unconformity. The sand sheet thins inland without a decrease of grain-size. The deposits consist frequently of three layers. The lower one may be cross-bedded with foresets dipping landward and indicating deposition during run-up. The overlying two sand layers are graded or parallel-laminated without indicators of current directions. Thus, it remains undecided whether they formed during run-up or return flow. Thin dark laminae rich in heavy minerals frequently mark the contacts between successive layers. Benthic foraminifera indicate an entrainment of sediment by the tsunami from water depths less than ca. 30 m water depth. On the Indian shelf these depths are present at distances of up to 5 km from the coast. In Kenya only one wave is recorded, which attained a run-up height of 3 m at a run-up distance of ca. 35 m from the tidal water line at the time of the tsunami impact. Only one layer of fine sand was deposited by the tsunami. It consists predominantly of heavy minerals supplied to the sea by a nearby river. The sand layer thins landward with a minor decrease in grain-size. Benthic foraminifera indicate an entrainment of sediment by the tsunami from water depths less than ca. 30 m water depth, reaching down potentially to ca. 80 m. The presence of only one tsunami-related sediment layer in Kenya, but three in India, reflects the impact of only one wave at the coast of Kenya, as opposed to several in India. Grain-size distributions in the Indian and Kenyan deposits are mostly normal to slightly positively skewed and indicate that the detritus was entrained by the tsunami from well sorted pre-tsunami deposits in nearshore, swash zone and beach environments.     

Tsunami-induced scour at coastal roadways: a laboratory study

Coastal roads are lifelines for bringing emergency personnel and equipment into affected areas after tsunamis, thus careful thought should be given to how to make roadways safer from tsunamis. Scouring at roadways is the primary damage caused by tsunamis; however, tsunami-induced scouring and beach erosion are less understood compared to tsunami runup and tsunami inundation. A set of laboratory experiments are reported in this study on tsunami-induced scour at a road model situated on a sandy beach. Our experiments showed that the distance between the shoreline and a roadway, which varies with tides, was a key factor affecting the scour depth at the road. Having the coastal road at about half of the inundation distance is not the most ideal location. The depth of road embedment did not affect the scour depth in our experiments. It was also found that for typical tsunamis, the scour depth is unlikely to reach its equilibrium stage. The information reported in this study is useful for local authorities to assess potential tsunami damage of roads and to have a better plan for tsunami disaster relief.     

Storegga tsunami deposits in a coastal lake on Suouroy, the Faroe Islands

Presumed deposits of the Storegga tsunami have been recognized in a coastal lake situated 4 m a.s.l. on the island of Suðuroy, the Faroe Islands. The stratigraphy in the lake reveals a major erosion and redepositional event. The deposited material ranges from sand and sandy gyttja, with marine shell fragments and foraminifera, to gyttja with rip-up clasts, wood fragments and thin sand layers. Diatom analysis indicates that the deposit contains 5-8% polyhalobous (full marine) species, decreasing to 1-2% in the undisturbed lacustrine gyttja above. The tsunami event was dated to some time between 7300 and 6400 14 C yr BP. Lithostratigraphic profiles in the lake suggest that at least two large waves inundated the basin. The first and largest wave eroded most or all of the sediments previously deposited in the basin. The next wave caused minor erosion of the redeposited material. The waves deposited two generations of sand overlain by organic conglomerates, after which followed a unit of suspension material and normal lacustrine gyttja.     

Modeling of inundation dynamics on Banda Aceh,Indonesia during the great Sumatra tsunamis December 26, 2004

The tsunami inundation flows on Banda Aceh, Indonesia reached 5 km inland during the December 26, 2004, event and devastated most of the houses, buildings, and infrastructure along the coast and killed more than 167,000 people. The overland flows from the northwest coast and the west coast collided at Lampisang village approximately 3.7 km from Ulee Lheue (northwest coast) and 6.8 km from Lhok Nga (west coast) as reported by survivors. Inundation modeling based on the nonlinear shallow-water wave equations reproduces the inundation pattern and demonstrates a colliding of the overland flows. The model suggests that wave characteristics on the northwest coast of Banda Aceh were different from those on the waves that impacted upon the west coast. The areas, which experienced higher inundation levels, did not always experience greatest overland flow speeds, and the damage areas mostly coincide with the flow speed distribution rather than the runup and inundation depth.     

Methodology and challenges of fire following earthquake analysis: an urban community study considering water and transportation networks

Typhoon Lionrock, also known as the national number 1610 in Japan, caused severe flooding in east Japan in August 28–31, 2016, leaving a death toll of 22. With a maximum sustained wind speed of ~?220 km/h from the Joint Typhoon Warning Center’s best track, Lionrock was classified as a category 4 hurricane in Saffir–Simpson Hurricane Wind Scale and as a typhoon in Japan Meteorological Agency’s scale. Lionrock was among unique typhoons as it started its landfall from north of Japan. Here, we studied the characteristics of this typhoon through tide gauge data analysis, field surveys and numerical modeling. Tide gauge analysis showed that the surges generated by Lionrock were in the ranges of 15–55 cm with surge duration of 0.8–3.1 days. Our field surveys revealed that the damage to coastal communities/structures was moderate although it caused severe flooding inland. We measured a maximum coastal wave runup of 4.3 m in Iwaisaki. Such a runup was smaller than that generated by other category 4 typhoons hitting Japan in the past. Our numerical model was able to reproduce the storm surge generated by the 2016 Typhoon Lionrock. This validated numerical model can be used in the future for typhoon-hazard studies along the coast of northeastern Japan. Despite relatively small surge/wave runups in coastal areas, Lionrock’s death toll was more than that of some other category 4 typhoons. We attribute this to various primary (e.g., flooding, surges, waves, strong winds) and secondary (e.g., landslides, coastal erosions, debris flows, wind-blown debris) mechanisms and their combinations and interactions that contribute to damage/death during a typhoon event.

     

Probabilistic Analysis of Tsunami Hazards*

Determining the likelihood of a disaster is a key component of any comprehensive hazard assessment. This is particularly true for tsunamis, even though most tsunami hazard assessments have in the past relied on scenario or deterministic type models. We discuss probabilistic tsunami hazard analysis (PTHA) from the standpoint of integrating computational methods with empirical analysis of past tsunami runup. PTHA is derived from probabilistic seismic hazard analysis (PSHA), with the main difference being that PTHA must account for far-field sources. The computational methods rely on numerical tsunami propagation models rather than empirical attenuation relationships as in PSHA in determining ground motions. Because a number of source parameters affect local tsunami runup height, PTHA can become complex and computationally intensive. Empirical analysis can function in one of two ways, depending on the length and completeness of the tsunami catalog. For site-specific studies where there is sufficient tsunami runup data available, hazard curves can primarily be derived from empirical analysis, with computational methods used to highlight deficiencies in the tsunami catalog. For region-wide analyses and sites where there are little to no tsunami data, a computationally based method such as Monte Carlo simulation is the primary method to establish tsunami hazards. Two case studies that describe how computational and empirical methods can be integrated are presented for Acapulco, Mexico (site-specific) and the U.S. Pacific Northwest coastline (region-wide analysis). * The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Behavior of historic tsunamis of volcanic origin as revealed by onshore tsunami deposits

According to the old documents, two historic tsunamis of volcanic origin attacked Hokkaido, northern Japan. They are the 1640 Komagatake event which killed more than 700 people and the 1741 Oshima-Ohshima event which killed 1467 people. In order to obtain more information of these old tsunami disasters, we studied onshore tsunami deposits associated with these events. Tsunami deposits are identified by their sedimentary structure and granulometric characteristics. We traced the 1640 and 1741 tsunami deposits showing similar features at outcrops, by making pits or trenches. Minimum runup heights of these historic tsunamis were revealed by these tsunami deposit distributions. Trench survey is one of the best way to find and study onshore paleo-tsunami deposit     

Evaluation of tsunami impacts on shallow marine sediments: An example from the tsunami caused by the 2003 Tokachi-oki earthquake, northern Japan

A combined approach of field geology and numerical simulation was conducted for evaluating the tsunami impacts on the shelf sediments. The 2003 Tokachi-oki earthquake, M 8.0, that occurred on 25 September 2003 off southeastern Hokkaido, northern Japan, generated a locally destructive tsunami. Maximum run-up height of the tsunami waves reached 4 m above sea level. In order to estimate the tsunami impacts on shallow marine sediments, we compared pre- and post-tsunami marine sediments in water depths of 38–112 m in terms of grain size, sedimentary structure, and microfossil content. Decreases of fine fractions, especially finer than very fine sand, which led to coarsen the mean grain size, were detected in the inner shelf of the northern part of the study area. Foraminiferal assemblages also changed in the coarsened sediments. On the other hand, the other shelf sediments largely unchanged or slightly fined. We also simulated the tsunami wave velocity and direction, and grain size entrained by the modeled tsunami. The numerical simulation resulted in that the 2003 tsunami could transport very fine sand in water depths shallower than 45–95 m at the northern part of the study area. This is comparable with the actual grain-size changes after the tsunami had passed. However, some storms and tidal currents might also be possible to stir the surface sediments after the pre-tsunami survey, so we could not conclude that the grain-size changes had been caused only by the tsunami. Nevertheless, a combined approach of sampling and modeling was powerful for estimating the tsunami impacts under the sea.