As a result of rapid societal development and urbanization, the pumping of groundwater has gradually increased. Land subsidence has thus become a common geological disaster, which can result in huge economic losses. Interferometric synthetic aperture radar (InSAR), with its large-scale and high-accuracy monitoring characteristics, can attain information on Earth surface deformation using the interferometric phase between couples of SAR images acquired at different times. Time-series results for the ground surface are the key information required to understand the deformation pattern and further study the reason for the subsidence. However, in recent research, most methods for resolving time-series deformation—like the Berardino method—that use residuals in functional model solving and distinguish high-pass displacement and the atmospheric component by filtering do not generally work well and functional models focusing on prior information in the time-series solution process are not always available. In this paper, to solve the above problems, 34 Sentinel-1A descending mode scenes of Mexico City captured between 2015/04/13 and 2016/09/10 are used as experimental data. Firstly, a new functional model is provided to obtain the deformation time-series. The nonlinear deformation and atmospheric phase are combined as an unknown parameter and the method of singular value decomposition (SVD) is used to solve this variable. The nonlinear displacement and atmospheric phase are then separated by the singular spectrum analysis (SSA) method. Finally, the total land subsidence time-series is obtained by adding together the linear displacement and nonlinear displacement. Two typical methods and the proposed method were compared using both unit weights and adaptive weights. The experimental results show that the proposed method can obtain a more accurate time-series deformation result. Moreover, the different weights do not result in significant differences and the solved atmospheric and nonlinear phases have good consistency with the interferogram phase.
Time-series analysis of surface deformation at Brady Hot Springs geothermal field (Nevada) using interferometric synthetic aperture radar
