Environmental managers and policymakers increasingly discuss trade-offs between ecosystem services (ESs). However, few studies have used nonlinear models to provide scenario-specific land-use planning. This study determined the effects of different future land use/land cover (LULC) scenarios on ESs in the Yili River Valley, China, and analyzed the trade-offs and synergistic response characteristics. We simulated land-use changes in the Yili River Valley during 2020–2030 under three different scenarios using a patch-generating land-use simulation (PLUS) model—business as usual (BAU), economic development (ED), and ecological conservation (EC). Subsequently, we evaluated the water yield (WY), carbon storage (CS), soil retention (SR), and nutrient export (NE) ESs by combining the PLUS and integrated valuation of ecosystem services and trade-offs (InVEST) models, thus exploring multiple trade-offs among these four ESs at a regional scale. For the BAU scenario, there are some synergistic effects between WY and SR in the Yili River Valley, in addition to significant trade-off effects between CS and NE. For the ED scenario, the rapid expansion of cropland and constructed land is at the expense of forested grassland, leading to a significant decline in ESs. For the EC scenario, the model predicted that the cumulative regional net future carbon storage, cumulative water retention, and cumulative soil conservation would all increase due to ecological engineering and the revegetation of riparian zones and that formerly steep agricultural land can be effective in improving ESs. Meanwhile, the trade-off effect would be significantly weakened between CS and NE. These results can inform decision makers on specific sites where ecological engineering is implemented. Our findings can enhance stakeholders’ understanding of the interactions between ESs indicators in different scenarios.

Pakistan is a flood-prone country and almost every year, it is hit by floods of varying magnitudes. This study was conducted to generate a flash flood map using analytical hierarchy process (AHP) and frequency ratio (FR) models in the ArcGIS 10.6 environment. Eight flash-flood-causing physical parameters were considered for this study. Five parameters were based on the digital elevation model (DEM), Advanced Land Observation Satellite (ALOS), and Sentinel-2 satellite, including distance from the river and drainage density slope, elevation, and land cover, respectively. Two other parameters were geology and soil, consisting of different rock and soil formations, respectively, where both layers were classified based on their resistance against water percolation. One parameter was rainfall. Rainfall observation data obtained from five meteorological stations exist close to the Chitral District, Pakistan. According to its significant importance in the occurrence of a flash flood, each criterion was allotted an estimated weight with the help of AHP and FR. In the end, all the parameters were integrated using weighted overlay analysis in which the influence value of the drainage density was given the highest value. This gave the output in terms of five flood risk zones: very high risk, high risk, moderate risk, low risk, and very low risk. According to the results, 1168 km2, that is, 8% of the total area, showed a very high risk of flood occurrence. Reshun, Mastuj, Booni, Colony, and some other villages were identified as high-risk zones of the study area, which have been drastically damaged many times by flash floods. This study is pioneering in its field and provides policy guidelines for risk managers, emergency and disaster response services, urban and infrastructure planners, hydrologists, and climate scientists.