Objectives In recent years, the climate change in Xinjiang region,China is complex. It is of great significance to study the terrestrial water storage (TWS) and its load deformation for the accurate maintenance of the regional reference frame.

Methods The satellite gravity data GRACE（gravity recovery and climate experiment）was used to invert the changes of land water storage in Xinjiang region from 2010 to 2014, and the independent component analysis was used to decompose it into spatiotemporal models, and the spatial characteristics of multiple time scales were analyzed. On this basis, the correction on the first-order of gravity potential coefficient and load Love number is carried out, the correction of the non-tidal atmospheric and ocean model (GAC) from GRACE is introduced to obtain the vertical change of the crust caused by the migration of the land water load. The GRACE inversion results are corrected by scale factor method derived from WaterGAP global hydrology model (WGHM). The influence of deformation is analyzed by introducing monthly precipitation data from global precipitation climatology project. And the deformation displacements of the 13 continuously operating reference system stations in the survey area are quantitatively compared. The correlation between TWS load signal of each station and the annual vertical displacement of global navigation satellite system (GNSS) is analyzed.

Results The results show that TWS in Xinjiang region presents multi-time scale characteristics, with obvious annual variation and long-term trend. The annual change is particularly significant near the western Pamirs Plateau, and the trend term is reflected in the long-term decrease year by year. The signal is relatively strong in the western Urumqi and Tianshan Mountains. In general, the change of land water load causes the vertical displacement of the crust in the study area, and the time series shows a small fluctuation with an amplitude of -1.5-1.5 mm, which is consistent with the deformation results of WGHM model. In addition to the obvious time lag, the time series of TWS load deformation and precipitation change are consistent with the overall trend. After GAC correction, the consistency between GRACE results and GNSS vertical displacement was enhanced, and the correlation coefficient and weighted root mean square（WRMS）were significantly improved. The correlation coefficient of each station was 0.63-0.91, and WRMS was 16.18%-58.97%.

Conclusions The research method in this paper can reveal more detailed hydrological signal characteristics in Xinjiang region, and the research results can provide important technical support for accurate maintenance of regional reference frame.