Near real-time GPS applications for tsunami early warning systems

GPS (Global Positioning System) technology is widely used for positioning applications. Many of them have high requirements with respect to precision, reliability or fast product delivery, but usually not all at the same time as it is the case for early warning applications. The tasks for the GPS-based components within the GITEWS project (German Indonesian Tsunami Early Warning System, Rudloff et al., 2009) are to support the determination of sea levels (measured onshore and offshore) and to detect co-seismic land mass displacements with the lowest possible latency (design goal: first reliable results after 5 min). The completed system was designed to fulfil these tasks in near real-time, rather than for scientific research requirements. The obtained data products (movements of GPS antennas) are supporting the warning process in different ways. The measurements from GPS instruments on buoys allow the earliest possible detection or confirmation of tsunami waves on the ocean. Onshore GPS measurements are made collocated with tide gauges or seismological stations and give information about co-seismic land mass movements as recorded, e.g., during the great Sumatra-Andaman earthquake of 2004 (Subarya et al., 2006). This information is important to separate tsunami-caused sea height movements from apparent sea height changes at tide gauge locations (sensor station movement) and also as additional information about earthquakes' mechanisms, as this is an essential information to predict a tsunami (Sobolev et al., 2007). <br><br> This article gives an end-to-end overview of the GITEWS GPS-component system, from the GPS sensors (GPS receiver with GPS antenna and auxiliary systems, either onshore or offshore) to the early warning centre displays. We describe how the GPS sensors have been installed, how they are operated and the methods used to collect, transfer and process the GPS data in near real-time. This includes the sensor system design, the communication system layout with real-time data streaming, the data processing strategy and the final products of the GPS-based early warning system components

Partial rupture of a locked patch of the Sumatra megathrust during the 2007 earthquake sequence

The great Sumatra–Andaman earthquake and tsunami of 2004 was a dramatic reminder of the importance of understanding the seismic and tsunami hazards of subduction zones [1,2,3,4]. In March 2005, the Sunda megathrust ruptured again, producing an event [5] of moment magnitude (Mw) 8.6 south of the 2004 rupture area, which was the site of a similar event in 1861 (ref. 6). Concern was then focused on the Mentawai area, where large earthquakes had occurred in 1797 (Mw = 8.8) and 1833 (Mw = 9.0) [6,7]. Two earthquakes, one of Mw = 8.4 and, twelve hours later, one of Mw = 7.9, indeed occurred there on 12 September 2007. Here we show that these earthquakes ruptured only a fraction of the area ruptured in 1833 and consist of distinct asperities within a patch of the megathrust that had remained locked in the interseismic period. This indicates that the same portion of a megathrust can rupture in different patterns depending on whether asperities break as isolated seismic events or cooperate to produce a larger rupture. This variability probably arises from the influence of non-permanent barriers, zones with locally lower pre-stress due to the past earthquakes. The stress state of the portion of the Sunda megathrust that had ruptured in 1833 and 1797 was probably not adequate for the development of a single large rupture in 2007. The moment released in 2007 amounts to only a fraction both of that released in 1833 and of the deficit of moment that had accumulated as a result of interseismic strain since 1833. The potential for a large megathrust event in the Mentawai area thus remains large

Plate-boundary deformation associated with the great Sumatra–Andaman earthquake

The Sumatra–Andaman earthquake of 26 December 2004 is the first giant earthquake (moment magnitude M_w > 9.0) to have occurred since the advent of modern space-based geodesy and broadband seismology. It therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events. Here we report estimates of the ground displacement associated with this event, using near-field Global Positioning System (GPS) surveys in northwestern Sumatra combined with in situ and remote observations of the vertical motion of coral reefs. These data show that the earthquake was generated by rupture of the Sunda subduction megathrust over a distance of >1,500 kilometres and a width of <150 kilometres. Megathrust slip exceeded 20 metres offshore northern Sumatra, mostly at depths shallower than 30 kilometres. Comparison of the geodetically and seismically inferred slip distribution indicates that ~30 per cent additional fault slip accrued in the 1.5 months following the 500-second-long seismic rupture. Both seismic and aseismic slip before our re-occupation of GPS sites occurred on the shallow portion of the megathrust, where the large Aceh tsunami originated. Slip tapers off abruptly along strike beneath Simeulue Island at the southeastern edge of the rupture, where the earthquake nucleated and where an M_w = 7.2 earthquake occurred in late 2002. This edge also abuts the northern limit of slip in the 28 March 2005 M_w = 8.7 Nias–Simeulue earthquake

Receiver function study in northern Sumatra and the Malaysian peninsula

International audienceIn this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities () near the surface in the range of 3.4-3.6 km s attributed to crystalline rocks and 3.6-4.0 km s in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30-34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low of 2.6-2.9 km s indicate sediment layers. High of 4.2 km s are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in between 10 and 20 km depth. Within the subducting slab ≈ 4.7 km s. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found

Tsunami inundation modeling for western Sumatra

A long section of the Sunda megathrust south of the great tsunamigenic earthquakes of 2004 and 2005 is well advanced in its seismic cycle and a plausible candidate for rupture in the next few decades. Our computations of tsunami propagation and inundation yield model flow depths and inundations consistent with sparse historical accounts for the last great earthquakes there, in 1797 and 1833. Numerical model results from plausible future ruptures produce flow depths of several meters and inundation up to several kilometers inland near the most populous coastal cities. Our models of historical and future tsunamis confirm a substantial exposure of coastal Sumatran communities to tsunami surges. Potential losses could be as great as those that occurred in Aceh in 2004

Insight into the 2004 Sumatra-Andaman earthquake from GPS measurements in southeast Asia

Data collected at ~60 Global Positioning System (GPS) sites in southeast Asia show the crustal deformation caused by the 26 December 2004 Sumatra-Andaman earthquake at an unprecedented large scale. Small but significant co-seismic jumps are clearly detected more than 3,000 km from the earthquake epicentre. The nearest sites, still more than 400 km away, show displacements of 10 cm or more. Here we show that the rupture plane for this earthquake must have been at least 1,000 km long and that non-homogeneous slip is required to fit the large displacement gradients revealed by the GPS measurements. Our kinematic analysis of the GPS recordings indicates that the centroid of released deformation is located at least 200 km north of the seismological epicentre. It also provides evidence that the rupture propagated northward sufficiently fast for stations in northern Thailand to have reached their final positions less than 10 min after the earthquake, hence ruling out the hypothesis of a silent slow aseismic rupture