This study combined emergy theory with water ecological footprint to analyze and evaluate the sustainable development status of water resources in Yunnan Province from 2012 to 2021. The water pollutant account and the water product account were supplemented to the water ecological footprint. The results showed that the water resources in Yunnan Province were generally in an ecological deficit during the study period, but their absolute values were very small. The ecological pressure index of water resources was high, while the sustainable utilization index of water resources was low, indicating that the water resources were in a state of unsustainable utilization. However, the emergy ecological footprint of water resources and the emergy ecological footprint of water resources of 10,000 Yuan GDP showed a decreasing trend, and the eco-economic coordination index of water resources was relatively high, indicating that the water resources will show a sustainable development trend in the future. The sustainable development of water resources was influenced by both natural ecological factors and social factors. Therefore, suggestions such as improving the utilization rate of water resources development, developing a green low-carbon economy, preventing and controlling water pollution, and protecting the water environment are conducive to promoting the sustainable development of water resources.

  • Provided a new research method for the sustainable development of water resources in karst cities.

  • Supplemented the water pollutant account and the aquatic product account.

  • Calculated energy density and the average value of regional water resources.

  • Considering natural ecological and social factors, and proposing optimal paths.

The sustainable use of water resources is an important part of the construction of ecological civilization and is also the main content of the study of resources and environment (Oliveira 2017; Gleeson et al. 2020; Tian et al. 2021; Meng & Wang 2022). Sustainable development and utilization of water resources is an important prerequisite for regional economic development and is closely related to regional sustainable development, and has important strategic significance and special functions in the region (Chen et al. 2016; Larsen et al. 2016; Ali et al. 2019; Song et al. 2022). Behind the rapid economic development and population growth, the demand for water resources from industry and society has gradually increased (Flörke et al. 2018; Wang et al. 2021a; Darko et al. 2022; Mann & Gupta 2022). The uncontrolled abuse of water resources has caused ecological and environmental problems such as water depletion, water environment pollution, soil erosion, and water shortage (Liu et al. 2021; Michelle et al. 2021; Faye 2022; Hatami et al. 2022; Hojjati-Najafabadi et al. 2022; Karaman et al. 2022; Nguyen et al. 2022). The ‘2023 United Nations World Water Development Report’ jointly published by UNESCO and UN-Water pointed out that 2 billion people (approximately 25% of the world's population) currently lack access to safe drinking water, and 2–3 billion people experience water scarcity for at least one month each year. If international cooperation is not strengthened, the problem of water shortage will become more and more serious in the coming decades, with the global urban population facing water shortage expected to grow from 930 million in 2016 to 1.7 billion to 2.4 billion in 2050 (Liu 2023). The United Nations Water Conference 2023 adopted the ‘Water Action Agenda’. In the face of the increasingly severe water resources situation in the world, the participants discussed solutions and called for strengthened international cooperation to jointly promote the sustainable development of water resources (Li & Zheng 2023). China's per capita water resources are at a backward level in the world, making it one of the world's water-poor countries, and more than 60% of large- and medium-sized cities in China suffer from water shortage problems (Wang & Li 2021). With economic development and population growth, the demand for water resources and aquatic products is increasing and the discharge of wastewater is increasing, so the sustainable development of water resources is affected. How to improve the current situation of water resources, maintain the sustainable utilization of water resources, and ensure the sustainable development of water resources have become an important link to promote the coordinated development of population, resources, economy, and ecological environment.

The sustainable development of water resources is a current research hotspot (Ishaque et al. 2023; Islam et al. 2023; Mikulčić et al. 2023; Rahmani et al. 2023; Zhang et al. 2023). Many scholars have carried out research work on the sustainable development of water resources. Farrokhzadeh et al. (2020) used a multi-objective optimization model linked with the Water Evaluation and Planning (WEAP) software to optimize water allocation decisions over multiple years and make optimal planning in water resources necessary for sustainable socio-economic development in the Sistan region and its Hamoun wetland, located in southeast Iran. Ezzeldin et al. (2022) put forward a rainwater harvesting technique that combines geographic information systems (GIS), remote sensing (RS), multi-criteria analysis, and hydrological modeling in a case study in Wadi Watir in the Sinai Peninsula, Egypt. This research addresses the challenges of water scarcity along with socio-economic and environmental pressures while achieving sustainable development goals. Abedzadeh et al. (2020) used Fuzzy Fault Tree Analysis to assess risk-based assessment of water resources development plans (WRDPs) in Homozgan province in the South of Iran, under the sustainable development framework. The proposed approach not only addresses the risk of WRDPs in compliance with sustainable development objectives but also facilitates decision-making for risk management by prioritizing the factors in the failure of plans. Ifediegwu (2022) used an integrated method to assess the groundwater potential zones in the Lafia district utilizing RS, GIS, and the analytic hierarchy process (AHP). The results can aid sustainable groundwater resource management in this area. Shams & Muhammad (2022) proposed a novel framework for a critical assessment of the implementation of integrated water resources management and water–energy–food nexus in Afghanistan toward sustainable water resources management. The above-mentioned studies have carried out the sustainable development of water resources in different regions from different perspectives, and certain research results have been achieved. However, water resources not only have ecological attributes but also have socio-economic attributes. The evaluation of sustainable development of regional water resources is a composite system that organically combines socio-economic, resource, and environmental factors. The above methods cannot better reflect the complexity and comprehensiveness of the regional water resources sustainable development system.

The ecological footprint method proposed by Rees (1992) and Wackernagel et al. (1999) can make up for the shortcomings of the above methods. However, the information covered by the data accounts collected by the ecological footprint model of water resources is not comprehensive, and there is a lack of water pollutant accounts and aquatic product accounts. Most of the selected parameters are global equilibrium factors and yield factors, without considering regional and time differences, which can cause errors between the calculation results and the actual situation, and cannot accurately reflect the utilization of water resources in the region. In view of the above reasons, firstly, the water ecological footprint account is supplemented with a water pollution account and a water product account in addition to a water consumption account (including agricultural water account, industrial water account, domestic water account, and ecological environment water account), which further improves the water ecological footprint account and more accurately reflects the impact of regional water pollution on the ecosystem and the pressure of economic development on the ecosystem. Secondly, this paper improves it with the help of emergy theory (Odum 1996), converting different kinds of energy in each account into unified solar emergy through solar emergy transformity. After unification, various forms of energy can be quantitatively analyzed and studied, which truly reflects the relationship between the water environment and economic and social development. Then, in the emergy ecological footprint model of water resources, the emergy ecological footprint of water resources is calculated by using the regional emergy density of water resources in the current year, and the emergy ecological carrying capacity of water resources is calculated by using the regional average emergy density of water resources. The calculation results are more in line with the actual situation of regional water resources ecological sustainable development, and the evaluation results are more reasonable. The water ecological footprint assessment method based on the emergy theory has the characteristics of more scientific, perfect, operational, and practical value, and it comprehensively considers the influence of social and environmental factors. It is one of the most important ways to study the status of regional water resource utilization, judge the carrying capacity of water resources, and evaluate the sustainable development capacity of regional water resources.

This paper combines the emergy theory and water ecological footprint to take advantage of the strengths of these two theoretical models. The emergy ecological footprint model of water resources is composed of the emergy ecological footprint of water consumption, the emergy ecological footprint of water pollution, and the emergy ecological footprint of aquatic products. The emergy ecological carrying capacity of water resources is composed of surface water chemical energy, groundwater chemical energy, and rainwater potential energy. Then, according to the regional water resources' emergy ecological carrying capacity, emergy ecological footprint, and other related data, five water resources sustainable development evaluation indexes were calculated, namely emergy ecological deficit/surplus of per capita water resources, ecological pressure index of water resources, sustainable utilization index of water resources, eco-economic coordination index of water resources, and emergy ecological footprint of water resource of 10,000 Yuan GDP. The current situation of sustainable utilization of water resources in Yunnan Province from 2012 to 2021 was analyzed and evaluated, and the internal driving factors were analyzed from natural ecological factors and social factors, and the optimization paths of water resources sustainable development were proposed to improve the current situation. The research results can provide new ideas for the evaluation and analysis of sustainable development of water resources in the karst region, which is valuable for alleviating the contradiction between water supply and demand of water resources and promoting the sustainable development of regional society and economy and the virtuous cycle of water ecological environment.

Emergy ecological footprint theory of water resources

Canadian ecological economist Rees (1992) pioneered the concept of ecological footprint and proposed the calculation method of the ecological footprint model in 1996 in collaboration with his doctoral student Wackernagel et al. (1999) to determine whether the area was sustainable by comparing regional ecological footprint and ecological carrying capacity. Odum (1996) proposed a new theory and method of emergy analysis, which can convert different types of energy into a unified unit of solar emergy through the solar emergy transformity, and then quantitatively analyze and study it to truly reflect the relationship between the ecological environment and economic and social development. Later, some scholars introduced the emergy theory into the ecological footprint model, and through continuous practice and application, it developed into a more systematic emergy ecological footprint model to determine the sustainable development of the region (Zadgaonkar & Mandavgane 2020; Liu et al. 2022). According to the above theory, this paper tried to combine the solar emergy in the emergy theory with the ecological carrying capacity and ecological footprint in the ecological footprint model to maximize the advantages of these two theoretical models. Ecosystems provide renewable and non-renewable resources. Non-renewable resources are gradually consumed, while renewable resources can be used continuously. Therefore, this paper selects water resources as the renewable resources of Yunnan Province in the study area. The emergy ecological carrying capacity indicators of water resources include surface water chemical energy, groundwater chemical energy, and rainwater potential energy. The emergy ecological footprint indicators of water resources include water consumption account, water pollutant account, and aquatic product account. Meanwhile, with the solar emergy transformity as a parameter, different types of energy in each account can be converted into the same form of solar emergy. Combined with five evaluation indexes of sustainable development of water resources, the current situation of sustainable utilization of water resources in the region can be analyzed and evaluated, which will more accurately simulate the dynamic change trend of sustainable development of regional water resources.

Emergy ecological carrying capacity model of water resources

Emergy theory

Emergy is the amount of energy in another form contained in the flowing or stored energy of existing materials on earth (Wang et al. 2021b). The material energy in the earth exists in different forms, and the energy cannot be directly added or subtracted from each other, but the energy values can be converted into solar emergy, and the solar emergy transformity is a way to convert them to each other (Lv 2009). The emergy is calculated as follows:
(1)
where E is the emergy; B is the available energy of the resources or material; and T is the solar emergy transformity. The main solar emergy transformity in this article is shown in Table 1.
Table 1

Main solar emergy transformity

AccountSolar emergy transformityReferences
Agricultural water 1.60 × 1012 sej/m3 Li et al. (2022a, 2022b
Industrial water 2.32 × 1012 sej/m3 Li et al. (2022a, 2022b
Domestic water 8.80 × 1011 sej/m3 Li et al. (2022a, 2022b
Eco-environment water 1.263 × 1012 sej/m3 Huang et al. (2019)  
Aquatic product 2.00 × 106 sej/t Lian (2014)  
Water pollution 4.94 × 1012 sej/t Yang (2016)  
Surface water chemical energy 1.54 × 104 sej/J Wang et al. (2023)  
Groundwater chemical energy 1.54 × 104 sej/J Wang et al. (2023)  
Rainwater potential energy 8.89 × 103 sej/J Wang et al. (2023)  
AccountSolar emergy transformityReferences
Agricultural water 1.60 × 1012 sej/m3 Li et al. (2022a, 2022b
Industrial water 2.32 × 1012 sej/m3 Li et al. (2022a, 2022b
Domestic water 8.80 × 1011 sej/m3 Li et al. (2022a, 2022b
Eco-environment water 1.263 × 1012 sej/m3 Huang et al. (2019)  
Aquatic product 2.00 × 106 sej/t Lian (2014)  
Water pollution 4.94 × 1012 sej/t Yang (2016)  
Surface water chemical energy 1.54 × 104 sej/J Wang et al. (2023)  
Groundwater chemical energy 1.54 × 104 sej/J Wang et al. (2023)  
Rainwater potential energy 8.89 × 103 sej/J Wang et al. (2023)  
Emergy density is the size of solar energy consumed by a country or region in a year, which mainly reflects the ecological and economic development level of the region. The better the ecological environment and the more active the economic development, the higher the emergy density is. The calculation formula is as follows:
(2)
(3)
(4)
where P is the regional emergy density of water resources (sej/hm2); Eg, Eu, and Er are the chemical energy of surface water, chemical energy of groundwater, and geopotential energy of rain (sej); S is the area of the region (3.94 × 107 hm2); Mg is surface water (m3); Mu is groundwater (m3); G is the Gibbs free energy, surface water Gibbs free energy is 4.94 J/g, groundwater Gibbs free energy is 4.90 J/g; ρ is the density of water (1,000 kg/m3); Rain is the average annual precipitation (mm); h is the average altitude (2,000 m); g is the acceleration of gravity (9.8 m/s2); Tg is the solar emergy transformity of surface water chemical energy; Tu is the solar emergy transformity of groundwater chemical energy; and Tr is the solar emergy transformity of rainwater potential energy.

Emergy ecological carrying capacity of water resources

The ecological carrying capacity of water resources is the maximum economic and social development pressure that can be carried by the water ecosystem under certain conditions in a certain period in the region (Huang et al. 2008; Wu et al. 2021; Shukui et al. 2022; Wang et al. 2022). According to the emergy theory, the water resources emergy in Yunnan Province is mainly composed of rainwater potential energy, surface water chemical energy, and groundwater chemical energy. Also considering the water resource utilization rate, at least 60% of water should be deducted to maintain biodiversity and ecological environment development (Liu & Yang 2021). The calculation formula is given as follows:
(5)
where Pw is the regional average emergy density of water resources (sej/hm2), which is the average value of regional emergy density of water resources in Yunnan Province in recent 10 years, and is 2.56 × 1015 sej/hm2. WEC is the emergy ecological carrying capacity of water resources (hm2); and Wec is the emergy ecological carrying capacity of per capita water resources (hm2/cap). The rest of the symbols have the same meaning as above.

Emergy ecological footprint model of water resources

Water ecological footprint is the conversion of the water consumed by a country or region to meet the sustainable development of regional society and economy into a biologically productive land area with special properties (Zhang & Zhang 2013). The emergy ecological footprint model of water resources is constructed by integrating emergy theory and the ecological footprint model. According to the actual situation in Yunnan Province, the primary account of the emergy ecological footprint model of water resources is divided into three secondary accounts: emergy ecological footprint of water consumption, emergy ecological footprint of water pollution, and emergy ecological footprint of aquatic products. The calculation formula is as follows:
(6)
where WEF, WEFC, WEFP, and WEFA are the emergy ecological footprint of water resources, emergy ecological footprint of water consumption, emergy ecological footprint of water pollution, and emergy ecological footprint of aquatic product (hm2), respectively; N is the total population in the study area (10,000 people); and Wef is the emergy ecological footprint of per capita water resources (hm2/cap).

Emergy ecological footprint of water consumption

The emergy ecological footprint of water consumption is the footprint of freshwater resources consumed by urban and rural residents in production and life (mainly including agricultural water consumption, industrial water consumption, domestic water consumption, and ecological environment water consumption) (Li et al. 2022a, 2022b). The agricultural water account includes farmland irrigation water, forest and fruit land irrigation water, grassland irrigation water, fishing pond replenishment water, and large-scale livestock and poultry breeding water. The industrial water account refers to the water used by industrial and mining enterprises for manufacturing, processing, cooling, air conditioning, purification, and washing in the production process, excluding the reuse of water within the enterprises. The domestic water account refers to urban domestic water and rural domestic water. Urban domestic water use includes residential water use and public water use, and rural domestic water use refers to rural residents' domestic water use. The ecological environment water account includes water for the urban environment supplied by artificial measures and water consumption for replenishment of some rivers, lakes, and wetlands, excluding the amount of water naturally satisfied by precipitation and runoff. The calculation formula is as follows:
(7)
where WRi represents the emergy ecological footprint of water consumption of the i-th category of water accounts (hm2), i = 1, 2, 3, 4 represent agricultural, industrial, domestic, and eco-environment water consumption accounts, respectively. Wri is the emergy ecological footprint of per capita water consumption of the i-th category of water accounts (hm2/cap); Ci is the total water consumption of the i-th category (m3); and Ti is the solar emergy transformity of water consumption of the i-th category (sej/m3). The rest of the symbols have the same meaning as above.

Emergy ecological footprint of water pollution

The emergy ecological footprint of water pollution is mainly the footprint of wastewater discharged by urban and rural residents' production and living (Li et al. 2022b). The total amount of wastewater discharge consists of domestic wastewater discharge and industrial wastewater discharge. Because the solar emergy transformity of domestic sewage and industrial wastewater are the same, the total amount of wastewater discharge is selected for calculation. The calculation formula is as follows:
(8)
where Q is the total amount of wastewater discharge (t); T1 is the solar emergy transformity of wastewater (sej/t); and Wefp is the emergy ecological footprint of per capita water pollution (hm2/cap). The rest of the symbols have the same meaning as above.

Emergy ecological footprint of aquatic products

The emergy ecological footprint of aquatic products is the footprint of producing various types of aquatic products (Li et al. 2022b). The calculation formula is as follows:
(9)
where A is the yield of aquatic products (t); λ is the energy conversion coefficient of the aquatic product (5.59 × 109 j/t) (Lian 2014); T2 is the solar emergy transformity of the aquatic product (sej/J); and Wefa is the emergy ecological footprint of per capita aquatic product (hm2/cap). The rest of the symbols have the same meaning as above.

Comprehensive evaluation indicators for sustainable development of water resources

Since there are differences among regions in terms of water resources, climate conditions, land area, population size, and economic level, it is not comprehensive enough to evaluate the sustainable development status of regions by only relying on two basic indicators (emergy ecological carrying capacity and emergy ecological footprint). In order to more objectively evaluate the regional sustainable development, on the basis of the original, five sustainable evaluation indicators (emergy ecological deficit/surplus of per capita water resources, ecological pressure index of water resources, sustainable utilization index of water resources, eco-economic coordination index of water resources, and emergy ecological footprint of water resource of 10,000 Yuan GDP) are comprehensively considered in order to make a more accurate evaluation of the ecological status of the region, so as to put forward more reasonable suggestions.

Emergy ecological deficit/surplus of per capita water resources

The emergy ecological deficit/surplus of per capita water resources is the difference between the emergy ecological carrying capacity of per capita water resources and the emergy ecological footprint of per capita water resources, which is used to judge whether the regional water resources utilization is in ecological surplus or shortage, and can clearly indicate the degree of sustainable utilization and the supply and demand situation of regional water resources. The calculation formula is:
(10)

When Wec = Wef, the supply and demand of water resources are in balance, indicating that the emergy ecological carrying capacity of regional water resources can meet the demand of urban and rural residents for production and living water. When Wef < Wec, there is an ecological surplus of water resources, indicating that the regional water resources supply can meet the demand, which is conducive to the sustainable development of the region. When Wef > Wec, there is an ecological deficit of water resources, indicating that the regional socio-economic development exceeds its ecological load on water resource utilization, which is not conducive to the sustainable development of the region.

Ecological pressure index of water resources

The ecological pressure index of water resources (WEPI) can be used to judge the degree of regional water resources development and utilization, and its value is the ratio of regional water resources' ecological footprint to water resources ecological carrying capacity. The calculation formula is:
(11)

When WEPI is close to 1, it indicates that the regional water resources supply and consumption are basically equal and are in the critical state of sustainable development of water resources. When 0 < WEPI < 1, it indicates that the regional water resources supply can support its consumption status and level. When WEPI > 1, it indicates that the regional water resources are in an unsustainable state, and the larger the value, the greater the regional water resources' ecological pressure. At this time, measures should be taken to improve the water's ecological carrying capacity.

Sustainable utilization index of water resources

The sustainable utilization index of water resources (WSUI) is expressed as the ratio of the difference between WEC and WEF to WEC. WSUI refers to the degree to which the total water resources that can be sustainably supplied in the region meet the demand for water resources from human socio-economic activities within a specific period of time. The calculation formula is as follows (Li et al. 2021):
(12)
  • (1)

    When 0 < WSUI ≤ 1, it means that the emergy ecological carrying capacity of water resources is surplus and can support the future emergy ecological footprint of water resources, which does not cause damage to the ability of sustainable use of water resources in the future, and the water resources are in a sustainable utilization state; when the value of WSUI is large, the water resources are in a strong sustainable utilization state; when the value of WSUI is small, the water resources are in a weakly sustainable utilization state.

  • (2)

    When WSUI < 0, it means that the ecological carrying capacity of water resources is not enough to support the ecological footprint of water resources, which has already caused damage to the ability of sustainable use of water resources in the future, and the water resources are in the state of unsustainable utilization.

  • (3)

    When WSUI = 0, the water resources are at the critical point of sustainable and unsustainable utilization.

Eco-economic coordination index of water resources

The eco-economic coordination index of water resources (WECI) can reflect the coordination between the water resources ecosystem and social economy, whether the regional socio-economic development is within the ecological carrying range, and the relationship between the emergy ecological footprint and the regional resource endowment. The calculation formula is:
(13)

The value range of WECI is (1, 1.414). The closer the WECI is to 1, the worse the ecological and economic coordination of water resources; the closer the WECI is to 1.414, the better the ecological and economic coordination of water resources. When WECI = 1.414, it means that the water resources are in the best ecological and economic coordination state.

Emergy ecological footprint of water resources of 10,000 Yuan GDP

The utilization efficiency of water resources can be expressed by the emergy ecological footprint of water resources of 10,000 Yuan GDP, which is the ratio of the emergy ecological footprint of regional water resources to GDP. The smaller the result, the higher the utilization efficiency of water resources is. The calculation formula is:
(14)
where WGDP is the emergy ecological footprint of water resources of 10,000 Yuan GDP (hm2/10,000); and GDP is the gross domestic product in the study area (10,000 Yuan).

Study area overview

Yunnan Province is located in the southwest of China, with geographical coordinates of 21°8′–29°15′ north latitude and 97°31′–106°11′ east longitude. It occupies a very advantageous geographical position. It is bordered by Guizhou Province and Guangxi Province in the east, Sichuan Province in the north, Tibet in the northwest, Myanmar in the west, Laos in the south, and Vietnam in the southeast. The total area of the province is 394,100 km2. The geographical location is shown in Figure 1. From the perspective of topography, Yunnan Province is mainly mountainous terrain, and its mountainous area accounts for more than 80% of the province's area. The terrain of Yunnan Province is inclined from northwest to southeast, showing a step-by-step decline, and the altitude difference is extremely different. Most of the areas are in the range of 1,000–3,500 m, and the average altitude is about 2,000 m. The study area is a typical mountainous plateau topography, and the unique natural geographic and regional conditions have created a variety of landform types, with rolling low mountains and hills, developing a variety of distinctive karst landforms. From the perspective of climatic conditions, Yunnan's climate is basically of the subtropical plateau monsoon type, with significant three-dimensional climatic characteristics, numerous types, small annual temperature differences, large daily temperature differences, distinct dry and wet seasons, and unusual vertical variations in temperature with topographic height. In addition, the precipitation in the province has the characteristics of uneven distribution in season and region. From the perspective of natural resources, Yunnan Province has the advantages of abundant mineral resources, flora and fauna resources, tourism resources, water resources, and energy resources. Among the energy resources, water resources and coal resources are relatively abundant, and the province has many rivers and lakes, which are characterized by abundant water resources and sufficient precipitation.
Figure 1

Geographical location of Yunnan Province.

Figure 1

Geographical location of Yunnan Province.

Close modal

The distribution area of karst in the world is nearly 2.2 × 107 km2, and the karst area in China is 3.4 × 106 km2, accounting for about one-third of the land area, mostly distributed in the southern regions with a warm and humid climate, involving many provinces such as Sichuan, Chongqing, Hunan, Guangdong, Guizhou, Guangxi, and Yunnan. The southwestern karst area with Guizhou, Guangxi, and Yunnan as the main body is one of the largest contiguous exposed carbonate rock distribution areas in the world. The area of Yunnan Province is 3.941 × 105 km2, of which 1.1 × 105 km2 is covered by karst landforms, accounting for 28.15% of the total area. Counties with karst landform area, which accounts for over 30% of the total area, are called karst-developed counties, while counties with karst area accounting for over 70% of the total area are called severe karst counties. There are 116 karst-developed counties and 65 severe karst counties in Yunnan Province (Chang 2023).

Data sources

This paper collates and collects relevant data on Yunnan Province from 2012 to 2021. The data on annual rainfall, agricultural water consumption, industrial water consumption, domestic water consumption, ecological and environmental water consumption, water consumption of 10,000 Yuan GDP, surface water resources, and underground water resources are obtained from ‘Yunnan Water Resources Bulletin’ (2012–2021), and the data of GDP, sewage discharge, aquatic products production, sewage discharge, aquatic products production, and year-end resident population are obtained from ‘Yunnan Statistical Yearbook’ (2013–2022). See Supplementary material for the specific collated categories.

Analysis of sustainable development of water resources

Analysis of emergy ecological carrying capacity of water resources

According to the average annual precipitation data from 2012 to 2021 in Supplementary Table S1 and the emergy ecological carrying capacity of water resources, rainwater, surface water, and groundwater in Supplementary Table S2, the emergy ecological carrying capacity of water resources and the average annual precipitation change curve are shown in Figure 2.
Figure 2

Changes in emergy ecological carrying capacity of water resources and average annual precipitation in Yunnan Province from 2012 to 2021.

Figure 2

Changes in emergy ecological carrying capacity of water resources and average annual precipitation in Yunnan Province from 2012 to 2021.

Close modal

As a renewable resource, water resources mainly come from rainwater, surface water, and groundwater. In 2012–2021, the fluctuation range of water resources emergy in Yunnan Province was 8.51 × 1022 ∼ 1.15 × 1023 sej, and the emergy ecological carrying capacity of water resources showed a fluctuating upward trend, with an overall increase rate of about 2.08%. Among them, the emergy ecological carrying capacity of water resources was the highest in 2017, which was 1.80 × 107 hm2, corresponding to an average annual precipitation of 1,351.5 mm. The lowest value was 1.33 × 107 hm2 in 2019, corresponding to an average annual precipitation of 1,008 mm. The emergy ecological carrying capacity of rainwater was the main contributor to the emergy ecological carrying capacity of water resources, accounting for 81.36% on average. The emergy ecological carrying capacity of surface water was 1.83 × 106 ∼ 2.63 × 106 hm2, accounting for 13.93% on average. The emergy ecological carrying capacity of groundwater was the smallest, and the annual variation range was not large. The average annual value was 7.44 × 105 hm2, accounting for 4.72% on average. The annual variation of the emergy ecological carrying capacity of water resources in Yunnan Province was closely related to the annual variation of the average annual precipitation. The emergy ecological carrying capacity of water resources increased with the increase of rainfall. It showed that the average annual precipitation in Yunnan Province played a decisive role in the emergy ecological carrying capacity of water resources. On the one hand, this reflects the weak self-regulation ability of water resources carrying capacity in Yunnan Province, which relies heavily on atmospheric precipitation. On the other hand, the intensity of human activities has increased with the continuous development of society, such as the rapid development of tourism, accelerated urbanization, and rapid economic development, which have caused excessive use of water resources and water pollution, resulting in the reduction of human available water resources and further affecting the emergy ecological carrying capacity of water resources.

Analysis of emergy ecological footprint of water resources

According to the data on the emergy ecological footprint of agricultural water consumption, industrial water consumption, domestic water consumption, ecological environment water consumption, emergy ecological footprint of water pollution, and emergy ecological footprint of aquatic products from 2011 to 2020 in Supplementary Table S3, the ecological footprint changes of water resources accounts in Yunnan Province from 2012 to 2021 are shown in Figure 3.
Figure 3

Changes in emergy ecological footprint of various water resources accounts in Yunnan Province from 2012 to 2021.

Figure 3

Changes in emergy ecological footprint of various water resources accounts in Yunnan Province from 2012 to 2021.

Close modal

During the study period, the emergy ecological footprint of water resources decreased from 1.72 × 107 hm2 to 1.68 × 107 hm2, with an overall decrease rate of about 2.31%. Among the ecological footprint accounts, the water consumption account had the greatest impact on the ecological footprint of water resources, with an average contribution rate of 59.59%. The aquatic product account and the water pollution account had a relatively small impact on the ecological footprint of water resources, with an average contribution rate of 20.29 and 20.12%, respectively. The water consumption account was divided into four types of water accounts: agricultural, industrial, domestic, and eco-environment. The agricultural water consumption account (WR1) accounted for 41.46%, the industrial water consumption account (WR2) accounted for 12.46%, the domestic water consumption account (WR3) accounted for 4.73%, and the ecological environment water consumption account (WR4) accounted for 0.94%. The emergy ecological footprint of agricultural water consumption accounted for the largest proportion, indicating that the water resources in Yunnan Province were mainly invested in agriculture, and the agricultural water use footprint showed a fluctuating upward trend in general, with an increased rate of 6.20%. The overall emergy ecological footprint of industrial water consumption showed a fluctuating downward trend, with a decline rate of 47.05%. The reason was that large enterprises or factories attach importance to the use of reused water and the recycling of recycled water, and the phenomenon of wasteful water resources was decreasing day by day. The overall emergy ecological footprint of domestic water consumption showed a fluctuating upward trend, with an increase rate of 50.79%. The reason for this was that the population in Yunnan Province is increasing, human activities need to consume a large amount of surface water or groundwater resources, and the domestic water consumption of residents has increased significantly. The emergy ecological footprint of ecological environment water consumption had the smallest proportion, but increased year by year, from 5.61 × 104 hm2 in 2012 to 2.67 × 105 hm2 in 2021. It reflected that Yunnan Province was striving to improve the urban and rural water ecological environment and green sustainable development. The emergy ecological footprint of water pollution increased from 3.25 × 106 hm2 in 2012 to 3.39 × 106 hm2 in 2021, with a growth rate of 4.03%, indicating that the discharge of wastewater in Yunnan Province showed an increasing trend, and water pollution management should be paid high attention. The emergy ecological footprint of aquatic products showed a fluctuating upward trend, with an increase rate of 4.96%. This was mainly due to the increase in regional population and the demand for aquatic products.

Analysis of sustainable utilization of water resources

According to the ecological deficit/surplus of per capita water resources, ecological pressure index of water resources, sustainable utilization index of water resources in Supplementary Table S4, water resources sustainable utilization analysis for 2012–2021 in Yunnan Province was produced as shown in Figure 4.
Figure 4

Analysis of sustainable utilization of water resources in Yunnan Province from 2012 to 2021.

Figure 4

Analysis of sustainable utilization of water resources in Yunnan Province from 2012 to 2021.

Close modal

From 2012 to 2021, the sustainable utilization index of water resources increased in fluctuation, and the change trend was negatively correlated with the trend of the ecological pressure index of water resources. The water resources in 2016, 2017, 2018, and 2020 were in surplus and weakly sustainable utilization. The ecological pressure index of water resources was 0.9164, 0.7837, 0.7946, and 0.9764, respectively. The sustainable utilization index of water resources was 0.0836, 0.2163, 0.2054, and 0.0236, respectively, indicating that the regional water resources supply in that year can support its consumption status and level, and the emergy ecological carrying capacity of water resources can support the future emergy ecological footprint of water resource growth, and there is no damage to the ability to sustain water resources use in the future. For the rest of the years, the water resources were in an unsustainable state of utilization, and the water resources used were in a dangerous state of ecological deficit. The ecological pressure index of water resources was greater than 1, indicating that the emergy ecological carrying capacity of water resources in that year was not enough to support the emergy ecological footprint of water resources, which has caused damage to the ability to continuously utilize water resources in the future. Among them, 2019 was the year with the largest deficit, with a value of − 0.1201, and it was also the year with the lowest sustainable utilization index, with a value of −0.4387. Because the annual average rainfall in this year was the smallest, the total amount of water resources was significantly smaller than in other years, and the emergy ecological carrying capacity of water resources decreased, but the emergy ecological footprint of water resources was the largest, leading to an increase in ecological pressure on the water resource. As a side reflection of the ecological pressure index of water resources, the sustainable utilization index of water resources and the emergy ecological carrying capacity of water resources are affected by the average annual rainfall. If the average annual rainfall becomes smaller, there may be a water ecological deficit or a low ecological surplus of water resources, which will lead to an increase in the ecological pressure of water resources, thus affecting the sustainable use of water resources in the future.

Analysis of eco-economic coordination of water resources

According to the data of the 2012–2021 eco-economic coordination index of water resources in Supplementary Table S4, the change in the eco-economic coordination index of water resources in Yunnan Province from 2012 to 2021 is shown in Figure 5.
Figure 5

Change in water resources eco-economic coordination index in Yunnan Province from 2012 to 2021.

Figure 5

Change in water resources eco-economic coordination index in Yunnan Province from 2012 to 2021.

Close modal

The eco-economic coordination index of water resources showed a fluctuating up-and-down trend. The fluctuation range of WECI was 1.392–1.414, and the multi-year average was 1.409, which indicated that the coordination between the water resources ecosystem and the socio-economic system was relatively stable. The WECI values for 2013, 2015, and 2020 were 1.414, indicating that the water resources were in the best ecological and economic coordination. The year 2019 had the lowest value of 1.392, which meant that the ecological and economic coordination of water resources was relatively poor. Due to the industrial restructuring and technological progress in Yunnan Province in recent years, the utilization efficiency of water resources has been improved. However, it is still necessary to pay attention to the sustainable utilization state of water resources, improve the ecological and social coordination level of water resources, and avoid the decline of the ecological coordination index.

Analysis of utilization efficiency of water resources

According to the emergy ecological footprint of water resources of 10,000 Yuan GDP in Supplementary Table S4 and the data of GDP and water consumption per 10,000 Yuan of GDP in Yunnan Province from 2012 to 2021 in Supplementary Table S1, the change in water resources utilization and development efficiency in Yunnan Province from 2012 to 2021 is shown in Figure 6.
Figure 6

Changes in water resources utilization and development efficiency in Yunnan Province from 2012 to 2021.

Figure 6

Changes in water resources utilization and development efficiency in Yunnan Province from 2012 to 2021.

Close modal

The emergy ecological footprint of water resources of 10,000 Yuan GDP showed a decreasing trend, from 0.1528 hm2/10,000 Yuan in 2012 to 0.0587 hm2/10,000 Yuan in 2021, with a decrease of about 61.60%. The reason for the decline was that, on the one hand, the GDP of Yunnan Province has been increasing year by year with rapid economic development, rising from 1.03 × 108 million Yuan in 2012 to 2.71 × 108 million Yuan in 2021. On the other hand, with technological progress, the water consumption per 10,000 Yuan of GDP has decreased year by year, from 144 m3/10,000 Yuan in 2012 to 59 m3/10,000 Yuan in 2021. During the period from 2012 to 2019, the emergy ecological footprint of water resources of 10,000 Yuan GDP decreased rapidly. Because of the rapid economic development and population growth in Yunnan Province, the demand for water consumption in each account increased. The stable stage from 2020 to 2021 was due to the progress of science and technology, the efficiency of water resources development and utilization was improved, resulting in a decrease in water consumption per 10,000 Yuan of added value produced. It also benefits from the government's development of a green economy and people's increased awareness of water conservation, thus improving the efficiency of water resources development and utilization. From the side, it reflects that water resources are constrained by internal factors and have limited space for sustainable utilization.

Driving factors analysis and optimization paths of sustainable development of water resources

Through the analysis of the current situation of sustainable development of water resources in Yunnan Province during the study period, it can be seen that water resources are generally in an ecological deficit and unsustainable utilization state. However, the changing trend of the comprehensive evaluation index of water resources sustainability shows that the sustainability of water resources is constrained by rainfall and internal factors, and the future ecological carrying capacity of water resources may not be able to support the social development needs. Therefore, this paper analyzes the internal driving factors of sustainable development of water resources from two aspects, namely natural ecological factors and social factors, and puts forward relevant optimization paths.

Natural ecological factors

Natural ecology mainly affects the sustainability of water resources in three main ways. From the perspective of natural geographical conditions, Yunnan is a mountainous province, which belongs to the mountainous plateau terrain. The landform is complex and diverse. The terrain is characterized by high peaks in the north and low peaks in the south, ranging from the highest point of 6,740 m, the main peak of Meili Snow Mountain, Kawagbo Peak, fell to the lowest point of 76.4 m in Hekou County, showing its significant changes in height. The transition from plateau to mountain gradually becomes gentle, the river valleys are deeply embedded under the plateau surface and developed, and the basins are intertwined and coexist among the mountains, which together constitute a complex and diverse landform. On the one hand, the high mountains and steep gorges are deep, the construction of reservoirs is extremely difficult, and the development and utilization of water resources is difficult, resulting in engineering water shortage. On the other hand, under the effect of gravity, the soil matrix of the mountain system loses a lot, the soil water content decreases rapidly, and the ecological environment is extremely fragile, which affects the carrying capacity of water resources. From the perspective of geological structure conditions, karst vulnerability in karst areas consists of mountain system vulnerability and topographic vulnerability. Due to the extremely rapid loss of surface water in karst areas, it is easy to cause frequent flooding and drought disasters in karst mountainous areas, resulting in seasonal water shortages. Therefore, the unique geological structure also affects the sustainable development of water resources.

From the perspective of the hydrological system, Yunnan Province is a typical karst mountainous terrain, which has a dual hydrogeological structure of surface and underground due to special geological and climatic conditions. Surface water infiltration is severe, making it difficult to form a sustainable water surface and having poor water storage capacity. Yunnan has a large difference in terrain height, which makes the exploitation of surface water resources difficult, and the underground water is deeply buried and extremely unevenly distributed in space. This is an important reason for the difficulty of obtaining water and surface drought in karst mountains. In summary, natural ecology mainly controls the sustainable development of water resources by affecting the strength of water resources carrying capacity.

Social factors

The influence of social factors is mainly reflected in increasing the ecological footprint of water resources, reducing the ecological carrying capacity of water resources during the study period, leading to the decrease of water resources' ecological surplus and the increase of water resources' ecological pressure. Firstly, from 2012 to 2021, Yunnan Province was in a period of rapid urbanization, with the population increasing from 46.59 million to 46.9 million, and GDP increasing from 1,030.947 billion to 2,714.676 billion. The demand for water for living and ecological environment is constantly increasing, and the water consumption of agriculture, aquatic products, and wastewater discharge has been at a high level. Secondly, in terms of land development and utilization, Yunnan Province is located in the karst mountainous area, and the increase in population makes the phenomenon of sloping land reclamation very common, which leads to long-term indiscriminate logging and soil erosion, and the water ecological environment is deteriorating, resulting in the decrease of water resources carrying capacity.

Optimal paths for sustainable development of water resources

In order to guarantee the sustainable development of water resources in Yunnan Province, based on the analysis of the driving factors and the current situation of water resources development and utilization, the following optimization paths are proposed:

  • (1)

    Strengthen policy guidance, change the development concept, and improve water resource development and utilization efficiency

The emergy ecological footprint of water resources increases with socio-economic development, and economic and population growth will definitely consume a large amount of water resources and water products. A good economic development model can improve the efficiency of water resource utilization, reduce water losses, and alleviate the ecological pressure on water resources. Yunnan Province should adhere to the policy of ‘ecological priority and green development’, improve agricultural irrigation technology to reduce agricultural water use, improve wastewater treatment technology to reduce wastewater discharge; reduce the input of water resources in aquatic products production; and improve the recycling rate of water resources to form a sustainable development model of water resources.

  • (2)

    Adjust the supply and demand structure of water resources, optimize the allocation of water resources, and develop a green low-carbon economy

The amount of industrial water and industrial sewage discharge plays an important role in the emergy ecological footprint of water resources. The government should adjust the structure of water resources supply and demand at the macro level, and optimize the allocation of water resources in time and space. Developing low-water consumption industries will help economic growth, upgrading or eliminating high-water consumption industries, adjusting and optimizing industrial structure, and forming a green carbon circular economy model. Through the development and utilization of surface water and groundwater, high-quality development of water conservancy projects to regulate the supply of water resources. At the same time, we should further promote the concept of water saving, implement step water price, improve the consumption structure of water resources, improve the utilization efficiency of water resources, and ensure the sustainable development of urban water resources.

  • (3)

    Prevention and control of water pollution, protection of the ecological environment, and improvement of the water resources systems and policies

Government departments should speed up the establishment and improvement of urban and rural sewage and enterprise sewage treatment project implementation progress, and effectively improve the water environment, to solve the enterprise excessive discharge of sewage and urban and rural domestic sewage excessive discharge, establish and improve the accountability system for water pollution incidents, and resolutely close down enterprises or production lines that seriously exceed the discharge standards. The environmental protection department should establish and improve the legal system of water resources protection, strengthen the comprehensive management of water areas, increase the efforts of pollution control and water environmental protection of important water areas, and protect the good state of surface water and groundwater quality such as rivers, lakes, and reservoirs. We should take the biological-ecological pollution control development path, improve the water environment, and improve the ecological environment of water resources' self-cleaning and repair capabilities.

Through a series of measures and policies, we can improve the ecological environment of water resources, enhance the restoration capacity and cleaning capacity of water resources' ecological environment, and accelerate the formation of a green development pattern that matches the ecological footprint of water resources' energy value with the ecological carrying capacity of water resources, and form a situation of harmonious coexistence between human and nature.

Water resources in Yunnan Province were generally in an ecological deficit and unsustainable utilization state during the study period. 2019 was the year with the largest deficit with its value of −0.1201 hm2/cap; 2019 was the year with the highest ecological pressure index of water resources with its value of 1.4387; 2019 was the year with the lowest sustainable utilization index of water resources with its value of −0.4387. In 2019, the emergy ecological carrying capacity of water resources was the lowest, the emergy ecological footprint of water resources was the highest, and the corresponding average annual rainfall was the smallest, with a value of 1,008 mm. During the study period, the ecological pressure index of water resources was relatively high, and the sustainable utilization index of water resources was relatively low. The emergy ecological carrying capacity of water resources was not enough to support the emergy ecological footprint of water resources, which has damaged the ability of sustainable utilization of water resources in the future, and the water resources were in an unsustainable state. Natural ecological factors control the level of sustainable development of water resources mainly by influencing the intensity of the carrying capacity of water resources. The influence of social factors is mainly reflected in the increase of the emergy ecological footprint of water resources in Yunnan Province during the study period, which reduces the emergy ecological carrying capacity of water resources, leading to the decrease of water resources ecological surplus and the increase of water resources ecological pressure. The sustainability of water resources is constrained by rainfall and internal factors, and the future ecological carrying capacity of water resources may not be able to support the social development needs.

During the research period, the emergy ecological carrying capacity of water resources showed a fluctuating increase, with an overall increase rate of approximately 2.08%. Overall, it was still at a relatively low level. The emergy ecological carrying capacity of rainwater was the main contributor to the emergy ecological carrying capacity of water resources, accounting for 81.36% on average. The emergy ecological carrying capacity of surface water and the emergy ecological carrying capacity of groundwater accounted for 13.93 and 4.72%, respectively. The annual change pattern of the emergy ecological carrying capacity of water resources in Yunnan Province is positively correlated with the annual change pattern of the average annual rainfall. The emergy ecological footprint of water resources decreased from 1.72 × 107 to 1.68 × 107 hm2 during the study period, with an overall decrease rate of about 2.31%. Among the ecological footprint accounts, the water consumption account had the greatest impact on the emergy ecological footprint of water resources, with an average contribution rate of 59.59%. The aquatic products account and the water pollution account had relatively small impacts on the emergy ecological footprint of water resources, with average contribution rates of 20.29 and 20.12%, respectively. However, the emergy ecological footprint of water pollution and the emergy ecological footprint of aquatic products showed an upward change, with an upward rate of 4.03 and 4.96%, respectively. During the study period, the emergy ecological carrying capacity of water resources showed an upward trend, the emergy ecological footprint of water resources emergy showed a downward trend, the ecological pressure index of water resources decreased in fluctuation, the sustainable utilization index of water resources increased in fluctuation, and the water resources showed a sustainable development trend. The emergy ecological footprint of water resources of 10,000 Yuan GDP decreased from 0.1528 hm2/10,000 Yuan in 2012 to 0.0587 hm2/10,000 Yuan in 2021, with a decrease of about 61.60%. The efficiency of water resources development and utilization has improved. The eco-economic coordination index of water resources ranged from 1.392 to 1.414, with a multi-year average of 1.409, indicating that the coordination of water resources ecosystem and socio-economic system was relatively stable. In recent years, industrial structure adjustment and scientific and technological progress have continuously improved the utilization efficiency of water resources, and water resources have shown a sustainable development trend.

The sustainable development of water resources is influenced by both natural ecological factors and social factors. The natural ecological factors are mainly manifested as serious soil erosion caused by the fragility of the ecological environment in karst mountainous areas, and the highly permeable hydrogeological structure resulting in weak surface water storage capacity and lack of surface water, thus affecting the level of the ecological carrying capacity of the regional water resources. The social factors are mainly manifested by the increase in demand for water products due to the increase in population and the increase in demand for water resources due to the rapid economic development, as well as the deterioration of the water environment due to the overload discharge of wastewater, which makes the ecological footprint of water resources increase and the ecological carrying capacity of water resources decrease. The optimal path for sustainable development of water resources mainly includes strengthening policy guidance, changing development concepts, improving water resources development and utilization rate; adjusting water resources supply and demand structure, optimizing water resources allocation; developing a green low-carbon economy; preventing and control water pollution, protect ecological environment, and improve water resources system policy.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Abedzadeh
S.
,
Roozbahani
A.
&
Heidari
A.
2020
Risk assessment of water resources development plans using fuzzy fault tree analysis
.
Water Resources Management
34
(
8
),
2549
2569
.
Chang
C.
2023
Geochemical Characteristics of Typical Red Weathering Crust in Karst Area Shilin, Yunnan Province
.
University of Science and Technology
,
Beijing
, p.
136
.
Chen
Y.
,
Zhang
S.
,
Zhang
Y.
,
Xu
L.
,
Qu
Z.
,
Song
G.
&
Zhang
J.
2016
Comprehensive assessment and hierarchical management of the sustainable utilization of urban water resources based on catastrophe theory
.
Journal of the Taiwan Institute of Chemical Engineers
60
,
430
437
.
Darko
G.
,
Obiri-Yeboah
S.
,
Takyi
S. A.
,
Amponsah
O.
,
Borquaye
L. S.
,
Amponsah
L. O.
&
Fosu-Mensah
B. Y.
2022
Urbanizing with or without nature: Pollution effects of human activities on water quality of major rivers that drain the Kumasi Metropolis of Ghana
.
Environmental Monitoring and Assessment
194
(
1
),
38
.
Farrokhzadeh
S.
,
Hashemi Monfared
S.
,
Azizyan
G.
,
Sardar Shahraki
A.
,
Ertsen
M.
&
Abraham
E.
2020
Sustainable water resources management in an arid area using a coupled optimization-simulation modeling
.
Water
12
(
3
),
885
.
Flörke
M.
,
Schneider
C.
&
McDonald
R. I.
2018
Water competition between cities and agriculture driven by climate change and urban growth
.
Nature Sustainability
1
(
1
),
51
58
.
Gleeson
T.
,
Cuthbert
M.
,
Ferguson
G.
&
Perrone
D.
2020
Global groundwater sustainability, resources, and systems in the Anthropocene
.
Annual Review of Earth and Planetary Sciences
48
(
1
),
431
463
.
Hojjati-Najafabadi
A.
,
Mansoorianfar
M.
,
Liang
T.
,
Shahin
K.
&
Karimi-Maleh
H.
2022
A review on magnetic sensors for monitoring of hazardous pollutants in water resources
.
Science of the Total Environment
824
,
153844
.
Huang
L.
,
Zhang
W.
,
Jiang
C.
&
Fan
X.
2008
Ecological footprint method in water resources assessment
.
Acta Ecologica Sinica
28 (
3
),
1279
1286
.
Huang
X.
,
Zhou
Y.
,
Yan
W.
&
Fang
G.
2019
Quantifying benefits of ecological water supply based on emergy analysis
.
Advances in Science and Technology of Water Resources
39
(
2
),
12
15, 36
.
Ishaque
W.
,
Mukhtar
M.
&
Tanvir
R.
2023
Pakistan's water resource management: Ensuring water security for sustainable development
.
Frontiers in Environmental Science
11
, 1096747.
Larsen
T. A.
,
Hoffmann
S.
,
Lüthi
C.
,
Truffer
B.
&
Maurer
M.
2016
Emerging solutions to the water challenges of an urbanizing world
.
Science
352
(
6288
),
928
933
.
Li
Z.
&
Zheng
X.
2023
Jointly Promoting Sustainable Development of Water Resources
.
People's Daily
.
Li
H.
,
Wang
Y.
,
Liu
T.
,
Jiang
X.
,
Sun
Y.
,
Chen
L.
&
Zheng
Y.
2021
Assessment on sustainable utilization of marine aquatic resources based on ecological footprint in China
.
Journal of Fisheries Research
43
(
5
),
509
516
.
Li
X.
,
Zhang
Y.
&
Shan
Y.
2022a
Sustainable utilization of water resources in Yulin city based on an emergy ecological footprint model
.
Arid Zone Research
39
(
4
),
1066
1075
.
Li
Y.
,
Yuan
Y.
,
Li
Z.
&
Guo
X.
2022b
Evaluation of sustainable utilization of water resources in Henan province based on spatial-temporal variation of energy water ecological footprint
.
Yellow River
44
(
6
),
100
104, 162
.
Lian
Y.
2014
Study on Modified Energy-Ecological Footprint of Minjiang River Basin and Its Socioeconomic Causes
.
Fujian Agriculture and Forestry University, Fujian
.
Liu
Y.
2023
The United Nations adopts the ‘Water Action Agenda’
.
Ecological Economy
39
(
6
),
1
4
.
Liu
K.
&
Yang
L.
2021
Characteristics of water resources ecological footprint based on emergy theory: Taking Beijing as an example
.
Research of Soil and Water Conservation
28
(
03
),
406
414
.
Lv
C.
2009
Research on Ecological Economic Value of Regional Water Resources Based on Emergy Theory
.
Zhengzhou University, Zhenghzou
.
Meng
F.
&
Wang
Y.
2022
Ecological wisdom contained in the belief in water god along the ancient silk road
.
Applied Ecology and Environmental Research
20
(
4
),
3155
3171
.
Michelle
T. H. V.
,
Jones
E. R.
,
Flörke
M.
,
Franssen
W. H. P.
,
Hanasaki
N.
,
Wada
Y.
&
Yearsley
J. R.
2021
Global water scarcity including surface water quality and expansions of clean water technologies
.
Environmental Research Letters
16
(
2
),
24020
.
Mikulčić
H.
,
Baleta
J.
,
Zhang
Z.
&
Klemeš
J. J.
2023
Sustainable development of energy, water and environmental systems in the changing world
.
Journal of Cleaner Production
390
,
135945
.
Nguyen
T. G.
,
Phan
K. A.
&
Huynh
T. H. N.
2022
Application of integrated-weight water quality index in groundwater quality evaluation
.
Civil Engineering Journal
8
(
11
),
2661
2674
.
Odum
H.
1996
Environmental Accounting: Emergy and Environmental Decision Making
.
John Wiley and Sons
,
New York
.
Shukui
T.
,
Qi
L.
&
Siyu
H.
2022
Spatial-temporal evolution of coupling relationship between land development intensity and resources environment carrying capacity in China
.
Journal of Environmental Management
301
,
113778
.
Tian
P.
,
Wu
H.
,
Yang
T.
,
Jiang
F.
,
Zhang
W.
,
Zhu
Z.
,
Yue
Q.
,
Liu
M.
&
Xu
X.
2021
Evaluation of urban water ecological civilization: A case study of three urban agglomerations in the Yangtze River Economic Belt, China
.
Ecological Indicators
123
,
107351
.
Wackernagel
M.
,
Onisto
L.
,
Bello
P.
,
Callejas Linares
A.
,
Susana López Falfán
I.
,
Méndez García
J.
,
Isabel Suárez Guerrero
A.
&
Guadalupe Suárez Guerrero
M.
1999
National natural capital accounting with the ecological footprint concept
.
Ecological Economics
29
(
3
),
375
390
.
Wang
Y.
&
Li
B.
2021
Sustainable development and utilization of water resources from the perspective of ecological environment protection
. In
2021 Academic Annual Meeting of Liaoning Provincial Water Conservancy Society
,
Shenyang, Liaoning, China
, p.
3
.
Wang
H.
,
Shi
Q.
,
Li
H.
,
Di
D.
,
Li
Z.
&
Jiang
M.
2023
Spatiotemporal evolution of water ecological footprint based on the emergy-spatial autocorrelation method
.
Environmental Science and Pollution Research
30
(
16
),
47844
47860
.
Wu
J.
,
Wang
Z.
&
Dong
L.
2021
Prediction and analysis of water resources demand in Taiyuan city based on principal component analysis and BP neural network
.
Journal of Water Supply: Research and Technology-Aqua
70
(
8
),
1272
1286
.
Yang
J.
2016
Urban Sustainability Assessment Based on the Emergy Ecological Footprint: A Case Study of Chongqing City
.
Chongqing University
.
Zadgaonkar
L. A.
&
Mandavgane
S. A.
2020
Framework for calculating ecological footprint of process industries in local hectares using emergy and LCA approach
.
Clean Technologies and Environmental Policy
22
(
10
),
2207
2221
.
Zhang
Y.
&
Zhang
H.
2013
Analysis of water ecological footprint in Guangxi based on ecosystem services
.
Acta Ecologica Sinica
33
(
13
),
4111
4124
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data