Abstract: Although coal has made a huge contribution to the development of the economy and socie-ty and its economic benefits have often attracted much attention, little research has focused on the ecosystem services of coalfields. Based on remote sensing data, meteorological data, and soil data in Shanxi coalfields during 1986, 2000, and 2015, we estimated soil conservation and water yield using the InVEST model, assessed the net primary productivity of vegetation using the CASA mo-del, and estimated sand fixation using the RWEQ model. Further, we simulated the spatial patterns of ecosystem services (ESs) using the k-means cluster analysis method and analyzed the influence factors of ESs using the Geodetector model in Shanxi coalfield areas. The results showed that soil conservation service, water yield service, and sand fixation service increased continuously. The high-value area of soil conservation service was mainly concentrated in the north of Hedong coalfield and the northeast of Qinshui coalfield, while the low-value area was distributed in the southwestern edge of Datong coalfield. The high-value area of water yield service was mainly concentrated in the northeast of Qinshui coalfield, while the low-value area was distributed in the northeast of Qinshui coalfield, Xishan coalfield and northwestern Qinshui coalfield. The high-value area for vegetation production service was mainly concentrated in the southeast of Qinshui coalfield, while the low-value area was distributed in Datong coalfield, Ningwu coalfield, Xishan coalfield, and northern Hedong coalfield. The distribution of low- and high-value areas of sand fixation service was unfixed. Ecosystem service bundles could be divided into four categories. The first category belonged to soil conservation service bundle, mainly distributed in the northern Ningwu coalfield, the northern Hedong coalfield, and the northern Qinshui coalfield. The second was water yield service bundle, mainly distributed in Huoxi coalfield and southern Qinshui coalfield. The third category belonged to vegetation production service bundle, mainly distributed in parts of Qinshui coalfield. The fourth category belonged to sand fixation service bundle, mainly distributed in the southern part of Hedong coalfield and Qinshui coalfield. Soil conservation service was greatly affected by temperature, digital elevation model (DEM), and industrial output value, with q values of 0.5, 0.3, and 0.2, respectively. Water yield service was greatly affected by precipitation, temperature, and DEM, with q values of 0.8, 0.3, and 0.2, respectively. The industrial output value, precipitation, and temperature q values of vegetation production service were 0.7, 0.6, and 0.2, respectively. The main influen-cing factors of sand fixation service were precipitation, temperature, and DEM, while the q values were 0.7, 0.3, and 0.3, respectively. The spatial distribution of coalfields ESs and the relationship between multiple ESs were closely related to natural and human factors. Therefore, maintaining the coordination relationship between natural-human factors and ecological services would be helpful to the management of the land reclamation, ecological reconstruction, and the sustainable development of coalfields ecosystem.

Key words: ecosystem service, bundle, driving force, spatio-temporal pattern, coalfield.