The July 15, 2009 Fiordland,New Zealand Tsunami: Real-Time Assessment

On 15 July 2009, a Mw 7.8 earthquake occurred off the New Zealand coast, which by serendipitous coincidence occurred while the International Tsunami Symposium was in session in Novosibirsk, Russia. The earthquake generated a tsunami that propagated across the Tasman Sea and was detected in New Zealand, Australia and as far away as the US West coast. Small boats close to the epicenter were placed in jeopardy, but no significant damage was observed despite a measured run-up height of 2.3 m in one of the Sounds in close proximity to the source (Wilson in GNS Science Report 46:62 2009). Peak-to-trough tsunami heights of 55 cm were measured at Southport, Tasmania and a height of 1 m was measured in Jackson Bay, New Zealand. The International Tsunami Symposium provided an ideal venue for illustration of the value of immediate real-time assessment and provided an opportunity to further validate the real time forecasting capabilities with the scientific community in attendance. A number of agencies with responsibility for tsunami forecast and/or warning, such as the NOAA Center for Tsunami Research, the Pacific Tsunami Warning Center, GNS Science in New Zealand, the Australian Bureau of Meteorology and the European Commission Joint Research Centre were all represented at the meeting and were able to demonstrate the use of state of the art numerical models to assess the tsunami potential and provide warning as appropriate.     

Twin Tsunamis Triggered by the 12 January 2010 Haiti Earthquake

On 12 January 2010, a magnitude M _w 7.0 earthquake occurred 25 km west–southwest of Haiti’s capital Port-au-Prince causing an estimated 316,000 fatalities, thereby exceeding any previous loss of life from a similar size earthquake. In addition, tsunami waves triggered by the earthquake caused at least three fatalities at Petit Paradis due to a complete lack of tsunami awareness. The International Tsunami Survey Team (ITST) was deployed within weeks of the event and covered the greater Bay of Port-au-Prince and more than 100 km of Hispaniola’s southern coastline. The collected survey data include more than 21 tsunami heights along with observations of coastal land level change. Maximum tsunami heights of 3 m have been measured for two independently triggered tsunamis.     

Effects of the Tohoku tsunami of March 11, 2011 on the Pacific coast of Russia

The records of the Tohoku tsunami of March 11, 2011, obtained at the nearest Deep-Ocean Assessment and Reporting of Tsunamis (DART) stations and the coastal telemetric recorders of the Russian Tsunami Warning System, are analyzed. Such parameters as tsunami arrival times, heights of the first and maximal waves, and predominant periods are presented. The eyewitness accounts and photographs of tsunami effects are presented. The tsunami height distribution along the coast of Kuril Islands is discussed.     

Slip Distribution and Seismic Moment of the 2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms and Coastal Geodetic Data

The slip distribution and seismic moment of the 2010 and 1960 Chilean earthquakes were estimated from tsunami and coastal geodetic data. These two earthquakes generated transoceanic tsunamis, and the waveforms were recorded around the Pacific Ocean. In addition, coseismic coastal uplift and subsidence were measured around the source areas. For the 27 February 2010 Maule earthquake, inversion of the tsunami waveforms recorded at nearby coastal tide gauge and Deep Ocean Assessment and Reporting of Tsunamis (DART) stations combined with coastal geodetic data suggest two asperities: a northern one beneath the coast of Constitucion and a southern one around the Arauco Peninsula. The total fault length is approximately 400 km with seismic moment of 1.7 × 10²² Nm (Mw 8.8). The offshore DART tsunami waveforms require fault slips beneath the coasts, but the exact locations are better estimated by coastal geodetic data. The 22 May 1960 earthquake produced very large, ~30 m, slip off Valdivia. Joint inversion of tsunami waveforms, at tide gauge stations in South America, with coastal geodetic and leveling data shows total fault length of ~800 km and seismic moment of 7.2 × 10²² Nm (Mw 9.2). The seismic moment estimated from tsunami or joint inversion is similar to previous estimates from geodetic data, but much smaller than the results from seismic data analysis.     

W Phase Inversion and Tsunami Inundation Modeling for Tsunami Early Warning: Case Study for the 2011 Tohoku Event

Centroid moment tensor solutions for the 2011 Tohoku earthquake are determined by W phase inversions using 5 and 10 min data recorded by the Full Range Seismograph Network of Japan (F-net). By a scaling relation of moment magnitude to rupture area and an assumption of rigidity of 4 × 10¹⁰ N m^?2, simple rectangular earthquake fault models are estimated from the solutions. Tsunami inundations in the Sendai Plain, Minamisanriku, Rikuzentakata, and Taro are simulated using the estimated fault models. Then the simulated tsunami inundation area and heights are compared with the observations. Even the simulated tsunami heights and inundations from the W phase solution that used only 5 min data are considerably similar to the observations. The results are improved when using 10 min of W phase data. These show that the W phase solutions are reliable to be used for tsunami inundation modeling. Furthermore, the technique that combines W phase inversion and tsunami inundation modeling can produce results that have sufficient accuracy for tsunami early warning purposes.     

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.     

国际海啸预警系统(ITWS)

介绍了国际海啸预警系统的构成、地震与海啸信息的检测、海啸预警信息的发布，并介绍了太平洋海啸预警中心和阿拉斯加海啸预警中心。     

Field Survey and Geological Effects of the 15 November 2006 Kuril Tsunami in the Middle Kuril Islands

The near-field expression of the tsunami produced by the 15 November 2006 Kuril earthquake (M_w 8.1–8.4) in the middle Kuril Islands, Russia, including runup of up to 20 m, remained unknown until we conducted a post-tsunami survey in the summer of 2007. Because the earthquake occurred between summer field expeditions in 2006 and 2007, we have observations, topographic profiles, and photographs from three months before and nine months after the tsunami. We thoroughly surveyed portions of the islands of Simushir and Matua, and also did surveys on parts of Ketoi, Yankicha, Ryponkicha, and Rasshua. Tsunami runup in the near-field of the middle Kuril Islands, over a distance of about 200 km, averaged 10 m over 130 locations surveyed and was typically between 5 and 15 m. Local topography strongly affected inundation and somewhat affected runup. Higher runup generally occurred along steep, protruding headlands, whereas longer inundation distances occurred on lower, flatter coastal plains. Sediment transport was ubiquitous where sediment was available—deposit grain size was typically sand, but ranged from mud to large boulders. Wherever there were sandy beaches, a more or less continuous sand sheet was present on the coastal plain. Erosion was extensive, often more extensive than deposition in both space and volume, especially in areas with runup of more than 10 m. The tsunami eroded the beach landward, stripped vegetation, created scours and trim lines, cut through ridges, and plucked rocks out of the coastal plain.     

Observations and Impacts from the 2010 Chilean and 2011 Japanese Tsunamis in California (USA)

The coast of California was significantly impacted by two recent teletsunami events, one originating off the coast of Chile on February 27, 2010 and the other off Japan on March 11, 2011. These tsunamis caused extensive inundation and damage along the coast of their respective source regions. For the 2010 tsunami, the NOAA West Coast/Alaska Tsunami Warning Center issued a state-wide Tsunami Advisory based on forecasted tsunami amplitudes ranging from 0.18 to 1.43 m with the highest amplitudes predicted for central and southern California. For the 2011 tsunami, a Tsunami Warning was issued north of Point Conception and a Tsunami Advisory south of that location, with forecasted amplitudes ranging from 0.3 to 2.5 m, the highest expected for Crescent City. Because both teletsunamis arrived during low tide, the potential for significant inundation of dry land was greatly reduced during both events. However, both events created rapid water-level fluctuations and strong currents within harbors and along beaches, causing extensive damage in a number of harbors and challenging emergency managers in coastal jurisdictions. Field personnel were deployed prior to each tsunami to observe and measure physical effects at the coast. Post-event survey teams and questionnaires were used to gather information from both a physical effects and emergency response perspective. During the 2010 tsunami, a maximum tsunami amplitude of 1.2 m was observed at Pismo Beach, and over $3-million worth of damage to boats and docks occurred in nearly a dozen harbors, most significantly in Santa Cruz, Ventura, Mission Bay, and northern Shelter Island in San Diego Bay. During the 2011 tsunami, the maximum amplitude was measured at 2.47 m in Crescent City Harbor with over $50-million in damage to two dozen harbors. Those most significantly affected were Crescent City, Noyo River, Santa Cruz, Moss Landing, and southern Shelter Island. During both events, people on docks and near the ocean became at risk to injury with one fatality occurring during the 2011 tsunami at the mouth of the Klamath River. Evaluations of maximum forecasted tsunami amplitudes indicate that the average percent error was 38 and 28 % for the 2010 and 2011 events, respectively. Due to these recent events, the California tsunami program is developing products that will help: (1) the maritime community better understand tsunami hazards within their harbors, as well as if and where boats should go offshore to be safe, and (2) emergency managers develop evacuation plans for relatively small “Warning” level events where extensive evacuation is not required. Because tsunami-induced currents were responsible for most of the damage in these two events, modeled current velocity estimates should be incorporated into future forecast products from the warning centers.     

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.     

The 11 March 2011 Tohoku Tsunami Survey in Rikuzentakata and Comparison with Historical Events

On 11 March 2011, a moment magnitude M _w = 9.0 earthquake occurred off the Japan Tohoku coast causing catastrophic damage and loss of human lives. In the immediate aftermath of the earthquake, we conducted the reconnaissance survey in the city of Rikuzentakata, Japan. In comparison with three previous historical tsunamis impacting the same region, the 2011 event presented the largest values with respect to the tsunami height, the inundation area and the inundation distance. A representative tsunami height of 15 m was recorded in Rikuzentakata, with increased heights of 20 m around rocky headlands. In terms of the inundation area, the 2011 Tohoku tsunami exceeded by almost 2.6 times the area flooded by the 1960 Chilean tsunami, which ranks second among the four events compared. The maximum tsunami inundation distance was 8.1 km along the Kesen River, exceeding the 1933 Showa and 1960 Chilean tsunami inundations by factors of 6.2 and 2.7, respectively. The overland tsunami inundation distance was less than 2 km. The tsunami inundation height linearly decreased along the Kesen River at a rate of approximately 1 m/km. Nevertheless, the measured inland tsunami heights exhibit significant variations on local and regional scales. A designated “tsunami control forest” planted with a cross-shore width of about 200 m along a 2 km stretch of Rikuzentakata coastline was completely overrun and failed to protect the local community during this extreme event. Similarly, many designated tsunami shelters were too low and were overwashed by tsunami waves, thereby failing to provide shelter for evacuees—a risk that had been underestimated.     

Coastal Sedimentation Associated with the Tohoku Tsunami of 11 March 2011 in South Kuril Islands, NW Pacific Ocean

Sediment deposited by the Tohoku tsunami of March 11, 2011 in the Southern Kurils (Kunashir, Shikotan, Zeleniy, Yuri, Tanfiliev islands) was radically different from sedimentation during local strong storms and from tsunamis with larger runup at the same location. Sediments from the 2011 Tohoku tsunami were surveyed in the field, immediately and 6 months after the event, and analyzed in the laboratory for sediment granulometry, benthos Foraminifa assemblages, and diatom algae. Run-up elevation and inundation distance were calculated from the wrackline (accumulations of driftwood, woody debris, grass, and seaweed) marking the distal edge of tsunami inundation. Run-up of the tsunami was 5 m at maximum, and 3–4 m on average. Maximum distance of inundation was recorded in river mouths (up to 630 m), but was generally in the range of 50–80 m. Although similar to the local strong storms in runup height, the tsunami generally did not erode the coast, nor leave a deposit. However, deposits uncharacteristic of tsunami, described as brown aleuropelitic (silty and clayey) mud rich in organic matter, were found in closed bays facing the South Kuril Strait. These closed bays were covered with sea ice at the time of tsunami. As the tsunami waves broke the ice, the ice floes enhanced the bottom erosion on shoals and destruction of low-lying coastal peatland even at modest ranges of runup. In the muddy tsunami deposits, silt comprised up to 64 % and clay up to 41.5 %. The Foraminifera assemblages displayed features characteristic of benthic microfauna in the near-shore zone. Deep-sea diatoms recovered from tsunami deposits in two closely situated bays, namely Krabovaya and Otradnaya bays, had different requirements for environmental temperature, suggesting these different diatoms were brought to the bays by the tsunami wave entraining various water masses when skirting the island from the north and from the south.     

Relationship of Tsunami Intensity to Source Earthquake Magnitude as Retrieved from Historical Data

Operational prediction of near-field tsunamis in all existing Tsunami Warning Systems (TWSs) is based on fast determination of the position and size of submarine earthquakes. Exceedance of earthquake magnitude above some established threshold value, which can vary over different tsunamigenic zones, results in issuing a warning signal. Usually, a warning message has several (from 2 to 5) grades reflecting the degree of tsunami danger and sometimes contains expected wave heights at the coast. Current operational methodology is based on two main assumptions: (1) submarine earthquakes above some threshold magnitude can generate dangerous tsunamis and (2) the height of a resultant tsunami is, in general, proportional to the earthquake magnitude. While both assumptions are physically reasonable and generally correct, statistics of issued warnings are far from being satisfactory. For the last 55 years, up to 75% of warnings for regional tsunamis have turned out to be false, while each TWS has had at least a few cases of missing dangerous tsunamis. This paper presents the results of investigating the actual dependence of tsunami intensity on earthquake magnitude as it can be retrieved from historical observations and discusses the degree of correspondence of the above assumptions to real observations. Tsunami intensity, based on the Soloviev-Imamura scale is used as a measure of tsunami “size”. Its correlation with the M _s and M _w magnitudes is investigated based on historical data available for the instrumental period of observations (from 1900 to present).     

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock     

南中国海区域海啸预警中心地震监测系统

简要介绍了南中国海区域海啸预警与减灾系统的建设和发展历程,同时重点阐述了地震监测系统构成及其基本功能。作为重要组成部分,地震监测系统通过地震数据的实时汇集、存储、自动处理和分析,并结合人机交互方式实现了地震定位、震源机制解和有限断层模型反演。实际应用表明,地震监测系统对全球6.0级以上地震定位时间不超过8 min,在震后10—15 min内完成W震相方法快速反演海底强震震源机制解,在震后短时间内完成有限断层模型反演,为海啸预警提供快速、准确、可靠的地震基本参数和震源特征参数。      

Study on the relation of river morphology and tsunami propagation in rivers

The relations of river morphology and tsunami propagation in rivers were studied at several rivers in the Tohoku region during The Great Chilean Tsunami of 2010 and The Great East Japan Tsunami of 2011. It was found that river mouth morphological features play an important role in the intrusion of low magnitude tsunamis in which the geological and geographical conditions are an important factor. Nevertheless, the effects of these features were not found in the case of an extreme tsunami wave. As the wave enters the river, the propagation depends on other factors. It was found that the intrusion distance correlates well to the riverbed slope. The measurements of water level and riverbed slope were analyzed to propose an empirical method for estimating the damping coefficient for the tsunami propagation in rivers based on the tsunami of 2011. The proposed empirical method was used to approximate the length of the tsunami intrusion into a river by assuming that the furthest distance is given for the ratio of local tsunami wave height to the tsunami wave height at the river entrance of 0.05 (5 %). The estimated intrusion length from the proposed method in this study shows a good comparison with measurement data.     

葡萄牙破坏性地震和海啸预警系统(DETWS)

本文介绍了葡萄牙破坏性地震和海啸预警系统(Destructive Earthquakes and Tsunami Warning System)的构成、地震与海啸信息的检测、海啸预警信息的发布。     

An Evaluation of Tsunami Forecasts from the T2 Scenario Database

A tsunami scenario database (T2) has recently been developed for use within the Joint Australian Tsunami Warning Centre (JATWC). This scenario database has proven to be a very useful tool for forecast guidance, issuing of tsunami warnings and general event analysis. In this paper, the T2 scenarios are described, and evaluated by comparing them with observations of sea level from tsunameters for a number of recent tsunami events. In general, the T2 scenario database performs very well in terms of predicting the arrival time of the tsunami and the wave amplitudes at tsunameter locations.     

Earthquake parameters and spatial distribution of coseismic effects in Southern Siberia and Mongolia

In this work we review earthquakes that happened in Southern Siberia and Mongolia within the coordinates of 42°–62° N and 80°–124° E and first propose relationships between earthquake parameters (a surface-wave earthquake magnitude M _s and an epicentral intensity(I ₀) based on the MSK-64 scale) and maximal distances from an earthquake epicenter (R _{e max}), hypocenter (R _{h max}), and a seismogenic fault (R _{f max}) to the localities of secondary coseismic effects. Special attention was paid to the study of these relationships for the effects of soil liquefaction. Hence, it was shown that secondary deformations from an earthquake were distributed in space away from an earthquake epicenter, than from an associating seismogenic fault. The effects of soil liquefaction are manifested by several times closer to a seismogenic fault, than all other effects, regardless of the type of tectonic movement in a seismic focus. Within the 40 km zone from an earthquake epicenter 44% of the known manifestations of liquefaction process occurred; within the 40 km zone from a seismogenic fault—90%. We propose the next relationship for effects of soil liquefaction: M _s = 0.007 × R _{e max} + 5.168 that increases the limits of the maximum epicentral distance at an earthquake magnitude of 5.2 ≤ M _s ≤ 8.1 as compared to the corresponding relationships for different regions of the world.     

Probability of Peak Ground Horizontal and Peak Ground Vertical Accelerations at Tehran and Surrounding Areas

The active faults near Tehran are capable of producing Mw magnitudes of 6.62?C7.23; at epicentral distances of 25?km from the active faults, and Mw?=?7.23, the peak ground horizontal acceleration, PGH, is between 386 and 730?cm/s² and peak ground vertical acceleration, PGV, is between 192 and 261?cm/s²; the historic earthquake of the fourth century BC, Mw?=?7.16 produced the highest estimated PGH acceleration in Tehran, between 254 and 479?cm/s² and PGV acceleration between 127 and 173?cm/s². Earthquakes from 1909 to 2008, within 300?km from Tehran, are used for calculation of magnitude frequency relation, and results are applied to estimate PGH and PGV accelerations as a function of magnitudes for a set of fixed epicentral distance and site conditions; also as a function of epicentral distances for a set of fixed magnitudes and site conditions. Poisson??s distribution is used for probability calculation of PGH and PGV accelerations for several exposure times, site conditions and epicentral distances; accelerations with very high probability, near 1, are presented. At an epicentral distance of 10?km and exposure time of 450?years, in the northern part of Tehran, close to Mosha and the Northern Tehran faults, PGH acceleration is 800?C420?cm/s² and PGV is 400?C220?cm/s² with high probability. At an epicentral distance of 25?km and 1,000?years exposure time, PGH is 610?C320?cm/s² and PGV is 310?C160?cm/s² with high probability, where larger values are for soft soil and smaller values are for hard rock.