## 3 Data, Methods and Results  
### SAR Data Processing  
The processing of InSAR time series analysis was completed using the LiCSBAS, an open-source package that integrates with LiCSAR products [53; 54; 55]. The processing chain and schematic description of data/methods used are shown in the flow chart of Figure S2. LiCSAR produced interferograms (wrapped and unwrapped) from SLC (single look complex) data of Sentinel-1 acquisitions. The interferograms were multilooked with a factor of 20 \(\times\) 4 in range and azimuth (46 \(\times\) 56 m spacing) and spatially filtered by a GAMMA adaptive power spectrum filter with an alpha value of 1.0 to decrease noise [56; 57]. With the SNAPHU software, LiCSAR unwraps the phase in two dimensions using a statistical cost technique [58]. We employed unwrapped interferograms and coherence images that were geocoded with a pixel spacing of 0.001 degree (100 m) and converted to GeoTIFF format for our time series analysis. We used the frames with ID: 080D_05196_131104 for descending track and 175A_05184_121313 for ascending track, respectively. For ascending and descending track (see Figure S1), we used 346 and 317 available interferograms from the LiCSAR portal, respectively (out of 99 and 120 SAR images, respectively; see Table S1 for image dates). The areas of interest are in the center-east of the frames ID 175A and 080D (Figure S1). The time span for the ascending track ranges was from 2015 to 2021 and for the descending from 2016 to 2021, respectively (see Table S1). Before the main processing, a tropospheric correction was applied on the unwrapped products using the generic atmospheric correction online service for InSAR (GACOS) data ([59]; pixel size 10 \(\times\) 10 km). GACOS utilizes the iterative tropospheric decomposition (ITD) model to separate stratified and turbulent signals from tropospheric total delays. The tropospheric corrections were used in GeoTIFF format. To evaluate the success of the correction, the calculation of standard deviation (STD) of unwrapped phase before and after the GACOS correction and the reduction rate of the STD was performed through the LiCSBAS software [53]. The STD of unwrapped phases for each interferogram was typically lower (see correlation diagrams in Figures 2 and 3; see also reduction of STD through the interferograms in Figures S3 and S4), indicating that this tropospheric correction approach greatly reduced tropospheric noise. We see in the diagrams (Figures 2 and 3) that the reduction of STD for the dataset of ascending orbit is lower (reduction rate 29.57%) than the reduction of STD for the descending orbit (36.29%). Figure 2: Correlation diagram showing standard deviation: (**a**), of unwrapped phases before and after the GACOS corrections and their reduction rates; (**b**), for ascending track 175. The gray line denotes a 1:1. The reduction of STD for the ascending orbit is moderate. The STD decreased from 6.35 rad to 4.41 rad on average with a mean reduction rate 29.57%. After the tropospheric correction, we performed the quality check function of LiCSBAS. Statistics such as average phase coherence (we applied a 0.3 threshold) and percentage of genuine unwrapped pixels are used to identify bad data. We also performed the network refinement by interferometric loop closure [60] and removed 100 interferograms (95 ascending and 5 descending orbit) from further processing. Unwrapped data may contain unwrapping mistakes, which can result in serious errors in the derived time series and should be eliminated or addressed prior to use. Moreover, after the quality checks, we applied the small baseline (SB) inversion on the network of interferograms (Figure 4). The least-squares approach was then used to get the mean displacement velocity (LOS) from the cumulative displacements. We used the NSBAS approach [61] to produce a more realistic time series of displacement even with a disconnected network (Figure 4).