Equilibrium Relationship between Ecosystem Service Supply and Consumption Driven by Economic Development and Ecological Restoration

: Economic development and ecological restoration guided by policies has had a profound impact on the relationship between ecosystem service supply and consumption in the past few decades. Exploring the relationship in speciﬁc study areas is necessary to support regional planning for sustainability. In this study, we analyzed the equilibrium relationship between ecosystem service supply and consumption driven by economic development and ecological restoration from 1995 to 2018 in Northeast China. We quantiﬁed the supply and consumption of three critical ecosystem services, food production, water yield and atmospheric regulation, and identiﬁed the main areas for supply and consumption of ecosystem services. We measured the quantity and spatial equilibrium relationship between ecosystem service supply and consumption driven by economic development and ecological restoration at economic megaregion scale. We found the new industrial structure shaped by regional economic development policies formed the new ecosystem service consumer demand. The development of agricultural industry not only resulted in the increase of water consumption, but also the scope of water consumption was expanded in Northeast China. The development of tertiary industry causes the change of atmospheric regulation consumption in urban agglomeration. What’s more, economic development exacerbated imbalances in the supply and consumption of ecosystem services. The unbalanced areas where food supply was less than consumption were in the clusters of Harbin, Changchun, Shenyang, and coastal cities. Since food products can circulate freely through the market, this imbalance will not fundamentally affect regional food security, but will improve the level of agricultural industrialization and scale. The implementation of the grain economic policy caused the rapid expansion of unbalanced areas where water supply was less than consumption. The incredible space mismatch between the supply and consumption of ecosystem services still exists for a long time. However, the implementation of the ecological restoration didn’t effectively solve the problem of insatiable ecosystem service consumption. The ecosystem service consumption brought by economic development need to be given the utmost attention when formulating ecological restoration policies and economic development policies. The results can provide insights for enhancing economic sustainability in northeast china, as well other economic megaregion with similar characteristics throughout the world.


Introduction
Ecosystems are natural resources that contribute to societal wellbeing in multiple ways [1]. The past few decades have witnessed the implementation of China's regional economic development policy, with the GDP increasing from US$97.988 billion in 1990 to US$888.74 billion in 2018 in Northeast China [2]. China has successively implemented regional economic development strategies such as "The rise of the eastern coast", "The development of the western region in China", "The strategy of revitalizing the old industrial base in Northeast China" and "The rise of Central China". Regional economic development strategies make a significant contribution to economic and social development, while posing a grave pressure to the stable operation of ecosystems. During the process of regional economic development, forest, grassland, wetland and other ecosystems have been converted into human-dominated ecosystems at an unprecedented scale and pace [3]. Because of differentiated ecological service output of different ecosystems and differentiated consumption of ecosystem services by different industries, this process inevitably changes the structure of the antecedent ecosystems and regional industry [4], resulting in the large-scale changes in the supply and consumption structure of ecosystem services. At the same time, since the 1990s, China has successively implemented ecological restoration such as the first and second phases of "The Natural Forest Conservation Program (NFCP)", "3-North Shelter Forest Program", grassland protection, wetland protection, and construction of a number of ecological function zones [5]. Economic sustainability has become one of the most pressing issues facing humanity today. Understanding the change of the relationship between the supply and consumption of ecosystem services driven by economic development and ecological restoration is crucial for achieving the balance between economic development and ecological protection [6].
Considerable research has demonstrated that the economic development is the main driving force of change in the supply and consumption structure of ecosystem services (ESs) [7]. ES consumption is the sum of all ecosystem services currently consumed by society in a given period of time. However, the effects of economic development and ecological restoration on the supply and consumption of ecosystem services are not consistent among different studies, and in some cases, research results are contradictory. For example, a number of studies have revealed the supply of ecosystem services were negatively correlated with economic development [8]. The main reason is that the implementation of regional economic development strategies is often accompanied by large-scale land use changes, which destroy the original ecosystem's circular function and weaken the supply capacity of ecosystem services. However, some studies have found the opposite results-ecosystem service supply would increase during economic development [9]. The main reason is that economic development has brought about the agglomeration and flow of population, which is conducive to the restoration of the ecosystem on a large area of idle land and enhances the supply capacity of ecosystem services. To sum up, the relationship between economic development and the supply and consumption of ecosystem services is complex. Extending the obtained results directly to other study area might lead to erroneous conclusions [10]. It is necessary to detect the relationship between economic development and the structure evolution of ecosystem service supply and consumption in specific research area in order to provide scientific basis for making reasonable regional economic development policies and ecological restoration policies.
Additionally, ecological restoration policies are also an important driving force for changing the ecosystem service supply structure. However, the impact of ecological restoration on the supply of ecosystem services is not consistent in different regions. Many studies believe that ecological restoration have improved the stability of ecosystems such as forests, grasslands, and wetlands, and curbed the deterioration of the ecosystem itself [11], along with improving the supply capacity of ecosystem services. However, other studies have found that the dilemma of insufficient ecosystem services in economic and industrial functional zones has not been improved, although forests, grasslands, wetlands and other ecosystems themselves have been restored [12]. Although a series of ecological restoration policies have been implemented in Northeast China, many industrial functional areas are still facing problems such as water shortages and environmental pollution.
Studying ecosystem services from the point of view of the ecosystem itself or the economic zone itself cannot effectively solve the current ecological problems such as water shortages and environmental pollution. Therefore, under the framework of a megaregion [13], a full understanding of the relationship between the supply and consumption of ecosystem services driven by economic development and ecological restoration, and the relationship between the ecosystem service supply area and the industrial function area is of great significance to minimize the negative impact of regional economic development on the ecosystem and to enhance the effect of ecological restoration policies. The regional economic development policy has caused serious degradation of the regional ecosystem, which not only poses a threat to regional ecological security, but also restricts future economic development, although promoting the rapid economic development [14]. In addition, the changes in the industrial structure caused by economic development policies have created new demand for ecosystem services, which are growing rapidly. A large number of studies mainly focus on the supply and value of ecosystem services [1,15]. However, few articles have studied the changes in the relationship between ecosystem service supply and consumption driven by economic development and ecological restoration [16]. Without an understanding of the linkages between the major consumer and supply areas of ecosystem services, the economic and industrial overdevelopment may lead to the degradation of the supply capacity of key ecosystem services in the major supply areas. In addition, the large amount of ecosystem service consumption formed by the excessive development of industrial function zones lacks the support of the corresponding ecological function zones, which will still cause an imbalance in the supply and consumption of ecosystem services and restrict the sustainable development of the economy.
Based on the above review of the literature, this study was designed to explore the equilibrium relationship between ecosystem service supply and consumption driven by economic development and ecological restoration. The study had two components. First, three types of ecosystem services, including water, air conditioning and food production, were selected to quantitatively assess the evolution of the spatio-temporal pattern of ecosystem service supply and consumption from 1995 to 2018. Second, we determine the main supply and consumption areas of ecosystem services and quantify the equilibrium relationship between ecosystem service supply and consumption driven by economic development and ecological restoration in economic megaregion. Finally, the impact of economic development and ecological restoration on the equilibrium relationship between ecosystem service supply and consumption was analyzed. Our results are expected to provide quantitative information for realizing coordinated development of economy and environment and enhancing economic region sustainability in northeastern China and other economic megaregion with similar characteristics throughout the world.

Data Sources
Spatial data and statistical data were used in this study to analyze the relationship between regional ecosystem service supply and consumption under the dual factor of economic development and ecological restoration. The Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) with spatial resolution of 1 km and Normalized Difference Vegetation Index (NDVI) was obtained from the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www. gscloud.cn). The radiation, precipitation, and other meteorological data were obtained from the Chinese meteorological data network (http://data.cma.cn/). The spatial meteorological data were interpolated using the kriging method. The soil data were obtained from harmonized world soil database constructed by the Food and Agriculture Organization of the United Nations and the Vienna International Institute for Applied Systems (http: //webarchive.iiasa.ac.at/Research/LUC/External-World-soil-database). The population, GDP and Land use data at a spatial resolution of 1 km were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) (http: //www.resdc.cn). The statistical data, such as food in output and food consumption per capita, obtained from China Statistical Yearbook and China Household Survey Yearbook. Water consumption in industry, agriculture and daily life is derived from China Water Resources Bulletin and China Statistical Yearbook. Air monitoring data obtained Bulletin on the Ecological Environment Ecological Environment Status Bulletin of Inner Mongolia Autonomous Region, Heilongjiang Province, Jilin Province and Liaoning Province. Air monitoring data were obtained from the PM 2.5 Historical data Center (https://www. aqistudy.cn) and the Ecological environment status bulletin of Inner Mongolia Autonomous Region, Heilongjiang Province, Jilin Province and Liaoning Province. The spatial air monitoring data were interpolated using the inverse distance weighted method. A small amount of missing data was replaced by data for neighboring years.

Food Supply (FS) and Food Consumption (FC)
Food production is an important part of regional economic and industrial development policies. What's more, food supply is crucial for food security and regional sustainable development [17]. A large number of studies have shown that there is a significant linear relationship between food yield, such as crops and livestock, and NDVI. According to the land use data, the total rice output is allocated according to the ratio of the grid NDVI value to the total NDVI value of the paddy field. The total output of foods such as corn, wheat, oilseeds, grains and vegetables are distributed according to the ratio of the grid NDVI value to the total NDVI value of dry farmland. The total production of meat and milk is allocated according to the ratio of the grid NDVI value to the total NDVI value of grassland. The total production of aquatic products is allocated according to the ratio of the grid NDVI value to the total NDVI value of the water area. The food supply capacity of each grid is then quantified. At the same time, the total food supply of northeast China is calculated by summation (Table 1). Main food consumption per capita includes the type and weight of food ( Table 2). The equation for Food supply and food consumption is as follows: where FS i is the supply of food in the i pixel, including rice, grains other than rice, meat, milk, and aquatic products, with a spatial resolution of 1 km. G sum is the total output of food (rice, grains other than rice, meat, milk, and aquatic products). NDVI i is the normalized vegetation index of i pixel. NDVI sum is the sum of NDVI values of cultivated land, grassland or water area respectively.
where FC i is the consumption of food in the i grid, with a spatial resolution of 1 km. FC pc is food consumption per capita. L pop stands for population density.

Water Supply (WS) and Water Consumption (WC)
Water production refers to the ability of ecosystems to intercept or store water resources from rainfall and to reduce surface runoff. In Northeast China, industry and agriculture consume a lot of water resources. Therefore, evaluating the supply and consumption of water production is critical to the rational use and protection of water resources. The InVEST model was adopted to quantify the supply of water production services. The equation for water supply and water consumption is as follows: where WS represents the annual water supply (mm), Aet i is the actual annual evapotranspiration of the grid unit (mm), Pet i is the potential evapotranspiration of the grid unit, Kc i is the crop evapotranspiration coefficient, Eto i is the reference evapotranspiration of vegetation, R i represents annual precipitation, Awc i is the available water content of plants, W is the empirical parameter, Z is the Zhang coefficient: (12) where WC refers to the water consumption in northeast China, WC pf is the water consumption in the paddy field, WC df is the water consumption in the dry farmland, WC il is the water consumption in industry and life, D sp is the total water consumption in the paddy field, D sd is the total water consumption in the dry farmland, and D sr is the total water consumption of industry and life, NDVI pi is the normalized leaf area index of the i-th grid in the paddy field, NDVI di is the normalized leaf area index of the i-th grid in the dry farmland, and NDVI psum is the sum of the normalized leaf area index of the paddy field, NDVI dsum is the sum of the normalized leaf area index of dry farmland, P i is the population of the i-th grid, and P sum is the total population in Northeast China. D a is the total water consumption of cultivated land. According to research on agricultural water consumption [18], the average water consumption per unit area of paddy fields is 6.135 times that of dry land in Northeast China.

Atmospheric Regulation Supply (AS) and Atmospheric Regulation Consumption (AC)
Atmospheric regulation means that the ecosystem maintains the balance of atmospheric chemical components, such as absorbing SO 2 , absorbing nitrogen oxides, dust, etc. Atmospheric regulation has been considered to be an important ecological problem in northeast China. The amount of atmospheric regulation for each pixel on the landscape was obtained by the equation described as follows: where AS i represents the annual amount of atmospheric regulation in the i-th land use type, A SO2 is the annual absorption of SO 2 , A Nx is the annual absorption of nitrogen compounds, and Apm10 is the annual absorption of inhalable particulate matter. With reference to a large number of documents on monitoring data [19,20], the absorption amount per unit area of each land use type SO 2 , nitride, pm10 was calculated by the average method ( Table 3): where AC represents the annual consumption of atmospheric regulation services per square kilometer (kg·km −2 ·a −1 ), C SO2 is the annual average SO 2 concentration (µg/m 3 ), C Nx is the annual average nitrogen oxide concentration (µg/m 3 ), C pm10 is the annual average inhalable particulate matter concentration (µg/m 3 ). According to the relevant research on the transit height of dust and haze in Northeast China [21,22], the pollutant suspension height (Ph) is set to 500. The number of days to consume atmospheric regulation services per year is set to 365. The unit conversion factor (β) is 0.001. Table 3. Amount of pollutants absorbed by each landscape type per unit area (kg·km −2 ·a −1 ). The supply and consumption ratio of ecosystem services (SCRE) was adopted to analyze the relationship between ecosystem service supply and consumption and to reveal the equilibrium level of ecosystem service supply and consumption. SCRE for each pixel was computed as follows [23]:

Value Absorption of SO 2 Absorption of NO X Absorption of PM10
where SCRE represents equilibrium relationship between ecosystem service supply and consumption, S is the supply of ecosystem services, C is consumption of ecosystem services, S max is the maximum supply of ESs within the space of the assessment area, C max is the maximum consumption of ESs within the space of the assessment area. SCRE > 0 indicates that the supply of ecosystem services is higher than consumption, SCRE = 0 indicates that the supply of ecosystem services is balanced with consumption, and SCRE < 0 indicates that the supply of ecosystem services is lower than consumption.

Spatial Equilibrium between Ecosystem Services Supply and Consumption
Local Indicators of Spatial Association (LISA) was adopted to reflect the clustering and matching relationship between ecosystem service supply and consumption, and to reveal the spatial equilibrium relationship between ecosystem service supply and consumption. The GEODA software is used to calculate LISA with the following equation: where w ij represents the spatial weight matrix between unit i and unit j, y i is the attribute value of unit i, y is the average of all attribute values, n is the total number of cells in the region. LISA > 0 indicates that the service balance of the space unit is a high-high value or low-low value spatial aggregation. LISA < 0 indicates that the service balance of the space unit is a high-low value or low-high value spatial aggregation.

Impact of Economic Development on the Equilibrium Relationship between Ecosystem Service Supply and Consumption
The impact of economic development on the equilibrium relationship between ecosystem service supply and consumption is complex, and linearity, "inverted U" shape, "S" shape and other nonlinear relationship may occur. We applied curve estimation to identify the relationship between economic development and the equilibrium relationship. Curve estimation in SPSS 19.0 for Windows (IBM, Armonk, New York, USA), was used to analyze the relationship between GDP and SCRE. The correlation between GDP and equilibrium coefficient (SCRE) of ecosystem service supply and consumption was calculated using SPSS 19.0 software from 1995 to 2018. The basic unit for curve estimation was the county. Due to the different size and unit for the indicators of SCRE and GDP, GDP were centralized before conducting curve estimation.

The Spatial Evolution of Food Supply
The spatial distribution of ecosystem service supply from 1995 to 2018 is shown in Figures 2 and 3. The map indicates a very high spatial variation of ecosystem service supply across the study area from 1995 to 2018. In 1995, the high food supply areas were mainly distributed in Jilin and Liaohe Plain, while the food supply capacity per unit area was slightly weaker in Heilongjiang and Inner Mongolia. According to Table 4, the food supply of Liaoning and Jilin was 14.60 and 20.47 million tons, which was 50.17% of the total food supply of Northeast China. The food supply was just 8.19 million tons in Inner Mongolia. The food supply of Heilongjiang was 26.65 million tons, which was just 38.12% of the total food supply of Northeast China. However, the main food supply areas were gradually shifting to Heilongjiang Province and Sanjiang Plain By 2010 in northeast china. The food supply of Heilongjiang was 54.53 million tons, which was 43.85% of the total food supply of Northeast China. The food supply of Liaoning was 19.41 million tons, which was just 15.61% of the total food supply of Northeast China in 2010. The high food supply areas were concentrated in the Songnen Plain and the Sanjiang Plain with the implementation of the grain production policy in 2018, while the food supply status of the Liaohe Plain was gradually declining. The food supply of Heilongjiang was 78.72 million tons, which was 46.43% of the total food supply of Northeast China in 2018. The food supply of Liaoning and Jilin was 14.60 and 20.47 million tons, which was just 36.32% of the total food supply of Northeast China. Table 4. Total supply and total consumption of ecosystem services in the four provinces of Inner Mongolia, Liaoning, Jilin and Heilongjiang. (FS-food supply (million tons); FC-food consumption (million tons); WS-water supply (million cubic meters); WC-water consumption (million cubic meters); AS-atmospheric regulation supply (million kilograms); AC-atmospheric regulation consumption (million kilograms)).

Spatial Evolution of Water Supply
The results obtained from Figures 2 and 3 indicated that the spatial distribution of water supply also changed greatly from 1995 to 2018. The high water supply areas were mainly concentrated in the east of Jilin Province and the southeast of Liaoning Province in 1995 and 2010. However, not only the high water supply areas were distributed in the east of Jilin Province and the southeast of Liaoning Province, but also the vast majority of Heilongjiang Province became high water supply areas with the implementation of the ecological restoration in 2018. This shows a certain trend that the high water supply area extends to the north and west. According to Table 4

Spatial Evolution of Atmospheric Regulation Supply
The results obtained from Figures 2 and 3 indicated that the spatial distribution of atmospheric regulation service supply has changed to some extent from 1995 to 2018. In general, the high supply areas of atmospheric regulation services are mainly distributed in the northeast of Inner Mongolia, the northeast and west of Heilongjiang province, the east of Jilin Province, and the east of Liaoning Province. However, the supply of atmospheric regulation service in Chifeng, Xingan, Chaoyang, Huludao, and Harbin showed a significant increase from 1995 to 2018. The atmospheric regulation supply in Inner Mongolia increased from 183,759.67 to 193,053.53 million kilograms, and the atmospheric regulation supply in Liaoning increased from 63,107.64 to 69,949.39 million kilograms, from 1995 to 2018. This showed that the supply capacity of atmospheric regulation services was improving in these areas from 1995 to 2018, with the implementation of the ecological restoration such as 3-North Shelter Forest Program, Grass protection and NFCP.

Spatial Evolution of Food Consumption
The spatial distribution of the food consumption shown in Figures 2 and 3 shows a different pattern: the spatial concentration trend of food consumption. The spatial distribution of food consumption moved from the surrounding areas to three major urban agglomerations of Harbin, Changchun and Liaoning from 1995 to 2010. With accelerating population concentration in the city, economic development has promoted the process of concentration of food consumption in major cities from 1995 to 2010. The spatial distribution of food consumption moved from the surrounding areas to three major urban agglomerations of Harbin, Changchun and Shenyang from 2010 to 2018. Benefiting from the excellent development of coastal cities, the centralized pattern of food consumption from 2010 to 2018 has shown a trend of gradual strengthening to coastal city clusters such as Shenyang and Dalian.

Spatial Evolution of Water Consumption
The spatial distribution of the water consumption changed remarkably from 1995 to 2018. The water consumption in Huludao, Chaoyang, Chifeng, Jinzhou and Fuxin showed a decreasing trend, while the water consumption shows an increasing and expanding trend in Sanjiang Plain and Songnen Plain of Heilongjiang Province (Table 4, Figures 2 and 3). According to Table 4, Heilongjiang's share of water consumption in northeast china increased from 40.68% (24,448.46 million cubic meters) to 47.63% (34,356.03 million cubic meters). According to the China Statistical Yearbook, the share of agricultural water in the total water consumption increased from 69.72% to 88.63% in Heilongjiang from 2004 to 2018. However, the share of agricultural water in the total water consumption decreased from 87.55% to 73.03% in Inner Mongolia from 2004 to 2018. The share of agricultural water in the total water consumption decreased from 65.11% to 61.78% in Liaoning from 2004 to 2018. These evolutions show that Heilongjiang Province, as the most important commercial grain base, not only increased the area of dry farmland and paddy fields, but also increased the consumption of agricultural water from 1995 to 2018 (Table 4, Figures 2 and 3). At the same time, due to the decline in food production status in southwestern Liaoning, agricultural water consumption has also experienced a major reduction. From 1995 to 2010, the maximum water consumption per square kilometer for domestic and industrial water in Shenyang, Changchun, and Harbin rose from 1.95 to 13.3 million cubic meters. Between 2010 and 2018, the maximum water consumption per square kilometer was reduced from 13,300,000 to 440,000 cubic meters (Figure 3), with the industrial structure shifted from high water-consuming industries to tertiary industry (Table 5).

Spatial Evolution of Atmospheric Regulation Consumption
The distribution of atmospheric regulation service consumption shows obvious variation in spatial structure from 1995 to 2018. The consumption areas for atmospheric regulation service were mainly distributed in the central and southwestern regions of Heilongjiang Province and the Shenyang urban agglomeration in 1995. The demand area for atmospheric regulation services showed a trend of extending to the central, southern and western parts of the study area from 1995 to 2018. According to Table 4 To some extent, this evolution was affected by the transformation of the economic structure from industrial industries to service industries. According to Table 5, the proportion of the secondary industry dropped significantly, while the proportion of the tertiary industry rose rapidly in northeast china.

Equilibrium Relationship between Food Supply and Consumption
The results obtained from Figure 4 indicate that the equilibrium relationship between food supply and consumption changed greatly from 1995 to 2018. The large-scale unbalanced areas where food supply was lower than consumption were mainly distributed in the Songnen Plain, the Liaohe Plain and the eastern part of the Inner Mongolia Autonomous Region in 1995. What's more, the most serious unbalanced areas were mainly concentrated in Harbin, Changchun, Shenyang and other areas. However, with the construction of Commodity Grain Base in northeast china in northeast China, the situation of insufficient food supply was changed in most areas of Songnen plain, Liaohe Plain and eastern Inner Mongolia in 2010. The scope of insufficient food supply was narrowed largely, mainly in Harbin urban agglomeration, Changchun urban agglomeration, Shenyang urban agglomeration, Dalian urban agglomeration, etc. from 1995 to 2010 (Figure 4). However, the minimum value of SCRE has changed from −1.6172 to −1.9130, showing a trend of increasing imbalance in food supply in urban agglomerations from 1995 to 2010. What's more, the scope of insufficient food supply expanded in southern Liaoning and coastal urban agglomerations from 2010 to 2018.

Equilibrium Relationship between Water Supply and Consumption
The equilibrium relationship between water supply and consumption shows obvious variation in spatial structure from 1995 to 2018 (Figure 4). The unbalanced areas where water supply was lower than consumption were mainly distributed in Harbin, Daqing, Suihua, Shenyang, Liaoyang, Anshan, Panjin, Yingkou in 1995. These distribution shows that the industrial structure of over-reliance on heavy industry exacerbated the imbalance between water supply and consumption. The unbalanced areas where water supply was lower than consumption were mainly distributed in western Liaoning Province, eastern Inner Mongolia Autonomous Region, Shenyang and Harbin. However, the highly unbalanced areas where water supply was higher than consumption (SCRE > 0.1002) remarkably decreased from 1995 to 2010. The unbalanced areas where water supply was lower than consumption expanded to central Liaoning and the Sanjiang Plain from 2010 to 2018. In addition, the highly unbalanced areas where water supply was higher than consumption (SCRE > 0.1002) remarkably increased in Heilongjiang Province, eastern Liaoning Province, eastern Jilin Province, and northeastern Inner Mongolia from 2010 to 2018, with the construction of ecological function zone for water conservation (Greater Khingan Mountains, Lesser Khingan Range and Changbai Mountain).

Equilibrium Relationship between Atmospheric Regulation Supply and Consumption
The unbalanced areas where atmospheric regulation supply was lower than consumption were mainly distributed in Harbin, Jiamusi, Shenyang, Liaoyang and Anshan in 1995. The unbalanced areas where atmospheric regulation supply was lower than consumption extended to the south from 1995 to 2010, mainly in Liaoning Province, Jilin Province, and Harbin. According to Table 5 (Table 4). These results show an expanding trend in the area of insufficient atmospheric regulation, with rapid economic (the secondary industry) development. What's more, the minimum value of SCRE has changed from −0.0782 to −0.0832, showing a trend of increasing imbalance in atmospheric regulation from 1995 to 2010. According to Tables 4 and 5, the unbalanced areas where the supply of air conditioning was lower than consumption had been greatly reduced from 2010 to 2018, with the proportion of the secondary industry dropped significantly. In reality, there were still large areas of potential imbalance in Heilong Province, Jilin Province and Liaoning Province in 2018.

Spatial Relationship between Food Supply and Consumption
The results of LISA obtained from Figure 5 indicates that the dominant matching type of food in Northeast China is low supply and low consumption. The space range of the two space matching types of low supply-high consumption and high supply-high consumption changed little from 1995 to 2010. The spatial range of matching types of high food supply and high consumption changed the most in 2010, which mainly extended in Heilongjiang Province and Sanjiang Plain. What's more, The spatial range of matching types of high food supply and high consumption extended in Songnen Plain and Sanjiang Plain from 2010 to 2018, with the development of the agricultural industry in Heilongjiang as the most important commodity grain base in China. However, the scope between low supply-high consumption and high supply-high consumption was narrowed from 2010 to 2018. These evolutions show that economic development have not only caused a trend of centralization of food demand, but also weakened the food supply capacity of areas with high urbanization levels. This type of specialization combined with trade between regions is an result of economic development.

Spatial Relationship between Water Supply and Consumption
The LISA results obtained from Figure 5 shows that the matching types of water supply and consumption changed significantly from 1995 to 2018. The areas with low supply and high consumption of water were mainly distributed in the south-central Songnen Plain and Liaohe Plain. The spatial scope of low supply and high consumption of water was gradually concentrated in the southern part of the study area, with a certain degree of expansion in Liaoning from 1995 to 2018. The spatial scope of high supply-high consumption of water gradually decreased from 1995 to 2018, mainly in the southeast of Liaoning Province. However, due to the implementation of the ecological restoration policy by the Chinese government since 1998, the spatial scope of high supply and low demand gradually expanded in Northeast China from 1995 to 2018, mainly in the east of Jilin Province and the northwest of Heilongjiang Province.

Spatial Relationship between Atmospheric Regulation Supply and Consumption
The LISA results obtained from Figure 5 showed that the areas with high supply and high consumption of atmospheric regulation were mainly located in the east of Jilin Province and the east of Heilongjiang Province in 1995. The areas with high supply and high consumption of atmospheric regulation gradually shifted to the south from 1995 to 2018, mainly in eastern Jilin Province, eastern Liaoning Province, Chaoyang, Chifeng, and Fuxin. The evolutions show that both the supply and consumption for atmospheric regulation were increasing in eastern Liaoning Province, Chaoyang, Chifeng, and Fuxin from 1995 to 2018. The areas with low supply and high consumption of atmospheric regulation increased substantially from 1995 to 2010, mainly in Inner Mongolia Autonomous Region, Liaoning Province and Jilin Province. The areas with low supply and high consumption of atmospheric regulation decreased substantially from 2010 to 2018, mainly in Liaoning Province, eastern Inner Mongolia Autonomous Region and southern Jilin Province. The evolutions show that the driving effect of economic development on atmospheric regulation shows an upward trend and then a downward trend. With the implementation of the ecological restoration policy, the areas with high supply and low consumption of atmospheric regulation were extending from 1995 to 2018.

Impact of Economic Development on the Equilibrium Relationship between Ecosystem Service Supply and Consumption
The results obtained from Table 6 showed that the equilibrium relationship between ecosystem service supply and consumption was negatively correlated with economy developing level. Economic development exacerbated the imbalance between food supply and consumption from 1995 to 2018. The main reason was that economic development intensified the agglomeration around the population, resulting in increased imbalances between food supply and consumption. What's more, there was a significant negative correlation between the level of economic development and the imbalanced level of water supply and consumption in 1995 and 2018 (Table 6). In addition, economic development exacerbated the imbalance between atmospheric regulation supply and consumption from 1995 to 2018 (Table 6). Table 6. Correlation between GDP and equilibrium coefficient (SCRE) of ecosystem service supply and consumption.

Ecosystem Service Supply and Consumption Evolution Driven by Economic Development and Ecological Restoration
The sustainability of economic development and the effectiveness of ecological restoration have received widespread attention in recent years [24]. Regional economic development policies not only promoted the development of agriculture, industry, and service industries in Northeast China, but also their structural relationships continued to dynamically evolve from 1995 to 2018. However, the new industrial structure formed a new ecosystem service consumer demand. The development of agricultural industry not only resulted in the increase of water consumption [25], but also the scope of water consumption was expanded in northeast China. The development of tertiary industry causes the change of atmospheric regulation consumption in urban agglomeration. Regional economic development not only promoted the concentration of population in coastal cities, Shenyang urban agglomeration, Harbin urban agglomeration, and Changchun urban agglomeration, but also resulted in the concentration of food consumption, industrial and domestic water consumption in cities. The ecological restoration policies improved the ecological services of water, food, and atmosphere regulation in northeast china in the past 25 years. Such improvement in ecosystem services was consistent with that from the China Ecosystem Assessment at the national scale [26], and from other studies [27]. However, there was a strong spatial difference in the improvement of ecological service supply.

Significance of Analyzing the Relationship between Ecosystem Service Supply and Consumption
Separating the supply of ecosystem services from consumption to conduct research cannot effectively solve the problem of insatiable ecosystem service consumption [28]. In the past 25 years, economic development strategies such as The strategy of revitalizing the old industrial base in Northeast China and The development of the western region in China were implemented, which greatly promoted the development of the regional economy. The output value of the primary, secondary, and tertiary industries greatly increased. However, the economic development had a profound impact on the relationship between the supply and consumption of ecosystem services, and there were considerable differences in multiple service types [29]. At the same time, in order to solve the ecological crisis, the government successively implemented ecological restoration policies such as The Natural Forest Conservation Program (NFCP), the construction of ecological function area, and the construction of ecological civilization since 1998 [30]. Although this restoration policy increased the value of ecosystem services in northeast china [26], the situation of water supply and air pollution in the region became more severe. The equilibrium relations and spatial matching relations were introduced into this study to examine changes in the relationship between ESs supply and consumption driven by economic development and ecological restoration. First of all, food production increased considerably in northeast China, with the development of the economy in the past 25 years. However, the expansion of urban land use types and the concentration of population caused by economic development exacerbated the imbalance of food supply. In fact, this imbalance is the result of specialization and regional division. Since food products can circulate freely through the market, this imbalance will not fundamentally affect regional food security, but will improve the level of agricultural industrialization and scale. Secondly, the economic development policy exacerbated the imbalance of water supply, while promoting agricultural and industrial development. Not only did the minimum value of the imbalance decrease from −0.9520 to −1.5524, but the range of the imbalance also expanded. The degree of spatial mismatch of water was also shocking, and it did not show a trend that the range of mismatch will shrink. Thirdly, economic development also affected the imbalance of atmospheric regulation in northeast China. The imbalance of atmospheric regulation from 1995 to 2010 became more severe in scope and extent. Although the imbalance of atmospheric regulation was moderating from 2010 to 2018, there were still quite extensive areas where supply and consumption didn't match.
The results show that the implementation of the ecological restoration didn't effectively solved the problem of unsatisfied consumption of ecosystem services due to the existence of imbalance and spatial mismatch, although improving the value of ecological services [26]. Due to the limited radiation range of water and atmosphere regulation [31], the free circulation of space in the market environment cannot be realized, so it didn't significantly improve social welfare like food. In order to improve the efficiency of ecological investment, the consumption-oriented factor should be taken into account in future ecological restoration policies. It is necessary to consider the potential ecosystem service demand types of future industry types, and establish a matching ecosystem service supply functional area [32], while formulating economic development policies.

Integrating Ecological and Economic Approaches
The new demand of ecosystem services triggered by regional economic development is not usually considered because obtaining this type of data is very time-consuming. However, in this study, we have demonstrated that not only the value of ecological services can be overestimated or underestimated [16], but there is also impossible to solve the ecological service supply and consumption gap caused by new ecological service consumption if we consider only the ecosystem service supply. What's more, the imbalance between supply and consumption weakens the supporting role of ecosystem to social and economic system. In fact, only a small percentage of the territory provided efficient supply that can meet the ecosystem service consumption. Solving the sustainability of economic development is a very important point of departure for ecosystem services research [33], as they facilitate a more inclusive ecological-economic approach by using economic methods to estimate changes in ecosystem service demand caused by economic development or changes in industrial structure, and exploring the interaction between economic and ecological processes. Incorporating the ecological service demand brought about by future industrial reform into the ecological restoration policy considerations will improve the efficiency of the ecological restoration policy and solve the sustainability problem of economic development.
This study is focused on the analysis of the relationship between ecosystem service supply and consumption, although to support decision-making it is necessary to analyses different ES because synergies and trade-offs between them may occur [34], and to analyses delivery method of ecosystem services. We considered that the trend of ecosystem service consumption brought by economic development would need to be given the utmost attention when analyzing ES for effective policy implementation and management, as the realization of the ecosystem service flow strongly depends on economic development process. Understanding which ecosystem services are more valued by economic social system is useful for targeting resources efficiently. In addition, under the framework of megaregions [12], studying the relationship between the supply and demand of ecosystem services can better understand the interaction process between ecosystems and economic systems. This type of study has important reference implications for the formulation of ecological restoration policies and economic development policies, as we identify the area with the highest or lowest demand based on changes of economic industry and show the equilibrium relationship between ecosystem service supply and consumption.

Conclusions
Our work presented an interdisciplinary framework by quantifying the spatial equilibrium relationship between the supply and consumption of ecosystem services to optimize ecosystem management, and used the framework to conduct a study on the relationship between the supply and consumption of ecosystem services driven by economic development and ecological restoration in Northeast China, which contributed to the formulation of economic development policies and ecological restoration policies.
The results show that the study of ecosystem service supply and consumption separately cannot effectively solve the problem that ecosystem service consumption cannot be satisfied. The economic development policies aggravated the spatial imbalance between the supply and consumption of ecosystem services in northeast China in the past 30 years, while promoting economic development. The imbalance between the supply and consumption of food, water and air in Northeast China became more severe. The unbalanced areas where food supply was less than consumption were in the clusters of Harbin, Changchun, Shenyang, and coastal cities, and there were a serious spatial mismatch between supply and consumption. Since food products can circulate freely through the market, this imbalance will not fundamentally affect regional food security, but will improve the level of agricultural industrialization and scale. The implementation of the grain economic policy not only caused the rapid expansion of unbalanced areas where water supply was less than consumption, but also aggravated the spatial mismatch between water supply and consumption in Northeast China. With economic development, the unbalanced situation in which the supply of atmospheric regulation was lower than the consumption became more severe (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010), and then experienced a reduction (2010-2018). However, the mismatch between atmospheric regulation supply and consumption was still relatively widespread from 2010 to 2018. On the other hand, the implementation of ecological restoration policies improved the supply of ecosystem services to varying degrees in Northeast China.
We concluded that economic development exacerbated spatial imbalances in the supply and consumption of ecosystem services over the past 30 years in Northeast China. What's more, the incredible space mismatch between the supply and consumption of ecosystem services still exists for a long time. The implementation of the ecological restoration didn't effectively satisfy the consumption of ecosystem services, although improving the quality and value of ecological services in Northeast China. Economic development policies which promote the change of industrial structure will form the new consumption demand of ecosystem services. The ESs consumption direction brought by economic development need to be given the utmost attention when formulating ecological restoration policies and economic development policies. The analysis of the relationship between the supply and consumption of ecosystem services can help understand the process of interaction between ecosystems and economic systems. Future ecological restoration policies should include the consumption of ecosystem services in the scope of decision-making to improve the efficiency of ecological restoration policies.

Data Availability Statement:
The data presented in this study are available on request from thecorresponding author. The data are not publicly available due to confidentiality assurance of eachparticipant's information.