Identifying the impacts of land use landscape pattern and climate changes on streamflow from past to future

Identifying the individual and combined hydrological response of land use landscape pattern and climate changes is key to effectively managing the ecohydrological balance of regions. However, their nonlinearity, effect size, and multiple causalities limit causal investigations. Therefore, this study aimed to establish a comprehensive methodological framework to quantify changes in the landscape pattern and climate, evaluate trends in streamflow response, and analyze the attribution of streamflow events in five basins in Beijing from the past to the future. Future climate projections were based on three general circulation models (GCMs) under two shared socioeconomic pathways (SSPs). Additionally, the landscape pattern in 2035 under a natural development scenario was simulated by the patch-generating land use simulation (PLUS). The Soil and Water Assessment Tool (SWAT) was applied to evaluate the streamflow spatial and temporal dynamics over the period 2005-2035 with multiple scenarios. A bootstrapping nonlinear regression analysis and boosted regression tree (BRT) model were used to analyze the individual and combined attribution of landscape pattern and climate changes on streamflow, respectively. The results indicated that in the future, the overall streamflow in the Beijing basin would decrease, with a slightly reduced peak streamflow in most basins in the summer and a significant increase in the autumn and winter. The nonlinear quadratic regression more effectively explained the impact of landscape pattern and climate changes on streamflow. The trends in the streamflow change depended on where the relationship curve was in relation to the threshold. In addition, the impacts of landscape pattern and climate changes on streamflow were not isolated but were joint. They presented a nonlinear, non-uniform, and coupled relationship. Except for the YongDing River Basin, the annual streamflow change was influenced more by the landscape pattern. The dominant factors and the critical pair interactions varied from basin to basin. Our findings have implications for city planners and managers for optimizing ecohydrological functions and promoting sustainable development.