Spatiotemporal Characteristics of Horizontal Crustal Deformation in the Sichuan–Yunnan Region Using GPS Data

Based on various velocity fields from Global Positioning System (GPS) data over nearly 20 years in the Sichuan–Yunnan region, this paper calculated the strain rate field and its spatiotemporal characteristics by using an improved least squares collocation method. We evaluated the calculated strain field by extensively discussing the impact of non-tectonic factors on the calculation. Subsequently, we described the present-day strain rate features and their spatiotemporal variations. The results indicate the necessity of considering the influence of non-tectonic factors when calculating the strain rate field by using GPS velocity data. Widespread strain accumulation is observed in the Sichuan–Yunnan region, and significant second strain rate invariant with an average value of 33.1 nanostrain/yr primarily occurs along the eastern boundary faults of the Sichuan–Yunnan rhomboid block, specifically the Xianshuihe–Anninghe–Zemuhe–Xiaojiang fault systems. These fault systems also demarcate the zones with negative and positive dilation strain. According to the spatiotemporal variations of strain rate fields, the northern and southern segments of the Xianshuihe Fault, the Anninghe Fault and its eastern adjacent faults, and the Xiaojiang Fault are undergoing intensifying strain. Consequently, these zones should be paid more attention due to their relatively higher seismic risk in the Sichuan–Yunnan region.