RESUMO
The effects of global climate change and human activities are anticipated to significantly impact ecosystem services (ESs), particularly in urban agglomerations of arid regions. This paper proposes a framework integrating the dynamic Bayesian network (DBN), system dynamics (SD) model, patch generation land use simulation (PLUS) model, and the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) model for predicting land use change and optimizing ESs spatial patterns that is built upon the SSP-RCP scenarios from CMIP6. This framework is applied to the oasis urban agglomeration on the northern slope of the Tianshan Mountains in Xinjiang (UANSTM), China. The findings indicate that both the SD model and PLUS model can accurately forecast the distribution of future land use. The SD model shows a relative error of less than 2.32%, while the PLUS model demonstrates a Kappa coefficient of 0.89. The land use pattern displays obvious spatial heterogeneity under different climate scenarios. The expansion of cultivated land and construction land is the main form of land use change in UANSTM in the future. The DBN model proficiently simulates the interactive relationships between ESs and diverse factors. The classification error rates for net primary productivity (NPP), habitat quality (HQ), water yield (WY), and soil retention (SR) are 20.04%, 3.47%, 4.45%, and 13.42%, respectively. The prediction and diagnosis of DBN determine the optimal ESs development scenario and the optimal ESs region in the study area. It is found that the majority of ESs in UANSTM are predominantly influenced by natural factors with the exception of HQ. The socio-economic development plays a minor role in such urban agglomerations. This study offers significant insights that can contribute to the fields of ecological protection and land use planning in arid urban agglomerations worldwide.