Spatiotemporal changes and driving factors of ecosystem services between karst and non-karst World Heritage sites in Southwest China

Understanding the spatiotemporal variation and drivers of ecosystem services is fundamental to optimal management and sustainable development of World Heritage (WH) sites. Although WH sites face multiple natural and anthropogenic threats, our understanding of their ecosystem services is still limited, especially for karst WH sites. In this study, we assessed habitat quality (HQ), carbon storage (CS), soil retention (SR), water conservation (WC), and the combined ecosystem service (CES) of karst and non-karst WH sites in Southwest China from 2000 to 2020 using the InVEST model. We also assessed trade-offs/synergies among ecosystem services using the spatial overlay method, and identified driving factors of variation in ecosystem services using geographical detector and structural equation models. The results showed that ecosystem services of the WH sites exhibited high spatiotemporal variation. In particular, there were higher values in the property zone than in the buffer zone, and an increasing trend in SR but a decreasing trend in HQ and CES over time. Compared to non-karst sites, karst WH sites had significantly lower values of HQ, CS, SR, and CES, but higher spatial heterogeneity in CS, WC, and CES. Weak trade-offs among ecosystem services dominated the WH sites, with the proportion of weak synergies increasing over time. Compared to non-karst sites, karst WH sites had a significantly lower proportion of strong synergies and a significantly higher proportion of weak synergies. The provision of ecosystem services was primarily influenced by natural factors (e.g., landscape division index and normalized difference vegetation index), followed by anthropogenic factors (e.g., distance from road and population density). Overall, these findings may have important implications for decision-making aimed at protecting the outstanding universal value, authenticity, and integrity of WH with different attributes.