Description: Satellite Interferometric Synthetic Aperture Radar (InSAR) is widely used for topographic, geological and natural resource investigations. However, most of the existing InSAR studies of ground deformation are based on relatively short periods and single sensors. This paper introduces a new multi-sensor InSAR time series data fusion method for time-overlapping and time-interval datasets, to address cases when partial overlaps and/or temporal gaps exist. A new Power Exponential Knothe Model (PEKM) fits and fuses overlaps in the deformation curves, while a Long Short-Term Memory (LSTM) neural network predicts and fuses any temporal gaps in the series. Taking the city of Wuhan (China) as experiment area, COSMO-SkyMed (2011–2015), TerraSAR-X (2015–2019) and Sentinel-1 (2019–2021) SAR datasets were fused to map long-term surface deformation over the last decade. An independent 2011–2020 InSAR time series analysis based on 230 COSMO-SkyMed scenes was also used as reference for comparison. The correlation coefficient between the results of the fusion algorithm and the reference data is 0.87 in the time overlapping region and 0.97 in the time-interval dataset. The correlation coefficient of the overall results is 0.78, which fully demonstrates that the algorithm proposed in our paper achieves a similar trend as the reference deformation curve. The experimental results are consistent with existing studies of surface deformation at Wuhan, demonstrating the accuracy of the proposed new fusion method to provide robust time series for the analysis of long-term land subsidence mechanisms.