Abstract: The Shendong mining area is an important coal production base in China,and the coal mining has caused strong disturbance to its fragile ecological environment.It is necessary to use scientific methods to reveal the evolution of the ecosystem service function of the Shendong mining area,and improve the management system of green mining and sustainable development in the mining area.To quantify the changes in ecosystem service functions in Shendong mining area,a Comprehensive Ecosystem Service (CES) evaluation model was constructed based on the InVEST model,which integrates three ecosystem services:water supply,soil conservation,and carbon storage,and measures the characteristics and patterns of changes in CES in the Shendong mining area from 1990 to 2018 in both the scale of mining area and mine.The results show that at the mining area scale,① the SCE in Shendong mining area decreased from 0.448 6 (1990) to 0.382 5 (2000),and then continued to increase to 0.471 6 (2015) and then decreased to 0.453 2 (2018),showing an overall increasing trend at a rate of 0.009 a-1.Spatially,it was the weakest in the northwest and gradually strengthened from the northwest to the southeast.② The spatial variability of the CES in the Shendong mining area was obvious,and the contribution of the change in the low-variable area to the change in the CES was higher than that of the medium-variable area and the abrupt change area,and the CES spatial change was closely related to vegetation cover and land use types.③ The spatial clustering of CES in the Shendong mining area was obvious,showing a pattern of “weakest in the northwest,strongest along the rivers,and insignificant in the southeast”.④ The driving factors of CES change in Shendong mining area were land use type,slope,elevation,rainfall and vegetation coverage in order.The synergy of two factors was higher than that of single factor,and the synergy between land use and other factors was the strongest.At the mine scale,① the CES in the Shendong mine area was highest in the medium-intensity mining area (SCE=0.502 8) and lowest in the very high-intensity mining area (SCE=0.430 8).② The CES in the Bulianta mine area was positively influenced by ecological management measures,the CES in the Daliuta and Huojitu mine area had been negatively affected by mining since 2010,and the CES in the Shigetai,Wulanmulun,Halagou and Yujialiang mine area had continued to suffer negative impacts since mining.③ The CES in the Daliuta rec-lamation experimental area had developed in the direction of becoming stronger.The model’s comprehensive as-sessment results can quantitatively reflect the changing pattern of the CES in the Shendong mining area.It is suggested that Shendong mining area should implement ecological management with the northwest as the key area,and focus on optimizing the land use structure by improving vegetation cover,and use scientific methods such as microbial reclamation in the mining area.In addition,the development and use of multi-scale integrated model for the optimization of multiple ecosystem services should be committed,with a view to providing more accurate and reliable decision support for the sustainable management of regional ecosystem services.