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.     

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

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

Finite element modeling of the July 12, 1993 Hokkaido Nansei-Oki tsunami

A fault plane model and a finite element hydrodynamic model are applied to the simulation of the Hokkaido Nansei-Oki tsunami of July 12, 1993. The joint performance of the models is assessed based on the overall ability to reproduce observed tsunami waveforms and to preserve mass and energy during tsunami propagation. While a number of observed characteristics of the waveforms are satisfactorily reproduced (in particular, amplitudes and arrival times at tidal gauges relatively close to the source, and general patterns of energy concentration), others are only marginally so (notably, wave periods at the same gauges, and wave heights along Okushiri); differences between observations and simulations are traceable to both the fault plane and the hydrodynamic models. Nonnegligible losses of energy occur throughout the simulated tsunami propagation. These losses seem to be due to a combination of factors, including numerical damping and possible deficiencies of the shallow water equations in preserving energy.     

Stratigraphic records of tsunamis along the Japan Sea,southwest Hokkaido,northern Japan

The stratigraphy of tsunami deposits along the Japan Sea, southwest Hokkaido, northern Japan, reveals tsunami recurrences in this particular area. Sandy tsunami deposits are preserved in small valley plains, whereas gravelly deposits of possible tsunami origin are identified in surficial soils covering a Holocene marine terrace and a slope talus. At least five horizons of tsunami events can be defined in the Okushiri Island, the youngest of which immediately overlies the Ko‐d tephra layer (1640 AD) and was likely formed by the historical Oshima‐Ohshima tsunami in 1741 AD. The four older tsunami deposits, dated using accelerator mass spectrometry ¹⁴C, were formed at around the 12th century, 1.5–1.6, 2.4–2.6, and 2.8–3.1 ka, respectively. Tsunami sand beds of the 1741 AD and circa 12th century events are recognized in the Hiyama District of Hokkaido Island, but the older tsunami deposits are missing. The deposits of these two tsunamis are found together at the same sites and distributed in regions where wave heights of the 1993 tsunami (Hokkaido Nansei‐oki earthquake, Mw = 7.7) were less than 3 m. Thus, the 12th century tsunami waves were possibly generated near the south of Okushiri Island, whereas the 1993 tsunami was generated towards the north of the island. The estimated recurrence intervals of paleotsunamis, 200–1100 years with an average of 500 years, likely represents the recurrence interval of large earthquakes which would have occurred along several active faults offshore of southwest Hokkaido.     

Tsunami deposits from the 1993 Southwest Hokkaido earthquake and the 1640 Hokkaido Komagatake eruption,northern Japan

The southwest Hokkaido tsunami of July 12th, 1993, left continuous onshore sand deposits along the west coast of Oshima Peninsuka, Hokkaido, northern Japan. We investigated spatial distribution and lithofacies of the new tsunami deposits for its identification of ancient tsunami deposits. An eyewitness acount and bent plants helped our interpretation of the onshore tsunami behavior. We regard the following properties as typical of the coastal tsunami sand deposits: (1) The deposits cover the surface almost continuously on gentle topography. (2) Deposit thicknesses and mean grain sizes descrease with distance from the sea. (3) Deposit thicknesses and lithofacies vary greatly across local surface undulation. (4) Graded bedding reflecting tsunami runup and backwash is present in thick deposits. (5) The deposits are widely distributed along the coast and extend inland several tens of meters to 100 m. We examined a candidate for the paleo-tsunami deposits associated with the 1640 Komagatake eruption, and confirmed that the similar patterns are typical of ancient tsunami deposits.     

An Overview of the 28 October 2012 M<Subscript>w</Subscript> 7.7 Earthquake in Haida Gwaii,Canada: A Tsunamigenic Thrust Event Along a Predominantly Strike-Slip Margin

The boundary between the Pacific and North America plates along Canada’s west coast is one of the most seismically active regions of Canada, and is where Canada’s two largest instrumentally recorded earthquakes have occurred. Although this is a predominantly strike-slip transform fault boundary, there is a component of oblique convergence between the Pacific and North America plates off Haida Gwaii. The 2012 M_w 7.7 Haida Gwaii earthquake was a thrust event that generated a tsunami with significant run up of over 7 m in several inlets on the west coast of Moresby Island (several over 6 m, with a maximum of 13 m). Damage from this earthquake and tsunami was minor due to the lack of population and vulnerable structures on this coast.     

Waveform and Spectral Analyses of the 2011 Japan Tsunami Records on Tide Gauge and DART Stations Across the Pacific Ocean

We studied the 11 March 2011 Tohoku tsunami through analysis of the sea level records from 21 tide gauge and 16 DART (Deep-ocean Assessment and Reporting of Tsunamis) stations from across the Pacific Ocean. The extreme power of this trans-oceanic tsunami was indicated by the trough-to-crest heights of 3.03 m at Arena Cove on the western coast of the USA and 3.94 m at Coquimbo on the southern coast of Chile. The average value of the maximum amplitude was 163.9 cm for the examined tide gauge records. At many coastal tide gauge stations the largest wave arrived several hours after the first arrival of the tsunami wave, and the tsunami lasted for a long time with an average duration of 4 days. On the contrary, at most of the DART stations in the deep ocean, the first wave was the largest, the tsunami amplitudes were smaller with an average maximum of 51.2 cm, and the durations were shorter with an average of 2 days. The two dominant tsunami periods on the DART records were 37 and 67.4 min, which are possibly attributed to the width and length of the tsunami source fault, respectively. The dimensions of the tsunami source was estimated as 233 km × 424 km. Wavelet analyses of tide gauge and DART records showed that most of the tsunami energy was distributed at the wide period band of around 10–80 min during the first hour after the tsunami arrival, then it was concentrated in a relatively narrower band. The frequency-time plots showed the switches and lapses of tsunami energy at the 35- and 65-min period bands.     

The Effect of the Great Barrier Reef on the Propagation of the 2007 Solomon Islands Tsunami Recorded in Northeastern Australia

The effect of offshore coral reefs on the impact from a tsunami remains controversial. For example, field surveys after the 2004 Indian Ocean tsunami indicate that the energy of the tsunami was reduced by natural coral reef barriers in Sri Lanka, but there was no indication that coral reefs off Banda Aceh, Indonesia had any effect on the tsunami. In this paper, we investigate whether the Great Barrier Reef (GBR) offshore Queensland, Australia, may have weakened the tsunami impact from the 2007 Solomon Islands earthquake. The fault slip distribution of the 2007 Solomon Islands earthquake was firstly obtained by teleseismic inversion. The tsunami was then propagated to shallow water just offshore the coast by solving the linear shallow water equations using a staggered grid finite-difference method. We used a relatively high resolution (approximately 250 m) bathymetric grid for the region just off the coast containing the reef. The tsunami waveforms recorded at tide gauge stations along the Australian coast were then compared to the results from the tsunami simulation when using both the realistic 250 m resolution bathymetry and with two grids having fictitious bathymetry: One in which the the GBR has been replaced by a smooth interpolation from depths outside the GBR to the coast (the “No GBR” grid), and one in which the GBR has been replaced by a flat plane at a depth equal to the mean water depth of the GBR (the “Average GBR” grid). From the comparison between the synthetic waveforms both with and without the Great Barrier Reef, we found that the Great Barrier Reef significantly weakened the tsunami impact. According to our model, the coral reefs delayed the tsunami arrival time by 5–10 minutes, decreased the amplitude of the first tsunami pulse to half or less, and lengthened the period of the tsunami.     

Tide Gauge Observations of 2004–2007 Indian Ocean Tsunamis from Sri Lanka and Western Australia

Tide gauge data collected from Sri Lanka (three stations) and Western Australia (eleven stations) during the Indian Ocean tsunamis, which occurred in December 2004, March 2005, July 2006, and September 2007, and incorporated five tsunamis, were examined to determine tsunami behaviour during these events. During the December 2004 tsunami, maximum wave heights of 3.87 m and 1.75 m were recorded at Colombo (Sri Lanka) and Bunbury (Western Australia), respectively. The results indicated that although the relative magnitudes of the tsunamis varied, the tsunami behaviour at each station was similar. This was due to the effect of the local and regional topography. At all tide gauges, the spectral energy corresponding to periods between 20 and 85 minutes increased during the tsunami. The sea-level data obtained from the west and south coasts of Sri Lanka (Colombo and Kirinda) indicated the importance of wave reflections from the Maldives Island chain, which produced the maximum wave two to three hours after the arrival of the first wave. In contrast, Trincomalee on the east coast did not show evidence of a reflected wave. Similarly, along the west coast of Australia, the highest waves occurred 15 hours after the arrival of the first wave. Here, based on travel times, we postulated that the waves were reflected from the Mascarene Ridge and/or the Island of Madagascar. Reflected waves were not present in the 2006 tsunami, where the primary waves propagated away from topographic features. One of the main influences of the tsunami was to set up oscillations at the local resonance frequency. Because Sri Lanka and Western Australia have relatively straight coastlines, these oscillations were related to the fundamental period of the shelf oscillation. For Colombo, this corresponded to 75-minute period, whereas in Geraldton and Busselton (Australia), the four-hour period was most prominent; at Jurien Bay and Fremantle, the resonance period was 2.7 hours.     

Evidence on the Koseda coast of Yakushima Island of a tsunami during the 7.3 ka Kikai caldera eruption

Previous research indicates that Yakushima Island, southwestern Japan, may have been struck by a huge tsunami before or soon after the arrival of the Koya pyroclastic flow during the 7.3 ka caldera‐forming Kikai eruption, but this has not yet been confirmed. This paper describes sedimentological and chronostratigraphic evidence showing that Unit TG, one of three gravel beds exposed on the Koseda coast of northeast Yakushima Island and investigated here, is a tsunami deposit. Unit TG is a poorly sorted, 30 cm thick gravel bed overlying a wave‐cut bench and underlying a Koya pyroclastic flow deposit. Sparse wood fragments in Unit TG were dated at 7 416–7 167 cal year BP. The constituent gravel clasts of Unit TG are similar in composition to those of modern beach and river deposits along the Koseda coast. Unit TG also contains pumice clasts whose chemistry is identical to that of pumice derived from the 7.3 ka eruption at Kikai caldera. The long‐axis orientations and composition of gravel clasts in Unit TG suggest that they were transported by a landward‐travelling high‐particle‐concentration flow, which suggests that Unit TG was deposited by a tsunami run‐up flow during the 7.3 ka Kikai caldera eruption, just before the arrival of the major Koya pyroclastic flow at the Koseda coast. Whether the 7.3 ka tsunami was caused by a volcanic eruption or an earthquake remains unclear, but Unit TG demonstrates that a tsunami arrived immediately before emplacement of a Koya pyroclastic flow.     

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.     

Sumatra tsunami: lessons from modeling

The need for the combination of seismic data with real-time wave height information for an effective prediction of tsunami impact is emphasized in the paper. A preliminary, but comprehensive study of arrival times, wave heights and run-up values at a number of locations and tide gage stations throughout the Indian Ocean seaboard is presented. Open ocean wave height data from satellite observations are analyzed and used in the reconstruction of a tsunami source mechanism for the December 26, 2004 event. The reconstructed source is then used to numerically estimate tsunami impact along the Indian Ocean seaboard, including wave height, and arrival times at 12 tide gage stations, and inundation at 3 locations on the coast of India. The December 2004, as well as the March 28, 2005 tsunamis are investigated and their differences in terms of tsunami generation are analyzed and presented as a clear example of the need for both, seismic and real-time tsunami data for a reliable tsunami warning system in the Indian Ocean.     

东海海域潜在地震海啸的数值模拟初步研究

针对东海冲绳海槽地区的地震地质背景,对东海海域潜在的地震海啸进行了预研究. 假设了冲绳海槽在发生8.5级大地震,断层错距高达15 m的极端地震情况引发的海啸对中国东部沿海地区的影响. 初步数值模拟结果表明,该地震引发的海啸的最大初始波高为4.3 m,4小时左右传至浙江沿岸,近岸各处波高为1——2 m,其中局部地区波高为2.4 m;约7——8小时靠近上海海岸线（若震源在中冲绳海槽地区,海啸传到上海最快大约7小时）,近岸波高约为1 m. 近岸区域地形变化复杂,海岛密布,局部地形条件可能会很大地影响实际各地点海啸波高,加上海啸在岸边爬高及港湾效应,估计波高还会升高. 给出了冲绳海槽南、中、北部发生潜在地震海啸的传播等时图. 笔者在东海设置了3个地震及海啸监测站,基于海啸模拟结果绘制了监测站处的海啸随时间演化曲线,分析了预研究成果对海啸预警可能发挥的作用.      

Numerical Simulations of Tsunami Waves and Currents for Southern Vancouver Island from a Cascadia Megathrust Earthquake

The 1700 great Cascadia earthquake (M = 9) generated widespread tsunami waves that affected the entire Pacific Ocean and caused damage as distant as Japan. Similar catastrophic waves may be generated by a future Cascadia megathrust earthquake. We use three rupture scenarios for this earthquake in numerical experiments to study propagation of tsunami waves off the west coast of North America and to predict tsunami heights and currents in several bays and harbours on southern Vancouver Island, British Columbia, including Ucluelet, located on the west coast of the island, and Victoria and Esquimalt harbours inside Juan de Fuca Strait. The earthquake scenarios are: an 1100-km long rupture over the entire length of the subduction zone and separate ruptures of its northern or southern segments. As expected, the southern earthquake scenario has a limited effect over most of the Vancouver Island coast, with waves in the harbours not exceeding 1 m. The other two scenarios produce large tsunami waves, higher than 16 m at one location near Ucluelet and over 4 m inside Esquimalt and Victoria harbours, and very strong currents that reach 17 m/s in narrow channels and near headlands. Because the assumed rupture scenarios are based on a previous earthquake, direct use of the model results to estimate the effect of a future earthquake requires appropriate qualification.     

Field survey of the 1994 Mindoro Island,Philippines tsunami

This is a report of the field survey of the November 15, 1994 Mindoro Island, Philippines, tsunami generated by an earthquake (M=7.0) with a strike-slip motion. We will report runup heights from 54 locations on Luzon, Mindoro and other smaller islands in the Cape Verde passage between Mindoro and Luzon. Most of the damage was concentrated along the northern coast of Mindoro. Runup height distribution ranged 3–4 m at the most severely damaged areas and 2–4 in neighboring areas. The tsunami-affected area was limited to within 10 km of the epicenter. The largest recorded runup value of 7.3 m was measured on the southwestern coast of Baco Island while a runup of 6.1 m was detected on its northern coastline. The earthquake and tsunami killed 62 people, injured 248 and destroyed 800 houses. As observed in other recent tsunami disasters, most of the casualties were children. Nearly all eyewitnesses interviewed described the first wave as a leading-depression wave. Eyewitnesses reported that the main direction of tsunami propagation was SW in Subaang Bay, SE in Wawa and Calapan, NE on Baco Island and N on Verde Island, suggesting that the tsunami source area was in the southern Pass of Verde Island and that the wave propagated rapidly in all directions. The fault plane extended offshore to the N of Mindoro Island, with its rupture originating S of Verde Island and propagating almost directly south to the inland of Mindoro, thereby accounting for the relatively limited damage area observed on the N of Mindoro.     

The 2010 Chilean Tsunami Off the West Coast of Canada and the Northwest Coast of the United States

The major (M _w = 8.8) Chilean earthquake of 27 February 2010 generated a trans-oceanic tsunami that was observed throughout the Pacific Ocean. Waves associated with this event had features similar to those of the 1960 tsunami generated in the same region by the Great (M _w = 9.5) 1960 Chilean Earthquake. Both tsunamis were clearly observed on the coast of British Columbia. The 1960 tsunami was measured by 17 analog pen-and-paper tide gauges, while the 2010 tsunami was measured by 11 modern digital coastal tide gauges, four NEPTUNE-Canada bottom pressure recorders located offshore from southern Vancouver Island, and two nearby open-ocean DART stations. The 2010 records were augmented by data from seven NOAA tide gauges on the coast of Washington State. This study examines the principal characteristics of the waves from the 2010 event (height, period, duration, and arrival and travel times) and compares these properties for the west coast of Canada with corresponding properties of the 1960 tsunami. Results show that the 2010 waves were approximately 3.5 times smaller than the 1960 waves and reached the British Columbia coast 1 h earlier. The maximum 2010 wave heights were observed at Port Alberni (98.4 cm) and Winter Harbour (68.3 cm); the observed periods ranged from 12 min at Port Hardy to 110–120 min at Prince Rupert and Port Alberni and 150 min at Bamfield. The open-ocean records had maximum wave heights of 6–11 cm and typical periods of 7 and 15 min. Coastal and open-ocean tsunami records revealed persistent oscillations that “rang” for 3–4 days. Tsunami energy occupied a broad band of periods from 3 to 300 min. Estimation of the inverse celerity vectors from cross-correlation analysis of the deep-sea tsunami records shows that the tsunami waves underwent refraction as they approached the coast of Vancouver Island with the direction of the incoming waves changing from an initial direction of 340° True to a direction of 15° True for the second train of waves that arrived 7 h later after possible reflection from the Marquesas and Hawaiian islands.     

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.     

Preliminary numerical simulation of potential earthquake-induced tsunami in East China Sea

In this paper,we present a numerical simulation of the propagation of a tsunami in the East China Sea,which might be induced by a hypothetical M8.5 earthquake in Okinawa Trough. Our results show that the initial maxi-mum wave height of tsunami could reach as high as 4.3 m for the hypothetical earthquake. It would take 3.5～4 hours for the tsunami to propagate to the coast of Zhejiang Province,and 7～8 hours to the near-shore of Shanghai. The peak tsunami height could be up to about 2 m in the coast of Zhejiang Province. Based on the numerical ex-periments,we plot the arrival time contours of tsunami in East China Sea and time history curves on the three ob-servational stations,and discussed the significance of the pre-analysis.     

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.     

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.