Contamination of tsunami sediments in a coastal zone inundated by the 26 December 2004 tsunami in Thailand

Tsunami sediments deposited in a coastal zone of Thailand by the 26 December 2004 tsunami wave were sampled within 50 days after the event. All surface and ground waters in tsunami- inundated zone revealed significant salinity at that time. The tsunami sediments, composed mainly of fine to medium sand, contain significantly elevated contents of salts (Na⁺, K⁺, Ca⁺², Mg⁺², Cl and SO ₄ ⁻² ) in water-soluble fraction, and of Cd, Cu, Zn, Pb in the bioavailable fraction and As in the exchangeable fraction in relation to the reference sample. The origin of contaminants is marine, as well as litho- and anthropogenic. The salts and Pb, Zn and Cu reveal high correlation to each other and to the mean grain size (pore size and porosity). Serious environmental hazard exists in that region because, due to gentle morphology, there is a risk of migration of the contaminants into ground waters and food chain.     

Use of tsunami waveforms for earthquake source study

Tsunami waveforms recorded on tide gauges, like seismic waves recorded on seismograms, can be used to study earthquake source processes. The tsunami propagation can be accurately evaluated, since bathymetry is much better known than seismic velocity structure in the Earth. Using waveform inversion techniques, we can estimate the spatial distribution of coseismic slip on the fault plane from tsunami waveforms. This method has been applied to several earthquakes around Japan. Two recent earthquakes, the 1968 Tokachi-oki and 1983 Japan Sea earthquakes, are examined for calibration purposes. Both events show nonuniform slip distributions very similar to those obtained from seismic wave analyses. The use of tsunami waveforms is more useful for the study of unusual or old earthquakes. The 1984 Torishima earthquake caused unusually large tsunamis for its earthquake size. Waveform modeling of this event shows that part of the abnormal size of this tsunami is due to the propagation effect along the shallow ridge system. For old earthquakes, many tide gauge records exist with quality comparable to modern records, while there are only a few good quality seismic records. The 1944 Tonankai and 1946 Nankaido earthquakes are examined as examples of old events, and slip distributions are obtained. Such estimates are possible only using tsunami records. Since tide-gauge records are available as far back as the 1850s, use of them will provide unique and important information on long-term global seismicity.     

Diatom biostratigraphy of tsunami deposits: Examples from the 1998 Papua New Guinea tsunami

Distinctive diatom assemblages may be associated with tsunami sediments and may often contrast with the assemblages found within sediments underlying the tsunami deposit as well as those associated with the modern coastal environment. Sediments associated with the 1998 tsunami that destroyed much of the Sissano lagoon area in northern Papua New Guinea have been investigated. Surface sediments from three transects across the sediment spit near Warapu have been examined for diatom content and preservation. The preservation is variable, and the data show an, often chaotic, assemblage that can be attributed to the tsunami waves incorporating and depositing diatoms from distinctive habitat zones during their runup and subsequent backwash. The diatoms identified within the Warapu sediments indicate an origin from within the inter-tidal and offshore area rather than from the beach–sand spit complex. The sand deposits disclose a high percentage, in excess of 75%, of broken diatom valves, and a predominance of centric (circular) species due to preferential preservation. The study demonstrates that the application of diatom biostratigraphy to modern tsunami deposits can be used in conjunction with other stratigraphical lines of evidence to interpret the source and provenance of historical and palaeo-tsunami deposits.     

Quality and durability assessments of the armourstones for two rubble mound breakwaters (Mersin,Turkey)

Due to economic reasons, natural stones (armourstones) of various sizes and qualities are very frequently used for the constructions of the breakwaters to protect coastal engineering structures from wave actions. Deterioration of the armourstones with time in the form of abrasion and disintegration may end up with the damage of the structures. Therefore, it is necessary to investigate the long-term performance and quality of the armourstones, which should be sound and durable. Mersin and Kumkuyu harbours were constructed using four different limestones obtained from two quarries. The limestones have different characters and site performances. In this study, the material and mass properties of the limestones taken from the quarries with known site performances as armourstones are investigated. The site performances and durability of the limestones are compared with the field measurements and laboratory works. Thus, the information obtained is used to assess long-term durability of the armourstones. The long-term performance of the Degirmencayi and Tirtar upper level limestones are observed to be good whereas it is rather poor for the Tirtar middle and lower level limestones. Comparison between the predicted and observed durabilities of the armourstones indicated that CIRIA/CUR, RDI_d, RERS, and wet to dry strength ratio give better results based on their field performances. However, the prediction of the durability of the limestones is poor in case RDIs, average pore diameter, and saturation coefficient are used.     

Hydraulic behavior of tsunami backflows: insights from their modern and ancient deposits

Tsunamis are unpredictable, catastrophic events, and so present enormous difficulties for direct studies in the field or laboratory. However, their sedimentary deposits yield evidence of a wide variety of hydrodynamic conditions caused by flow transformations on a spatial and temporal scale. Tsunami deposits ranging from the Miocene to modern times identified at different localities along the Chilean coast are described to provide a database of their characteristics. Among the typical features associated with tsunami deposits are well-rounded megaclasts eroded from coastal alluvial fans or beaches by very dense, competent flows. Sand injections from the base of these flows into the substrate indicate very high dynamic pressures, whereas basal shear carpets suggest hyperconcentrated, highly sheared flows. Turbulence develops in front of advancing debris flows, as indicated by megaflutes at the base of scoured channels.     

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.     

Identification of tsunami deposits considering the tsunami waveform: An example of subaqueous tsunami deposits in Holocene shallow bay on southern Boso Peninsula, Central Japan

This study proposes a tsunami depositional model based on observations of emerged Holocene tsunami deposits in outcrops located in eastern Japan. The model is also applicable to the identification of other deposits, such as those laid down by storms. The tsunami deposits described were formed in a small bay of 10–20-m water depth, and are mainly composed of sand and gravel. They show various sedimentary structures, including hummocky cross-stratification (HCS) and inverse and normal grading. Although, individually, the sedimentary structures are similar to those commonly found in storm deposits, the combination of vertical stacking in the tsunami deposits makes a unique pattern. This vertical stacking of internal structures is due to the waveform of the source tsunamis, reflecting: 1) extremely long wavelengths and wave period, and 2) temporal changes of wave sizes from the beginning to end of the tsunamis.

The tsunami deposits display many sub-layers with scoured and graded structures. Each sub-layer, especially in sandy facies, is characterized by HCS and inverse and normal grading that are the result of deposition from prolonged high-energy sediment flows. The vertical stack of sub-layers shows incremental deposition from the repeated sediment flows. Mud drapes cover the sub-layers and indicate the existence of flow-velocity stagnant stages between each sediment flow. Current reversals within the sub-layers indicate the repeated occurrence of the up- and return-flows.

The tsunami deposits are vertically divided into four depositional units, Tna to Tnd in ascending order, reflecting the temporal change of wave sizes in the tsunami wave trains. Unit Tna is relatively fine-grained and indicative of small tsunami waves during the early stage of the tsunami. Unit Tnb is a protruding coarse-grained and thickest-stratified division and is the result of a relatively large wave group during the middle stage of the tsunami. Unit Tnc is a fine alternation of thin sand sheets and mud drapes, deposited from waning waves during the later stage of the tsunami. Unit Tnd is deposited during the final stage of the tsunami and is composed mainly of suspension fallout. Cyclic build up of these sub-layers and depositional units cannot be explained by storm waves with short wave periods of several to ten seconds common in small bays.     

Exact analytical solutions of nonlinear problems of tsunami wave run-up on slopes with different profiles

A review of papers investigating tsunami wave run-up on a beach is given and the control parameters of the problem are revealed. There are two such parameters in the case of ideal fluid: the bottom sloping angle and the breaking parameter. A stage-by-stage approach for finding run-up characteristics is formulated: the linear calculation of shoreline oscillations and the subsequent non-linear transformation of the solution according to the Riemann method. Solution of the nononedimensional problems of wave run-up on a beach in the linear formulation is obtained.     

Effects of rainy season on mobilization of contaminants from tsunami deposits left in a coastal zone of Thailand by the 26 December 2004 tsunami

Study on contamination of tsunami sediments deposited on 26 December 2004 conducted shortly after the tsunami in coastal zone of Thailand revealed elevated contents of salts in water-soluble and some heavy metals and arsenic in bioavailable fractions (Szczuciński et al. in Env Geol 49:321–331, 2005). Few months later rainy season started and effected in total rainfall of over 3,300 mm. This paper presents results of survey repeated 1 year after the tsunami. To assess the effects of rainy season on mobilization of previously determined potential contaminants, the same locations were sampled again and analysed with the same methods. The tsunami deposit layer was well preserved but in many locations the sediments were coarser than just after tsunami due to washing out of finer fractions. The water-soluble salts contents were strongly reduced after the rainy season. However, the concentrations of acid leachable heavy metals and metalloids were still elevated in comparison to reference sample from an area not impacted by tsunami. It is possible that the metals and metalloids are successively moved to more bioavailable fraction from forms which were more resistant to mobilization.     

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.     

Three-dimensional slope reliability and risk assessment using auxiliary random finite element method

This paper aims to propose an auxiliary random finite element method (ARFEM) for efficient three-dimensional (3-D) slope reliability analysis and risk assessment considering spatial variability of soil properties. The ARFEM mainly consists of two steps: (1) preliminary analysis using a relatively coarse finite-element model and Subset Simulation, and (2) target analysis using a detailed finite-element model and response conditioning method. The 3-D spatial variability of soil properties is explicitly modeled using the expansion optimal linear estimation approach. A 3-D soil slope example is presented to demonstrate the validity of ARFEM. Finally, a sensitivity study is carried out to explore the effect of horizontal spatial variability. The results indicate that the proposed ARFEM not only provides reasonably accurate estimates of slope failure probability and risk, but also significantly reduces the computational effort at small probability levels. 3-D slope probabilistic analysis (including both 3-D slope stability analysis and 3-D spatial variability modeling) can reflect slope failure mechanism more realistically in terms of the shape, location and length of slip surface. Horizontal spatial variability can significantly influence the failure mode, reliability and risk of 3-D slopes, especially for long slopes with relatively strong horizontal spatial variability. These effects can be properly incorporated into 3-D slope reliability analysis and risk assessment using ARFEM.     

Hazard assessment of the Tidal Inlet landslide and potential subsequent tsunami, Glacier Bay National Park, Alaska

An unstable rock slump, estimated at 5 to 10 × 10⁶ m³, lies perched above the northern shore of Tidal Inlet in Glacier Bay National Park, Alaska. This landslide mass has the potential to rapidly move into Tidal Inlet and generate large, long-period-impulse tsunami waves. Field and photographic examination revealed that the landslide moved between 1892 and 1919 after the retreat of the Little Ice Age glaciers from Tidal Inlet in 1890. Global positioning system measurements over a 2-year period show that the perched mass is presently moving at 3–4 cm annually indicating the landslide remains unstable. Numerical simulations of landslide-generated waves suggest that in the western arm of Glacier Bay, wave amplitudes would be greatest near the mouth of Tidal Inlet and slightly decrease with water depth according to Green’s law. As a function of time, wave amplitude would be greatest within approximately 40 min of the landslide entering water, with significant wave activity continuing for potentially several hours.     

The tsunami of 426 BC in the Maliakos gulf, eastern Greece

In 427 BC, a major earthquake occurred in ancient Greece. In particular, Attica, Boeotia, and the island of Euboea were the areas where seismic activity was most frequent. The fact that these events happened in conjuction with the Peloponnesian war provides us with an account made by historians of the war. Such an account is the one made by Thucydides.During the spring and summer of 426 BC shocks continued to take place. This time, the sea area between the island of Euboea and the mainland (Maliakos gulf) was also affected and as a result, a seismic sea-wave of considerable size formed. The tsunami, as it is better known, swept the surrounding coastal area. Major topographic alteration of the area occurred, resulting in a huge loss of life and the destruction of cities.In this paper, the author attempts to describe this event and to explain scientifically how it happened, and how this affected the shape of the area and human life.All the evidence used in this paper has been cross-referenced with at least one other historic or scientific source. Although it was extremely difficult to uncover hidden detail about an event so far in the past, any facts that could not be verified have not been included.     

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

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

Estimation of tsunami hazard from volcanic activity — Suggested methodology with Augustine volcano,Alaska as an example

A general approach for the estimation of tsunami height and hazard in the vicinity of active volcanoes has been developed. An empirical relationship has been developed to estimate the height of the tsunami generated for an eruption of a given size. This relationship can be used to estimate the tsunami hazard based on the frequency of eruptive activity of a particular volcano. This technique is then applied to the estimation of tsunami hazard from the eruption of the Augustine volcano in Alaska. Modification of this approach to account for a less than satisfactory data base and differing volcanic characteristics is also discussed with the case of the Augustine volcano as an example. This approach can be used elsewhere with only slight modifications and, for the first time, provides a technique to estimate tsunami hazard from volcanic activity, similar to a well-established approach for the estimation of tsunami hazard from earthquake activity.     

Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

Damage formation associated with bending under a constant moment

Deformation within the Earth's lithosphere is largely controlled by the rheology of the rock. Fracture and faulting are characterized by elastic rheologies with brittle mechanisms, while folding and flow are characterized by plastic and/or viscous rheologies due to ductile mechanisms. However, it has been recognized that deformation that resembles ductile behavior can be produced within the confines of the brittle lithosphere. Specific examples are folds that form in the shallow crust, steep hinges at subduction zones that are accompanied by seismicity, and large-scale deformation at plate boundaries. In these cases, the brittle lithosphere behaves elastically with fracture and faulting yet produces ductile behavior. In this paper, we attempt to simulate such ductile behavior in elastic materials using continuum damage mechanics. Engineers utilize damage mechanics to model the continuum deformation of brittle materials. We utilize a modified form of damage mechanics that represents a reduction in frictional strength of preexisting fractures and faults. We use this empirical approach to simulate the bending of the lithosphere under the application of a constant moment.We use numerical simulations to obtain elastostatic solutions for plate bending and where the longitudinal stress at a particular node exceeds a yield stress, we apply damage to reduce Young's modulus at the node. Damage is calculated at each time step by a power-law relationship of the ratio of the yield stress to the longitudinal stress and the yield strain to the longitudinal strain. This results in the relaxation of the material due to increasing damage. To test our method, we apply our damage rheology to an infinite plate deforming under a constant bending moment. We simulate a wide range of behaviors from slow relaxation to instantaneous failure, over timescales that span six orders of magnitude. Using this method, stress relaxation produces elastic-perfectly plastic behavior in cases where failure does not occur. For cases of failure, we observe a rapid increase in damage leading to failure, analogous to the acceleration of microcrack formation and acoustic emissions prior to failure. The changes in the rate of damage accumulation in failure cases are similar to the changes in b-values of acoustic emissions observed in triaxial compression tests of fractured rock and b-value changes prior to some large earthquakes. Thus continuum damage mechanics can simulate the phenomenon of ductile behavior due to brittle mechanisms as well as observations of laboratory experiments and seismicity.     

Investigation of long waves in the tsunami frequency band on the southwestern shelf of Kamchatka

Two inexpensive cable bottom pressure stations were installed on the southwestern shelf of Kamchatka (Okhotsk Sea) in 1987 and two more in 1988 to provide longwave measurements in the tsunami frequency band, to investigate the generating mechanism of these waves, and to test the instrumentation. Microfluctuations of atmospheric pressure were recorded simultaneously. Two cable lines were torn off by ship anchors in March 1989 but others are still working in spite of highly dynamic activity on beaches and in hard ice regimes. Careful data analysis of two months of observations (September–October, 1987) showed that: (1) the atmospheric spectra were very stable and monotonic in the period range 2–50 min and corresponded to a power law of ^–2.3, (2) the direct generation of long waves by atmospheric pressure fluctuations was negligible, (3) there was high correlation between background longwave oscillations and sea state, (4) the structure of the offshore longwave field was in good agreement with theoretical estimates of standing waves for a linear slope.     

Boulders of MIS 5 age deposited by a tsunami on the coast of Otago, New Zealand

A series of elevated imbricated boulders were investigated on the Otago coastline, southeast New Zealand, through field surveying and optical luminescence dating. By using established hydrodynamic relationships of sediment transport the energy required to move the clasts was calculated and compared to the historic record of marine inundations of that coast. The boulders are platy in shape and are over 2 m long in some cases, and are sourced from a locally outcropping conglomerate unit which appears to be the only lithology on this section of coast that erodes to produce clasts of this size. It is estimated that the boulders were deposited by a tsunami between 2 and 3 m high during the latter part of Marine Isotope Stage 5. They therefore represent the first pre-Holocene tsunami deposit and one composed of large boulders described on the New Zealand coastline.     

Evaluation of geochemical impact of tsunami on Pichavaram mangrove ecosystem,southeast coast of India

The 26 December 2004-tsunami has deposited sediments in the Pichavaram mangrove ecosystem, east coast of India. Ten surface and three core sediment samples were collected within 30 days of the event and analyzed for nutrients. Water samples were also analyzed to see the impact of tsunami on the geochemical behavior of nutrients. An increase in the concentration of various nutrients namely nitrate and phosphate was observed. The geochemistry of the mangrove forest was observed to be influenced by a number of factors like rapid increase of aquaculture farms, agricultural practices and the anthropogenic discharge from the nearby-inhabited areas. Further the sediment column was disturbed due to energetic tsunami waves, which has caused a sheer increase in the dissolved oxygen in water. As a result, the change in the redox potential has resulted in change in the nutrients absorbed/associated with the sediments. In addition, role of retreating water after tsunami on the nutrient geochemistry was also evaluated.