Spatial–Temporal Multivariate Correlation Analysis of Ecosystem Services and Ecological Risk in Areas of Overlapped Cropland and Coal Resources in the Eastern Plains, China

The overlapped areas of cropland and coal resources play a fundamental role in promoting economic and social progress. However, intensive mining operations in high water-level areas have brought significant spatial–temporal heterogeneity and ecological problems. From the dual dimensions of the ecosystem service value (ESV) and ecological risk (ER), it is of great significance to explore the influence characteristics of underground mining on the landscape, such as above-ground cultivated land, which is valuable to achieving regional governance and coordinated development. In this study, taking Peixian as the research area, a multiple-dimensional correlation framework was constructed based on the revised ESV and ER, integrating the grey relational degree, spatial–temporal heterogeneity, disequilibrium, and inconsistency index to explore the ESV and ER assessment and correlation characteristics from 2010 to 2020. The results show that (1) the ESV showed a high agglomerated distribution pattern in the east, with a net decrease of 13.61%. (2) The ER decreased by 78.18 and was concentrated in the western and southern regions, with overall contiguous and local agglomeration characteristics. This indicates that the ecological security of the region has improved. (3) The comprehensive grey correlation between the cultural service value and the ecological risk index was the highest. Furthermore, the spatial–temporal heterogeneity of the ESV and ER weakened, and the disequilibrium rose and then fell, indicating that the ecosystem gradually tended to be stable. The study is crucial for overlapped cropland and coal resource areas to maintain stability and sustainable development. The multivariate correlation framework provides practical value for ecosystem management and risk control.