Articles | Volume 26, issue 8
Research article
26 Apr 2022
Research article |  | 26 Apr 2022

Decreased virtual water outflows from the Yellow River basin are increasingly critical to China

Shuang Song, Shuai Wang, Xutong Wu, Yongyuan Huang, and Bojie Fu

Water scarcity is an emerging threat to food security and socio-economic prosperity, and it is crucial to assess crop production response to water scarcity in large river basins. The water footprint, which considers water use in supply chains, provides a powerful tool for assessing the contributions of water resources within a certain region by tracking the volume and structure of virtual water flows. In this study of the structure of the water footprint network from a complexity perspective, we reassessed the significance of water resources for crop services in a large river basin with a severe water shortage – the Yellow River basin (YRB) of China. The temporal increase of the complexity index indicated that the virtual water outflows (VWFs) from the YRB were becoming increasingly critical to China; i.e. the ability of YRB to produce crops boosted the difficulty of its water being replaced by water exporting from other basins. Decomposition of complexity suggested that during the 1980s to 2000s, the temporally increased complexity was due mainly to the lack of competitors and the increasing uniqueness of crops supporting VWFs. This complexity deeply embedded the YRB into the footprints of a water network that facilitated further development with constrained water resources. Still, it also reinforced reliance from other regions on YRB's scarce water. Based on this analysis, we suggest that resource regulation should be carried out appropriately to ensure ecological sustainability and high-quality development of river basins.

1 Introduction

Water scarcity is an emerging threat to food security and socio-economic prosperity (Zhou et al.2020; Liu et al.2017; Dolan et al.2021; Mekonnen and Hoekstra2016). Water resources play an essential role in producing crops within a river basin. The export of those crops to other river basins can greatly exacerbate the negative consequences of water shortages. Although China's per capita water resources are only 25 % of the global average, China is trying to decrease its dependence on the imported crop, become self-sufficient in food production, and ensure that its food supply is secure. A reasonable assessment of the contribution of the water resources in a river basin to domestic crops supplies will be the first step in balancing the water–food nexus (Wang et al.2020). The Yellow River basin (YRB) is China's most critical agricultural production area. Although the YRB uses only 2.66 % of the water resources in China, the provinces it flows through account for over 35.63 % of the national grain production (data from, official statistics, last access: 14 April 2022), with 41 % of its grain production consumed outside the basin (Zhuo et al.2016b, 2020). Because water consumption for agricultural purposes once accounted for 80 % of the natural runoff of water in the YRB and caused the Yellow River to dry, the supply of water in the YRB is now strictly controlled (Wang et al.2019). Although the availability of water resources constrains agricultural production in the YRB, the national crop yield per unit of water resources has been increasing rapidly (Zhou et al.2020). An assessment of the changing status of agricultural production in the YRB may therefore provide a valuable example to guide the development of water resources in other resource-deficient basins.

The water footprint and virtual water (Conference on Priorities for Water Resources Allocation and Management1993), a geographically explicit indicator that involves water use in supply chains, has provided a powerful tool for assessing the contribution of water resources to a basin and tracking the transfer of water resources across regions (Jaramillo and Destouni2015; Oki and Kanae2004). Virtual water can be associated with specific products that transfer through complex trade relationships between geographic units (Hoekstra2014). Virtual water trade networks consist of virtual water outflows (VWFs) that are embedded in production and consumption trajectories of crops (Oki and Kanae2004; Chini et al.2018; Bae and Dall'erba2018). For example, the VWF from China's water-rich south to its water-scarce north has been quantified, and the network represented by crop trade between provinces has been mapped (Zhai et al.2019; Zhuo et al.2016a). However, because these volume-based studies have ignored the complexity of crop supply and demands (linked with trade) and the uniqueness of each basin (determined by regional characteristics), they have failed to consider the structure of water footprint networks. The structure of a water footprint network reflects the inherent heterogeneity of the distribution of the resources and the pattern of production and consumption. This heterogeneity is consistent with the fact that a water resource cannot be simply replaced by the same volume of water in another basin (Yu and Ding2021; Zhuo et al.2016a; Li et al.2020). This replacement problem reflects that water is not the only resource needed for widely circulated agricultural products. Because of path dependence, others resources, such as the unique hydrothermal conditions and the status of infrastructure within a basin, all determine the position of the basin (or a region) in a water footprint network (Best2019). Complexity, a rapidly developing concept in fields of study such as complexity science, network science, and development economics, represents the “capacity” that is embedded in a certain region, based on measurement of its position in structural terms (Hidalgo2021; Arthur2021; Meng et al.2020; Hidalgo et al.2007). For water resources that are unevenly distributed geographically, complexity can characterize a basin's overall capacity to provide water services required for national crop production. Because complexity-based metrics have been extensively studied in the empirical assessment of economic vitality through a bipartite-networks approach, the concept of complexity should be taken into consideration in the context of water footprint networks to upgrade the toolbox used for integrated water resources management (Hidalgo2021; Hidalgo and Hausmann2009; Liu et al.2017). With such a toolbox, the regional significance of the water supply services provided by the basin through agricultural production can be comprehensively reassessed in a way that takes into account structural factors (Fig. 1).

Recently, the growing idea is that the heart of China's crop production should cease to be the water-poor basin (e.g. the YRB) and that other water-rich basins (southern areas of China) should contribute more (Liu et al.2021; Zhuo et al.2016a). Another call is to ease the water shortage in the YRB by transferring water across the basins as soon as possible (Liu et al.2021). These strategies rely on a comprehensive assessment of the significance of water resources. Thus it would be an oversight to ignore the water footprint network in a structural context because each basin differs in terms of its capacity to use water resources to provide services for crop production (Li et al.2020; Mekonnen and Hoekstra2020, 2011). In this study, we took a complexity perspective concerning water footprint networks and assessed the significance of the YRB to China's crop supply. Our results showed that although the YRB had reduced its virtual water outflow because of resource constraints, its importance to crop production in China increased when considering water footprint networks. The complexity of the water footprint network enabled our approach to provide a new perspective for understanding the changes in the status of a basin with a severe water shortage concerning national crop production.

2 Methods

An overall method route is illustrated here (Fig. 1).

Figure 1Virtual water transfers between regions through multiple dominant crops. When the proportion of VWF from a region is large enough when considering a specific crop, we establish a connection between the region and the crop. In that way, we abstract region–crop bipartite networks for analysis. Here, we give a more straightforward illustration of our method: (1) illustrating a nation comprised of three regions (A, B, and C) where three types of crops (a, b, and c) can be the productions for supporting virtual water flows (VWFs). Region B does not produce crop a, and region C produces a negligible crop a and b. (2) When considering VWF volume, it is hard to compare regions A and B, as they support 40 % national total VWF. (3) However, when the structure is involved, their position in region–crop bipartite networks are different, typically in their differences of diversification (N=1), uniqueness (N=2), and competitiveness (N=3) (the parameter N refers to levels of decomposition, detailed interpretation in Sect. 2.4). (4) Therefore, considering both volume (by ignoring eligible VWF crops) and structure, the GENEPY index is a method for distilling information (see Sect. 2.3). Our main results derive from comparisons between indexes of YRB, China, and random networks (as a benchmark). We generated three different null models by shuffling links in different ways, and Table 2 gives them mathematical descriptions.

2.1 Dataset and VWF estimation

We produced province–crop associations by using a rich dataset on the water footprint of crop production and consumption in China, which is openly accessible (, last access: 14 April 2022). The dataset contains annual statistics for the period 1978–2008 for 22 individual crops (wheat, maize, rice, sorghum, barley, millet, potato, sweet potato, soybean, groundnuts, sunflower, rapeseed, sugar beet, sugar cane, cotton, spinach, tomato, cabbage, apple, grapes, tea, tobacco). For each crop, the annual water footprint of both production and consumption are available at the province level of China (n=31, no data in Taiwan, Hong Kong, and Macau are available). The dataset was widely used when evaluating Chinese virtual water flows (Xie et al.2020; Sun et al.2021; Zhuo et al.2016b). As national water footprint in total production and total consumption is equal for each crop j, the difference between the water footprint of production and consumption indicates the surplus partial for province i, considered to be virtual water outflows (VWFs):

(1) D i j = F i , j production - F i , j consuption , if F i , j production - F i , j consuption > 0 D i j = 0 , if F i , j production - F i , j consuption 0 .

In this way, we obtained the annual volume matrix of VWF Dij in year y, where i and j indicate a specific province and a particular crop; Fproduction and Fconsumption are virtual water embedded in crops production and consumption directly from the original dataset. We focus on outflows, so virtual water inputs (i.e. if consumption is more extensive than its production) are not considered (set to zero).

2.2 Construction of network

To construct a bipartite network, we reduce the VWF volume matrix Dijy to a binary matrix to indicate whether there is a linking edge between a province and a crop (for capturing the network topology). The relative comparative advantage (RCA) procedure was used to construct the bipartite network between provinces and crops. The RCA of production referred to a region's (provinces in this study) export of a particular product (crops, here) in terms of its proportion of the total trade of that product. An analogy can be made to the VWF (Dij) of province i as a proportion of the national total volume of VWF (iDij), embedded in a specific crop j. A higher RCA indicates a larger part of shares among the total national VWF, i.e. more virtual water exported from the province through a specific crop. Firstly, for a specific year y, RCA matrix of province i and crop j are calculated according to (Balassa1965; Dolan et al.2021; Sciarra et al.2020)

(2) RCA i j = D i j / j D i j i D i j / i j D i j .

Then, for y[1978,2008], we constructed such bipartite network RCAij year by year to capture topology changes of Chinese VWF by time. Then, the network matrix M is given by Mij=1 if RCAij≥1, and 0 if RCAij<1.

2.3 Quantitative metrics of complexity

We used a GENEPY index to distil information on the networks in reference to economic complexity (Sciarra et al.2020):

(3) GENEPY i = x = 1 2 λ x X i , x 2 2 + 2 x = 1 2 λ x 2 X i , x 2 ,

where Xi,1 and Xi,2 are the normalized eigenvectors of province i corresponding to the first two largest eigenvalues λ1 and λ2 of the proximity matrix Nii*:

(4) N i i * = j M i j M i * j k i k i * ( k j ) 2 ,  if  i i *  where,  k i = j M i j , k j = i M i j / k i N i i * = 0 ,  if  i = i * .

Here, M is the constructed network matrix, ki is the degree (how many types of crops connected in the network) of the province i, and kj represents the degree of a crop corrected by how easy it is found within the network. The redundant information of the self-proximity (i.e. when i=i*) is deleted by setting corresponding values to zero. In addition, the symmetric square matrix N is interpreted as the mathematical description of the weighted topology of the network – such that the provinces are the nodes and the similarities between the VWF-supporting crops are the links connecting them. Then, the eigenvector centrality of the nodes (i.e. the matrix N) can be a helpful tool to interpret the complexity of the network. In a practical sense, two eigenvalues vectors Xi,1 and Xi,2 of each province i can be combined into unique metrics to distil its ability in producing highly competitive crops. A higher GENEPY index indicates a potential superiority (VWF embedded are more accessible and irreplaceable) in the current VWF networks. For more details in math, we refer the reader to Sciarra et al. (2020). In this way, based on the idea of dimensionality reduction, we could use the average GENEPY index of the YRB to simply assess its importance to the virtual water of China:

(5) GENEPY YRB = 1 n i YRB n ( GENEPY i ) GENEPY China = 1 N i China N ( GENEPY i ) .

Here, iYRB indicates that i is one of the significant provinces that heavily rely on water resources in the YRB, and GENEPYi is the complexity index of the province i according to the Eq. (3).

2.4 Decomposition of complexity

Decomposition where “reflection” applied, of the complexity index “GENEPY” explains the main factors of network changes. The reflection method can describe a structure of bipartite network (Hidalgo and Hausmann2009; Hidalgo2021). To a certain province i, the number of existed connections ki,N varies according to different reflecting times N:

(6) k i , N = 1 k i , 0 j M i j k j , N - 1 , k j , N = 1 k j , 0 i M i j k i , N - 1 ,

where j represents a certain crop and i a certain province, Mij defines the network, and ki,0 represents the observed levels of diversification of a province (the number of products exported by that province). When reflecting time N takes different integer values, it adds further potential connections to the total number in the previous layer (i.e. N−1) of the network. We therefore characterized each province i through the vector ki(ki,1,ki,2,ki,3) in its different dimensions. For example, N=1, with initial conditions given by the degree (i.e. the number of links of provinces and crops): ki,0=ki, same as ki=jMij in Eq. (4).

Table 1Interpretations of the first three pairs of variables describing the province–crop network through the method of reflections.

Download Print Version | Download XLSX

With reference to existing complexity studies, the first three major dimensions can therefore be intuitively explained when N=1, N=2, and N=3 as summarized in Table 1. When the reflection method is used in this way, it reflects the crop diversification, crop uniqueness, and regional competitiveness of the water footprint outflow of the YRB. Although we could have used the reflection approach to continue iterating for more complex explanations, decomposing complexity into three steps helped explain the changes in complexity more clearly and intuitively.

2.5 Null models and sensitivity tests

We randomly created provincial–crop bipartite networks for a sensitivity test. We calculated the same metrics as a comparable reference value to decide whether the influence of networks structure was trivial (Fig. 1). The idea behind the randomization procedure is that we can create a null model starting from the data but shuffling the links of the network while conserving some of its statistical properties. We randomly generated (executed by Python 3.9 and NumPy 1.2) three scenarios where the bipartite networks have the same (1) number of edges, (2) edge sequences on provinces, and (3) edge sequences on crops, with the original dataset, respectively (Table 2 and Fig. 1).

Table 2How different null models were generated.

Download Print Version | Download XLSX

3 Results

3.1 Increasing complexity with decreasing VWF volume

The total volume of virtual water outflows from the provinces in the YRB decreased continuously during the study period (Fig. 2a), and its proportion of the national total decreased by an even more significant percentage (Fig. 2b). In 1978, there were five provinces in the YRB whose virtual water outflows (VWF) exceeded the national average, but there were only three of them in 2008, and their overall ranking had decreased significantly. Though the total volume, share, and ranking of VWF across provinces of the YRB were all decreasing, the average complexity index of the YRB was holistically higher than that of China (Fig. 3a). The gap between the two increased rapidly after 1985. It reached its widest point in 1993 when the average complexity index of the YRB was about 1.4 times the national average (Fig. 3b). After then, the difference between the two decreased with some fluctuations, but the complexity of the YRB remained about 1.2 times the national average (Fig. 3b).

Figure 2Virtual water outflows (VWFs) of the Yellow River basin (YRB) changes from 1978 to 2008. (a) Total VWF in the YRB and China. (b) Proportion of total VWF in the YRB to the national volume. (c) Ranking of VWF in Chinese provinces in 1978. Blue bars are provinces in the YRB. (d) Ranking of VWF in Chinese provinces in 2008. Blue bars are provinces in the YRB.


Figure 3(a) Average complexity index of YRB and of China from 1978 to 2008. (b) Ratios of average complexity index of YRB to that of China from 1978 to 2008. The solid line was fitted with generalized additive models (GAMs). The grey shaded area indicates a 90 % confidence interval. The red dashed line indicates the baseline where the average complexity index of the YRB is equal to that of China.


3.2 Decomposition of changing complexity

Firstly, an indication of crop diversification (N=1, Fig. 4a) was that almost every region was transferring virtual water through 7–8 dominant crops on average. The peak of diversification occurred before about 2000 when the major crop species supporting virtual water transfer were most abundant. Still, there is no significant difference between the YRB, China, and the random network (null models).

When N=2 (Fig. 4b), however, both the YRB and the national average had a significantly (p<0.1) lower number of competitors that transferred virtual water through their similar dominant crops than the random network. The implication is that the crops supporting VWF of the YRB were not all the same to other regions – and are usually more unique than other regions. The opposite nature of the trends of the YRB before the 2000s indicated that its uniqueness was rising and significantly higher (i.e. reflecting index was lower, p<0.1) than the national average. For N=3 (Fig. 4c), whereas the national average of competitiveness was relatively similar to that of the random networks, the total number of types of crops that are competitors with the YRB were significantly lower between 1987 and 2000. The implication is that the competitiveness of the YRB exceeded the average level. In general, since the 1980s, the uniqueness and competitiveness of the YRB have been very different from that of the whole country. However, that difference gradually disappeared in the 2000s.

Figure 4Trends of different dimensions based on the reflecting method decomposed from complexity (see Sect. 2.3 and Table 1). Blue and grey colours show the average levels of the YRB and China, respectively. (a) N=1, number of major crops that support virtual water outflow (VWF), an indication of diversification of a certain region. (b) N=2, the total number of provinces where similar dominant crops supporting VWF are included, a smaller number indicating the greater uniqueness of the particular region. (c) N=3, the total number of dominant crop types supporting VWF from all competitors (who have similar dominant crops). A smaller number of competitors means higher competitiveness of the YRB.


4 Discussion

The use of regional VWF, especially within a water footprint network, has been a common approach to assess the importance of water resources, but complex structures have not been comprehensively evaluated in these assessments (Mekonnen and Hoekstra2020, 2011; Fang et al.2014). The Yellow River, an agriculture-oriented basin with substantial heterogeneity in the upper, middle, and lower reaches, has scarce water resources that are the cornerstone of its crop production (Wang et al.2019). With “The Reform and Opening” of China since 1978, domestic crop trades have gradually increased both in diversity and the volume, and that increase may have caused the total VWF volume to keep increasing. At the same time, however, the YRB was the first basin to apply a water allocation scheme because extreme water shortages caused the river to frequently dry after the 1970s (Wang et al.2019). Therefore, during the time that the national VWF was transferred mainly from the north (where the YRB is located) to the south (Zhuo et al.2016b), the contribution of the YRB, which is deficient in water resources, decreased in volume (Fig. 2). However, considering the structure of the crop–water footprint bipartite network, that did not indicate the YRB had been decreasing in importance to China. On the contrary, our complexity-based analysis revealed the differences in complexity between the YRB and China since 1978. As a result, the YRB's water footprints have been complex for other basins to replace crop production. Our decomposition results suggested that this difficulty was due mainly to a low number of competitors and the increasing uniqueness of crops supporting VWF from the YRB. The Yellow River is a large river that crosses an arid, semi-arid, and monsoon climate zone. The small number of competitors and the relatively high uniqueness of the YRB depend partially on the unique geographical conditions in the YRB (Fu et al.2017). For example, high-quality apples produced in the fertile but arid Loess Plateau in the middle reaches of the YRB and the wet areas in the lower reaches together account for more than 90 % of the total apple production in China regions, which can compete with the YRB. The increasing uniqueness of the YRB means that its ability to improve the quality of its agricultural products and push competitors out of the crop supply network is growing. Those regions that can compete with a particular crop from the YRB (such as Xinjiang, which also exports many high-quality apples) are limited by an extreme shortage of water resources and cannot increase diversification, hence lacking overall competitiveness. Assessing the value of natural resources has been a constant problem because of the trade-off between economic efficiency, social equity, and resource availability (Dalin et al.2015; Grafton et al.2018; Yoon et al.2021). The “Porter hypothesis” of the environment has proposed that environmental regulation can stimulate technological innovation. Similarly, the stimulus of resource scarcity can improve optimization of a market allocation for efficient use of resources (Wagner2004; Luptáčik2010). The complexity of the YRB had risen above the national average since about 1987, when the river basin was regulated. Because of water scarcity, the Chinese government required provinces to adhere to strict resource quotas (Wang et al.2018). During this period, although the total VWF of the YRB and its proportion decreased, the increasing complexity indicated that the resource-constrained YRB was gaining market advantages. In different ways, this advantage has manifested itself as an increase of the position of the YRB in the virtual water network and by an increase of both competitiveness and uniqueness (Fig. 4) (Fang and Chen2015; Fang et al.2014; Yang et al.2012). After the 1990s, the complexity index of China evidenced a similar trend of improvement, which was concomitant with the period when the strict control of water resources and a transformation to water conservation were implemented throughout the whole country (Zhou et al.2020; Liu and Yang2012). According to the literature (Zhou et al.2020), the slowdown of the growth of China's water consumption can be divided into two stages by the 1990s. While growth during the earlier period occurred mainly in arid and semi-arid areas (e.g. the YRB), growth during the latter period affected the whole country. Therefore, the increase of the national average level of complexity lagged behind that of the YRB. This lag was probably the result of the rest of China pursuing structural advantages later than the YRB. Traditional development economic theory points out that the use of resources with comparative advantage is a prerequisite for producing an economically efficient division of labour and a marketing network (Hidalgo and Hausmann2009). In this case, however, before water resources are severely restricted, it may not be more cost-effective to intensify a comparative advantage than to expand resource investment. This conclusion is consistent with the theories related to the “peak water use” and the “efficiency paradox” (Gleick and Palaniappan2010; Grafton et al.2018). Thus, the complexity index may depict the process by which the YRB and then China pursued this structural advantage within the water footprint network as constrained by resources.

In gradually being deeply embedded in the water footprint network, the quality of crops in the YRB has constantly been improved. That improvement has become the structural driving force behind basin development. According to our analysis, the Yellow River basin (YRB) has already developed competitiveness. As a result, the VWF of the YRB was well-structured into national networks. Therefore, it is not easy to replace crop productions with other regions. As a result, water diverted to the Yellow River basin is an option to balance water shortage and crop production with food security (Long et al.2020). However, it must be pointed out that the YRB is a basin with a severe water shortage, and the proposed water resource regulations were also promulgated during a crisis period when the river was drying up (Wang et al.2019, 2018). Like a double-edged sword, the increased complexity facilitated the continued development of the YRB with constrained water resources and embedded it more deeply into the water footprints networks where scarce water could not be easily replaced. From this structural perspective, complexity provided a reminder to help guide China as quickly as possible from resource-dependent development to high-quality development that would enhance regional competitiveness. However, when there is a crisis in the availability of resources, regulations must be implemented. Those regulations may deepen structural issues and make it more challenging to solve the ecological crisis exposed completely.

Our research shows that introducing the concept of complexity into virtual water research is an exciting attempt with new perspectives, but there are still many limitations to this initial exploration. Complexity also depends on other factors in analysed regions, such as policies, markets, infrastructure, and labour resources. However, these are not the focus of this paper and are most similar within provinces. Moreover, that no suitable water footprint data for an extended period to watershed boundaries are available is another challenge for focusing on water resources' influence on complexity. That is why we construct networks by regarding provinces as regional nodes. However, combinations between provincial (Zhuo et al.2016b) and watershed datasets (Xie et al.2020) can be better for accurate assessments. This study focuses more on changing trends rather than the mechanisms. Developing a more sophisticated dataset for further analysis can be one of the future directions. Furthermore, there is also a gap between VWF complexity and its influence on water resources management. If complexity (levels of a region embedded into VWF networks) is good or dangerous to water-poor and water-rich basins is still on an open discussion.

5 Conclusions

This paper took the structure of the water footprint network from a complexity perspective. It assessed the significance of water resources for crop services in a large river basin with a severe water shortage – the Yellow River basin (YRB). From 1978 to 2008, the number of virtual water flows (VWFs) from the YRB and its percentage of the total VWF from the national total decreased significantly. The fact that the YRB has lagged behind the national VWF trend is probably related to restricted water use policies. However, our results showed that the complexity of the YRB increased and was significantly higher than the national average during the period from the 1980s to the 2000s. Decomposition of complexity suggested that this pattern was due mainly to few competitors and the increasing uniqueness of supporting VWF crops for the YRB. Based on an assessment of water conservation policies in China, we suggested that the initial promulgation of resource regulations was the key to a competitiveness-oriented transformation for crop production in the YRB. Subsequently, the increased complexity enabled the YRB to develop with limited water resources. However, it also more deeply embedded the YRB into water footprint networks where scarce water cannot be easily replaced. From the complexity analysis, we point out that resource regulation should be carried out at an appropriate stage to ensure ecological sustainability and high-quality development of river basins.

Code availability

All code is available in a GitHub repository: (SongshGeo2022).

Data availability

Published data are available from (waterfootprint network2022). We are very grateful to Zhuo et al. (2016a) for their excellent work in providing a rich dataset of virtual water studies on Chinese crops.

Author contributions

SW and BF designed this research. SS performed the research and analysed data. SS and YH wrote the paper. XW and YH revised and polished the manuscript and gave major advice.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank Yinan Xiao (School of Artificial Intelligence, Peking University), who verified our application of the mathematical method.

Financial support

This research has been supported by the National Natural Science Foundation of China (grant no. 42041007) and the Fundamental Research Funds for the Central Universities.

Review statement

This paper was edited by Giuliano Di Baldassarre and reviewed by two anonymous referees.


Arthur, W. B.: Foundations of Complexity Economics, Nature Reviews Physics, 3, 136–145,, 2021. a

Bae, J. and Dall'erba, S.: Crop Production, Export of Virtual Water and Water-saving Strategies in Arizona, Ecol. Econ., 146, 148–156,, 2018. a

Balassa, B.: Trade Liberalisation and “Revealed” Comparative Advantage, Man. Sch., 33, 99–123,, 1965. a

Best, J.: Anthropogenic Stresses on the World's Big Rivers, Nat. Geosci., 12, 7–21,, 2019. a

Chini, C. M., Djehdian, L. A., Lubega, W. N., and Stillwell, A. S.: Virtual Water Transfers of the US Electric Grid, Nature Energy, 3, 1115–1123,, 2018. a

Conference on Priorities for Water Resources Allocation and Management: Proceedings of the Conference on Priorities for Water Resources Allocation and Management: Natural Resources and Engineering Advisers Conference, Southampton, July 1992, Overseas Development Administration, London, (last access: 14 April 2022), 1993. a

Dalin, C., Qiu, H., Hanasaki, N., Mauzerall, D. L., and Rodriguez-Iturbe, I.: Balancing Water Resource Conservation and Food Security in China, P. Natl. Acad. Sci. USA, 112, 4588–4593,, 2015. a

Dolan, F., Lamontagne, J., Link, R., Hejazi, M., Reed, P., and Edmonds, J.: Evaluating the Economic Impact of Water Scarcity in a Changing World, Nat. Commun., 12, 1915,, 2021. a, b

Fang, D. and Chen, B.: Ecological Network Analysis for a Virtual Water Network, Environ. Sci. Technol., 49, 6722–6730,, 2015. a

Fang, D., Fath, B. D., Chen, B., and Scharler, U. M.: Network Environ Analysis for Socio-Economic Water System, Ecol. Indic., 47, 80–88,, 2014. a, b

Fu, B., Wang, S., Liu, Y., Liu, J., Liang, W., and Miao, C.: Hydrogeomorphic Ecosystem Responses to Natural and Anthropogenic Changes in the Loess Plateau of China, Annu. Rev. Earth Pl. Sc., 45, 223–243,, 2017. a

Gleick, P. H. and Palaniappan, M.: Peak Water Limits to Freshwater Withdrawal and Use, P. Natl. Acad. Sci. USA, 107, 11155–11162,, 2010. a

Grafton, R. Q., Williams, J., Perry, C. J., Molle, F., Ringler, C., Steduto, P., Udall, B., Wheeler, S. A., Wang, Y., Garrick, D., and Allen, R. G.: The Paradox of Irrigation Efficiency, Science, 361, 748–750,, 2018. a, b

Hidalgo, C. A.: Economic Complexity Theory and Applications, Nature Reviews Physics, 3, 92–113,, 2021. a, b, c

Hidalgo, C. A. and Hausmann, R.: The Building Blocks of Economic Complexity, P. Natl. Acad. Sci. USA, 106, 10570–10575,, 2009. a, b, c

Hidalgo, C. A., Klinger, B., Barabási, A.-L., and Hausmann, R.: The Product Space Conditions the Development of Nations, Science, 317, 482–487,, 2007. a

Hoekstra, A. Y.: Water Scarcity Challenges to Business, Nature Climate Change, 4, 318–320,, 2014. a

Jaramillo, F. and Destouni, G.: Local Flow Regulation and Irrigation Raise Global Human Water Consumption and Footprint, Science, 350, 1248–1251,, 2015. a

Li, M., Wiedmann, T., Liu, J., Wang, Y., Hu, Y., Zhang, Z., and Hadjikakou, M.: Exploring Consumption-Based Planetary Boundary Indicators: An Absolute Water Footprinting Assessment of Chinese Provinces and Cities, Water Res., 184, 116163,, 2020. a, b

Liu, J. and Yang, W.: Water Sustainability for China and Beyond, Science, 337, 649–650,, 2012. a

Liu, J., Yang, H., Gosling, S. N., Kummu, M., Flörke, M., Pfister, S., Hanasaki, N., Wada, Y., Zhang, X., Zheng, C., Alcamo, J., and Oki, T.: Water Scarcity Assessments in the Past, Present, and Future, Earth's Future, 5, 545–559,, 2017. a, b

Liu, Y., Shi, M., Shen, D., Shao, S., Deng, H., and Liao, Y.: Water use and coordinated regional development, Regional Economic Review, 20–31, 2021. a, b

Long, D., Yang, W., Scanlon, B. R., Zhao, J., Liu, D., Burek, P., Pan, Y., You, L., and Wada, Y.: South-to-North Water Diversion Stabilizing Beijing's Groundwater Levels, Nat. Commun., 11, 3665,, 2020. a

Luptáčik, M.: Scarcity and Efficiency, in: Mathematical Optimization and Economic Analysis, edited by: Luptácik, M., Springer Optimization and Its Applications, Springer, New York, NY, 3–24,, 2010. a

Mekonnen, M. M. and Hoekstra, A. Y.: The green, blue and grey water footprint of crops and derived crop products, Hydrol. Earth Syst. Sci., 15, 1577–1600,, 2011. a, b

Mekonnen, M. M. and Hoekstra, A. Y.: Four Billion People Facing Severe Water Scarcity, Science Advances, 2, e1500323,, 2016. a

Mekonnen, M. M. and Hoekstra, A. Y.: Blue Water Footprint Linked to National Consumption and International Trade Is Unsustainable, Nature Food, 1, 792–800,, 2020. a, b

Meng, J., Fan, J., Ludescher, J., Ankit, A., Chen, X., Bunde, A., Kurths, J., and Schellnhuber, H. J.: Complexity Based Approach for El Nino Magnitude Forecasting before the “Spring Predictability Barrier”, P. Natl. Acad. Sci. USA, 117, 177–183,, 2020. a

Oki, T. and Kanae, S.: Virtual Water Trade and World Water Resources, Water Sci. Technol., 49, 203–209,, 2004. a, b

Sciarra, C., Chiarotti, G., Ridolfi, L., and Laio, F.: Reconciling Contrasting Views on Economic Complexity, Nat. Commun., 11, 3352,, 2020. a, b, c

SongshGeo: complexity_yrb, GitHub,, last access: 14 April 2022. a

Sun, J., Yin, Y., Sun, S., Wang, Y., Yu, X., and Yan, K.: Review on Research Status of Virtual Water: The Perspective of Accounting Methods, Impact Assessment and Limitations, Agr. Water Manage., 243, 106407,, 2021. a

Wagner, M.: The Porter Hypothesis Revisited: A Literature Review of Theoretical Models and Empirical Tests, Tech. Rep. 0407014, University Library of Munich, Germany, (last access: 14 April 2022), 2004. a

Wang, Y., Peng, S., Jiang, G., and Fang, H.: Thirty Years of the Yellow River Water Allocation Scheme and Future Prospect, MATEC Web Conf., 246, 01083,, 2018. a, b

Wang, Y., Zhao, W., Wang, S., Feng, X., and Liu, Y.: Yellow River Water Rebalanced by Human Regulation, Sci. Rep.-UK, 9, 9707,, 2019. a, b, c, d

Wang, Z., Xia, J., Zhou, M., Deng, S., and Li, T.: Modelling Hyperconcentrated Floods in the Middle Yellow River Using an Improved River Network Model, Catena, 190, 104544,, 2020. a

waterfootprint network:, last access: 14 April 2022. a

Xie, P., Zhuo, L., Yang, X., Huang, H., Gao, X., and Wu, P.: Spatial-Temporal Variations in Blue and Green Water Resources, Water Footprints and Water Scarcities in a Large River Basin: A Case for the Yellow River Basin, J. Hydrol., 590, 125222,, 2020. a, b

Yang, Z., Mao, X., Zhao, X., and Chen, B.: Ecological Network Analysis on Global Virtual Water Trade, Environ. Sci. Technol., 46, 1796–1803,, 2012. a

Yoon, J., Klassert, C., Selby, P., Lachaut, T., Knox, S., Avisse, N., Harou, J., Tilmant, A., Klauer, B., Mustafa, D., Sigel, K., Talozi, S., Gawel, E., Medellín-Azuara, J., Bataineh, B., Zhang, H., and Gorelick, S. M.: A Coupled Human–Natural System Analysis of Freshwater Security under Climate and Population Change, P. Natl. Acad. Sci. USA, 118, e2020431118,, 2021. a

Yu, D. and Ding, T.: Assessment on the Flow and Vulnerability of Water Footprint Network of Beijing City, China, Journal of Cleaner Production, 293, 126126,, 2021. a

Zhai, M., Huang, G., Liu, L., Xu, X., and Li, J.: Transfer of Virtual Water Embodied in Food: A New Perspective, Sci. Total Environ., 659, 872–883,, 2019.  a

Zhou, F., Bo, Y., Ciais, P., Dumas, P., Tang, Q., Wang, X., Liu, J., Zheng, C., Polcher, J., Yin, Z., Guimberteau, M., Peng, S., Ottle, C., Zhao, X., Zhao, J., Tan, Q., Chen, L., Shen, H., Yang, H., Piao, S., Wang, H., and Wada, Y.: Deceleration of China's Human Water Use and Its Key Drivers, P. Natl. Acad. Sci. USA, 117, 7702–7711,, 2020. a, b, c, d

Zhuo, L., Mekonnen, M. M., and Hoekstra, A. Y.: The Effect of Inter-Annual Variability of Consumption, Production, Trade and Climate on Crop-Related Green and Blue Water Footprints and Inter-Regional Virtual Water Trade: A Study for China (1978–2008), Water Res., 94, 73–85,, 2016a. a, b, c, d

Zhuo, L., Mekonnen, M. M., Hoekstra, A. Y., and Wada, Y.: Inter- and Intra-Annual Variation of Water Footprint of Crops and Blue Water Scarcity in the Yellow River Basin (1961–2009), Adv. Water Resour., 87, 29–41,, 2016b. a, b, c, d

Zhuo, L., Li, M., Wu, P., Huang, H., and Liu, Y.: Assessment of crop related physical-virtual water coupling flows and driving forces in Yellow River basin, J. Hydraul. Eng., 51, 1059–1069,, 2020. a

Short summary
A reasonable assessment of the contribution of the water resources in a river basin to domestic crops supplies will be the first step in balancing the water–food nexus. Our results showed that although the Yellow River basin had reduced its virtual water outflow, its importance to crop production in China had been increasing when water footprint networks were considered. Our complexity-based approach provides a new perspective for understanding changes in a basin with a severe water shortage.