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

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

Tsunamis in the New Zealand archaeological record

Historical and geological records both indicate tsunami inundation of New Zealand in the 700 years since the first human settlement. In addition, Maori oral traditions refer to unusual waves that might have been tsunami waves, although the accounts are open to other interpretations. Tsunami evidence has rarely been proposed from archaeological sites, primarily because of a limited understanding of the requisite evidence and environmental context. We list a criteria suggesting possible tsunami inundation of archaeological sites based upon geoarchaeological data, and use them in a case study from the Archaic Maori occupation site at Wairau Bar. The list is possibly incomplete, but indicates that archaeological investigations can gain from assessments of changing environmental conditions through time at any individual site. Our intention is not to prove tsunami inundation; rather, it is to point to archaeological sites as possible sources of information. We highlight the potential of the Wairau Bar site for further investigation.     

Global tsunami deposits database

The National Geophysical Data Center and co-located World Data Center for Geophysics and Marine Geology provide integrated access to historical tsunami event, deposit, and proxy data. Historical events are important for understanding the frequency and intensity of relatively recent tsunamis. Deposit data collected during post-tsunami field surveys provide information on tsunami erosion, sedimentation, flow depths, inundation, and run-up. Deposit data from prehistoric tsunami events extend the record to pre-recorded times, constrain tsunami recurrence intervals, and estimate the minimum magnitude of tsunami inundation. Proxies indicate that an event capable of producing a tsunami occurred, but are not direct evidence of a tsunami. All of these data are used to develop tsunami hazard assessments, provide guidance to warning centers, validate models, inform community preparedness efforts, and educate the public about tsunami risks.     

Tsunami sediments and their foraminiferal assemblages

Tsunami hazard assessment begins with a compilation of past events that have affected a specific location. Given the inherent limitations of historical archives, the geological record has the potential to provide an independent dataset useful for establishing a richer, chronologically deeper time series of past events. Recent geological studies of tsunami are helping to improve our understanding of the nature and character of tsunami sediments. Wherever possible, geologists should be working to improve the research ‘tool kit’ available to identify past tsunami events. Marine foraminifera (single celled heterotrophic protists) have often been reported as present within tsunami-deposited sediments but in reality, little information about environmental conditions, and by analogy, the tsunami that deposited them, has been reported even though foraminifera have an enormous capacity to provide meaningful palaeo-environmental data. Here, we review what foraminifera are, describe their basic form and significance, summarise where they have been reported in tsunami sediments and identify what can be learnt from them. We review the gaps in our understanding and make recommendations to assist researchers who examine foraminiferal assemblages in order to enhance their use within tsunami geology.     

Numerical simulation of tsunamis — Its present and near future

Hindcasting of a tsunami by numerical simulations is a process of lengthy and complicated deductions, knowing only the final results such as run-up heights and tide records, both of which are possibly biased due to an insufficient number of records and due to hydraulic and mechanical limitation of tide gauges. There are many sources of error. The initial profile, determined with seismic data, can even be different from the actual tsunami profile. The numerical scheme introduces errors. Nonlinearity near and on land requires an appropriate selection of equations. Taking these facts into account, it should be noted that numerical simulations produce satisfactory information for practical use, because the final error is usually within 15% as far as the maximum run-up height is concerned.The state-of-the-art of tsunami numerical simulations is critically summarized from generation to run-up. Problems in the near future are also stated. Fruitful application of computer graphics is suggested.     

Physical criteria for distinguishing sandy tsunami and storm deposits using modern examples

Modern subaerial sand beds deposited by major tsunamis and hurricanes were compared at trench, transect, and sub-regional spatial scales to evaluate which attributes are most useful for distinguishing the two types of deposits. Physical criteria that may be diagnostic include: sediment composition, textures and grading, types and organization of stratification, thickness, geometry, and landscape conformity.

Published reports of Pacific Ocean tsunami impacts and our field observations suggest that sandy tsunami deposits are generally < 25 cm thick, extend hundreds of meters inland from the beach, and fill microtopography but generally conform to the antecedent landscape. They commonly are a single homogeneous bed that is normally graded overall, or that consists of only a few thin layers. Mud intraclasts and mud laminae within the deposit are strong evidence of tsunami deposition. Twig orientation or other indicators of return flow during bed aggradation are also diagnostic of tsunami deposits. Sandy storm deposits tend to be > 30 cm thick, generally extend < 300 m from the beach, and will not advance beyond the antecedent macrotopography they are able to fill. They typically are composed of numerous subhorizontal planar laminae organized into multiple laminasets that are normally or inversely graded, they do not contain internal mud laminae and rarely contain mud intraclasts. Application of these distinguishing characteristics depends on their preservation potential and any deposit modifications that accompany burial.

The distinctions between tsunami and storm deposits are related to differences in the hydrodynamics and sediment-sorting processes during transport. Tsunami deposition results from a few high-velocity, long-period waves that entrain sediment from the shoreface, beach, and landward erosion zone. Tsunamis can have flow depths greater than 10 m, transport sediment primarily in suspension, and distribute the load over a broad region where sediment falls out of suspension when flow decelerates. In contrast, storm inundation generally is gradual and prolonged, consisting of many waves that erode beaches and dunes with no significant overland return flow until after the main flooding. Storm flow depths are commonly < 3 m, sediment is transported primarily as bed load by traction, and the load is deposited within a zone relatively close to the beach.     

Impact of the 2004 Indian Ocean tsunami along the Tamil Nadu coastline: field survey review and numerical simulations

The 2004 Indian Ocean tsunami had a significant impact on the Tamil Nadu coast in India. In this paper, a range of field survey data collected by different survey teams available in literature have been analyzed and compiled to serve as a basis for validation of numerical tsunami simulations. The individual field surveys reveal a significant scatter in the run-up data between the different teams, which point out that the uncertainty in these data must be taken into account when using them for validation. The inundation of the 2004 Indian Ocean tsunami is simulated for the coastlines of Chennai and Nagapattinam based on high-resolution topography. Different spatially uniform Manning friction as well as heterogeneous friction maps is used. Overall, the simulation results showed a good agreement with the field observations, but there are also some observed spatial variability in the goodness of the fit between the data and simulations. In some areas, clear discrepancies are found. The results obtained using detailed land use maps including spatially variable friction are not significantly more accurate than those employing spatially constant values. For most areas, parameters indicating relatively low friction provided best match with the observations. This may also suggest that the inundation is often strongly governed by local variations in topography.     

Advancing tsunami risk assessment by improving spatio-temporal population exposure and evacuation modeling

Tsunamis are among the most destructive and lethal of coastal hazards. These are time-specific events, and despite directly affecting a narrow strip of coastline, a single occurrence can have devastating effects and cause massive loss of life, especially in urbanized coastal areas. In this work, in order to consider the time dependence of population exposure to tsunami threat, the variation of spatio-temporal population distribution in the daily cycle is mapped and analyzed in the Lisbon Metropolitan Area. High-resolution daytime and nighttime population distribution maps are developed using ‘intelligent dasymetric mapping,’ that is, applying areal interpolation to combine best-available census data and statistics with land use and land cover data. Workplace information and mobility statistics are considered for mapping daytime distribution. In combination with a tsunami hazard map, information on infrastructure, land use and terrain slope, the modeled population distribution is used to assess people’s evacuation speed, applying a geospatial evacuation modeling approach to the city of Lisbon. The detailed dynamic population exposure assessment allows producing both daytime and nighttime evacuation time maps, which provide valuable input for evacuation planning and management. Results show that a significant amount of population is at risk, and its numbers increase dramatically from nighttime to daytime, especially in the zones of high tsunami flooding susceptibility. Also, full evacuation can be problematic in the daytime period, even if initiated immediately after a major tsunami-triggering earthquake. The presented approach greatly improves tsunami risk assessment and can benefit all phases of the disaster management process.     

Analysis of the global tsunami data for vulnerability and risk assessment

For the assessment of tsunami risk and vulnerability, one has to make use of past tsunami observations. The most comprehensive tsunami databases for the world have been prepared by the National Geophysical Data Center of USA which are listed on their website for all the four oceans as well as the following marginal seas: Caribbean Sea, Mediterranean Sea, Black Sea, Red Sea and Gulf of Mexico. The dataset goes back as far as the first century AD and lists the events on a confidence rating scale of 0–4; 0 being an erroneous entry and 4 being a definite tsunami. Based on these various datasets for different geographical areas, a comprehensive global dataset was prepared in this study, which included only tsunami events with confidence rating of 3 and 4, meaning either probable or definite. In this composite and abridged global tsunami database there is no distinction either according to geography or tsunami strength as implied by its impact on the coast. A simple and straightforward statistical analysis suggests an almost complete randomness and no patterns that can be used for future tsunami predictions with a few minor exceptions.     

An assessment of educational tsunami evacuation map designs in Washington and Oregon

Educational tsunami evacuation map brochures in Washington and Oregon have been developed locally, resulting in significant differences between the types of tsunami hazard information they include. This paper identifies six tsunami hazard information types present in 38 brochures in Washington and Oregon: (1) tsunami hazard zone, (2) road network, (3) assembly areas, (4) evacuation guidance, (5) infrastructure, and (6) terrain. It compares and contrasts these information types in the maps and text of six of the brochures, including a proposed design standard in Oregon. Design differences of all 38 brochure maps are then organized using principles of cartographic abstraction, which describe mapmaker decisions about selection, generalization, and symbolization of information. We further use this framework to situate the information content of a new interactive Google Maps tool in Oregon. Our assessment identifies limitations of current tsunami hazard information that may be relevant to improving tsunami education. In theory, more advanced evacuation map tools can play an important role in reducing the limitations of tsunami hazard information relevant to the public. The new Google Maps tool addresses few of these limitations. Recognizing how map-making decisions define the underlying information content of evacuation maps can facilitate much needed future evaluations and developments in evacuation map design.     

Mapping of seawater inundation along Nagapattinam based on field observations

Nagapattinam, in the east coast of India, was severely affected during the deadliest Indian Ocean tsunami of December 26, 2004. The tsunami caused heavy damage to life and property, and the death toll was about 3,378 in Nagapattinam taluk. Certain villages along the coast witnessed large inundation while adjacent villages were protected from the fury of the tsunami waves. This study was carried out to examine the underlying causes for the vulnerability along Nagapattinam coast with the help of field observations, remote sensing, and geographical information system as tools. Coastal areas with high sand dunes have been protected from tsunami, and areas adjacent to backwaters were inundated. Realtime Kinematic Global Positioning System and high-resolution satellite data were used to map the topographic information and maximum extent of inundation. Thematic maps on land use, land cover, and coastal geomorphology were generated using remote sensing and field data. Using field data as the primary source of information, tsunami hazard maps have been generated for Nagapattinam.     

Distinguishing tsunami and storm deposits: An example from Martinhal, SW Portugal

Geological identification of past tsunamis is important for risk assessment studies, especially in areas where the historical record is limited or absent. The main problem when using the geological evidence is to distinguish between tsunami and storm deposits. Both are high-energy events that may leave marine traces in coastal stratigraphic sequences. At Martinhal, SW Portugal both storm surge and tsunami deposits are present at the same site within a single stratigraphic sequence, which makes it suitable to study the differences between them, excluding variations caused by local factors.

The tsunami associated with the Lisbon earthquake of November 1st 1755 AD, had a major impact on the geomorphology and sedimentology of Martinhal. It breached the barrier and laid down an extensive sheet of sand, as described in eyewitness reports. Besides the tsunami deposit the stratigraphy of Martinhal also displays evidence for storm surges that have breached and overtopped the barrier, flooding the lowland and leaving sand layers. Both marine-derived flood deposits show similar grain size characteristics and distinctive marine foraminifera. The most important differences are the rip-up clasts and boulders exclusively found in the tsunami deposit and the landward extent of the tsunami deposit that everywhere exceeds that of the storm deposits. Identification of both depositional units was only possible using a collection of different data and extensive stratigraphical information from cores as well as trenches.     

A method for estimating casualties due to the tsunami inundation flow

This study develops a method for estimating the number of casualties that may occur while people evacuate from an inundation zone when a tsunami has inundated an area. The method is based on a simple model of hydrodynamic forces as they affect the human body. The method uses a Tsunami casualty index (TCI) computed at each grid point of a numerical tsunami model to determine locations and times within the tsunami inundation zone where evacuation during the tsunami inundation is not possible and therefore where casualties are likely to occur. The locations and times can be combined with information about population density to compute the potential number of casualties. This information is useful in developing tsunami evacuation routes that avoid such locations. To illustrate the method, it is applied to the Seattle waterfront in Washington State, USA, that is under the threat of possible tsunami disasters due to Seattle Fault earthquakes. Preliminary results suggest that the tsunami casualties may occur within the Seattle waterfront for 15 min, during the time interval from 3 to 18 min after a large Seattle Fault tsunami is generated when the background tide level is mean high water.     

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.     

Tsunamigenic slope failures: the Pacific Islands ‘blind spot’?

We discuss issues related to a recognised shortcoming in existing tsunami hazard assessments for Pacific Island Countries and Territories (PICTs), that of tsunamigenic slope failures (TSFs). Currently, TSFs are most likely underrepresented as sources in existing tsunami databases for two key reasons. First, relatively low magnitude earthquakes associated with subduction zones are generally assigned as the tsunamigenic source, as opposed to the TSFs they generate. A reassessment of such ‘anomalous tsunamis’ may yield clues that serve to reassign their tsunamigenic source. Second, there are thousands of oceanic islands and seamounts scattered across the Pacific and flank collapse of volcanic edifices such as these is a largely unquantified tsunamigenic threat. However, while it is now possible to model such TSFs, this is unlikely to happen in the near future because of the lack of detailed bathymetry and landslide mass data. Recent developments in the identification of past tsunamis in the Pacific Islands have developed a unique range of indicators that can be used for identifying such events. These are geological, oral tradition and archaeological components that include, but are not limited to, a modified Darwinian model of atoll formation, coastal megaclasts, oral traditions of vanished islands and giant waves, and the abandonment of prehistoric coastal sites. As such, the most logical way forward is to use the multiple indicators available to us to identify evidence of past tsunamis.     

Tsunami hazard risk of a future volcanic eruption of Kolumbo submarine volcano,NE of Santorini Caldera,Greece

Kolumbo submarine volcano, located NE of Santorini caldera in the Aegean Sea, has only had one recorded eruption during historic times (1650 AD). Tsunamis from this event severely impacted the east coast of Santorini with extensive flooding and loss of buildings. Recent seismic studies in the area indicate a highly active region beneath Kolumbo suggesting the potential for future eruptive activity. Multibeam mapping and remotely operated vehicle explorations of Kolumbo have led to new insights into the eruptive processes of the 1650 AD eruption and improved assessments of the mechanisms by which tsunamis were generated and how they may be produced in future events. Principal mechanisms for tsunami generation at Kolumbo include shallow submarine explosions, entrance of pyroclastic flows into the sea, collapse of rapidly accumulated pyroclastic material, and intense eruption-related seismicity that may trigger submarine slope collapse. Compared with Santorini, the magnitude of explosive eruptions from Kolumbo is likely to be much smaller but the proximity of the volcano to the eastern coast of Santorini presents significant risks even for lower magnitude events.     

A simple model for calculating tsunami flow speed from tsunami deposits

This paper presents a simple model for tsunami sedimentation that can be applied to calculate tsunami flow speed from the thickness and grain size of a tsunami deposit (the inverse problem). For sandy tsunami deposits where grain size and thickness vary gradually in the direction of transport, tsunami sediment transport is modeled as a steady, spatially uniform process. The amount of sediment in suspension is assumed to be in equilibrium with the steady portion of the long period, slowing varying uprush portion of the tsunami. Spatial flow deceleration is assumed to be small and not to contribute significantly to the tsunami deposit. Tsunami deposits are formed from sediment settling from the water column when flow speeds on land go to zero everywhere at the time of maximum tsunami inundation. There is little erosion of the deposit by return flow because it is a slow flow and is concentrated in topographic lows. Variations in grain size of the deposit are found to have more effect on calculated tsunami flow speed than deposit thickness. The model is tested using field data collected at Arop, Papua New Guinea soon after the 1998 tsunami. Speed estimates of 14 m/s at 200 m inland from the shoreline compare favorably with those from a 1-D inundation model and from application of Bernoulli's principle to water levels on buildings left standing after the tsunami. As evidence that the model is applicable to some sandy tsunami deposits, the model reproduces the observed normal grading and vertical variation in sorting and skewness of a deposit formed by the 1998 tsunami.     

Historical tsunamis in the Adriatic Sea

Tsunami activity in the Adriatic Sea from the sixteenth century until the present has been analysed with the ultimate goal to improve the European tsunami catalogue and provide data for a new geo-database of tsunami events in the European-Mediterranean region. The study encompasses twenty-seven events, nine on the western and eighteen on the eastern coast of the Adriatic, with special attention being devoted to contemporary sources and to local journals and newspapers. For all the analysed events, the path of information from coeval sources, through the nineteenth century and up to modern tsunami catalogues, has been constructed. Tsunamis on the western coast have already been studied, but to obtain a coherent picture of tsunamigenic activity in the Adriatic Sea, they have been included in this work. Furthermore, the study was extended to see whether they had propagated to the opposite coast. Most of the events on the eastern coast have now been systematically analysed for the first time. The search of bibliographical sources revealed three new reports on tsunamis on the eastern coast that had not been previously recorded in international publications. The study established that, out of the eighteen eastern Adriatic events, twelve can be considered false, while six were true tsunamis. In the last 600 years, fifteen true tsunami events occurred in the Adriatic. One was very strong, six were strong or rather strong, and eight were light tsunamis. As a final result of this analysis, carried out according to standardised criteria, fifteen Adriatic tsunami events will be inserted in the TRANSFER (Tsunami Risk ANd Strategies For the European Region) database for the European-Mediterranean region.     

Impact of the National Tsunami Hazard Mitigation Program on Operations of the Richard H. Hagemeyer Pacific Tsunami Warning Center

The first 7 years of the National Tsunami Hazard Mitigation Program (NTHMP) have had a significant positive impact on operations of the Richard H. Hagemeyer Pacific Tsunami Warning Center (PTWC). As a result of its seismic project, the amount and quality of real-time seismic data flowing into PTWC has increased dramatically, enabling more rapid, accurate, and detailed analyses of seismic events with tsunamigenic potential. Its tsunameter project is now providing real-time tsunameter data from seven strategic locations in the deep ocean to more accurately measure tsunami waves as they propagate from likely source regions toward shorelines at risk. These data have already been used operationally to help evaluate potential tsunami threats. A new type of tsunami run-up gauge has been deployed in Hawaii to more rapidly assess local tsunamis. Lastly, numerical modeling of tsunamis done with support from the NTHMP is beginning to provide tools for real-time tsunami forecasting that should reduce the incidence of unnecessary warnings and provide more accurate forecasts for destructive tsunamis.     

Reassessing the tsunami risk in major ports and harbors of California I: San Diego

We assess tsunami hazards in San Diego Bay, California, using newly identified offshore tsunami sources and recently available high resolution bathymetric/topographic data. Using MOST (Titov and Synolakis, J Waterways Port Coastal Ocean Eng ASCE 124(4):57–171, 1998), we simulate locally, regionally and distant-generated tsunamis. Local tsunami source models use more realistic fault and landslide data than previous efforts. With the exception of the Alaska-Aleutian Trench, modeling results suggest that local sources are responsible for the largest waves within the San Diego Bay and Mission Bay. Because San Diego Bay is relatively well protected by North Island and the Silver Strand, the wave heights predicted are consistently smaller inside the harbor than outside. However, historical accounts, recent tsunamis and our predictions show that San Diego Bay is vulnerable to strong tsunami induced currents. More specifically, large currents are expected inside the harbor for various distant and local tsunami sources with estimated flow velocities exceeding 100 cm/s. Such currents have been damaging to harbor facilities, such as wharves and piers, and may cause boats to break from moorings and ram into adjacent harbor structures, as observed in recent historic tsunamis. More recently, following the M _w 8.8 February 27, 2010 Chile earthquake, tsunami-currents damaged docks/piers in Shelter Island confirming our findings. We note that the first generation of inundation maps in use in San Diego County by emergency management was based on much larger “worst case but realistic scenarios” (Synolakis et al. 2002a), which reflected the understanding of offshore hazards pervasive ten years ago. Large inundation and overland flow depths were observed primarily in local tsunami source simulations. In particular, locally induced tsunamis appear capable to overtop the Silver Strand. The results suggest that further work needs to be carried out with respect to local tsunami sources as they seem to have worse impact in the San Diego region than previously thought but probably low probability of occurrence. We also predict that a coastal community can be devastated simultaneously by large waves inundating shores and large currents in locations with small flow depths.