Secular crustal deformation and interplate coupling of the Japanese Islands as deduced from continuous GPS array, 1996-2001

Data from the nation-wide GPS continuous tracking network that has been operated by the Geographical Survey Institute of Japan since April 1996 were used to study crustal deformation in the Japanese Islands. We first extracted site coordinate from daily SINEX files for the period from April 1, 1996 to February 24, 2001. Since raw time series of station coordinates include coseismic and postseismic displacements as well as seasonal variation, we model each time series as a combination of linear and trigonometric functions and jumps for episodic events. Estimated velocities were converted into a kinematic reference frame [Heki, K., 1996. Horizontal and vertical crustal movements from three-dimensional very long baseline interferometry kinematic reference frame: implication for reversal timescale revision. J. Geophys. Res., 101: 3187-3198.] to discuss the crustal deformation relative to the stable interior of the Eurasian plate. A Least-Squares Prediction technique has been used to segregate the signal and noise in horizontal as well as vertical velocities. Estimated horizontal signals (horizontal displacement rates) were then differentiated in space to calculate principal components of strain. Dilatations, maximum shear strains, and principal axes of strain clearly portray tectonic environments of the Japanese Islands. On the other hand, the interseismic vertical deformation field of the Japanese islands is derived for the same GPS data interval. The GPS vertical velocities are combined with 31 year tide gage records to estimate absolute vertical velocity. The results of vertical deformation show that (1) the existence of clear uplift of about 6 mm/yr in Shikoku and Kii Peninsula, whereas pattern of subsidence is observed in the coast of Kyushu district. This might reflect strong coupling between the Philippine Sea plate and overriding plate at the Nankai Trough and weak coupling off Kyshu, (2) no clear vertical deformation pattern exists along the Pacific coast of northeastern Japan. This might be due to the long distance between the plate boundary (Japan trench) and overriding plate where GPS sites are located, (3) significant uplift is observed in the southwestern part of Hokkaido and in northeastern Tohoku along the Japan Sea coast. This is possibly due to the viscoelastic rebound of the 1983 Japan Sea (M-w 7.7) and the 1993 Hokkaido-Nansei-Oki (M-w 7.8) earthquakes and/or associated with distributed compression of incipient subduction there. We then estimate the elastic deformation of the Japanese Islands caused by interseismic loading of the Pacific and Philippine Sea subduction plates. The elastic models account for most of the observed horizontal velocity field if the subduction movement of the Philippine Sea Plate is 100% locked and if that of the Pacific Plate is 70% locked. However, the best fit for vertical velocity ranges from 80% to 100% coupling factor in southwestern Japan and only 50% in northeastern Japan. Since horizontal data does not permit the separation of rigid plate motion and interplate coupling because horizontal velocities include both contributions, we used the vertical velocities to discriminate between them. So, we can say there is strong interplate coupling (80%-100%) over the Nankaido subduction zone, whereas it is about 50% only over the Kurile-Japan trench. (c) 2006 Elsevier B.V. All rights reserved.