Simulating changes in the value of ecosystem services caused by land use changes in large cities under multiple scenarios is of great significance for cities to formulate land use policies and improve ecosystem services. Take Shijiazhuang, which is in the process of rapid urbanization, as an example. Based on the remote sensing image data and statistical yearbook of 1988, 1998, 2008, and 2018 as the basic data to analyze and estimate the 30 years of land use and ecosystem service value changes in Shijiazhuang. According to this, the CA-Markov model was used to simulate the land use change in Shijiazhuang under three scenarios in 2030 and estimate the value of ecosystem services under each scenario, using grid tools to visually express the spatial distribution of ecosystem service values and the degree of agglomeration under three scenarios. The results indicate that the most obvious feature of land use change in Shijiazhuang from 1988 to 2018 was that the farmland area decreased year by year, the built-up expanded rapidly, the farmland area decreased by 86,874.75 hm2 in 30 years, and the built-up increased by 154,711.90 hm2. In 1988, 1998, 2008, and 2018, the ecosystem service value of Shijiazhuang was 32.578 billion yuan, 32.799 billion yuan, 29.944 billion yuan, and 31.251 billion yuan respectively. In 2030, under three scenarios of natural development, farmland protection, and ecological protection, the value of ecosystem services is 331.111 billion yuan, 33.670 billion yuan, and 33.891 billion yuan in order. The hot spots are mainly concentrated in the northwest and southwest of Shijiazhuang, and cold spots are concentrated in the eastern cities, counties, and districts. Based on changes in land use brought about by urban expansion, simulating the value of ecosystem services under multiple scenarios in the future, providing scientific guidance for building urban ecological networks, and realizing sustainable urban ecological development.
Keywords: CA-Markov; Change of land use; Ecosystem service value; Hot spot analysis.