Next Article in Journal
Livelihood Capital Effects on Famers’ Strategy Choices in Flood-Prone Areas—A Study in Rural China
Next Article in Special Issue
Construction and Optimization of an Ecological Network in Zhengzhou Metropolitan Area, China
Previous Article in Journal
COVID-19 Pandemic Psychological Impact and Volunteering Experience Perceptions of Medical Students after 2 Years
Previous Article in Special Issue
Research on Village Type Identification and Development Strategy under the Background of Rural Revitalization: A Case of Gaochun District in Nanjing, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Variation of Ecosystem Services Value under Land Use/Cover Change in the Black Soil Region of Northeastern China

1
School of Public Administration and Law, Northeast Agricultural University, Harbin 150030, China
2
Land Remote Sensing Big Data Technology Innovation Center, Harbin 150030, China
3
School of Economics and Management, Northeast Agricultural University, Harbin 150030, China
4
School of Public Administration, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(12), 7533; https://doi.org/10.3390/ijerph19127533
Submission received: 21 April 2022 / Revised: 15 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022
(This article belongs to the Special Issue Land Use and Urban-Rural Integration)

Abstract

:
A better understanding of the dynamic variation in the ecosystem service value (ESV) under land use/cover change (LUCC) is conductive to improving ecosystem services and environmental protection. The present study took Landsat TM/ETM remote sensing images and socio-economic statistic data as data sources and extracted land-use data using RS and GIS technology at 5-year intervals from 1990 to 2020. Then, we interpreted the spatio-temporal characteristics of LUCC and analyzed ESV changes using the value equivalence method in the black soil region of northeastern China (BSRNC). The main results showed that land use changed significantly during the study period. Cultivated land continued to expand, especially paddy areas, which increased by 1.72 × 106 ha, with a relative change of 60.9% over 30 years. However, grassland decreased by 2.47 × 106 ha, with a relative change of −60.6% over 30 years. The ESV showed a declining trend, which decreased by CNY 607.96 million during 1990–2020. The decline in forest and grassland caused a significant decline in the ESV. Furthermore, the ESV sensitivity coefficients were less than one for all of the different categories of ecosystem services. LUCC has a considerable impact on ESV in the BSRNC, resulting in ecosystem function degradation. As a result, future policies must emphasize the relationship between food security and environmental protection in situations of significant land-use change.

1. Introduction

Ecosystem services refer to services that affect human well-being and that are derived from the ecosystem on the Earth’s terrestrial surface [1,2]. Hence, the links between man and nature are closely related and may even be inseparable. As such, this interconnection implies repercussions for ecosystem functioning. Previous studies have highlighted that human manipulation and the degradation of the environment have reached the point in which the cumulative global impacts affect our food, fiber, and water systems and threaten people’s health and well-being [3]. However, the population growth issues that are connected with socioeconomic development may reflect the contradiction between human activities and human well-being [4]. Under this context, ecosystem service value (ESV) is indispensable when balancing the need for food production and ecological protection because it can clarify the ecological assets and values provided by a region [5].
Additionally, changes in the functioning of the Earth’s ecological systems can pose a significant threat to human survival. In this situation, implementing suitable land and ecosystem services are essential carriers of terrestrial ecosystems. Land-use change also affects regional ecosystem service functions by changing the distribution of regional resources and ecosystem types. The landscapes of ecological systems have been considerably affected by ongoing land use/cover change (LUCC), which represents the biggest element of global change. In addition to the background noted earlier, another study shed light on that fact that since the mid-20th century, 60% of global ecosystem services have been degraded due to irrational land utilization [6]. This mode of land use has appeared in different ways in different regions. This scenario has resulted in desertification, deforestation, climatic change, and decreased carbon sequestration [7,8]. This situation is visible in large agricultural countries [9,10], especially China. Within this region, the ecosystem service value (ESV) is a significant indicator of the ecosystem function changes caused by LUCC [11,12,13]. Therefore, analyzing the dynamic variation in ESV under LUCC has become a widespread issue. A detailed comprehension of the variability in the ESV may potentially improve the performance of global environmental policies.
ESV approaches have been used to evaluate ecosystem services on a global and regional scale [14,15,16]. In that sense, ESV assessment may provide a momentous reference for ecosystem protection [17]. Costanza et al. [14] first proposed analyzing ESV changes on a global scale. Due to the tempo-spatial heterogeneity, this study cannot be directly applied to the national scale, especially in China. Xie et al. [16] revised the equivalent factor to establish the national scale spatio-temporal dynamic assessment method, resulting in different regional conditions showing unique ecosystem services and functions [18,19]. Nonetheless, this method has been improved. The evaluation method and ESV equivalence table have been further revised based on existing studies and factors [20].
Many studies have revealed that land-use dynamics can cause ESV change [9,21]. LUCC may increase ecosystem supply or cause ecological degradation [22]. For instance, population growth and food demands have driven the transfer of cultivated land to grassland [23,24]. The rapid shrinkage of forest and grassland areas has led to the degradation of ecosystem functions. In this situation, it is essential to analyze the relationship between ecosystem functions and land use to better understand how LUCC affects ecosystem functions. However, existing land and ecological policies ignore assessing and reflecting the hidden ecological value and cost of land use and rarely consider the relationship between land use and ecosystem services [25,26].
In China, the government has launched ecological restoration and protection policies to improve ecosystem function [27,28,29]. Previous research has shown that establishing ecological restoration and conservation policies enhances ecosystem services and functions [30,31]. Some studies have also neglected to cover the loss in the ESV caused by the lack of land-use policies implemented at the regional level [32]. Hence, the sustainable development of ecological functions and systems requires a comprehensive analysis of the variations in LUCC and ESV on a large spatial scale. With the development of remote sensing and GIS technologies, the quantitative analysis of LUCC and ESV has become more accessible [33]. In this work, we estimated the ESV and evaluated the relationship between LUCC and the ESV, providing data reference that can be used by regional governments to coordinate the interaction between land-use policies and ecosystem services. Some researchers have shown that the ESV often varies in terms of its temporal and spatial characteristics [34,35]. For example, Shang et al. found that the interaction between the NDVI and soil types had the most substantial impact on the spatial differentiation of the IESV in the study area, with an explanatory power of 21.95% [5]. Therefore, it is essential to carry out some studies in different regions to analyze changes in the ESV.
The black soil region of northeastern China (BSRNC) is one of the four black soil regions in the world. This area represents the area with the largest amount of grain production in China, accounting for 1/4 of the country’s total grain output [36,37]. The BSRNC has evolved from natural forest and grassland ecosystems to artificial farmland ecosystems [38]. About 70% of the land categories have shown changes that have resulted in deforestation, soil, and water loss [39]. Irrational land use has led to decreases in the fertility of black land and the degradation of ecosystem services. Primarily, cultivated land expansion has placed tremendous pressure on the ecosystem. The government has established ecological restoration and protection measures to enhance ecosystem functions, such as supporting the sustainable development of prominent grain-producing regions and land conservation tillage planning [40]. Although these measures have been implemented to restore and protect ecosystems, the ESV has rarely been assessed. Furthermore, the variation in the ESV under LUCC is not clear in this vital region. To some extent, the lack of data and information on the spatio-temporal changes resulting from LUCC hinders black land protection and ecological restoration processes. Therefore, an accurate study outlining the rules of ESV variation under LUCC plays a strategic role in promoting environmental conservation and natural resource utilization in the BSRNC.
This paper takes the whole BSRNC as its study area to research ESV changes caused by LUCC from 1990 to 2020. The objectives of this study are as follows: (1) analyze the temporal and spatial variation caused by LUCC in the BSRNC; (2) revise the equivalence factor to calculate the ESV; And (3) evaluate the credibility of the research results. To address the objectives, we first introduce the study area and describe the source of the data and process them. Then, we use the dynamic degree index of land utilization to quantitatively analyze the variation in the temporal and spatial characteristics of LUCC. Furthermore, the study used an equivalence factor to determine the ESV related to the physical geography of the BSRNC and then assessed the results. Following this, we discuss the relationship between the ESV and LUCC and present a policy to protect ecosystem function and the ecological environment. Finally, we put forward conclusions and directions for future studies.

2. Materials and Methods

2.1. Study Area

The BSRNC is located between the latitudes of 41°04′–51°37′ N and the longitudes of 120°30′–134°06′ E (Figure 1). The topography conditions in this area are complex, surrounding the Greater Khingan, Lesser Khingan, and Changbai Mountains. The regional landform mainly comprises hills and mountains distributed in the Sanjiang, Songnen, and Liaohe Plains. The BSRNC belongs to the temperate-continental monsoon climate zone and has average annual precipitation ranging from 500–700 mm/year. Furthermore, rainfall is mainly concentrated in the crop-growing season, from April to September, accounting for 80% to 90% of the annual precipitation. Suitable and superior natural geographical conditions promote the stable development of regional agriculture, making the BSRNC a substantial grain base [38]. The BSRNC plays a key role in China’s ecological conservation and food security [41]. However, the growing human population has significantly prompted LUCC in the region, resulting in a reclamation rate in the BSRNC of about 75%. Additionally, unsustainable human land-use practices contribute to the degradation of cultivated land quality in black soil significantly, threatening national food security. Biodiversity, nutrient maintenance, hydrological networks, and other ecological services all suffer as a result of long-term and spatial reclamation. Therefore, we chose the BSRNC as the research area because it is significant for national food security and we want to strengthen ecological protection in this area.

2.2. Data Sources

This study acquired the 1:100 million spatial distribution data of Chinese soil types from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn, accessed on 16 April 2022). We extracted the spatial distribution data of black soil from the spatial distribution data on soil types. This paper determined the location of the BSRNC according to the list outlining the protection and utilization of black land and the principle of centralized contiguous land (Figure 1). Landsat TM/ETM remote sensing images were the source of the LUCC data used in this study. As such, the spatial resolution was 30 × 30 m, and the comprehensive accuracy was higher than 93% [42,43].
This paper obtained land-use type data from the Chinese remote sensing monitoring database of land-use status. In the BSRNC, cultivated land (paddy and dry land) accounts for the highest area. The mixed coniferous broad-leaved, broad-leaved, and bush forest types are the main components of the forest ecosystem. The grassland ecosystem is dominated by high grassland coverage. The water area classification includes rivers, lakes, tidal flats, and beaches. According to the land source characteristics in the BSRNC and the results of a secondary ecosystem classification by Xie et al. [44], we divided the land-use types into six categories: paddy, dry land, forest, grassland, water area, and unused land. The grain output and area data were from the 2010 and 2020 editions of the Heilongjiang Statistical Yearbook, Liaoning Statistical Yearbook, and Jilin Statistical Yearbook. The grain price data were sourced from China’s 2020 Agricultural Price Survey Yearbook [45].

2.3. Methods

LUCC dynamics alter the variations in the ESV, which can cause the degradation of ecosystem services [9,21]. This paper takes LUCC as the breakthrough point to study the changes in the ESV in the BSRNC from 1990 to 2020. We interpreted the land use dynamic index and land use dynamic degree to quantitatively analyze the temporal and spatial variation characteristics of LUCC. We quoted the unit price of food production service and value coefficient (VC) of farmland to calculate the ESV and quantitatively analyze the variation in its spatio-temporal characteristics. Then, the sensitivity coefficients of various land-use types were determined to verify the feasibility of the evaluation results. Further, we revealed the regularities in LUCC based on physical geographical location to lay an essential foundation to further explore ESV changes in the BSRNC. Then, according to Xie et al. [15,16] and the actual situation in the BSRNC, this paper does not consider the changes in construction land and does not involve the calculation of its ESV.

2.3.1. Quantitative Calculation of LUCC

This study analyzed LUCC by comparing the areas occupied by each land-use type during the study period. We calculated total area gains and losses across each land-use type. In this article, we chose two indices, the land use dynamic index and land use dynamic degree [25], to analyze LUCC in the BSRNC using the following equation [46,47]:
D = U y U x U x × 1 T
where D refers to the dynamic degree of a certain land-use type; Uy is the area of a certain land-use type at the end of the study period; Ux is the area at the beginning of the study period; T is the year during the research period.

2.3.2. Measurement of the ESV

a. Equivalent factor correction
In this article, the equivalent factors for different land-use types were determined according to Xie et al.’s 2017 study [16] and the geographical characteristics of the BSRNC. Scholars have also used Xie’s method to revise the equivalent factor when conducting studies in northeast China [48]. The study used the data collected by other scholars to correct the equivalence factor. Because the farmland comprising the study area is basically paddy fields and dry land, the equivalence of the paddy field and dry land factors correspond to those for paddy and dry land in Xie et al. [16]. The forests in the study area include the coniferous, mixed coniferous broad-leaved, broad-leaved, and bush forest types. Therefore, the forest factor is the average value of coniferous, mixed coniferous broad-leaved, broad-leaved, and bush forest factors. Grassland basically is prairie, so the grassland factor adopts the value of the prairie factor. The water area factor is the value of water. Most of the unused land is bare land, and the unused land factor adopts the value of the bare land factor. This research examined whether physical geographical factors influence changes in ecological land use. As a result, the changes in the construction land area and ESV are not considered in this paper. The equivalence factors used in this study are shown in Table 1.
b. Unit Price of Food Production Service on Cultivated Land
This article revised the ESV of the cultivated land ecosystem per unit area by referring to the research results of Xie et al. [16]. Meanwhile, the equivalence factor of the ESV per unit area was defined as being equal to 1/7 of the average economic value of the grain production of cultivated land. The formula is as follows:
E a = 1 / 7 i = 1 n m i p i q i M , i = 1 , 2 n
where Ea is the economic value of grain production per unit area of cultivated land; i is crop species; pi is the average price of the crop i (RMB/kg); qi is the yield per unit area of the crop i (kg/ha); mi is the sowing area of the crop i (ha); and M is the sown area of all crops(ha).
This study chose rice, maize, and soybean as the main crops and counted their yield, sowing area, yield per unit area, and price in 2020, which was finally calculated to be CNY 2433.86/ha.
c. Value coefficient
To determine the equivalent value coefficient, we referred to Costanza et al. [1,14] and Xie et al. [15]. We adopted the value coefficient to calculate the ESV per unit area of different land-use types. It was calculated as follows:
V C i j = e i j × E a , i , j = 1 , 2 n
where VCij refers to the unit price of the ecological service i of ecosystem j (CNY/ha); i refers to the type of ecological service function; j refers to the ecosystem type; eij provides the equivalence factor of the unit price of the ecological service i in ecosystem j.
d. Regional ESV
Based on the land-use data from 1990 to 2020 and the VC of different land types, we calculated the ESV in the BSRNC to provide a data reference for ecosystem restoration. The calculation formula of ESV is as follows [14]:
E S V = A k × V C k , E S V = A k × V C k f
In the formula, ESV represents the ecosystem service value of non-construction land. Ak is the distribution area of land-use type K in the study area. VCk is the ecosystem service value coefficient per unit area of land-use type K. VCkf is the value of the ecological service function of land-use type K.

2.3.3. Sensitivity Analysis

Sensitivity analysis was used to certify the elasticity between the total ESV and VC of different land-use types, representing a significant portion of ESV evaluation [9]. The economic concept of the elasticity coefficient was used to calculate the sensitivity of the ESV to determine the dependence of the ESV. The specific calculation equation is as follows [9,49,50]:
C S = ( E S V n E S V m ) / E S V m ( V C n l V C m l ) / V C m l
where CS is the sensitivity coefficient. ESVm and ESVn are the initial and adjusted ESV, respectively. VCml and VCnl are the initial and adjusted VC, respectively.
If CS is greater than 1, it indicates that the ESV is elastic to the VC and that the research results have poor accuracy and low credibility. On the contrary, if CS is less than one, it means that ESV is inelastic to VC and that the research results are credible [51,52,53].

3. Results

3.1. Land-Use Change and Dynamics

3.1.1. Temporal Land-Use Changes

Land use in the BSRNC changed significantly from 1990 to 2020. Cultivated land (paddy and dry land) and forestland were the dominant land-use types in the BSRNC. These two land-use types represented more than 80% of the total area. Out of all of the land-use types, cultivated land accounted for the largest proportion of the total area coverage (>48%), followed by forestland (35%), while water area, grassland, and unused land account for a smaller proportion of area coverage.
Figure 2 shows that the changes that took place between land-use types in the study area were complex. From 1990 to 2020, the paddy land area increased continuously from 2.80 × 106 ha to 4.52 × 106 ha. The grassland land area exhibited a continuous decline during the study period, from 4.01 × 106 ha to 1.54 × 106 ha, decreasing by 2.47 × 106 ha. In particular, grassland area decreased sharply in 2005–2020 and 2015–2020, from 3.36 × 106 ha to 2.35 × 106 ha and from 2.30 × 106 ha to 1.54 × 106 ha, respectively. During the whole study period, dry land significantly increased from 16.58 × 106 ha to 17.68 × 106 ha, increasing by 1.10 × 106 ha. In the same sense, from 1990 to 1995, dry land increased from 16.58 × 106 ha to 17.72 × 106 ha, increasing by 1.14 × 106 ha. Following this trend, we noted that dry land increased from 17.11 × 106 ha to 17.68 × 106 ha during the period of 2015–2020, increasing by 0.57 × 106 ha.
From 1990 to 2020, the forest land area declined significantly, from 17.48 × 106 ha to 16.40 × 106 ha, decreasing by 1.08 × 106 ha. As for water, land area showed a slowly declining trend during the study period. The overall change was not obvious, and its area decreased by 0.32 × 106 ha. Unused land showed a slowly increasing trend from 1990 to 2020. With the exception of the period from 2005 to 2010, the changes in unused land were not more obvious than they were in other years. Among all of the land-use types, unused land accounts for a very small proportion, representing only 0.08% of the total.
Overall, essential transitions happen between cultivated land, forestland, and grassland. Due to large-scale reclamation, more forestland and grassland are converted into cultivated land. Thus, the increase in cultivated land is the largest in the BSRNC. The percentage of cultivated land increased from 41.6% in 1990 to 48.7% in 2020, increasing by 7.1%. The forest and grassland decreased significantly. The percentage of forestland decreased from 37.5% in 1990 to 35.2% in 2020.

3.1.2. Land-Use Spatial Change

Notable land use/cover conversions were identified between 1990 and 2020 (Figure 3). As depicted in Figure 3, cultivated land (paddy and dry land) expansion was clearly observed across all of the BSRNC during the study period. However, extensive deforestation was observed in the five geographic regions. Most of the forest types were lost in the Greater and Lesser Khingan Mountains and in the Songnen Plain. In contrast, cultivated land expansion occurred in the Liaohe Plain, Sanjiang Plain, and Songnen Plain. Moreover, during the research period, the spatial distribution of grassland changed dramatically due to grassland expansion in the Sanjiang Plain, which replaced cultivated land. Unused land became distributed throughout the study area.
Additionally, a large amount of forest was converted to cultivated land in the Songnen Plain from 1990 to 1995. Then, the conversion of dry land to paddy is significant. This conversion was mainly concentrated in the Sanjiang Plain from 1995 to 2000. The main drivers for this were to ensure food security policies and respond to climate change. From 2005 to 2010, a direct conversion of forest and unused land to cultivated land was noted in the Sanjiang Plain. From 2015 to 2020, dry land and paddy conservation were found in the Songnen Plain. During this period, paddy field expansion was the result of increasing demands for food due to population growth.

3.1.3. Land-Use Dynamic Changes

Several significant changes in the dynamic degree of land-use can be observed in the whole study area (Table 2). An overall increasing trend can be observed in paddy, dry land, and unused land, with paddy land increasing the most, up to 1.72 × 106 ha, with the rate of change being 60.9% and a dynamic degree of 2.03%. However, forest, grassland, and water area show a decreasing trend. Grassland decreased the most, dropping by −2.47 × 106 ha, with the rate of change being −60.6% and the dynamic degree being −2.02%. Forest decreased by 1.08 × 106 ha, with the rate of change being −6.18% and the dynamic degree being −0.20%.
Overall, the paddy land area increased by 60.9% from 1990 to 2020, in contrast to a minor decrease between 2000 to 2005 (−1.99%). The highest paddy change dynamic (18.03%) was observed between 2010 to 2015. This increase could have been attributed to the country’s policies targeted at enhancing food security. Grassland decreased significantly from1990 to 2020 (−60.6%), while a slight increase was witnessed from 2000 to 2005 (0.21%). In the forest, there seems to be a slight decreasing trend of −6.18% from 1990 to 2020.
Additionally, the change trends in the water land area and unused land are fundamentally opposed to one another. The water land area decreased by −24.65%, and the unused land area increased by 24.08% during the study period. From 2015 to 2020, the relative change in the grassland and water land area was the most significant, changing by −33.12% and −25.25%, respectively. This situation showed that because of the rising prices of agricultural products and the government’s agricultural development policies, local governments and agricultural producers continue to impact ecological land in order to develop cultivated land due to the multiple incentives of economic benefits.

3.2. Variations of ESV

3.2.1. Temporal Variation Analysis of ESV

From 1990 to 2020, the total ESV showed negative growth in the BSRNC. As shown in Table 3, the total ESV decreased from CNY 14,792.86 million to CNY 14,184.90 million. In fact, during the study period, forest land had the highest ESV (>57%). However, unused land contributed to the total ESV the least, contributing about 0.001%. The ESV provided by cultivated land has increased continuously over the past 30 years, from CNY 1883.59 million to CNY 2154.35 million. The proportion of arable land increased from 12.7% to 16.3%.
However, the ESV of forest and grassland decreased significantly, from CNY 8407.49 million to CNY 7887.65 million and from CNY 494.61 million to CNY 189.45 million, respectively. The results showed that the cultivated land area experienced an obvious increase during the study period. The ecosystem’s unit area value was moderate, and it could not counterbalance the decrease in the region’s ESV induced by forest and grassland losses. The ESV of water area and unused land fluctuated from 1990 to 2020. The ESV of water area fell by CNY 57.15 million, whereas the ESV of unused land grew by CNY 3.43 million.
Only paddy, dry land, and unused land showed a growing trend in the ESV among the six land categories that were investigated. The ESVs of other land-use types fell, particularly those of grassland and water area, indicating that the ecosystems in the BSRNC are under tremendous downward pressure and require further adjustments.

3.2.2. ESV Spatial Variation Analysis

From 1990 to 2020, the amount and distribution of land types in the BSRNC changed considerably, which had a significant impact on the spatial change in the ESV. In particular, cultivated land expansion impacted the spatial distribution changes in the ESV.
As shown in Figure 4, we found that the area with the most significant decline in the ESV was located in the Greater and Lesser Khingan Mountains and Changbai Mountains, and this was mainly because of the massive forest land loss in these mountains. The expansion of cultivated land led to a decline in forest area in the region—numerous forest land types with high-value equivalence coefficients result in a high ESV. However, the cultivated land low-value equivalence coefficients result in a low ESV. The increase (CNY 163.35 million RMB) in the ESV caused by the expansion of arable land cannot compensate for the decrease (-CNY 519.84 million) in the ESV caused by the decrease in the forest area. As a result, the ESV in these mountain regions has been declining, seriously affecting the ecosystem services and ecosystem functions in these regions.
From 1990 to 2005, the ESV continuously increased in the Sanjiang Plain and Liaohe Plain, which was mainly due to cultivated land. In 2005, the implementation of returning farmland to forest and grassland restricted the conversion of forest into cultivated land, which enabled the ESV of the region to experience negative growth. After 2015, implementing the national food security strategy and establishing high-standard farmland resulted in the positive growth of the cultivated land ecosystem.

3.2.3. Value Analysis of Various Land Ecosystem Service Functions

Table 4 shows the changes in the ecosystem services for various land types. Regulation services contributed the most to the proportion of the ESV, while cultural services accounted for the lowest proportion of the total ESV from 1990 to 2020. In addition, a continuously increasing trend was documented in regulation services, while provision, support, and culture services showed fluctuating changes in the past 30 years.
The highest to lowest proportion of the service value of the land ecosystem for each year comprises hydrological regulation, climate regulation, gas regulation, soil formation and retention, biodiversity, environmental purification, food production, aesthetic landscape, raw material production, water supply, and nutrient cycling. The ESV of hydrological and climate regulation accounted for more than 50% of the total ESV. From 1990 to 2020, continuous increases were found in food production, increasing from CNY 577.35 million to CNY 638.36 million. The value of other ecosystem service functions decreased at different rates. Hydrological regulation experienced a significant drop of CNY 836.46 million during the study period. From 2015 to 2020, there was an obvious decrease in the water supply that was determined to be closely related to the rapid decline of the forest and grassland areas.
Meanwhile, the expansion of cultivated land in the BSRNC is essential for food production, but this expansion leads to ecosystem destruction. A continuous decline was observed in the water supply, gas regulation, climate regulation, environmental purification, and hydrological regulation. Climate regulation dropped from CNY 2817.18 million in 1990 to CNY 2601.94 million in 2020. Environmental purification dropped from CNY 988.97 million in 1990 to CNY 886.60 million in 2020. A series of data shows that the expansion of cultivated land leads to a decline in the ESV.

3.3. Sensitivity Coefficient Analysis

The Coefficient of Sensitivity (CS) was used to evaluate the sensitivity and consistency of our ESV estimations. As such, the assessed ecosystem service value for 1990 to 2020 is presented in Table 5. This result shows different aspects, and there are various land-use types.
The sensitivity coefficients of the ESV to the VC for diverse land-use types were less than one. In fact, from 1990 to 2020, water area recorded its highest CS of 0.82% in 2020. During the same period, the VC of unused land was 0.01%. Land-use types such as paddy and dry land, which represent high-land-use area, incurred a moderate value. Paddy land recorded a maximum value of 0.35% in 1990. Alternatively, the value VC of dry land remains the same over the years at 0.06%.
Consequently, if we consider the order of the degree of sensibility, we can note the following order: water area > grassland > paddy > forest > dry land > unused land. This context demonstrated that the ESV in the BSRNC lacked elasticity to VC and that the results were reliable. In conclusion, the analysis revealed that the paddy VC declined significantly from 0.36% in 1990 to 0.22% in 2020, a differential of −0.14%. Alternatively, the VC of the water area increased from 0.77% in 1990 to 0.82% in 2020.

4. Discussion

4.1. Interpretation of LUCC in the BSRNC

Over the past 30 years, a large number of forest and grassland areas have been replaced with cultivated land. The area of cultivated land increased from 19.38 × 106 ha in 1990 to 22.21 × 106 ha in 2020, increasing by 3.83 × 106 ha. There have changes between dry land and paddy areas, most of which involve the conversion of dry land to paddy. The paddy land area in the BSRNC continued to expand from 2.80 × 106 ha to 4.52 × 106 ha, increasing by 1.72 × 106 ha during the study period. Mao et al. [40] have shown that China has also implemented a series of policies to promote agricultural development and grain production. For example, China launched a food security project to increase grain production and ensure food security in the face of a rapidly growing population [23].
Additionally, the implementation of conversion from dry land to paddy has caused wetland and grassland loss and degradation, which may lead to the problem of water shortage [54]. Natural resources, such as water, may have a significant influence on this land transition, notably paddy expansion for increased food production. We found that there were abundant water resources to increase the paddy area in the BSRNC, which has an average surface runoff of 2634.3 × 108 m3/year and precipitation of 567.5 × 108 m3/year [36]. In this context, it is evident that water supply resources provide favorable conditions for agricultural production, which can also explain the expansion of paddy land area during the study period. Meanwhile, climate data from meteorological stations in the colder areas of China displayed that the air temperature warming rate was 0.1 °C/decade during 1986–2015, with intensified warming (0.2 °C/decade) after the 2000s. Climatic warming increased the accumulated temperature during the rice growing season and stimulated the climate adaptability of crops in this region [55]. Under these factors, grain production has increased rapidly, but the sustainable utilization of resources and the effective management of ecosystems face serious challenges.
Intense LUCC has a negative impact on the ecosystem and ecological environment. China developed larger-scale ecosystem restoration policies, particularly in the BSRNC. These policies include the Grain for Green Policy and the Natural Forests Conservation Policy [30,56]. Although many ecosystem restoration policies have been implemented, the improvement of ecosystem services still needs more time [57]. Hence, analyzing the relationship between the ESV and LUCC is indispensable for ecologically sustainable development.

4.2. The Relationship between ESV and LUCC

Many studies have shown that LUCC significantly affects changes in the ESV [39,58]. In particular, irrational land-use change can lead to ecosystem degradation. Many studies have used the CS to evaluate the sensitivity and consistency of their ESV estimations. As such, this study used the CS to assess the relationship between the ESV and LUCC. Based on the results of the sensitivity analysis, the sensitivity coefficients of ESV and VC for different land-use types were less than one. These results demonstrated that LUCC is closely related to ecosystem services in the BSRNC.
The intrinsic relationship between humans and natural systems was effectively shown by examining the LUCC process and its associated consequences on ecosystem services. Furthermore, the expansion of the cultivated land area to promote food production has led to a significant decrease in hydrological, climate, and gas regulation and as well as a decrease in the total ESV. Our results indicate that between 1990 and 2020, the ecosystem providing the ESV decreased by CNY 607.96 million. The significant increase in food production during the study period caused paddy expansion. Because of the high demand for water in paddy areas, the water supply has shrunk, decreasing by CNY 187.56 million. Projects to increase the grain production capacity and to improve rice yield in order to promote food grain have also led to forest and grassland degradation, leading to a significant decrease in climate regulation, environmental purification, hydrological regulation, and soil formation and retention. In the BSRNC, forest areas decreased substantially, especially due to the conversion to cultivated land, and global warming benefits from the reclamation of cultivated land [59]. There is also a need to feed the increasing population [41]. The conversion from forest land to other land-use types has contributed to local ESV loss the most. Analyzing these trends raises the question of whether it is possible to compare ecosystem services over time. Since ecosystem services are related to the value society assigns to the goods and services produced by nature, the exact delivery of a service might be valued quite differently over time. To avoid these problems, we focused on comparing the potential of the land-use pattern to provide ecosystem functions.
In this study, the decrease in the ESV is mainly due to the expansion of cultivated land, which resulted in forest and grassland losses. The two land types are considered significant providers of ecosystem services. Large-scale deforestation and reclamation have led to the deterioration of the ecosystem, which has a seriously impact on the ecological environment in the Greater Khingan Mountains and Lesser Khingan Mountains. With the Songnen Plain and Sanjiang Plain at the core of grain output, the region has focused on developing modern agriculture and expanding cultivated land. Notably, after 2005, agricultural producers have still continued to impact ecological land to develop cultivated land due to the government’s agricultural development policy in the Sanjiang Plain. Hence, the deterioration of ecological land leads to reductions in regional natural resources and the degradation of ecosystem services such as climate regulation, gas regulation, water supply, and biodiversity.
LUCC significantly drives variations in the ESV at the national scale [60,61,62]. Ouyang et al. [30] researched improved ecosystem services from 2000 to 2010 at the country scale and found that ecosystem services were distinctive in different regions due to spatial heterogeneity and diverse policies [63,64]. Rapid socioeconomic development and policy changes may accelerate LUCC or transfer land into a new type. Concurrently, many developed and developing regions are experiencing a decline in cultivated land and an increase in urban land [51,65,66]. Therefore, many scholars focus on the relationship between urban land expansion and ESV [60,66]. However, BSRNC, a major grain-producing area in China, is mainly manifested in the relationship between the disorderly expansion of cultivated land and ESV. In terms of the transition of forest and grassland to cultivated land, this transfer facilitates the rapid decline in the ESV, which may have a negative impact on the regional ecological environment and system.

4.3. Regional Land-Use Policy Implications for Ecological Protection

The BSRNC represents important forest resources and a grain production base, with large-scale reclaimed forest accounting for over 1000 ha since 1990 [33]. Land consolidation projects and deforestation activities have significantly increased the amount of potential cultivated land in the BSRNC. As a result of these measures, cultivated land has rapidly expanded while forest and grassland have reduced significantly, resulting in a substantial impact on ecosystems. This is related to the government’s implementation of several major programs after the year 2000, including the Grain for Green project and the natural forest protection project [30,56]. The implementation of these ecological projects has significantly improved the ecological environment but has also led to declines in the cultivated land area. From 2000 to 2005, the decreased arable land accounted for 63.5% of the total land in China [67]. In this context, food security has once again become a core concern of the Chinese government. In 2010, China implemented the grain production capacity construction project to increase grain production and ensure food security. The local government recognized the Songnen and Sanjiang Plains as being vital to modern agriculture development. Even though the government has adopted several legislations to safeguard ecosystems, ecological conservation still faces significant challenges.
The regional functional orientation determines regional land-use change. Therefore, the conflict between food production and ecological protection leads to changes in the relative importance of cultivated land and ecological land resources in different periods and then changes the relative scarcity of different types of land resources. As the prominent grain-producing area in China, the BSRNC should emphasize the security of grain yield and how to increase it and should pay attention to regional ecological security as well as strengthen the ecological support for food security, rationally plan the ecological resources of different land types, protect the ecological land, and establish a dual assessment mechanism for food security and ecological security.
The decrease in the ESV varies in different regions, so the government should optimize the spatial patterns of the main functional zones for food and ecological security according to the actual regional conditions. In the ecological zones for the Greater Khingan Mountains and Lesser Khingan Mountains, land resources should be well protected, and the ecological balance should be effectively maintained. Additionally, comprehensive management should be carried out to ensure the coordinated development of agricultural production and the environmental environment. This will result in improved grain production in the Sanjiang Plain and Songnen Plain and leads down the road to sustainable agricultural development.
In summary, the synchronization of ecosystem conservation strategies is hindered by our lack of understanding regarding variation in the ecosystem service value and fundamental laws, spatiotemporal changes in biophysical conditions, and ecosystem structures. As such, for above this scenario, future LULC policies should pay more attention to balancing the relationship between agricultural land expansion and ecosystem services.

5. Conclusions

In this article, we analyzed the variations in the ESV under LUCC in the BSRNC from 1990 to 2020 using a value equivalence method. The results showed that land use in the BSRNC changed significantly from 1990 to 2020. Cultivated land continued to increase from 19.38 × 106 ha to 22.21 × 106 ha. However, grassland in the BSRNC continued to decrease, decreasing from 4.01 × 106 ha to 1.54 × 106 ha. Cultivated land such as paddy areas continued to increase, from 2.80 × 106 ha to 4.52 × 106 ha, increasing by 1.72 × 106 ha. Policy implementation is a crucial factor affecting the expansion of paddy areas and of the conversion of dry land into paddy areas. The ecosystem resulted in the ESV decreasing by CNY 607.96 million during the study period. Cultivated land, forest, and water body areas were the main land-use types that made the greatest contributions to the ESV. The change trends observed in the ESV for land use in the BSRNC are consistent with the changes in its main components. Then, the evaluation results of the ESV proved to be reliable due to the sensitivity indexes being less than one (<1). The sensitivity coefficients followed water body > grassland > paddy > forest > dry land > unused land in 2020, which indicates that the land-use changes lack flexibility in the changes in the ESV.
This study can provide more effective policy-making and data support for ecosystem restoration and environmental protection in China’s black land region. However, the research still has some limitations. We did not study the temporal and spatial changes of construction land or calculate its ESV. We also only discussed the impact of LUCC on the ESV and did not consider the impact of other factors on the ESV, such as urbanization. Consequently, to improve the accuracy of the assessment results, future studies should focus on the spatial-temporal changes in the ESV that are driven by both natural and economic factors. Moreover, the driving factors resulting in ESV changes in the BSRNC should also be considered in future research.

Author Contributions

Conceptualization, Q.L. and G.D.; methodology, Q.L.; validation, Q.L., L.W. and G.D.; formal analysis, L.W.; investigation, Q.L., L.W., and J.L.; resources, Q.L. and G.D.; data curation, Y.L. and W.L.; writing—original draft preparation, L.W.; writing—review and editing, Q.L., G.D. and B.F.; visualization, L.W.; supervision, Q.L. and G.D.; project administration, Q.L.; funding acquisition, Q.L. and S.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant numbers 41901208, National Natural Science Foundation of China, grant number 42101217, China Postdoctoral Science Foundation, grant number 2021M700738.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy or other restrictions.

Acknowledgments

We thank Yunkai Li and Bonoua Faye for their detailed and helpful reviews of an early draft of the manuscript. We also express our gratitude to the professionals at Northeast Agricultural University who participated in the research and survey.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Costanza, R.; de Groot, R.; Sutton, P.; van der Ploeg, S.; Anderson, S.J.; Kubiszewski, I.; Farber, S.; Turner, R.K. Changes in the global value of ecosystem services. Glob. Environ. Chang. 2014, 26, 152–158. [Google Scholar] [CrossRef]
  2. Schägner, J.P.; Maes, J.; Brander, L.M.; Hartje, V. Mapping Ecosystem Services’ Values: Current Practice and Future Prospects. SSRN J. 2013, 4, 33–46. [Google Scholar] [CrossRef] [Green Version]
  3. Whitmee, S.; Haines, A.; Beyrer, C.; Boltz, F.; Capon, A.G.; de Souza Dias, B.F.; Ezeh, A.; Frumkin, H.; Gong, P.; Head, P.; et al. Safeguarding human health in the Anthropocene epoch: Report of The Rockefeller Foundation–Lancet Commission on planetary health. Lancet 2015, 386, 1973–2028. [Google Scholar] [CrossRef]
  4. Cumming, G.S.; Buerkert, A.; Hoffmann, E.M.; Schlecht, E.; von Cramon-Taubadel, S.; Tscharntke, T. Implications of agricultural transitions and urbanization for ecosystem services. Nature 2014, 515, 50–57. [Google Scholar] [CrossRef] [PubMed]
  5. Shang, Y.; Wang, D.; Liu, S.; Li, H. Spatial-Temporal Variation and Mechanisms Causing Spatial Differentiation of Ecosystem Services in Ecologically Fragile Regions Based on Value Evaluation: A Case Study of Western Jilin, China. Land 2022, 11, 629. [Google Scholar] [CrossRef]
  6. Valle Junior, R.F.; Varandas, S.G.P.; Sanches Fernandes, L.F.; Pacheco, F.A.L. Environmental land use conflicts: A threat to soil conservation. Land Use Policy 2014, 41, 172–185. [Google Scholar] [CrossRef]
  7. Sandifer, P.A.; Sutton-Grier, A.E.; Ward, B.P. Exploring connections among nature, biodiversity, ecosystem services, and human health and well-being: Opportunities to enhance health and biodiversity conservation. Ecosyst. Serv. 2015, 12, 1–15. [Google Scholar] [CrossRef] [Green Version]
  8. Loomes, R.; O’Neill, K. Nature’s Services: Societal Dependence on Natural Ecosystems. Pac. Conserv. Biol. 2000, 6, 274. [Google Scholar] [CrossRef] [Green Version]
  9. Kreuter, U.P.; Harris, H.G.; Matlock, M.D.; Lacey, R.E. Change Ecosyst. Serv. Values San Antonio Area Texas. Ecol. Econ. 2001, 39, 333–346. [Google Scholar] [CrossRef]
  10. Morri, E.; Santolini, R. Ecosystem Services Valuation for the Sustainable Land Use Management by Nature-Based Solution (NbS) in the Common Agricultural Policy Actions: A Case Study on the Foglia River Basin (Marche Region, Italy). Land 2022, 11, 57. [Google Scholar] [CrossRef]
  11. Don, A.; Schumacher, J.; Freibauer, A. Impact of tropical land-use change on soil organic carbon stocks—A meta-analysis: Soil organic carbon and land-use change. Glob. Chang. Biol. 2011, 17, 1658–1670. [Google Scholar] [CrossRef] [Green Version]
  12. Lambin, E.F.; Geist, H. (Eds.) Land-Use and Land-Cover Change: Local Processes and Global Impacts; Global Change—The IGBP Series; Springer: Berlin/Heidelberg, Germany, 2006; ISBN 978-3-540-32201-6. [Google Scholar]
  13. Taylor, M.J.; Aguilar-Støen, M.; Castellanos, E.; Moran-Taylor, M.J.; Gerkin, K. International migration, land use change and the environment in Ixcán, Guatemala. Land Use Policy 2016, 54, 290–301. [Google Scholar] [CrossRef]
  14. Costanza, R.; d’Arge, R.; de Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Ecol. Econ. 1998, 25, 3–15. [Google Scholar] [CrossRef]
  15. Xie, G.; Zhen, L.; Lu, C.; Xiao, Y.; Chen, C. Expert knowledgy based valuation method of ecosystem services in China. Nat. Res. 2008, 23, 911–919. [Google Scholar]
  16. Xie, G.; Zhang, C.; Zhen, L.; Zhang, L. Dynamic changes in the value of China’s ecosystem services. Ecosyst. Serv. 2017, 26, 146–154. [Google Scholar] [CrossRef]
  17. Tong, X.; Brandt, M.; Yue, Y.; Horion, S.; Wang, K.; Keersmaecker, W.D.; Tian, F.; Schurgers, G.; Xiao, X.; Luo, Y.; et al. Increased vegetation growth and carbon stock in China karst via ecological engineering. Nat. Sustain. 2018, 1, 44–50. [Google Scholar] [CrossRef]
  18. Xing, L.; Xue, M.; Wang, X. Spatial correction of ecosystem service value and the evaluation of eco-efficiency: A case for China’s provincial level. Ecol. Indic. 2018, 95, 841–850. [Google Scholar] [CrossRef]
  19. Wang, H.; Zhou, S.; Li, X.; Liu, H.; Chi, D.; Xu, K. The influence of climate change and human activities on ecosystem service value. Ecol. Eng. 2016, 87, 224–239. [Google Scholar] [CrossRef]
  20. Fei, L.; Shuwen, Z.; Jiuchun, Y.; Liping, C.; Haijuan, Y.; Kun, B. Effects of land use change on ecosystem services value in West Jilin since the reform and opening of China. Ecosyst. Serv. 2018, 31, 12–20. [Google Scholar] [CrossRef]
  21. Hu, H.; Liu, W.; Cao, M. Impact of land use and land cover changes on ecosystem services in Menglun, Xishuangbanna, Southwest China. Environ. Monit. Assess. 2008, 146, 147–156. [Google Scholar] [CrossRef]
  22. Polasky, S.; Nelson, E.; Pennington, D.; Johnson, K.A. The Impact of Land-Use Change on Ecosystem Services, Biodiversity and Returns to Landowners: A Case Study in the State of Minnesota. Environ. Resour. Econ. 2011, 48, 219–242. [Google Scholar] [CrossRef]
  23. Mao, D.; He, X.; Wang, Z.; Tian, Y.; Xiang, H.; Yu, H.; Man, W.; Jia, M.; Ren, C.; Zheng, H. Diverse policies leading to contrasting impacts on land cover and ecosystem services in Northeast China. J. Clean. Prod. 2019, 240, 117961. [Google Scholar] [CrossRef]
  24. Song, F.; Su, F.; Zhu, D.; Li, L.; Li, H.; Sun, D. Evaluation and driving factors of sustainable development of the wetland ecosystem in Northeast China: An emergy approach. J. Clean. Prod. 2020, 248, 119236. [Google Scholar] [CrossRef]
  25. Su, K.; Wei, D.; Lin, W. Evaluation of ecosystem services value and its implications for policy making in China—A case study of Fujian province. Ecol. Indic. 2020, 108, 105752. [Google Scholar] [CrossRef]
  26. Sonter, L.J.; Johnson, J.A.; Nicholson, C.C.; Richardson, L.L.; Watson, K.B.; Ricketts, T.H. Multi-site interactions: Understanding the offsite impacts of land use change on the use and supply of ecosystem services. Ecosyst. Serv. 2017, 23, 158–164. [Google Scholar] [CrossRef] [Green Version]
  27. Benayas, J.M.R.; Newton, A.C.; Diaz, A.; Bullock, J.M. Enhancement of Biodiversity and Ecosystem Services by Ecological Restoration: A Meta-Analysis. Science 2009, 325, 1121–1124. [Google Scholar] [CrossRef]
  28. Ma, S.; Qiao, Y.-P.; Jiang, J.; Wang, L.-J.; Zhang, J.-C. Incorporating the implementation intensity of returning farmland to lakes into policymaking and ecosystem management: A case study of the Jianghuai Ecological Economic Zone, China. J. Clean. Prod. 2021, 306, 127284. [Google Scholar] [CrossRef]
  29. Barral, M.P.; Rey Benayas, J.M.; Meli, P.; Maceira, N.O. Quantifying the impacts of ecological restoration on biodiversity and ecosystem services in agroecosystems: A global meta-analysis. Agric. Ecosyst. Environ. 2015, 202, 223–231. [Google Scholar] [CrossRef] [Green Version]
  30. Ouyang, Z.; Zheng, H.; Xiao, Y.; Polasky, S.; Liu, J.; Xu, W.; Wang, Q.; Zhang, L.; Xiao, Y.; Rao, E.; et al. Improvements in ecosystem services from investments in natural capital. Science 2016, 352, 1455–1459. [Google Scholar] [CrossRef]
  31. Feng, Q.; Zhao, W.; Hu, X.; Liu, Y.; Daryanto, S.; Cherubini, F. Trading-off ecosystem services for better ecological restoration: A case study in the Loess Plateau of China. J. Clean. Prod. 2020, 257, 120469. [Google Scholar] [CrossRef]
  32. Wu, C.; Chen, B.; Huang, X.; Dennis Wei, Y.H. Effect of land-use change and optimization on the ecosystem service values of Jiangsu province, China. Ecol. Indic. 2020, 117, 106507. [Google Scholar] [CrossRef]
  33. Zhang, G.; Liu, J.; Zhang, Z. Analysis of Deforested-Land Reclamation during the Late 20th Century in China Based on Remote Sensing and GIS. Geogr. Res. 2003, 22, 221–228. [Google Scholar]
  34. Yin, D.; Li, X.; Li, G.; Zhang, J.; Yu, H. Spatio-Temporal Evolution of Land Use Transition and Its Eco-Environmental Effects: A Case Study of the Yellow River Basin, China. Land 2020, 9, 514. [Google Scholar] [CrossRef]
  35. Kc, A.; Wagle, N.; Acharya, T.D. Spatiotemporal Analysis of Land Cover and the Effects on Ecosystem Service Values in Rupandehi, Nepal from 2005 to 2020. ISPRS Int. J. Geo-Inf. 2021, 10, 635. [Google Scholar] [CrossRef]
  36. Pan, T.; Zhang, C.; Kuang, W.; De Maeyer, P.; Kurban, A.; Hamdi, R.; Du, G. Time Tracking of Different Cropping Patterns Using Landsat Images under Different Agricultural Systems during 1990–2050 in Cold China. Remote Sens. 2018, 10, 2011. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, Z.; Song, K.; Ma, W.; Ren, C.; Zhang, B.; Liu, D.; Chen, J.M.; Song, C. Loss and Fragmentation of Marshes in the Sanjiang Plain, Northeast China, 1954–2005. Wetlands 2011, 31, 945–954. [Google Scholar] [CrossRef]
  38. Pan, T.; Zhang, C.; Kuang, W.; Luo, G.; Du, G.; Yin, Z. Large-scale rain-fed to paddy farmland conversion modified land-surface thermal properties in Cold China. Sci. Total Environ. 2020, 722, 137917. [Google Scholar] [CrossRef]
  39. Yan, H.; Liu, J.; Huang, H.Q.; Tao, B.; Cao, M. Assessing the consequence of land use change on agricultural productivity in China. Glob. Planet. Chang. 2009, 67, 13–19. [Google Scholar] [CrossRef]
  40. Mao, D.; Luo, L.; Wang, Z.; Wilson, M.C.; Zeng, Y.; Wu, B.; Wu, J. Conversions between natural wetlands and farmland in China: A multiscale geospatial analysis. Sci. Total Environ. 2018, 634, 550–560. [Google Scholar] [CrossRef]
  41. Mao, D.; Wang, Z.; Wu, C.; Song, K.; Ren, C. Examining forest net primary productivity dynamics and driving forces in northeastern China during 1982–2010. Chin. Geogr. Sci. 2014, 24, 631–646. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, J.; Kuang, W.; Zhang, Z.; Xu, X.; Qin, Y.; Ning, J.; Zhou, W.; Zhang, S.; Li, R.; Yan, C.; et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
  43. Ning, J.; Liu, J.; Kuang, W.; Xu, X.; Zhang, S.; Yan, C.; Li, R.; Wu, S.; Hu, Y.; Du, G.; et al. Spatiotemporal patterns and characteristics of land-use change in China during 2010–2015. J. Geogr. Sci. 2018, 28, 547–562. [Google Scholar] [CrossRef] [Green Version]
  44. Xie, G.; Zhang, C.; Zhang, L.; Chen, W.; Li, S. lmprovement of the Evaluation Method for Ecosystem Service Value Based on Per Unit Area. Nat. Res. 2015, 30, 1243–1254. [Google Scholar]
  45. China Yearbook of Agricultural Price Survey. Available online: https://navi.cnki.net/knavi/yearbooks/YBAXS/detail (accessed on 21 April 2022).
  46. Wang, W.; Guo, H.; Chuai, X.; Dai, C.; Lai, L.; Zhang, M. The impact of land use change on the temporospatial variations of ecosystems services value in China and an optimized land use solution. Environ. Sci. Policy 2014, 44, 62–72. [Google Scholar] [CrossRef]
  47. Tan, Z.; Guan, Q.; Lin, J.; Yang, L.; Luo, H.; Ma, Y.; Tian, J.; Wang, Q.; Wang, N. The response and simulation of ecosystem services value to land use/land cover in an oasis, Northwest China. Ecol. Indic. 2020, 118, 106711. [Google Scholar] [CrossRef]
  48. Wang, L.-J.; Ma, S.; Zhao, Y.-G.; Zhang, J.-C. Ecological restoration projects did not increase the value of all ecosystem services in Northeast China. For. Ecol. Manag. 2021, 495, 119340. [Google Scholar] [CrossRef]
  49. Qiu, L.; Pan, Y.; Zhu, J.; Amable, G.S.; Xu, B. Integrated analysis of urbanization-triggered land use change trajectory and implications for ecological land management: A case study in Fuyang, China. Sci. Total Environ. 2019, 660, 209–217. [Google Scholar] [CrossRef]
  50. Yang, Y.; Wang, K.; Liu, D.; Zhao, X.; Fan, J. Effects of land-use conversions on the ecosystem services in the agro-pastoral ecotone of northern China. J. Clean. Prod. 2020, 249, 119360. [Google Scholar] [CrossRef]
  51. Tianhong, L.; Wenkai, L.; Zhenghan, Q. Variations in ecosystem service value in response to land use changes in Shenzhen. Ecol. Econ. 2010, 69, 1427–1435. [Google Scholar] [CrossRef]
  52. Liu, Y.; Li, J.; Zhang, H. An ecosystem service valuation of land use change in Taiyuan City, China. Ecol. Model. 2012, 225, 127–132. [Google Scholar] [CrossRef]
  53. Wang, S.; Wu, B.; Yang, P. Assessing the changes in land use and ecosystem services in an oasis agricultural region of Yanqi Basin, Northwest China. Environ. Monit. Assess. 2014, 186, 8343–8357. [Google Scholar] [CrossRef] [PubMed]
  54. Chen, H.; Wang, G.; Lu, X.; Jiang, M.; Mendelssohn, I.A. Balancing the Needs of China’s Wetland Conservation and Rice Production. Environ. Sci. Technol. 2015, 49, 6385–6393. [Google Scholar] [CrossRef] [PubMed]
  55. Li, Z.; Liu, Z.; Anderson, W.; Yang, P.; Wu, W.; Tang, H.; You, L. Chinese Rice Production Area Adaptations to Climate Changes, 1949–2010. Environ. Sci. Technol. 2015, 49, 2032–2037. [Google Scholar] [CrossRef] [PubMed]
  56. Wen, X.; Théau, J. Spatiotemporal analysis of water-related ecosystem services under ecological restoration scenarios: A case study in northern Shaanxi, China. Sci. Total Environ. 2020, 720, 137477. [Google Scholar] [CrossRef]
  57. Chazdon, R.L. Beyond Deforestation: Restoring Forests and Ecosystem Services on Degraded Lands. Science 2008, 320, 1458–1460. [Google Scholar] [CrossRef] [Green Version]
  58. Liu, B.; Pan, L.; Qi, Y.; Guan, X.; Li, J. Land Use and Land Cover Change in the Yellow River Basin from 1980 to 2015 and Its Impact on the Ecosystem Services. Land 2021, 10, 1080. [Google Scholar] [CrossRef]
  59. Li, X.; Zhao, G.; Li, B. Analysis of Changing Situation of Land Resources in Northeast China. Res. Soil Water Conserv. 2010, 17, 68–74. [Google Scholar]
  60. Zhang, J.; Qu, M.; Wang, C.; Zhao, J.; Cao, Y. Quantifying landscape pattern and ecosystem service value changes: A case study at the county level in the Chinese Loess Plateau. Glob. Ecol. Conserv. 2020, 23, e01110. [Google Scholar] [CrossRef]
  61. Talukdar, S.; Pal, S. Modeling flood plain wetland transformation in consequences of flow alteration in Punarbhaba river in India and Bangladesh. J. Clean. Prod. 2020, 261, 120767. [Google Scholar] [CrossRef]
  62. Jin, G.; Chen, K.; Liao, T.; Zhang, L.; Najmuddin, O. Measuring ecosystem services based on government intentions for future land use in Hubei Province: Implications for sustainable landscape management. Landsc. Ecol. 2021, 36, 2025–2042. [Google Scholar] [CrossRef]
  63. Wang, Y.; Dai, E.; Yin, L.; Ma, L. Land use/land cover change and the effects on ecosystem services in the Hengduan Mountain region, China. Ecosyst. Serv. 2018, 34, 55–67. [Google Scholar] [CrossRef]
  64. Shi, G.; Ye, P.; Ding, L.; Quinones, A.; Li, Y.; Jiang, N. Spatio-Temporal Patterns of Land Use and Cover Change from 1990 to 2010: A Case Study of Jiangsu Province, China. Int. J. Environ. Res. Public Health 2019, 16, 907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Prishchepov, A.V.; Müller, D.; Dubinin, M.; Baumann, M.; Radeloff, V.C. Determinants of agricultural land abandonment in post-Soviet European Russia. Land Use Policy 2013, 30, 873–884. [Google Scholar] [CrossRef] [Green Version]
  66. Qu, S.; Hu, S.; Li, W.; Wang, H.; Zhang, C.; Li, Q. Interaction between urban land expansion and land use policy: An analysis using the DPSIR framework. Land Use Policy 2020, 99, 104856. [Google Scholar] [CrossRef]
  67. Liu, X.; Zhao, C.; Song, W. Review of the evolution of cultivated land protection policies in the period following China’s reform and liberalization. Land Use Policy 2017, 67, 660–669. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Ijerph 19 07533 g001
Figure 2. Temporal variations in land-use in the BSRNC.
Figure 2. Temporal variations in land-use in the BSRNC.
Ijerph 19 07533 g002
Figure 3. Spatial distribution of land-use transfer in the BSRNC from 1990 to 2020. Notes: 1: Paddy; 2: dry land; 3: forest; 4: grassland; 5: water area; 6: unused land. (a) is the spatial distribution of land-use transfer in the BSRNC from 1990 to 1995; (b) is the spatial distribution of land-use transfer in the BSRNC from 1995 to 2000; (c) is the spatial distribution of land-use transfer in the BSRNC from 2000 to 2005; (d) is the spatial distribution of land-use transfer in the BSRNC from 2005 to 2010; (e) is the spatial distribution of land-use transfer in the BSRNC from 2010 to 2015; (f) is the spatial distribution of land-use transfer in the BSRNC from 2015 to 2020; (g) is the spatial distribution of land-use transfer in the BSRNC from 1990 to 2020.
Figure 3. Spatial distribution of land-use transfer in the BSRNC from 1990 to 2020. Notes: 1: Paddy; 2: dry land; 3: forest; 4: grassland; 5: water area; 6: unused land. (a) is the spatial distribution of land-use transfer in the BSRNC from 1990 to 1995; (b) is the spatial distribution of land-use transfer in the BSRNC from 1995 to 2000; (c) is the spatial distribution of land-use transfer in the BSRNC from 2000 to 2005; (d) is the spatial distribution of land-use transfer in the BSRNC from 2005 to 2010; (e) is the spatial distribution of land-use transfer in the BSRNC from 2010 to 2015; (f) is the spatial distribution of land-use transfer in the BSRNC from 2015 to 2020; (g) is the spatial distribution of land-use transfer in the BSRNC from 1990 to 2020.
Ijerph 19 07533 g003
Figure 4. Spatial variation of the ESV in the BSRNC. Notes: The subfigures represent the change in ESV for different land use types.
Figure 4. Spatial variation of the ESV in the BSRNC. Notes: The subfigures represent the change in ESV for different land use types.
Ijerph 19 07533 g004
Table 1. Equivalence value per unit area of the ecological service of non-construction land in the BSRNC.
Table 1. Equivalence value per unit area of the ecological service of non-construction land in the BSRNC.
Primary TypesSecondary TypesPaddyDry LandForestGrasslandWater AreaUnused Land
ProvisionFood production1.360.850.250.100.800.00
Raw material production0.090.400.580.140.230.00
Water supply−2.630.020.300.088.290.00
RegulationGas regulation1.110.671.910.510.770.02
Climate regulation0.570.365.711.342.290.00
Environmental purification0.170.101.670.445.550.10
Hydrological regulation2.720.273.740.98102.240.03
SupportSoil formation and retention0.011.032.320.620.930.02
Nutrient cycling0.190.120.180.050.070.00
Biodiversity0.210.132.120.562.550.02
CultureAesthetic landscape0.090.060.980.251.890.01
Total3.894.0119.765.07125.610.20
Table 2. Evolution of dynamic land-use types in the period of 1990–2020 (%).
Table 2. Evolution of dynamic land-use types in the period of 1990–2020 (%).
Periods Land-Use Types
PaddyDry LandForestGrasslandWater AreaUnused Land
1990–1995Relative change0.716.85−2.85−10.91−10.18−5.28
Dynamic degree0.141.37−0.57−2.18−2.04−1.06
1995–2000Relative change16.920.84−1.81−6.007.99−6.04
Dynamic degree3.380.17−0.36−1.201.60−1.21
2000–2005Relative change−1.990.72−0.060.21−0.77−2.59
Dynamic degree−0.400.14−0.010.04−0.15−0.52
2005–2010Relative change8.16−1.01−1.25−30.272.3839.49
Dynamic degree1.63−0.20−0.25−6.050.487.90
2010–2015Relative change18.03−3.92−0.01−2.132.30−2.53
Dynamic degree3.61−0.790.00−0.430.46−0.51
2015–2020Relative change9.663.33−0.34−33.12−25.255.28
Dynamic degree1.930.67−0.07−6.62−5.051.06
1990–2020Relative change60.906.64−6.18−60.60−24.6524.08
Dynamic degree2.030.22−0.20−2.02−0.810.79
Table 3. Ecosystem service values of different land-use types from 1990 to 2020 (CNY 100 million).
Table 3. Ecosystem service values of different land-use types from 1990 to 2020 (CNY 100 million).
Year
Land-Use Types
QuantityPaddyDry LandForestGrasslandWater AreaUnused Land
1990ESV265.381618.218407.49494.613992.9414.23
1995ESV267.261729.068167.85440.653786.2813.48
2000ESV312.471743.648020.02414.213672.9112.66
2005ESV306.241756.108015.32415.063842.9112.34
2010ESV331.221738.337915.14289.413934.3417.21
2015ESV390.941670.077914.26283.254024.9116.77
2020ESV428.731725.627887.65189.453935.7917.66
1990–1995Variation1.88110.85−239.64−53.96−206.66−12.25
1995–2000Variation45.2114.58−147.83−26.44−113.37−0.82
2000–2005Variation−6.2312.46−4.700.85170−0.32
2005–2010Variation24.98−17.77−100.18−125.6591.434.87
2010–2015Variation59.72−68.26−0.88−6.1690.57−0.44
2015–2020Variation37.7955.55−26.61−93.85−89.120.89
1990–2020Variation163.35107.41−519.84−305.21−57.153.43
Table 4. ESV of different categories of ecosystem services from 1990 to 2020 (CNY 100 million).
Table 4. ESV of different categories of ecosystem services from 1990 to 2020 (CNY 100 million).
Primary TypesSecondary Types1990199520002005201020152020
ProvisionFood production577.35594.82613.15613.38615.18622.04638.36
Raw material production435.31437.14435.09436.02428.58423.13424.31
Water supply227.62195.58181.34183.57169.13134.2840.06
RegulationGas regulation1234.421222.321222.381222.101204.981210.431213.34
Climate regulation2817.182736.482699.932698.462640.032642.422601.94
Environmental purification988.97948.54948.34946.57934.32938.41886.60
Hydrological regulation5233.604855.475088.165059.465107.625217.084397.14
SupportSoil formation and retention1494.921485.571470.891473.371442.901425.301417.63
Nutrient cycling145.08145.54146.74146.76145.35146.20147.98
Biodiversity1105.921069.611059.481058.501036.991039.021009.11
CultureAesthetic landscape532.50513.50510.40509.79500.55501.91482.42
Total14,792.8714,204.5714,375.9014,347.9814,225.6314,300.2213,258.89
Table 5. Sensitivity coefficient resulting from the adjustment of the equivalent value coefficient.
Table 5. Sensitivity coefficient resulting from the adjustment of the equivalent value coefficient.
Land-Use Types1990199520002005201020152020
Paddy0.360.350.300.310.290.240.22
Dry land0.060.060.060.060.060.060.06
Forest0.060.060.060.060.060.060.06
Grassland0.250.280.300.300.430.440.65
Water area0.770.810.790.800.780.760.82
Unused land0.010.010.010.010.010.010.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Q.; Wang, L.; Du, G.; Faye, B.; Li, Y.; Li, J.; Liu, W.; Qu, S. Dynamic Variation of Ecosystem Services Value under Land Use/Cover Change in the Black Soil Region of Northeastern China. Int. J. Environ. Res. Public Health 2022, 19, 7533. https://doi.org/10.3390/ijerph19127533

AMA Style

Li Q, Wang L, Du G, Faye B, Li Y, Li J, Liu W, Qu S. Dynamic Variation of Ecosystem Services Value under Land Use/Cover Change in the Black Soil Region of Northeastern China. International Journal of Environmental Research and Public Health. 2022; 19(12):7533. https://doi.org/10.3390/ijerph19127533

Chicago/Turabian Style

Li, Quanfeng, Lu Wang, Guoming Du, Bonoua Faye, Yunkai Li, Jicheng Li, Wei Liu, and Shijin Qu. 2022. "Dynamic Variation of Ecosystem Services Value under Land Use/Cover Change in the Black Soil Region of Northeastern China" International Journal of Environmental Research and Public Health 19, no. 12: 7533. https://doi.org/10.3390/ijerph19127533

APA Style

Li, Q., Wang, L., Du, G., Faye, B., Li, Y., Li, J., Liu, W., & Qu, S. (2022). Dynamic Variation of Ecosystem Services Value under Land Use/Cover Change in the Black Soil Region of Northeastern China. International Journal of Environmental Research and Public Health, 19(12), 7533. https://doi.org/10.3390/ijerph19127533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop