Recognizing changes in ecosystem services (ES) and their relationships is the basis of achieving sustainable regional development. Regional collaborative development has become the core strategy of the development of the Beijing–Tianjin–Hebei (BTH) region. However, sub regions have different ecological changes and relationships. Here, we quantify and map ES, including water yield, sediment retention, carbon sequestration and grain productive capacity in 2000, 2005, 2010 and 2015, using several biophysical models and explore the relationships of spatial correction, trade-offs and synergies among multiple ES in different spatial scales. Results across the four years show that the quality and variation tendency of ES from each region are spatially heterogeneous. The relationship between ES that are not significant in the entire region shows different correlations in individual ecological–functional zones. From the perspective of regional disparity, the effect of land use factor and correlative mechanisms among ES are analyzed. To observe the spatiotemporal variations and relationships of ES in individual regions, land use management policies are proposed on the basis of the results of the relationships among ES.