No abstract

After a large subduction earthquake, crustal deformation continues to occur, with a complex pattern of evolution. This postseismic deformation is due primarily to viscoelastic relaxation of stresses induced by the earthquake rupture and continuing slip (afterslip) or relocking of different parts of the fault. When postseismic geodetic observations are used to study Earth's rheology and fault behaviour, it is commonly assumed that short-term (a few years) deformation near the rupture zone is caused mainly by afterslip, and that viscoelasticity is important only for longer-term deformation. However, it is difficult to test the validity of this assumption against conventional geodetic data. Here we show that new seafloor GPS (Global Positioning System) observations immediately after the great Tohoku-oki earthquake provide unambiguous evidence for the dominant role of viscoelastic relaxation in short-term postseismic deformation. These data reveal fast landward motion of the trench area, opposing the seaward motion of GPS sites on land. Using numerical models of transient viscoelastic mantle rheology, we demonstrate that the landward motion is a consequence of relaxation of stresses induced by the asymmetric rupture of the thrust earthquake, a process previously unknown because of the lack of near-field observations. Our findings indicate that previous models assuming an elastic Earth will have substantially overestimated afterslip downdip of the rupture zone, and underestimated afterslip updip of the rupture zone; our knowledge of fault friction based on these estimates therefore needs to be revised.

[1] On 11 March 2011, the devastating M9.0 Tohoku Earthquake occurred on the interface of the subducting Pacific plate, and was followed by a huge tsunami that killed about 20,000 people. Several geophysical studies have already suggested that the very shallow portion of the plate interface might have played an important role in producing such a large earthquake and tsunami. However, the sparsity of seafloor observations leads to insufficient spatial resolution of the fault slip on such a shallow plate interface. For this reason, the location and degree of the slip has not yet been estimated accurately enough to assess future seismic risks. Thus, we estimated the coseismic slip distribution based on terrestrial GPS observations and all available seafloor geodetic data that significantly improve the spatial resolution at the shallow portion of the plate interface. The results reveal that an extremely large (greater than 50 m) slip occurred in a small (about 40 km in width and 120 km in length) area near the Japan Trench and generated the huge tsunami. The estimated slip distribution and a comparison of it with the coupling coefficient distribution deduced from the analysis of the small repeating earthquakes suggest that the 2011 Tohoku Earthquake released strain energy that had accumulated over the past 1000 years, probably since the Jogan Earthquake in 869. The accurate assessments of seismic risks on very shallow plate interfaces in subduction zones throughout the world can be obtained by improving the quality and quantity of seafloor geodetic observations.

The northeastern Japan arc is located in a typical subduction zone and is a seismically active region where large interplate earthquakes have occurred repeatedly. The nationwide GPS network has made it possible to investigate the crustal deformation in unprecedented detail; however, vertical displacements are less accurate than the horizontal ones and have not been used to constrain interplate coupling models. In this study, GPS carrier phase data are analyzed in order to estimate three components of site displacements. The result indicates uplift along the Japan Sea coast and subsidence along the Pacific Ocean coast. Using the vertical site velocities together with the horizontal ones as the data for a geodetic inversion method, a new model of interplate coupling is constructed to explain both the horizontal and vertical velocities for the interseismic period from 1997 to 2001. The model demonstrates strong coupling in the area off Miyagi and in the area off Aomori through the area off Tokachi, including the areas of asperities. The downdip limit of the locked fault zone on the plate boundary has been estimated seismologically to be at a depth of 50–60 km in northeastern Japan, and the distribution of slip deficits derived in this study indicates that the transition zone extends down to a depth of about 100 km.