Next Article in Journal
Different Effect of Cloud Seeding on Three Dam Basins, Korea
Previous Article in Journal
Numerical Simulations of Swirling Water Jet Atomization: A Mesh Convergence Study
Previous Article in Special Issue
Sustainability Analysis of the Water Environment Carrying Capacity of Harbin City Based on an Optimized Set Pair Analysis Posture-Deviation Coefficient Method Evaluation Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of the Evolutionary Characteristics and Influencing Factors of Land and Water Resource Systems in Major Grain-Producing Areas

1
College of Science, Northeast Agricultural University, Harbin 150030, China
2
School of Water Conservancy & Civil Engineering, Northeast Agricultural University, Harbin 150030, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(14), 2553; https://doi.org/10.3390/w15142553
Submission received: 15 June 2023 / Revised: 6 July 2023 / Accepted: 11 July 2023 / Published: 12 July 2023

Abstract

:
In the process of rapid advancement of agricultural production, the dynamic evolution characteristics of land and water resources in the main grain-producing areas and the influencing factors are less studied. This study takes Heilongjiang Province, the main grain-producing area in China, as an example, constructs an index system from three functions: production, life, and ecology, uses information entropy to determine the weights and importance of each index, uses the rate of change of the index to determine the basic data of the dynamic development of the system, combines the weights to determine the dynamic evolution characteristics, and compares and analyzes them with the static evolution characteristics determined by the actual data. The results showed that there were differences in the important indicators under different conditions, and the important indicators under static and dynamic conditions were the proportion of the tertiary industry to GDP(A7) and GDP per unit area(A8), with importance weights of 7.45% and 8.0%, respectively. The static evolution index of the land and water resource system increased slowly from 0.16 to 0.91, while the dynamic evolution index fluctuated and declined from 0.58 to 0.34, indicating that the ability of the land and water resource system in the study area to maximize comprehensive benefits is constantly weakening. Managers can pre-control the development speed of the important indicators under dynamic conditions and promote the sustainable development of the land and water resource system.

1. Introduction

Among the 17 sustainable development goals (SDGs) proposed by the United Nations in 2015, agriculture and food systems occupy eight [1]. It is evident that agriculture has become a key link in the sustainable development of human society. With the substantial increase in grain production in China, problems such as low utilization of water resources [2], serious pollution of land and water resources [3], and food safety caused by agricultural environmental pollution have arisen in agricultural production [4], so there is a need for rational planning and management of soil and water resource systems in major grain-producing areas.
As a collection of natural resources, land and water resources are essential resources for agricultural production. To improve the efficiency of agricultural land and water resource utilization to guarantee food security, scholars are currently conducting scientific research on land and water resource development and utilization [5], optimal management of land and water resources [6], and risk prevention [7], mainly from the actual data of land and water resource systems.
In terms of the development and utilization of land and water resources, the construction of an indicator system is still the main focus, and research is being conducted on the effective evaluation of land and water resources. Cheng et al. objectively analyzed the carrying capacity of land and water resources as well as the main influencing factors in the main grain-producing areas by using a cloud model that converts between qualitative and quantitative analyses [8,9]. Tan et al. assessed the carrying capacity of land and water resources in all Chinese provinces from 2009–2017 by constructing a DPSIR model [10]. Wen et al. conducted quantitative research on the spatiotemporal dynamic changes of the carrying capacity of land and water resources in grain-producing areas from four aspects based on the entropy weight-TOPSIS model [11]. Wang et al. conducted a comprehensive assessment of regional resource sustainability development by applying the newly established relative resource carrying capacity method, which considers population size and economic scale as resource carrying capacity [12]. He and Wang, based on ecological footprint theory, constructed the evaluation framework of the trend of soil and water resource carrying capacity changes from two perspectives on soil and water resource support and pressure [13].
In terms of optimizing the management of land and water resources, which started in the 1940s, Masse’s reservoir optimization scheduling theory was developed [14]. Singh et al.’s linear programming model was used to manage land and water resources to achieve the maximization of output and revenue [15]. Subsequently, due to the multiple uncertainties in agricultural production, mathematical models, such as dynamic programming, stochastic programming, and fuzzy mathematics, have been widely applied to the management of land and water resources [16]. After the continuous improvement of information technology and remote sensing technology [17], these technologies have also been widely applied to the optimization management of land and water resources. With the development of computer technology and in-depth interdisciplinary research, Koven et al. combined the atmospheric simulation model and the ecological simulation model to explore the impact of climate change on water quantity, water quality, and ecology [18]. Yan et al. explored the matching pattern of crop and agricultural land and water resources under future climate change conditions [19]. Li et al. conducted a study on the optimal allocation of agricultural land and water resources in terms of agricultural water resources, energy, and food production under changing environmental conditions [20]. These studies have gradually brought the optimal management of land and water resources to a more mature stage.
As for the risk analysis of land and water resources, since the development and utilization of land and water resources have negative impacts on the environment and the vulnerability of land and water resources themselves [21], the risk analysis of land and water resources aims to quantify and evaluate the possible impacts and losses of the utilization of land and water resources on the environment [22]. In the assessment of water environments, the health-risk assessment model recommended by the US Environmental Protection Agency (EPA) is mainly adopted [23]. The sampling distribution theory and Bayesian inference theory are often used to address the uncertainty issues in risk assessment [24,25]. In terms of the risk assessment of the water quantity population, mathematical statistics and simulation methods, such as system dynamics, are mainly used [26,27]. In the risk assessment of land use, cloud models and landscape models combining fuzziness and uncertainty are mainly used [28]. In the risk analysis of land pollution, quantitative evaluation methods based on ecological footprint, mathematical models, and value accounting [29,30] are used with remote sensing, MapGIS, and other technologies for comprehensive evaluation and management [31,32]. Therefore, the rapidly changing scientific and technological advancements provide favorable tools for scientific research on land and water resources.
Overall, in the evaluation and management of land and water resources, research on environmental changes such as global population expansion, temperature rise, and ecological degradation has become the main direction at present, and research methods and perspectives are continuously enriched and improved. However, there is a lack of research on the dynamic evolution of land and water resources under changing environments and the influencing factors, which makes it difficult to realize the advanced planning of optimal management of land and water resources. In this paper, Heilongjiang province, a major grain-producing region in China, was selected as the study area, and the period from 2004 to 2020 was used as the study period, during which agricultural water use and arable land area continued to grow, and grain production increased rapidly. It constructs an indicator system based on the “Production–Life–Ecology” function of land and water resources, determines the indicator weights by using the information entropy theory, which has the ability to measure uncertain information, and determines the dynamic development index of each indicator by the rate of change. We compare the static development indices determined from the raw data of each indicator, and analyze the evolution trends and influencing factors of land and water resource systems in Heilongjiang Province under dynamic and static conditions, and provide new ideas and reference bases for the advance planning and management of regional land and water resources.

2. Research Area

The study area, Heilongjiang Province, is located in an alpine region with a long winter time span and geographical coordinates: 121°11′~135°05′ E and 43°25′~53°33′ N, see Figure 1. As one of the important commodity grain bases in China, its grain production has been growing continuously since 2003 and will increase nearly threefold to 7.54×107 tons in 2020. At the same time, agricultural water consumption and arable land area also continued to grow, increasing by 1.85 times and 1.74 times, respectively, at the peak. Although the water consumption of CNY 10,000 output value fell by 3.99 × 102 m3, it was still nearly four times higher than the national data of the same period, and the water use efficiency was relatively low. The per capita cultivated land area of the province increased from 0.26 hm2 in 2003 to 0.54 hm2. Although the converted fertilizer use started to show a downward trend in 2015, it increased from 1.26 × 106 tons in 2003 to 2.24 × 106 tons in 2020, with an increase of 77.78%. Therefore, it is necessary to conduct research on the evolutionary trends and influencing factors of regional soil and water resource systems.

3. Research Methods

3.1. Constructing Indicator System

The establishment of the indicator system of the land and water resource system is the basis and prerequisite for analyzing the evolution trend, and also the main basis for determining the influencing factors [33]. In this paper, we observe the basic principles of data accessibility, objectivity, and overall operability of the indicator system [34]. The index system is constructed with the objective of analyzing the impact of the land and water resource system on the sustainable development of the region, starting from the “Production–Life–Ecology” functions of reflecting the production capacity, livelihood protection, and ecological maintenance of land and water resources [35,36] (See Table 1). Among them, the production capacity is the ability to support agricultural and non-agricultural production and economic development, which is the fundamental driving force for the change of livelihood and ecological functions; the livelihood protection is the ability to guarantee regional life and social development, which is the fundamental purpose of production and ecological functions; the ecological maintenance is the ability to maintain and protect the ecological environment, which is the important environmental resources guarantee for the development of production and life functions. The disorderly or vicious development of one type of function among the “Production–Life–Ecology” functions will directly affect the coordinated development of the other functions. Only the coordinated and common development of the “Production–Life–Ecology” functions can promote the maximization of the comprehensive benefits of regional land and water resources and the sustainable development of the land and water resource system [37].

3.2. Determination of Indicator Weights

For this paper, the selected data are statistical data, and there are many uncertainties in the process of obtaining statistical data; the advanced theoretical information entropy [38,39], which quantifies the uncertainty information, is used to determine the weights of evaluation indices. Based on the degree of influence of the relative changes of index data on the system as a whole, the key indicators in determining the sustainable development of the land and water resource system are searched for, so that the evaluation of the evolution characteristics of the land and water resource system and the analysis of the main influencing factors can be carried out. The calculation steps are as follows:
Step 1: Standardization of metrics
Suppose the sample set is { x * i j i = 1 , 2 , , q ; j = 1 , 2 , p } , where x * i j is the value of the i indicator in j year, and q , p are the number of indicators and the length of the time series, respectively.
Positive   indicators :   x i , j = x * i j x min ( i ) x max ( i ) x min ( i )
Negative   indicators :   x i , j = x max ( i ) x * i j x max ( i ) x min ( i )
In the above equation, x max ( i ) and x min ( i ) denote the maximum and minimum values of the indicator i , respectively, and x i , j is the standardized index value.
Step 2: Calculate the entropy value of each indicator
According to the information entropy theory, the entropy value of each index is determined.
E i = 1 log p j = 1 p x i j log x i j
x i j = 1 + x i , j j = 1 p ( 1 + x i , j )
In the above equation, E i is the entropy value of indicator i .
Step 3: Calculate the indicator weights
The greater the information entropy value of each indicator, the smaller the effect on the soil and water resource system, based on this to determine the weights.
d i = 1 E i q i = 1 q E i
In the above equation, d i is the weight of indicator i .

3.3. Identification of Trends in the Evolution of Land and Water Resource Systems

During the study period, there were large differences in the development rate of each indicator, which had a large impact on the development trend of the soil and water resource system in the study area. In this paper, the actual data of each indicator are used to determine the static evolution index of the soil and water resource system, which reflects the current status of the development and utilization of the soil and water resource system under a certain time section. The rate of change derived from the actual data determines the dynamic evolution index of the soil and water resource system, reflecting the development rate of each indicator of the soil and water resource system under changing conditions. The details are as follows:
The static evolution index is determined by normalizing the actual data of each indicator and combining the weights of the indicators (see Equation (6)).
The dynamic evolution index is based on the actual discrete data of each indicator, and the development function of the indicator ( f i ( x ) ) is derived using the Fourier series and Gaussian series to fitting functions (see Equation (9)), and the rate of change of the development function is used to determine the development rate function of each indicator ( V i ( x ) ). The time series of the study period 2005–2019 was used as the independent variable to obtain the rate of change of each indicator at different times as the base data under dynamic conditions. Combined with the calculation method of indicator weights, the weights of the indicators under dynamic conditions ( d i ) are obtained, and then the dynamic evolution index of the soil and water resource system is determined (see Equation (7)).
Static
Evolution Index:
S T j = i = 1 q d i x i , j (6)
Dynamic
Evolution Index:
D T j = i = 1 q d i       V i ( x ) x = j (7)
V i ( x ) = d f i ( x ) / d x (8)
Fourier series
fitting:
f i ( x ) = a 0 2 + n = 1 ( a n cos ( n ω x i ( j ) ) + b n sin ( n ω x i ( j ) ) ) (9)
Gaussian series
fitting:
f i ( x ) = n = 1 a n e x i ( j ) b n c n 2 (10)
R 2 = 1 j = 2004 2020 ( x i ( j ) f i ( x ) x = j ) 2 j = 2004 2020 ( x i ( j ) x i ( j ) ¯ ) 2 (11)
R M S E = 1 p j = 2004 2020 ( x i ( j ) f i ( x ) x = j ) 2 (12)
In the above equation, f i ( x ) is the development function of the indicator i , V i ( x ) is the rate of change of indicator i , d   i are the dynamic indicator weights, V i ( x ) x = j  is the base datum of the indicator i in the year j under dynamic conditions, n in the fit function is determined by R 2 , R M S E .

3.4. Data Sources

Based on the actual data of each indicator in the study area from 2004 to 2020, this paper determines the weights of each indicator to analyze the main influencing factors and conducts a comprehensive evaluation of the evolution trend. The data involved in the transportation land are the total length of roads and railway lines, the total investment in environmental management and the total amount of emissions are from the China Statistical Yearbook on Environment (2004–2020), the forest area is from the China Statistical Yearbook (2004–2020), the urban land area is from the China Urban Construction Statistical Yearbook (2005–2020), and other original data are from the Heilongjiang Statistical Yearbook (2004–2020), with some missing data supplemented by interpolation.

4. Results Analysis and Discussion

4.1. Comparative Analysis of Factors Influencing Land and Water Resource Systems under Different Conditions

Determine the importance of each indicator by weight calculation method and screen out the main factors, which will affect the comprehensive benefits of land and water resources. The actual discrete data of each indicator are brought directly into the Equations (1) to (5) to determine the static importance of each indicator ( d i ). The basic data of each indicator under dynamic conditions are brought into Equations (1) to (5) to obtain the dynamic importance of each indicator ( d i ), see Figure 2.
As shown in Figure 2, there are five indicators of the production function of land and water resources that vary greatly in importance under different conditions. They are effective irrigation rate (A1), water consumption per unit of output value (A2), secondary industry water consumption (A4), the proportion of tertiary industry in GDP (A7) and GDP per unit area (A8), and the differences in importance are above 2%, while the differences in importance of other indicators are below 1%. Under different conditions, the proportion of tertiary industry in GDP (A7) is a statically important indicator (indicator importance of 7.45%), the GDP per unit area (A8) is a dynamically important indicator (indicator importance of 8.00%), and under dynamic conditions, the indicator importance of water consumption per unit of output value (A2) also reaches over 7%, which is a very important indicator. It means that the proportion of tertiary industry in GDP (A7) and GDP per unit area (A8) are the fundamental drivers of changes in livelihood and ecological functions of land and water resources in static and dynamic development, and their indicator meanings are the level of non-agricultural production and the ability of the land to support economic development, respectively. Here, there are differences in the production functions of the important indicators under different conditions.
The importance differences of seven indicators in the life function are small under different conditions, and only the importance difference of traffic density (B3) reaches 3.32%. In addition, the difference in importance of population density (B7) is 1.15% under different conditions, and the difference in importance of other indicators is less than 1%. Under different conditions, urban land density (B5) is a statically important indicator (index importance of 6.26%) and traffic density (B3) is a dynamically important indicator (index importance of 6.74%). It indicates that regardless of the development in that condition, concerning the production and ecological functions of land and water resources, their main purpose is to promote the carrying capacity of land and water resources for social development and life, and there is very little difference in important indicators under different conditions.
Among the ecological functions, the proportion of water for the ecological environment (C2) has the largest importance difference of 2.02% under different conditions. In addition, the forest coverage rate (C1), sewage treatment rate (C4), and proportion of environmental investment (C6) differed by more than 1% under different conditions, while the importance of other indicators differed by less than 1%. Under different conditions, forest coverage rate (C1) is a static important indicator (indicator importance is 6.78%), and SO2 emission per unit area (C5) has an importance of 5.12%, which is the main indicator under static conditions. Meanwhile, SO2 emission per unit area (C5) is an important indicator in dynamic conditions (indicator importance of 5.60%), while forest coverage rate (C1) is another major indicator in dynamic development (indicator importance of 4.95%). It indicates that forest coverage rate and SO2 emission per unit area are important indicators of ecology functions regardless of development under any conditions. The meanings of these two indicators, which are the ability to maintain the ecological environment and the degree of ecological risk, directly affect the ability of land and water resources to provide security for production and life.

4.2. Characterization of the Evolution of Land and Water Resource Systems under Different Conditions

The actual data of each indicator in the study area from 2004 to 2020 were substituted into Equations (1) to (6) to determine the static evolution index of the “Production–Life–Ecology” functions of land and water resources (see Figure 3). Substitute the actual data of each indicator into Equations (9) to (12); the fitting functions of each index were selected by R 2 1 and R M S E 0 and implemented by Python, in which water consumption per unit of arable land (A3), reclamation rate (A5), water resources per capita (B1), fertilizer and pesticide film use per land (C3), and sewage treatment rate (C4) were fitted by Gaussian series, and all other indices were fitted by Fourier series. Considering that 2004 and 2020 are the time boundaries of this study, which has uncertainty, this paper only selects the change rate of each indicator in the study area from 2005 to 2019 as the initial data under dynamic conditions, and in combination with Equations (1) to (5) and Equations (7) to (8), the dynamic evolution index of the “Production–Life–Ecology” functions of land and water resources is determined (see Figure 4).

4.2.1. Static Evolution Index of “Production–Life–Ecology” Functions of Land and Water Resources

The static evolution time series diagram of the “Production–Life–Ecology” functions of land and water resources shows that the “Production–Life–Ecology” functions basically show a slow growth trend. The static evolution index of the production function increased continuously from 0.034 in 2004 to 0.37 in 2019, with a short decline in 2020, which is the largest index change among the land and water resource criterion layer, with a range of 0.34. The static evolution index of the life function increased continuously from 0.072 in 2004 to 0.36 in 2020, which is the largest index change among the land and water resource criterion layer, with a range of 0.29. The static evolution index of ecology function has a minimum of 0.025 in 2007 and a maximum of 0.20 in 2018 with some fluctuation, but it is the least variable index among the land and water resource criterion layer, with a range of 0.16.
The share of tertiary sector in GDP (A7), an important indicator of production function, changed by 0.67 during the study period, while the GDP per unit area (A8), as the fourth most important indicator, changed by 2.28. The urban land density (B5), as an important indicator of the life function, changed by 0.32 during the study period, while the amount of water resources per capita (B1), the fourth most important indicator, changed by 2.84. Forest coverage rate (C1) and SO2 emission per unit area (C5), as important indicators of ecology functions, varied by 0.11 and 0.85, respectively, during the study period, while sewage treatment rate (C4), which ranked fifth in terms of importance, varied by 30.79. It can be seen that under static conditions, there is a large difference between indicators that actually varies a lot among those that are important.

4.2.2. Dynamic Evolution Index of the “Production–Life–Ecology” Functions of Land and Water Resources

From the dynamic evolution of the “Production–Life–Ecology” functions of land and water resources, we can see that the “Production–Life–Ecology” functions basically show a fluctuating downward trend. The dynamic evolution index of the production function fluctuates up from 2005 to 0.32 in 2008 and continues to decline to 0.13 in 2018, which is the largest change in the index of the land and water resource criterion layer, with a range of 0.19. The dynamic evolution index of life function fluctuates down from 2005 to 0.15 in 2007, and continues to rise to 0.20 in 2013, and fluctuates down to 0.15 in 2019. It fluctuates more frequently, but it is the least variable of the indices of the land and water resource guideline layer, which is only 0.05. The dynamic evolution index of ecological function fluctuates and decreases from a maximum of 0.17 in 2005 to 0.04 in 2019, showing an overall fluctuating downward trend. It is the larger index change among the land and water resource guideline layers, with a range of 0.12.
The GDP per unit area (A8) and water consumption per unit of output value (A2), important indicators of production function, changed by 4.29 and 4.46, respectively, during the study period. The indicator of secondary industry water consumption (A4) ranked third in importance, and maximum variation reached 7.85. Traffic density (B3), an important indicator of life function, changed by 4.36 during the study period, also the greatest variation indicator. Forest coverage rate (C1) and SO2 emission per unit area (C5), as important indicators of ecological functions, varied by 0.26 and 1.31, respectively, during the study period, while sewage treatment rate (C4) ranked third in importance, and maximum variation reached 16.86. It can be seen that under dynamic conditions, the difference between the actual indicators that change a lot and the important indicators is small, and the important indicators determined by the rate of change can be used to reflect the development trend of the soil and water resource system, as well as the planning and management of land and water resources by regulating the important indicators.

4.3. Comparative Analysis of the Evolutionary Characteristics of the “Production–Life–Ecology” Functions of Land and Water Resources under Different Conditions

Under static conditions, the evolution index of the “Production–Life–Ecology” functions of land and water resources reflects the current situation of the development and utilization of land and water resources under a certain time section, so it basically shows a slow growth trend. The large differences between the indicators with large changes in the actual situation and the important indicators assessed in this paper may be due to the fact that the static evolution index of land and water resources is directly influenced by policies and objective conditions, and it is difficult to reflect the differences brought about by the evolution of each indicator and the intensity of action over time, and it is difficult to reflect the differences in the indicators and the intensity of action over time. Under the dynamic conditions, the dynamic evolution index of the “Production–Life–Ecology” functions of land and water resources has obvious fluctuation characteristics, and the overall decline of the dynamic evolution index indicates that the dynamic evolution index, which is determined by the rate of change of indicators, can better reflect the uncertainty in the development and utilization of land and water resources. The overall decreasing trend of the dynamic evolution index also indicates that the ability of the study area to obtain the maximum comprehensive benefits of land and water resources is gradually weakening. Moreover, the difference between the indicators with large data changes and the important indicators during the study period is small, which indicates that managers can use important indicators to make timely adjustments to the planning and management programs of land and water resources.
In determining the important indicators, the differences between the important indicators of livelihood and ecological functions under different conditions are small, and the important indicators with large differences mainly come from the indicators of production functions. And it is clear from the data of each indicator during the study period that the study area has good conditions of agricultural production, abundant water, and forest resources. Therefore, managers should pay attention to the development of non-agricultural industries and promote regional economic development mainly by controlling the rate of development of GDP per unit area (A8), of which GDP per unit area (A8) is an important indicator of the production function. And rational planning of urban land and sewage treatment should be carried out to ensure high-quality living conditions and a livable living environment for residents [40]. While raising awareness of ecological protection, we strengthen ecological security [41], promote the coordinated development of soil and water resource systems, and obtain the best comprehensive benefits of land and water resources in the study area.

5. Conclusions

This paper takes the main grain-producing areas as an example, establishes the indicator system from the “Production–Life–Ecology” functions of land and water resources, uses information entropy to determine the importance of the indicators, conducts a comparative analysis of the actual discrete data and dynamic data of each indicator for the evolution trend, and obtains the following main conclusions:
(1)
Under different conditions, the important indices of “Production–Life–Ecology” functions are different. Among production and ecological functions, there are many indicators with a 2% or more difference in indicator importance, while there is only one indicator with great difference in the index importance of life function. The managers should pay attention to the ability of production and economic development ability of land and water resources, as well as the ecological environmental protection ability of land and water resources.
(2)
The static evolution indices of the “Production–Life–Ecology” functions of land and water resources basically show a slow growth trend. During the study period, there are differences between the indicators with large data changes and important indicators, so managers can promote the sustainable development of regional land and water resources from four important indicators: the proportion of tertiary industry to GDP, urban land density, forest coverage, and SO2 emissions per unit area.
(3)
The dynamic evolution index of the “Production–Life–Ecology” functions of land and water resources basically shows a fluctuating downward trend. The difference between indicators with large changes and important indicators during the study period is small, so that managers can adjust the evolution mechanism of the land and water resource system in terms of important indicators to produce the best benefit to regional land and water resources.

Author Contributions

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

Funding

This research was funded by the National Natural Science Foundation of China under grant No: 52179007, and the Natural Science Foundation of Heilongjiang Province under grant No. LH2022E009.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Acknowledgments

The authors thank the National Natural Science Foundation of China under grant No: 52179007 and the Natural Science Foundation of Heilongjiang Province under grant No. LH2022E009.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Z.Q.; Guo, J.Y.; You, X.Y. A study on spatio-temporal coordination and driving forces of urban land and water resources utilization efficiency in the Yangtze River Economic Belt. J. Water Clim. Chang. 2023, 14, 272–288. [Google Scholar] [CrossRef]
  2. Liu, W.; Ma, Z.; Lei, B. Spatiotemporal Distribution of Irrigation Water Use Efficiency from the Perspective of Water Footprints in Heilongjiang Province. Water 2022, 14, 1232. [Google Scholar] [CrossRef]
  3. Wang, Z.W.; Feng, J.; Fang, M.; Ye, X.Y. Sources, Influencing Factors, and Pollution Process of Inorganic Nitrogen in Shallow Groundwater of a Typical Agricultural Area in Northeast China. Water 2020, 12, 3292. [Google Scholar] [CrossRef]
  4. Sun, F.F.; Dai, Y.; Yu, X.H. Air pollution, food production and food security: A review from the perspective of food system. J. Inst. Agric. 2017, 16, 2945–2962. (In Chinese) [Google Scholar] [CrossRef]
  5. Li, M.; Singh, V.P.; Fu, Q.; Liu, D.; Yang, G.Q. Efficient allocation of agricultural land and water resources for soil environment protection using a mixed optimization-simulation approach under uncertainty. Geoderma 2019, 353, 55–69. [Google Scholar] [CrossRef]
  6. Jiang, Q.X.; He, X.L.; Wang, Z.L.; Wu, Y.X.; Liao, H.Y. Water resources optimal allocation model based on interval multi-stage stochastic planning and its application. Adv. Sci. Tech. Water 2022, 42, 1–7. (In Chinese) [Google Scholar]
  7. Zhao, Z.Y.; Li, W.C.; Zhang, Y.Z.; Wang, X. Evaluation of agricultural land and water resources security in Ningxia based on DPSIR model. J. Zhejiang Agric. 2017, 29, 1336–1346. (In Chinese) [Google Scholar]
  8. Fu, Q.; Liu, Y.; Li, T.X.; Cui, S.; Liu, D.; Cheng, K. Analysis of Water Utilization in Grain Production from Water Footprint Perspective in Heilongjiang Province. Trans. Chin. Soc. Agric. Mach. 2017, 48, 184–192. (In Chinese) [Google Scholar]
  9. Ren, Y.T.; Xu, D.Y.; Cheng, K. Risk Assessment of Ecological Security of Regional Land and water resources Based on Dynamic Perspective. North. Hortic. 2018, 36, 172–178. (In Chinese) [Google Scholar] [CrossRef]
  10. Tan, C.C.; Peng, Q.H.; Ding, T.; Zhou, Z.X. Regional Assessment of Land and Water Carrying Capacity and Utilization Efficiency in China. Sustainability 2021, 13, 9183. [Google Scholar] [CrossRef]
  11. Wen, Q.; Sun, J.T.; Fan, L.Y.; Li, X.W.; Li, Q.S. Evaluation of Agricultural Soil-Water Resources Carrying Capacity(ASWCC)and Relation Analysis Based on Entropy Weight TOPSIS in Henan Province. Soil Water Conserv. 2022, 29, 333–338. (In Chinese) [Google Scholar]
  12. Wang, F.; Wang, C.; Ye, Y.; Wen, B. An Improved Method for Evaluating Regional Resource Carrying Capacities: A Case Study of the Tarim River Basinin Arid China. Pol. J. Environ. Stud. 2019, 28, 2415–2428. [Google Scholar] [CrossRef] [PubMed]
  13. He, Y.H.; Wang, Z.R. Water-land resource carrying capacity in China: Changing trends, main driving forces, and implications. J. Clean. Prod. 2022, 331, 130003. [Google Scholar] [CrossRef]
  14. You, X.Y.; Xie, X.M.; Sun, S.J.; Wang, H. Review on the present situation and future prospect of water resources deployment models in China. J. China Inst. Water Resour. Hydropower 2004, 2, 131–140. (In Chinese) [Google Scholar]
  15. Singh, A. Optimal allocation of water and land resources for maximizing the farm income and minimizing the irrigation-induced environmental problems. Stoch. Environ. Res. Risk Assess. 2017, 31, 1147–1154. [Google Scholar] [CrossRef]
  16. Baalousha, H.M.; Tawabini, B.; Seers, T.D. Fuzzy or Non-Fuzzy? A Comparison between Fuzzy Logic-Based Vulnerability Mapping and DRASTIC Approach Using a Numerical Model. A Case Study from Qatar. Water 2021, 13, 1288. [Google Scholar] [CrossRef]
  17. Xu, X.Y.; Li, J.; Tolson, B.A. Progress in integrating remote sensing data and hydrologic modeling. Prog. Phys. Geog. 2014, 38, 464–498. [Google Scholar] [CrossRef] [Green Version]
  18. Koven, C.D.; Lawrence, D.M.; Riley, W.J. Permafrost carbon-climate feedback is sensitive to deep soil carbon decomposability but not deep soil nitrogen dynamics. Proc. Natl. Acad. Sci. USA 2015, 112, 3752–3757. [Google Scholar] [CrossRef]
  19. Yan, Y.J.; Zhou, H.F.; Zhu, W.; Yao, L.L. Matching patterns of land and water resources in rainfed agricultural areas of Central Asia under future climate change scenarios. Agric. Res. Arid. Reg. 2021, 39, 184–193. (In Chinese) [Google Scholar]
  20. Li, M.; Zhao, L.; Zhang, C.L.; Liu, Y.D.C.; Fu, Q. Optimization of agricultural resources in water-energy-food nexus in complex environment: A perspective on multienergy coordination. Energy Convers. Manag. 2022, 258, 115537. [Google Scholar] [CrossRef]
  21. Wang, H.R.; Qian, L.X.; Zhao, Z.Y.; Wang, Y. Theory and assessment method of water resources risk. J. Hydraul. Eng.-ASCE 2019, 50, 980–989. (In Chinese) [Google Scholar]
  22. Jiang, Q.X.; Zhou, Z.M.; Wang, Z.L.; Fu, Q.; Wang, T.; Zhao, Y.Z. Risk assessment and optimization of water resources shortage based on water and land resources coupling. Trans. Chin. Soc. Agric. Eng. 2017, 33, 136–143. (In Chinese) [Google Scholar]
  23. Li, H.Z.; Pei, Y.Y.; You, J. Ecological risk assessment of combined pollution in watersheds. Chin. Sci. Bull. 2019, 64, 3412–3428. (In Chinese) [Google Scholar]
  24. Sun, L.; Sun, J.; Li, Y.P.; Suo, C.; Liu, J.; Gao, P.P. Development of an ensemble Bayesian inference-based copula approach for bivariate risk evaluation of extreme precipitation under climate change. Int. J. Bioclimatol. 2022, 42, 8755–8776. [Google Scholar] [CrossRef]
  25. Dana, L.K.; Curtis, L.S. Bayesian inference in probabilistic risk assessment—The current state of the art. Reliab. Eng. Syst. Saf. 2008, 94, 628–643. [Google Scholar] [CrossRef]
  26. Liu, B.; Jeanne, H.J.; McBean, E.; Li, Y. Risk assessment of hybrid rain harvesting system and other small drinking water supply systems by game theory and fuzzy logic modeling. Sci. Total Environ. 2020, 708, 134436. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, L.B.; Hu, Y.N.; Jin, J.L.; Wu, C.J.; Zhou, Y.L.; Cui, Y. Dynamic prediction of water resources carrying capacity of Chaohu Basin and system optimization regulation based on system dynamics simulation. J. Lake Sci. 2021, 33, 242–254. (In Chinese) [Google Scholar]
  28. Zhu, Z.Y.; Mei, Z.K.; Xu, X.Y.; Feng, Y.Z.; Ren, G.X. Landscape Ecological Risk Assessment Based on Land Use Change in the Yellow River Basin of Shaanxi, China. Int. J. Environ. Res. Public Health 2022, 19, 9547. [Google Scholar] [CrossRef]
  29. Yi, S.Y.; Li, X.N.; Chen, W.P. A Classification System for the Sustainable Management of Contaminated Sites Coupled with Risk Identification and Value Accounting. Int. J. Environ. Res. Public Health 2023, 20, 1470. [Google Scholar] [CrossRef] [PubMed]
  30. Georgios, K.; Dionisios, G.; Nymphodora, P.; Loannis, M. Topsoil pollution as ecological footprint of historical mining activities in Greece. Land. Degrad. Dev. 2018, 29, 2025–2035. [Google Scholar] [CrossRef]
  31. Xiao, H.X.; Pu, S.Y.; He, F.K.; Liu, S.B. Progress of the Application of Remote Sensing Technology in the Soil Pollution Research. Earth Environ. 2020, 48, 622–630. (In Chinese) [Google Scholar]
  32. Ye, M.L.; Yang, M.L.; Liu, C.Y.; Ma, Y.H.; Wang, Q. Application of Hyperspectral Remote Sensing in Monitoring Heavy Metals in Soil. Environ. Monit. Tech. Adm. 2018, 30, 1–5. (In Chinese) [Google Scholar]
  33. Hu, Z.B.; Pang, Y.; Xu, R.C.; Yu, H.; Niu, Y.; Wu, C.G.; Liu, Y. Systematic Evaluation and Influencing Factors Analysis of Water Environmental Carrying Capacity in Taihu Basin, China. Water 2023, 15, 1213. [Google Scholar] [CrossRef]
  34. Cheng, K.; He, K.X.; Sun, N.; Fu, Q. Comprehensive evaluation of eco-environmental resources in the main grain-producing areas of China. Ecol. Inform. 2023, 75, 102059. [Google Scholar] [CrossRef]
  35. Xie, X.T.; Li, X.S. Spatio-temporal evolution characteristics and influencing factors of “production-living-ecological” functions in Henan Province, China. Trans. Chin. Soc. Agric. Eng. 2021, 37, 243–252. (In Chinese) [Google Scholar]
  36. Liu, Q.F.; Wang, X.K.; Zhu, Q.; Song, J.P.; Yan, Q.Y.; Zhao, Y. Coupling relationship of water resources carrying capacity system in Tibet Autonomous Region based on “production-living-ecological” function. J. Nat. Resour. 2023, 38, 1618–1631. (In Chinese) [Google Scholar] [CrossRef]
  37. Nahidul Karim, M.; Daher, B. Evaluating the Potential of a Water-Energy-Food Nexus Approach toward the Sustainable Development of Bangladesh. Water 2021, 13, 366. [Google Scholar] [CrossRef]
  38. Shannon, C.E. The Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, Y.J.; Wang, Y.Z.; Wang, J.P. Objective Attributes Weights Determining Based on Shannon Information Entropy in Hesitant Fuzzy Multiple Attribute Decision Making. Math. Probl. Eng. 2014, 2014, 1–7. [Google Scholar] [CrossRef] [Green Version]
  40. Maria, M.B.; Martinez-Del-Rio, J.; Antolin-Lopez, R. Trade-offs among urban sustainability, pollution and livability in European cities. J. Clean. Prod. 2019, 224, 651–660. [Google Scholar] [CrossRef]
  41. Wang, C.; Yu, C.; Chen, T.; Feng, Z.; Wu, K. Can the establishment of ecological security patterns improve ecological protection? An example of Nanchang, China. Sci. Total Environ. 2020, 740, 140051. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Geographical location of the study area.
Figure 1. Geographical location of the study area.
Water 15 02553 g001
Figure 2. Comparison of the importance of each index under different conditions.
Figure 2. Comparison of the importance of each index under different conditions.
Water 15 02553 g002
Figure 3. Static evolution trend of the “Production–Life–Ecology” functions of land and water resources.
Figure 3. Static evolution trend of the “Production–Life–Ecology” functions of land and water resources.
Water 15 02553 g003
Figure 4. Dynamic evolution trend of the “Production–Life–Ecology” functions of land and water resources.
Figure 4. Dynamic evolution trend of the “Production–Life–Ecology” functions of land and water resources.
Water 15 02553 g004
Table 1. Index system of water and land resources’ “Production–Life–Ecology” functions.
Table 1. Index system of water and land resources’ “Production–Life–Ecology” functions.
Guideline LayerIndicator LayerUnitIndicator DescriptionIndicator MeaningProperty
Production functionEffective irrigation rate (A1)%Effective irrigated area/Cultivated areaAgricultural production level+
Water consumption per unit of output value (A2) m3/yuanTotal water consumption/Total outputThe level of water resources to support the economy
Water consumption per unit of arable land (A3) 104 m3/hm2Agricultural water consumption/Arable land areaThe level of water resources to support agriculture production
Secondary industry water consumption (A4) m3/yuanIndustrial water consumption/Industrial outputLevel of water resources support for non-agricultural production
Reclamation rate (A5) %Arable land area/Total land areaAgricultural production capacity carried by the land+
Grain yield per unit area (A6)t/hm2Total food production/Arable land areaLevel of agricultural production made of land+
The proportion of tertiary industry in GDP (A7) %Tertiary industry output/Total outputNon-agricultural production levels+
GDP per unit area (A8) 104 yuan/hm2Total production value/Total land areaLand supports the level of economic development+
Life functionWater resources per capita (B1) 104 m3/peopleTotal water resources/Total populationThe ability of water resources to support livelihood standards+
Domestic water per capita (B2) 104 m3/peopleDomestic water consumption/Total populationLivelihood standards guaranteed by water resources+
Traffic density (B3)%Transportation land/Total land areaThe ability of land resources to support livelihood standards+
Food production per capita (B4)t/peopleTotal food production/Total populationSocial level guaranteed by land resources +
Urban land density (B5) %Urban land area/Total land areaThe level of social development supported by land+
Urbanization level (B6)%Urban population/Total populationSocial development Level+
Population density (B7) people/hm2Total population/Total land areaThe ability of land to secure a standard of livelihood+
Ecology functionForest coverage rate (C1)%Forest area/Total land areaThe ability of land and water resources to maintain the ecosystem+
Proportion of water for ecological environment (C2) %Ecological water use/Total water useThe ability of water resources to maintain the ecosystem+
Unit area of chemical fertilizer and pesticide film use (C3)t/hm2(Fertilizer use + Pesticide use + Film use)/Cultivated land areaThe level of land resources carrying ecological risks
Sewage treatment rate (C4)%Wastewater treatment volume/Total discharge volumeThe level of water resources carrying ecological risks+
SO2 emissions per unit area (C5)t/km2SO2 Emissions/Town areaEcological risk level
Proportion of environmental investment (C6) %Total investment in environmental governance/GDPThe level of protection of the ecological environment+
Notes: “+” are positive indicators; “−” are negative indicators; they indicate, respectively, relationship between numerical size and promotion of coordinated and sustainable development of land and water resources.
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

Cheng, K.; Fu, Q.; Sun, N.; Wang, Z.; Zhao, Y. Comparative Analysis of the Evolutionary Characteristics and Influencing Factors of Land and Water Resource Systems in Major Grain-Producing Areas. Water 2023, 15, 2553. https://doi.org/10.3390/w15142553

AMA Style

Cheng K, Fu Q, Sun N, Wang Z, Zhao Y. Comparative Analysis of the Evolutionary Characteristics and Influencing Factors of Land and Water Resource Systems in Major Grain-Producing Areas. Water. 2023; 15(14):2553. https://doi.org/10.3390/w15142553

Chicago/Turabian Style

Cheng, Kun, Qiang Fu, Nan Sun, Zixin Wang, and Yuxin Zhao. 2023. "Comparative Analysis of the Evolutionary Characteristics and Influencing Factors of Land and Water Resource Systems in Major Grain-Producing Areas" Water 15, no. 14: 2553. https://doi.org/10.3390/w15142553

APA Style

Cheng, K., Fu, Q., Sun, N., Wang, Z., & Zhao, Y. (2023). Comparative Analysis of the Evolutionary Characteristics and Influencing Factors of Land and Water Resource Systems in Major Grain-Producing Areas. Water, 15(14), 2553. https://doi.org/10.3390/w15142553

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