Use of InSAR data for measuring land subsidence induced by groundwater withdrawal and climate change in Ardabil Plain, Iran

The Ardabil plain, with an approximate area of 1097.2 km2 in northwestern Iran, has experienced land subsidence due to intensive groundwater withdrawal and long seasons of drought in recent years. Different techniques have been used to investigate and evaluate subsidence in this region including: Global Positioning Systems (GPS), Levelling, and Geotechnical methods. These methods are typically expensive, time-consuming, and identify only a small fraction of the areas prone to subsidence. This study employs an Interferometric Synthetic Aperture Radar (InSAR) technique to measure the long-term subsidence of the plain. An open-source SAR interferometry time series analysis package, LiCSBAS, that integrates with the automated Sentinel-1 InSAR processor (COMET-LiCSAR) is used to analyze Sentinel-1 satellite images from October 2014 to January 2021. Processing of Sentinel-1 images shows that the Ardabil plain has been facing rapid subsidence due to groundwater pumping and reduced rainfall, especially between May 2018 to January 2019. The maximum subsidence rate was 45 mm/yr, measured at the southeastern part of the plain. While providing significant advantages (less processing time and disk space) over other InSAR processing packages, implementation of the LiCSBAS processing package and its accuracy for land subsidence measurements at different scales needs further evaluation. This study provides a procedure for evaluating its efficiency and accuracy for land subsidence measurements by comparing its measurements with the results of the GMTSAR and geotechnical numerical modeling. The results of geotechnical numerical modeling showed land subsidence with an average annual rate of 38 mm between 2006 and 2020, which was close to measurements using the InSAR technique. Comparison of the subsidence measurements of the Ardabil plain using the LiCSBAS package with results obtained from other techniques shows that LiCSBAS is able to accurately detect land deformation at large scales (~ km). However, they may not be optimized for more local deformations such as infrastructure monitoring.