A novel ecological risk assessment model of ecosystem services was developed.
DNDC, SWAP, Biome-BGC and AEZ models were used for assimilation.
Integrated degradation risk of four key ecosystem services was assessed.
Land-use change posed a great degradation risk on ecosystem services in 1985–2020.
Severe deterioration in 2020 would be centered in some small and less developed cities.