Tsunami Observations Research Articles

Validation of linearity assumptions for using tsunami waveforms in joint inversion of kinematic rupture models: Application to the 2010 Mentawai Mw 7.8 tsunami earthquake

AbstractTsunami observations have particular importance for resolving shallow offshore slip in finite‐fault rupture model inversions for large subduction zone earthquakes. However, validations of amplitude linearity and choice of subfault discretization of tsunami Green's functions are essential when inverting tsunami waveforms. We explore such validations using four tsunami recordings of the 25 October 2010 Mentawai Mw 7.8 tsunami earthquake, jointly inverted with teleseismic body waves and 1 Hz GPS (high‐rate GPS) observations. The tsunami observations include near‐field and far‐field deep water recordings, as well as coastal and island tide gauge recordings. A nonlinear, dispersive modeling code, NEOWAVE, is used to construct tsunami Green's functions from seafloor excitation for the linear inversions, along with performing full‐scale calculations of the tsunami for the inverted models. We explore linearity and finiteness effects with respect to slip magnitude, variable rake determination, and subfault dimensions. The linearity assumption is generally robust for the deep water recordings, and wave dispersion from seafloor excitation is important for accurate description of near‐field Green's functions. Breakdown of linearity produces substantial misfits for short‐wavelength signals in tide gauge recordings with large wave heights. Including the tsunami observations in joint inversions provides improved resolution of near‐trench slip compared with inversions of only seismic and geodetic data. Two rupture models, with fine‐grid (15 km) and coarse‐grid (30 km) spacing, are inverted for the Mentawai event. Stronger regularization is required for the fine model representation. Both models indicate a shallow concentration of large slip near the trench with peak slip of ~15 m. Fully nonlinear forward modeling of tsunami waveforms confirms the validity of these two models for matching the tsunami recordings along with the other data.

Parallel Implementation of Dispersive Tsunami Wave Modeling with a Nesting Algorithm for the 2011 Tohoku Tsunami

Because of improvements in offshore tsunami observation technology, dispersion phenomena during tsunami propagation have often been observed in recent tsunamis, for example the 2004 Indian Ocean and 2011 Tohoku tsunamis. The dispersive propagation of tsunamis can be simulated by use of the Boussinesq model, but the model demands many computational resources. However, rapid progress has been made in parallel computing technology. In this study, we investigated a parallelized approach for dispersive tsunami wave modeling. Our new parallel software solves the nonlinear Boussinesq dispersive equations in spherical coordinates. A variable nested algorithm was used to increase spatial resolution in the target region. The software can also be used to predict tsunami inundation on land. We used the dispersive tsunami model to simulate the 2011 Tohoku earthquake on the Supercomputer K. Good agreement was apparent between the dispersive wave model results and the tsunami waveforms observed offshore. The finest bathymetric grid interval was 2/9 arcsec (approx. 5 m) along longitude and latitude lines. Use of this grid simulated tsunami soliton fission near the Sendai coast. Incorporating the three-dimensional shape of buildings and structures led to improved modeling of tsunami inundation.

Tsunami Observations in Rivers from a Perspective of Tsunami Interaction with Tide and Riverine Flow

The observations of the 2011 Tohoku tsunami and the 2010 Chilean tsunami in several rivers in Japan and in the Columbia River in the USA are analyzed for patterns of tsunami behavior in river environments. Tsunamis in rivers exhibit actions very different from those observed on an open coast, but very similar among different rivers, though the action scale in different rivers varies greatly. We describe two tsunami effects in rivers as observed in field data. First, the river tide modulates the tsunami wave in a very specific way common to all rivers. Second, a strong near-field tsunami can cause significant prolonged water accumulation in lower river reaches. Both effects are inherent in tidal river environments, and have been reproduced numerically in a simplified 1-D river using a non-linear, shallow-water model with bottom friction. The numerical experiments highlight the indispensable role of a tsunami’s interaction with tide and riverine flow.

Observations and Numerical Modeling of the 2012 Haida Gwaii Tsunami off the Coast of British Columbia

A major (M w 7.7) earthquake occurred on October 28, 2012 along the Queen Charlotte Fault Zone off the west coast of Haida Gwaii (formerly the Queen Charlotte Islands). The earthquake was the second strongest instrumentally recorded earthquake in Canadian history and generated the largest local tsunami ever recorded on the coast of British Columbia. A field survey on the Pacific side of Haida Gwaii revealed maximum runup heights of up to 7.6 m at sites sheltered from storm waves and 13 m in a small inlet that is less sheltered from storms (Leonard and Bednarski 2014). The tsunami was recorded by tide gauges along the coast of British Columbia, by open-ocean bottom pressure sensors of the NEPTUNE facility at Ocean Networks Canada’s cabled observatory located seaward of southwestern Vancouver Island, and by several DART stations located in the northeast Pacific. The tsunami observations, in combination with rigorous numerical modeling, enabled us to determine the physical properties of this event and to correct the location of the tsunami source with respect to the initial geophysical estimates. The initial model results were used to specify sites of particular interest for post-tsunami field surveys on the coast of Moresby Island (Haida Gwaii), while field survey observations (Leonard and Bednarski 2014) were used, in turn, to verify the numerical simulations based on the corrected source region.

Tsunami Simulation Using Sources Inferred from Various Measurement Data: Implications for the Model Forecast

Model forecast applications use various models of tsunami sources inferred from different measurement data. Even the same type of observation data can produce substantially different tsunami source models during a real-time forecast when more data are obtained during the real-time analysis. Improved tsunami observations enable investigation of the influence of such model source variability on the final forecast using different source data sets of several events. The 2010 Maule, Chile and 2011 Tohoku, Japan tsunamis were two recent events that provide ample observations throughout the Pacific and were, thus, used here to study the sensitivity of different model inputs for forecasting. The sources for these events were derived using the following three different methods: (1) real time or post event inversion of tsunameter water level data; (2) prediction of sea floor deformations via analysis of seismic wave forms and application of a finite fault model; and (3) prediction of sea floor deformation using real-time GPS data. For the March 11, 2011 Tohoku tsunami, two examples of each method are used, while for the February 27, 2010 Maule event, only one tsunameter inversion and one finite fault model method were used due to a much more limited data set. Observed data from the Deep-ocean Assessment and Reporting for Tsunamis (DART) network, Japan GPS buoys, and select tide gauges across the Pacific were compared with forecasts to assess the sensitivity of these three methods using root-mean-square error analysis. We divided the analysis by the type of data and the distance from the source. This sensitivity analysis showed that increasing the resolution of a tsunami source model does not necessarily improve tsunami forecast quality, even in the near-field. Instead, the findings suggest that when forecasting coastal impact, defining the overall energy characteristic of a tsunami source may be more important than refining small source details. Source models based on direct tsunami observations are better at reproducing a tsunami signal: this finding is not very surprising but has implications for tsunami forecasting and warning operations.

Advances in earthquake and tsunami sciences and disaster risk reduction since the 2004 Indian ocean tsunami

The December 2004 Indian Ocean tsunami was the worst tsunami disaster in the world’s history with more than 200,000 casualties. This disaster was attributed to giant size (magnitude M ~ 9, source length >1000 km) of the earthquake, lacks of expectation of such an earthquake, tsunami warning system, knowledge and preparedness for tsunamis in the Indian Ocean countries. In the last ten years, seismology and tsunami sciences as well as tsunami disaster risk reduction have significantly developed. Progress in seismology includes implementation of earthquake early warning, real-time estimation of earthquake source parameters and tsunami potential, paleoseismological studies on past earthquakes and tsunamis, studies of probable maximum size, recurrence variability, and long-term forecast of large earthquakes in subduction zones. Progress in tsunami science includes accurate modeling of tsunami source such as contribution of horizontal components or “tsunami earthquakes”, development of new types of offshore and deep ocean tsunami observation systems such as GPS buoys or bottom pressure gauges, deployments of DART gauges in the Pacific and other oceans, improvements in tsunami propagation modeling, and real-time inversion or data assimilation for the tsunami warning. These developments have been utilized for tsunami disaster reduction in the forms of tsunami early warning systems, tsunami hazard maps, and probabilistic tsunami hazard assessments. Some of the above scientific developments helped to reveal the source characteristics of the 2011 Tohoku earthquake, which caused devastating tsunami damage in Japan and Fukushima Dai-ichi Nuclear Power Station accident. Toward tsunami disaster risk reduction, interdisciplinary and trans-disciplinary approaches are needed for scientists with other stakeholders.

Localized fault slip to the trench in the 2010 Maule, Chile Mw = 8.8 earthquake from joint inversion of high‐rate GPS, teleseismic body waves, InSAR, campaign GPS, and tsunami observations

AbstractThe 27 February 2010, Mw 8.8 Maule earthquake ruptured ~500 km along the plate boundary offshore central Chile between 34°S and 38.5°S. Establishing whether coseismic fault offset extended to the trench is important for interpreting both shallow frictional behavior and potential for tsunami earthquakes in the region. Joint inversion of high‐rate GPS, teleseismic body waves, interferometric synthetic aperture radar (InSAR), campaign GPS, and tsunami observations yields a kinematic rupture model with improved resolution of slip near the trench. Bilateral rupture expansion is resolved in our model with relatively uniform slip of 5–10 m downdip beneath the coast and two near‐trench high‐slip patches with >12 m displacements. The peak slip is ~17 m at a depth of ~15 km on the central megathrust, located ~200 km north from the hypocenter and overlapping the rupture zone of the 1928 M ~8 event. The updip slip is ~16 m near the trench. Another shallow near‐trench patch is located ~150 km southwest of the hypocenter, with a peak slip of 12 m. Checkerboard resolution tests demonstrate that correctly modeled tsunami data are critical to resolution of slip near the trench, with other data sets allowing, but not requiring slip far offshore. Large interplate aftershocks have a complementary distribution to the coseismic slip pattern, filling in gaps or outlining edges of large‐slip zones. Two clusters of normal faulting events locate seaward along the plate motion direction from the localized regions of large near‐trench slip, suggesting that proximity of slip to the trench enhanced extensional faulting in the underthrusting plate.

Tsunami observations in the open ocean

Deep-sea tsunami measurements play a major role in understanding the physics of tsunami wave generation and propagation, and in the creation of an effective tsunami warning system. The paper provides an overview of the history of tsunami recording in the open ocean from the beginning (about 50 years ago) to the present day. It describes modern tsunami monitoring systems, including the Deep-ocean Assessment and Reporting of Tsunamis (DART), innovative Japanese bottom cable projects, and the NEPTUNE-Canada geophysical bottom observatory. The specific peculiarities of seafloor longwave observations in the deep ocean are discussed and compared with those recorded in coastal regions. Tsunami detection in bottom presure observations is exemplified based on analysis of distant (22000 km) records of the 2004 Sumatra tsunami in the northeastern Pacific.

Simulation of Tsunami Hazards Affecting the East Cape Region, New Zealand

The east coast of New Zealand is known for being exposed to a variety of tsunami sources, both those arising from the nearby Hikurangi subduction zone and its associated crustal faults, and those arising from more distant parts of the Pacific. Using numerical simulations with a parallelized computer model, we assess the tsunami inundation hazard posed to the most populated coastal communities on the east coast of the East Cape region, New Zealand, which had not been evaluated before. Our tsunami inundation hazard study is based on severe but realistic scenario events from a selection of local and distant earthquake tsunami sources. Such modelling covers a gap in knowledge caused by the short historical record of tsunamis in this region (covering <200 years), and the sparse observations of historical tsunamis due to a largely rural population. We identify that the worst flooding is often not associated with the first wave arrivals, and that coastal oscillations can last a long time in distant events. The modelling results allow us to determine typical characteristics of the pattern of flooding, specific to an area. This information is important for emergency planning and preparedness.

Rupture process of the 2010 Mw 7.8 Mentawai tsunami earthquake from joint inversion of near‐field hr‐GPS and teleseismic body wave recordings constrained by tsunami observations

AbstractThe 25 October 2010 Mentawai tsunami earthquake (Mw 7.8) ruptured the shallow portion of the Sunda megathrust seaward of the Mentawai Islands, offshore of Sumatra, Indonesia, generating a strong tsunami that took 509 lives. The rupture zone was updip of those of the 12 September 2007 Mw 8.5 and 7.9 underthrusting earthquakes. High‐rate (1 s sampling) GPS instruments of the Sumatra GPS Array network deployed on the Mentawai Islands and Sumatra mainland recorded time‐varying and static ground displacements at epicentral distances from 49 to 322 km. Azimuthally distributed tsunami recordings from two deepwater sensors and two tide gauges that have local high‐resolution bathymetric information provide additional constraints on the source process. Finite‐fault rupture models, obtained by joint inversion of the high‐rate (hr)‐GPS time series and numerous teleseismic broadband P and S wave seismograms together with iterative forward modeling of the tsunami recordings, indicate rupture propagation ~50 km up dip and ~100 km northwest along strike from the hypocenter, with a rupture velocity of ~1.8 km/s. Subregions with large slip extend from 7 to 10 km depth ~80 km northwest from the hypocenter with a maximum slip of 8 m and from ~5 km depth to beneath thin horizontal sedimentary layers beyond the prism deformation front for ~100 km along strike, with a localized region having &gt;15 m of slip. The seismic moment is 7.2 × 1020 N m. The rupture model indicates that local heterogeneities in the shallow megathrust can accumulate strain that allows some regions near the toe of accretionary prisms to fail in tsunami earthquakes.

Observations of meteorological tsunamis along the south-west Australian coast

Tide gauge data were used to identify the occurrence, characteristics, and cause of tsunamis of meteorological origin (termed ‘meteotsunamis’) along the Western Australian coast. This is the first study to identify meteotsunamis in this region, and the results indicated that they occur frequently. Although meteotsunamis are not catastrophic to the extent of major seismically induced basin-scale events, the wave heights of meteotsunamis examined at some local stations in this study were higher than those recorded through seismic tsunamis. In June 2012, a meteotsunami contributed to an extreme water-level event at Fremantle, which recorded the highest water level in over 115 years. Meteotsunamis (wave heights >0.4 m, when the mean tidal range in the region is ~0.5 m) were found to coincide with thunderstorms in summer and the passage of low-pressure systems during winter. Spectral analysis of tide gauge time series records showed that existing continental seiche oscillations (periods between 30 min and 5 h) were enhanced during the meteotsunamis, with a high proportion of energy transferred to the continental shelf oscillation period. Three recent meteotsunami events (22 March 2010, 10 June 2012, and 7 January 2013) two due to summer thunderstorms and one due to a winter frontal system were chosen for detailed analysis. The meteotsunami amplitudes were up to a factor 2 larger than the local tidal range and sometimes contributed up to 85 % of the non-tidal water signal. A single meteorological event was found to generate several meteotsunamis along the coast, up to 500 km apart, as the air pressure disturbance propagated over the continental shelf; however, the topography and local bathymetry of the continental shelf defined the local sea-level resonance characteristics at each location. With the available data (sea level and meteorological), the exact mechanisms for the generation of the meteotsunamis could not be isolated.

Special Issue on Challenges of Earthquake Forecast Research Illuminated by the 2011 Tohoku-Oki Earthquake

The 2011 off the Pacific coast of Tohoku Earthquake, hereafter referred to as Tohoku-Oki earthquake, occurring off northeastern Japan’s Pacific coast on March 11, 2011 had a moment magnitude of 9.0 and generated a tsunami responsible for most of the deaths of the event’s 19,000 victims. Identifying scientifically what happened before, on, and after March 11 is one starting point for a discussion on how to reduce casualties and mitigate the impact of such natural disasters. The 14 papers in this special issue cover incidents related to pre-, co- and post-seismic phenomena, including volcanoes. Three papers discuss why and how such a large quake occurred. Three more papers go into the implications of short- and long-term crustal deformations seen in northeastern Japan. Four papers detail short- and long-term phenomena leading to the Tohoku-Oki quake. Two papers discuss real-time tsunami forecasting based on off-shore and on-shore geodetic, seismic and tsunami observation data. The last two papers explore the effects of the 2011 temblor on volcanic phenomena. The magnitude 9.0 produced in the 2011 event is the largest historically recorded in Japan and may not necessarily have been anticipated beforehand, and the generation mechanism behind such a gigantic occurrence is not yet completely understood. Even so, preparations should be made for such earthquakes in other parts of Japan and in other countries. The Nankai trough is an example of areas that require our attention. A national project for observation and study for earthquake prediction is now being integrated into a new program, Earthquake and Volcano Hazards Observation and Research Program (2014-2019). Studies presented in this special issue are also being supported in part by this program. We are certain that readers will find that this special issue will contribute much to our understanding of gigantic earthquakes and at least some of the measure to be taken in preparation for such natural phenomena. Finally, we extend our sincere thanks to all of the contributors and reviewers involved with these articles.

階層ベイズモデルによる津波波源逆解析解の信用区間と観測誤差の同時定量推定

階層ベイズモデルによる津波波源逆解析解の信用区間と観測誤差の同時定量推定

OWC波力発電ブイのエネルギー吸収性能に関する基礎的研究

The omni-directional OWC wave power generation device is proposed towards installation of the tsunami observation buoy which can autonomously generate electricity and employ. The applicability and wave energy absorption performance of the OWC buoy are examined by a water tank experiment. The possibility of the energy absorption performance enhancement of the OWC wave power generation buoy are examined by making an experiment on two kinds of OWC buoy models of the bottom opening type and the side opening type. This study shows the basic energy absorption performance and the possibility of the energy absorption performance enhancement of the OWC buoy.

High-frequency radar observations of the June 2013 US East Coast meteotsunami

We report here on the observation and offline detection of the meteotsunami off the New Jersey coast on June 13, 2013, using coastal radar systems and tide gauges. This work extends the previous observations of tsunamis originating in Japan and Indonesia. The radars observed the meteotsunami 23 km offshore, 47 min before it arrived at the coast. Subsequent observations showed it moving onshore. The neighboring tide gauge height reading provides confirmation of the radar observations near the shore.

Forecasting tsunamis in Poverty Bay, New Zealand, with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.

Impact of Near-Field, Deep-Ocean Tsunami Observations on Forecasting the 7 December 2012 Japanese Tsunami

Following the devastating 11 March 2011 tsunami, two deep-ocean assessment and reporting of tsunamis (DART®)(DART® and the DART® logo are registered trademarks of the National Oceanic and Atmospheric Administration, used with permission) stations were deployed in Japanese waters by the Japanese Meteorological Agency. Two weeks after deployment, on 7 December 2012, a M w 7.3 earthquake off Japan’s Pacific coastline generated a tsunami. The tsunami was recorded at the two Japanese DARTs as early as 11 min after the earthquake origin time, which set a record as the fastest tsunami detecting time at a DART station. These data, along with those recorded at other DARTs, were used to derive a tsunami source using the National Oceanic and Atmospheric Administration tsunami forecast system. The results of our analysis show that data provided by the two near-field Japanese DARTs can not only improve the forecast speed but also the forecast accuracy at the Japanese tide gauge stations. This study provides important guidelines for early detection and forecasting of local tsunamis.

Reverberations on the Watery Element: A Significant, Tsunamigenic Historical Earthquake Offshore the Carolina Coast

ABSTRACT We investigate an early nineteenth‐century earthquake that has been previously cataloged but not previously investigated in detail or recognized as a significant event. The earthquake struck at approximately 4:30 a.m. LT on 8 January 1817 and was widely felt throughout the southeastern and mid‐Atlantic United States. Around 11:00 a.m. the same day, an eyewitness described a 12‐inch tide that rose abruptly and agitated boats on the Delaware River near Philadelphia. We show that the timing of this tide is consistent with the predicted travel time for a tsunami generated by an offshore earthquake 6–7 hours earlier. By combining constraints provided by the shaking intensity distribution and the tsunami observation, we conclude that the 1817 earthquake had a magnitude of low‐ to mid‐ M 7 and a location 800–1000 km offshore of South Carolina. Our results suggest that poorly understood offshore source zones might represent a previously unrecognized hazard to the southern and mid‐Atlantic coast. Both observational and modeling results indicate that potential tsunami hazard within Delaware Bay merits consideration: the simple geometry of the bay appears to catch and focus tsunami waves. Our preferred location for the 1817 earthquake is along a diffuse northeast‐trending zone defined by instrumentally recorded and historical earthquakes. The seismotectonic framework for this region remains enigmatic.

Relationship Between Maximum Tsunami Amplitude and Duration of Signal

All available tsunami observations at tide gauges situated along the North American coast were examined to determine if there is any clear relationship between maximum amplitude and signal duration. In total, 89 historical tsunami recordings generated by 13 major earthquakes between 1952 and 2011 were investigated. Tidal variations were filtered out of the signal and the duration between the arrival time and the time at which the signals drops and stays below 0.3 m amplitude was computed. The processed tsunami time series were evaluated and a linear least-squares fit with a 95 % confidence interval was examined to compare tsunami durations with maximum tsunami amplitude in the study region. The confidence interval is roughly 20 h over the range of maximum tsunami amplitudes in which we are interested. This relatively large confidence interval likely results from variations in local resonance effects, late-arriving reflections, and other effects.

Recent Progress of Seismic Observation Networks in Japan

Before the occurrence of disastrous Kobe earthquake in 1995, the number of high sensitivity seismograph stations operated in Japan was nearly 550 and was concentrated in the Kanto and Tokai districts, central Japan. In the wake of the Kobe earthquake, Japanese government has newly established the Headquarters for Earthquake Research Promotion and started the reconstruction of seismic networks to evenly cover the whole Japan. The basic network is composed of three seismographs, i.e. high sensitivity seismograph (Hi-net), broadband seismograph (F-net), and strong motion seismograph (K-NET). A large majority of Hi-net stations are also equipped with a pair of strong motion sensors at the bottom of borehole and the ground surface (KiK-net). A plenty of high quality data obtained from these networks are circulated at once and is producing several new seismological findings as well as providing the basis for the Earthquake Early Warning system. In March 11, 2011, "Off the Pacific coast of Tohoku Earthquake" was generated with magnitude 9.0, which records the largest in the history of seismic observation in Japan. The greatest disaster on record was brought by huge tsunami with nearly 20 thousand killed or missing people. We are again noticed that seismic observation system is quite poor in the oceanic region compared to the richness of it in the inland region. In 2012, NIED has started the construction of ocean bottom seismic and tsunami observation network along the Japan Trench. It is planned to layout 154 stations with an average spacing of 30km, each of which is equipped with an accelerometer for seismic observation and a water pressure gauge for tsunami observation. We are expecting that more rapid and accurate warning of earthquake and tsunami becomes possible by this observing network.