Next Article in Journal
Study of Urban Flooding Response under Superstandard Conditions
Previous Article in Journal
Modeling of Soil Loss by Water Erosion and Its Impacts on the Cantareira System, Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Synergetic Development of Water and Land Resources Based on a Coupling Coordination Degree Model

1
Yellow River Laboratory, Zhengzhou University, Zhengzhou 450001, China
2
School of Water Conservancy and Civil Engineering, Zhengzhou University, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(8), 1491; https://doi.org/10.3390/w15081491
Submission received: 19 January 2023 / Revised: 18 March 2023 / Accepted: 6 April 2023 / Published: 11 April 2023
(This article belongs to the Section Soil and Water)

Abstract

:
The interaction between water and land resources (WALRs) has been further enhanced with the development of human production activities. Evaluating the synergetic development (SD) level of WALRs is conducive to discovering the weakness of comprehensive utilization of resources and promoting sustainable development. However, previous studies did not clearly elucidate the effects of the synergetic development between WALRs (SD-WALRs). For evaluation methods, the impact of various subsystem development levels on the whole system is often ignored due to its unclear definition. Therefore, in this research, the concept of the SD-WALRs was defined based on synergetic theory. By using the “Driving Force–Pressure–State–Impact–Response” (DPSIR) model, comprehensive evaluation index systems of WALR development were established. The index systems were regarded as efficacy functions of the coupling coordination degree (CCD) model, and the evaluation model of the SD-WALR level was constructed based on it. Taking Luoyang City as an example, using the latest 10 years (2010–2019) with available data as the study period, the results showed that the value of the SD degree increased from 0.609 to 0.789 during the study period, which reached the level of intermediate synergetic development. In general, the development of WALRs showed a positive trend from “high coupling but low synergetic development” to “high synergetic development”.

1. Introduction

Water and land resources (WALRs) are the basic components of nature and the material basis that is essential to human survival and development, which play an irreplaceable role in maintaining ecosystem and society stability [1,2,3]. In recent years, rapid urbanization, population growth and improvement in people’s living standards have made the tension between human demand and resource supply more serious. How to realize sustainable development and the utilization of WALRs has become a key issue that needs to be discussed in depth [4,5]. However, WALRs are not independent; they influence and restrict each other [6]. On the one hand, there is a natural mutual relationship, where the hydrological cycle is an important sub-cycle of the geological cycle, which plays an important role in shaping landforms, weathering rocks and material removal. Land resources are also the main carriers of water storage and change. On the other hand, water is required for almost all land-use activities, such as cultivation, mining and industrial production [7,8]. Meanwhile, the quantity, quality and distribution of water are affected by land-use structure and vegetative cover. Anthropogenic activities based on WALRs have further strengthened their mutual relationship. It is obvious that the quantity and quality of WALRs can directly affect the development and utilization of resources, which may cause the eco-economic benefits generated by production for humans to increase or decrease. Furthermore, the quantity and quality matching degree of WALRs and reciprocity are also remarkably influential factors for eco-economic benefits generation, where the unbalanced or uncoordinated utilization of WALRs will impede sustainable development and even cause ecological degradation [9]. Therefore, realizing a high-level SD-WALRs is the foundation of regional sustainable development and environmental protection, while evaluating the level of the SD-WALRs can regulate the cooperation effect of regional WALRs, discover weaknesses in the comprehensive utilization of resources and promote sustainable development.
The complex interactions of WALRs were investigated in a large number of studies, which can be divided into two categories. One involves taking water resources or land resources as a separate research subject to discuss the impact of changes on the other one. Bronstert et al. [10] summarized how the changes in land utilization, such as agricultural activities, vegetation changes, road construction and urbanization, may affect hydrological processes. The studies of Weber et al. [11] and Farley et al. [12] showed that runoff yield would be increased when grassland expanded or forests shrunk. Xia et al. [13] and Capra et al. [14] investigated the impacts of slope erosion and underlying surface changes on runoff production, respectively. Zhang [15] pointed out that the optimal allocation of agricultural water can ensure the production capacity of cultivated land. These studies focused on the impact of changes in either water or land resources, such as physical and chemical properties and their evolution laws, that were caused by the other resource. Although they had strong pertinence, the mutual relationship between water and land resources could not be analyzed from an overall perspective. The feedback between changes in WALR and the society, economy and ecology could not be fully interpreted either.
The other kind of studies investigated water resources and land resources as two distinct but closely connected subsystems based on systematic theory. These studies analyzed the interaction of WALRs in terms of coupling and coordinated relationships and took the effects of the whole system on the economy, society and eco-environment into account. For example, Tong and Chen [16] coupled land utilization with surface water quality to analyze their bidirectional relationship. Dong et al. [17] analyzed the coupling effects and spatiotemporal differentiation between the urbanization process and WALR exploitation. Lv [18] evaluated the comprehensive bearing capacity of two subsystems through the coupling coordination degree of WALRs. Tan [19], on the basis of Lv’s research, proposed corresponding development modes according to various comprehensive bearing capacities. It can be seen that systematic theory could be used to analyze the interaction of WALRs overall and reflect its intensity by calculating the value of the coupling coordination degree. However, water or land resources are relative subsystems. For their collaboration effects, the mutual matching and orderly development of subsystem elements are crucial. Unbalanced, disorganized development of any subsystem will impede its efficient utilization and sustainable development. In previous studies, the effects analysis of the SD-WALRs was incomplete, which may cause people to not know what state a WALR system is in and how it can achieve high synergetic development. Furthermore, the effects of the subsystems’ development level on the evolution trend of the overall system were neglected.
The effects of the SD-WALRs were defined based on synergetic theory in this study. With full consideration of the driving force and pressure from society and the economy, the state of WALR utilization under their mutual relationship and the measures of environmental protection taken by human beings, the comprehensive evaluation index systems of WALRs were constructed using a DPSIR model, which could assess the value of the subsystems development level. Then, the SD-WALR evaluation model could be constructed by taking the comprehensive evaluation index systems as efficacy functions of the CCD model. The overall research idea is shown in Figure 1.

2. Material and Methods

2.1. Study Area

Luoyang City (112°16′ E–112°37′ E, 34°32′ N–34°45′ N) is located in the west of Henan Province, China. It is a sub-center of the Central Plains Urban Agglomeration and has a high strategic position. The city covers 15,230 km2 of land and has complex terrain. Mountainous and hilly regions account for a large proportion of its area, while plains only account for 13.8%. Land resources that are easy to utilize are still scarce. Meanwhile, water resources in Luoyang are insufficient. Since water resources in this city are mainly supplied by precipitation and groundwater, its distribution is remarkably spatiotemporally uneven and affected by the terrain and climate. Rapid population growth and social economic development have increased the pressure on the exploitation and utilization of WALRs. Moreover, the uneven distribution of water has also greatly restricted the utilization of land resources. The shortage of WALRs and difficult natural conditions for utilization limit the further development of Luoyang City. It is obvious that Luoyang must realize the SD-WALRs to overcome the shortage of resources and break the bottleneck of development. The geographical location of Luoyang and some spatial characteristics of its WALRs are shown in Figure 2. The original data used for this research were from the Henan Statistical Bureau (HSB 2010–2019), Luoyang Statistical Bureau (LSB 2010–2019), Luoyang Water Conservancy Bureau (LWCB 2010–2019) and “Luoyang City land use overall planning”.

2.2. Effects of the SD-WALRs

Synergetic theory was first proposed by the German physicist Haken in 1969. It refers to the process in which various parts of a system or different systems cooperate to generate a new structure [20]. Moreover, the SD of systems emphasizes that the structure of subsystems would be more stable, the function would be more effective and the relationship of subsystems would be highly correlated through interaction. Therefore, the effects of the SD could be summarized into three parts based on the studies of Takayuki [21] and other scholars. First, the elements of subsystems fit each other and develop orderly. Second, the coupling relationship between subsystems is strong and the cooperation is efficient. Last but not least, the whole system adapts to the external eco-environment. It is notable that these three parts are only distinguished for convenience of understanding. In interpreting the effects of the SD, they are inseparable. The effects are shown schematically in Figure 3.
Then, based on the above analysis and the relationship between resources and production activities, the effects of the SD-WALRs can be expressed as follows:
(1)
The elements of the WALR subsystem match each other and develop in an orderly manner. The main goal is that the quantity and quality of WALRs can meet human needs, ensure production needs and maintain ecological stability. In addition, there are no adverse consequences, such as environmental deterioration, resource scarcity or a decline in people’s living standards under the current utilization state. The structure and functions of the whole system are stable and developing.
(2)
The second part says that the interaction between the two subsystems of WALRs is strong and coordinative. The essential feature is that the quantity and distribution of WALRs match. Furthermore, the SD-WALRs is mainly reflected in society and economy, which means that utilization and production of WALRs are in close combination, water use for various types of land can be guaranteed and the water use of various departments is efficient.
(3)
The last part emphasizes the adaptability of WALRs to the eco-economic system’s development and evolution. It refers to people being able to obtain higher eco-economic benefits in the current utilization state. Resources have great potential for exploration to support future development.
In summary, the evaluation method must meet several requirements: it can not only evaluate the subsystems’ development levels and strengths of their mutual relationships accurately but also reflect the adaptability of WALRs to the external eco-environment.

2.3. Evaluation Index of the SD-WALR Subsystems Development State Based on the DPSIR Model

The DPSIR model is an improved version of the PSR (pressure–state–response) conceptual model proposed by the OECD (Organisation for Economic Co-operation and Development). It was first adopted by the European Environment Agency in 1993. Now, the model is often used to analyze environmental issues, sustainable development and other relevant problems because it can provide clear ideas for the structure or function analysis of complex systems and avoid missing decisive information [22,23]. In addition, it can assist in selecting indexes of the evaluation model. Therefore, the DPSIR model can simplify the analysis process regarding complex interactions between various elements of the WALRs system, and scientifically select representative indexes that are pertinent and maneuverable [24]. The model is composed of five parts, namely, D (driving force), P (pressure), S (state), I (impact) and R (response). Considering the role of each part in the evaluation of the SD-WALRs on the basis of its effects given above, the relationship of each part is shown in Figure 4.
Based on the effects of the SD-WALRs, the index systems for evaluating the development level of WALR subsystems can be selected using the DPSIR model.
For the water or land resources subsystem, the main driving factors of system development are the natural and socio-economic driving forces. The occurrence of extreme natural conditions is random and not sequential [22], and synergistic development of WALRs is more concerned with the matching and coordination of water and land resources under the conditions of human activities. Therefore, in the short term, the socio-economic driving force should be the main driving factor used to evaluate the development level of WALRs. Population growth and social and economic development will affect the distribution and utilization pattern of regional WALRs. Therefore, the resident population, gross domestic product (GDP) and urbanization rate were selected to reflect the effect of the social and economy on the subsystems of WALRs. In addition, water pollution and land degradation indirectly affect the availability of regional resources. Water whose quality is worse than grade V will lose its usable value, where the higher the proportion of water worse than grade V, the lower the amount of water resources that will be available. Meanwhile, bare land cannot be directly utilized or requires a large amount of money to be invested in governance before developed. The higher the vegetation coverage, the smaller the degree of land degradation. Therefore, it is necessary to add the proportion of water worse than grade V and the ratio of vegetation coverage as evaluation indexes.
As a direct factor of subsystem changes, pressure comes from two aspects: one is the demand of social and economic development for the quantity and quality of resources, while the other is the competing demands for water and land resources in utilization. The change in land-use structure will reallocate water resources. Meanwhile, water resources also constrain various land-based production activities. For the water resources subsystem, the water resources per capita, total water consumption in the region and volume of sewage were selected to reflect the pressure of its development. For the land resources, the land resources per capita, proportion of regional development and ecological index of cultivated land were chosen to reflect it. Since the competing demands from production activities for WALRs are mainly reflected in the structure of land and water use, corresponding indexes were added in the pressure part of the DPSIR model.
The state is the realistic reflection of the WALR subsystem under the driving force and pressure. The matching degree of resources is directly related to whether humans can realize sustainable utilization, which must be taken into account in the evaluation [25]. Moreover, effective development and utilization of resources can help humans to obtain greater eco-economic benefits, which are also the reflection of the function and structure stability. Therefore, water consumption per unit of industrial added value, irrigation water per unit area, and an urban or rural domestic comprehensive water quota can reflect the efficiency of industrial, agricultural and domestic water use. The proportion of effective irrigation area and per capita housing area can also reflect the land resource utilization status of various departments. The GDP per unit industrial area can also reflect the intensive use of the industrial land area.
The impact is the final result of the system’s development and evolution. The orderly development of WALR subsystems causes the utilization of resources to be more efficient, which causes unit resources to produce greater value and conspicuously improves the living standards of human beings. Obviously, GDP per unit resource should be an important index to reflect the impact. Engel’s coefficient is usually used to reflect the living standards of residents, which was also an indicator used in the model. Moreover, the water resources exploitation ratio and land diversity index can reflect the health degree of system development. They were also selected as evaluation indexes.
To achieve orderly development of the systems, humans must adjust their production intensity and mode, which is the response of human society regarding WALR subsystems [22]. People adapting land for cultivation or reclaiming land from the lake may cause soil erosion. Production activities will also generate sewage and solid wastes. In order to prevent the system’s structure and functions from degrading, people have increased investment to control soil erosion and dispose of sewage and other wastes. Therefore, investment in a water conservancy project, the area of soil erosion controlled, unit land assets investment, and the proportion of sewage and solid wastes disposal were selected to measure the model’s response part.
Based on the above analysis, the evaluation indexes of the development level of WALR subsystems were selected and shown in Table 1 and Table 2.

2.4. Evaluation Model of the SD-WALRs

The level of the SD-WALRs was evaluated using the CCD model, which is based on the concept of the capacity coupling degree in physics [26,27]. If the synergistic development of the whole system is regarded as a multi-objective problem, the contribution of each element in a subsystem to the development is a sub-goal. The effect of each element on the synergistic development can be expressed by setting efficacy coefficients that are evaluated using the DPSIR model, and the efficacy functions can reflect the relationship of various efficacy coefficients. The CCD model sets efficacy functions to reflect the contribution of subsystems to the whole system development and evaluates the positive coupling effects of the subsystems’ interaction using the coupling coordination degree. Then, the results can reflect whether the system has good coupling, as well as harmonious and orderly evolution [28]. The model is very suitable for analyzing the development trend of the WALR system and is very intuitive.
According to the effects, the main determinants of the SD-WALRs under their mutual relationship are the development level of subsystems. This can be evaluated by setting comprehensive evaluation index systems. At the same time, close integration, good linkage and efficient interactions are also characteristics of the SD-WALRs. Therefore, the comprehensive evaluation index constructed above can be used as the efficacy functions of the CCD model. Then, the cooperation and interaction strength of the WALR subsystems in their respective development level can be represented using the coupling coordination degree. In this way, the level of the SD-WALRs can be evaluated. The model-constructing steps are mainly divided into the following two parts:
(1)
A comprehensive evaluation of the development level of WALR subsystems.
Suppose U1 and U2 are the development levels of the two WALR subsystems, then:
U 1 2 = z i j λ i
where zij is the normalized value of index i of the water or land resource subsystem in year j, which can be calculated via the min-max normalization (MMN) method using the original data, and λi is the weight of index i.
Common weighting methods can be generally divided into two categories: one is subjective weighting methods represented by the expert grading (EG) method and the analytic hierarchy process (AHP) model, while the other is objective weighting methods, such as the entropy weight (EW) method and the criteria importance through intercriteria correlation (CRITIC) method. However, it is inevitable that the subjective weighting method could be affected by human cognition and experience, while the objective weighting method sometimes cannot take the importance of indexes in special situations into account because it is completely based on objective data. Therefore, this research made a combination and used the EW method and attribute hierarchical model (AHM) to weight the indexes. The AHM is an improved version of the AHP model, which has the advantages of the AHP model and omits the consistency check when using it [29]. Meanwhile, the EW method weights indicators based on the information entropy of the data [30]. Combining the AHM and EW methods to calculate weights can make the best use of the advantages and bypass the disadvantages of these two. The main steps are as follows:
Suppose λi is the combined weight of the index i of WALRs, which can be calculated using Formula (2):
λ i = ω A H M ω E W i = 1 n ω A H M ω E W
where ωAHM and ωEW are the weights found using the AHM and EW methods, respectively;
ωAHM-i can be obtained using Formulas (3) and (4):
ω A H M i = 2 n ( n 1 ) j = 1 n μ i j
μ i j     = a k a k + 1 ( a i j = k ) 1 a k + 1 ( a i j = 1 k ) 0.5     ( a i j = 1 , i j ) 0   ( a i j = 1 , i = j )
where ωAHM-i is the weight of index i calculated using the AHM; μij is the relative attribute measure; aij is the relative importance of factor i to factor j, and its reciprocal is found according to Satty’s original proposal; and the transformation relationship between aij and μij is shown in Formula (4). k is the ratio of importance pairs determined using Satty’s method [31]. ɑ is the attribute measure conversion parameter, which usually takes the value 1 or 2.
ωEW can be obtained using Formulas (5)–(7):
ω E W i = 1 H j n i = 1 m H j
H j = 1 ln n i = 1 n y i j ln y i j
y i j = ( 1 + z i j ) i = 1 n ( 1 + z i j )
where ωEW-i is the weight of index i found using the EW method, Hj is the entropy value of the indicator i in year j and yij represents the proportion of the index i in year j.
(2)
Evaluation model of the SD-WALRs
Taking U1 and U2 as the efficacy function of the CCD model, the evaluation model of the SD-WALRs is as follows:
C = 2 × U 1 × U 2 U 1 + U 2
T = α U 1 + β U 2
S D = C × T
where C is the coupling degree of the WALR subsystems and ranges from 0 to 1. When C = 0, there is no interaction between the two subsystems, and the extent of coupling between subsystems increases with C. T is the coordination degree of the WALRs, which is used to show the positive coupling extent of the subsystems. ɑ and β are undetermined coefficients of the subsystems’ contributions; because water and land resources are equally important for human development, ɑ and β were both set equal to 0.5. SD is the index of the SD-WALR level, which represents the degree of coupling and the harmonious and orderly evolution of the WALRs.
With the study of system coupling or synergy by [28,32,33], combined with the effects of the SD-WALRs in this paper, the level of the SD-WALRs can be classified according to Table 3.

3. Results

Since the required data has not been fully released since 2020, the latest 10 years that are available (2010~2019) were selected as the research period, which could not only avoid the loss of timeliness of the research but also reflect the changing trend of the development level. The original data of the water and land resources subsystem indexes used in the research were obtained from the website (http://www.henan.gov.cn/ and http://shuili.ly.gov.cn/?m=home&c=Lists&a=index&tid=87) (accessed on 10 May 2022) The weights of indicators found using three methods based on Formulas (2)–(7) are shown in Table 4.
The combined weights had the advantages of subjective and objective weighting methods. They started from the entropy value and modified the weights of some indexes according to the special situation in the study area, such as the structure of the tertiary, industry and policy management. Then, the annual data of the study area were normalized using the MMN method. Based on Formula (1), the results of the development level of the WALRs in Luoyang from 2010 to 2019 are shown in Table 5, and the value change trends of the various DPSIR parts are shown in Figure 5 and Figure 6.
C, T and SD were calculated by Formulas (8)–(10) based on the results of Table 4. According to the classification in Table 3, the evaluation results from 2010 to 2019 are displayed in Table 6 and the trends of the SD-WALRs are given in Figure 7.

4. Discussion

4.1. Evaluation Results of the WALR Subsystems Development

It can be seen from Table 5 and Figure 5 and Figure 6 that the value of the WALR subsystem development level showed a fluctuating upward trend during the study period, which indicated that the internal elements of the WALR subsystem gradually developed from “disorder” to “order”. The values of pressure and state changed greatly. Compared with 2010, the pressure values of the water subsystems decreased from 0.197 to 0.094, and the pressure values of land subsystems decreased from 0.237 to 0.088, respectively, while the value of state increased from 0.041 and 0.011 to 0.179 and 0.114, respectively. The values of impact and response also greatly improved compared with the beginning of the study period. In order to analyze the reasons for their changes, the radar figures of normalized data of the indexes in 2011, 2014 and 2019 are shown in Figure 8, and the trends of some indexes with remarkable changes from 2010 to 2019 are shown in Figure 9.
Combined with the original data and Figure 8 and Figure 9, it can be seen that the increase in the population and urbanization process enhanced the driving force of the system. However, the pressure of the system decreased and the state value increased significantly. This was mainly helped by the effective implementation of various measures for WALR protection and efficient utilization in Luoyang in recent years, such as the control of soil erosion and vigorously promoting the popularization of water-saving production technology and water-saving appliances. The government continuously increased investment in the development, utilization and treatment of WALRs. From 2010 to 2019, the accumulated soil erosion in Luoyang reached 2100 km2 and the ratio of solid waste and sewage treatment was maintained above 95%. The changes in the state indexes, such as the increase in the proportion of effective irrigation area, the decrease in the average irrigation water consumption per area and water consumption industrial added value per CNY 10 thousand, reflected the fact that the utilization efficiency of the WALRs in Luoyang gradually improved, and the development of the WALRs was more in line with the current situation of social and economic development. Figure 9c shows that the area of the industrial land and green land significantly increased compared with the beginning of the study period, where many unused and uncultivated areas were converted to industrial land and green land. The total regional water consumption changed little, the proportion of industrial water consumption decreased moderately, and the proportion of ecological water consumption increased remarkably, which not only reflected the improvement in industrial water-use efficiency but also reflected the fact that the water-use structure of Luoyang was more reasonable.
It is worth noting that the response values of the water resources subsystem fluctuated greatly from 2014 to 2015, and for the land resources subsystem, this phenomenon appeared in 2011, 2012 and 2014. From Figure 9e, the fluctuation was mainly related to the reduction in the soil-erosion-controlled area and the ratio of sewage and waste treatment in these years. Compared with 2013, the soil-erosion-controlled area in 2014 decreased by 550 km2, and the sewage and solid waste treatment ratio was the lowest in the study period. In addition, the improvement of the state value was highly related to the development and utilization of the WALRs, and the matching degree of the total amount and distribution of the WALRs was stable for ten years. In this study, the method used for calculating the matching degree of the WALRs was the same as in Li’s research [9]. As shown in Figure 9d, the Gini coefficient reflecting the spatial matching degree of the WALRs in Luoyang was stable at about 0.1 and the spatial matching degree of the WALRs was classified as good. However, from 2010 to 2014, the system was not in a high-level coordinated state. The spatial and temporal matching of water resources is an important factor that restricts regional sustainable development. However, we cannot judge whether the WALRs reached a high level of SD only from the degree of spatial matching.

4.2. Evaluation Results of the SD-WALRs

Table 6 and Figure 7 show that the coupling degree of the WALRs in Luoyang was stable and close to 1 during the study period. The value of the comprehensive coordination index was low in the early years but it rose greatly later, showing a fluctuating upward trend as a whole. This indicated that there was a strong interaction between water and land resources in Luoyang, and the cooperative relationship between the water and land resources subsystems was constantly strengthened. From 2010 to 2019, the index of the SD-WALRs in Luoyang increased from 0.603 to 0.789, which meant an increase from the primary synergetic development level to intermediate synergetic development. The overall trend showed that the subsystem developed from “high coupling but low synergetic” to “high coupling and synergetic”. Some of the indexes in the impact part, such as the GDP per water unit, GDP per unit land area, land diversity index, development and utilization ratio of WALRs, and Engel’s coefficient improved greatly. This reflected the fact that the WALRs had a positive impact on the outside of the system at this synergetic level.
Lv [18] evaluated the interaction between WALRs in Guiyang City. Because his research did not elucidate the effects of the SD-WALRs, the interaction between the WALR subsystems was evaluated using the coupling coordination degree. According to his method, the indexes of evaluation were selected from the utilization and protection states of the WALR subsystems; then, the level of interaction could be evaluated using the CCD model. By using his method, the trends from his research and this study are shown in Figure 10. It can be clearly seen that the values of the coupling degree and coupling coordination degree evaluated using Lv’s method were smaller than from this research, the coordinated degree values were close and the trends of these three indexes were almost consistent. The reason for the difference is that impact of the interaction of the WALRs on the respective development of subsystems was not considered in the index-selection process. In the analysis of this research, the development level of water or land resources subsystems under their mutual relationship is vital for the SD-WALRs. For a more scientific and accurate evaluation, indexes that can reflect the matching degree of WALRs, the structure of development and utilization, and the impact on socio-economy can be selected for use in the evaluation systems.

5. Conclusions

Based on the synergetic theory, this study comprehensively analyzed the effects of the SD-WALRs. The evaluation index systems reflecting the development level of the WALR subsystem were selected using the DPSIR model. Then, the weights of indexes were assigned through the combination of subjective and objective weighting methods. Finally, the evaluation model of the SD-WALRs was constructed with the help of the CCD model. Taking Luoyang City as an example, the case verification was completed. The evaluation results showed that Luoyang City maintained a strong interaction between the WALRs for many years. In recent years, Luoyang City scientifically formulated land-use planning and continuously strengthened water resource management. With the increase in the area of soil erosion controlled, the promotion of efficient water-saving technologies and the improvement of waste treatment, the level of subsystems development and the cooperation effect were improved. This showed that the development and utilization of WALRs will be more scientific and reasonable. Meanwhile, it was found that the growth in the development level of the land resource subsystem was higher than the water resource subsystem, where water resources will restrict the further development and utilization of land resources. To achieve a higher level of synergetic development of the WALRs, Luoyang must further strengthen the optimization and management of water resources subsystems by optimizing the water-use structure, improving water-use efficiency and increasing investment in water conservancy project construction. However, there are still some problems. The DPSIR model is a simplification of the interaction relationship of various elements in the system, which can describe the overall coordinated evolution trend of the system in a macro way. However, there is insufficient explanation of the interaction relationship, evolution mode and feedback mechanism between various elements in the system. In addition, the evaluation model in this study was only an annual evaluation of the level and change trend. The detailed process of system evolution was not depicted, which still needs further research.

Author Contributions

All authors contributed to this study’s conception and design. Conceptualization: C.L.; methodology: C.L. and W.X.; formal analysis and investigation: M.L.; writing—original draft preparation: C.L. and W.X.; writing—review and editing: C.L., W.X., S.W. and Y.H.; supervision: S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Technologies Research and Development Program of China (2021YFC3000204) and the National Natural Science Foundation of China (NSCF-52079125).

Data Availability Statement

Data of Luoyang City used in research were from the statistical yearbook of Henan Province, the water resources bulletin of Henan Province and Luoyang, etc. The specific data source can be accessed at http://www.henan.gov.cn/ (accessed on 10 May 2022) and http://shuili.ly.gov.cn/?m=home&c=Lists&a=index&tid=87 (accessed on 10 May 2022).

Acknowledgments

This research was funded by the Key Technologies Research and Development Program of China (2021YFC3000204) and the National Natural Science Foundation of China (NSCF-52079125). The authors are grateful to colleagues and friends who shared their meteorological and hydrological data with us. We also thank the reviewers for insightful comments that improved an earlier version of this manuscript.

Conflicts of Interest

The authors declare that they have no competing interests or other interests that might be perceived to influence the results and/or discussion reported in this paper.

References

  1. Liu, Y.; Wu, C. Situation of Land-Water Resources and Analysis of Sustainable Food Security in China. J. Nat. Resour. 2002, 17, 270–275. [Google Scholar]
  2. Dobriyal, P.; Qureshi, A.; Badola, R.; Hussain, S.A. A review of the methods available for estimating soil moisture and its implications for water resource management. J. Hydrol. 2012, 458, 110–117. [Google Scholar] [CrossRef]
  3. Fitton, N.; Alexander, P.; Arnell, N.; Bajzelj, B.; Calvin, K.; Doelman, J.; Gerber, J.S.; Havlik, P.; Hasegawa, T.; Herrero, M.; et al. The vulnerabilities of agricultural land and food production to future water scarcity. Glob. Environ. Chang. Hum. Policy Dimens. 2019, 58, 10. [Google Scholar] [CrossRef]
  4. Qin, L. A Study on the Sustainable Use of Water Resources and Environmental Economy in the Huaihe River Basin. Ph.D. Thesis, Hohai University, Nanjing, China, 2001. [Google Scholar]
  5. Chen, B.; Hanb, M.Y.; Peng, K.; Zhou, S.L.; Shao, L.; Wu, X.F.; Wei, W.D.; Liu, S.Y.; Li, Z.; Li, J.S.; et al. Global land-water nexus: Agricultural land and freshwater use embodied in worldwide supply chains. Sci. Total Environ. 2018, 613, 931–943. [Google Scholar] [CrossRef]
  6. Jiang, Q. Study on Carrying Capacity Evaluation and Dynamic Simulation of Sustainble Utilization of Water and Land Resources in Sanjiang Plain. Ph.D. Thesis, Northeast Agricultural University, Harbin, China, 2011. [Google Scholar]
  7. Lv, Y.; Chen, J.; Zhao, H. Evaluation on Water and Land Economy Coupling of Countries in Yarlungzangbo-Brahmaputra River Basin and Nepal. Water Resour. Prot. 2018, 34, 61–66. [Google Scholar]
  8. Zhang, Y.; Yan, Z.X.; Song, J.X.; Wei, A.L.; Sun, H.T.; Cheng, D.D. Analysis for spatial-temporal matching pattern between water and land resources in Central Asia. Hydrol. Res. 2020, 51, 994–1008. [Google Scholar] [CrossRef]
  9. Li, H.; Zhou, W.; Zhuang, Y.; Ma, C.; Liu, B. Agricultural Water and Soil Resouces Matching Patterns and Carrying Capacity in Yan’ an City. Trans. Chin. Soc. Agric. Eng. 2016, 32, 156–162. [Google Scholar]
  10. Bronstert, A.; Niehoff, D.; Burger, G. Effects of climate and land-use change on storm runoff generation: Present knowledge and modelling capabilities. Hydrol. Process. 2002, 16, 509–529. [Google Scholar] [CrossRef]
  11. Weber, A.; Fohrer, N.; Moller, D. Long-term land use changes in a mesoscale watershed due to socio-economic factors–effects on landscape structures and functions. Ecol. Model. 2001, 140, 125–140. [Google Scholar] [CrossRef]
  12. Farley, K.A.; Jobbagy, E.G.; Jackson, R.B. Effects of afforestation on water yield: A global synthesis with implications for policy. Glob. Chang. Biol. 2005, 11, 1565–1576. [Google Scholar] [CrossRef]
  13. Xia, J.; Qiao, Y.; Song, X.; Ye, A.; Zhang, X. Analysis about Effect Rules of Underlying Surface Change to the Relatiopnship between Rainfall and Runoff in the Chabagou Catchment. Resour. Sci. 2007, 29, 70–76. [Google Scholar]
  14. Capra, A.; Porto, P.; Scicolone, B. Relationships between rainfall characteristics and ephemeral gully erosion in a cultivated catchment in Sicily (Italy). Soil Tillage Res. 2009, 105, 77–87. [Google Scholar] [CrossRef]
  15. Zhang, Y. Research on Optimal Allocation of Land and Water Resources of Cultivated Land in Naoli River Basin. Ph.D. Thesis, Northeastern University, Shenyang, China, 2019. [Google Scholar]
  16. Tong, S.T.Y.; Chen, W.L. Modeling the relationship between land use and surface water quality. J. Environ. Manag. 2002, 66, 377–393. [Google Scholar] [CrossRef]
  17. Dong, W.; Yang, Y.; Zhang, Y. Coupling Effect and Spatio-Temporal Differentiation between Oasis City Development and Water-Land Resources. Resour. Sci. 2013, 35, 1355–1362. [Google Scholar]
  18. Lv, T. Research on Comprehensive Evaluation of Land-Water Resources Carrying Capacity and Its Regulations Mechanism—A Case Study of Guiyang City. Ph.D. Thesis, Zhejiang University, Hangzhou, China, 2015. [Google Scholar]
  19. Tan, K. Study on Coupling and Regulation of Water and Land Resources in Wenshan City of Karst Mountain Area. Master’s Thesis, Yunnan University, Kunming, China, 2020. [Google Scholar]
  20. Chen, S. Research on land-ocean-air synergistic principle in the development of water resources in China for the 21st century. J. Hydraul. Eng. 2006, 37, 1403–1406. [Google Scholar]
  21. Qiao, J. Research on Synergetic Development Mechanism and Planning Strategies of Urban River System Landscape. Ph.D. Thesis, Northeast Forestry University, Harbin, China, 2017. [Google Scholar]
  22. Yu, B.; Lv, C. Application of DPSIR Framework for Analyses of Sustainable Agricultural Development. China Popul. Resour. Environ. 2004, 14, 70–74. [Google Scholar]
  23. Vannevel, R. Using DPSIR and Balances to Support Water Governance. Water 2018, 10, 118. [Google Scholar] [CrossRef] [Green Version]
  24. Gari, S.R.; Newton, A.; Icely, J.D. A review of the application and evolution of the DPSIR framework with an emphasis on coastal social-ecological systems. Ocean Coast. Manag. 2015, 103, 63–77. [Google Scholar] [CrossRef] [Green Version]
  25. Weatherhead, E.K.; Howden, N.J.K. The relationship between land use and surface water resources in the UK. Land Use Policy 2009, 26, S243–S250. [Google Scholar] [CrossRef]
  26. Li, Y.F.; Li, Y.; Zhou, Y.; Shi, Y.L.; Zhu, X.D. Investigation of a coupling model of coordination between urbanization and the environment. J. Environ. Manag. 2012, 98, 127–133. [Google Scholar] [CrossRef]
  27. Fan, Y.P.; Fang, C.L.; Zhang, Q. Coupling coordinated development between social economy and ecological environment in Chinese provincial capital cities-assessment and policy implications. J. Clean. Prod. 2019, 229, 289–298. [Google Scholar] [CrossRef]
  28. Ariken, M.; Zhang, F.; Chan, N.W.; Kung, H.T. Coupling coordination analysis and spatio-temporal heterogeneity between urbanization and eco-environment along the Silk Road Economic Belt in China. Ecol. Indic. 2021, 121, 12. [Google Scholar] [CrossRef]
  29. Chen, Q. Attribute Hierarchical Model-A New Method of Unstructured Decision Making. Acta Sci. Nat. Univ. Pekin. 1998, 34, 12–16. [Google Scholar]
  30. Chen, P.Y. Effects of the entropy weight on TOPSIS. Expert Syst. Appl. 2021, 168, 13. [Google Scholar] [CrossRef]
  31. Agarwal, R.; Garg, P.K. Remote Sensing and GIS Based Groundwater Potential & Recharge Zones Mapping Using Multi-Criteria Decision Making Technique. Water Resour. Manag. 2016, 30, 243–260. [Google Scholar]
  32. Zameer, H.; Yasmeen, H.; Wang, R.; Tao, J.; Malik, M.N. An empirical investigation of the coordinated development of natural resources, financial development and ecological efficiency in China. Resour. Policy 2020, 65, 13. [Google Scholar] [CrossRef]
  33. Jin, C. Study on High Quality Coordinated Development of Ecological Protection and Economy in Central Plains Urban Agglomeration and Its Influencing Factors. Ph.D. Thesis, Liaoning University, Shenyang, China, 2021. [Google Scholar]
Figure 1. Research ideas and steps on the evaluation of the SD-WALRs.
Figure 1. Research ideas and steps on the evaluation of the SD-WALRs.
Water 15 01491 g001
Figure 2. Geographical location of Luoyang (a), land type and distribution of Luoyang (b), isohyet of Luoyang in 2019 (c) and groundwater depth of Luoyang in 2019 (d).
Figure 2. Geographical location of Luoyang (a), land type and distribution of Luoyang (b), isohyet of Luoyang in 2019 (c) and groundwater depth of Luoyang in 2019 (d).
Water 15 01491 g002
Figure 3. The effects of synergistic development.
Figure 3. The effects of synergistic development.
Water 15 01491 g003
Figure 4. Relationship of each part of the DPSIR model in the SD-WALRs.
Figure 4. Relationship of each part of the DPSIR model in the SD-WALRs.
Water 15 01491 g004
Figure 5. Trends of the development level of the water resources.
Figure 5. Trends of the development level of the water resources.
Water 15 01491 g005
Figure 6. Trends of the development level of the land resources.
Figure 6. Trends of the development level of the land resources.
Water 15 01491 g006
Figure 7. Trends of the C, T and SD and the evaluation classifications of the SD.
Figure 7. Trends of the C, T and SD and the evaluation classifications of the SD.
Water 15 01491 g007
Figure 8. Normalized data of the WALR indexes: (a1,a2,a3) normalized data of the water resources subsystem indexes in 2010, 2014 and 2019, respectively; (b1,b2,b3) normalized data of the land resources subsystem indexes in the corresponding years.
Figure 8. Normalized data of the WALR indexes: (a1,a2,a3) normalized data of the water resources subsystem indexes in 2010, 2014 and 2019, respectively; (b1,b2,b3) normalized data of the land resources subsystem indexes in the corresponding years.
Water 15 01491 g008
Figure 9. Part of the indexes with remarkable changes from 2010 to 2019: (a,b) main driving force indexes of water and land respectively, (c) trends of two representative resources utilization states, (d) trends of the matching degree and GDP per unit resources, and (e) measures taken to avoid degradation by humans.
Figure 9. Part of the indexes with remarkable changes from 2010 to 2019: (a,b) main driving force indexes of water and land respectively, (c) trends of two representative resources utilization states, (d) trends of the matching degree and GDP per unit resources, and (e) measures taken to avoid degradation by humans.
Water 15 01491 g009
Figure 10. Evaluation results of this research and Lv’s method.
Figure 10. Evaluation results of this research and Lv’s method.
Water 15 01491 g010
Table 1. Evaluation index system of the development level of water resources subsystems.
Table 1. Evaluation index system of the development level of water resources subsystems.
ItemIndexUnitNotes or Calculation MethodItemIndexUnitNotes or Calculation Method
Driving forceGDP (+) (D1)CNY 108The drive for socio-economic evolutionStateGini coefficient (−) (S2)/Spatial matching degree of WALRs, calculation method from Li [9]
Resident population (+) (D2)104 p
Urbanization rate (+) (D3)%Urban/total populationIrrigation water per unit area (+) (S3)Current water-use efficiency
Proportion of water worse than grade V (−) (D4)%Length of river worse than grade V/total river lengthWater consumption per unit of industrial added value (−) (S4)m3/CNY 104
PressureLand-use structureCultivated area (−) (P1)km2Pressure of land-use structure on water resources subsystemUrban domestic comprehensive water quota
(−) (S5)
m3/d
Plantation area (−) (P2)km2Rural domestic comprehensive water quota
(−) (S6)
m3/d
Industrial area (−) (P3)km2ImpactGDP per unit water (+) (I1)CNY/m³GDP/water consumption
Residential building area (−) (P4)km2Water resources exploitation ratio (−) (I2)%Total water consumption/total regional water resources
Greenland aera (−) (P5)km2Engel coefficient (−) (I3)%Food expenditure/total personal expenditure
Water resources per capita (+) (P6)m3/pPressure of regional resources on water resources subsystemResponseInvestment of water conservancy project (+) (R1)CNY 104Total investment in water conservancy project construction and management
Total water consumption (−) (P7)108 m3Ratio of sewage treatment (+) (R2)%Measures taken to avoid system function degradation
Volume of sewage (−) (P8)108 m3Area of soil erosion controlled (+) (RW)km2
StateMatching coefficient of WALRs (+) (S1)/Quantity matching degree of WALRs, calculation method from Li [9]
Table 2. Evaluation index system of synergistic development of land resources.
Table 2. Evaluation index system of synergistic development of land resources.
ItemIndexUnitNotes or Calculation MethodItemIndexUnitNotes or Calculation Method
Driving forceGDP (+) (D1)CNY 108The drive for socio-economic evolutionStateMatching coefficient of WALRs (+) (S1)/The same as Table 1
Resident population (+) (D2)104 pGini coefficient (−) (S2)/
Urbanization rate (+) (D3)%Proportion of effective irrigation area (+) (S3)%Present situation of land use
Ratio of vegetation coverage (+) (D4)%The drive for eco-environment evolutionGDP per unit industrial area (+) (S4)CNY/m2
PressureWater-use struct-ureProportion of agricultural water use (−) (P1)%Pressure of water -use structure on land resources subsystemPer capita housing area (+) (S5)m3
Proportion of industrial water use (−) (P2)%ImpactGDP per unit area (+) (I1)CNY/m3GDP/land development area
Proportion of domestic water for residents (−) (P3)%Land diversity index (+) (I2)/ H = i = 1 n p i ln p i
H: land diversity index;
pi: proportion of land type i;
n: number of land-use types
Proportion of ecological water use (−) (P4)%
Land resources per capita (+) (P5)m2/pPressure of regional resources on land resources subsystemEngel coefficient (−) (I3)%The same as Table 1
ResponseUnit land assets investment (+) (R1)CNY/m2Measures taken to avoid system function degradation
Proportion of regional development (−) (P6)%Construction land area/total areaArea of soil erosion controlled (+) (R2)km2
Ecological index of cultivated land (−) (P7)kg/m2(Total amount of pesticide, chemical fertilizer and agricultural film)/cultivated areaHarmless treatment ratio of solid waste (+) (R3)%
Note: The raw data of all indicators are from Henan Statistical Bureau (HSB 2010–2019), Luoyang Statistical Bureau (LSB 2010–2019), Luoyang Water Conservancy Bureau (LWCB 2010–2019) and “Luoyang City land use overall planning”.
Table 3. Classification of the synergistic development level.
Table 3. Classification of the synergistic development level.
ClassesRange of IndexTypes of SD-WALRs
Unbalanced recession0 < SD ≤ 0.1Extremely unbalanced recession
0.1 < SD ≤ 0.2Severely unbalanced recession
0.2 < SD ≤ 0.3Intermediate unbalanced recession
0.3 < SD ≤ 0.4Mildly unbalanced recession
Transition stage0.4 < SD ≤ 0.5On the verge of recession
0.5 < SD ≤ 0.6Reluctantly synergetic development
Synergetic development0.6 < SD ≤ 0.7Primary synergetic development
0.7 < SD ≤ 0.8Intermediate synergetic development
0.8 < SD ≤ 0.9Good synergetic development
0.9 < SD ≤ 1High-quality synergetic development
Table 4. Weights of indicators found using three methods.
Table 4. Weights of indicators found using three methods.
Index Weight of the Water Resources SubsystemIndex Weights of the Land Resources Subsystem
ItemsωAHMωEWλiItemsωAHMωEWλi
D10.02780.03580.0236D10.02780.03520.0210
D20.03690.05060.0444D20.03690.04970.0393
D30.01530.03860.0140D30.01530.03790.0124
D40.02000.02550.0121D40.02000.04900.0210
P10.03230.04500.0345P10.02120.03480.0158
P20.01670.04490.0178P20.04590.05910.0581
P30.04480.03960.0422P30.01810.04130.0160
P40.01920.03010.0137P40.01170.05870.0147
P50.01880.03610.0161P50.08810.04480.0845
P60.05420.04180.0538P60.07700.06240.1029
P70.06580.04140.0647P70.03800.03640.0296
P80.04830.05670.0651S10.04610.04690.0463
S10.04420.04770.0501S20.04610.03480.0344
S20.04420.03540.0372S30.09050.06800.1318
S30.07310.04820.0837S40.08380.03090.0555
S40.08690.02810.0580S50.03350.03650.0262
S50.02580.03300.0202I10.07550.03420.0553
S60.02580.05890.0361I20.03750.03800.0305
I10.07550.03640.0653I30.03750.03850.0309
I20.03750.03170.0282R10.04950.04600.0488
I30.03750.03920.0349R20.05000.04970.0532
R10.04950.05740.0675R30.05000.06690.0717
R20.05000.04760.0565
R30.05000.05050.0600
Table 5. Evaluation results of the development level of the WALR subsystems. Unit: 10−2.
Table 5. Evaluation results of the development level of the WALR subsystems. Unit: 10−2.
ItemYear2010201120122013201420152016201720182019
Water resources subsystemD0.001.592.313.153.905.146.336.718.519.42
P19.7017.5011.427.2813.0117.3616.4813.5311.819.41
S10.4615.5115.8117.717.5312.948.259.6612.4315.61
I4.075.034.081.834.046.307.769.4210.5310.79
R2.166.506.1410.951.114.8112.3211.7711.5812.00
U36.3946.1239.7740.9239.6046.5551.1351.0854.8557.22
Land resources subsystemD0.020.591.242.082.994.105.216.518.139.37
P23.6819.6519.7419.9222.0516.0110.0811.929.688.84
S4.9110.979.344.468.4119.3317.6420.9125.3626.03
I1.112.513.453.174.626.547.809.5010.6711.43
R6.560.401.416.331.576.218.1910.7111.2812.04
U36.2834.1435.1835.9539.6452.1948.9159.5465.1267.7
Table 6. Evaluation results of the SD of water and land resources.
Table 6. Evaluation results of the SD of water and land resources.
YearCTSDLevel of the SD-WALRs
20100.9990.3630.603Primary synergetic development
20110.9890.4010.630Primary synergetic development
20120.9980.3750.612Primary synergetic development
20130.9980.3840.619Primary synergetic development
20140.9990.3960.629Primary synergetic development
20150.9980.4940.702Intermediate synergetic development
20160.9990.5000.707Intermediate synergetic development
20170.9970.5530.743Intermediate synergetic development
20180.9960.6000.773Intermediate synergetic development
20190.9960.6250.789Intermediate synergetic development
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

Lv, C.; Xu, W.; Ling, M.; Wang, S.; Hu, Y. Evaluation of Synergetic Development of Water and Land Resources Based on a Coupling Coordination Degree Model. Water 2023, 15, 1491. https://doi.org/10.3390/w15081491

AMA Style

Lv C, Xu W, Ling M, Wang S, Hu Y. Evaluation of Synergetic Development of Water and Land Resources Based on a Coupling Coordination Degree Model. Water. 2023; 15(8):1491. https://doi.org/10.3390/w15081491

Chicago/Turabian Style

Lv, Cuimei, Wenchao Xu, Minhua Ling, Sensen Wang, and Yuguang Hu. 2023. "Evaluation of Synergetic Development of Water and Land Resources Based on a Coupling Coordination Degree Model" Water 15, no. 8: 1491. https://doi.org/10.3390/w15081491

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