Integrated assessment of the supply–demand relationship of ecosystem services in the Loess Plateau during 1992–2015

ABSTRACT Landscape alterations and ecosystem services (ES) are crucial elements that affect the socio-ecological development of ecologically fragile regions. To provide scientific support for the land-use planning and regional ecological restoration to achieve sustainable development goals (SDGs) in the Loess Plateau from 1992 to 2015, this study proposes an integrated assessment framework for assessing the supply – demand match and equilibrium of ES combining potential relationships between landscape alterations, ES, and sustainable development at the county and regional scales. Results show that more than 85% of local counties have a supply – demand mismatch in terms of ecosystem services. Supply – demand distributions for provisioning and supporting services were relatively balanced compared to cultural services, which were imbalanced. Although the overall supply – demand relationship is relatively balanced because of the significant influence of ecological restoration, it shows an increasing trend toward imbalance as human – land conflict is prominent in some regions. SDGs emphasize on specific ES information such as ecosystem conservation, sustainable agriculture, and urban construction. Furtherly, positive impacts from landscape dynamics can improve the supply capacity of the ES and contribute to regional sustainability.


Introduction
Supply -demand equilibrium of ecosystem services (ES) is a central issue in socio-ecological management, which effectively indicates that ES dynamics are affected by landscape patterns, decision-making, and human well-being (Dai et al. 2015;Wei et al. 2017). ES involve various human benefits obtained from ecosystems, thereby reflecting environmental conditions and utility (Ramel et al. 2020). The Millennium Ecosystem Assessment (MEA) classifies ES into provisioning, regulating, cultural, and supporting services, and is one of the most influential classifications and evaluation frameworks for ES (Marques et al. 2013). The term supply refers to the capacity of ecosystems to provide various products and services, whereas the term demand refers to the willingness to acquire or consume multiple products and services (Castillo-Eguskitza, Martín-López, and Onaindia 2018). Generally, the match between supply and demand of ES is analyzed in terms of quantity, quality, time, space, and stakeholders . Due to mismatches between supply and demand of ES, there are quantitative differences and spatial misalignments (Geijzendorffer, Martin-Lopez, and Roche 2015), indicating changes in profit and loss between the provision of ecosystem functions and social needs (Lin et al. 2022). The imbalance or mismatches in the supply and demand of ES severely constrains the sustainability of the socio-ecological system at global and regional scales, indicating that the supply -demand assessment of ES is becoming a key challenge for interdisciplinary research and landscape management decision-making (Schirpke et al. 2019).
Supply -demand studies on ES are rapidly evolving from conceptual and integrated assessment methods to quantitative analysis , and mainstream methods include valuation based on ecological economics (e.g., monetary estimation and ecological footprint) (Costanza et al. 2014;Mancini et al. 2018) and processes-based models (Marques et al. 2013;Othoniel et al. 2019). However, these methods mainly focus on the supply side of the service, while the theoretical and assessment methods for ES demand are not yet fully developed. Assessment methods for ES demand mainly reflect two aspects. First, for a certain ES, supply services mostly use various forms of consumption as demand indicators, such as per capita water and food consumptions in each administrative unit (Yu et al. 2021). For regulation and support services, they use emissions (Xia et al. 2022), importance or vulnerability (Larondelle and Lauf 2016;Feurer et al. 2021), and accessibility (Ala-Hulkko et al. 2019) as demand indicators. For cultural services, they mostly use public preference and willingness as demand indicators. Second, for ES as a whole, socioeconomic indicators (e.g., land use intensity, population density, and GDP) are usually chosen to characterize ES demand . The above methods have been frequently applied in supply -demand assessment studies of ES at global and regional scales (Torres, Tiwari, and Atkinson 2021). For ecological economics and processes-based models, complex multisource data are required from biophysical and socioeconomic fields. Many studies related to supply and demand of ES have been conducted, including assessing supply and demand ratios of ES for water supply in Mediterranean basins (Boithias et al. 2014) and spatial patterns of supply and demand of ES in the Jing-Jin-Ji region in China (Chen et al. 2022). Most of these studies used methods such as supply and demand ratios of ES and coupled coordination models to determine supply and demand imbalances of ES (Cai et al. 2021;Tao et al. 2022). Previous studies tend to aggregate to a single ecosystem type at the administrative unit or watershed scale, with less attention to vulnerable areas (Cao et al. 2021;Wei et al. 2021). Usually, these studies spatially cover one or many typical ES indicators, and an integrated supply -demand relationship assessment of ES is lacking (Deng, Li, and Gibson 2016). The coupling coordination model is still inadequate as it only involves ecological variables without exploring the coupling between land use/cover changes (LUCC) and supply and demand of ES. The matrix approach is becoming a spatially explicit tool to evaluate supplydemand changes of ES by integrating LUCC data (Burkhard et al. 2009;Ramel et al. 2020), which has obvious advantages and high efficiency in regional scenarios (Lee et al. 2016;Luiza Petroni, Siqueira-Gay, and Lucia Casteli Figueiredo Gallardo 2022). The matrix method considers the changes in supply and demand of ES caused by the dynamics of LUCC. Datasets of LUCC are used as measurement proxies to facilitate the estimation process of ES and their changes in relation to changes in LUCC, with obvious advantages and high efficiency in regional scenarios (Arowolo et al. 2018). Preliminary applications of the matrix assessment focused on the identification of spatiotemporal dynamics of supply -demand of ES under land use scenarios (Sun et al. 2019). The Gini coefficient is an emerging measure of the degree of equilibrium in the distribution of evaluation indicators (Cai et al. 2017;Benra and Nahuelhual 2019), such as inequalities between urban green space area and ES (Chen et al. 2022) and inequalities between ES and well-being . Overall, a comprehensive assessment of supply -demand relationship of ES is still scarce due to the limitation of its data availability and suitable methods (Sun et al. 2020).
Sustainable ES management is an emerging issue (Yuan and Lo 2020). The Intergovernmental Science -Policy Platform on Biodiversity and Ecosystem Services focuses on the sustainable assessment of ES and biodiversity at the global scale. The United Nations has proposed 17 Sustainable Development Goals (SDGs) with 169 nested targets as an inclusive solution ). Among them, SDG 13 (climate action), SDG 14 (life underwater), and SDG 15 (life on land) emphasize improving the value of ES and livelihoods to increase the resilience of socio-ecological systems (Lee et al. 2016). The abovementioned studies show that changes in the supply -demand relationships of ES due to anthropogenic disturbances could have cascading effects and affect the implementation of SDGs and human well-being (Linders et al. 2020). Furthermore, LUCC processes under anthropogenic stress can also affect the implementation of SDGs by changing the supply -demand relationship of ES. Thus, a quantitative assessment of the association between ES and SDGs is imperative . Therefore, the analysis of the LUCC -ES-SDGs nexus can lead decision-makers to develop rational ecological restoration projects and socio-ecological regulation schemes, and thus it is an emerging hotspot in sustainable socio-ecological systems management (Peng et al. 2021). However, the level of ES changes required to cause impacts that could be used to monitor SDGs and reflect LUCC processes remains unclear .
The Chinese Loess Plateau (LP) is an ecologically fragile area with various ecological goods and services that are spatially imbalanced. For example, the regional supply of clean water cannot meet the demand due to shortages (Feng et al. 2020), and there is an imbalance in the supply -demand of carbon sequestration in forests and grasslands . China has implemented a series of socio-ecological system management measures in the LP, such as the Grain for Green (GfG) program (Fu et al. 2011 and the Loess Plateau Comprehensive Management Outline (2010-2030) in 2011. The "greenness" of the LP is recovering due to these efforts ); however, rapid vegetation expansion may lead to unsustainable water demand by the ecosystem (Wu, Liu et al. 2019) and also affect related ES (e.g., food security) (Chen et al. 2015). To address the scientific issues and policy needs mentioned above, we (1) analyzed the LUCC spatiotemporal dynamics during 1992-2015 using the spatial mapping method; (2) assessed the multi-scale spatiotemporal variation characteristics, matching patterns, and equilibrium of supply and demand of ES integrating LUCC, matrix -index assessment, and Gini methods; (3) identified the potential impact of supply and demand of ES driven by the LUCC processes on the SDGs using a semiquantitative approach. The comprehensive assessment of the nexus between LUCC, ES, and SDGs under the supply -demand dynamics of ES will provide crucial scientific support for the management of socioecological systems and optimization of decisionmaking in the LP.

Study area
The LP (33.68° N-41.82° N, 100.85° E-114.55° E) accounts for 6.76% (6.4 × 10 5 km 2 ) of the total area of China ( Figure 1). It covers seven provinces/autonomous regions (including 341 counties). The elevation in this region is 75-5,149 m, the annual mean temperature is 3.6-14.3°C, the annual precipitation is 150-700 mm, and the annual evaporation is 1,400-2,000 mm . From 1999 to 2015, 88.20% of the vegetation in the LP was significantly restored, and the ecological management was effective (Lü et al. 2021). Socioeconomically, this region is a crucial part of the national policy for the large-scale development of western China (LSDWC). In 2015, the gross domestic product reached 5.38 trillion yuan, and the total and urban populations in the study area reached approximately 114 and 62.06 million, respectively. Regional landscape alterations have led to complex spatiotemporal heterogeneity in the supply and demand of ES (Yurui et al. 2021). For example, the increase in grain production in some regions has led to water shortages for supporting human activities, whereas 21.3% of counties in the LP showed a severe decrease in grain output (Feng et al. 2016;Wu, Liu et al. 2019), which seriously affects the implementation of SDGs in the region.

LUCC analysis
The spatial change of landscape in the LP was analyzed using LUCC maps following the kernel density method using the following equation: The kernel density function k() is a non-parametric spatial analysis method to calculate element densities in their neighborhood. The closer the value is to 1, the greater the amount of change in the element area . We converted the yearly changes in each land type into point data and calculated the kernel density values k() using ArcGIS 10.5, h is the bandwidth, x-x i is the distance from the estimated point x to the sample x i , and n is the number of samples.

Supply -demand matrices of ES
The ES indicators were collected from different landscape patterns based on the environmental and socioeconomic characteristics of the LP (Liu, Liang, and Hashimoto 2020). Following the high salience principle of supply and demand of ES, we selected services of high regional interest, including eight provisioning services, six regulating services, four cultural services, and three supporting services. We used the modified LUCC-based score matrix, which follows a three-step assessment procedure to assess the potential supplydemand relationships in ES (Burkhard et al. 2009;Ramel et al. 2020). First, the 44 land types in Burkhard's matrices were merged into six based on similar land classes. Second, we searched in Web of Science™ for literature on Burkhard's matrix to calibrate the supply and demand matrix based on the matrix values of case studies (Tao et al. 2018;Wu, Liu et al. 2019;Sun et al. 2020;Li et al. 2022). Third, the final scores were discussed and determined by four professors who worked together. Therefore, a threshold scoring (0-2) was used to develop the supply -demand matrix of ES integrating 21 ES indicators with six landscape types ( Figure 2). Specifically, the supply matrix reflected the capacity of the landscape to supply ecological products or services and divided it into low (<1), medium (1-2), and high ( = 2) supply. The supply capacity of forests and waterbodies in provisioning services was 1.25, and for others it was ≤1. Forests, grasslands, and waterbodies had higher supply capacities in regulating services compared to other land types. The supply capacities of major land cover types were >1 in cultural and supporting services (Figure 2a). The demand matrix assessed the level of human demand for ecological products or services from a specific landscape, including low (<1), medium (1-2), and high ( = 2) demand. The demand level was higher in the human-dominated built-up landscape than in the other landscapes. Natural ecosystem types have low demand levels because of their low population distribution and low product consumption. Cultivated land has high demand values (e.g., pollination and photosynthesis) in regulating and supporting services ( Figure 2b). Second, an integrated index I AESC was used to indicate the supply -demand levels of the LUCC system at a specific spatial scale (Vihervaara et al. 2010;Liu, Liang, and Hashimoto 2020).
ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi where y j is the average supply capacity or demand level of ES, A j and A are the area of LUCC type j and the total area of this region (unit: km 2 ), respectively. After using Equation 2, all the values of I AESC were normalized using Z-scores using Equations 3-5 (Cao et al. 2021). Z is the normalized supply -demand levels of ES at a specific administrative unit, μ f is the supply values or demand values of the f-th unit, � μ and S are the mean and standard deviation of μ f in all the units, respectively. To break down the characteristics of the changes in the overall supply -demand patterns of ES in the LP and to identify the mismatch between the supply and demand of ES at spatiotemporal scales, the supply (x-axis) and demand (y-axis) values were plotted to generate four quadrant diagrams and obtain high supply -high demand (H -H), low supplyhigh demand (L -H), low supply -low demand (L -L), and high supply -low demand (H -L) patterns.

Equilibrium assessment of ecosystem services
An integrated Gini coefficient, which allows for more targeted ecological management, was used to reflect the equilibrium of supply and demand of ES in the LP (Wu, Liu et al. 2020). First, the normalized values of I AESC were arranged in an ascending order, and the Lorenz curve was plotted. A single Gini coefficient was derived using the cumulative percentage of I AESC from demand (x-axis) and supply (y-axis) in a specific ES type. Finally, the Gini coefficients were aggregated to obtain the integrated Gini coefficient, and its value reflected the equilibrium of supply -demand distribution of ES in the LP (Cai et al. 2017).
where B 1 is the area between the 45° diagonal and the Lorenz curve, and B 2 is the area under the curve. Based on previous studies (Teng et al. 2011;Wu, Liu et al. 2020), the types of Gini values are divided as: high equilibrium (<0.2), medium equilibrium (0.2-0.3), basic equilibrium (0.3-0.4), disequilibrium (0.4-0.6), and extreme disequilibrium (>0.6) of supply -demand distribution of ES.
In addition, LUCC -ES-SDGs nexus was analyzed using frequency statistics (Geijzendorffer et al. 2017;Wood et al. 2018), including the influencing intensity of LUCC -ES and the frequency of ES mentioned in the 2030 Agenda for Sustainable Development policy document (United Nations 2015). Thus, a multi-step and integrated assessment framework was developed to analyze the supply -demand equilibrium of ES at the county and regional scales in the LP during 1992-2015 ( Figure 3).

LUCC analysis
Area analysis of land-use showed that the major landscapes in LP in 2015 were grassland (43.21%), cultivated land (41.8%), and forests (10.74%) (Figure 1). Except for 2003, the area of cultivated and bare land from 1992 to 2015 showed a decrease of 2.45 and 2.3%, respectively. The area of built-up land, grasslands, and forests increased by 1.95, 1.74, and 1.01%, respectively, whereas the area of waterbodies increased slightly (0.05%). Affected by the national ecological restoration policies (e.g., GfG), the decreased cultivated land area was transferred to new grassland and forests. The policies (e.g., LSDWC) also drove the rapid economic development and urbanization, resulting in a significantly increased area of built-up land (2.05% in 2015).
The values of kernel density (>0.5) showed that the major changing types are built-up land, grassland, and cultivated land during the period. Cultivated land decreased significantly in the central LP, for example, a decrease of 496.83 km 2 in Yan'an city (especially in the Wuqi, Zhidan, and Ansai counties). Increase in grassland was mainly observed at the borders of Shaanxi, Shanxi, and Inner Mongolia provinces (e.g., a 1792.92 km 2 increase in Yulin City). In Jinzhong and Changzhi cities, forests had the most significant increase of 1164.94 and 1417.14 km 2 , respectively. The expansion of built-up land was primarily centered around Xi'an, Yuncheng, Taiyuan, and Luoyang cities, where the area of built-up land increased by 1090.23, 893.79, 460.49, and 670.97 km 2 , respectively. The most significant reduction in bare land area (3644.34 km 2 ) was observed in Baotou city, as well as in the other regions of the Inner Mongolia province (Figure 4).

Supply -demand patterns of ecosystem services
The supply -demand patterns of ES were divided into four types using the quadrant diagrams ( Figure 5). The overall trend in the pattern of matching supply and demand for ES at the county scale was the same from The spatial ES distribution showed that the supply capacity was generally low in the southwestern LP, and  the demand level was high in the central region. A coexistent pattern with the H -L type was detected in the northeastern regions of the LP in 1992. After 1992, the distribution of supply -demand showed a similar pattern, but with significant differences in the proportions ( Figure 6). For example, the H -L pattern accounted for 58.24% in 2015, mainly located in the forests and grassland landscapes, indicating that the ecological environment had improved significantly through long-term water -soil conservation and ecorestoration. The region proposes to adopt an ecological priority strategy to facilitate the movement of typical ES  across the region while safeguarding the stability of the socioeconomic system. The L -H pattern was mainly concentrated in the provincial capital cities (e.g., Xi'an, Zhengzhou) and their surrounding areas, with high population density, human-dominated land use, and excessive ES consumption. It is necessary to improve the land-use efficiency and vegetation coverage in these area to mitigate the potential conflicts between low supply and high demand of ES. The L -L pattern was mainly distributed in the northern and southern parts of the LP, including 13 counties (50,274 km 2 ), a dominant landscape in the desert area with weak service supply capacity and slow economic development. The area with the H -H pattern was the smallest (39,410 km 2 ), and it has typical historical and cultural cities in the LP, such as Lintong, Helan, Zichang, Fugu, Shouyang, and Yuanping.

Supply -demand equilibrium of ecosystem services
The Gini coefficient of provisioning service showed a high equilibrium (<0.2) of supply -demand in the LP from 1992 to 2000. Further, it changed into a basic equilibrium (0.3246) affected by a decrease in provisioning of crops and timber for the implementation of the GfG policy from 2000 to 2015 (Figure 7a). The Gini coefficient of regulating services showed a disequilibrium during 1992-2015, indicating a large gap between the supply and demand distribution of regulating services. Particularly, the Gini coefficient decreased from 0.5683 in 2008 to 0.5516 in 2015, showing a slight decrease in the mismatch between the supply and demand of regulating services affected by the improvement in vegetation restoration and soil conservation capacity (Figure 7b). The Gini coefficient of cultural services showed an extreme disequilibrium (>0.6) during 1992-2015. There are many historical, cultural, and natural heritage sites in the study area. However, the development of recreation and ecotourism was inadequate, resulting in a significant imbalance between the supply and demand of cultural services (Figure 7c). The Gini coefficient of supporting services increased from 0.2139 in 1992 to 0.241 in 2015, indicating a medium equilibrium of supply -demand in supporting services. For example, a synergistic interaction between carbon and nitrogen cycles that was affected by the restoration of forest and grassland occurred during the period (Figure 7d). The synthesized Gini values (0.3-0.4) showed that the total supply -demand of ES in the LP exhibited equilibrium during 1992-2015. However, the trend toward a balanced supply and demand for synthesized ES has weakened (Figure 7e). Demand control and targeted ecological management program enhancements must be implemented to ensure ecological integrity and consistency.

Assessment of potential nexus of LUCC -ES-SDGs
The term ecosystem services has only been mentioned explicitly in SDG 15.1 (sustainable use of terrestrial ecosystem and services) in the document Transforming our World: The 2030 Agenda for Sustainable Development (https://sdgs.un.org/ 2030agenda), but all the SDG categories have described ES. The frequencies of four ES types in the description of 17 SDGs was provisioning (40), regulating (54), cultural (21), and supporting services (18) (Figure 8). Among these SDGs, the most legible description for the nexus between LUCC, ES, and SDGs are Goal 15 (life on land), Goal 11 (sustainable cities and communities), and Goal 2 (zero hunger), which were mentioned 58, 18, and 17 times, respectively. For example, the LP region is a vital crop producing area in China, where long-term severe land degradation has caused a decline in soil fertility. The LUCC management of suitable cultivated land can improve the ecological environment of local crop production and create a sustainable food production system (SDG 2) (Figure 9). In recent years, China has promoted water-efficient agricultural techniques for dryland farming, watercourse management projects in the valley plains, and the construction of checkdams in the LP, which have improved agricultural industrial structure, water and soil conservation, and farmers' income simultaneously, as well as promoted SDG 2 and SDG 6. In addition, the typical land use forms such as the restoration of vegetation and optimization of livestock structure (e.g., adjusting the grass planting area and developing forage farming and grass-livestock industry) resulted in the rapid restoration of degraded land, which promoted the achievement of SDG 15. Overall, the restoration of degraded land improved the supply of ES, such as soil conservation, flood protection, nutrient cycling, and climate regulation. The diverse LUCC processes described above contributed positively to the supply -demand equilibrium of ES and they can also promote the achievement of key SDGs (e.g., SDG 15) ( Figure 9). (P-provisioning services, R-regulating services, C-cultural services, and S-supporting services).

Applicability of LUCC matrices in the Loess Plateau
The integrated matrix method used in this study can assess spatial changes in the supply and demand of ES. The advantage of this method is that the expert scores can be corrected quickly to obtain a higher accuracy based on local ecological features. The coefficient of sensitivity (CS) was calculated to validate the reasonableness of the assessment method and the validity of the results (Wu, Liu et al. 2019). If CS < 1, ES are insensitive to the quantified value, indicating that the quantified method is reasonable. From the calculations, we found that the values for cultivated land and grassland were slightly greater than 1 while the CS values for all other land use types were <1, indicating that the results of the matrix model based on land use types were generally plausible (Table 1). However, accurate measurement of the realistic supply situation and demand level of ES still needs an in-depth exploration. It is essential to pay a great deal of attention to the connotation of supply and demand of ES to ensure the consistency of analysis and accuracies of quantitative methods.

Link age of LUCC -ES-SDGs under ecological restoration
The implementation of ecological restoration projects has led to significant landscape alterations of the LP. For example, we found that 12.57 million km 2 of cultivated land was converted to other land types during the research period, which is consistent with the previous LUCC studies in the LP region (Zhang et al. 2020;Liang, Hashimoto, and Liu 2021). LUCC affects ES mainly through different vegetation cover types. Factors such as topography and meteorology depend largely on physical geographic conditions, which do not change quickly . Therefore, spatial differences in the supply and demand of ES in different areas were closely related to the local spatial landscape patterns and socioeconomic conditions. Particularly, counties containing higher proportions of forests and grassland tended to be stronger in the provision of ES, while those with higher built-up land had lower supply capacities of ES or deficit in terms of supply -demand. The assessment also indicated that the enhancement of ecological projects and improvement of ecosystem management would alleviate the supply -demand imbalance of ES in the LP, which is generally consistent with conclusions in earlier studies (Wu, Liu et al. 2019;Jiang et al. 2021). In addition, a series of large-scale ecological measures since 1999 have led to profound changes in the ecological environment of the LP. For example, the vegetation cover in the LP increased by 17.06% and soil conservation efficacy increased by an average of 22% from 2000 to  2019 (Chinese Academy of Sciences 2020), and the desertification process in the northern LP began to reverse during the same period. These land use patterns have significantly improved the supply capacity of ES and promoted the achievement of local SDG15. Supply -demand relationships in provisioning services were balanced in the study area; for example, food services can meet the local market demand in accordance with SDG2 (zero hunger) that focuses on sustainable food production systems (Chinese Academy of Sciences 2020). The supply process of ES is not immune to spillover; the shift in the center of gravity of food was toward the northwest in China in recent years (Zhong, Liu, and Wang 2022). The cultivated land area was reduced, regional vegetation recovery was accompanied by soil drying (Lü et al. 2021), and agricultural cultivation is highly dependent on regulating services (e.g., climate regulation and pollination). The imbalance between the supply and demand of LP regulation services is not conducive to SDG realization. Natural and resource-rich habitats should be created or protected by developing multifunctional forested landscapes combined with diverse crop rotations. Promotion of multifunctional sustainable crop production is required to ensure the sustainable development of the LP.

Implications of the assessment of sustainable development in the Loess Plateau
The LP occupies an important position in ensuring the overall ecological security and resource supply of China. We found that the high supply area of LP, which accounts for more than half of its total area, mainly belongs to the economically underdeveloped regions, and the urbanization level of this region lags behind in terms of ecological and environmental improvement. Therefore, achieving a win -win situation of socioeconomic development and environmental protection is still a key challenge for this region. In addition, the ecological barrier of LP, which mainly consists of forest, grassland, and cultivated land, is mainly distributed in the Guanzhong Plain and the eastern part of the Yellow River Basin (Zhang et al., 2020), and it plays an important role in maintaining regional ecological security. However, most areas in this region have low supply and high demand pattern of ES, and the growth rate of construction land is more significant than that of forest and grassland. Demand and consumption by human populations in these areas are increasing, and therefore, these are key areas in terms of ecological protection along with the strengthening of the supply and demand flow of ES with neighboring areas. In conclusion, good land use management is needed to enhance the synergy between multiple ES (Wu, Wang et al. 2019), balance between supply and demand of ES, and realization of SDGs.

Research limitations and prospects
The CCI LUCC dataset is the first global-scale land cover map with a high spatiotemporal resolution, and our study showed that this data could reveal changing trends of land use types at a regional scale over a long time series, which is consistent with the study by Yang et al. (Yang et al. 2017). However, the spatial resolution of the CCI product (300 m) does not imply that all image cells are pure, and there is uncertainty in identifying features such as small and non-contiguous cultivated land and narrow rivers in the study area. Improving the LUCC data resolution to solve the mixed image element problem is necessary in the future. Although the 2030 Agenda for Sustainable Development provides relevant suggestions for eliminating social inequalities and protecting the ecosystem, it neither explains how to use natural resources equitably nor establishes indicators to measure progress on achieving ecological goals. Currently, provisioning and regulating services are the most studied, and cultural and supporting services are much less represented; however, in spite of this lacuna, more importance is provided for global sustainability policies. Therefore, research on theories and assessment methods must be increased for cultural and supporting services, particularly on methodology-related studies that can develop indicators.
In this study, the integrated assessment framework was used to assess the mismatch and equilibrium phenomena of the supply and demand of ES at multiple scales, and the results are useful to improve the understanding of socio-ecological conflict mechanisms. It is recommended that results of the assessment of supply and demand of ES should be integrated into landscape management and environmental decision-making to reduce the loss of socio-ecological capital in the LP. Data monitoring of socio-ecological systems and nexus assessment of LUCC -ES-SDGs should be the key domains in the future.

Conclusions
Our study assessed ecosystem services and regional implementation of SDGs using an integrated assessment framework. We identified spatiotemporal changes in landscape patterns, matrices, and composite indices of supply and demand of ES in the LP from 1992 to 2015 and analyzed the patterns of change and equilibrium characteristics of supply -demand patterns of ES at the county and regional scales. The potential impacts of LUCC -ES processes on SDGs were also analyzed. We found that land use strategies such as intensive land use and ecological restoration have driven landscape alterations in the LP. The reduction in cultivated land and bare land resulted in an increase in ecological land, indicating that the regional ecological environment has been restored. Although local provisioning and supporting services were generally in a more balanced state, there was a severe imbalance between supply and demand for regulating and cultural services. In addition, corresponding land use strategies and socioeconomic and ecological management measures should be rationally developed and implemented for different supply -demand matching patterns of ES to mitigate environmental deficit and promote sustainable socio-ecological system development. Finally, the LUCC -ES-SDGs nexus is complex, and the SDGs related to ES mainly emphasize ecosystem conservation, sustainable agriculture, and urban development. It is suggested that information on LUCC -ES should be integrated into the monitoring and evaluation programs of SDGs to provide scientific support for decision-making of sustainable development in the LP.