Abstract
Understanding the dynamic changes of relationships between ecosystem services (ESs) and their dominant factors can effectively adjust human activities to adapt proactively to global climate change. In this study, the Huang-Huai-Hai Plain (HHHP) was selected to assess the dynamics of four key ESs (NPP, net primary productivity; WY, water yield; SC, soil conservation; FP, food production) from 2000 to 2020. The constraint lines of interactions among ESs were extracted based on a segmented quantile regression model. On this basis, the effects of both human activities and natural factors on the key features of the interactions between ESs were quantified with the help of automatic linear model. The results indicated that two types of constraint relationships, including exponential and humped-shaped, existed among the six pairs of ESs. In the past two decades, small changes in NPP thresholds would lead to large variations in other ESs thresholds. Precipitation and normalized difference vegetation index were the key factors to determine the constraint strength of ESs in the HHHP. The potential maximum value of WY in the HHHP could be increased by adjusting landscape shape to make it more complicated. This study helps to improve the potential of target ESs and provides a decision-making basis for promoting regional sustainable development.