Analysis of Crustal Movement and Deformation in Mainland China Based on CMONOC Baseline Time Series

In this paper, we propose a method for the analysis of tectonic movement and crustal deformation by using GNSS baseline length change rates or baseline linear strain rates. The method is applied to daily coordinate solutions of continuous GNSS stations of the Crustal Movement Observation Network of China (CMONOC). The results show that: (a) The baseline linear strain rates are uneven in space, which is prominent in the Tianshan, Sichuan-Yunnan, Qinghai-Tibet Plateau, and Yanjing areas, with a maximum value of 1 × 10−7 a−1, and about two orders smaller in the South China block, the Northeast block, and the inner area of the Tarim basin, where the average baseline linear strain rates are 1.471 × 10−9 a−1, 2.242 × 10−9 a−1, and 3.056 × 10−9 a−1, respectively; (b) Active crustal deformation and strong earthquakes in the Xinjiang area are mainly located in the north and south sides of the Tianshan block; the compression deformations both inside the Tarim block and in the southern Tianshan fault zone are all increasing from east to west, and the Tarim block is not a completely “rigid block”, with the shrinkage rate in the west part at about 1~2 mm/a; (c) The principal directions of crustal deformation in the Xinjiang, Tibet, and Sichuan-Yunnan regions are generally in the north—south compression and east—west extension, indicating that the collision and wedging between the Indian and Eurasian plates are still the main source of tectonic movements in mainland China.

Independent Component Extraction from the Incomplete Coordinate Time Series of Regional GNSS Networks

Independent component analysis (ICA) is one of the most effective approaches in extracting independent signals from a global navigation satellite system (GNSS) regional station network. However, ICA requires the involved time series to be complete, thereby the missing data of incomplete time series should be interpolated beforehand. In this contribution, a modified ICA is proposed, by which the missing data are first recovered based on the reversible property between the original time series and decomposed principal components, then the complete time series are further processed with FastICA. To evaluate the performance of the modified ICA for extracting independent components, 24 regional GNSS network stations located in North China from 2011 to 2019 were selected. After the trend, annual and semiannual terms were removed from the GNSS time series, the first two independent components captured 17.42, 18.44 and 17.38% of the total energy for the North, East and Up coordinate components, more than those derived by the iterative ICA that accounted for 16.21%, 17.72% and 16.93%, respectively. Therefore, modified ICA can extract more independent signals than iterative ICA. Subsequently, selecting the 7 stations with less missing data from the network, we repeatedly process the time series after randomly deleting parts of the data and compute the root mean square error (RMSE) from the differences of reconstructed signals before and after deleting data. All RMSEs of modified ICA are smaller than those of iterative ICA, indicating that modified ICA can extract more exact signals than iterative ICA.

Transient Deformation in California From Two Decades of GPS Displacements: Implications for a Three-Dimensional Kinematic Reference Frame

Our understanding of plate boundary deformation has been enhanced by transient signals observed against the backdrop of time-independent secular motions. We make use of a new analysis of displacement time series from about 1,000 continuous Global Positioning System (GPS) stations in California from 1999 to 2018 to distinguish tectonic and nontectonic transients from secular motion. A primary objective is to define a high-resolution three-dimensional reference frame (datum) for California that can be rapidly maintained with geodetic data to accommodate both secular and time-dependent motions. To this end, we compare the displacements to those predicted by a horizontal secular fault slip model for the region and construct displacement and strain rate fields. Over the past 19 years, California has experienced 19 geodetically detectable earthquakes and widespread postseismic deformation. We observe postseismic strain rate variations as large as 1,000 nstrain/year with moment releases equivalent up to an Mw6.8 earthquake. We find significant secular differences up to 10 mm/year with the fault slip model, from the Mendocino Triple Junction to the southern Cascadia subduction zone, the northern Basin and Range, and the Santa Barbara channel. Secular vertical uplift is observed across the Transverse Ranges, Coastal Ranges, Sierra Nevada, as well as large-scale postseismic uplift after the 1999 Mw7.1 Hector Mine and 2010 Mw7.2 El Mayor-Cucapah earthquakes. We also identify areas of vertical land motions due to anthropogenic, natural, and magmatic processes. Finally, we demonstrate the utility of the kinematic datum by improving the accuracy of high-spatial-resolution 12-day repeat-cycle Sentinel-1 Interferometric Synthetic Aperture Radar displacement and velocity maps.