Abstract: Land subsidence is a natural geological phenomenon in which the surface elevation drops. It can severely destroy urban infrastructure and threaten urban safety if it occurs in densely populated cities with a high social development degree. The analysis of the evolution characteristics of land subsidence can reflect the degree of the influence of land subsidence on the ground infrastructures, and building an efficient land subsidence prediction model is of great significance for preventing and controlling land subsidence and protecting urban safety. This study obtained the spatial-temporal information on land subsidence using the persistent scatterer interferometric synthetic aperture Radar (PS-InSAR) method first and then verified the information using leveling to get high precision. Then, this study analyzed the general spatial-temporal characteristics of the land subsidence field using an empirical orthogonal function. The analysis results are as follows. Spatial modal No. 1 had a high variance contribution rate, almost representing the general spatial evolution of the study area. Its corresponding time coefficient showed a significant linear trend. By contrast, spatial mode No. 2 had a low variance contribution rate and a seasonally significant time coefficient. Finally, the time series of the regional land subsidence were predicted using both long short-term memory (LSTM) and Attention-LSTM models. The prediction results indicate that the Attention-LSTM model was superior to the LSTM model, with the mean square error loss (MSE-loss) of as low as 0.01. This prediction method expands the application of deep learning in the study of land subsidence.