Land Subsidence Susceptibility Mapping Using Interferometric Synthetic Aperture Radar (InSAR) and Machine Learning Models in a Semiarid Region of Iran

Abstract

Most published studies identify groundwater extraction as the leading cause of land subsidence (LS). However, the causes of LS are not only attributable to groundwater extraction. Other land-use practices can also affect the occurrence of LS. In this study, radar interferometric techniques and machine learning (ML) models were used for the prediction, susceptibility zoning, and prioritization of influential variables in the occurrence of LS in the Bakhtegan basin. The LS rate was characterized by applying an interferometric synthetic aperture radar (InSAR). The recursive feature elimination (RFE) method was used to detect and select the dominant combination of indicators to prepare an LS susceptibility map. Three ML models, including random forest (RF), k-nearest neighbors (KNN), and classification and regression trees (CART), were used to develop predictive models. All three models had acceptable performance. Among the ML models, the RF model performed the best (i.e., Nash–Sutcliffe efficiency, Kling–Gupta efficiency, correlation coefficient, and percent bias metrics of 0.76, 0.78, 0.88, and 0.70 for validating phase, respectively). The analysis conducted on all three ML model outputs showed that high and very high LS susceptibility classes were located on or near irrigated agricultural land. The results indicate that the leading cause of land LS in the study region is not due to groundwater withdrawals. Instead, the distance from dams and the proximity to anticlines, faults, and mines are the most important identifiers of LS susceptibility. Additionally, the highest probability of LS susceptibility was found at distances less than 18 km from synclines, 6 to 13 km from anticlines, 23 km from dams, and distances less than 20 to more than 144 km from mines. The validated methods presented in this study are reproducible, transferrable, and recommended for mapping LS susceptibility in semiarid and arid climate zones with similar environmental conditions.