Abstract:Ecosystem service bundle is the combination of multiple ecosystem services, and it could serve as the dominant function characterization of ecosystem multifunctionality. The identification of the spatial and functional changes and the driving factors of ecosystem service bundles can provide reference for dominant functions based zoning management strategies, land versatility coordination, and integrative promotion of regional ecosystem services. Taking the Daqing River Basin as an example, this study selected six ecosystem services as criteria based on the actual situation, including two types of supply services (water yield and crop production) and four kinds of regulation services (water retention, water purification, soil retention, and carbon sequestration). On the township scale, the characteristics of ecosystem service bundle dynamics were discussed, and the main driving factors were investigated. The results indicated that:(1) According to the dominant ecosystem service, the Daqing River Basin can be divided into three ecosystem service bundles, namely, ecological regulation bundle (B1), agricultural product supply bundle (B2) and human settlement bundle (B3); (2) From 2000 to 2015, the spatial pattern of ecological regulation bundle changed the most greatly, with a change rate of 19.6%, while the agricultural product supply bundle and human settlement bundle had relatively strong spatial stability, with a change rate lower than 5.0%. The changes of ecosystem service bundle types mainly happened at the junction of different bundles; (3) The ecosystem services of all the three bundles showed overall increasing trends simultaneously, with an average increase rate of 8.62%; (4) Initial status of natural background and ecological protection policies were the main factors driving the change of ecological regulation bundle. Meanwhile, the temporal and spatial changes of agricultural product supply bundle and human settlement bundle were mainly driven by agricultural products supply and built-up land area, respectively. The results of this study could provide effective approach to manage the multi-functional properties of the land and enhance the overall function of territory space.