Underground mining has caused drastic disturbances to regional ecosystems and soil nutrients. Understanding the three-dimensional (3D) spatial distribution of soil nutrients in mining area farmland is crucial for agricultural production and environmental management. However, few studies have reported the 3D spatial distribution of soil organic matter (SOM) in coal mining subsidence area. In our study, a sequential Gaussian simulation (SGS) algorithm was used to analyse the spatial distribution of SOM based on observations of 180 soil samples in the Zhaogu mine in China. The results showed that the SOM content had considerable variation in spatial distribution at different soil depths (0–20, 20–40, 40–60 cm) and decreased with the increase in soil depth. The spatial variability of surface organic matter was the largest, and the coefficient of variation was 29.38%, which was moderately mutated. The spatial distribution of SOM also varied among slope locations. The SOM content was higher upslope and downslope than on the middle slope. In addition, given a threshold, SGS can be used to calculate the probability that the organic matter content at any position is lower or higher than the given value. The research results provide a reference for land reclamation and precision agriculture.

Soil organic matter (SOM) scale effects are critical for crop growth and food security, especially in coal–grain complexes. However, few studies describe the spatial variation in SOM and its influencing factors at different sampling scales. Here, geostatistical theory and mathematical statistical methods were adopted to analyze the spatial variation characteristics of and structural differences in SOM in the coal mining subsidence area at Zhaogu No. 2 Mine at different sampling scales. The results showed that SOM varied spatially at large, medium, and small scales, and the coefficients of variation were 28.07%, 14.93%, and 14.31%, respectively, which are moderate values. The characteristic functions of the SOM content at different sampling scales differed, and the spatial structure scale effect was obvious. The spatial distribution of the SOM content fitted by the multiscale fitting model method was generally the same as the spatial distribution law of the SOM content fitted by the single scale kriging interpolation method; however, in terms of the detailed expression and spatial distribution of small-scale SOM content, the fitting model method was more accurate, and the accuracy increased by 36%. At the different sampling scales, sample size and soil type had specific effects on the SOM spatial distribution. These results provide research concepts and technical countermeasures for improving food security and the ecological environment in the coal–grain complex and help ensure sustainable agricultural lands.

Urban impervious surface is considered one of main factors affecting urban heat island and urban waterlogging. It is commonly extracted utilizing the original linear spectral mixture analysis (LSMA) model. However, due to the deficiencies of this method, many improvements and modifications have been proposed. In this paper, a modified dynamic endmember linear spectral mixture analysis (DELSMA) model was introduced and tested in Zhengzhou, China, using different images of Landsat series satellites. The accuracy and performance of DELSMA model was evaluated in terms of $RMSE$ , $r$ and $R^2$ . Results show that (1) the DELSMA model performed equally well for Landsat-5 Thematic Mapper (TM) and Landsat-7 Enhanced Thematic Mapper (ETM+) images, and obtained better accuracy by using Landsat-8 Operational Land Imager (OLI) than Landsat TM/ETM+; (2) the DELSMA model achieved a better performance than the original LSMA model consistently, using images of Landsat from different sensors. Based exclusively on the overall accuracy reports, the DELSMA model proved to be a more efficient method for extracting impervious surface. Our study will provide a reliable method of impervious surface estimation for the urban planner and management in monitoring urban expansion, revealing urban heat island, and estimating urban surface runoff, using time-series Landsat imagery.

Soil nitrogen is very important for crop growth and development. However, the factors affecting the three-dimensional spatial distribution of soil total nitrogen (TN), particularly in coal mining subsidence areas, are unclear. In this study, Markov geostatistics was used to analyse the three-dimensional spatial distribution characteristics and influential factors of TN by examining 180 soil samples from the Zhaogu mine in China. The results showed that the TN content was significantly different at different soil depths (0–20, 20–40, 40–60 cm) and decreased with increasing soil depth. The variation coefficient of the TN content decreased gradually from top to bottom, ranging from 18.18 to 25.62%. In addition, the TN content was greatly affected by mining subsidence, rainfall, irrigation, fertilization and management mode. The factors that influenced the TN content also varied across different slope positions. The TN content of the upslope was the highest, and the TN content of the middle slope was the lowest. These results can provide research ideas and technical countermeasures for ecological environment improvement and sustainable land development in coal mining subsidence areas.