Next Article in Journal
Recent Progress of Electrochemical Energy Devices: Metal Oxide–Carbon Nanocomposites as Materials for Next-Generation Chemical Storage for Renewable Energy
Next Article in Special Issue
Comparative Study of Environmental Assessment Methods in the Evaluation of Resources and Environmental Carrying Capacity—A Case Study in Xinjiang, China
Previous Article in Journal
Exploring Socio-Cognitive Mindfulness in the Context of Sustainable Consumption
Previous Article in Special Issue
Quantitative Assessment of Soil Erosion Based on CSLE and the 2010 National Soil Erosion Survey at Regional Scale in Yunnan Province of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment and Pressure Response Analysis of the Water Footprint of Agriculture and Livestock: A Case Study of the Beijing–Tianjin–Hebei Region in China

1
State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(13), 3693; https://0-doi-org.brum.beds.ac.uk/10.3390/su11133693
Submission received: 18 June 2019 / Revised: 2 July 2019 / Accepted: 4 July 2019 / Published: 5 July 2019
(This article belongs to the Collection Risk Assessment and Management)

Abstract

:
Excessive water consumption, associated with regional agriculture and livestock development and rapid urbanization, has caused significant stress to the ecological health and sustainable use of water resources. We used the water footprint theory to quantify the spatiotemporal characteristics and variation in the water footprint of agriculture and livestock (WF-AL) in the Beijing–Tianjin–Hebei region of China (2000–2016). We predicted the spatial distribution and sustainability of regional water resources at different levels of annual precipitation. Results showed that the average county WF-AL rose from 8.03 × 108 m3 in 2000 to 10.89 × 108 m3 in 2016. There was spatial heterogeneity compared to the average city WF-AL. The WF-AL varied between the mountains and the plains. The scale of the WF-AL was one of the main reasons for differences in the consumption and distribution of water resources. The development of regional water resources deteriorated from a stable state to an unstable state from 2000 to 2016. Only 5.8% of the areas maintained a stable state of water resources. Even in the predicted wet years, no improvements were found in the instability of water resources in four areas centered on the counties of Xinji, Daming, Luannan, and Weichang. To achieve a medium and long-term balance between WF-AL development and water resource recovery, the WF-AL should be limited and combined with reservoir and cross-regional water transfer.

1. Introduction

The consumption of freshwater resources is an integral part of modern society but it has gradually become a limiting factor for society [1]. According to the 2017 United Nations World Water Development Report, 60% of the world’s population live in water-deficient areas, and by 2030 the global freshwater demand will grow by 50% compared with current demand [2,3,4]. China is one of the countries suffering the most severe water shortages worldwide, and problems related to the uneven spatial distribution of water resources and excessive water consumption are extremely serious in this country [5,6,7]. The Beijing-Tianjin-Hebei (BTH) region in China’s semi-arid and semi-humid areas typifies these problems. In 2017, the BTH region had per capita water resources of 286 m3/person, less than 30% of the international minimum standard [8,9]. Urbanization of the BTH region has reached 64.9%, and 2.35% of the land area has to bear 7.24% of the population’s living security and 10.36% of the socio-economic aggregate. The siphon effect of urbanization will increase the scale of regional agriculture and livestock development in the non-core urban areas of the BTH region, further increasing the demand for water resources [10,11]. However, 65% of the surface rivers in the BTH region are facing a shortage of water resources due to seasonal flows, as well as water shortages because of agricultural non-point source pollution and secondary pollution from transferred ecological water [12,13]. Regional agriculture and livestock development has caused significant stress on the ecological health of river water and the sustainable use of water resources in the BTH region [14].
With regard to the scarcity of regional water resources, the governments in the BTH region have adopted water-saving control measures to deal with the impact of water shortages [15,16,17]. However, the control measures focus on cross-regional water resource regulation in the core urban areas, or small-scale agricultural water-saving demonstrations and experiments in the non-core urban areas; no effective measures have been taken to restrict the over-exploitation of water resources in regional agriculture and livestock production [18,19]. Given the current climatic conditions, which are relatively stable, it will be difficult to change the trend of regional water shortages in the short to medium term [20]. It is, therefore, necessary to maintain a dynamic balance between the development of agriculture and livestock production and the conservation and recovery of water resources [21].
The water footprint theory was proposed in 2002. Because of its simple and easy application, this theory has been extensively used in research into the consumption and evaluation of water resources at macro scales, such as national and regional levels [22,23,24,25,26,27]. Researchers have explored the impacts of land use, typical crop production, pollution risk, demographic change, and socioeconomic fluctuations on regional water resources at province and city levels in the BTH region [15,18,28,29,30,31]. However, existing studies rarely involved small-scale horizontal comparisons between counties, which weakens the analysis of spatial differences in the water footprint and affects the feasibility of implementing recommendations. Temporally, previous studies have not analyzed the long-term cumulative effects of changes in the water footprint, and have ignored the analysis of regional water resources under different precipitation conditions. Regional evaluation studies also reduce the effect of natural differences in the distribution of agriculture and livestock development caused by factors such as terrain and urbanization.
In the present study, we estimated the water footprint of agriculture and livestock (WF-AL) in 155 counties (including districts) in the non-core urban areas of the BTH region from 2000 to 2016 based on typical agriculture and livestock production. We then analyzed the distribution of the WF-AL based on precipitation and terrain characteristics across different areas. The main objectives of this study were: (1) to quantitatively evaluate the spatiotemporal variation of the WF-AL in the BTH region in 2000–2016; (2) to estimate the response pressure of agricultural and livestock development on water resources in the BTH region across different levels of annual precipitation; and (3) to provide recommendations and plans for adjusting the distribution of agriculture and livestock industries and improving the allocation of water resources. The results can help to optimize water resource allocation and implement a sustainable medium and long-term development plan for agriculture and livestock in the BTH region. The study also provides a reference for research into the conservation of water resources in similar semi-arid and semi-humid areas.

2. Materials and Methods

2.1. Study Region

The BTH region is located in the Haihe River Basin in northern China. This region is comprised of the Beijing municipality, Tianjin municipality, and Hebei Province, with a total area of approximately 216,000 km2. The terrain is complex, generally tilting from the northwest to the southeast, consisting of a mountain area, a buffer zone (mountain–plain transition zone), a plains area, and a coastal plain (Figure 1). The BTH region has a temperate continental monsoon climate. The average annual precipitation range is between 400 and 600 mm from north to south and the precipitation is mostly concentrated in June–September, which is typical for a semi-arid and semi-humid area.
The BTH region covers 13 cities, from which the specific areas studied here were selected according to the “Beijing-Tianjin-Hebei Collaborative Development Plan” and the “13th Five-Year Development Plan” of each city, as well as the development trends of counties and their differences in industrial and agricultural development (as of 2018). After excluding the core urban areas of each city, a total of 155 county-level administrative areas were selected to form the study region for WF-AL in the BTH region. These included 31 counties in the mountain area, 37 counties in the buffer zone, 80 counties in the plains area, and 7 counties on the coastal plain.
Yield data for major agricultural and livestock products and area data for administrative divisions in the 155 county-level administrative areas from 2000 to 2016 were derived from the relevant statistical yearbooks and socioeconomic statistical bulletins of Beijing, Tianjin, and Hebei (2000–2017) (http://www.nlc.cn/). The data on precipitation and water resource use were derived from province/city water resource bulletins and city statistical yearbooks (2000–2017). Classification of the terrain and watersheds was based on the 30-m digital elevation model of China (2010).

2.2. WF-AL Analysis

2.2.1. Calculation Methods

The water footprint refers to the sum of water resources contained in direct or indirect goods and services consumed by humans per year within a certain area [32]. Agriculture and livestock production and subsequent product consumption are an integral part of human production and livelihoods, as well as one of the most important forms of water consumption [33]. The water footprint of agriculture and livestock (WF-AL) refers to the total water consumption during the life cycle of agricultural production and livestock breeding, in which the life cycle of agricultural and livestock products is generally one year [34]. Currently, there are two main methods for the calculation of WF-AL: the top-down evaluation method and the province/city-level bottom-up accumulation method [35,36].
A third method, the county-level bottom-up accumulation method, is based on county-level sub-units, while the data are dependent on county statistical yearbooks. The county-level method is superior to the two methods mentioned above in terms of its use for the quantification, evaluation, and optimal management of small and medium-sized watershed units across administrative areas. Therefore, the county-level bottom-up accumulation method was adopted to calculate WF-AL and its spatial density in our study region. The calculation methods are expressed as follows:
W F A L = P i × V W
where WF-AL represents the total WF-AL in a region, which quantifies the water consumption of regional agriculture and livestock; Pi represents the mass of an agriculture and livestock product consumed in the corresponding region; and VW represents the virtual water content per unit mass of the agricultural and livestock products in the corresponding region.
ρ [ W F A L ] = W F A L ρ H 2 O × A rea
where ρ[WF-AL] represents the spatial density of WF-AL in a region, which measures the spatial variation of WF-AL per unit area in the region; ρH2O represents the density of water (ρH2O = 1000 kg/m3); and Area represents the area corresponding to the region.
L [ W F A L ] = ρ [ W F A L ] P = W F A L ρ H 2 O × A rea P
where L[WF-AL] represents the difference between the spatial density of WF-AL and the annual precipitation in a region, which measures the development intensity and surplus or loss of water resources for agriculture and livestock in the region and P refers to the annual precipitation in the region. L[WF-AL] < -200mm indicates that the development of regional water resources is in a healthy state, while L[WF-AL] ∈ (−100,100) mm indicates a relatively balanced state and L[WF-AL] > 400mm indicates an unstable state [37].
K [ W F A L ] = ρ [ W F A L ] × 1 P = W F A L ρ H 2 O × A rea × 1 P
where K[WF-AL] represents the ratio between the spatial density of the WF-AL and the annual precipitation in a region, which measures the response relationship between regional agriculture and livestock development and precipitation.

2.2.2. Virtual Water Content of Agriculture and Livestock Products

The WF-AL per unit mass is numerically equivalent to its virtual water content, which is determined by water consumption during the production process of crop planting/livestock rearing. In the non-urban industrial areas, agriculture and livestock products form the most important part of regional water consumption. There are two main methods for calculating the virtual water content of agricultural products (VWc): estimation based on the measured evapotranspiration of crops, and estimation based on the inversion of products. For estimation of a long-term series of virtual water content for agriculture and livestock, the former method is applied to single small watersheds with station observation experiments, while the latter method is applied to watersheds with limited empirical data and a large area.
Limited to the measured meteorological data in the BTH region, the actual evapotranspiration of agricultural products is difficult to obtain. Here, we adopted the inversion of products method to estimate the VWc of agricultural products. We selected agricultural products with a total yield > 5% throughout the study region. These products included eight food crops and cash crops. Their VWc values are listed in Table 1 [38,39,40,41]. The virtual water content of livestock products (VWa) refers to the water consumption in the life cycle of livestock production during feeding. Because large-scale farming is the dominant method for the main livestock products in the study region, we selected six typical livestock products based on the results of similar studies in this region. Their VWa values are given in Table 1 [14,32,42,43,44,45,46,47,48].

3. Results

3.1. Structure of the WF-AL

3.1.1. County WF-AL

From 2000 to 2016, the average annual county WF-AL (WF-ALcounty) increased from 8.03 × 108 m3 to 10.89 × 108 m3. The average annual growth rate of WF-ALcounty was 1.9%, indicating a slow upward trend. However, for individual counties, there were spatiotemporal differences in WF-ALcounty. From 2000 to 2002, the Chongli District of Zhangjiakou had the lowest value (0.39 × 108 m3), while the maximum value was found in the Gaogcheng District of Shijiazhuang (31.69 × 108 m3). In 2015–2016, the Mentougou District of Beijing had the lowest value (0.25 × 108 m3), while the maximum value occurred in the Xinji County of Shijiazhuang (34.65 × 108 m3; Figure 2). The average annual growth rate of the Max(WF-ALcounty) was 0.6%, while the Min(WF-ALcounty) had an average annual growth rate of 4.9%. The ratio of Max(WF-ALcounty) to Min(WF-ALcounty) increased from 81.3 in the early stage to 138.6 in the later stage. In regard to spatial distribution, the top 10% of areas in terms of Max(WF-ALcounty) showed an agglomeration effect. Most were located in the plain areas of Shijiazhuang and Handan, or the coastal plain of Tangshan. By contrast, the top 10% of areas in terms of Min(WF-ALcounty) were scattered in the northwest and central east of the BTH region.
Based on Figure 2, we selected 30% and 70% of the average WF-ALcounty as the reference values [49]. Thus, the WF-ALcounty values 7.13 × 108 m3 and 11.26 × 108 m3 were used as the thresholds for regional low–medium, and medium–high water consumption in the BTH region. Table 2 shows that the proportion of counties with low water consumption gradually decreased from 51.6% to 32.9%, while the proportion of counties with high water consumption doubled from 19.4% to 38.1%. The proportion of counties with medium water consumption showed periodic fluctuations, but generally remained unchanged. Spatially, the low water consumption areas which were located in the mountain area and buffer zone became medium water consumption areas, while the medium water consumption areas located in the buffer zone and plains area became high water consumption areas. The water consumption in the BTH region progressively increased from the northwest to the southeast.

3.1.2. City WF-AL

The average city WF-AL (WF-ALcity) of Tangshan, Shijiazhuang, Handan, and Qinhuangdao was higher than 11.26 × 108 m3, which is in the high water consumption category (Figure 3). By contrast, the average WF-ALcity of Zhangjiakou was lower than 7.13 × 108 m3, which is in the low water consumption category. The remaining eight cities were in the medium water consumption category, although Hengshui and Cangzhou were close to the high water consumption category (i.e., WF-ALcity was generally within the range of medium–high water consumption). From 2000 to 2016, the average WF-ALcity increased from 8.00 × 108 m3 to 11.05 × 108 m3, with an average annual growth rate of 2.0%. The average annual WF-ALcity of the 13 cities fluctuated in the range of 2.92 × 108 m3–18.99 × 108 m3.
After subdividing the water consumption of each city into counties, we found that the WF-ALcounty showed fluctuations consistent with the corresponding WF-ALcity in Tangshan, Hengshui, and Cangzhou (Figure 4). However, there was an increase of WF-ALcounty in local areas or an overall decreasing trend of WF-ALcity found in Zhangjiakou, Baoding, Shijiazhuang, and Handan. Most counties in these four cities were located in the mountain area and buffer zone of the BTH region. This demonstrates that the spatial differences in the WF-ALcity were substantially reduced compared with those in the WF-ALcounty.

3.1.3. Terrain WF-AL

Because the development of agriculture and livestock varies with terrain, there were considerable differences in the planting structure of crops and the type of livestock. The BTH region can be divided into a mountain area, a buffer zone, a plains area, and a coastal plain (Figure 1). The terrain WF-AL (WF-ALterrain) of the buffer zone was generally close to WF-ALcounty, while the corresponding values of the mountain area and the plains area respectively decreased by 21.9% and 7.2%, and the value of the coastal plain markedly increased by 22.4% (Table 3).
Under similar climatic conditions, the land area of agricultural planting in the mountain area and the buffer zone was smaller than that in the plain area and the coastal plain due to limitations imposed by the terrain. In contrast, livestock farming was developed on a larger scale in the mountain area and buffer zone compared with the other two terrains. Thus, there were important differences in the development pattern of agriculture and livestock production between these terrains, and their WF-AL ratios were not comparable. The WF-AL ratios of agriculture and livestock products in the mountain area and plain area were much higher than those in the buffer zone and coastal plain (Figure 5). This is because, compared to agricultural products, livestock products require secondary processing with more complex processes of water consumption. The climatic conditions of the semi-arid and semi-humid regions make them unsuitable for large-scale farming of livestock because of limited water resources.

3.2. Density Variation of the WF-AL

3.2.1. Surplus or Loss of Water Footprint Density

In the northwest and southeast areas of the BTH region (Figure 6), the L[WF-AL] overall remained < 0 from 2000 to 2008, indicating a healthy and balanced state of regional water resource development. However, a fragmented and patchy trend was found in these areas from 2009 to 2016, with L[WF-AL] > 400 mm in some counties (e.g., Zhuolu and Huailai of Zhangjiakou). This indicates that regional water resource development shifted into an unstable state. In the west and central areas of the BTH region, the L[WF-AL] gradually increased from 2000, with periodic fluctuations in local areas. There was no obvious center of surplus or loss for regional water resources, and regional water resource development varied between stable and unstable states.
Long-term significant losses of water resources were found in the north, central southwest, and central east areas of the BTH region. The L[WF-AL] was between 100 and 200 mm in local areas only before 2008, while it was > 400 mm in the remaining areas with the Max[L[WF-AL]] close to 2000 mm. Four areas with water resources in an unstable state were formed between 2015 and 2016, which respectively were centered in Xinji County, Shijiazhuang (L[WF-AL] = 1981 mm) and Daming County, Handan (L[WF-AL] = 1145 mm) in the plain area, Luannan County, Tangshan in the coastal plain (L[WF-AL] = 1847 mm), and Weichang County, Chengde in the mountain area (L[WF-AL] = 618 mm).
Between 2000 and 2016, L[WF-AL] showed small fluctuations within the six intervals between –200 mm and 400 mm, while negatively correlated fluctuations were found within the two intervals of L[WF-AL] < –200 mm and > 400 mm (Table 4). The healthy development areas of water resources sharply decreased during 2000–2014, whereas the unstable development areas rapidly increased in 2003–2005 and 2009–2014. The L[WF-AL] was relatively stable between 2015 and 2016.

3.2.2. Ratio of Water Footprint Density

Owing to the huge terrain differences, microclimate variations are found in local areas of the BTH region, and the precipitation fluctuates substantially at the county level. Therefore, we used the ratio of water footprint density to reduce the spatial differences across the region and quantify the level of regional water resource development and use over time and space (Figure 7).
The proportion of counties with K[WF-AL] < 40% decreased from 12.3% in 2000 to 5.8% in 2008, and then remained almost constant until 2016 (Table 5). The proportion of counties with K[WF-AL] in the range of 40%–70% decreased from 22.6% in 2008 to 9.0% in 2016. No significant fluctuations were found in the proportion of counties with K[WF-AL] in the range of 70%–90% or 90%–110%. However, the proportion of counties with K[WF-AL] > 110% significantly increased from 41.3% in 2000 to 60.6% in 2016, and the increase mostly occurred in the range of >200%.
Taking the development and use rate of river water resources at 40% as the standard threshold for classifying the carrying capacity of river development, only 5.8% of counties in the BTH region consistently met this standard [50,51]. Most of these areas were distributed in the mountain area (Shangyi County of Zhangjiakou, and Laiyuan and Fuping County of Baoding) and the buffer zone (Mentougou and Changping District of Beijing). This is because the counties of Shangyi, Laiyuan, and Fuping have natural groundwater supplies, with nearly 80% forest and grass coverage, complete surface river systems, and stable wetland, lake, and reservoir ecosystems. By contrast, in the Mentougou and Changping Districts, due to the policy restrictions of Beijing, the development of livestock is under strict control, and the areas suitable for agriculture have been used for urban development, which compensates for the limitations of an insufficient surface water supply and leads to a relatively low development and use rate of water resources. The plain area is concentrated in the Jinnan, Xiqing, and Dongli Districts of Tianjin. These three districts have a higher urbanization rate, with limited areas suitable for agricultural development and production. The livestock production in these counties is mainly based on fisheries, which has little impact on water consumption. In addition, these three counties are close to the estuary of Bohai Bay and receive high precipitation. The recent water supply by the Middle Line Project of South-to-North Water Division also contributes to abundant regional water resources in these counties.

3.3. Prediction and Early-Warning of WF-AL

3.3.1. Early-Warning in Different Precipitation Years

We fitted the scenarios for four types of precipitation years (wet, relatively wet, relatively dry, and dry) with the development status of agriculture and livestock in the BTH region between 2015 and 2016 as the baseline (Figure 8) [52]. The simulation was performed without consideration of future control measures such as the transformation of agriculture and livestock or the promotion of water-saving technologies.
In the different precipitation years, the proportions of different terrains with K[WF-AL] < 0.4 were less than 15%. In particular, the proportions of the buffer zone and plains area were less than 5%. Thus, within the range of normal precipitation fluctuations, the development of agriculture and livestock approached or exceeded the threshold for the natural bearing capacity in most of the areas in the BTH region. However, with the K[WF-AL] of 0.4–0.7 and under the precipitation conditions of relatively wet and wet years, the proportions of the mountain area, buffer zone, plains area, and coastal plain that were stable were 19.4%, 10.8%–24.3%, 12.5%–15.0%, and 0.0%, respectively. This suggests that the water ecosystem health in the mountain area and buffer zone could recover in the short term if the precipitation conditions are suitable and the scale of agriculture and livestock is partially limited.
The proportions of mountain area, buffer zone, and plains area with a of K[WF-AL] of 0.7–0.9 increased by 10.0%–13.5% during the conversion from dry years to wet years, while increase of the proportion of the coastal plain with this K[WF-AL] range reached 28.6%. For these areas, the variations in natural precipitation and the scale of agriculture and livestock development have comparable effects on the regional water ecosystems. Thus, the planning of sustainable production of agriculture and livestock should consider different precipitation conditions.
In the wet years, the proportions of mountain area, buffer zone, plains area, and coastal plain with a K[WF-AL] > 0.9 reached 48.4%, 56.8%, 68.8%, and 57.1%, respectively. The mountain area was concentrated in Weichang County of Chengde and Zhuolu County of Zhangjiakou. The buffer zone counties included Luannan County of Tangshan, Dingzhou County of Baoding, and Wu’an County of Handan. The plains area was centered on Gaocheng County of Shijiazhuang and Daming County of Handan. The coastal plain was concentrated in Jinghai District of Tianjin and Changli County of Qinhuangdao. In the abovementioned areas, agriculture and livestock were over-developed, which had considerable negative effects on the enrichment of regional natural water resources. Therefore, while emphasizing the limited production of agriculture and livestock plus the protection of the region’s inherent interests, we should consider using the Middle Line Project of South-to-North Water Division as the main trunk to build and connect the branch channels of cross-regional water transfer in the BTH region, thereby constructing a high-speed water network in this water-deficient region.

3.3.2. Pressure Response Control of Water Resources

Cross-regional water transfer channels are similar to highway networks, and require comprehensive consideration of land resources and construction costs. Large, medium, and small-sized reservoir systems are well-established in the mountain area and buffer zone of the BTH region. However, the overall water conservation function of natural rivers, which are the links between reservoirs, is still weak in these areas. Only some rivers (Jumahe, Baihe, and Luanhe) can function as water transfer channels. The geological conditions are complex, and it is difficult to construct new artificial water transfer channels in a short period of time. Therefore, the adjustment plan for regional development should be based on limiting the scale of agricultural and livestock development and reservoir water transfer.
Compared to the mountain area and buffer zone, it is obviously advantageous to construct a network of new water transfer channels in the plains and coastal plain areas of the BTH region. In addition to relying on natural rivers, the construction of a new artificial water transfer network is feasible in terms of site selection, construction difficulty, and terrain and geological constraints. Especially in severely water-deficient areas, cross-regional water transfer should be combined with limited agriculture and livestock production to mitigate the regional water shortage.
Therefore, taking into account the current situation of K[WF-AL] in different areas of the BTH region, we propose the following plan (Table 6, Figure 9) to establish the medium and long-term balance of agricultural and livestock development and water resource recovery.

4. Discussion

4.1. Merits and Limitations of the Multi-Scale Study

Based on the classification of WF-AL at different scales, the county-level data can reflect the spatiotemporal dynamics of regional WF-AL more directly than the province/city-level data. Counties also serve as grass-root administrative units in China, which are largely autonomous in the implementation of policies and the adjustment and control of agricultural and livestock locations. This system has great advantages for the development and co-ordination of agriculture and livestock production.
This study identified the following limitations: (1) some county-level administrative units were subjected to slight jurisdictional adjustments or mergers, and in these areas, the WF-AL data showed abnormal fluctuations before and after the year of adjustment, and the interval stability of the data was lower compared with the province/city-level data; and (2) it was difficult to obtain the data for virtual water content per unit mass of agricultural and livestock products, so the use of province/city-level average data to some extent affected the accuracy of the study results.

4.2. Optimization of the Ecological Nodes

In the overlap area at the watershed scale, the county-level administrative units were commonly associated with problems such as abnormal WF-AL compared to adjacent areas or a WF-AL much higher than the water resource capacity threshold. Therefore, it was necessary to adopt composite measurements of water resources according to local conditions: (1) counties with water shortages can be supplied by cross-regional water transfers from the Middle Line Project of South-to-North Water Division, the Yellow River-to-Baiyangdian Water Transfer Project, and the Water Diversion Project from Luanhe River to Tianjin City; and (2) in counties where the development intensity of agriculture and livestock is too high, the scale of agriculture and livestock development can be limited or reduced through detailed classification at the village and town levels and delineation of ecologically sensitive areas for water resources within counties.

5. Conclusions

Urbanization is the main direction of development in the urban agglomeration areas of the BTH region. The non-core urban counties of the 13 cities have therefore taken on part of the existing agricultural and livestock production from the core urban areas. The allocation and restoration of regional water resources has become the focus of regional agriculture and livestock development. Based on 14 typical agricultural and livestock products, we calculated the WF-AL and its density in 155 counties of the BTH region over the period 2000–2016. We also predicted the spatial distribution of the WF-AL across different precipitation years. The results showed that in areas with complex terrain conditions, the spatial differences in the WF-AL at the city level were not significant compared with those at county and terrain levels, thus making it difficult to explain the response relationship of water resources to the development of agriculture and livestock. We found that the WF-AL was clearly affected by terrain. Both the ratio of agriculture and livestock and the average value of the WF-AL were lower in the mountain area and plains area than in the buffer zone and coastal plain. Agriculture and livestock production developed rapidly in the BTH region in 2003–2005 and 2012–2014, along with significant expansion of unstable water resource development areas. In 2015–2016, a desirable balance of water resources was maintained only in Shangyi County, Zhangjiakou, Laiyuan and Fuping County, Baoding, Mentougou and Changping District, Beijing, and Jinnan, Xiqing, and Dongli District, Tianjin. Furthermore, the simulation results in different precipitation years showed that in predicted wet years, the regional water resources generally recovered, but there were still four cross-city regional water-deficient areas in the mountain area, buffer zone, plains area, and coastal plain. Therefore, taking into account the current development for agriculture and livestock in 2015–2016, we constructed a medium and long-term plan for sustainable agriculture and livestock production in the BTH region through a combination of limited development of agriculture and livestock, reservoir water transfer, and cross-regional water transfer. The goal is to achieve a balance between the optimal development of agriculture and livestock and the recovery of water resources in the BTH region. This study is of significance in that it implements quantitative evaluation and constructs an optimal plan for agriculture and livestock development in the BTH region at the county level based on the water footprint theory. It is hoped that the results of this study can effectively influence the optimal development of water resource allocation and provide a reference for the management of the city and county-level agriculture and livestock departments in the BTH region.

Author Contributions

Conceptualization, C.C., X.L. (Xiaohan Lu) and X.L. (Xuyong Li); Formal analysis, C.C., X.L. (Xiaohan Lu) and X.L. (Xuyong Li); Investigation, C.C.; Writing—original draft, C.C., X.L. (Xiaohan Lu); Writing—review & editing, C.C., X.L. (Xiaohan Lu) and X.L. (Xuyong Li).

Funding

This research was funded by National Key R&D Program of China, grant number (2016YFC0503007, 2017YFC0505803), National Natural Science Foundation of China, grant number 41771531 and the National Water Pollution Control and Treatment Science and Technology Major Project of China, grant number 2015ZX07203-005.

Acknowledgments

We appreciate the constructive comments of the reviewers and editors on this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoekstra, A.Y.; Mekonnen, M.M. The water footprint of humanity. Proc. Natl. Acad. Sci. USA 2012, 109, 3232–3237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mekonnen, M.M.; Hoekstra, A.Y. Four billion people facing severe water scarcity. Sci. Adv. 2016, 2, e1500323–e1500328. [Google Scholar] [CrossRef]
  3. Ridoutt, B.G.; Pfister, S. A revised approach to water footprinting to make transparent the impacts of consumption and production on global freshwater scarcity. Glob. Environ. Chang. Hum. Policy Dimens. 2010, 20, 113–120. [Google Scholar] [CrossRef]
  4. Miglietta, P.P.; De Leo, F.; Toma, P. Environmental Kuznets curve and the water footprint: An empirical analysis. Water Environ. J. 2017, 31, 20–30. [Google Scholar] [CrossRef]
  5. Zhang, C.; Anadon, L.D. A multi-regional input-output analysis of domestic virtual water trade and provincial water footprint in China. Ecol. Econ. 2014, 100, 159–172. [Google Scholar] [CrossRef]
  6. Xiong, W.; Holman, I.; Lin, E.; Conway, D.; Jiang, J.; Xu, Y.; Li, Y. Climate change, water availability and future cereal production in China. Agric. Ecosyst. Environ. 2010, 135, 58–69. [Google Scholar] [CrossRef]
  7. Jiang, Y. China’s water scarcity. J. Environ. Manag. 2009, 90, 3185–3196. [Google Scholar] [CrossRef]
  8. He, S.; Shao, X. Spatial Clustering and Coupling Coordination of Population-Land-Economic Urbanization in Beijing-Tianjin-Hebei Region. Econ. Geogr. 2018, 38, 95–102. [Google Scholar]
  9. Wu, J.; Xu, D.; Xie, W.; Peng, J. Spatialization of Demographic Data at Medium Scale Based on Remote Sensing Images: Regarding Beijing-Tianjin-Hebei as an Example. Acta Sci. Nat. Univ. Pekin. 2015, 51, 707–717. [Google Scholar]
  10. Han, Y.; Li, X.; Huang, H.; Jia, D. Spatial and temporal distribution of water footprint of main crops and its influencing factors in Beijing-Tianjin-Hebei region. South—North Water Transf. Water Sci. Technol. 2018, 16, 26–34. [Google Scholar]
  11. Li, C.; Xu, M.; Wang, X.; Tan, Q. Spatial analysis of dual-scale water stresses based on water footprint accounting in the Haihe River Basin, China. Ecol. Indic. 2018, 92, 254–267. [Google Scholar] [CrossRef]
  12. Ren, L.L.; Wang, M.R.; Li, C.H.; Zhang, W. Impacts of human activity on river runoff in the northern area of China. J. Hydrol. 2002, 261, 204–217. [Google Scholar] [CrossRef]
  13. Yang, Y.; Tian, F. Abrupt change of runoff and its major driving factors in Haihe River Catchment, China. J. Hydrol. 2009, 374, 373–383. [Google Scholar] [CrossRef]
  14. Mekonnen, M.M.; Hoekstra, A.Y. A Global Assessment of the Water Footprint of Farm Animal Products. Ecosystems 2012, 15, 401–415. [Google Scholar] [CrossRef] [Green Version]
  15. Sun, S. Water footprints in Beijing, Tianjin and Hebei: A perspective from comparisons between urban and rural consumptions in different regions. Sci. Total Environ. 2019, 647, 507–515. [Google Scholar] [CrossRef] [PubMed]
  16. Chu, Y.; Shen, Y.; Yuan, Z. Water footprint of crop production for different crop structures in the Hebei southern plain, North China. Hydrol. Earth Syst. Sci. 2017, 21, 3061–3069. [Google Scholar] [CrossRef] [Green Version]
  17. White, D.J.; Feng, K.; Sun, L.; Hubacek, K. A hydro-economic MRIO analysis of the Haihe River Basin’s water footprint and water stress. Ecol. Model. 2015, 318, 157–167. [Google Scholar] [CrossRef]
  18. Wang, Z.; Liang, L.; Sun, Z.; Wang, X. Spatiotemporal differentiation and the factors influencing urbanization and ecological environment synergistic effects within the Beijing-Tianjin-Hebei urban agglomeration. J. Environ. Manag. 2019, 243, 227–239. [Google Scholar] [CrossRef]
  19. Liu, M.; Xu, X.; Wang, H.; Wang, F. Water Footprint and Spatial-temporal Analysis of Hebei Province Based on Virtual Water Theory. J. Nat. Resour. 2012, 27, 1022–1034. [Google Scholar]
  20. Li, Y.; Zhang, Z.; Shi, M. Restrictive Effects of Water Scarcity on Urban Economic Development in the Beijing-Tianjin-Hebei City Region. Sustainability 2019, 11, 2452. [Google Scholar] [CrossRef]
  21. Fan, L.; Wang, H.; Liu, Z.; Li, N. Quantifying the Relationship between Drought and Water Scarcity Using Copulas: Case Study of Beijing-Tianjin-Hebei Metropolitan Areas in China. Water 2018, 10, 1622. [Google Scholar] [CrossRef]
  22. Serio, F.; Miglietta, P.P.; Lamastra, L.; Ficocelli, S.; Intini, F.; De Leo, F.; De Donno, A. Ground water nitrate contamination and agricultural land use: A grey water footprint perspective in Southern Apulia Region (Italy). Sci. Total Environ. 2018, 645, 1425–1431. [Google Scholar] [CrossRef] [PubMed]
  23. Chapagain, A.K.; Hoekstra, A.Y.; Savenije, H.H.G.; Gautam, R. The water footprint of cotton consumption: An assessment of the impact of worldwide consumption of cotton products on the water resources in the cotton producing countries. Ecol. Econ. 2006, 60, 186–203. [Google Scholar] [CrossRef]
  24. Yano, S.; Hanasaki, N.; Itsubo, N.; Oki, T. Water Scarcity Footprints by Considering the Differences in Water Sources. Sustainability 2015, 7, 9753–9772. [Google Scholar] [CrossRef] [Green Version]
  25. Paterson, W.; Rushforth, R.; Ruddell, B.L.; Konar, M.; Ahams, I.C.; Gironas, J.; Mijic, A.; Mejia, A. Water Footprint of Cities: A Review and Suggestions for Future Research. Sustainability 2015, 7, 8461–8490. [Google Scholar] [CrossRef] [Green Version]
  26. Novoa, V.; Ahumada-Rudolph, R.; Rojas, O.; Saez, K.; de la Barrera, F.; Luis Arumi, J. Understanding agricultural water footprint variability to improve water management in Chile. Sci. Total Environ. 2019, 670, 188–199. [Google Scholar] [CrossRef] [PubMed]
  27. Li, J. Scenario analysis of tourism’s water footprint for China’s Beijing-Tianjin-Hebei region in 2020: Implications for water policy. J. Sustain. Tour. 2018, 26, 127–145. [Google Scholar] [CrossRef]
  28. Li, Y.; Wang, X.; Zhu, C.; Qiao, R. Research for Calculation and Prediction of Water Resource Ecological Footprint in Xingtai City, Hebei Province. Res. Soil Water Conserv. 2014, 21, 227–230. [Google Scholar]
  29. Wang, X.; Li, X.; Fischer, G.; Sun, L.; Tan, M.; Xin, L.; Liang, Z. Impact of the changing area sown to winter wheat on crop water footprint in the North China Plain. Ecol. Indic. 2015, 57, 100–109. [Google Scholar] [CrossRef]
  30. Yongqiang, C.A.O.; Jing, M.A. Empirical Study on Water Footprint of Water Resources Management in Hebei Province. J. Yangtze River Sci. Res. Inst. 2011, 28, 18–21. [Google Scholar]
  31. Han, Y.; Jia, D.; Zhuo, L.; Sauvage, S.; Sanchez-Perez, J.-M.; Huang, H.; Wang, C. Assessing the Water Footprint of Wheat and Maize in Haihe River Basin, Northern China (1956–2015). Water 2018, 10, 867. [Google Scholar] [CrossRef]
  32. Cheng, X.; Sun, R.; Chen, L.; Kong, P. Spatial and temporal patterns of the water footprint in the Beijing-Tianjin-Hebei region with consideration of crop and animal products and domestic water. Acta Ecol. Sin. 2018, 38, 4461–4472. [Google Scholar]
  33. Chapagain, A.K.; Hoekstra, A.Y. The water footprint of coffee and tea consumption in the Netherlands. Ecol. Econ. 2007, 64, 109–118. [Google Scholar] [CrossRef]
  34. Philipp, D.; Putman, B.; Thoma, G. Asas-Csas Annual Meeting Symposium on Water Use Efficiency at The Forage-Animal Interface: Life cycle assessment of forage-based livestock production systems. J. Anim. Sci. 2019, 97, 1865–1873. [Google Scholar] [CrossRef] [PubMed]
  35. Allan, J.A. Virtual water-the water, food, and trade nexus useful concept or misleading metaphor? Water Int. 2003, 28, 106–113. [Google Scholar] [CrossRef]
  36. Wang, Y.Y.; Wang, H.X.; Cai, Y. Calculation and analysis of water footprint in Beijing City. Chin. J. Eco-Agric. 2011, 19, 954–960. [Google Scholar] [CrossRef]
  37. Bao, C.; He, D. Spatiotemporal characteristics of water resources exploitation and policy implications in the Beijing-Tianjin-Hebei Urban Agglomeration. Prog. Geogr. 2017, 36, 58–67. [Google Scholar]
  38. Zhang, Y.; Li, Y.; Ouyang, Z.; Liu, J. The grey water footprint of the winter wheat-summer maize crop rotation system of the North China Plain. Acta Ecol. Sin. 2015, 35, 6647–6654. [Google Scholar]
  39. Wei, Y.; Tang, D.; Ding, Y.; Agoramoorthy, G. Incorporating water consumption into crop water footprint: A case study of China’s South-North Water Diversion Project. Sci. Total Environ. 2016, 545, 601–608. [Google Scholar] [CrossRef]
  40. Mekonnen, M.M.; Hoekstra, A.Y. The green, blue and grey water footprint of crops and derived crop products. Hydrol. Earth Syst. Sci. 2011, 15, 1577–1600. [Google Scholar] [CrossRef] [Green Version]
  41. Wang, Y.; Liu, J.; Zhao, D. Assessing Water Resources Based on Theory of Water Footprint—A Case Study in Xuanhua District, Zhangjiakou City, Hebei Province. Bull. Soil Water Conserv. 2018, 38, 213–219. [Google Scholar]
  42. Kayatz, B.; Harris, F.; Hillier, J.; Adhya, T.; Dalin, C.; Nayak, D.; Green, R.F.; Smith, P.; Dangour, A.D. “More crop per drop”: Exploring India’s cereal water use since 2005. Sci. Total Environ. 2019, 673, 207–217. [Google Scholar] [CrossRef] [PubMed]
  43. Miglietta, P.P.; De Leo, F.; Ruberti, M.; Massari, S. Mealworms for Food: A Water Footprint Perspective. Water 2015, 7, 6190–6203. [Google Scholar] [CrossRef]
  44. Ridoutt, B.G.; Sanguansri, P.; Harper, G.S. Comparing Carbon and Water Footprints for Beef Cattle Production in Southern Australia. Sustainability 2011, 3, 2443–2455. [Google Scholar] [CrossRef] [Green Version]
  45. Ma, D.; Xian, C.; Zhang, J.; Zhang, R.; Ouyang, Z. The Evaluation of Water Footprints and Sustainable Water Utilization in Beijing. Sustainability 2015, 7, 13206–13221. [Google Scholar] [CrossRef] [Green Version]
  46. Ibidhi, R.; Ben Salem, H. Water footprint and economic water productivity of sheep meat at farm scale in humid and semi-arid agro-ecological zones. Small Rumin. Res. 2018, 166, 101–108. [Google Scholar] [CrossRef]
  47. Lamastra, L.; Miglietta, P.P.; Toma, P.; De Leo, F.; Massari, S. Virtual water trade of agri-food products: Evidence from italian-chinese relations. Sci. Total Environ. 2017, 599, 474–482. [Google Scholar] [CrossRef] [PubMed]
  48. Yu, Y.; Zhang, H.; Hu, H. Water Resources Carrying Capacity of Livestock Husbandry in China Based on Water Footprint Theory. Resour. Sci. 2012, 34, 394–400. [Google Scholar]
  49. Lu, Y.; Payen, S.; Ledgard, S.; Luo, J.F.; Ma, L.; Zhang, X.Y. Components of feed affecting water footprint of feedlot dairy farm systems in Northern China. J. Clean Prod. 2018, 183, 208–219. [Google Scholar] [CrossRef]
  50. Shi-feng, Z.; Jun-xu, C. Research on the Risk of Water Resources Shortage in North China. J. Nat. Resour. 2009, 24, 1192–1199. [Google Scholar]
  51. Yuan, W.; Lianxi, S.; Ke, L.I.; Hongyan, S.U.N. Analysis of present situation of water resources and countermeasures for sustainble development in China. J. Water Resour. Water Eng. 2008, 19, 10–14. [Google Scholar]
  52. Xiqin, W.; Yuan, Z. The Allowable Exploitation Rate of Rivers Water Resources of the Seven Major Rivers in China. J. Nat. Resour. 2008, 23, 500–506. [Google Scholar]
Figure 1. Study area and terrain types.
Figure 1. Study area and terrain types.
Sustainability 11 03693 g001
Figure 2. County water footprint of agriculture and livestock production (WF-ALcounty) of the BTH region, 2000–2016.
Figure 2. County water footprint of agriculture and livestock production (WF-ALcounty) of the BTH region, 2000–2016.
Sustainability 11 03693 g002aSustainability 11 03693 g002b
Figure 3. City water footprint of agriculture and livestock (WF-ALcity) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Figure 3. City water footprint of agriculture and livestock (WF-ALcity) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Sustainability 11 03693 g003
Figure 4. Spatial differences within the city water footprint of agriculture and livestock (WF-ALcity) in the Beijing–Tianjin–Hebei (BTH) region.
Figure 4. Spatial differences within the city water footprint of agriculture and livestock (WF-ALcity) in the Beijing–Tianjin–Hebei (BTH) region.
Sustainability 11 03693 g004
Figure 5. Water footprint ratios of agriculture and livestock products in the Beijing–Tianjin–Hebei (BTH) region.
Figure 5. Water footprint ratios of agriculture and livestock products in the Beijing–Tianjin–Hebei (BTH) region.
Sustainability 11 03693 g005
Figure 6. Differences in the water footprint density of agriculture and livestock (L[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Figure 6. Differences in the water footprint density of agriculture and livestock (L[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Sustainability 11 03693 g006
Figure 7. Water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Figure 7. Water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Sustainability 11 03693 g007aSustainability 11 03693 g007b
Figure 8. Water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region under different precipitation years.
Figure 8. Water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region under different precipitation years.
Sustainability 11 03693 g008
Figure 9. A schematic for limited production of agriculture and livestock in the Beijing–Tianjin–Hebei (BTH) region.
Figure 9. A schematic for limited production of agriculture and livestock in the Beijing–Tianjin–Hebei (BTH) region.
Sustainability 11 03693 g009
Table 1. Virtual water consumption of agriculture and livestock products in the Beijing–Tianjin–Hebei (BTH) region (m3/kg).
Table 1. Virtual water consumption of agriculture and livestock products in the Beijing–Tianjin–Hebei (BTH) region (m3/kg).
RegionAgricultural and Livestock Products
Grain CropsCash CropsLivestock Products
WheatMaizePaddyTubersPeanutCottonVegetableFruitPorkBeefMuttonPoultryEggMilk
Hebei1.381.192.191.21.25.50.10.683.618.119.983.508.652.2
Beijing1.230.841.40.71.54.40.20.58
Tianjin1.250.761.190.881.55.220.10.48
Table 2. Proportional distribution of the county water footprint of agriculture and livestock (WF-ALcounty) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Table 2. Proportional distribution of the county water footprint of agriculture and livestock (WF-ALcounty) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
2000–20022003–20052006–20082009–20112012–20142015–2016
Low consumption51.6%44.5%43.2%36.8%30.3%32.9%
Medium consumption29.0%27.1%32.9%28.4%34.2%29.0%
High consumption19.4%28.4%23.9%34.8%35.5%38.1%
Table 3. Distribution of the terrain water footprint of agriculture and livestock (WF-ALterrain) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016 (108 m3).
Table 3. Distribution of the terrain water footprint of agriculture and livestock (WF-ALterrain) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016 (108 m3).
2000–20022003–20052006–20082009–20112012–20142015–2016County Number
Mountain Area3.794.905.446.016.827.5331
Buffer Zone8.639.539.3310.2810.9510.8937
Plain Area9.3610.6210.4411.3611.8911.9180
Coastal Plain8.219.8110.2512.6413.6713.807
Table 4. Proportional distribution of the water footprint density of agriculture and livestock (L[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Table 4. Proportional distribution of the water footprint density of agriculture and livestock (L[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
L[WF-AL] (mm)<−200(−200,−100)(−100,0)(0,100)(100,200)(200,300)(300,400)>400
2000–200227.7%9.7%14.2%14.8%9.7%4.5%4.5%14.8%
2003–200523.9%11.0%8.4%11.6%10.3%9.7%4.5%20.6%
2006–200821.3%10.3%12.3%13.5%8.4%8.4%7.1%18.7%
2009–201116.1%11.6%10.3%12.3%10.3%6.5%7.7%25.2%
2012–20149.0%14.8%10.3%12.3%10.3%9.0%5.2%29.0%
2015–20169.7%12.3%10.3%16.1%9.0%6.5%5.2%31.0%
Table 5. Proportional distribution of the water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
Table 5. Proportional distribution of the water footprint density ratio of agriculture and livestock (K[WF-AL]) in the Beijing–Tianjin–Hebei (BTH) region, 2000–2016.
K[WF-AL]<40%40%–70%70%–90%90%–110%>110%>200%
2000–200212.3%21.3%10.3%14.8%41.3%12.9%
2003–20057.1%21.3%8.4%12.9%50.3%16.8%
2006–20085.8%22.6%9.7%15.5%46.5%17.4%
2009–20115.8%12.9%12.9%11.6%56.8%19.4%
2012–20145.2%11.0%12.9%11.6%59.4%23.9%
2015–20165.8%9.0%11.0%13.5%60.6%25.8%
Table 6. The plan for limited production of agriculture and livestock in the Beijing–Tianjin–Hebei (BTH) region.
Table 6. The plan for limited production of agriculture and livestock in the Beijing–Tianjin–Hebei (BTH) region.
K[WF-AL], NowK[WF-AL], FutureLimited Production of AgricultureLimited Production of LivestockWater Transfer
0–0.4Maintain---
0.4–0.5drop to 0.4-10%-
0.5–0.7Drop to 0.410%30%-
0.7–0.8Drop to 0.7-10%Reservoir water transfer
0.8–0.9Drop to 0.710%10%Reservoir water transfer
0.9–1.1Drop to 0.910%15%Cross-region water transfer
1.1–1.3Drop to 0.910%30%Cross-region water transfer
>1.3Drop to 0.930%30%Cross-region water transfer

Share and Cite

MDPI and ACS Style

Cao, C.; Lu, X.; Li, X. Risk Assessment and Pressure Response Analysis of the Water Footprint of Agriculture and Livestock: A Case Study of the Beijing–Tianjin–Hebei Region in China. Sustainability 2019, 11, 3693. https://0-doi-org.brum.beds.ac.uk/10.3390/su11133693

AMA Style

Cao C, Lu X, Li X. Risk Assessment and Pressure Response Analysis of the Water Footprint of Agriculture and Livestock: A Case Study of the Beijing–Tianjin–Hebei Region in China. Sustainability. 2019; 11(13):3693. https://0-doi-org.brum.beds.ac.uk/10.3390/su11133693

Chicago/Turabian Style

Cao, Chen, Xiaohan Lu, and Xuyong Li. 2019. "Risk Assessment and Pressure Response Analysis of the Water Footprint of Agriculture and Livestock: A Case Study of the Beijing–Tianjin–Hebei Region in China" Sustainability 11, no. 13: 3693. https://0-doi-org.brum.beds.ac.uk/10.3390/su11133693

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop