Tsunami early warning using GPS‐Shield arrays

[1] The 2004 catastrophic Indian Ocean tsunami has strongly emphasized the need for reliable tsunami early warning systems. Another giant tsunamigenic earthquake may occur west of Sumatra, close to the large city of Padang. We demonstrate that the presence of islands between the trench and the Sumatran coast makes earthquake-induced tsunamis especially sensitive to slip distribution on the rupture plane as wave heights at Padang may differ by more than a factor of 5 for earthquakes having the same seismic moment (magnitude) and rupture zone geometry but different slip distribution. Hence reliable prediction of tsunami wave heights for Padang cannot be provided using traditional, earthquake-magnitude-based methods. We show, however, that such a prediction can be issued within 10 minutes of an earthquake by incorporating special types of near-field GPS arrays (“GPS-Shield”). These arrays measure both vertical and horizontal displacements and can resolve higher order features of the slip distribution on the fault than the seismic moment if placed above the rupture zone or are less than 100 km away of the rupture zone. Stations in the arrays are located as close as possible to the trench and are aligned perpendicular to the trench, i.e., parallel to the expected gradient of surface coseismic displacement. In the case of Sumatra and Java, the GPS-Shield arrays should be placed at Mentawai Islands, located between the trench and Sumatra and directly at the Sumatra and Java western coasts. We demonstrate that the “GPS-Shield” can also be applied to northern Chile, where giant earthquakes may also occur in the near future. Moreover, this concept may be applied globally to many other tsunamigenic active margins where the land is located above or close to seismogenic zones.

[1]  Emile A. Okal,et al.  Seismology: Speed and size of the Sumatra earthquake , 2005, Nature.

[2]  Eric L. Geist,et al.  Local Tsunamis and Earthquake Source Parameters , 1998 .

[3]  K. Sieh,et al.  Deformation and Slip Along the Sunda Megathrust in the Great 2005 Nias-Simeulue Earthquake , 2006, Science.

[4]  Chen Ji,et al.  Source Description of the 1999 Hector Mine, California, Earthquake, Part I: Wavelet Domain Inversion Theory and Resolution Analysis , 2002 .

[5]  Hans-Peter Plag,et al.  Rapid determination of earthquake magnitude using GPS for tsunami warning systems , 2006 .

[6]  Yehuda Bock,et al.  Detection of arbitrarily large dynamic ground motions with a dense high‐rate GPS network , 2004 .

[7]  J. Brune Tectonic stress and the spectra of seismic shear waves from earthquakes , 1970 .

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

[9]  J. Zumberge,et al.  Precise point positioning for the efficient and robust analysis of GPS data from large networks , 1997 .

[10]  R. Dmowska,et al.  Local Tsunamis and Distributed Slip at the Source , 1999 .

[11]  F. Pollitz,et al.  Stress changes along the Sunda trench following the 26 December 2004 Sumatra‐Andaman and 28 March 2005 Nias earthquakes , 2006 .

[12]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[13]  Kenneth W. Hudnut,et al.  Uplift and subsidence associated with the great Aceh‐Andaman earthquake of 2004 , 2006 .

[14]  Steven N. Ward,et al.  The 1960 Chile earthquake: inversion for slip distribution from surface deformation , 1990 .

[15]  Peter M. Shearer,et al.  Extent, duration and speed of the 2004 Sumatra–Andaman earthquake imaged by the Hi-Net array , 2005, Nature.

[16]  Michael R. Craymer,et al.  Current tectonics of northern Cascadia from a decade of GPS measurements , 2003 .

[17]  Yehuda Bock,et al.  Instantaneous geodetic positioning at medium distances with the Global Positioning System , 2000 .

[18]  R. Lawrence Edwards,et al.  Submergence and uplift associated with the giant 1833 Sumatran subduction earthquake: Evidence from coral microatolls , 1999 .

[19]  E. Okal Seismic parameters controlling far-field tsunami amplitudes: A review , 1988 .

[20]  H. Kanamori,et al.  A slow earthquake , 1979 .

[21]  Ragnar Stefánsson,et al.  Nucleation Points of Recent Mainshocks in Southern Iceland, Mapped by b-Values , 2006 .

[22]  Rongjiang Wang,et al.  PSGRN/PSCMP - a new code for calculating co- and post-seismic deformation, geoid and gravity changes based on the viscoelastic-gravitational dislocation theory , 2006, Comput. Geosci..

[23]  Chalermchon Satirapod,et al.  Insight into the 2004 Sumatra–Andaman earthquake from GPS measurements in southeast Asia , 2005, Nature.

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

[25]  Costas E. Synolakis,et al.  Tsunami inundation modeling for western Sumatra , 2006, Proceedings of the National Academy of Sciences.

[26]  Rongjiang Wang,et al.  A simple orthonormalization method for stable and efficient computation of Green's functions , 1999, Bulletin of the Seismological Society of America.

[27]  F. Pollitz,et al.  The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets , 2005, Science.

[28]  Chen Ji,et al.  Coseismic Slip and Afterslip of the Great Mw 9.15 Sumatra–Andaman Earthquake of 2004 , 2007 .

[29]  Eric Deleersnijder,et al.  An efficient Eulerian finite element method for the shallow water equations , 2005 .

[30]  John McCloskey,et al.  Indonesian earthquake: Earthquake risk from co-seismic stress , 2005, Nature.

[31]  L. B. Freund,et al.  A two-dimensional analysis of surface deformation due to dip-slip faulting , 1976 .

[32]  D. Boore Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra , 1983 .

[33]  Kurt L. Feigl,et al.  Crustal deformation and fault slip during the seismic cycle in the North Chile subduction zone, from GPS and InSAR observations , 2004 .

[34]  Yehuda Bock,et al.  Plate-boundary deformation associated with the great Sumatra–Andaman earthquake , 2006, Nature.

[35]  S. Jade,et al.  Co-seismic and post-seismic displacements in Andaman and Nicobar Islands from GPS measurements , 2005 .

[36]  Chen Ji,et al.  Rupture Process of the 2004 Sumatra-Andaman Earthquake , 2005, Science.

[37]  V. Gahalaut,et al.  Constraints on 2004 Sumatra Andaman earthquake rupture from GPS measurements in Andaman Nicobar Islands , 2006 .

[38]  M. Rothacher,et al.  Towards real‐time tsunami amplitude prediction , 2006 .

[39]  R. Dmowska,et al.  Seismic activity in the Sumatra–Java region prior to the December 26, 2004 (Mw = 9.0–9.3) and March 28, 2005 (Mw = 8.7) earthquakes , 2006 .

[40]  J. Clague,et al.  Giant Earthquakes Beneath Canada's West Coast , 1996 .

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

[42]  V. Gahalaut,et al.  Constraints on rupture of the December 26, 2004, Sumatra earthquake from far-field GPS observations , 2005 .

[43]  S. Tinti,et al.  A smoothing algorithm to enhance finite-element tsunami modelling: An application to the 5 February 1783 Calabrian case, Italy , 1995 .

[44]  Y. Okada Surface deformation due to shear and tensile faults in a half-space , 1985 .

[45]  Kojiro Irikura,et al.  Semi-Empirical Estimation of Strong Ground Motions During Large Earthquakes , 1983 .

[46]  Matthias Ohrnberger,et al.  Tracking the rupture of the Mw = 9.3 Sumatra earthquake over 1,150 km at teleseismic distance , 2005, Nature.