Characterization of the 26 December 2004 tsunami deposits in Andaman Islands (Bay of Bengal, India)

The large tsunami, which was generated by an earthquake on 26 December 2004, affected most of the countries around the Indian Ocean. A total of 48 tsunamigenic surface sediments and nine core samples have been collected from various coastal geomorphological features such as beaches, estuaries/creeks and mangrove areas in the Andaman Islands. These samples were analysed for textural analysis and geochemical studies to evaluate effects of the tsunami on sediment contamination. The studied sediments, deposited by the 26 December 2004 tsunami in Andaman group of islands, belong to poorly sorted, coarse to medium sands. Generally the concentration of heavy metals in the tsunamigenic surface sediments is mainly in the order of Cu > Mn > Fe > Zn > Pb during the post-tsunami (2005) and Cu > Fe > Mn > Zn > Pb during the post-monsoon (2008). The analysed core samples show that tsunami sediments have been preserved at certain depths from the sampling locations and indicate that they were derived from shallow littoral to neritic depths. The approximate width of deposits deposited by the 26 December 2004 Tsunami in Diglipur and Mayabandar areas (North Andaman) is ~10 cm, in Rangat and Baratang (Middle Andaman) the thickness of the deposits is ~15 cm. In Chidiyatapu, Junglighat, Rutland Islands and Havelock Island (South Andaman) the thickness of the deposits is ~30, ~8, ~25 and ~5 cm, respectively, and in Hut Bay (Little Andaman) the thickness of the deposits is about ~15 cm.     

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

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

Modeling the 2004 Sumatra tsunami at Seychelles Islands: site-effect analysis and comparison with observations

The M _w = 9.1 mega-thrust Sumatra–Andaman earthquake that occurred on December 26, 2004, was followed by a devastating tsunami that killed hundreds of thousands of people and caused catastrophic effects on human settlements and environments along many coasts of the Indian Ocean, where even countries very far from the source were affected. One of these cases is represented by the Republic of Seychelles, where the tsunami reached the region about 7 h after the earthquake and produced relevant damages, despite the country was more than 4,500 km far from the seismic source. In the present work, we present and discuss a study of the 2004 Sumatra tsunami by means of numerical simulations with the attention focused on the effects observed at the Seychelles Archipelago, a region never previously investigated with this approach. The case is interesting since these islands lay on a very shallow oceanic platform with steep slopes so as the ocean depth changes from thousands to few tens of meters over short distances, with significant effects on the tsunami propagation features: the waves are strongly refracted by the oceanic platform and the tsunami signal is modified by the introduction of additional frequencies. The study is used also to validate the UBO-TSUFD numerical code on a real tsunami event in the far field, and the results are compared with the available observations, i.e., the sea level time series recorded at the Pointe La Rue station, Mahé Island, and run-up measurements and inundation lines surveyed few weeks after the tsunami at Praslin Island, where the tsunami hit during low tide. Synthetic results are found in good agreement with observations, even though some of the observations remain not fully solved. Moreover, simulations have been run in high-tide condition since the 2004 Sumatra tsunami hitting at high tide can be taken as the worst-case scenario for the Seychelles islands and used for tsunami hazard and risk assessments.     

Sheltered coastal environments as archives of paleo-tsunami deposits: Observations from the 2004 Indian Ocean tsunami

The 2004 earthquake left several traces of coseismic land deformation and tsunami deposits, both on the islands along the plate boundary and distant shores of the Indian Ocean rim countries. Researchers are now exploring these sites to develop a chronology of past events. Where the coastal regions are also inundated by storm surges, there is an additional challenge to discriminate between the deposits formed by these two processes. Paleo-tsunami research relies largely on finding deposits where preservation potential is high and storm surge origin can be excluded. During the past decade of our work along the Andaman and Nicobar Islands and the east coast of India, we have observed that the 2004 tsunami deposits are best preserved in lagoons, inland streams and also on elevated terraces. Chronological evidence for older events obtained from such sites is better correlated with those from Thailand, Sri Lanka and Indonesia, reiterating their usefulness in tsunami geology studies.     

Numerical simulation and field survey of inundation due to 2004 Indian Ocean tsunami in Trincomalee,Sri Lanka

This paper outlines the field measurements and numerical modelling carried out to develop a high-resolution tsunami inundation map, as a case study, for the city of Trincomalee on the east coast of Sri Lanka, which was devastated by the 2004 tsunami. We employ the deterministic approach together with numerical simulations based on the probable worst-case scenario to derive the inundation map. Linear and non-linear versions of shallow-water equations have been utilized to simulate tsunami propagation and onshore inundation, respectively. The field data considered in the present paper comprise the extent of inundation, the tsunami heights and the arrival times whilst the model results include the spatial distribution of the flow depth, the peak current speeds and the momentum flux. The computed extent of onshore inundation reproduces the observed overall pattern of inundation in most areas barring the south-eastern part of the city. Further, the model simulations suggest maximum flow depths up to about 2 m in most areas of the city whilst patches of flow depths exceeding 2 m can be seen in a narrow strip along the coastline. The computed current speeds also exceed 3 m/s at some locations adjacent to the shoreline.     

Tectonic geomorphology,tsunamis and environmental hazards: reference to Andaman-Nicobar Islands

Geomorphic features associated with earthquakes and tsunamis have received wide attention in estimating uplift and subsidence after the tectonic event. Although various techniques are in vogue in estimating the uplift and subsidence after the 2004 Andaman-Sumatra earthquake and subsequent tsunami, remote sensing techniques have been proved to be quite handy to study the geomorphic changes. In the present study, geomorphic changes associated with the destructive event of December 2004 have been analyzed. The emergent and subsident coasts around the smaller islands in the Andaman region have been identified. The coral reef area that has been subjected to uplift or subsidence in some of the islands of the Andaman and Nicobar region is delineated, and the net areal extents of these coral beds have been computed. Of the six islands studied in Andaman region, coral reef of four islands was subjected to uplift, and around two islands the area was subsided. The uplifted area varied from 0.10 to 11 km², and subsidence was about 0.50 km². In Nicobar region, the subsidence of coral reefs was recorded. This study helps to monitor the coastal environments and the destruction due to natural hazards.     

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

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

Socio-economic implications based on interviews with fishermen following the Indian Ocean tsunami

Following the devastating tsunami of 26 December 2004 in the Indian Ocean, there was a need to give a voice to the affected population. Hence a survey was conducted in the tsunami-affected regions of India. The tsunami mainly affected the states of Tamil Nadu, Kerala and Andhra Pradesh and the Union Territory of Pondicherry, all in south India, as well as the Andaman and Nicobar Islands of India in the Bay of Bengal. For various logistical reasons, no survey was conducted in the Andaman and Nicobar Islands. The survey was conducted during 21 January to 19 February 2005 and from 1 March to 8 March 2005. A total of eight people, arranged into four teams, simultaneously conducted the survey based upon a prepared questionnaire comprising a total of 16 questions. The total number of villages surveyed was 161, and the overall results of the survey are reported here. Among many observations, capacity building during the construction process, relocation and housing issues and tsunami education and awareness were prominent.     

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.     

Mega-Tsunami of 26th December, 2004： Indian initiative for early warning system and mitigation of oceanogenic hazards

The 26th December 2004 earthquake of Mw 9.3 is the second largest earthquake ever to have been recorded.This generated a tsunami which affected several Asian countries. In India, the Andaman & Nicobar group of islands, and coastal states of Tamil Nadu, Andhra Pradesh and Kerala were severely affected. Here, we briefly provide an outline of the approach taken by India for an early warning system for mitigation of oceanogenic disasters.     

Tsunamis versus storm deposits from Thailand

Along the Andaman (west) coast of Thailand, the 2004 tsunami depositional features associated with the 2004 tsunami were used to describe the characteristics of tsunamis in a place far away from the effect of both recent and ancient storms. The current challenge is that a lack of precise sedimentological characteristics have been described that will differentiate tsunami deposits from storm deposits. Here, in sedimentological senses, we reviewed the imprints of the sedimentological characteristics of the 2004 tsunami and older deposits and then compared them with storm deposits, as analyzed from the deposits found along the eastern (Gulf of Thailand; GOT) coast of Thailand. We discuss the hydraulic conditions of the 2004 tsunami and its predecessors, on the Andaman coast, and compare them to storm flows found on the coast of the GOT. Similar to an extensive tsunami inflow deposit, a storm flow overwash has very similar sedimentary structures. Well-preserved sedimentary structures recognized in sand sheets from both tsunami and storms include single and multiple normal gradings, reverse grading, parallel, incline and foreset lamina, rip-up clasts, and mud drapes. All these sedimentary structures verify the similarity of tsunami and storm inflow behavior as both types of high-energy flow start to scour the beach zone. Antidunes are likely to be the only unique internal sedimentary structures observed in the 2004 tsunami deposit. Rip-up clasts are rare within storm deposits compared to tsunami deposits. We found that the deposition during the outflow from both tsunami and storms was rarely preserved, suggesting that it does not persist for very long in the geological record.     

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.     

Method of calculating tsunami travel times in the Andaman Sea region

A new model to calculate tsunami travel times in the Andaman Sea region has been developed. The model specifically provides more accurate travel time estimates for tsunamis propagating to Patong Beach on the west coast of Phuket, Thailand. More generally, the model provides better understanding of the influence of the accuracy and resolution of bathymetry data on the accuracy of travel time calculations. The dynamic model is based on solitary wave theory, and a lookup function is used to perform bilinear interpolation of bathymetry along the ray trajectory. The model was calibrated and verified using data from an echosounder record, tsunami photographs, satellite altimetry records, and eyewitness accounts of the tsunami on 26 December 2004. Time differences for 12 representative targets in the Andaman Sea and the Indian Ocean regions were calculated. The model demonstrated satisfactory time differences (<2 min/h), despite the use of low resolution bathymetry (ETOPO2v2). To improve accuracy, the dynamics of wave elevation and a velocity correction term must be considered, particularly for calculations in the nearshore region.     

The post-depositional changes of the onshore 2004 tsunami deposits on the Andaman Sea coast of Thailand

The Indian Ocean tsunami flooded the coastal zone of the Andaman Sea and left tsunami deposits with a thickness of a few millimetres to tens of centimetres over a roughly one-kilometre-wide tsunami inundation zone. The preservation potential and the post-depositional changes of the onshore tsunami deposits in the coastal plain setting, under conditions of a tropical climate with high seasonal rainfall, were assessed by reinvestigating trenches located along 13 shore-perpendicular transects; the trenches were documented shortly after the tsunami and after 1, 2, 3 and 4 years. The tsunami deposits were found preserved after 4 years at only half of the studied sites. In about 30% of the sites, the tsunami deposits were not preserved due to human activity; in a further 20% of the sites, the thin tsunami deposits were eroded or not recognised due to new soil formation. The most significant changes took place during the first rainy season when the relief of the tsunami deposits was levelled; moderate sediment redeposition took place, and fine surface sediments were washed away, which frequently left a residual layer of coarse sand and gravel. The fast recovery of new plant cover stabilised the tsunami deposits and protected them against further remobilisation during the subsequent years. After five rainy seasons, tsunami deposits with a thickness of at least a few centimetres were relatively well preserved; however, their internal structures were often significantly blurred by roots and animal bioturbation. Moreover, soil formation within the deposits caused alterations, and in the case of thin layers, it was not possible to recognise them anymore. Tsunami boulders were only slightly weathered but not moved. Among the various factors influencing the preservation potential, the thickness of the original tsunami deposits is the most important. A comparison between the first post-tsunami survey and the preserved record suggests that tsunamis with a run-up smaller than three metres are not likely to be preserved; for larger tsunamis, only about 50% of their inundation area is likely to be presented by the preserved extent of the tsunami deposits. Any modelling of paleotsunamis from their deposits must take into account post-depositional changes.     

The assessment of health impact caused by energy use in urban areas of China: an intake fraction–based analysis

The Indian Ocean Tsunami of December 2004 caused inundation of seawater along the Northern coast of Tamil Nadu, India, resulting in loss of 8,000 people with extensive damage to properties. The paper describes the inundation of seawater in two northern districts, namely Kancheepuram and Villupuram districts, which showed distinct patterns of inundation of seawater and run-up levels due to variations in geomorphic features. TUNAMI N2 model was used to predict the seawater inundation for earthquakes occurred in 1881 at Car Nicobar, Sumatra 2004 and a worst-case scenario. The coastal areas with beaches having gentle slope showed more inundation compared with coastal areas having varied slope and habited by sand dunes and coastal vegetation. Appreciable inundation of seawater with tsunami simulated for 1881 Car Nicobar indicated that proximity to the source plays a major role besides earthquake parameters in causing inundation. The worst-case scenario generated from subduction zone of Car Nicobar using Sumatra 2004 earthquake parameters revealed extreme vulnerability of coasts of both the districts to giant tsunamis.     

A geological record of multiple Pleistocene tsunami inundations in an oceanic island: The case of Maio,Cape Verde

In the Central Atlantic archipelagos – the Canaries, Cape Verde, Madeira and the Azores – tsunami hazard is often regarded as low, when compared with other extreme wave events such as hurricanes and storms. The geological record of many of these islands, however, suggests that tsunami hazard may be underestimated, notwithstanding being lower than in areas adjacent to subduction zones, such as the margins of the Pacific and Indian oceans. Moreover, tsunamis in oceanic islands are generally triggered by local large-scale volcanic flank collapses, for which little is known about their frequency, making it difficult to estimate the probability of a new occurrence. Part of the problem lies in the fact that tsunami deposits are usually difficult to date, and few islands in the world exhibit evidence for repeated tsunami inundation on a protracted timescale. This study reports on the presence of abundant tsunami deposits (conglomerates and sandstones) on Maio Island (Cape Verde) and discusses their stratigraphy, sedimentological characteristics, probable age and tsunamigenic source. Observations indicate that four distinct inundation events of variable magnitude took place during the Pleistocene. One of the tsunami deposits yielded a high-confidence U/Th age of 78·8 ± 0·9 ka, which overlaps within error with the 73 ± 7 ka age proposed for Fogo volcano's flank collapse, an event known to have had a significant tsunami impact on nearby Santiago Island. This shows that the Fogo tsunami also impacted Maio, resulting in runups in excess of 60 m above coeval sea-level at ca 120 km from the source. Two older deposits, possibly linked to recurrent flank collapses of the Tope de Coroa volcano in Santo Antão Island, yielded lower-confidence ages of 479 to 390 ka and 360 to 304 ka. A younger deposit (<78 ka) remains undated. In summary, the geological record of Maio exhibits well-preserved evidence of repeated tsunami inundation, reinforcing the notion that tsunami hazard is not so low at volcanic archipelagos featuring prominent and highly-active volcanoes such as in Cape Verde.     

Chronology of major terrace forming events in the Andaman islands during the last 40 kyr

Major earthquakes that trigger tsunamis are great natural hazards. The devastations caused by the December 26, 2004 Sumatran earthquake, and the March 11, 2011 Japan earthquake, and associated tsunamis will remain in our memories for a long time. Such events reaffirm the need for studying the cause and effects of large earthquakes of the past and to prepare the world better for the future. In such an effort, to understand the pattern of earthquakes and their effects on the geomorphic evolution, we have studied deformation history in the Andaman and Nicobar Islands, located in one of the most active convergent margins of the world. Focusing on tectonically formed coastal terraces and determining the timing of their formation from the exposed dead corals, we have been able to reconstruct the history of major earthquakes in these islands for the last 40 kyr. Our results in conjunction with the existing radiocarbon age data from coastal terraces of these islands appear to suggest that the frequency of major earthquakes (M > 7) in the region has increased during the last 9 kyr. In confirmation with some earlier work, we find evidences for a major earthquake and a tsunami between 500-600 cal yr BP and possibly 4 others during 6–9 cal kyr BP. Our results also indicate that there has been a continuous subsidence of the south Andaman Islands.     

Tracing tsunami signatures of the ad 551 and ad 1303 tsunamis at the Gulf of Kyparissia (Peloponnese,Greece) using direct push in situ sensing techniques combined with geophysical studies

The western Peloponnese was repeatedly hit by major tsunami impacts during historical times as reported by historical accounts and recorded in earthquake and tsunami catalogues. Geological signatures of past tsunami impacts have also been found in many coastal geological archives. During the past years, abundant geomorphological and sedimentary evidence of repeated Holocene tsunami landfall was found between Cape Katakolo and the city of Kyparissia. Moreover, neotectonic studies revealed strong crust uplift along regional faults with amounts of uplift between 13 m and 30 m since the mid-Holocene. This study focuses on the potential of direct push in situ sensing techniques to detect tsunami sediments along the Gulf of Kyparissia. Direct push measurements were conducted on the landward shores of the Kaiafa Lagoon and the former Mouria Lagoon from which sedimentary and microfaunal evidence for tsunami landfall are already known. Direct push methods helped to decipher in situ high-resolution stratigraphic records of allochthonous sand sheets that are used to document different kinds of sedimentological and geomorphological characteristics of high-energy inundation, such as abrupt increases in grain size, integration of muddy rip-up clasts and fining upward sequences which are representative of different tsunami inundation pulses. These investigations were completed by sediment coring as a base for local calibration of geophysical direct push parameters. Surface-based electrical resistivity tomography and seismic data with highly resolved vertical direct push datasets and sediment core data were all coupled in order to improve the quality of the geophysical models. Details of this methodological approach, new in palaeotsunami research, are presented and discussed, especially with respect to the question of how the obtained results may help to facilitate tracing tsunami signatures in the sedimentary record and deciphering geomorphological characteristics of past tsunami inundation. Using direct push techniques and based on sedimentary data, sedimentary signatures of two young tsunami impacts that hit the Kaiafa Lagoon were detected. Radiocarbon age control allowed the identification of these tsunami layers as candidates for the ad 551 and ad 1303 earthquake and tsunami events. For these events, there is reliable historical data on major damage on infrastructure in western Greece and on the Peloponnese. At the former Mouria Lagoon, corroborating tsunami traces were found; however, in this case it is difficult to decide whether these signatures were caused by the ad 551 or the ad 1303 event.     

A 3-tier tsunami vulnerability assessment technique for the north-west coast of Peninsular Malaysia

The 2004 tsunami that struck the Sumatra coast gave a warning sign to Malaysia that it is no longer regarded as safe from a future tsunami attack. Since the event, the Malaysian Government has formulated its plan of action by developing an integrated tsunami vulnerability assessment technique to determine the vulnerability levels of each sector along the 520-km-long coastline of the north-west coast of Peninsular Malaysia. The scope of assessment is focused on the vulnerability of the physical characteristics of the coastal area, and the vulnerability of the built environment in the area that includes building structures and infrastructures. The assessment was conducted in three distinct stages which stretched across from a macro-scale assessment to several local-scale and finally a micro-scale assessment. On a macro-scale assessment, Tsunami Impact Classification Maps were constructed based on the results of the tsunami propagation modelling of the various tsunami source scenarios. At this stage, highly impacted areas were selected for an assessment of the local hazards in the form of local flood maps based on the inundation modelling output. Tsunami heights and flood depths obtained from these maps were then used to produce the Tsunami Physical Vulnerability Index (PVI) maps. These maps recognize sectors within the selected areas that are highly vulnerable to a maximum tsunami run-up and flood event. The final stage is the development of the Structural Vulnerability Index (SVI) maps, which may qualitatively and quantitatively capture the physical and economic resources that are in the tsunami inundation zone during the worst-case scenario event. The results of the assessment in the form of GIS-based Tsunami-prone Vulnerability Index (PVI and SVI) maps are able to differentiate between the various levels of vulnerability, based on the tsunami height and inundation, the various levels of impact severity towards existing building structures, property and land use, and also indicate the resources and human settlements within the study area. Most importantly, the maps could help planners to establish a zoning scheme for potential coastline development based on its sensitivity to tsunami. As a result, some recommendations on evacuation routes and tsunami shelters in the potentially affected areas were also proposed to the Government as a tool for relief agencies to plan for safe evacuation.