Tsunami Data Assimilation Without a Dense Observation Network

The tsunami data assimilation method enables tsunami forecasting directly from observations, without the need of estimating tsunami sources. However, it requires a dense observation network to produce desirable results. Here we propose a modified method of tsunami data assimilation for regions with a sparse observation network. The method utilizes interpolated waveforms at virtual stations. The tsunami waveforms at the virtual stations between two existing observation stations are estimated by shifting arrival times with the linear interpolation of observed arrival times and by correcting the amplitudes for their water depths. In our new data assimilation approach, we employ the Optimal Interpolation algorithm to both the real observations and virtual stations, in order to construct a complete wavefront of tsunami propagation. The application to the 2004 Sumatra‐Andaman earthquake and the 2009 Dusky Sound, New Zealand, earthquake reveals that addition of virtual stations greatly helps improve the tsunami forecasting accuracy. Plain Language Summary Data assimilation is a method to combine observation and numerical simulation and is widely used in weather forecast. The data assimilation methods have been recently applied for tsunami forecast in North America and Japan where dense observation networks exist. In this study, we proposed a data assimilation method by introducing virtual observation data from neighboring real observations. We applied the method for the Indian Ocean with the 2004 Sumatra‐Andaman earthquake tsunami and offshore New Zealand with the 2009 Dusky Sound earthquake tsunami. We found that the method greatly improved the forecasting accuracy and the method could be used for the regions with sparse observation network.

[1]  Fumihiko Imamura,et al.  Near-field tsunami forecasting from cabled ocean bottom pressure data , 2009 .

[2]  M. Leonard Analysis of tide gauge records from the December 2004 Indian Ocean tsunami , 2006 .

[3]  Kenji Satake,et al.  Tsunami Source of the 2004 Sumatra–Andaman Earthquake Inferred from Tide Gauge and Satellite Data , 2007 .

[4]  Takuto Maeda,et al.  Successive estimation of a tsunami wavefield without earthquake source data: A data assimilation approach toward real‐time tsunami forecasting , 2015 .

[5]  S. Samsonov,et al.  Oblique slip on the Puysegur subduction interface in the 2009 July MW 7.8 Dusky Sound earthquake from GPS and InSAR observations: implications for the tectonics of southwestern New Zealand , 2010 .

[6]  K. Furlong,et al.  Intraplate deformation adjacent to the Macquarie Ridge south of New Zealand—The tectonic evolution of a complex plate boundary , 2007 .

[7]  J. Tromp,et al.  The Great Sumatra-Andaman Earthquake , 2005 .

[8]  T. Ishibe,et al.  Comparative study of two tsunamigenic earthquakes in the Solomon Islands: 2015 Mw 7.0 normal‐fault and 2013 Santa Cruz Mw 8.0 megathrust earthquakes , 2016 .

[9]  John McCloskey,et al.  Seismology: Earthquake risk on the Sunda trench , 2005, Nature.

[10]  Toshitaka Baba,et al.  Tsunami magnitudes determined from ocean‐bottom pressure gauge data around Japan , 2004 .

[11]  N. Litchfield,et al.  The 2009 Mw 7.8 earthquake on the Puysegur subduction zone produced minimal geological effects around Dusky Sound, New Zealand , 2011 .

[12]  Array Observations of the 2012 Haida Gwaii Tsunami Using Cascadia Initiative Absolute and Differential Seafloor Pressure Gauges , 2015 .

[13]  H. Kanamori,et al.  The Great Sumatra-Andaman Earthquake of 26 December 2004 , 2005, Science.

[14]  Daisuke Matsuoka,et al.  Parallel Implementation of Dispersive Tsunami Wave Modeling with a Nesting Algorithm for the 2011 Tohoku Tsunami , 2015, Pure and Applied Geophysics.

[15]  Iyan E. Mulia,et al.  Tsunami data assimilation of Cascadia seafloor pressure gauge records from the 2012 Haida Gwaii earthquake , 2016 .

[16]  G. Laske,et al.  The character of seafloor ambient noise recorded offshore New Zealand: Results from the MOANA ocean bottom seismic experiment , 2012 .

[17]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[18]  Harold O. Mofjeld,et al.  The NTHMP Tsunameter Network , 2005 .

[19]  Isamu Aida,et al.  RELIABILITY OF A TSUNAMI SOURCE MODEL DERIVED FROM FAULT PARAMETERS , 1978 .

[20]  M. Eblé,et al.  Deep-Ocean Measurements of Tsunami Waves , 2015, Pure and Applied Geophysics.

[21]  Philip L.-F. Liu,et al.  Tsunami hazard and early warning system in South China Sea , 2009 .

[22]  B. Bradley,et al.  THE MW 7.6 DUSKY SOUND EARTHQUAKE OF 2009: PRELIMINARY REPORT , 2010 .

[23]  Vasily V. Titov,et al.  Real-Time Tsunami Forecasting: Challenges and Solutions , 2003 .

[24]  K. Ford,et al.  Tools and Techniques: Radiometric Methods , 2015 .

[25]  Kenji Satake,et al.  Green's Function‐Based Tsunami Data Assimilation: A Fast Data Assimilation Approach Toward Tsunami Early Warning , 2017 .

[26]  Gavin P. Hayes,et al.  Quantifying potential tsunami hazard in the Puysegur subduction zone, south of New Zealand , 2010 .

[27]  J. Dettmer,et al.  Time reverse imaging for far‐field tsunami forecasting: 2011 Tohoku earthquake case study , 2015 .

[28]  Kenji Satake,et al.  Linear and nonlinear computations of the 1992 Nicaragua earthquake tsunami , 1995 .

[29]  T. Waseda,et al.  Preparing for the Future Nankai Trough Tsunami: A Data Assimilation and Inversion Analysis From Various Observational Systems , 2017 .

[30]  Aditya Riadi Gusman,et al.  A methodology for near‐field tsunami inundation forecasting: Application to the 2011 Tohoku tsunami , 2014 .

[31]  Le principe de Huygens , 1924 .

[32]  A. Gusman,et al.  Improving Forecast Accuracy With Tsunami Data Assimilation: The 2009 Dusky Sound, New Zealand, Tsunami , 2019, Journal of Geophysical Research: Solid Earth.

[33]  Eric M. Dunham,et al.  Tsunami Wavefield Reconstruction and Forecasting Using the Ensemble Kalman Filter , 2019, Geophysical Research Letters.