Landslide susceptibility mapping is the first and most important step involved in landslide hazard assessment. The purpose of the present study is to compare three nonlinear approaches for landslide susceptibility mapping and test whether coal mining has a significant impact on landslide occurrence in coal mine areas. Landslide data collected by the Bureau of Land and Resources are represented by the X, Y coordinates of its central point; causative factors were calculated from topographic and geologic maps, as well as satellite imagery. The five-fold cross-validation method was adopted and the landslide/non-landslide datasets were randomly split into a ratio of 80:20. From this, five subsets for 20 times were acquired for training and validating models by GIS Geostatistical analysis methods, and all of the subsets were employed in a spatially balanced sample design. Three landslide models were built using support vector machine (SVM), logistic regression (LR), and artificial neural network (ANN) models by selecting the median of the performance measures. Then, the three fitted models were compared using the area under the receiver operating characteristics (ROC) curves (AUC) and the performance measures. The results show that the prediction accuracies are between 73.43% and 87.45% in the training stage, and 67.16% to 73.13% in the validating stage for the three models. AUCs vary from 0.807 to 0.906 and 0.753 to 0.944 in the two stages, respectively. Additionally, three landslide susceptibility maps were obtained by classifying the range of landslide probabilities into four classes representing low (0–0.02), medium (0.02–0.1), high (0.1–0.85), and very high (0.85–1) probabilities of landslides. For the distributions of landslide and area percentages under different susceptibility standards, the SVM model has more relative balance in the four classes compared to the LR and the ANN models. The result reveals that the SVM model possesses better prediction efficiency than the other two models. Furthermore, the five factors, including lithology, distance from the road, slope angle, elevation, and land-use types, are the most suitable conditioning factors for landslide susceptibility mapping in the study area. The mining disturbance factor has little contribution to all models, because the mining method in this area is underground mining, so the mining depth is too deep to affect the stability of the slopes.

It is of great significance for the monitoring and protection of the original ecological environment in coal mining areas to identify the ground subsidence and quantify its influence on the surface vegetation. The surface deformation and vegetation information were obtained by using spaceborne SAR and Landsat OLI images in the Xishan Coalfield. The relative change rate, coefficient of variation, and trend analysis methods were used to compare the vegetation growth trends in the subsidence center, subsidence edge, and non-subsidence zones; and the vegetation coverage was predicted by the pixel dichotomy and grey model from 2021 to 2025. The results indicated that the proportions of vegetation with high fluctuation and serious degradation were 6.60% and 5.64% in the subsidence center, and its NDVI values were about 10% lower than that in the subsidence edge and non-subsidence zones. In addition, vegetation coverage showed a wedge ascending trend from 2013 to 2020, and the prediction values of vegetation coverage obtained by GM (1,1) model also revealed this trend. The residuals of the predicted values were 0.047, 0.047, and 0.019 compared with the vegetation coverage in 2021, and the vegetation coverage was the lowest in the subsidence center, which was consistent with the law obtained by using NDVI. Research suggested that ground subsidence caused by mining activities had a certain impact on the surface vegetation in the mining areas; the closer to the subsidence center, the greater the fluctuation of NDVI, and the stronger the vegetation degradation trend; conversely, the smaller the fluctuation, and the more stable the vegetation growth.

The Qilian Mountains, located on the northeastern edge of the Qinghai–Tibet Plateau, are characterized by unique high-altitude and cold-climate terrain, where permafrost and seasonally frozen ground are extensively distributed. In recent years, with global warming and increasing precipitation on the Qinghai–Tibet Plateau, permafrost degradation has become severe, further exacerbating the fragility of the ecological environment. Therefore, timely research on surface deformation and the freeze–thaw patterns of alpine permafrost in the Qilian Mountains is imperative. This study employs Sentinel-1A SAR data and the SBAS-InSAR technique to monitor surface deformation in the alpine permafrost regions of the Qilian Mountains from 2017 to 2023. A method for spatiotemporal interpolation of ascending and descending orbit results is proposed to calculate two-dimensional surface deformation fields further. Moreover, by constructing a dynamic periodic deformation model, the study more accurately summarizes the regular changes in permafrost freeze–thaw and the trends in seasonal deformation amplitudes. The results indicate that the surface deformation time series in both vertical and east–west directions obtained using this method show significant improvements in accuracy over the initial data, allowing for a more precise reflection of the dynamic processes of surface deformation in the study area. Subsidence is predominant in permafrost areas, while uplift mainly occurs in seasonally frozen ground areas near lakes and streams. The average vertical deformation rate is 1.56 mm/a, with seasonal amplitudes reaching 35 mm. Topographical (elevation; slope gradient; aspect) and climatic factors (temperature; soil moisture; precipitation) play key roles in deformation patterns. The deformation of permafrost follows five distinct phases: summer thawing; warm-season stability; frost heave; winter cooling; and spring thawing. This study enhances our understanding of permafrost deformation characteristics in high-latitude and high-altitude regions, providing a reference for preventing geological disasters in the Qinghai–Tibet Plateau area and offering theoretical guidance for regional ecological environmental protection and infrastructure safety.

Based on the Chinese Geomorphologic Database at 1 : 1 000 000 scales, the distribution of the typical loess geomorphologic types (such as the loess tableland, loess ridge and loess knoll) is acquired in the Loess Plateau of China. Then, based on the SRTM (Shuttle Radar Topography Mission) digital elevation model (DEM) data and topographic analysis methods, the slope spectrums are computed for the typical loess geomorphologic types and their subtypes. Through achieving the tendency line of the slope spectrum and analysing the slope spectrum characteristics of the loess typical geomorphologic types, the transitional relationships are explored: (1) the general rule is that loess tableland transitions to loess ridge, and then to loess knoll. (2) The specific relationships for the subtypes are as follows: in loess tableland, the transition is from loess terrace to complete tableland, then to residual tableland, and finally to beam tableland. In the loess ridge, the transition is from oblique ridge to knoll ridge, and the final stage is the loess knoll.

Based on the results of remote sensing data interpretation, this paper aims to simulate and predict the mountain permafrost distribution changes affected by the mean decadal air temperature (MDAT), from the 1990s to the 2040s, in the Qilian Mountains. A bench-mark map is visually interpreted to acquire a mountain permafrost distribution from the 1990s, based on remote sensing images. Through comparison and estimation, a logistical regression model (LRM) is constructed using the bench-mark map, topographic and land coverage factors and MDAT data from the 1990s. MDAT data from the 2010s to the 2040s are predicted according to survey data from meteorological stations. Using the LRM, MDAT data and the factors, the probabilities (p) of decadal mountain permafrost distribution from the 1990s to the 2040s are simulated and predicted. According to the p value, the permafrost distribution statuses are classified as ‘permafrost probable’ (p > 0.7), ‘permafrost possible’ (0.7 ≥ p ≥ 0.3) and ‘permafrost improbable’ (p < 0.3). From the 1990s to the 2040s, the ‘permafrost probable’ type mainly degrades to that of ‘permafrost possible’, with the total area degenerating from 73.5 × 10³ km² to 66.5 × 10³ km². The ‘permafrost possible’ type mainly degrades to that of ‘permafrost impossible’, with a degradation area of 6.5 × 10³ km², which accounts for 21.3% of the total area. Meanwhile, the accuracy of the simulation results can reach about 90%, which was determined by the validation of the simulation results for the 1990s, 2000s and 2010s based on remote sensing data interpretation results. This research provides a way of understanding the mountain permafrost distribution changes affected by the rising air temperature rising over a long time, and can be used in studies of other mountains with similar topographic and climatic conditions.