Search Results (445)

Sea-level rise has important implications for the economic and infrastructure security of coastal cities. Land subsidence further exacerbates relative sea-level rise. The Beijing–Tianjin–Hebei region (BTHR) along the Bohai Bay is one of the areas most severely affected by ground subsidence in the world. This study applies the Small Baseline Subset Interferometric Synthetic Aperture Radar (SBAS InSAR) method to analyze 47 ALOS PALSAR-2 images with five frames, mapping subsidence across 21,677.7 km² and revealing spatial patterns and trends over time from 2015 to 2021. This is one of the few published research studies for large-scale and long-term analysis of its kind using ALOS-2 data in this region. The results reveal the existence of six major areas affected by severe subsidence in the study area, with the most pronounced in Jinzhan Town, Beijing, with the maximum subsiding velocity of −94.42 mm/y. Except for the two subsidence areas located in Chaoyang District of Beijing and Guangyang District of Langfang City, the other areas with serious subsidence detected are all located in suburban areas; this means that the strict regulations of controlling urban subsidence for downtown areas in the BTHR have worked. The accumulated subsidence is highly correlated with the time in the time series. Moreover, the subsidence of 161.4 km of the Beijing–Tianjin Inter-City High-Speed Railway (HSR) and 194.5 km of the Beijing–Shanghai HSR (out of a total length of 1318 km) were analyzed. It is the first time that PALSAR-2 data have been used to simultaneously investigate the subsidence along two important HSR lines in China and to analyze relatively long sections of the routes. The above two railways intersect five and seven subsiding areas, respectively. Within the range of the monitored railway line, the percentage of the section with subsidence velocity below −10 mm/y in the monitoring length range is 11.2% and 27.9%; this indicates that the Beijing–Shanghai HSR has suffered more serious subsidence than the Beijing–Tianjin Inter-City HSR within the monitoring period. This research is also beneficial for assessing the subsidence risk associated with different railways. In addition, this study further analyzed the potential reasons for the serious land subsidence of the identified areas. The results of the geological interpretation still indicate that the main cause of subsidence in the area is due to hydrogeological characteristics and underground water withdrawal. Full article

Remote sensing provides an efficient roadmap in geological analysis and interpretation. However, some challenges arise when remote sensing techniques are integrated with machine learning in geological surveys. Factors including irregular spatial distribution, sample imbalance, interclass resemblances, regolith, and geochemical similarities impede geological feature diagnosis, interpretation, and identification across varied remote sensing datasets. To address these limitations, a hybrid-attention-integrated long short-term memory (LSTM) network is employed to diagnose, interpret, and identify lithological feature representations in a remote sensing-based geological analysis using multisource data fusion. The experimental design integrates varied datasets including Sentinel-2A, Landsat-9, ASTER, ALOS PALSAR DEM, and Bouguer anomaly gravity data. The proposed model incorporates a hybrid attention mechanism (HAM) comprising channel and spatial attention submodules. HAM utilizes an adaptive technique that merges global-average-pooled features with max-pooled features, enhancing the model’s accuracy in identifying lithological units. Additionally, a channel separation operation is employed to allot refined channel features into clusters based on channel attention maps along the channel dimension. The comprehensive analysis of results from comparative extensive experiments demonstrates HAM-LSTM’s state-of-the-art performance, outperforming existing attention modules and attention-based models (ViT, SE-LSTM, and CBAM-LSTM). Comparing HAM-LSTM to baseline LSTM, the HAM module’s integrated configurations equip the proposed model to better diagnose and identify lithological units, thereby increasing the accuracy by 3.69%. Full article

The North Java coastal area, known as the Pantura region, is experiencing significant land subsidence, with certain areas sinking up to 10 cm per year. Pekalongan is among the most affected, with subsidence rates between 10 and 19 cm annually, mainly due to groundwater extraction, sediment compaction, and coastal erosion. Other coastal cities, like Semarang and Demak, show rates averaging 4 to 10 cm per year. This rapid subsidence is due to favorable geological conditions and ongoing urban development. This study investigates land subsidence in Pekalongan using the PS-InSAR method and dynamic visualization of time-series land cover data. PS-InSAR was applied to 45 scenes from ALOS-2 PALSAR-2 to monitor subsidence from 2014 to 2022. The results were validated with in situ subsidence benchmarks. Urban development dynamics were analyzed through land cover and land use change (LULC) and population density over the same period, using the GLC_FCS30D dataset in the GEE to detect non-natural LULC. The PS-InSAR results indicated that over 60.9% of investigation points experienced subsidence, up to 100 cm between 2014 and 2022. Ground validation showed an 83% agreement with PS-InSAR results. A statistical analysis of LULC from 2014 to 2022 did not show significant built-up area development, but the extension of salt marshes and water bodies indicated subsidence expansion. The population density reached 6873 people per square km by 2022, causing extensive groundwater use for domestic and industrial purposes, further aggravating the subsidence. Full article

Estimation of forest above-ground biomass (AGB) using multi-source remote sensing data is an important method to improve the accuracy of the estimate. However, selecting remote sensing factors that can effectively improve the accuracy of forest AGB estimation from a large amount of data is a challenge when the sample size is small. In this regard, the Least Absolute Shrinkage and Selection Operator (Lasso) has advantages for extensive redundant variables but still has some drawbacks. To address this, the study introduces two Least Absolute Shrinkage and Selection Operator Lasso-based variable selection methods: Least Absolute Shrinkage and Selection Operator Genetic Algorithm (Lasso-GA) and Variance Inflation Factor Least Absolute Shrinkage and Selection Operator (VIF-Lasso). Sentinel 2, Sentinel 1, Landsat 8 OLI, ALOS-2 PALSAR-2, Light Detection and Ranging, and Digital Elevation Model (DEM) data were used in this study. In order to explore the variable selection capabilities of Lasso-GA and VIF-Lasso for remote sensing estimation of forest AGB. It compares Lasso-GA and VIF-Lasso with Boruta, Random Forest Importance Selection, Pearson Correlation, and Lasso for selecting remote sensing factors. Additionally, it employs eight machine learning models—Random Forest (RF), Extreme Gradient Boosting (XGBoost), Support Vector Machine (SVM), Bayesian Regression Neural Network (BRNN), Elastic Net (EN), K-Nearest Neighbors (KNN), Extremely Randomized Trees (ETR), and Stochastic Gradient Boosting (SGBoost)—to estimate forest AGB in Wuyi Village, Zhenyuan County. The results showed that the optimized Lasso variable selection could improve the accuracy of forest biomass estimation. The VIF-Lasso method results in a BRNN model with an R² of 0.75 and an RMSE of 16.48 Mg/ha. The Lasso-GA method results in an ETR model with an R² of 0.73 and an RMSE of 16.70 Mg/ha. Compared to the optimal SGBoost model with the Lasso variable selection method (R² of 0.69, RMSE of 18.63 Mg/ha), the VIF-Lasso method improves R² by 0.06 and reduces RMSE by 2.15 Mg/ha, while the Lasso-GA method improves R² by 0.04 and reduces RMSE by 1.93 Mg/ha. From another perspective, they also demonstrated that the RX sample count and sensitivity provided by LiDAR, as well as the Horizontal Transmit, Vertical Receive provided by Microwave Radar, along with the feature variables (Mean, Contrast, and Correlation) calculated from the Green, Red, and NIR bands of optical remote sensing in 7 × 7 and 5 × 5 windows, play an important role in forest AGB estimation. Therefore, the optimized Lasso variable selection method shows strong potential for forest AGB estimation using multi-source remote sensing data. Full article

Grazing lands are crucial for agricultural productivity, ecological stability, and carbon sequestration, underscoring the importance of monitoring vegetation biomass for the effective management of these ecosystems. Remote sensing data, including optical vegetation indices (VIs) like the Normalized Difference Vegetation Index (NDVI), are widely used to monitor vegetation dynamics due to their simplicity and high sensitivity. In contrast, radar-based VIs, such as the Polarimetric Radar Vegetation Index (PRVI), offer additional advantages, including all-weather imaging capabilities, a wider saturation range, and sensitivity to the vegetation structure information. This study introduces an enhanced form of the PRVI, termed the Normalized PRVI (NPRVI), which is calibrated to a 0 to 1 range, constraining the minimum value to reduce the background effects. The calibration and range factor were derived from statistical analysis of PRVI components across vegetated regions in the Contiguous United States (CONUS), using dual-polarization C-band Sentinel-1 and L-band ALOS-PALSAR data on the Google Earth Engine (GEE) platform. Machine learning models using NPRVI and NDVI demonstrated their complementarity with annual herbaceous biomass data from the Rangeland Analysis Platform. The results showed that the Random Forest Model outperformed the other machine learning models tested, achieving R² ≈ 0.51 and MAE ≈ 498 kg/ha (relative MAE ≈ 32.1%). Integrating NPRVI with NDVI improved biomass estimation accuracy by approximately 10% compared to using NDVI alone, highlighting the added value of incorporating radar-based vegetation indices. NPRVI may enhance the monitoring of grazing lands with relatively low biomass compared to other vegetation types, while also demonstrating applicability across a broad range of biomass levels and in diverse vegetation covers. Full article

Flash floods are one of the most hazardous natural events globally, characterized by their rapid onset and unpredictability, often overwhelming emergency preparedness and response systems. In the arid environment of Saudi Arabia, Wadi Al-Lith watershed is particularly prone to flash floods, exacerbated by sudden storms and the region’s distinct topographical features. This study focuses on the morphometric characterization and comparative analysis of flash flood risk within the Wadi Al-Lith basin. To assess flood susceptibility, two widely adopted methodologies were employed: the morphometric ranking approach and El-Shamy’s method. A 12.5-m resolution ALOS PALSAR digital elevation model (DEM) was used to delineate the watershed and generate a detailed drainage network via Arc-Hydro tools in the ArcGIS 10.4 software. Fifteen morphometric parameters were analyzed to determine their influence on flood potential and hazard prioritization. The findings of this study provide crucial insights for regional flood risk management, offering an improved understanding of flash flood dynamics and assisting in developing effective mitigation strategies for Wadi Al-Lith and similar environments. The findings reveal that Wadi Al-Lith comprises multiple sub-catchments with varying degrees of vulnerability to flash flooding. According to the morphometric hazard analysis (MHA), certain sub-catchments, including sc-2, sc-4, sc-5, sc-6, sc-10, sc-12, sc-13, and sc-15, emerge as highly susceptible to flood hazards, while others (sc-1 and sc-9) fall into moderate risk categories. In contrast, the application of El-Shamy’s method provides a different ranking of flood risks across the watershed’s sub-catchments, offering a comparative view of flood susceptibility. The insights gained from this dual-analysis approach are expected to support the development of targeted flood prevention and mitigation strategies, which are essential for minimizing the future impacts of flash flooding in the Wadi Al-Lith watershed and ensuring better preparedness for local communities. Full article

Ground subsidence in urban areas is mainly due to natural or anthropogenic activities, and it seriously threatens the healthy and sustainable development of the city and the security of individuals’ lives and assets. Shanghai is a megacity of China, and it has a long history of ground subsidence due to the overexploitation of groundwater and urban expansion. Time Series Synthetic Aperture Radar Interferometry (TS-InSAR) is a highly effective and widely used approach for monitoring urban ground deformation. However, it is difficult to obtain long-term (such as over 10 years) deformation results using single-platform SAR satellite in general. To acquire long-term surface deformation monitoring results, it is necessary to integrate data from multi-platform SAR satellites. Furthermore, the deformations are the result of multiple factors that are superimposed, and relevant studies that quantitatively separate the contributions from different driving factors to subsidence are rare. Moreover, the time series cumulative deformation results of massive measurement points also bring difficulties to the deformation interpretation. In this study, we have proposed a long-term surface deformation monitoring and quantitative interpretation method that integrates multi-platform TS-InSAR, PCA, and K-means clustering. SAR images from three SAR datasets, i.e., 19 L-band ALOS-1 PALSAR, 22 C-band ENVISAT ASAR, and 20 C-band Sentinel-1A, were used to retrieve annual deformation rates and time series deformations in Shanghai from 2007 to 2018. The monitoring results indicate that there is serious uneven settlement in Shanghai, with a spatial pattern of stability in the northwest and settlement in the southeast of the study area. Then, we selected Pudong International Airport as the area of interest and quantitatively analyzed the driving factors of land subsidence in this area by using PCA results, combining groundwater exploitation and groundwater level change, precipitation, temperature, and engineering geological and human activities. Finally, the study area was divided into four sub-regions with similar time series deformation patterns using the K-means clustering. This study helps to understand the spatiotemporal evolution of surface deformation and its driving factors in Shanghai, and provides a scientific basis for the formulation and implementation of precise prevention and control strategies for land subsidence disasters, and it can also provide reference for monitoring in other urban areas. Full article

Satellite-based estimation of forest variables including forest biomass relies on model-based approaches since forest biomass cannot be directly measured from space. Such models require ground reference data to adapt to the local forest structure and acquired satellite data. For wide-area mapping, such reference data are too sparse to train the biomass retrieval model and approaches for calibrating that are independent from training data are sought. In this study, we compare the performance of one such calibration approach with the traditional regression modelling using reference measurements. The performance was evaluated at four sites representative of the major forest biomes in Europe focusing on growing stock volume (GSV) prediction from time series of C-band Sentinel-1 and Advanced Land Observing Satellite Phased Array L-band Synthetic Aperture Radar (ALOS-2 PALSAR-2) backscatter measurements. The retrieval model was based on a Water Cloud Model (WCM) and integrated two forest structural functions. The WCM trained with plot inventory GSV values or calibrated with the aid of auxiliary data products correctly reproduced the trend between SAR backscatter and GSV measurements across all sites. The WCM-predicted backscatter was within the range of measurements for a given GSV level with average model residuals being smaller than the range of the observations. The accuracy of the GSV estimated with the calibrated WCM was close to the accuracy obtained with the trained WCM. The difference in terms of root mean square error (RMSE) was less than 5% units. This study demonstrates that it is possible to predict biomass without providing reference measurements for model training provided that the modelling scheme is physically based and the calibration is well set and understood. Full article

The Greater Amanzule Peatlands (GAP) in Ghana is an important biodiversity hotspot facing increasing pressure from anthropogenic land-use activities driven by rapid agricultural plantation expansion, urbanisation, and the burgeoning oil and gas industry. Accurate measurement of how these pressures alter land cover over time, along with the projection of future changes, is crucial for sustainable management. This study aims to analyse these changes from 2010 to 2020 and predict future scenarios up to 2040 using multi-source remote sensing and machine learning techniques. Optical, radar, and topographical remote sensing data from Landsat-7, Landsat-8, ALOS/PALSAR, and Shuttle Radar Topography Mission derived digital elevation models (DEMs) were integrated to perform land cover change analysis using Random Forest (RF), while Cellular Automata Artificial Neural Networks (CA-ANNs) were employed for predictive modelling. The classification model achieved overall accuracies of 93% in 2010 and 94% in both 2015 and 2020, with weighted F1 scores of 80.0%, 75.8%, and 75.7%, respectively. Validation of the predictive model yielded a Kappa value of 0.70, with an overall accuracy rate of 80%, ensuring reliable spatial predictions of future land cover dynamics. Findings reveal a 12% expansion in peatland cover, equivalent to approximately 6570 ± 308.59 hectares, despite declines in specific peatland types. Concurrently, anthropogenic land uses have increased, evidenced by an 85% rise in rubber plantations (from 30,530 ± 110.96 hectares to 56,617 ± 220.90 hectares) and a 6% reduction in natural forest cover (5965 ± 353.72 hectares). Sparse vegetation, including smallholder farms, decreased by 35% from 45,064 ± 163.79 hectares to 29,424 ± 114.81 hectares. Projections for 2030 and 2040 indicate minimal changes based on current trends; however, they do not consider potential impacts from climate change, large-scale development projects, and demographic shifts, necessitating cautious interpretation. The results highlight areas of stability and vulnerability within the understudied GAP region, offering critical insights for developing targeted conservation strategies. Additionally, the methodological framework, which combines optical, radar, and topographical data with machine learning, provides a robust approach for accurate and detailed landscape-scale monitoring of tropical peatlands that is applicable to other regions facing similar environmental challenges. Full article

The Vilcanota is the second-largest snow-capped mountain range in Peru, featuring 380 individual glaciers, each with its own unique characteristics that must be studied independently. However, few studies have been conducted in the Vilcanota range to monitor and track the area and volume changes of the Suyuparina and Quisoquipina glaciers. Notably, there are only a few studies that have approached this issue using LIDAR technology. Our methodology is based on a combination of optical, radar and LIDAR data sources, which allowed for constructing coherent temporal series for the both the perimeter and volume changes of the Suyuparina and Quisoquipina glaciers while accounting for the uncertainty in the perimeter detection procedure. Our results indicated that, from 1990 to 2013, there was a reduction in snow cover of 12,694.35 m² per year for Quisoquipina and 16,599.2 m² per year for Suyuparina. This represents a loss of 12.18% for Quisoquipina and 22.45% for Suyuparina. From 2006 to 2013, the volume of the Quisoquipina glacier decreased from 11.73 km³ in 2006 to 11.04 km³ in 2010, while the Suyuparina glacier decreased from 6.26 km³ to 5.93 km³. Likewise, when analyzing the correlation between glacier area and precipitation, a moderate inverse correlation (R = −0.52, p < 0.05) was found for Quisoquipina. In contrast, the correlation for Suyuparina was low and nonsignificant, showing inconsistency in the effect of precipitation. Additionally, the correlation between the snow cover area and the annual mean air temperature (R = −0.34, p > 0.05) and annual minimum air temperature (R = −0.36, p > 0.05) was low, inverse, and not significant for Quisoquipina. Meanwhile, snow cover on Suyuparina had a low nonsignificant correlation (R = −0.31, p > 0.05) with the annual maximum air temperature, indicating a minimal influence of the measured climatic variables near this glacier on its retreat. In general, it was possible to establish a reduction in both the area and volume of the Suyuparina and Quisoquipina glaciers based on freely accessible remote sensing data. Full article

Mapping forest canopy height is critical for climate modeling and forest management, and tropical forests present unique challenges for remote sensing due to their dense vegetation and complex structure. The advent of ICESat-2 and GEDI, two advanced lidar datasets, offers new opportunities for improving canopy height estimation. In this study, we used footprint-level canopy height products from ICESat-2 and GEDI, combined with features extracted from Landsat-8, PALSAR-2, and FABDEM products. The AutoGluon stacking ensemble learning algorithm was employed to construct inversion models, generating 30 m resolution continuous canopy height maps for the tropical forests of Puerto Rico. Accuracy validation was performed using the high-resolution G-LiHT airborne lidar products. Results show that tropical forest canopy height inversion remains challenging, with all models yielding relative root mean square errors (rRMSE) exceeding 0.30. The stacking ensemble model outperformed all base learners, and the GEDI-based map had slightly higher accuracy than the ICESat-2-based map, with RMSE values of 4.81 and 4.99 m, respectively. Both models showed systematic biases, but the GEDI-based model exhibited less underestimation for taller canopies, making it more suitable for biomass estimation. The proposed approach can be applied to other forest ecosystems, enabling fine-resolution canopy height mapping and enhancing forest conservation efforts. Full article

This study investigates soil loss in the Wadi Bin Abdullah watershed using the Revised Universal Soil Loss Equation (RUSLE) combined with advanced tools, such as remote sensing and the Geographic Information System (GIS). By leveraging the ALOS PALSAR Digital Elevation Model (DEM), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) rainfall data, and the Digital Soil Map of the World (DSMW), the research accurately evaluates soil loss loads. The methodology identifies significant variations in soil loss rates across the entire watershed, with values ranging from 1 to 1189 tons per hectare per year. The classification of soil loss into four stages—very low (0–15 t/ha/yr), low (15–45 t/ha/yr), moderate (45–75 t/ha/yr), and high (>75 t/ha/yr)—provides a nuanced perspective on soil loss dynamics. Notably, 20% of the basin exhibited a soil loss rate of 36 tons per hectare per year. These high rates of soil erosion are attributed to certain factors, such as steep slopes, sparse vegetation cover, and intense rainfall events. These results align with regional and global studies and highlight the impact of topography, land use, and soil properties on soil loss. Moreover, the research emphasizes the importance of integrating empirical soil loss models with modern technological approaches to identify soil loss-prone locations and precisely quantify soil loss rates. These findings provide valuable insights for developing environmental management strategies aimed at mitigating the impacts of soil loss, promoting sustainable land use practices, and supporting resource conservation efforts in arid and semi-arid regions. Full article

The scientific field of precision agriculture employs increasingly innovative techniques to optimize inputs, maximize profitability, and reduce environmental impact. However, obtaining a high number of soil samples is challenging in order to make precision agriculture viable. There is a trade-off between the amount of data needed and the time and resources spent to obtain these data compared to the accuracy of the maps produced with more or fewer points. In the present study, the research was based on an exhaustive dataset of apparent electrical conductivity (aEC) containing 3906 points distributed along 26 transects with spacing between each of up to 40 m, measured by the proximal soil sensor EM38-MK2, for a grain-producing area of 72 ha in São Paulo, Brazil. A second sparse dataset was simulated, showing only four transects with a 400 m distance and, in the end, only 162 aEC points. The aEC map via ordinary kriging (OK) from the grid with 26 transects was considered the reference, and two other mapping approaches were used to map aEC via sparse grid: kriging with external drift (KED) and geographically weighted regression (GWR). These last two methods allow the increment of auxiliary variables, such as those obtained by remote sensors that present spatial resolution compatible with the pivot scale, such as data from the Landsat-8, Aster, and Sentinel-2 satellites, as well as ten terrain covariates derived from the Alos Palsar digital elevation model. The KED method, when used with the sparse dataset, showed a relatively good fit to the aEC data (R² = 0.78), with moderate prediction accuracy (MAE = 1.26, RMSE = 1.62) and reasonable predictability (RPD = 1.76), outperforming the GWR method, which had the weakest performance (R² = 0.57, MAE = 1.78, RMSE = 2.30, RPD = 0.81). The reference aEC map using the exhaustive dataset and OK showed the highest accuracy with an R² of 0.97, no systematic bias (ME = 0), and excellent precision (RMSE = 0.56, RPD = 5.86). Management zones (MZs) derived from these maps were validated using soil texture data from clay samples measured at 0–10 cm depth in a grid of 72 points. The KED method demonstrated the highest potential for accurately defining MZs for irrigation, producing a map that closely resembled the reference MZ map, thereby providing reliable guidance for irrigation management. Full article

The determination of swelling/shrinking phenomena, from natural and anthropogenic activity, is examined in this study through the synergy of various remote sensing methodologies. For the period of 2016–2022, a time-series InSAR analysis of Sentinel-1 satellite images, with a Coherent Change Detection procedure, was conducted to calculate the Normalized Coherence Difference. These were combined with Sentinel-2 multispectral data by exploiting the Normalized Difference Vegetation Index to create multi-temporal image composites. In addition, ALOS-Palsar DEM derivatives highlighted the geomorphological characteristics, which, in conjunction with the satellite imagery outcomes and other auxiliary spatial datasets, were embedded within a Multi-Criteria Decision Analysis (MCDA) model. The synergy of the remote sensing and GIS techniques’ applicability within the MCDA model highlighted the zones undergoing seasonal swelling/shrinking processes in Pyrgos–Parekklisia and Moni regions in Cyprus. The accuracy assessment of the produced final MCDA outcome provided an overall accuracy of 72.4%, with the Kappa statistic being 0.66, indicating substantial agreement of the MCDA outcome with the results from a Persistent Scatterer Interferometry analysis and ground-truth observations. Thus, this study offers decision-makers a powerful procedure to monitor longer- and shorter-term swelling/shrinking phenomena. Full article

The realization of a harmonious relationship between the natural environment and economic development has always been the unremitting pursuit of traditional mineral resource-based cities. With rich reserves of iron and coal ore resources, Laiwu has become an important steel production base in Shandong Province in China, after several decades of industrial development. However, some serious environmental problems have occurred with the quick development of local steel industries, with ground subsidence and consequent secondary disasters as the most representative ones. To better evaluate possible ground collapse risk, comprehensive approaches incorporating the common deformation monitoring with small-baseline subset (SBAS)-synthetic aperture radar interferometry (InSAR) technique, environmental factors analysis, and risk evaluation are designed here with ALOS PALSAR and Sentinel-1 SAR observations. A retrospect on the ground deformation process indicates that ground deformation has largely decreased by around 51.57% in area but increased on average by around −5.4 mm/year in magnitude over the observation period of Sentinel-1 (30 July 2015 to 22 August 2022), compared to that of ALOS PALSAR (17 January 2007 to 28 October 2010). To better reveal the potential triggering mechanism, environmental factors are also utilized and conjointly analyzed with the ground deformation time series. These analysis results indicate that the ground deformation signals are highly correlated with human industrial activities, such underground mining, and the operation of manual infrastructures (landfill, tailing pond, and so on). In addition, the evaluation demonstrates that the area with potential collapse risk (levels of medium, high, and extremely high) occupies around 8.19 km², approximately 0.86% of the whole study region. This study sheds a bright light on the safety guarantee for the industrial operation and the ecologically friendly urban development of traditional steel production industrial cities in China. Full article