Groundwater resources, unlike surface water, are more vulnerable to disturbances and contaminations, as they take a very long time and significant cost to recover. So, predictive modeling and prevention strategies can empower policymakers for efficient groundwater governance through informed decisions and recommendations. Due to the importance of groundwater quality modeling, the hardness susceptibility mapping using machine learning (ML) models has not been explored. For the first time, the current research aimed to predict groundwater hardness susceptibility using the ML models. The performance of two ensemble models of boosted regression trees (BRT) and random forest (RF) is investigated through the arrangement of a comparative study with multivariate discriminant analysis (MDA). According to the hardness values in 135 groundwater quality monitoring wells, the hard and soft water are determined; then, 11 predictor variables including distance from the sea (DFS), land use, elevation, distance from the river (DFR), depth to groundwater (DTGW), pH, precipitation (PCP), evaporation (E), groundwater level (GWL), curvature, and lithology are used for predicting the groundwater hardness susceptibility map. Results indicated that the variables of DFR, DTGW, elevation, and DFS had a higher contribution to the modeling process. So, the high harness areas are mostly related to low elevations, low DTGW, and proximity to river and sea, which facilitate the percolation conditions for minerals containing calcium or magnesium into groundwater.

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.