Scrutinization of land subsidence rate using a supportive predictive model: Incorporating radar interferometry and ensemble soft-computing.

Land subsidence is a huge challenge that land and water resource managers are still facing. Radar datasets revolutionize the way and give us the ability to provide information about it, thanks to their low cost. But identifying the most important drivers need for the modeling process. Machine learning methods are especially top of mind amid the prediction studies of natural hazards and hit new heights over the last couple of years. Hence, putting an efficient approach like integrated radar-and-ensemble-based method into practice for land subsidence rate simulation is not available yet which is the main aim of this research. In this study, the number of 52 pairs of radar images were used to identify subsidence from 2014 to 2019. Then, using the simulated annealing (SA) algorithm the key variables affecting land subsidence were identified among the topographical parameters, aquifer information, land use, hydroclimatic variables, and geological and soil factors. Afterward, three individual machine learning models (including Support Vector Machine, SVM; Gaussian Process, GP; Bayesian Additive Regression Tree, BART) along with three ensemble learning approaches were considered for land subsidence rate modeling. The results indicated that the subsidence varies between 0 and 59 cm in this period. Comparing the Radar results with the permanent geodynamic station exhibited a very strong correlation between the ground station and the radar images (R = 0.99, RMSE = 0.008). Parsing the input data by the SA indicated that key drivers are precipitation, elevation, percentage of fine-grained materials in the saturated zone, groundwater withdrawal, distance to road, groundwater decline, and aquifer thickness. The performance comparison indicated that ensemble models perform better than individual models, and among ensemble models, the nonlinear ensemble approach (i.e., BART model combination) provided better performance (RMSE = 0.061, RSR = 0.42, R = 0.83, PBIAS = 2.2). Also, the distribution shape of the probability density function in the non-linear ensemble model is much closer to the observations. Results indicated that the presence of significant fine-grained materials in unconsolidated aquifer systems can clarify the response of the aquifer system to groundwater decline, low recharge, and subsequent land subsidence. Therefore, the interaction between these factors can be very dangerous and intensify subsidence.