Next Article in Journal
Soil Carbon Storage, Enzymatic Stoichiometry, and Ecosystem Functions in Indian Himalayan Legume-Diversified Pastures
Next Article in Special Issue
Evaluating Performance of Multiple Machine Learning Models for Drought Monitoring: A Case Study of Typical Grassland in Inner Mongolia
Previous Article in Journal
Integrating Urban Energy Resilience in Strategic Urban Planning: Sustainable Energy and Climate Action Plans and Urban Plans in Three Case Studies in Italy
Previous Article in Special Issue
Impact of Drought on Land Productivity and Degradation in the Brazilian Semiarid Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Distribution of the Ecosystem Provisioning Service and Its Correlation with Food Production in the Songhua River Basin, Northeastern China

1
Key Laboratory of Agricultural Water Resources, Hebei Key Laboratory of Agricultural Water-Saving, Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050001, China
2
College of Advanced Agricultural Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
*
Authors to whom correspondence should be addressed.
Land 2024, 13(4), 451; https://doi.org/10.3390/land13040451
Submission received: 1 March 2024 / Revised: 27 March 2024 / Accepted: 31 March 2024 / Published: 2 April 2024

Abstract

:
Provisioning services are essential components of ecosystem services. Food production is usually a driver of land use change, which has the effect on altering the provisioning services of ecosystems. As one of the main areas of food production in China, the provisioning services of the Songhua River Basin (SHRB) should be taken seriously. In view of this, it is urgent to carry out a study on the assessment of provisioning services in the SHRB to provide data support and scientific reference for the optimization of the spatial pattern of land use in the basin, the sustainable development of agriculture, and the formulation of differentiated protection policies. In this study, based on the equivalent factor method for the unit area value and spatial autocorrelation with the Moran’s I, we assessed the provisioning services values (PSV) of the SHRB every ten years during the period of 2000–2020 under different land use types and analyzed the relationships between different PSV and the production of four different food types, including rice, wheat, corn, and soja. The main conclusions are as follows: (1) From 2000 to 2020, the area of paddy fields in the SHRB increased and then decreased, while the area of dry lands continued to increase. The land use transfer matrix showed a significant expansion of paddy fields (+0.55 × 104 km2), shrinkage of grassland (−0.72 × 104 km2), and loss of water body (−0.43 × 104 km2) in the SHRB from 2000 to 2020; (2) The PSV in the SHRB showed an increasing trend from 2000 to 2020, growing by 16.73 × 1010 RMB, with the growth in 2010–2020 being greater than in 2000–2010. The order of increase in each type of PSV was: water supply > food supply > raw material supply; (3) Spatially, the increase in PSV per unit and total PSV in the SHRB was lesser in the center and greater in the east and west. Meanwhile, the spatial distribution of various PSV showed that the value of unit area food supply was higher in the central and eastern plains, while the raw material supply and water supply were higher in the western and eastern hilly areas. (4) In terms of spatial correlation, the distribution of soja production with the total PSV, food supply, raw materials supply, and water supply services values were positively spatially correlated. However, the production of rice, wheat, and corn with the total PSV, food supply, and raw materials supply services values were negatively spatially correlated. Cluster analysis revealed that changing the crop cultivation structure could protect the ecosystem and increase the value of ecosystem services.

1. Introduction

Ecosystem provisioning services, one of the four main categories of ecosystem services, are relevant to human survival and development because they provide ecosystem goods, including the supply of food, raw materials, and water [1,2,3,4]. Food and water are the lifeblood of human society [5]. Although the conversions of natural ecosystems into agricultural and developed areas have led to vast increases in the production of food [6], land use change can cause fragmentation of landscapes [7]. These increase the vulnerability of ecosystems to extreme events and put them at a high risk of disasters [8]. Water is widely acknowledged as a fundamental component of food security. With continuous increases in grain production, China has entered a stage of rapid agricultural development. At the same time, it has caused an increase in pressure on agricultural ecology, soil, water, and other natural resources [9,10]. Natural resources play a key role in human wellbeing, and provisioning services values (PSV) can serve as a useful tool for assessing their importance [11,12]. Quantifying the PSV not only helps to improve public recognition of natural capital and understanding of ecosystem quality and change, but also provides a scientific basis for the remunerative use of ecological resources, related policy formulation, and regional ecological planning [13,14].
Only 135 million hm2 of arable land is available in China, which is a serious shortage of resources. For the economic and social development of China as well as its national security, food security and self-sufficiency are strategic tasks [15]. China’s Black Soil Granary is one of the most famous fertile black soil areas in the world [16,17]. It is a major supplier of commercial grains such as corn and rice in China, accounting for about a fifth of the country’s annual grain production. The Songhua River Basin (SHRB) is an important agricultural region in Northeast China and is among the areas of China’s Black Soil Granary, using only 3.6% of the country’s land area and 3.5% of its total water resources to produce about 18% of the country’s cereal output. The problems of black soil erosion and future hydrological drought events in the SHRB are frequent and seriously endanger the social environment. With the economic and social development, the population in the SHRB is increasing rapidly and the problem of food production and water supply are becoming increasingly prominent. Due to the overexploitation of water resources by irrigated agriculture, the average water storage is significantly lower than pre-2000 levels, despite some recovery in the last 20 years [18]. According to the Statistical Yearbook and Water Resources Bulletin in the SHRB (www.slwr.gov.cn, (accessed on 1 August 2023)), the area of cereals sown increased by 7.75 × 104 km2 and the amount of irrigation water for agriculture increased by 7.47 × 109 m3 over the period 2000–2020. In 2020, the water supply in the SHRB is 44.91 × 109 m3, of which 37.27 × 109 m3 is for agriculture, accounting for over 80% of the water supply. The increase in the regional agricultural areas, agricultural water consumption, and food production have an impact on the amount of water resources available in the region. Food production therefore has a significant impact on regional provisioning services, particularly water supply, and there is a lack of research on the impact of food production on provisioning services in the SHRB. Simultaneously, due to the difference in water and land consumption per unit yield of different crops, while ensuring the growth of grain production, the adjustment of the internal grain planting structure is bound to affect the water and soil resource consumption of China’s grain production [19]. The value of ecosystem services changes with the changing land use patterns in the SHRB. The study of food security evaluation in the region mainly starts from the perspective of arable land resources. In addition, the valuation of ecosystem services in the region has mostly focused on the comprehensive ecosystem service value of a section of the basin or the whole region [20]. Fewer studies have examined the provisioning services and food production of the entire basin.
Therefore, in this study, with analyzing and calculating the temporal and spatial distribution characteristics of food, raw materials, and water supply services in the SHRB, Northern China, and their correlation with food production, it attempts to evaluate the influencing factors on food production in the SHRB from the perspective of the ecosystem provisioning service. Our targets were to (1) analyze the land use and PSV from 2000 to 2020 in the SHRB; (2) assess the spatio-temporal dynamics of unit PSV change and total PSV change from 2000 to 2020 in the municipalities of the SHRB; and (3) reveal the impact of different crops on urban PSV and propose effective measures to improve the provisioning services values. The results provide a theoretical basis for scientific research on land use management and food security in the SHRB, Northern China.

2. Materials and Methods

2.1. Study Area

The SHRB is located between 119°52′–132°31′ E, 41°42′–51°38′ N with a drainage area of 55.68 × 104 km2, comprising mainly 25 municipalities in Jilin and Heilongjiang provinces, China, a small part in Liaoning and Inner Mongolia Autonomous Regions, China (Figure 1). The basin is in the cold temperate zone, which is hot and rainy in summer and cold and dry in the winter. The average annual temperature is 3–5 °C and the average annual precipitation is about 500 mm [21]. The soil type in the basin is diverse, while chernozem and phaeozem are more widely distributed, accounting for 8% and 35% of the total area, respectively. The land is fertile, producing rice, corn, wheat, and soybeans. With 13% of the country’s land area and 8% of the population, it contributes a quarter of the country’s food production.

2.2. Data Sources

The land use data with a spatial resolution of 1 km are obtained from the Resource and Environment Science Data Centre of the Chinese Academy of Sciences (http://www.resdc.cn, (accessed on 1 December 2022)), including three periods of 2000, 2010, and 2020. The land use types were divided into seven types, namely, forestland, grassland, water body, construction land, bare land, paddy field, and dry land. The crop yield and sown area data for each city in the SHRB are from the 2000, 2010, and 2020 statistical yearbooks of provincial and municipal governments, provided by the local statistical bureaus of Jilin Province, Heilong Province, and Inner Mongolia. The 2020-unit price data of grain comes from the Agricultural Market Information System (AMIS), which is an inter-agency platform.

2.3. Calculation of the PSV

According to the national land use classification standard system, the types of land use in the SHRB are classified into seven primary types. An equivalent factor value (D) is the average annual economic value of grain production for 1 hectare of farmland, which is 1/7 of the national average value of grain in the same year, calculated as [22]:
D = 1 7 · S A
where S is the total economic value of annual crops and A is the annual crop sown area.
The amount of economic value per unit area of the ecosystem service value equivalent factor in the SHRB was calculated to be Chinese Yuan (RMB) 167,624.18/km2, 42,634.69/km2, and 933.82/km2 in 2020, 2010, and 2000, respectively.
The ecological service value (ESVij, in RMB/km2/a) of the land corresponding to each raster is accounted for by the formula:
E S V ij = Y ij · D
where Yij is the ESV equivalence factor for the ijth raster corresponding to the land use type, which was determined based on the table of secondary ecosystem service value equivalence factors in China as modified by Xie et al. [23] and Zhang et al. [20]. The value of PSV per unit area in the SHRB was shown in Table 1.

2.4. Moran Index

Spatial autocorrelation is divided into global spatial autocorrelation and local spatial autocorrelation. Global spatial autocorrelation can explore the correlation between variables in the whole study area, while bivariate global spatial autocorrelation analysis can reveal the spatial correlation and significance between different variables in the study area [24]. The bivariate Moran index allows for the assessment of correlation and the degree of difference between the two indicators. The Moran index can be derived from the analysis using GeoDa software 1.20.0 and is therefore used here to assess the spatial correlation between crop acreage and the PSV in the SHRB municipalities. The bivariate Moran index can be calculated between −1.00 and 1.00, with values greater than 0 indicating a good correlation between the two variables. If the calculated Moran index is 0, it indicates that there is no good correlation between the two variables, while values less than 0 indicate a large difference between the two variables.

2.5. Data Analysis

The land use data and municipal boundaries of the SHRB were overlaid using ArcGIS 10.2 to obtain the land use of the municipalities in the SHRB through spatial extraction. GeoDa software 1.20.0 was used for spatial correlation analysis. All data were statistically processed using IBM SPSS Statistics 25, Microsoft Excel 2023, and Origin Lab 2022.

3. Results

3.1. Land Use Change

The land use types in the SHRB are displayed in Table 2, respectively. The highest proportion of forestland is found in the SHRB in 2000, followed by dry land and grassland, while the proportion of water body, construction land, bare land, and paddy field are all below 10%. It is similar in 2020, with the highest proportion of forestland and dry land being at 71.49%. The trend in land use change from 2000 to 2020 is a moderate decrease in grassland and water body and an increase in construction land, bare land, and paddy field, with a decrease in grassland and water body by 1.29% and 0.78%, respectively, from approximately 6.66 × 104 km2 and 1.50 × 104 km2 in 2000 to only 5.94 × 104 km2 and 1.07 × 104 km2 in 2020. Over the same period, construction land and paddy field increased from nearly 1.37 × 104 km2 and 2.60 × 104 km2 in 2000 to 1.66 × 104 km2 and 3.15 × 104 km2 in 2020. In addition, forestland and dry land all showed small changes.
For land use transfer (Table 3), throughout the whole study period the area converted from grassland to forestland, dry land, and bare land in 2000 to 2020 is 163.55, 109.84, and 63.22 km2, respectively. Water body is mainly converted to bare land, with an area of 37.09 km2. Dry land in 2020 is mainly converted from forestland and grassland in 2000, with an area of 200.93 km2 and 109.84 km2, respectively, while the increase in construction land is mainly converted from dry land. This may be driven by population, social, and economic factors. Rural areas have increased the demand for cultivated land, resulting in the reverse ecological process of converting forestland and other ecological land into dry land. With the acceleration of urbanization and the expansion of urban area, the process of converting dry land into construction land has also occurred. The decrease in grassland and the increase in paddy field and dry land may be due to the degradation of grassland caused by overgrazing and anthropogenic mowing, resulting in the conversion of grassland to wetland, which is reclaimed and turned into farmland.

3.2. Variations in the PSV of the SHRB

The provisioning services and its three secondary categories are listed in Table 4. The total PSV of the SHRB increased significantly, from 0.10 × 1010 RMB in 2000 to 16.83 × 1010 RMB in 2020. We can see that the PSV had an increasing trend during 2000–2010 (4.39 × 1010 RMB) and 2010–2020 (12.34 × 1010 RMB). Regarding the three secondary categories, the values of services were all increased for food supply, raw materials supply, and water supply in the period of 2000–2020. Developments in science and technology have increased the productivity of the land and, as a result, crop yields have increased dramatically.
Figure 2 displays the spatial distribution of the PSV per unit in 2000, 2010, and 2020, as well as the total PSV changes in 2000–2010, 2010–2020, and 2000–2020 in the SHRB at the city scale. These distribution patterns varied among the 25 cities in the SHRB (there are 28 cities in total, but 3 of them have very little area within the Songhua River Basin, so they will not be analyzed here). In 2000, 2010, and 2020, the distribution patterns of low PSV per unit were the same, and were concentrated in the Songnen Plain, such as Qiqihar, Baicheng, and Changchun. The high PSV per unit were mainly distributed in the eastern and northwestern mountainous areas of the SHRB. During the three periods, the lowest PSV per unit were found in the Baicheng city, with 44,807.99 RMB/Km2 in 2010 and 167,366.47 RMB/Km2 in 2020, respectively. The highest PSV per unit were found in the Baishan city, with 2753.71 RMB/Km2, 125,155.95 RMB/Km2, and 506,017.25 RMB/Km2, respectively.
In terms of spatial variations during the entire period (2000–2020), the total PSV in all 25 cities showed an increasing trend. There was a lesser increase in the central plain area of the SHRB and a greater increase in the eastern mountains and western mountain–plain transition zone areas. Changes in the total PSV between 2000 and 2010 were relatively small, and while all cities saw increases in the total PSV, only one city in the eastern region of the SHRB saw an increase of more than 0.80 × 1010 RMB. The total PSV in the SHRB from 2010 to 2020 were generally consistent with the trend of increase from 2000 to 2020, with smaller increases in the total PSV in the central plains of the watershed, and larger increases in the eastern mountainous and western mountain–plain transition zones.

3.3. Spatial Distribution of Secondary Categories of the PSV and Their Change

Changes in the value of secondary categories of provisioning services by unit area of the SHRB are shown in Figure 3. The value of food supply per unit area increased significantly from 2000 to 2020, with the value of food supply in the central and eastern regions higher than that in the west. This is mainly due to the flat terrain and fertile soil in the central region, coupled with scientific and technological progress as well as national policy support in guaranteeing agricultural water security, combating soil erosion problems, and protecting black soil resources. During the period 2000–2010, the Harbin and Changchun cities, as the provincial capitals of the Heilongjiang and Jilin provinces, experienced rapid urbanization, which required a large amount of land, and part of the arable land was converted into construction land, so the value of food supply service has declined to some extent. In contrast, the spatial distribution pattern of the value of raw material supply and water supply services per unit area in the cities in the basin was opposite to the distribution of food supply, which was low in the central region and high in the east and west. It is worth noting that the value of water supply service of Qiqihar city in 2020 was −1014.42 RMB. The negative value of water supply service in 2020 may be since the water body of Qiqihar city has shrunk, and more than half of the transferred area of water body has been reclaimed as paddy field and dry land, and the ESV equivalent coefficient of the paddy field is negative. In addition, the application of fertilizers and pesticides on agricultural land will easily increase the pollution load of water sources.

3.4. Correlation between the PSV and Food Production

To clarify whether food production has significantly affected the PSV, a bivariate spatial autocorrelation analysis was conducted using the unit area yield of major crops in the SHRB in 2020 with provisioning services and three secondary categories. Table 5 shows the global bivariate Moran index between the PSV and farmland production of unit area in the SHRB in 2020. The results showed that the Moran index was mostly less than 0, indicating a weak negative spatial correlation. The index of soja unit area production with the total PSV, food supply, raw materials supply, and water supply services values were greater than 0, indicating a weak positive spatial correlation.
However, the global Moran index is unable to estimate the spatial autocorrelation structure. Anselin [25] proposes a local Moran index that can test for local area spatial correlation: bivariate local spatial autocorrelation LISA clustering plot based on a z-test (p = 0.05). Figure 4 shows the cluster analysis of food supply, raw materials supply, and water supply in relation to the unit area production of rice, wheat, corn, and soja for each region of the SHRB. There is no spatial autocorrelation between the production of the four crop types and food supply in the majority of the SHRB. However, in terms of raw material supply, the relationship between production and services for the four crop types in the northwestern part of the basin shows H-L, and the H-L for water supply occurs in the southwestern part of the basin. Therefore, from the point of view of food supply, the production of the four crops does not affect it in most areas of the SHRB. The northwestern part of the basin is not suitable for crop cultivation in terms of raw material supply. And for the southwestern part of the basin, the water supply is not sufficient to support food production. Based on the spatial relationship between the two variables, if the PSV of the ecosystem in some municipalities in the SHRB are to be increased, it can be tentatively concluded that changing the crop cultivation structure can increase the supply service value, e.g., the Hinggan League should reduce the cultivation of all four crops. As an important forest and agricultural area in China, the Hinggan League is rich in forest and grassland resources, and its natural resource assets are mainly ecological assets [26]. Improvements in technology and agricultural practices have led to increased crop yields, while crop selection is also important. Other studies have also shown that food production varied considerably with future land use choices and management, and that maintaining food and water quantities can be difficult. Future changes may lead to significant trade-offs and even lose-lose outcomes across the range of situations studied [27].

4. Discussion

4.1. Impacts of Climate Change on the PSV and Food Supply

Climate change affects ecosystem service provisioning mainly through direct or indirect. Climate change directly affects services such as water supply, food supply, carbon sinks, net primary productivity (NPP), and biodiversity distribution [28,29]. Millennium Ecosystem Assessment recognizes that climate change is causing a downward trend in the level of provisioning of most ecosystem services. It has been suggested that climate change affects global food production, and that global warming has made precipitation in the semiarid subtropical regions increasingly scarce, which has led to a decline in food production capacity in the region [30]. Briner [31] et al. found that climate change has a greater impact on the function of food production in high-altitude areas, while the middle- and low-altitude areas are mainly affected by land use changes driven by economic interests. One study analyzed the impact of climate change on crop production and income in Japan and found that climate change reduced rice production and farmers’ income in the eastern and northern parts of Japan, but increased rice production and farmers’ income in the western part [32].
It has been shown that climate change has a significant impact on ecosystem water conservation and soil retention functions and is a major driver of water yield and soil retention enhancement [33]. However, due to global warming, the frequency and intensity of extreme weather events are increasing globally and regionally, which will profoundly affect the structure and functioning of ecosystems in various places, thus negatively affecting ecosystem services [28]. Studies on the impacts of climate extremes on ecosystem services have mostly focused on the impacts of extreme precipitation, drought, and high temperatures. For example, although extreme precipitation leads to an increase in water supply services [34], it reduces food supply functions due to soil erosion and soil degradation. In 2022, the SHRB also encountered heavy precipitation events, and many farmlands were heavily waterlogged, severely affecting food production. To cope with extreme weather, research and management of soil erosion management can be carried out on a watershed basis, and research on soil erosion management techniques and mechanisms in black soil areas can be strengthened. On a watershed basis, dynamic monitoring of soil erosion and erosion gully management can be systematically carried out, the construction of small watersheds and irrigation areas in black soil areas can be accelerated, small-scale farmland water conservancy infrastructure can be improved, and the configuration of farmland protection forest networks can be optimized.

4.2. Ecological Effects of Different Land Use Types

The ecological effects of land use change in the SHRB are of two kinds: enhancing and reducing ecosystem services, thus showing that different types of land use change have different ecological effects [25,26]. As arable land is converted into forestland and grassland, the main ecological improvements are achieved, while forestland, grassland, and water bodies are negatively affected, mainly in the form of forestland and grassland into arable land and water bodies into bare land. It is worth noting that grassland degradation is severe in the SHRB due to the massive replacement of grassland by dry land and bare land. There is a growing conflict between the supply of agricultural products and the shortage of pasture, both of which result from the mismanagement of overgrazing and reclaimed grasslands [35]. Thus, it is necessary to improve the black soil conservation system integration program for the black soil farming system, soil-type characteristics, and environmental obstacle factors in different regions. For example, the eastern and central parts of the basin should pay full attention to the problem of grassland degradation caused by overgrazing and reclamation. Unplanned reclamation and encroachment of grasslands and water bodies to grow crops, especially rice, will not only reduce vegetation cover, but also negatively affect the value of ecosystem services in the central part of the Songhua River Basin. More action is needed to optimize crop cultivation structures and to protect the degraded grassland and water body. For the SHRB, different types of land use, especially cropland and forestland, need to be spatially optimized for food-friendly land with full consideration of their ecosystem service value. A focus can be put on the healthy evolution process of cropland systems in the SHRB, adopting the strategy of protecting the health of cropland systems through the principles of “fallow”, “rotation”, “conservation”, “retreat”, and “comprehensive management”, to provide an important guarantee for the sustainable utilization of cropland.

4.3. Impact of Cropping Systems on the PSV and Food Supply

Due to population growth and economic development, urbanization intensified, putting great pressure on land [36,37]. In the 13th Five-Year Plan, as part of its efforts to encourage the development of two urban agglomerations in the northeast, Harbin–Changchun and south-central Liaoning, the state is strongly supporting them [38], resulting in a serious occupation of arable land for construction around large- and medium-sized cities such as Harbin and Changchun. A combination of precipitation variability and vegetation degradation has been causing serious problems in recent years [39]. There is a decrease in grassland quality, a decrease in carrying capacity, and a decrease in soil and water conservation. Artificial grass cutting transforms grassland into wetland, which is then reclaimed as paddy fields or drylands. As marsh wetlands disappear and artificial wetlands increase, it is primarily due to food shortages.
In particular, the results show that the total PSV of the SHRB in the 2000–2020 period increased by 16.73 × 1010 RMB. This is mainly because the direct economic benefits are derived from the physical material goods created by the provisioning services. The centralized distribution of farmland in the central and eastern regions is the main area for food and raw material supply. The value of provisioning services in the mountainous areas in the east and west mainly derives from the raw material supply of forests. The value of food supply of farmland and the raw material supply of forests and grasslands in the western agro-pastoral ecotone are the main sources of the value of provisioning services. At the same time, we found during our fieldwork that human activities in the Songnen Plain strongly affect the degradation of grassland, mainly reclamation, grass cutting, overgrazing, etc. (Figure 5 ). And farmers’ participation in the conservation of biodiversity and its services is crucial for agricultural production. So, farmers can firstly actively participate in publicizing and educating rural grassroots environmental protection laws, regulations, and ecological knowledge, recognize rural ecology, understand the urgency and importance of ecological protection, and enhance their awareness of ecological protection. Secondly, they can protect the farming cultural heritage, by insisting on fewer demolitions and more changes, cautiously cutting down trees, prohibiting digging up mountains and not filling in lakes, respecting the rural texture, and forming a harmonious unity of mountains, water, forests, fields, and lakes. Then, they can use seed banks and traditional farming techniques to promote local crop diversity and guarantee food security. Finally, they can learn traditional techniques that are resource-saving, ecologically conserving, and environmentally friendly, and combine them with modern agronomy, incorporating mechanized, automated, informatized, and intelligent means.

4.4. Limitations

This paper used a land use map with a 1 km resolution remote sensing image. In subsequent studies, higher resolutions can be used to improve the accuracy of the findings. Secondly, in this paper we have chosen three years of data: 2000, 2010, and 2020. In future studies, the conclusions can be made more convincing by lengthening the time scale. Thirdly, while this study has explored the spatial and temporal variation in PSV and the relationship between PSV and food production, the mechanisms have not been studied in detail. Finally, this study only examines the relationship between food production and the PSV in the SHRB from the city scale. In the future, we will consider the relationship between food production and other ecosystem services.

5. Conclusions

This study shows the changes in land use and the PSV during the last 20 years and the correlation between food production and PSV in the SHRB. The results show a significant expansion of paddy fields (+0.55 × 104 km2), shrinkage of grassland (−0.72 × 104 km2), and loss of water body (−0.43 × 104 km2) in the SHRB from 2000 to 2020. At the same time, the total PSV show an increasing trend, increasing by 16.73 × 1010 RMB in 2000–2020. The high values are found in the western and eastern portions of the SHRB for both unit PSV and total PSV. From 2000 to 2020, the value of secondary categories of provisioning services by unit area increases in most cities in the SHRB. A specific analysis of the value of the three provisioning services, in terms of the value of food supply, finds that the value is higher in the central and eastern plains than in the western region. In contrast, the value of raw material supply and water supply are higher in the east and west and lower in the central region. In the spatial autocorrelation relation, from the point of view of food supply, the production of the four crops does not affect it in most areas of the SHRB. The northwestern part of the basin is not suitable for crop cultivation in terms of raw material supply. And for the southwestern part of the basin, the water supply is not sufficient to support food production. Although the method used in this study was relatively limited, it provided a viable reference for scholars studying the ES and the ESV.

Author Contributions

H.Y. and J.C. designed research; Y.Z. performed research; Y.Z. and C.Z. analyzed data; Y.Z. and H.Y. wrote this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Chinese Academy of Sciences Strategic Pioneer Science and Technology Special Project [XDA28020503].

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Acknowledgments

We gratefully acknowledge the anonymous reviewers for spending their valuable time and providing constructive comments on this manuscript.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Baude, M.; Meyer, B.C. Changes in landscape structure and ecosystem services since 1850 analyzed using landscape metrics in two German municipalities. Ecol. Indic. 2023, 152, 110365. [Google Scholar] [CrossRef]
  2. Keesstra, S.D.; Bouma, J.; Wallinga, J.; Tittonell, P.; Smith, P.; Cerdà, A.; Montanarella, L.; Quinton, J.N.; Pachepsky, Y.; van der Putten, W.H.; et al. The significance of soils and soil science towards realization of the United Nations Sustainable Development Goals. SOIL 2016, 2, 111–128. [Google Scholar] [CrossRef]
  3. Yuan, Y.; Wu, S.; Yu, Y.; Tong, G.; Mo, L.; Yan, D.; Li, F. Spatiotemporal interaction between ecosystem services and urbanization: Case study of Nanjing City, China. Ecol. Indic. 2018, 95, 917–929. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Liu, Y.; Zhang, Y.; Liu, Y.; Zhang, G.; Chen, Y. On the spatial relationship between ecosystem services and urbanization: A case study in Wuhan, China. Sci. Total Environ. 2018, 637–638, 780–790. [Google Scholar] [CrossRef] [PubMed]
  5. Li, J.; Yu, Y.; Wang, X.; Zhou, Z. System dynamic relationship between service water and food: Case study at Jinghe River Basin. J. Clean. Prod. 2022, 330, 129794. [Google Scholar] [CrossRef]
  6. Lawler, J.J.; Lewis, D.J.; Nelson, E.; Plantinga, A.J.; Polasky, S.; Withey, J.C.; Helmers, D.P.; Martinuzzi, S.; Pennington, D.; Radeloff, V.C. Projected land-use change impacts on ecosystem services in the United States. Proc. Natl. Acad. Sci. USA 2014, 111, 7492–7497. [Google Scholar] [CrossRef] [PubMed]
  7. Keshtkar, H.; Voigt, W. Potential impacts of climate and landscape fragmentation changes on plant distributions: Coupling multi-temporal satellite imagery with GIS-based cellular automata model. Ecol. Inform. 2016, 32, 145–155. [Google Scholar] [CrossRef]
  8. Zhou, T. New physical science behind climate change: What does IPCC AR6 tell us? Innovation 2021, 2, 100173. [Google Scholar] [CrossRef]
  9. Lu, Y.; Jenkins, A.; Ferrier, R.C.; Bailey, M.; Gordon, I.J.; Song, S.; Huang, J.; Jia, S.; Zhang, F.; Liu, X.; et al. Addressing China’s grand challenge of achieving food security while ensuring environmental sustainability. Sci. Adv. 2015, 1, e1400039. [Google Scholar] [CrossRef]
  10. Zhang, F.; Chen, X.; Vitousek, P. An experiment for the world. Nature 2013, 497, 33. [Google Scholar] [CrossRef]
  11. Liu, M.; Jia, Y.; Zhao, J.; Shen, Y.; Pei, H.; Zhang, H.; Li, Y. Revegetation projects significantly improved ecosystem service values in the agro-pastoral ecotone of northern China in recent 20 years. Sci. Total Environ. 2021, 788, 147756. [Google Scholar] [CrossRef] [PubMed]
  12. Small, N.; Munday, M.; Durance, I. The challenge of valuing ecosystem services that have no material benefits. Glob. Environ. Chang. 2017, 44, 57–67. [Google Scholar] [CrossRef]
  13. Chen, D.; Zhong, L. Review of the value evaluation and realization mechanism of ecosystem services. Chin. J. Agric. Resour. Reg. Plan. 2023, 44, 84–94. [Google Scholar]
  14. Yin, N.; Wang, S.; Liu, Y. Ecosystem service value assessment: Research progress and prospects. Chin. J. Ecol. 2021, 40, 233–244. [Google Scholar]
  15. 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]
  16. Cao, Y.; Hua, L.; Tang, Q.; Liu, L.; Cai, C. Evaluation of monthly-scale soil erosion spatio-temporal dynamics and identification of their driving factors in Northeast China. Ecol. Indic. 2023, 150, 110187. [Google Scholar] [CrossRef]
  17. Wen, D.; Liang, W. Soil Fertility Quality and Agricultural Sustainable Development in the Black Soil Region of Northeast China. Environ. Dev. Sustain. 2001, 3, 31–43. [Google Scholar] [CrossRef]
  18. Liu, M.; Guo, Y.; Zhang, X.; Shen, Y.-J.; Zhang, Y.; Pei, H.; Min, L.; Wang, S.; Shen, Y. China’s Black Soil Granary is increasingly facing extreme hydrological drought threats. Sci. Bull. 2023, 68, 481–484. [Google Scholar] [CrossRef] [PubMed]
  19. Yan, X. Study on Water-Land-Food (WLF) Nexus in Northeast China; Chinese Academy of Agricultural Sciences: Beijing, China, 2020. [Google Scholar]
  20. Zhang, Y.; Qu, J.; Wang, L. Spatial and temporal variation, and distribution uniformity of ESV. J. Fuzhou Univ. (Nat. Sci. Ed.) 2020, 48, 653–660. [Google Scholar]
  21. Wang, S.; Wang, Y.; Ran, L.; Su, T. Climatic and anthropogenic impacts on runoff changes in the Songhua River basin over the last 2015, 56 years (1955–2010), Northeastern China. CATENA 2015, 127, 258–269. [Google Scholar] [CrossRef]
  22. Liu, J.; Liu, M.; Tian, H.; Zhuang, D.; Zhang, Z.; Zhang, W.; Tang, X.; Deng, X. Spatial and temporal patterns of China’s cropland during 1990–2000: An analysis based on Landsat TM data. Remote Sens. Environ. 2005, 98, 442–456. [Google Scholar] [CrossRef]
  23. Xie, G.; Lu, C.; Leng, Y. Ecological assets valuation of the Tibetan Plateau. J. Nat. Resour. 2003, 18, 189–196. [Google Scholar]
  24. Wang, Y.; Ma, J. Effects of land use change on ecosystem services value in Guangxi section of the Pearl River-West River Economic Belt at the county scale. Acta Ecol. Sin. 2020, 40, 7826–7839. [Google Scholar]
  25. Anselin, L. Local Indicators of Spatial Association—LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  26. Bo, W.; Xiao, Y.; Wang, L.; Wang, X.; Ouyang, Z. Assessment of the status of ecological assets and variation of its characteristics: A case study of Hinggan League, Inner Mongolia. Acta Ecol. Sin. 2019, 39, 5425–5432. [Google Scholar]
  27. Qiu, J.; Carpenter, S.R.; Booth, E.G.; Motew, M.; Zipper, S.C.; Kucharik, C.J.; Chen, X.; Loheide, S.P., II; Seifert, J.; Turner, M.G. Scenarios reveal pathways to sustain future ecosystem services in an agricultural landscape. Ecol. Appl. 2018, 28, 119–134. [Google Scholar] [CrossRef] [PubMed]
  28. Grimm, N.B.; Groffman, P.; Staudinger, M. Climate change impacts on ecosystems and ecosystem services in the United States: Process and prospects for sustained assessment. Clim. Chang. 2016, 135, 97–109. [Google Scholar] [CrossRef]
  29. Mina, M.; Bugmann, H.; Cordonnier, T.; Irauschek, F.; Klopcic, M.; Pardos, M.; Cailleret, M. Future ecosystem services from European mountain forests under climate change. J. Appl. Ecol. 2016, 54, 389–401. [Google Scholar] [CrossRef]
  30. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since. Science 1980, 2011, 333. [Google Scholar]
  31. Briner, S.; Elkin, C.; Huber, R. Evaluating the relative impact of climate and economic changes on forest and agricultural ecosystem services in mountain regions. J. Environ. Manag. 2013, 129, 414–422. [Google Scholar] [CrossRef]
  32. Kunimitsu, Y. Regional impacts of long-term climate change on rice production and agricultural income: Evidence from computable general equilibrium analysis. J. Jpn. Soc. Civ. Eng. 2014, 70, 13–19. [Google Scholar] [CrossRef]
  33. Kuglerová, L.; Jyväsjärvi, J.; Ruffing, C.; Muotka, T.; Jonsson, A.; Andersson, E.; Richardson, J.S. Cutting edge: A comparison of contemporary practices of riparian buffer retention around small streams in Canada, Finland, and Sweden. Water Resour. Res. 2020, 56, e2019WR026381. [Google Scholar] [CrossRef]
  34. Chiang, L.; Lin, Y.; Huang, T.; Schmeller, D.; Verburg, P.; Liu, Y.; Ding, T. Simulation of ecosystem service responses to multiple disturbances from an earthquake and several typhoons. Landsc. Urban Plan. 2014, 122, 41–55. [Google Scholar] [CrossRef]
  35. Zhou, J.; Zhang, F.; Xu, Y.; Gao, Y.; Xie, Z. Evaluation of land reclamation and implications of ecological restoration for agro-pastoral ecotone: Case study of Horqin Left Back Banner in China. Chin. Geogr. Sci. 2017, 27, 772–783. [Google Scholar] [CrossRef]
  36. Song, W.; Deng, X. Land-use/land-cover change and ecosystem service provision in China. Sci. Total Environ. 2017, 576, 705–719. [Google Scholar] [CrossRef] [PubMed]
  37. Zhang, Z.; Xia, F.; Yang, D.; Huo, J.; Wang, G.; Chen, H. Spatiotemporal characteristics in ecosystem service value and its interaction with human activities in Xinjiang, China. Ecol. Indic. 2020, 110, 105826. [Google Scholar] [CrossRef]
  38. Mao, D.; Wang, Z.; Wu, J.; Wu, B.; Zeng, Y.; Song, K.; Yi, K.; Luo, L. China’s wetlands loss to urban expansion. Land Degrad. Dev. 2018, 29, 2644–2657. [Google Scholar] [CrossRef]
  39. Piao, S.; Ciais, P.; Huang, Y.; Shen, Z.; Peng, S.; Li, J.; Zhou, L.; Liu, H.; Ma, Y.; Ding, Y.; et al. The impacts of climate change on water resources and agriculture in China. Nature 2010, 467, 43–51. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Land 13 00451 g001
Figure 2. The PSV per unit and total PSV change in the SHRB from 2000 to 2020.
Figure 2. The PSV per unit and total PSV change in the SHRB from 2000 to 2020.
Land 13 00451 g002
Figure 3. The values of the unit area of food supply (ac), raw materials supply (df), and water supply (gi) in the SHRB from 2000 to 2020.
Figure 3. The values of the unit area of food supply (ac), raw materials supply (df), and water supply (gi) in the SHRB from 2000 to 2020.
Land 13 00451 g003
Figure 4. Bivariate LISA distribution map in the SHRB at the city scale.
Figure 4. Bivariate LISA distribution map in the SHRB at the city scale.
Land 13 00451 g004
Figure 5. Anthropogenic grass cutting in the SHRB has resulted in grassland converting to bare land, and bare land converting to farmland.
Figure 5. Anthropogenic grass cutting in the SHRB has resulted in grassland converting to bare land, and bare land converting to farmland.
Land 13 00451 g005
Table 1. The final value of ecosystem services per unit area in the SHRB.
Table 1. The final value of ecosystem services per unit area in the SHRB.
Land Use TypeForestlandGrasslandWater BodyConstruction LandBare LandPaddy FieldDry Land
Ecosystem provisioning serviceFood supply0.670.220.530.020.001.360.85
Raw materials supply0.350.330.150.030.000.090.40
Water supply2.200.186.250.080.00−2.630.02
Table 2. Land use change in the SHRB (2000–2020).
Table 2. Land use change in the SHRB (2000–2020).
Year200020102020
Land Use TypeArea (104 km2)Percent (%)Area (104 km2)Percent (%)Area (104 km2)Percent (%)
Forestland21.8839.4421.8539.4021.5438.90
Grassland6.6612.016.6612.005.9410.72
Water body1.502.711.492.691.071.93
Construction land1.372.471.402.521.663.00
Bare land3.626.533.566.423.987.18
Paddy field2.604.692.584.653.155.68
Dry land17.8232.1317.9232.3218.0532.59
Table 3. Land use transfer of the SHRB from 2000 to 2020 (unit: km2).
Table 3. Land use transfer of the SHRB from 2000 to 2020 (unit: km2).
20002020
ForestlandGrasslandWater BodyConstruction LandBare LandPaddy FieldDry Land
Forestland1737.32137.6412.0111.6067.5218.11200.93
Grassland163.55303.256.448.5163.2210.33109.84
Water body9.127.0849.323.3437.0913.9928.11
Construction land9.806.192.3627.613.2511.2776.63
Bare land45.7038.5411.796.82151.1522.8384.84
Paddy field18.133.505.0215.076.70135.6575.73
Dry land170.1797.5319.7893.0568.51102.321228.33
Table 4. PSV change in the SHRB from 2000 to 2020.
Table 4. PSV change in the SHRB from 2000 to 2020.
Ecosystem ServicesPSV (1010 RMB)PSV Change (1010 RMB)
2000201020202000–20102010–20202000–2020
Provisioning services0.104.4916.834.3912.3416.73
Food supply0.031.526.021.494.505.99
Raw materials supply0.020.742.870.722.132.85
Water supply0.052.237.942.185.717.89
Table 5. Bivariate Moran index for the SHRB in 2020.
Table 5. Bivariate Moran index for the SHRB in 2020.
Major CropsRiceWheatCornSoja
Total PSV−0.016−0.028−0.0860.091
Food supply−0.093−0.014−0.0870.024
Raw materials supply−0.032−0.097−0.0880.066
Water supply0.006−0.02−0.0640.084
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Yang, H.; Zhu, C.; Cao, J. Spatial and Temporal Distribution of the Ecosystem Provisioning Service and Its Correlation with Food Production in the Songhua River Basin, Northeastern China. Land 2024, 13, 451. https://doi.org/10.3390/land13040451

AMA Style

Zhao Y, Yang H, Zhu C, Cao J. Spatial and Temporal Distribution of the Ecosystem Provisioning Service and Its Correlation with Food Production in the Songhua River Basin, Northeastern China. Land. 2024; 13(4):451. https://doi.org/10.3390/land13040451

Chicago/Turabian Style

Zhao, Yuhan, Hui Yang, Chunyu Zhu, and Jiansheng Cao. 2024. "Spatial and Temporal Distribution of the Ecosystem Provisioning Service and Its Correlation with Food Production in the Songhua River Basin, Northeastern China" Land 13, no. 4: 451. https://doi.org/10.3390/land13040451

APA Style

Zhao, Y., Yang, H., Zhu, C., & Cao, J. (2024). Spatial and Temporal Distribution of the Ecosystem Provisioning Service and Its Correlation with Food Production in the Songhua River Basin, Northeastern China. Land, 13(4), 451. https://doi.org/10.3390/land13040451

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