Re-examination of the Source Mechanism of the 1998 Papua New Guinea Earthquake and Tsunami

— Simulation of tsunami propagation and runup of the 1998 Papua New Guinea (PNG) earthquake tsunami using the detailed bathymetry measured by JAMSTEC and adding bathymetric data at depths less than 60 m is carried out, reproducing the tsunami energy focus into Warapu and Arop along the Sissano Lagoon. However, the computed runup heights in the lagoon are still lower than those measured. Even if the error in estimating the fault parameters is taken into consideration, computational results are similar. Analysis by the wave ray method using several scenarios of the source size of the tsunami and location by the wave ray method suggests that a source characterized by small size in water 1,000-m deep approximately 25 km offshore the lagoon, best fits the arrival determined from the interviews with eyewitnesses. A two-layer numerical model simulating the interaction of the tsunami with a landslide is employed to study the behavior of a landslide-generated tsunami with different size sand depths of the initial slide just outside the lagoon. A landslide model with a volume of 4–8 × 10⁹ m³ is selected as the best in order to reproduce the distribution of the measured tsunami runup in the lagoon. The simulation of a tsunami generated in two stages, fault and landslide, could show good agreement with the runup heights and distribution of the arrival time, but a time gap of around 10 minutes remains, suggesting that a tsunami generated by the mainshock at 6:49 PM local time is too small for people to notice, and the following tsunami triggered by landslide or mass movement near the lagoon about ten minutes after the mainshock attacked the coast and caused the huge damage.     

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.     

Computation of Tsunamis in the East (Japan) Sea Using Dynamically Interfaced Nested Model

-- A study of 1983 and 1993 tsunami events in the East (Japan) Sea using multigrid dynamically-interfaced finite difference models was performed to produce detailed features of coastal inundations along the eastern coast of Korea. The computational domain is composed of several subregions with different grid sizes connected in parallel or inclined directions. The innermost subregion represents the coastal alignment reasonably well and has a grid size of about 30 meters. Numerical simulations have been performed in the framework of shallow-water equations (linear, as well as nonlinear) over the plane or spherical coordinate system, depending on the dimensions of the subregion. Results of simulations show the general agreements with the observed data of the runup height for both tsunamis. The evolution of the distribution function of tsunami heights is studied numerically and it shows a tendency to the log-normal curve for a long distance from the source.     

Tsunami Source of the 2010 Mentawai, Indonesia Earthquake Inferred from Tsunami Field Survey and Waveform Modeling

The 2010 Mentawai earthquake (magnitude 7.7) generated a destructive tsunami that caused more than 500 casualties in the Mentawai Islands, west of Sumatra, Indonesia. Seismological analyses indicate that this earthquake was an unusual “tsunami earthquake,” which produces much larger tsunamis than expected from the seismic magnitude. We carried out a field survey to measure tsunami heights and inundation distances, an inversion of tsunami waveforms to estimate the slip distribution on the fault, and inundation modeling to compare the measured and simulated tsunami heights. The measured tsunami heights at eight locations on the west coasts of North and South Pagai Island ranged from 2.5 to 9.3 m, but were mostly in the 4–7 m range. At three villages, the tsunami inundation extended more than 300 m. Interviews of local residents indicated that the earthquake ground shaking was less intense than during previous large earthquakes and did not cause any damage. Inversion of tsunami waveforms recorded at nine coastal tide gauges, a nearby GPS buoy, and a DART station indicated a large slip (maximum 6.1 m) on a shallower part of the fault near the trench axis, a distribution similar to other tsunami earthquakes. The total seismic moment estimated from tsunami waveform inversion was 1.0 × 10²¹ Nm, which corresponded to M_w 7.9. Computed coastal tsunami heights from this tsunami source model using linear equations are similar to the measured tsunami heights. The inundation heights computed by using detailed bathymetry and topography data and nonlinear equations including inundation were smaller than the measured ones. This may have been partly due to the limited resolution and accuracy of publically available bathymetry and topography data. One-dimensional run-up computations using our surveyed topography profiles showed that the computed heights were roughly similar to the measured ones.     

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.     

Runup Characteristics of Symmetrical Solitary Tsunami Waves of “Unknown” Shapes

The problem of tsunami wave runup on a beach is discussed in the framework of the rigorous solutions of the nonlinear shallow-water theory. We present an analysis of the runup characteristics for various shapes of the incoming symmetrical solitary tsunami waves. It will be demonstrated that the extreme (maximal) wave characteristics on a beach (runup and draw-down heights, runup and draw-down velocities and breaking parameter) are weakly dependent on the shape of incident wave if the definition of the “significant” wavelength determined on the 2/3 level of the maximum height is used. The universal analytical expressions for the extreme wave characteristics are derived for the runup of the solitary pulses. They can be directly applicable for tsunami warning because in many cases the shape of the incident tsunami wave is unknown.     

Tsunami generation of the 1993 Hokkaido Nansei-Oki earthquake

Heterogeneous fault motion of the 1993 Hokkaido Nansei-Oki earthquake is studied by using seismic, geodetic and tsunami data, and the tsunami generation from the fault model is examined. Seismological analyses indicate that the focal mechanism of the first 10 s, when about a third of the total moment was released, is different from the overall focal mechanism. A joint inversion of geodetic data on Okushiri Island and the tide gauge records in Japan and Korea indicates that the largest slip, about 6 m, occurred in a small area just south of the epicenter. This corresponds to the initial rupture on a fault plane dipping shallowly to the west. The slip on the northernmost subfault, which is dipping to the east, is about 2 m, while the slips on the southern subfaults, which are steeply dipping to the west, are more than 3 m. Tsunami heights around Okushiri Island are calculated from the heterogeneous fault model using different grid sizes. Computation on the smaller grids produces large tsunami height that are closer to the observed tsunami runup heights. Tsunami propagation in the nearly closed Japan Sea is examined as the free oscillation of the Japan Sea. The excitation of the free oscillation by this earthquake is smaller than that by the 1964 Niigata or 1983 Japan Sea earthquake.     

Combined Effects of Tectonic and Landslide-Generated Tsunami Runup at Seward, Alaska During the M W 9.2 1964 Earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the $M_{\rm W}$ 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5?min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211?million m³ (Haeussler et?al. in Submarine mass movements and their consequences, pp 269?C278, 2007). The first tectonic tsunami wave arrived in Resurrection Bay about 30?min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et?al. in Pure Appl Geophys 166:131?C152, 2009). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559?C572, 1994) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.     

Source Fault Model of the 1771 Yaeyama Tsunami, Southern Ryukyu Islands, Japan, Inferred from Numerical Simulation

The 1771 Yaeyama tsunami is successfully reproduced using a simple faulting model without submarine landslide. The Yaeyama tsunami (M 7.4), which struck the southern Ryukyu Islands of Japan, produced unusually high tsunami amplitudes on the southeastern coast of Ishigaki Island and caused significant damage, including 12,000 casualties. Previous tsunami source models for this event have included both seismological faults and submarine landslides. However, no evidence of landslides in the source has been obtained, despite marine surveying of the area. The seismological fault model proposed in this study, describing a fault to the east of Ishigaki Island, successfully reproduces the distribution of tsunami runup on the southern coast of the Ryukyu Islands. The unusual runup heights are found through the numerical simulation attributable to a concentration of tsunami energy toward the southeastern coast of Ishigaki Island by the effect of the shelf to the east. Thus, the unusual runup heights observed on the southeastern coast of Ishigaki Island can be adequately explained by a seismological fault model with wave-ray bending on the adjacent shelf.     

Satellite Data for a Rapid Assessment of Tsunami Inundation Areas after the 2011 Tohoku Tsunami

The M _w = 9.0 earthquake that occurred off the coast of Japan’s Tohoku region produced a great tsunami causing catastrophic damage and loss of life. Within hours of the tsunami event, satellite data were readily available and massive media coverage immediately circulated thousands of photographs and videos of the tsunami. Satellite data allow a rapid assessment of inundated areas where access can be difficult either as a result of damaged infrastructure (e.g., roads, bridges, ports, airports) or because of safety issues (e.g., the hazard at Nuclear Power Plant at Fukushima). In this study, we assessed in a day tsunami inundation distances and runup heights using satellite data (very high-resolution satellite images from the GeoEye1 satellite and from the DigitalGlobe worldview, SRTM and ASTER GDEM) of the Tohoku region, Northeast Japan. Field survey data by Japanese and other international scientists validated our results. This study focused on three different locations. Site selection was based on coastal morphologies and the distance to the tsunami source (epicenter). Study sites are Rikuzentakata, Oyagawahama, and Yagawahama in the Oshika Peninsula, and the Sendai coastal plain (Sendai City to Yamamoto City). Maximum inundation distance (6 km along the river) and maximum runup (39 m) at Rikuzentakata estimated from satellite data agree closely with the 39.7 m inundation reported in the field. Here the ria coastal morphology and horn shaped bay enhanced the tsunami runup and effects. The Sendai coastal plain shows large inundation distances (6 km) and lower runup heights. Natori City and Wakabayashi Ward, on the Sendai plain, have similar runup values (12 and 16 m, respectively) obtained from SRTM data; these are comparable to those obtained from field surveys (12 and 9.5 m). However, at Yagawahama and Oyagawahama, Miyagi Prefecture, both SRTM and ASTER data provided maximum runup heights (41 to 45 m and 33 to 34 m, respectively), which are higher than those measured in the field (about 27 m). This difference in DEM and field data is associated with ASTER and SRTM DEM’s pixel size and vertical accuracy, the latter being dependent on ground coverage, slope, aspect and elevation. Countries with less access to technology and infrastructure can benefit from the use of satellite imagery and freely available DEMs for an initial, pre-field surveys, rapid estimate of inundated areas, distances and runup, and for assisting in hazard management and mitigation after a natural disaster.     

Field Measurements and Numerical Simulations of the 2004 Tsunami Impact on the East Coast of Sri Lanka

This paper employs a numerical model of tsunami propagation together with documented observations and field measurements of the evidence left behind by the tsunami in December 2004, to identify and interpret the factors that have contributed to the significant spatial variability of the level of tsunami impact along the coastal belt of the eastern province of Sri Lanka. The model results considered in the present analysis include the distribution of the amplitude of the tsunami and the pattern of wave propagation over the continental shelf off the east coast, while the field data examined comprise the maximum water levels measured at or near the shoreline, the horizontal inundation distances and the number of housing and other buildings damaged. The computed maximum amplitude of the tsunami at water points nearest the shoreline along the east coast shows considerable variation ranging from 2.2 m to 11.4 m with a mean value of 5.7 m; moreover, the computed amplitudes agree well with the available field measurements. We also show that the shelf bathymetry off the east coast, particularly the submarine canyons at several locations, significantly influences the near-shore transformation of tsunami waves, and consequently, the spatial variation of the maximum water levels along the coastline. The measured values of inundation also show significant variation along the east coast and range from 70 m to 4560 m with a median value of 700 m. Our analyses of field data also show the dominant influence of the coastal topography and geomorphology on the extent of tsunami inundation. Furthermore, the measured inundation distances indicate no apparent correlation with the computed tsunami heights at the respective locations. We also show that both the computed tsunami heights and the measured inundation distances for the east coast closely follow the log-normal statistical distribution.     

Runup of Tsunami Waves in U-Shaped Bays

The problem of tsunami wave shoaling and runup in U-shaped bays (such as fjords) and underwater canyons is studied in the framework of 1D shallow water theory with the use of an assumption of the uniform current on the cross-section. The wave shoaling in bays, when the depth varies smoothly along the channel axis, is studied with the use of asymptotic approach. In this case a weak reflection provides significant shoaling effects. The existence of traveling (progressive) waves, propagating in bays, when the water depth changes significantly along the channel axis, is studied within rigorous solutions of the shallow water theory. It is shown that traveling waves do exist for certain bay bathymetry configurations and may propagate over large distances without reflection. The tsunami runup in such bays is significantly larger than for a plane beach.     

Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami

— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented.     

Reconstruction of Tsunami Inland Propagation on December 26, 2004 in Banda Aceh, Indonesia, through Field Investigations

This paper presents the results from an extensive field data collection effort following the December 26, 2004 earthquake and tsunami in Banda Aceh, Sumatra. The data were collected under the auspices of TSUNARISQUE, a joint French-Indonesian program dedicated to tsunami research and hazard mitigation, which has been active since before the 2004 event. In total, data from three months of field investigations are presented, which detail important aspects of the tsunami inundation dynamics in Banda Aceh. These include measurements of runup, tsunami wave heights, flow depths, flow directions, event chronology and building damage patterns. The result is a series of detailed inundation maps of the northern and western coasts of Sumatra including Banda Aceh and Lhok Nga. Among the more important findings, we obtained consistent accounts that approximately ten separate waves affected the region after the earthquake; this indicates a high-frequency component of the tsunami wave energy in the extreme near-field. The largest tsunami wave heights were on the order of 35 m with a maximum runup height of 51 m. This value is the highest runup value measured in human history for a seismically generated tsunami. In addition, our field investigations show a significant discontinuity in the tsunami wave heights and flow depths along a line approximately 3 km inland, which the authors interpret to be the location of the collapse of the main tsunami bore caused by sudden energy dissipation. The propagating bore looked like a breaking wave from the landward side although it has distinct characteristics. Patterns of building damage are related to the location of the propagating bore with overall less damage to buildings beyond the line where the bore collapsed. This data set was built to be of use to the tsunami community for the purposes of calibrating and improving existing tsunami inundation models, especially in the analysis of extreme near-field events.     

Numerical Modelling of the 26th December 2004 Indian Ocean Tsunami for the Southeastern Coast of India

A numerical simulation of the 26th December, 2004 Indian Ocean tsunami of the Tamil Nadu coastal zone is presented. The simulation approach is based on a fully nonlinear Boussinesq tsunami propagation model and included an accurate computational domain and a robust coseismic source. The simulation is first confronted to available tide gauge and runup observations. The agreement between observations and the predicted wave heights allowed a reasonable validation of the simulation. As a result, a full picture of the tsunami impact is provided over the entire coastal zone Tamil Nadu. The processes responsible for coastal vulnerability are discussed.     

Modeling of the 2011 Japan Tsunami: Lessons for Near-Field Forecast

During the devastating 11 March 2011 Japanese tsunami, data from two tsunami detectors were used to determine the tsunami source within 1.5 h of earthquake origin time. For the first time, multiple near-field tsunami measurements of the 2011 Japanese tsunami were used to demonstrate the accuracy of the National Oceanic and Atmospheric Administration (NOAA) real-time flooding forecast system in the far field. To test the accuracy of the same forecast system in the near field, a total of 11 numerical models with grids telescoped to 2 arcsec (~60 m) were developed to hindcast the propagation and coastal inundation of the 2011 Japanese tsunami along the entire east coastline of Japan. Using the NOAA tsunami source computed in near real-time, the model results of tsunami propagation are validated with tsunami time series measured at different water depths offshore and near shore along Japan’s coastline. The computed tsunami runup height and spatial distribution are highly consistent with post-tsunami survey data collected along the Japanese coastline. The computed inundation penetration also agrees well with survey data, giving a modeling accuracy of 85.5 % for the inundation areas along 800 km of coastline between Ibaraki Prefecture (north of Kashima) and Aomori Prefecture (south of Rokkasho). The inundation model results highlighted the variability of tsunami impact in response to different offshore bathymetry and flooded terrain. Comparison of tsunami sources inferred from different indirect methods shows the crucial importance of deep-ocean tsunami measurements for real-time tsunami forecasts. The agreement between model results and observations along Japan’s coastline demonstrate the ability and potential of NOAA’s methodology for real-time near-field tsunami flooding forecasts. An accurate tsunami flooding forecast within 30 min may now be possible using the NOAA forecast methodology with carefully placed tsunameters and large-scale high-resolution inundation models with powerful computing capabilities.     

Onshore tsunami deposits caused by the 1993 Southwest Hokkaido and 1983 Japan Sea earthquakes

Onshore tsunami deposits resulting from the 1993 Southwest Hokkaido and 1983 Japan Sea earthquakes were described to evaluate the feasibility of tsunami deposits for inferring paleoseismic events along submarine faults. Tsunami deposits were divided into three types, based on their composition and aerial distribution: (A) deposits consisting only of floating materials, (B) locally distributed siliclastic deposits, and (C) widespread siliclastic deposits. The most widely distributed tsunami deposits consist of the first two types. Type C deposits are mostly limited to areas where the higher tsunami runup was observed. The scale of tsunami represented by vertical tsunami runup is an important factor controlling the volume of tsunami deposits. The thickest deposits, about 10 cm, occur behind coastal dunes. To produce thick siliclastic tsunami deposits, a suitable source area, such as sand bar or dune, must be available in addition to sufficient vertical tsunami runup. Estimation of the amounts of erosion and deposition indicates that tsunami deposits were derived from both onshore and shoreface regions. The composition and grain size of the tsunami deposits strongly reflect the nature of the sedimentary materials of their source area. Sedimentary structures of the tsunami deposits suggest both low and high flow régimes. Consequently, it seems very difficult to identify tsunami deposits based only on grain size distribution or sedimentary structure of a single site in ancient successions.     

Lituya Bay Landslide Impact Generated Mega-Tsunami 50th Anniversary

On July 10, 1958, an earthquake M_w 8.3 along the Fairweather fault triggered a major subaerial landslide into Gilbert Inlet at the head of Lituya Bay on the southern coast of Alaska. The landslide impacted the water at high speed generating a giant tsunami and the highest wave runup in recorded history. The mega-tsunami runup to an elevation of 524 m caused total forest destruction and erosion down to bedrock on a spur ridge in direct prolongation of the slide axis. A cross section of Gilbert Inlet was rebuilt at 1:675 scale in a two-dimensional physical laboratory model based on the generalized Froude similarity. A pneumatic landslide tsunami generator was used to generate a high-speed granular slide with controlled impact characteristics. State-of-the-art laser measurement techniques such as particle image velocimetry (PIV) and laser distance sensors (LDS) were applied to the decisive initial phase with landslide impact and wave generation as well as the runup on the headland. PIV provided instantaneous velocity vector fields in a large area of interest and gave insight into kinematics of wave generation and runup. The entire process of a high-speed granular landslide impact may be subdivided into two main stages: (a) Landslide impact and penetration with flow separation, cavity formation and wave generation, and (b) air cavity collapse with landslide run-out and debris detrainment causing massive phase mixing. Formation of a large air cavity — similar to an asteroid impact — in the back of the landslide is highlighted. A three-dimenional pneumatic landslide tsunami generator was designed, constructed and successfully deployed in the tsunami wave basin at OSU. The Lituya Bay landslide was reproduced in a three-dimensional physical model at 1:400 scale. The landslide surface velocities distribution was measured with PIV. The measured tsunami amplitude and runup heights serve as benchmark for analytical and numerical models.     

Effects of bathymetry on tsunami propagation: Application of ray tracing to tsunamis

Ray tracing of seismic surface waves is applied to tsunami propagation to examine bathymetric effect along its propagation path. Computations are made for trans-Pacific tsunamis and for near-field tsunamis in the Japan Sea. For tsunamis across the Pacific Ocean, the comparison to a uniform ocean shows that focusing and defocusing, due to bathymetry, are significant for some combinations of source and receiver. For example, the refraction of rays is predominant at the East Pacific Rise for the tsunami from Chile. The tsunamis in the Japan Sea are strongly affected by the shallow Yamato Rise. The predicted arrival time and amplitude distribution generally agree with the observations from an actual tsunami. Since the computation can be made very quickly, the method is useful for preliminary analysis of tsunami propagation, such as in an operational warning system or in the determination of computational area for finite-difference computation.