Land subsidence induced by groundwater level decline has spatiotemporal variations. Taking the Interferometric Synthetic Aperture Radar (InSAR) results and the groundwater subsidence data acquired by the monitoring stations as the source material, this paper aims to reveal the spatiotemporal heterogeneity of groundwater-land subsidence in Beijing plain by using the Wind Rose Map (WRM) method and the Change Point Analysis (CPA) method. The WRM results show that the amount and variation in subsidence differs in different directions. This method detected the formation of new subsidence centers and the slowdown of land subsidence in 2008. The CPA results show that obvious changes are detected in subsidence development at the Wangsiying (WSY), Tianzhu (TZ) and Wangjing (WJ) stations. However, there is a relatively stable trend of groundwater decline and land subsidence at the Tianzhu (TZ) station. The stages of land subsidence development show a significant response to groundwater. Moreover, changes in land subsidence also show delayed response behind the changes in groundwater level. The time-lag could be affected by the variation in amplitude of the groundwater level.