Consistent estimation of strain-rate fields from GNSS velocity data using basis function expansion with ABIC

• Main Authors: Tomohisa Okazaki, Yukitoshi Fukahata, Takuya Nishimura
• Format: Article
• Language: English
• Published: SpringerOpen 2021-07-01
• Series: Earth, Planets and Space
• Subjects: Velocity field; Strain rate; GNSS; ABIC; Inversion analysis
• Online Access: https://doi.org/10.1186/s40623-021-01474-5

Abstract Present day crustal displacement rates can be accurately observed at stations of global navigation satellite system (GNSS), and crustal deformation has been investigated by estimating strain-rate fields from discrete GNSS data. For this purpose, a modified least-square inversion method was proposed by Shen et al. (J Geophys Res 101:27957–27980, 1996). This method offers a simple formulation for simultaneously estimating smooth velocity and strain-rate fields from GNSS data, and it has contributed to clarify crustal deformation fields in many regions all over the world. However, we notice three theoretical points to be examined when we apply the method: mathematical inconsistency between estimated velocity and strain-rate fields, difficulty in objectively determining the optimal value of a hyperparameter that controls smoothness, and inappropriate estimation of uncertainty. In this study, we propose a method of basis function expansion with Akaike’s Bayesian information criterion (ABIC), which overcomes the above difficulties. Application of the two methods to GNSS data in Japan reveals that the inconsistency in the method of Shen et al. is generally insignificant, but could be clear in regions with sparser observation stations such as in islet areas. The method of basis function expansion with ABIC shows a significantly better performance than the method of Shen et al. in terms of the trade-off curve between the residual of fitting and the roughness of velocity field. The estimated strain-rate field with the basis function expansion clearly exhibits a low strain-rate zone in the forearc from the southern Tohoku district to central Japan. We also find that the Ou Backbone Range has several contractive spots around active volcanoes and that these locations well correspond to the subsidence areas detected by InSAR after the 2011 Tohoku-oki earthquake. Thus, the method of basis function expansion with ABIC would serve as an effective tool for estimating strain-rate fields from GNSS data.