Search

We investigate the postseismic deformation of the 2004 Sumatra–Andaman earthquake (SAE) using 5 years of Global Positioning System (GPS) data located in northern Sumatra. Continuous GPS data from northern Sumatra suggest that the relaxation time in the vertical displacement is longer than horizontal displacements. This implies that there are multiple physical mechanisms that control the postseismic deformation, which refer to afterslip and viscoelastic relaxation. In this study, we introduce an analysis strategy of postseismic deformation to simultaneously calculate multiple mechanisms of afterslip and viscoelastic relaxation. The afterslip inversion results indicate that the distribution of the afterslip and the coseismic slip are compensatory of each other. Also, afterslip has a limited contribution to vertical deformation in northern Sumatra. In our rheology model, we use a gravitational Maxwell viscoelastic response and the result indicates that the elastic layer thickness is 65 ± 5 km and the Maxwell viscosity is 8.0 ± 1.0 × 1018 Pa s. We find that afterslip plus Maxwell viscoelastic relaxation are appropriate to explain the deformation in northern Sumatra. We also find that our rheology model reproduces the long-term features of the GPS time series in Thailand. Applying our rheology model to the data in Andaman Islands our afterslip estimation is located at the down-dip part of the plate boundary. Finally, we showed that our rheology model is applicable to the GPS datasets of postseismic deformation of the 2004 SAE located in northern Sumatra, Thailand, and Andaman–Nicobar, respectively.

© 2016 The Author(s). The concept of a weak asthenospheric layer underlying Earth's mobile tectonic plates is fundamental to our understanding of mantle convection and plate tectonics. However, little is known about the mechanical properties of the asthenosphere (the part of the upper mantle below the lithosphere) underlying the oceanic crust, which covers about 60 per cent of Earth's surface. Great earthquakes cause large coseismic crustal deformation in areas hundreds of kilometres away from and below the rupture area. Subsequent relaxation of the earthquake-induced stresses in the viscoelastic upper mantle leads to prolonged postseismic crustal deformation that may last several decades and can be recorded with geodetic methods. The observed postseismic deformation helps us to understand the rheological properties of the upper mantle, but so far such measurements have been limited to continental-plate boundary zones. Here we consider the postseismic deformation of the very large (moment magnitude 8.6) 2012 Indian Ocean earthquake to provide by far the most direct constraint on the structure of oceanic mantle rheology. In the first three years after the Indian Ocean earthquake, 37 continuous Global Navigation Satellite Systems stations in the region underwent horizontal northeastward displacements of up to 17 centimetres in a direction similar to that of the coseismic offsets. However, a few stations close to the rupture area that had experienced subsidence of up to about 4 centimetres during the earthquake rose by nearly 7 centimetres after the earthquake. Our three-dimensional viscoelastic finite-element models of the post-earthquake deformation show that a thin (30-200 kilometres), low-viscosity (having a steady-state Maxwell viscosity of (0.5-10) × 10 18 pascal seconds) asthenospheric layer beneath the elastic oceanic lithosphere is required to produce the observed postseismic uplift.

The Philippine Sea plate has subducted beneath southwest Japan at the Nankai Trough, causing megathrust earthquakes every 100-140 years. The subduction process shows lateral variations along the trough, such as a bending of the plate interface beneath the Kii Peninsula and rapid changes of strike and dip off the east coast of Kyushu. We invert three-dimensional GPS velocities using ABIC geodetic inversion technique to estimate interseismic slip de.cit distribution at the Nankai Trough. The studied area is roughly from Kii peninsula to eastern Kyushu (131.5-137E, 31-35N). In the ABIC geodetic inversion some initial constraints are needed to stabilize solutions, such that the slip at the outer margin of the fault is set to zero. For the estimation of interseismic slip de.cit on the plate interface, however, these constraints are unrealistic because slip de.cit distribution continues later.ally to the outside of the target area. In this study we propose a repeated inversion method that shifts target area laterally along the plate boundary. This method averages slip de.cit estimates at the same point obtained from more than three consecutive inversions. The estimates obtained near the lateral margin, which may be strongly affected by the non-slip constraints at the outer margin, are excluded from the averaging. As a result final solutions are stabilized regardless the setting of the region. The final slip de.cit distribution at the Nankai Trough shows high contrast of strong and weak plate locking zones. Down-dip variation of the slip de.cit is emphasized more clearly but the lateral variation becomes moderate compared with the single inversion result.

We have studied crustal deformation in the Nankai subduction zone, southwest Japan, based on three-dimensional GPS velocity fields. Oblique subduction of the Philippine Sea plate has caused two different modes of deformation of the overriding plate: interseismic crustal shortening in the direction of plate convergence, and permanent lateral movement of the forearc. The block boundary dividing the forearc is the Median Tectonic Line (MTL); however, we assumed that its shallower portion is fully or partially locked to a certain depth. The plate boundary and the MTL are represented by many rectangular faults. We carried out inversion analyses with a priori information to estimate simultaneously slip deficit rates at those rectangular faults, together with the rate of lateral movement of the forearc. The results show that the seismogenic subduction faults at a depth of 5–25 km are strongly locked. As for the transition zone at 25–35 km, the inversion analysis results in stronger locking than shown in previous studies, especially when the vertical velocity data are weighted. The rates of lateral forearc movement and slip deficit at the MTL are nearly comparable but in a reverse sense to each other. This shows that the shallower portion of the MTL is strongly locked but that stationary aseismic slip is occurring in the deeper part.

The interseismic strain accumulation process at the Nankai subduction zone, the plate-boundary region between the Philippine Sea and the Eurasian plates off Southwest Japan, has been revealed by GPS measurements spanning five years. Large strains characterized by a contraction of 2.2 ∼ 3.4 × 10−7/yr have accumulated in a narrow region extending parallel to the strike of the plate boundary. The contraction axes lie in a northwest-southeast direction nearly parallel to the direction of relative motion of the Philippine Sea plate with respect to the Eurasian plate. GPS-derived deformation data, including both rate and direction, are consistent with the expected earthquake deformation cycle at the Nankai subduction zone. Current interplate coupling in this zone is strong; much of the relative plate motion has been taken up by deformation of the marginal parts of the Eurasian plate.

Although the Kozu Island Volcano, one of the Izu Islands Volcanoes in the south part Central Japan, is an active volcano, there is no record of the eruption for about 1100 years since the last eruption in 833 A.D. Since 1988, frequent earthquake swarms are observed around the Kozu Island, and the uplift of 2–4 cm/yr is observed on the island by tidal observations. Station velocities detected by GPS measurements since 1989 show velocities that differ from the convergent velocity of the Philippine Sea plate calculated from plate motion models. A local GPS network with 12 stations is occupied around the volcano, and the GPS measurements are repeated every about six month since July 1996. Inflated deformation of 2–4 cm/yr are detected from the GPS measurements and the pressure source is estimated to be located in the northeastern part of the island at a depth of 2.1 km using Mogi solution. Negative gravity changes of more than 30 microgal are also measured above the pressure source in the period November 1998 to July 1999, consistent with uplift.