The Transoceanic 1755 Lisbon Tsunami in Martinique

On 1 November 1755, a major earthquake of estimated M _w=8.5/9.0 destroyed Lisbon (Portugal) and was felt in the whole of western Europe. It generated a huge transoceanic tsunami that ravaged the coasts of Morocco, Portugal and Spain. Local extreme run-up heights were reported in some places such as Cape St Vincent (Portugal). Great waves were reported in the Madeira Islands, the Azores and as far as the Antilles (Caribbean Islands). An accurate search for historical data allowed us to find new (unpublished) information concerning the tsunami arrival and its consequences in several islands of the Lesser Antilles Arc. In some places, especially Martinique and the Guadeloupe islands, 3?m wave heights, inundation of low lands, and destruction of buildings and boats were reported (in some specific locations probably more enclined to wave amplification). In this study, we present the results of tsunami modeling for the 1755 event on the French island of Martinique, located in the Lesser Antilles Arc. High resolution bathymetric grids were prepared, including topographic data for the first tens of meters from the coastline, in order to model inundations on several sites of Martinique Island. In order to reproduce as well as possible the wave coastal propagation and amplification, the final grid was prepared taking into account the main coastal features and harbour structures. Model results are checked against historical data in terms of wave arrival, polarity, amplitude and period and they correlate well for Martinique. This study is a contribution to the evaluation of the tele-tsunami impact in the Caribbean Islands due to a source located offshore of Iberia and shows that an 8.5 magnitude earthquake located in the northeastern Atlantic is able to generate a tsunami that could impact the Caribbean Islands. This fact must be taken into account in hazard and risk studies for this area.     

Numerical Modelling of Tsunamis Generated by Hypothetical Landslides in the Strait of Georgia,British Columbia

--A modified and corrected version of the viscous slide model of Jiang and LeBlond (1994) is used to assess the tsunami risk associated with hypothetical underwater slope failures in two coastal areas of British Columbia having potentially unstable sediment deposits: (a) Malaspina Strait, separating the mainland coast and Texada Island in the central Strait of Georgia; and (b) Roberts Bank on the foreslope of the Fraser River Delta in the southern Strait of Georgia. The intent of this study is to demonstrate the capability of the model for tsunami risk assessment and to improve upon previous studies of tsunami risk in the region based on reasonable submarine landslide scenarios. The potential risk from tsunamis associated with slide failures has been examined, but the likelihood of failure events themselves was not considered. For the Malaspina Strait scenarios, simulated tsunamis are generated by failure of a lobe of perched sediment situ ated on the slope of eastern Texada Island. Failure as a flow slide of the estimated 1,250,000 m³ of sediment generates wave troughs reaching ц.9 m and trough-to-crest heights of 6 to 8 m along the coast of Texada Island. At Cape Cockburn, on the opposite side of the strait, wave heights of 1.5 to 2.0 m are produced. For Roberts Bank, simulated waves are examined for two separate failure scenarios. The larger slide (Case 1) involves the failure of a sediment lobe with lateral dimensions of 7 2 3 km² and volume of 0.75 km³ while the smaller slide (Case 2) fails a sediment lobe with dimensions of 4 2 2.6 km² and volume of 0.23 km³. Computations were made both for high (+3 m) and low (х m) tide conditions. For both failure volumes, maximum wave amplitudes (up to 18 m for Case 1 and 8 m for Case 2) occur on the coasts of Mayne and Galiano Islands, opposite the source area. Wave amplitudes are much smaller (1 to 4 m) on the mainl and coast because of the reflection of the initial waves from Roberts Bank. Additional numerical experiments were conducted for both regions to estimate the sensitivity of the computed tsunami wave heights to input parameters, such as slide viscosity, bulk density, and slide position.     

The effects of the 2004 tsunami on a coastal aquifer in Sri Lanka

On December 26, 2004, the earthquake off the southern coast of Sumatra in the Indian Ocean generated far-reaching tsunami waves, resulting in severe disruption of the coastal aquifers in many countries of the region. The objective of this study was to examine the impact of the tsunami on groundwater in coastal areas. Field investigations on the east coast of Sri Lanka were carried out along a transect located perpendicular to the coastline on a 2.4 km wide sand stretch bounded by the sea and a lagoon. Measurements of groundwater table elevation and electrical conductivity (EC) of the groundwater were carried out monthly from October 2005 to August 2007. The aquifer system and tsunami saltwater intrusion were modeled using the variable-density flow and solute transport code HST3D to understand the tsunami plume behavior and estimate the aquifer recovery time. EC values reduced as a result of the monsoonal rainfall following the tsunami with a decline in reduction rate during the dry season. The upper part of the saturated zone (down to 2.5 m) returned to freshwater conditions (EC < 1000 μS/cm) 1 to 1.5 years after the tsunami, according to field observations. On the basis of model simulations, it may take more than 15 years for the entire aquifer (down to 28 m) to recover completely, although the top 6 m of the aquifer may become fresh in about 5 years.     

Tsunami Modeling to Validate Slip Models of the 2007 M w 8.0 Pisco Earthquake, Central Peru

Following the 2007, August 15th, M _w 8.0, Pisco earthquake in central Peru, Sladen et al. (J Geophys Res 115: B02405, 2010) have derived several slip models of this event. They inverted teleseismic data together with geodetic (InSAR) measurements to look for the co-seismic slip distribution on the fault plane, considering those data sets separately or jointly. But how close to the real slip distribution are those inverted slip models? To answer this crucial question, the authors generated some tsunami records based on their slip models and compared them to DART buoys, tsunami records, and available runup data. Such an approach requires a robust and accurate tsunami model (non-linear, dispersive, accurate bathymetry and topography, etc.) otherwise the differences between the data and the model may be attributed to the slip models themselves, though they arise from an incomplete tsunami simulation. The accuracy of a numerical tsunami simulation strongly depends, among others, on two important constraints: (i) A fine computational grid (and thus the bathymetry and topography data sets used) which is not always available, unfortunately, and (ii) a realistic tsunami propagation model including dispersion. Here, we extend Sladen’s work using newly available data, namely a tide gauge record at Callao (Lima harbor) and the Chilean DART buoy record, while considering a complete set of runup data along with a more realistic tsunami numerical that accounts for dispersion, and also considering a fine-resolution computational grid, which is essential. Through these accurate numerical simulations we infer that the InSAR-based model is in better agreement with the tsunami data, studying the case of the Pisco earthquake indicating that geodetic data seems essential to recover the final co-seismic slip distribution on the rupture plane. Slip models based on teleseismic data are unable to describe the observed tsunami, suggesting that a significant amount of co-seismic slip may have been aseismic. Finally, we compute the runup distribution along the central part of the Peruvian coast to better understand the wave amplification/attenuation processes of the tsunami generated by the Pisco earthquake.     

Best-fit distribution and log-normality for tsunami heights along coastal lines

The best-fit distribution of the tsunami height was investigated along the Eastern Coast of Korean Peninsula. Firstly, the tsunami heights corresponding to the nine probable undersea earthquakes were obtained along the coastline using the numerical simulation. The method of L-moment ratio diagram was used to identify the best-fit probability density function of the tsunami heights caused by each undersea earthquake. The result indicates the generalized Pareto distribution is the best-fit distribution representing the tsunami heights regardless of the characteristics of the undersea earthquakes. This is particularly because the area of high tsunami heights and its relative magnitude to the adjacent locations were similar for the most simulations cases. In addition, this study further investigated the reason why the tsunami height distribution is not represented by the log-normal (LN) distribution as suggested by the previous studies. Result of the investigation indicates that the log-normality of the tsunami heights can be preserved when the length of a coastal line is not long such that the homogeneity of the length of the wave propagation paths reaching at different locations of the coastal line is preserved. This subsequently secures the central limit theorem making the distribution of the tsunami heights have the LN distribution. As the length of the coastal line increases, the deviation of the tsunami height distribution from the log-normality increases.     

Linear and nonlinear computations of the 1992 Nicaragua earthquake tsunami

Numerical computations of tsunamis are made for the 1992 Nicaragua earthquake using different governing equations, bottom frictional values and bathymetry data. The results are compared with each other as well as with the observations, both tide gauge records and runup heights. Comparison of the observed and computed tsunami waveforms indicates that the use of detailed bathymetry data with a small grid size is more effective than to include nonlinear terms in tsunami computation. Linear computation overestimates the amplitude for the later phase than the first arrival, particularly when the amplitude becomes large. The computed amplitudes along the coast from nonlinear computation are much smaller than the observed tsunami runup heights; the average ratio, or the amplification factor, is estimated to be 3 in the present case when the grid size of 1 minute is used. The factor however may depend on the grid size for the computation.     

Slip Distribution and Seismic Moment of the 2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms and Coastal Geodetic Data

The slip distribution and seismic moment of the 2010 and 1960 Chilean earthquakes were estimated from tsunami and coastal geodetic data. These two earthquakes generated transoceanic tsunamis, and the waveforms were recorded around the Pacific Ocean. In addition, coseismic coastal uplift and subsidence were measured around the source areas. For the 27 February 2010 Maule earthquake, inversion of the tsunami waveforms recorded at nearby coastal tide gauge and Deep Ocean Assessment and Reporting of Tsunamis (DART) stations combined with coastal geodetic data suggest two asperities: a northern one beneath the coast of Constitucion and a southern one around the Arauco Peninsula. The total fault length is approximately 400 km with seismic moment of 1.7 × 10²² Nm (Mw 8.8). The offshore DART tsunami waveforms require fault slips beneath the coasts, but the exact locations are better estimated by coastal geodetic data. The 22 May 1960 earthquake produced very large, ~30 m, slip off Valdivia. Joint inversion of tsunami waveforms, at tide gauge stations in South America, with coastal geodetic and leveling data shows total fault length of ~800 km and seismic moment of 7.2 × 10²² Nm (Mw 9.2). The seismic moment estimated from tsunami or joint inversion is similar to previous estimates from geodetic data, but much smaller than the results from seismic data analysis.     

Modeling the seismic source and tsunami generation of the December 12, 1992 Flores Island,Indonesia, earthquake

On December 12, 1992 a large earthquake (M _s 7.5) occurred just north of Flores Island, Indonesia which, along with the tsunami it generated, killed more than 2,000 people. In this study, teleseismicP andSH waves, as well asPP waves from distances up to 123°, are inverted for the orientations and time histories of multiple point sources. By repeating the inversion for reasonable values of depth, time separation and spatial separation, a 2-fault model is developed. Next, the vertical deformation of the seafloor is estimated from this fault model. Using a detailed bathymetric model, linear and nonlinear tsunami propagation models are tested. The data consist of a single tide gauge record at Palopo (650 km to the north), as well as tsunami runup height measurements from Flores Island and nearby islands. Assuming a tsunami runup amplification factor of two, the two-fault model explains the tide gauge record and the tsunami runup heights on most of Flores Island. It cannot, however, explain the large tsunami runup heights observed near Leworahang (on Hading Bay) and Riangkroko (on the northeast peninsula). Massive coastal slumping was observed at both of these locations. A final model, which in addition to the two faults, includes point sources of large vertical displacement at these two locations explains the observations quite well.     

Tsunami Sediment Characteristics at the Thai Andaman Coast

This paper describes and summarizes the 2004 Indian Ocean tsunami sediment characteristics at the Thai Andaman coast. Field investigations have been made approximately 3 years after the 2004 Indian Ocean tsunami event. Seven transects have been examined at five locations. Sediment samples have been collected for grain-size analyses by wet-sieve method. Tsunami sediments are compared to three deposits from coastal sub-environments. The mean grain-size and standard deviation of deposits show that shoreface deposits are fine to very fine sand, poorly to moderately well sorted; swash zone deposits are coarse to fine sand, poorly to well sorted; berm/dune deposits are medium to fine sand, poorly to well sorted; and tsunami deposits are coarse to very fine sand, poorly to moderately well sorted. A plot of deposit mean grain-size versus sorting indicates that tsunami deposits are composed of shoreface deposits, swash zone deposits and berm/dune deposits as well. The tsunami sediment is a gray sand layer deposited with an erosional base on a pre-existing soil (rooted soil). The thickness of the tsunami sediment layer is variable. The best location for observation of the recent tsunami sediment is at about 50–200 m inland from the coastline. In most cases, the sediment layer is normally graded. In some cases, the sediment contains rip-up clasts of muddy soils and/or organic matter. The vertical variation of tsunami sediment texture shows that the mean grain-size is fining upward and landward. Break points of slope in a plot of standard deviation versus depth mark a break in turbulence associated with a transition to a lower or higher Reynolds number runup. This can be used to evaluate tsunami sediment main layer and tsunami sediment sub layers. The skewness of tsunami sediment indicates a grain size distribution with prominent finer-grain or coarse-grain particles. The kurtosis of tsunami sediment indicates grain-size distributions which are flat to peak distribution (or multi-modal to uni-modal distribution) upward. Generally, the major origins of tsunami sediment are swash zone and berm/dune zone sands where coarse to medium sands are the significant material at these locations. The minor origin of tsunami sediment is the shoreface where the significant materials are fine to very fine sands. However, for a coastal area where the shoreface slope is mild, the major origin of tsunami sediment is the shoreface. The interpretation of runup number from tsunami sediment characteristics gets three runups for the 2004 Indian Ocean tsunami at the Thai Andaman coast. It corresponds to field observations from local eyewitnesses. The 1st runup transported and deposited more coarse particles than the following runups. Overall, the pattern of onshore tsunami sediment transportation indicates erosion at swash zone and berm/dune zone, followed by dynamic equilibrium at an area behind the berm/dune zone and after that deposition at inland zone until the limit of sediment inundation. The total deposition is a major pattern in onshore tsunami sediment transportation at the deposition zone which the sediment must find in the direction of transport.     

Numerical Simulation of the 2011 Tohoku Tsunami Based on a New Transient FEM Co-seismic Source: Comparison to Far- and Near-Field Observations

In this work, we simulate the 2011 M9 Tohoku-Oki tsunami using new coseismic tsunami sources based on inverting onshore and offshore geodetic data, using 3D Finite Element Models (FEM). Such FEMs simulate elastic dislocations along the plate boundary interface separating the stiff subducting Pacific Plate from the relatively weak forearc and volcanic arc of the overriding Eurasian plate. Due in part to the simulated weak forearc materials, such sources produce significant shallow slip (several tens of meters) along the updip portion of the rupture near the trench. To assess the accuracy of the new approach, we compare observations and numerical simulations of the tsunami's far- and near-field coastal impact for: (i) one of the standard seismic inversion sources (UCSB; Shao et al. 2011); and (ii) the new FEM sources. Specifically, results of numerical simulations for both sources, performed using the fully nonlinear and dispersive Boussinesq wave model FUNWAVE-TVD, are compared to DART buoy, GPS tide gauge, and inundation/runup measurements. We use a series of nested model grids with varying resolution (down to 250 m nearshore) and size, and assess effects on model results of the latter and of model physics (such as when including dispersion or not). We also assess the effects of triggering the tsunami sources in the propagation model: (i) either at once as a hot start, or with the spatiotemporal sequence derived from seismic inversion; and (ii) as a specified surface elevation or as a more realistic time and space-varying bottom boundary condition (in the latter case, we compute the initial tsunami generation up to 300 s using the non-hydrostatic model NHWAVE). Although additional refinements are expected in the near future, results based on the current FEM sources better explain long wave near-field observations at DART and GPS buoys near Japan, and measured tsunami inundation, while they simulate observations at distant DART buoys as well or better than the UCSB source. None of the sources, however, are able to explain the largest runup and inundation measured between 39.5° and 40.25°N, which could be due to insufficient model resolution in this region (Sanriku/Ria) of complex bathymetry/topography, and/or to additional tsunami generation mechanisms not represented in the coseismic sources (e.g., splay faults, submarine mass failure). This will be the object of future work.     

Tsunami threat in the Indian Ocean from a future megathrust earthquake west of Sumatra

Several independent indicators imply a high probability of a great (M > 8) earthquake rupture of the subduction megathrust under the Mentawai Islands of West Sumatra. The human consequences of such an event depend crucially on its tsunamigenic potential, which in turn depends on unpredictable details of slip distribution on the megathrust and how resulting seafloor movements and the propagating tsunami waves interact with bathymetry. Here we address the forward problem by modelling about 1000 possible complex earthquake ruptures and calculating the seafloor displacements and tsunami wave height distributions that would result from the most likely 100 or so, as judged by reference to paleogeodetic data. Additionally we carry out a systematic study of the importance of the location of maximum slip with respect to the morphology of the fore-arc complex. Our results indicate a generally smaller regional tsunami hazard than was realised in Aceh during the December 2004 event, though more than 20% of simulations result in tsunami wave heights of more than 5 m for the southern Sumatran cities of Padang and Bengkulu. The extreme events in these simulations produce results which are consistent with recent deterministic studies. The study confirms the sensitivity of predicted wave heights to the distribution of slip even for events with similar moment and reproduces Plafker's rule of thumb. Additionally we show that the maximum wave height observed at a single location scales with the magnitude though data for all magnitudes exhibit extreme variability. Finally, we show that for any coastal location in the near field of the earthquake, despite the complexity of the earthquake rupture simulations and the large range of magnitudes modelled, the timing of inundation is constant to first order and the maximum height of the modelled waves is directly proportional to the vertical coseismic displacement experienced at that point. These results may assist in developing tsunami preparedness strategies around the Indian Ocean and in particular along the coasts of western Sumatra.     

Field Survey of the 27 February 2010 Chile Tsunami

On 27 February 2010, a magnitude M _w?=?8.8 earthquake occurred off the coast of Chile??s Maule region causing substantial damage and loss of life. Ancestral tsunami knowledge from the 1960 event combined with education and evacuation exercises prompted most coastal residents to spontaneously evacuate after the earthquake. Many of the tsunami victims were tourists in coastal campgrounds. The international tsunami survey team (ITST) was deployed within days of the event and surveyed 800?km of coastline from Quintero to Mehuín and the Pacific Islands of Santa María, Mocha, Juan Fernández Archipelago, and Rapa Nui (Easter). The collected survey data include more than 400 tsunami flow depth, runup and coastal uplift measurements. The tsunami peaked with a localized runup of 29?m on a coastal bluff at Constitución. The observed runup distributions exhibit significant variations on local and regional scales. Observations from the 2010 and 1960 Chile tsunamis are compared.     

Tsunamis on the Pacific Coast of Canada Recorded in 1994–2007

In the last 15 years there have been 16 tsunami events recorded at tide stations on the Pacific Coast of Canada. Eleven of these events were from distant sources covering almost all regions of the Pacific, as well as the December 26, 2004 Sumatra tsunami in the Indian Ocean. Three tsunamis were generated by local or regional earthquakes and two were meteorological tsunamis. The earliest four events, which occurred in the period 1994–1996, were recorded on analogue recorders; these tsunami records were recently re-examined, digitized and thoroughly analysed. The other 12 tsunami events were recorded using digital high-quality instruments, with 1-min sampling interval, installed on the coast of British Columbia (B.C.) in 1998. All 16 tsunami events were recorded at Tofino on the outer B.C. coast, and some of the tsunamis were recorded at eight or more stations. The tide station at Tofino has been in operation for 100 years and these recent observations add to the dataset of tsunami events compiled previously by S.O. Wigen (1983) for the period 1906–1980. For each of the tsunami records statistical analysis was carried out to determine essential tsunami characteristics for all events (arrival times, maximum amplitudes, frequencies and wave-train structure). The analysis of the records indicated that significant background noise at Langara, a key northern B.C. Tsunami Warning station located near the northern end of the Queen Charlotte Islands, creates serious problems in detecting tsunami waves. That station has now been moved to a new location with better tsunami response. The number of tsunami events observed in the past 15 years also justified re-establishing a tide gauge at Port Alberni, where large tsunami wave amplitudes were measured in March 1964. The two meteorological events are the first ever recorded on the B.C. coast. Also, there have been landslide generated tsunami events which, although not recorded on any coastal tide gauges, demonstrate, along with the recent investigation of a historical catastrophic event, the significant risk that landslide generated tsunami pose to coastal and inland regions of B.C.     

Forecasting tsunamis in Poverty Bay,New Zealand,with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.

     

A Revised Evaluation of Tsunami Hazards along the Chinese Coast in View of the Tohoku-Oki Earthquake

Japan’s 2011 Tohoku-Oki earthquake and the accompanying tsunami have reminded us of the potential tsunami hazards from the Manila and Ryukyu trenches to the South China and East China Seas. Statistics of historical seismic records from nearly the last 4 decades have shown that major earthquakes do not necessarily agree with the local Gutenberg-Richter relationship. The probability of a mega-earthquake may be higher than we have previously estimated. Furthermore, we noted that the percentages of tsunami-associated earthquakes are much higher in major events, and the earthquakes with magnitudes equal to or greater than 8.8 have all triggered tsunamis in the past approximately 100 years. We will emphasize the importance of a thorough study of possible tsunami scenarios for hazard mitigation. We focus on several hypothetical earthquake-induced tsunamis caused by M _w 8.8 events along the Manila and Ryukyu trenches. We carried out numerical simulations based on shallow-water equations (SWE) to predict the tsunami dynamics in the South China and East China Seas. By analyzing the computed results we found that the height of the potential surge in China’s coastal area caused by earthquake-induced tsunamis may reach a couple of meters high. Our preliminary results show that tsunamis generated in the Manila and Ryukyu trenches could pose a significant threat to Chinese coastal cities such as Shanghai, Hong Kong and Macao. However, we did not find the highest tsunami wave at Taiwan, partially because it lies right on the extension of an assumed fault line. Furthermore, we put forward a multi-scale model with higher resolution, which enabled us to investigate the edge waves diffracted around Taiwan Island with a closer view.     

Biochemical characteristics of alga-bacterial mats and invertebrates from shallow-water hydrothermal fields of the West Pacific Ocean

In this study we assessed the biochemical adaptations of biotic specimens from the hydrothermal vent fields from different geographical zones. Algabacterial and bacterial mats and specimens of bivalves (9 species) and sea urchins (2 species) from different coastal shallow-water fields of hydrothermal activity along the West Rift zones of the Pacific Ocean were sampled to compare biochemical activities to volcanogenic chemical characteristics. The algabacterial and bacterial mat cells and subcellular fractions of bivalve gills and sea urchin guts were measured for 5-aminolevulinic acid (ALA), cytochromes b, (c + c₁), (a + a₃) and P450 (CYP), and activities of ALA synthetase and CYP-dependent aminopyrine N-demethylase (APND). The algabacterial mats from the Kuril Islands and New Zealand displayed comparable ALA and cytochrome levels and enzyme activities. Similar ALA synthesis, APND activity, and levels of cytochromes b and (c + c₁) were observed in bacterial mats from hydrothermal fields from Kuriles, New Britain Island and New Zealand. More active ALA synthesis, APND activity and enhanced content of cytochromes b and (c + c₁) and lower levels of CYP were found in bacterial mats than in algabacterial mats. Further, all bivalves influenced by hydrothermal discharges had elevated CYP levels and APND activity as well as a significant increase in ALA content and ALA synthesis (P < 0.05). Statistically sufficient (P < 0.05) alterations in ALA synthesis, APND activity and total CYP level were found in sea urchins under hydrothermal influence. Our findings contribute to the global assessment of hydrothermal effluents on biota and indicate that living conditions near hydrothermal vents accelerate metabolism of bivalves and sea urchins in examined hydrothermal zones in the northern and southern West Pacific Ocean.     

The landslides and tsunamis of the 30th of December 2002 in Stromboli analysed through numerical simulations

On the 30th of December 2002 two tsunamis were generated only 7 min apart in Stromboli, southern Tyrrhenian Sea, Italy. They represented the peak of a volcanic crisis that started 2 days before with a large emission of lava flows from a lateral vent that opened some hundreds of meters below the summit craters. Both tsunamis were produced by landslides that detached from the Sciara del Fuoco. This is a morphological scar and is the result of the last collapse of the northwestern flank of the volcanic edifice, that occurred less than 5 ka b.p. The first tsunami was due to a submarine mass movement that started very close to the coastline and that involved about 20×10⁶ m³ of material. The second tsunami was engendered by a subaerial landslide that detached at about 500 m above sea level and that involved a volume estimated at 4–9×10⁶ m³. The latter landslide can be seen as the retrogressive continuation of the first failure. The tsunamis were not perceived as distinct events by most people. They attacked all the coasts of Stromboli within a few minutes and arrived at the neighbouring island of Panarea, 20 km SSW of Stromboli, in less than 5 min. The tsunamis caused severe damage at Stromboli.In this work, the two tsunamis are studied by means of numerical simulations that use two distinct models, one for the landslides and one for the water waves. The motion of the sliding bodies is computed by means of a Lagrangian approach that partitions the mass into a set of blocks: we use both one-dimensional and two-dimensional schemes. The landslide model calculates the instantaneous rate of the vertical displacement of the sea surface caused by the motion of the underwater slide. This is included in the governing equations of the tsunami, which are solved by means of a finite-element (FE) technique. The tsunami is computed on two different grids formed by triangular elements, one covering the near-field around Stromboli and the other also including the island of Panarea.The simulations show that the main tsunamigenic potential of the slides is restricted to the first tens of seconds of their motion when they interact with the shallow-water coastal area, and that it diminishes drastically in deep water. The simulations explain how the tsunamis that are generated in the Sciara del Fuoco area, are able to attack the entire coastline of Stromboli with larger effects on the northern coast than on the southern. Strong refraction and bending of the tsunami fronts is due to the large near-shore bathymetric gradient, which is also responsible for the trapping of the waves and for the persistence of the oscillations. Further, the first tsunami produces large waves and runup heights comparable with the observations. The simulated second tsunami is only slightly smaller, though it was induced by a mass that is approximately one third of the first. The arrival of the first tsunami is negative, in accordance with most eyewitness reports. Conversely, the leading wave of the second tsunami is positive.     

Boulder deposition during major tsunami events

A remarkable accumulation of marine boulders located above the present spring tide level has occurred in two coastal lowlands of the Algarve (Portugal). The size‐interval of the particles studied here is seldom reported in the literature in association with extreme events of coastal inundation, thus making this study of relevance to many other coasts worldwide. The spreads of boulders extend several hundred meters inland and well beyond the present landward limit of storm activity. The marine origin of the boulders is demonstrated by well‐developed macro‐bioerosion sculpturing and in situ skeletal remains of endolithic shallow marine bivalves. The good state preservation of the fossils within the boulders indicates that abrasion during transport and redeposition was not significant. We envisage boulder deposition as having taken place during the Lisbon tsunami of ad 1755 through the simultaneous landward entrainment of coarse particles from nearshore followed by rapid shoreward suspended‐dominated transport and non‐graded redeposition that excluded significant sorting by weight or boulder dimensions. We use numerical hydrodynamic modeling of tsunami (and storm) waves to test the observational data on boulder dimensions (density, size, distribution) on the most likely processes of sediment deposition. This work demonstrates the effectiveness of the study of boulder deposits in tsunami reconstruction. Copyright © 2011 John Wiley & Sons, Ltd.     

Tsunami damage to coastal defences and buildings in the March 11th 2011 M w 9.0 Great East Japan earthquake and tsunami

On March 11th 2011 a M _w 9.0 mega-thrust interface subduction earthquake, the Great East Japan Earthquake, occurred 130 km off the northeast coast of Japan in the Pacific Ocean at the Japan Trench, triggering tsunami which caused damage along 600 km of coastline. Observations of damage to buildings (including vertical evacuation facilities) and coastal defences in Tōhoku are presented following investigation by the Earthquake Engineering Field Investigation Team (EEFIT) at 10 locations in Iwate and Miyagi Prefectures. Observations are presented in the context of the coastal setting and tsunami characteristics experienced at each location. Damage surveys were carried out in Kamaishi City and Kesennuma City using a damage scale for reinforced concrete (RC), timber and steel frame buildings adapted from an earlier EEFIT tsunami damage scale. Observations show that many sea walls and breakwaters were overtopped, overturned, or broken up, but provided some degree of protection. We show the extreme variability of damage in a local area due to inundation depth, flow direction, velocity variations and sheltering. Survival of many RC shear wall structures shows their high potential to withstand local earthquake and significant tsunami inundation but further research is required into mitigation of scour, liquefaction, debris impact, and the prevention of overturning failure. Damage to steel and timber buildings are also discussed. These observations are intended to contribute to mitigation of future earthquake and tsunami damage by highlighting the key features which influence damage level and local variability of damage sustained by urban coastal infrastructure when subjected to extreme tsunami inundation depths.     

Field survey report on tsunami disasters caused by the 1993 Southwest Hokkaido earthquake

Detailed field work at Okushiri Island and along the southwest coast of Hokkaido has revealed quantitatively (1) the advancing direction of tsunami on land, (2) the true tsunami height (i.e., height of tsunami, excluding its splashes, as measured from the ground) and (3) the flow velocity of tsunami on land, in heavily damaged areas. When a Japanese wooden house is swept away by tsunami, bolts that tie the house to its concrete foundation resist until the last moment and become bent towards the direction of the house being carried away. The orientations of more than 850 of those bent bolts and iron pipes (all that can be measured, mostly at Okushiri Island) and fell-down direction of about 400 trees clearly display how tsunami behaved on land and caused serious damage at various places. The true tsunami height was estimated by using several indicators, such as broken tree twigs and a window pane. The flow velocity of tsunami on land was determined by estimating the hydrodynamic force exerted on a bent handrail and a bent-down guardrail by the tsunami throughin situ strength tests.Contrary to the wide-spread recognition after the tsunami hazard, our results clearly indicate that only a few residential areas (i.e., Monai, eastern Hamatsumae, and a small portion at northern Aonae, all on Okushiri Island) were hit by a huge tsunami, with true heights reaching 10 m. Southern Aonae was completely swept away by tsunami that came directly from the focal region immediately to the west. The true tsunami height over the western sea wall of southern Aonae was estimated as 3 to 4 m. Northern Aonae also suffered severe damage due to tsunami that invaded from the corner zone of the sand dune (8 m high) and tide embankment at the northern end of the Aonae Harbor. This corner apparently acted as a tsunami amplifier, and tide embankment or breakwater can be quite dangerous when tsunami advances towards the corner it makes with the coast. The nearly complete devastation of Inaho at the northern end of Okushiri Island underscored the danger of tsunami whose propagation direction is parallel to the coast, since such tsunami waves tend to be amplified and tide embankment or breakwater is constructed low towards the coast at many harbors or fishing ports. Tsunami waves mostly of 2 to 4 m in true height swept away Hamatsumae on the southeast site of Okushiri Island where there were no coastal structures. Coastal structures were effective in reducing tsunami hazard at many sites. The maximum flow velocity at northern Aonae was estimated as 10 to 18 m/s (Tsutsumi et al., 1994), and such a high on-land velocity of tsunami near shore is probably due to the rapid shallowing of the deep sea near the epicentral region towards Okushiri Island. If the advancing direction, true height, and flow velocity of tsunami can be predicted by future analyses of tsunami generation and progagation, the analyses will be a powerful tool for future assessment of tsunami disasters, including the identification of blind spots in the tsunami hazard reduction.