Spatiotemporal Heterogeneity of GNSS Vertical Displacements Driven by Environmental Loading Across the Complex Topography of Southwest China

Environmental loading is a major driver of nonlinear GNSS vertical displacements, yet its spatiotemporal heterogeneity remains insufficiently understood in regions with complex topography. In this study, we investigate the environmental loading effects on GNSS vertical motions across Southwest China using observations from a network of 66 stations. Singular Spectrum Analysis (SSA) and Empirical Orthogonal Function (EOF) analysis were applied to extract annual signals, while component-wise RMS reduction quantified hydrological and atmospheric loading contributions. Spatial statistical analysis, cross-wavelet transform, and k-means clustering examined correlation patterns and phase hysteresis between GNSS observations and modeled loads. Results show that hydrological loading dominates seasonal vertical oscillations, but crustal responses exhibit pronounced spatial heterogeneity controlled by regional topography and hydro-climatic gradients. EOF analysis reveals a dipole pattern induced by the Hengduan Mountains’moisture-blocking effect. Atmospheric loading anomalously dominates the eastern Sichuan Basin, whereas Yunnan displays strong amplitudes with high heterogeneity due to karst hydrogeology. Phase analysis identifies three distinct regimes: a rapid elastic response on the Tibetan Plateau, (with the lag of ~20 ± 5 days, correlation coefficient R ≈ 0.65), intermediate delays in Yunnan (~60 ± 5 days, R ≈ 0.58), and pronounced hysteresis in the Sichuan Basin (~105 ± 5 days, R ≈ 0.38) linked to slow groundwater diffusion and poroelastic processes. These findings highlight the critical role of local hydrogeological dynamics in modulating GNSS vertical deformation and provide new insights for improving environmental loading corrections in complex mountainous regions.