Assessment and Effect of Mining Subsidence on Farmland in Coal–Crop Overlapped Areas: A Case of Shandong Province, China

: Farmland protection and food security is highly focused on in China. However, coal mining has caused negative consequences to cropland in coal–crop overlapped areas (COA), especially in eastern China. Thus, revealing the spatiotemporal impact of coal mining on farmland on a large scale is crucial for coordinating coal mining and grainland protection. In this study, Shandong Province, a representative coal–grain composite area, was selected as a research case to evaluate the damage of mining subsidence on farmland. Firstly, the ﬁeld investigation and mining subsidence prediction revealed the current situation and trend of farmland damage caused by mining in 2021 and 2030. Then, we evaluated the impact of cropland damage on grain yield. Finally, farmland landscape patterns, ecological environment, and social stability due to mining subsidence were discussed. The results show that: (1) the damaged cropland in 2021 was 6.40 × 10 4 ha, of which 1.40 × 10 4 ha was non-yield. By 2030, the accumulative area of damaged cropland is estimated to reach 7.52 × 10 4 ha. (2) By 2025 and 2030, the farmland yield will be reduced by 16.44 × 10 4 t and 18.45 × 10 4 t in overlapped areas of Shandong. (3) The subsidence led to cropland fragmentation, and the terrestrial ecosystem became an aquatic ecosystem, further intensifying the contradiction between more people and less land. This study provides a reference for coordinating coal and grain production and formulating cropland protection strategies in similar regions. Meanwhile, it also provides a scientiﬁc basis for the government to formulate land reclamation indicators, technology, management, and acceptance standards and establish and implement the reclamation reward and punishment system.


Introduction
Coal is China's main energy, accounting for more than 55% of the total primary energy consumption, and it plays a strategic role in long-term development [1][2][3]. While coal mining has brought benefits, it has also led to land and ecological environmental damages [4]. For example, underground mining leads to numerous issues, such as ground subsidence [5,6], damage to land and vegetation [7,8], and the destruction of landscapes [9][10][11]. As 96% of China's coal mining is underground mining, this leads to a significant amount of land subsidence damage [12][13][14]. According to statistics, the accumulated area of coal mining subsidence in China has reached 1.35 × 10 6 ha, and it is still rising at an annual rate of 7.0 × 10 4 ha [15]. Most subsided land is distributed in plain areas with high groundwater levels in eastern China [16]. After ground subsidence, high-level groundwater quickly rises above the ground elevation, resulting in a large amount of surface water [17,18].
Coal subsidence areas with high groundwater levels in China are also primary coal and grain production areas, i.e., areas where the coal-crop overlapping area (COA) [19]. territory has a high phreatic level, abundant water sources, concentrated precipitation, fertile soil, and China's highest arable land rate. It is an important agricultural area for grain, cotton, and oil products. The province's land area is 15.80 × 10 4 km 2 , accounting for 1.64% of the country's total area, and the cultivated land area is 7.52 × 10 4 km 2 , accounting for 6.17% of the country's total cropland area. The agricultural conditions in the territory are excellent, and crops are grown in both summer and autumn. Wheat, corn, and sweet potatoes are the three major local food crops. for 1.64% of the country's total area, and the cultivated land area is 7.52 × 10 4 km 2 , accounting for 6.17% of the country's total cropland area. The agricultural conditions in the territory are excellent, and crops are grown in both summer and autumn. Wheat, corn, and sweet potatoes are the three major local food crops.
In addition, Shandong Province has rich mineral reserves, including 15 pillar minerals that the national economy depends on for development. It is one of China's 14 largescale coal production bases, the Luxi coal base. Its coal-bearing area is about 1.65 × 10 4 km 2 , accounting for 10.45% of the province's total land area. By the end of 2020, the province had 102 production mines and 3 mines under construction, mainly scattered around Jining, Tai'an, Zaozhuang, Jinan, Heze, Yantai, and Dezhou City (Figure 1). Shandong province is one of the earliest provinces in China to discover, exploit, and utilize coal resources. Due to long-term mining, many subsidence areas have been formed on the ground. In particular, Shandong Province is a high-water-level area. Large-scale coal mining has resulted in severe land damage, and the situation of cultivated land protection is even more serious.  In addition, Shandong Province has rich mineral reserves, including 15 pillar minerals that the national economy depends on for development. It is one of China's 14 large-scale coal production bases, the Luxi coal base. Its coal-bearing area is about 1.65 × 10 4 km 2 , accounting for 10.45% of the province's total land area. By the end of 2020, the province had 102 production mines and 3 mines under construction, mainly scattered around Jining, Tai'an, Zaozhuang, Jinan, Heze, Yantai, and Dezhou City ( Figure 1). Shandong province is one of the earliest provinces in China to discover, exploit, and utilize coal resources. Due to long-term mining, many subsidence areas have been formed on the ground. In particular, Shandong Province is a high-water-level area. Large-scale coal mining has resulted in severe land damage, and the situation of cultivated land protection is even more serious.

Distribution Characteristics of Coal-Crop Overlapped Areas
Using ArcGIS spatial analysis to superimpose the distribution map of coal resources and cropland resources, it was found that the overlapped areas of crop and coal production in Shandong Province are mainly around Heze City, Zaozhuang City, Jining City, Tai'an City, Liaocheng City, Dezhou City, Jinan City, Zibo City, Weifang City, and the northern part of Yantai City (Figure 2). The overlapped areas cover an area of 101.66 × 10 4 ha, which account for 6.44%, 13.53%, and 64.61% of the province's total land area, total cropland area, and coal-bearing area, respectively. The land use type in coal-bearing areas within Shandong Province was dominated by cultivated land. The established mining area in the province was 23.61 × 10 4 ha, accounting for 23.22% of its area of overlapped crop and coal production.
Using ArcGIS spatial analysis to superimpose the distribution map of coal resources and cropland resources, it was found that the overlapped areas of crop and coal production in Shandong Province are mainly around Heze City, Zaozhuang City, Jining City, Tai'an City, Liaocheng City, Dezhou City, Jinan City, Zibo City, Weifang City, and the northern part of Yantai City (Figure 2). The overlapped areas cover an area of 101.66 × 10 4 ha, which account for 6.44%, 13.53%, and 64.61% of the province's total land area, total cropland area, and coal-bearing area, respectively. The land use type in coal-bearing areas within Shandong Province was dominated by cultivated land. The established mining area in the province was 23.61 × 10 4 ha, accounting for 23.22% of its area of overlapped crop and coal production.

Current Situation of Damaged Farmland
Based on the subsidence area map provided by various mining areas in Shandong Province, this study developed a summary table for collecting various cities' current coal mining subsidence areas. Simultaneously, the current state of cultivated land destruction was analyzed in detail based on the provincial government and several cities' statistical planning data, in conjunction with several field studies and verifications.
The current area of coal mining subsidence, reduced-yield, and non-yield in the summary table, is provided by Shandong Provincial Energy Bureau, Natural Resources Bureau, and mining enterprises.

Mining Subsidence Prediction
The Mining Subsidence Prediction System (MSPS) can automatically calculate rock formations and surface deformation caused by coal mining based on mining conditions, prediction functions, and parameters. The mathematical models used by the software include the probability integral method, mountainous surface movement prediction model, and steep prediction model for coal seam mining. Among them, the probability integral

Current Situation of Damaged Farmland
Based on the subsidence area map provided by various mining areas in Shandong Province, this study developed a summary table for collecting various cities' current coal mining subsidence areas. Simultaneously, the current state of cultivated land destruction was analyzed in detail based on the provincial government and several cities' statistical planning data, in conjunction with several field studies and verifications.
The current area of coal mining subsidence, reduced-yield, and non-yield in the summary table, is provided by Shandong Provincial Energy Bureau, Natural Resources Bureau, and mining enterprises.

Mining Subsidence Prediction
The Mining Subsidence Prediction System (MSPS) can automatically calculate rock formations and surface deformation caused by coal mining based on mining conditions, prediction functions, and parameters. The mathematical models used by the software include the probability integral method, mountainous surface movement prediction model, and steep prediction model for coal seam mining. Among them, the probability integral method is commonly used for predicting mining subsidence in coal mining areas and is suitable for coal seams with an inclination angle of less than 45 • [42]. The coal seams in all mining areas in Shandong Province are all near-level coal seams, thus conforming to the scope of application of the probability integral method. When all the coal seams in the mining area are fully mined, the amount of subsidence at any point A (x, y) on the surface is W A [43]. The calculation formula is shown in Equation (1). where W 0 represents the maximum sink value, and it satisfies the following Equation (2); r represents the main influence radius, and variable r satisfies the following Equation (3).
where m denotes the mining thickness of the coal seam; q refers to the subsidence coefficient; α is the coal seam inclination.
where H stands for the elevation difference between the point A (x, y) to be calculated on the ground and the micro-element point (c, t) on the coal seam; tan β represents the main influence angle tangent; P refers to the mining area. Consequently, this study collected each mining area's overall development plan, geological report, annual mining plans, and mining engineering plan. It used MSPS software to simulate and predict the land subsidence affected by coal mining in each city in 2025 and 2030 based on the mining plan and mineral resource planning of each mine. The damage degree was further separated according to the groundwater level of each mine, the subsidence water accumulation, and the influence of subsidence on vegetation (Table 1). If there is no water accumulation on the surface, it is light damage; if there is seasonal water accumulation on the ground, it is minor damage; if there is perennial water accumulation, it is severe damage.

Damaged Farmland Forecast
Using the spatial analysis function of ArcGIS, the predicted coal mining subsidence data, water accumulation data, and coal-grain composite area were superimposed, and the composite factor was extracted to calculate the damaged cropland area by coal mining.

Construction of a Model for Farmland Yield Caused Coal Mining Subsidence
Severe subsidence caused by coal mining leads to cropland failure, while moderate and minor subsidence causes varying degrees of reduction in crop production. Moderate subsidence caused a decrease in grain production between 20% and 40%, whereas minor subsidence caused a decrease in grain production between 10% and 20% [44]. Therefore, combined with the study of cultivated land destruction in the coal-grain composite areas in eastern China [30], the reduction in cultivated land corresponding to severe, moderate, and minor subsidence was 100%, 30%, and 15%, respectively. Based on the above analysis, we constructed an estimation model (Equation (4)) to approximate the impact of cultivated land damage on grain output within in study area: where Q represents the quantity of grain reduction caused by subsidence in cropland; E represents the average grain yield (valued to 6329 kg/ha); S 1 , S 2 , and S 3 represent the area of minor, moderate, and severe cropland damage, respectively; P 1 , P 2 , and P 3 represent the proportion of grain reduction in minor, moderate, and severe subsided in the cropland, respectively.

Results
Based on the survey data analysis of farmland damage status and the prediction of mining subsidence according to the two stages of 2025 and 2030, the evolution law of farmland damage in composite areas of Shandong Province was intuitively reflected. Then, the grain reduction and other adverse effects were analyzed to support future cropland protection and subsidence land management planning.

Current Distribution of Damaged Farmland
As revealed by survey statistics, as of the end of 2021, the land subsidence area affected by underground coal mining in Shandong Province had reached 8.26 × 10 4 ha. Regarding land subsidence (Figure 3), cropland damage is distributed in Jining City, Tai'an City, Zaozhuang City, Jinan City, Heze City, Yantai City, and Dezhou City. Among them, Jinan City and Tai'an City have the most extensive damage range, and about 70% of the damaged farmland in Shandong Province focuses on these two cities.

Spatiotemporal Distribution of Coal Mining Subsidence Prediction
According to the mining plans of various mines in Shandong Province, MSPS is used to predict the impact of coal mining subsidence in 2025 and 2030 (Figures 4 and 5). It can be seen that from 2025 to 2030, the subsidence area caused by coal mining will gradually expand. Among them, the minor damage range gradually decreases, while the moderate and severe damage range continues to expand. As shown in Table 2, the damaged area of cultivated land in overlapped areas was as high as 6.40 × 10 4 ha, including 5.00 × 10 4 ha of reduced-yield farmland and 1.40 × 10 4 ha of non-productive farmland. Within the entire province's territory, disaster-affected cultivated land and non-yield farmland in Jining City are the largest, which are 2.90 × 10 4 ha and 0.84 × 10 4 ha, respectively. The damaged area accounts for 45.37% of the province's total area of damaged farmland. It is followed by the damaged cultivated land area of 1.53 × 10 4 ha in Tai'an City, accounting for 23.91% of the total area of damaged farmland in the province. The reduced and no-production farmland areas are 1.36 × 10 4 ha and 0.17 × 10 4 ha, respectively. The whole damaged area of Zaozhuang is 0.61 × 10 4 ha, of which the reduced production area is 0.46 × 10 4 ha, and the extinct production area is 0.15 × 10 4 ha. The whole broken area in Jinan City is 0.65 × 10 4 ha. Although it is not much different from the damaged area in Zaozhuang City, its reduced production area of 0.64 × 10 4 ha is much larger than that in Zaozhuang City. Heze City and Yantai City are slightly damaged, with the damaged area accounting for 6.12% and 3.96% of the total damaged area. However, the cropland damaged area, reduced-productive area, and non-productive area in Heze City are higher than those in Yantai City. These areas are 0.25 × 10 4 ha, 0.20 × 10 4 ha, and 0.05 × 10 4 ha in Yantai City. The area of cultivated land destruction in Dezhou City is the smallest, at 0.06 × 10 4 ha. Still, the non-productive farmland area in the city is as high as 0.04 × 10 4 ha, accounting for 64.96% of the total damaged farmland area in the city.

Spatiotemporal Distribution of Coal Mining Subsidence Prediction
According to the mining plans of various mines in Shandong Province, MSPS is used to predict the impact of coal mining subsidence in 2025 and 2030 (Figures 4 and 5). It can be seen that from 2025 to 2030, the subsidence area caused by coal mining will gradually expand. Among them, the minor damage range gradually decreases, while the moderate and severe damage range continues to expand.   As shown in Table 3, the damage to the land by coal mining subsidence gradually increases. By 2025, it is estimated that the subsidence area of the study area will reach 10.27 × 10 4 ha, accounting for 0.65% of the province's total land area. Among them, the severe damage area is 0.82 × 10 4 ha, accounting for 0.82% of the total cropland damaged area. By 2030, the subsidence area of the district will reach 12.26 × 10 4 ha, accounting for 0.78% of the total land area. As seen in Table 3, the severe damage area is 1.37 × 10 4 ha, accounting for 11.18% of the total cropland-damaged area. It follows that from 2021 to 2030, the subsidence area caused by coal mining will gradually expand. The proportion of minor damage area to total damaged area decreased, while the moderate and severe damage area increased.

Spatiotemporal Distribution of Farmland Damage Prediction
The subsidence prediction results and subsidence classification data in 2025 and 2030 were superimposed over the coal-grain composite area to obtain the data on cultivated land affected by coal mining, as shown in Figures 6 and 7 and Table 4. It can be concluded that, from 2025 to 2030, the total damaged cropland area shows an increasing trend, and the damaged area of each damage degree is also expanding. By 2025, the damaged cropland area will reach 6.97 × 10 4 ha, of which the severely damaged area will be 1.63 × 10 4 ha, accounting for 23.44% of the total damaged cropland area. By 2030, the damaged cropland area will reach 7.52 × 10 4 ha, of which the severely damaged area will be 1.89 × 10 4 ha, accounting for 25.12% of the total damaged cropland area. The results indicate that by 2030, about 7.39% of the cultivated land in Shandong Province will be damaged due to coal mining, and the damage will become more serious. Furthermore, about 1.86% of the arable land will be submerged in water, forming perennial stagnant water areas that are unable to produce food. damaged farmland area in Jining City was the largest, with an area of 3.17 × 10 4 ha. It is followed by Tai'an City, with an area of 1.51 × 10 4 ha, and the smallest is Dezhou City, which has an area of 0.08 × 10 4 ha. During 2025-2030, the damaged farmland area in Jining City is expected to increase the most, followed by Heze City, with Yantai City increasing the least. Compared with the period from 2021 to 2025, except for Jining City, Heze City, and Dezhou City, the area of each damage degree has changed significantly, and other cities have not changed much.     There are also differences in the spatial distribution of cropland damage in different periods (Tables 5 and 6). Except for Jinan City, which stopped coal mining in 2021, other cities will continue to expand coal mining from 2021 until at least 2030. Before 2025, the damaged farmland area in Jining City was the largest, with an area of 3.17 × 10 4 ha. It is followed by Tai'an City, with an area of 1.51 × 10 4 ha, and the smallest is Dezhou City, which has an area of 0.08 × 10 4 ha. During 2025-2030, the damaged farmland area in Jining City is expected to increase the most, followed by Heze City, with Yantai City increasing the least. Compared with the period from 2021 to 2025, except for Jining City, Heze City, and Dezhou City, the area of each damage degree has changed significantly, and other cities have not changed much.

Reduced Grain Production
Land subsidence has highly adverse effects on agricultural production. Coal mining destroys the surrounding soil layers, changes the infiltration conditions and evaporation paths of cultivated land, and seriously affects the water environment and nutrient environment of the cultivated soil, resulting in the decline of soil productivity and the failure of crops to grow normally. In areas with high groundwater levels, subsidence caused serious water accumulation in cultivated land and damaged agricultural supporting facilities. These disruptions severely affected food production and cash crop planting, leading to reduced agricultural output or, in some cases, total crop failure. Apart from some farmers being forced to adjust from their previous planting industry to aquaculture, many farmers have no growing land, and agricultural economic losses are profound.
According to the analysis of the forecast results and the estimation of the grain production reduction model (Equation (4)), by 2025 and 2030, grain production will be reduced by 16.44 × 104 t and 18.45 × 104 t, respectively. Figures 8 and 9 show that from 2025 to 2030, the grain yield reduction in Jining City will change the most. Except for Jinan City, which will see no change in grain production due to no longer mining after 2021, there will be no change in grain yield reduction in Yantai City and the extent of partial damage in other cities, mainly because the damaged area will change little. It is estimated that by 2030, minor, moderate, and severe damage will cause a reduction in grain output of 4.19 × 10 4 t, 2.30 × 10 4 t, and 11.95 × 10 4 t, accounting for 22.73%, 12.47%, and 64.80% of the total decline in output, respectively. In summary, local grain production will be greatly reduced due to continuous mining, which will significantly hinder the development of agricultural production.
will be no change in grain yield reduction in Yantai City and the extent of partial damage in other cities, mainly because the damaged area will change little. It is estimated that by 2030, minor, moderate, and severe damage will cause a reduction in grain output of 4.19 × 10 4 t, 2.30 × 10 4 t, and 11.95 × 10 4 t, accounting for 22.73%, 12.47%, and 64.80% of the total decline in output, respectively. In summary, local grain production will be greatly reduced due to continuous mining, which will significantly hinder the development of agricultural production.

Change of Farmland Landscape Pattern
With the disturbance of coal mining, the surface collapsed on a large scale ( Figure  10). Contiguous deep subsidence basins were formed in some areas such as Jining City, Zaozhuang City, and Heze City, resulting in the gradual evolution of flat farmland into sloping fields, perennial waterlogging areas, seasonal waterlogging areas, etc. Continuous coal production will inevitably lead to larger-scale surface subsidence, and the farmland landscape will change accordingly. By the end of 2025, the accumulated water area of cultivated land in the coal-grain compound zone in Shandong Province will be 2.71 × 10 4 ha, and by the end of 2030, the accumulated water area will increase to 3.10 × 10 4 ha. Most of the increased water is coming from cropland. Cropland has lost its landscape advantages, and water landscape advantages have risen significantly. Currently, Shandong Province is optimizing and adjusting the layout of coal resource development in the province and reducing the existing production capacity. In the future, surface subsidence will slowly increase, and the extent of cultivated land evolving into stagnant water and converging with sloping cropland will also increase.

Change of Farmland Landscape Pattern
With the disturbance of coal mining, the surface collapsed on a large scale (Figure 10). Contiguous deep subsidence basins were formed in some areas such as Jining City, Zaozhuang City, and Heze City, resulting in the gradual evolution of flat farmland into sloping fields, perennial waterlogging areas, seasonal waterlogging areas, etc. Continuous coal production will inevitably lead to larger-scale surface subsidence, and the farmland landscape will change accordingly. By the end of 2025, the accumulated water area of cultivated land in the coal-grain compound zone in Shandong Province will be 2.71 × 10 4 ha, and by the end of 2030, the accumulated water area will increase to 3.10 × 10 4 ha. Most of the increased water is coming from cropland. Cropland has lost its landscape advantages, and water landscape advantages have risen significantly. Currently, Shandong Province is optimizing and adjusting the layout of coal resource development in the province and reducing the existing production capacity. In the future, surface subsidence will slowly increase, and the extent of cultivated land evolving into stagnant water and converging with sloping cropland will also increase. ha, and by the end of 2030, the accumulated water area will increase to 3.10 × 10 4 ha. Most of the increased water is coming from cropland. Cropland has lost its landscape advantages, and water landscape advantages have risen significantly. Currently, Shandong Province is optimizing and adjusting the layout of coal resource development in the province and reducing the existing production capacity. In the future, surface subsidence will slowly increase, and the extent of cultivated land evolving into stagnant water and converging with sloping cropland will also increase.

Damage to Ecological Environment
As shown in Figure 11, during the coal mining process, it is easy to destroy the land structure of the mining area and its surrounding areas. The collapse of coal mining causes slopes in formerly flat farmland, disturbing the originally relatively stable soil structure and geological environment [45][46][47][48]. Slight neglect of coal mining may adversely affect

Damage to Ecological Environment
As shown in Figure 11, during the coal mining process, it is easy to destroy the land structure of the mining area and its surrounding areas. The collapse of coal mining causes slopes in formerly flat farmland, disturbing the originally relatively stable soil structure and geological environment [45][46][47][48]. Slight neglect of coal mining may adversely affect the groundwater system and destroy the balance of the groundwater system, resulting in surface water leakage in mining areas and causing extensive water accumulation in damaged cropland [49]. In addition to the combined effects of topography and climate, the vegetation in the goaf will also be severely damaged. The vegetation in the subsidence areas will take a long time to recover [50], causing the ecological environment to be forced into a vicious circle, bringing severe hidden dangers to mining area development.
re 2022, 12, x FOR PEER REVIEW 13 of 17 the groundwater system and destroy the balance of the groundwater system, resulting in surface water leakage in mining areas and causing extensive water accumulation in damaged cropland [49]. In addition to the combined effects of topography and climate, the vegetation in the goaf will also be severely damaged. The vegetation in the subsidence areas will take a long time to recover [50], causing the ecological environment to be forced into a vicious circle, bringing severe hidden dangers to mining area development. This study predicts that by 2025, an area of 2.71 × 10 4 ha of cropland ecosystem will be transformed into a wetland ecosystem in Shandong province due to coal mining if no relevant control measures are taken. By 2030, the study area will have an area of about 3.10 × 10 4 ha of cropland ecosystems transformed. The balance of the original cropland ecosystem will be disrupted, and the system energy flow and logistics will be blocked by the external circulation, resulting in a drastic reverse change in the cropland ecological environment. Living organisms will lose their habitat and niches, crops will not grow normally, and cropland ecological safety will decline.

Impact of Farmland Damage on Social Stability
Shandong Province, as a typical coal-grain complex area where farmers' income, agricultural product supply, mine output, and miner safety are the most important areas of the 'three rural' and 'four mines' issues. After coal mining, large cultivated land areas lie under accumulated water, topsoil resources are eroded, and soil organic matter is lost [51][52][53]. It is estimated that the soil loss will increase by about 20% for every 1° increase in the original grainland slope. At the same time, the tillage layer becomes shallow, generally from 25 cm to 15 cm. Compared with the normal cultivated land, the organic matter content of the subsidence-affected cultivated land decreases by about 30% [31]. Farmland damage, area reduction, and productivity decline lead to many farmers' income reduction and employment pressure increase. Surface subsidence and surface cracks caused by mining destroy buildings, traffic roads, and communication lines and greatly inconvenience farmers' production and life [54,55]. On issues such as compensation for land damage, land acquisition, land reclamation, village demolition, and personnel resettlement, local governments, mining enterprises, and farmers continue to have disputes and prominent disagreements, which can seriously affect the harmony and stability of local communities.
Because Shandong Province is also a high groundwater level area, damaged farmland will appear to be waterlogging and even become a perennial waterlogging area. As a result, the reclamation cost increases, and the reclamation rate decreases, which brings many inconveniences to the reclamation work.
In conclusion, considering the current scale of coal mining, the subsidence area will This study predicts that by 2025, an area of 2.71 × 10 4 ha of cropland ecosystem will be transformed into a wetland ecosystem in Shandong province due to coal mining if no relevant control measures are taken. By 2030, the study area will have an area of about 3.10 × 10 4 ha of cropland ecosystems transformed. The balance of the original cropland ecosystem will be disrupted, and the system energy flow and logistics will be blocked by the external circulation, resulting in a drastic reverse change in the cropland ecological environment. Living organisms will lose their habitat and niches, crops will not grow normally, and cropland ecological safety will decline.

Impact of Farmland Damage on Social Stability
Shandong Province, as a typical coal-grain complex area where farmers' income, agricultural product supply, mine output, and miner safety are the most important areas of the 'three rural' and 'four mines' issues. After coal mining, large cultivated land areas lie under accumulated water, topsoil resources are eroded, and soil organic matter is lost [51][52][53]. It is estimated that the soil loss will increase by about 20% for every 1 • increase in the original grainland slope. At the same time, the tillage layer becomes shallow, generally from 25 cm to 15 cm. Compared with the normal cultivated land, the organic matter content of the subsidence-affected cultivated land decreases by about 30% [31]. Farmland damage, area reduction, and productivity decline lead to many farmers' income reduction and employment pressure increase. Surface subsidence and surface cracks caused by mining destroy buildings, traffic roads, and communication lines and greatly inconvenience farmers' production and life [54,55]. On issues such as compensation for land damage, land acquisition, land reclamation, village demolition, and personnel resettlement, local governments, mining enterprises, and farmers continue to have disputes and prominent disagreements, which can seriously affect the harmony and stability of local communities.
Because Shandong Province is also a high groundwater level area, damaged farmland will appear to be waterlogging and even become a perennial waterlogging area. As a result, the reclamation cost increases, and the reclamation rate decreases, which brings many inconveniences to the reclamation work.
In conclusion, considering the current scale of coal mining, the subsidence area will be further increased in the future. The damage of subsidence to land resources, the negative impact on agriculture, ecology, and society, and the negative impact on the production, living, and living of residents in the subsidence area will be further increased. Mining areas should combine with the actual local situation and implement green mining technology to reduce these negative impacts caused by surface subsidence from the source [56,57]. According to the predicted damage to cropland, the mining areas can choose appropriate reclamation time, elevation, and layout to simulate the reclamation scenario while mining and change the traditional land reclamation method [58]. With the goal of restoring grainland, increasing reclamation rate, and reducing reclamation cost, the main measures for cultivated land protection and reclamation are proposed by combining natural restoration and engineering measures. Facing other problems such as water accumulation and facilities destruction in subsidence areas, the government and mining area should reasonably plan and design an ecological reconstruction scheme from the aspects of water accumulation degree, waste utilization, and farmers' demand to create an efficient ecological agriculture landscape [59].

Limitations and Outlooks
There is an interaction effect between coal mining and grain production in COA. The time continuity, spatial scalability, and substantial interference of coal mining have seriously threatened the high-quality cultivated land. What makes it worse is that such a threat will be intensified in time and space dimensions, thus threatening the sustainable development of the mining area. Therefore, it is key to coordinate coal mining and farmland protection to study cropland damage and its impact on the overlapped areas. In the past, the prediction of cropland destruction was mostly based on empirical models [20,31], or the probability integral method was used to predict cropland damage on a small scale [35]. Therefore, in this paper, the provincial scale is selected to predict the change of cropland destruction with the coal mining development in overlapped areas. The analysis of the current situation and future trend of cropland damage indicates that damage caused by long-term, large-scale, and high-intensity coal mining continues to increase, and the damage degree continues to change. This result is consistent with the characteristics of cropland damage caused by coal mining in the region [60]. In addition, previous studies revealed the impact of cropland damage using quantitative or qualitative analysis methods from the soil, vegetation, agriculture, ecological environment, and other aspects [61][62][63][64]. These provide the foundation for this study to analyze the impact of coal mining on cultivated land in Shandong Province. Through calculation, investigation, and literature analysis, it was found that the ever-changing scope of grainland damage will cause changes in the soil quality, landscape pattern, and ecological environment quality in the region and bring many social problems.
However, due to the complexity of the evolution process of cropland damage, certain errors will occur when predicting its evolution law. In addition, the change of cropland landscape patterns and ecological environment is a relatively complex process, which is affected not only by coal mining but also by geological conditions, climate, human activities, and other factors in the coal mining area. The factors of cultivated land destruction also involve many aspects. Therefore, it is worthy of further study to fully consider the influencing factors of cropland damage in the coal-grain composite areas and reduce the error between the predicted area and the actual damage. Meanwhile, it is necessary to analyze further the impact of cultivated land destruction on mining areas and its impact mechanism.

Conclusions
Coal and cropland resources in COA have a large overlapped area and a wide distribution range. Coal mining will inevitably cause a large amount of farmland damage. Based on the provincial scale, taking Shandong Province as an example, this paper analyzed the spatiotemporal impact of coal mining on cropland in coal-cropland overlapping areas. The research results suggest that the cropland area in overlapped regions will continue to decrease with the continuous exploitation of coal resources. The damage degree will continue to increase, seriously affecting local agriculture, the environment, and social stability. The research results lay a foundation for other coal-grain complex areas to clarify the current situation, trend, and mechanism of cropland damage. At the same time, grasping the temporal and spatial impact of mining on cropland in overlapped areas is also the premise and basis for formulating coordinated countermeasures for mining and cropland protection.
Therefore, it is recommended that relevant departments consider the following points when formulating grainland protection plans: (1) overlapping characteristics and evolution trends of damaged cropland should be fully considered. (2) Appropriate reclamation methods should be selected. Concurrent mining and reclamation technology can be promoted for areas with high groundwater levels to save more topsoil resources in advance. (3) It is necessary to rebuild the ecological environment according to local conditions. In addition, we should vigorously promote the green mining technology of coal mines to reduce the impact of coal mining subsidence from the source.
Given the limitations of this study, future research needs to pay more attention to how to reduce the error between the predicted damage and the actual damage and explore the impact mechanism of farmland damage on the mining area development to realize the sustainable development of coal mining and farmland protection.