Regional Water Footprint Assessment: a Case Study of Leshan City

This paper presents an assessment of urban water footprint in the period of 2001 to 2012 by taking Leshan City, China as a typical case study. The water footprint is calculated by the sum of the water footprints of various sectors, i.e., crop production, animal products, industrial processes, domestic waster, eco-environment, and virtual water trade. Results show that the water footprints of the various sectors rose by degrees varying from 19% to 55%, which gave rise to an increase of the total water footprint of 43.13% from 2001 to 2012. Crop production and animal products are identified as the major water intensive sectors, accounting for about 68.97% of the total water footprint. The water footprint in the Northeastern area of Leshan City is greater than that of the Southwestern area in the period 1992–2012, resulted in an expansion of water footprint in the Sha Wan and Wu Tongqiao Districts due to the development of urbanization. The application of water footprint assessment is expected to provide insight into the improvement of urban water efficiency, and thus aid in better water resources management.


Introduction
With the development of urbanization, water scarcity has become more and more prominent [1][2][3].In particular, the shortage of water resources has had a great impact on China's economic and social development, i.e., whether sufficient food can be provided to feed such a large and growing population [4][5][6][7][8].To promote water conservation, secure water efficiency corresponding to urban growth is an increasing challenge for China's policy makers [6,9].
The water footprint, as a comprehensive indicator for water resources management, is interpreted as the total amount of freshwater use to measure water consumption and pollution [10][11][12][13].Proposed by Hoekstra in 2003, assessment of water footprint has been widely applied to various categories in the past decade, e.g., products, business services, geographical areas, etc. [7,12,14].
In the field of product assessment, Chapagain et al. [15] assessed its associated impacts on the water resources based upon the water footprint of cotton consumption, which were transitional.Ridoutt et al. [16] applied a LCA-based methodology to water footprint accounting of agricultural food products, by which the farming stage was identified as the major sector of water consumption, and was strongly influenced by the downstream food manufacturers.Similar results were found by Ercin et al. [17], Gerbens-Leenes and Hoekstra [18], Ridoutt et al. [19], Van Oel and Hoekstra [20], Sustainability 2015, 7, 16532-16547 cities in China.In this context, the area with the largest water consumption as well as the time variation of the utilization of water resources of a specific area can be identified, in order to lay out reasonable policy mechanisms for water resource allocation.In addition, a computational tool has been devised to help policy makers better understand and employ the method of water footprint assessment.

Methodology and Data
The total water footprint in a geographically delineated area is seen as a sum of the water footprint accounting process [12].The actual water consumption of one region equals to the local demand of water resources and imported virtual water, in which the former is mainly divided into the demand for crop production, animal products, industrial processes, domestic water and eco-environment, as shown in Figure 1.Other sectors related to the water consumption are omitted in the system boundary.For example, water withdrawal and purification for domestic use are incorporated into the water footprint of domestic water, which is measured as direct water consumption.This is because the domestic sewage should be treated to discharge into nearby rivers, according to the National Standard of the second-class water quality [43].
Sustainability 2015, 7, page-page cities in China.In this context, the area with the largest water consumption as well as the time variation of the utilization of water resources of a specific area can be identified, in order to lay out reasonable policy mechanisms for water resource allocation.In addition, a computational tool has been devised to help policy makers better understand and employ the method of water footprint assessment.

Methodology and Data
The total water footprint in a geographically delineated area is seen as a sum of the water footprint accounting process [12].The actual water consumption of one region equals to the local demand of water resources and imported virtual water, in which the former is mainly divided into the demand for crop production, animal products, industrial processes, domestic water and ecoenvironment, as shown in Figure 1.Other sectors related to the water consumption are omitted in the system boundary.For example, water withdrawal and purification for domestic use are incorporated into the water footprint of domestic water, which is measured as direct water consumption.This is because the domestic sewage should be treated to discharge into nearby rivers, according to the National Standard of the second-class water quality [43].

Water Footprint of Crop Production
The average water demand of crop production mainly depends on the types of crops, the region of growth, and the mode of irrigation system.For a specific crop, the average water demand is calculated per district, given as follows [44]: where is the average water demand (m 3 per tonne) of crop in district of a specific urban area, the water requirement of crop production (m 3 per ha), and the crop yield (tonne per ha).
Here, CWR can be approximately substituted by the accumulated crop evapo-transpiration ET (in mm/day), given as follows [45]: where is the crop coefficient, and is the crop evapotranspiration in an ideal environment of growth ( / ), which is calculated by using the Food and Agriculture Organization (FAO) Penman-Monteith equation [44,46,47]: System boundary of regional water footprint assessment.

Water Footprint of Crop Production
The average water demand of crop production mainly depends on the types of crops, the region of growth, and the mode of irrigation system.For a specific crop, the average water demand is calculated per district, given as follows [44]: where AWD cd is the average water demand (m 3 per tonne) of crop c in district d of a specific urban area, CWR cd the water requirement of crop production (m 3 per ha), and A cd the crop yield (tonne per ha).
Here, CWR cd can be approximately substituted by the accumulated crop evapo-transpiration ET c (in mm/day), given as follows [45]: where K c is the crop coefficient, and ET 0 is the crop evapotranspiration in an ideal environment of growth (mm{d), which is calculated by using the Food and Agriculture Organization (FAO) Penman-Monteith equation [44,46,47]: where: R n = net radiation at the crop surface, MJ{ `m2 d ˘; G = flux of soil heat, MJ{ `m2 d ˘; T = average temperature, ˝C; u 2 = wind speed while measuring at 2 meters height, m{s; e s = pressure of saturation vapour, kPa; e a = actual vapour pressure, kPa; ∆ = slope of the curve between saturation vapour pressure and temperature, kPa{ ˝C; γ = psychrometric constant, kPa{ ˝C.

Water Footprint of Animal Products
The water footprint of animal products is composed of the virtual water content of a live animal during its entire lifespan, and the water consumption while distributing various meat products from the animals [45].The virtual water contains the water use of feeding, drinking water and servicing, calculated as follows [48]: where VWC feed , VWC drink , and VWC servicing denote the virtual water contents corresponding to feeding, drinking and servicing respectively, m 3 per tonne.VWC feed can be measured as the sum of the water requirement of the prepared feed mix and the virtual water of various feed ingredients contained, shown as follows: where: Q mix = water demand of mixing the feed, m 3 {d; VWC i = virtual water content of the ith feed crop, m 3 {t; C i = quantity of feed crop consumed by the animal daily, t{d; W = average live weight of the animal at the end of its lifespan, t; The virtual water content from drinking water VWC drink can be calculated by: VWC drink " where: Q d = the daily drinking water consumed by the animal, m 3 {d; W = average live weight of the animal at the end of its lifespan, t; The virtual water content from service water is deemed as the water consumed by the farmyard cleaning and maintaining, animal washing during the entire lifespan of the live animal, which is calculated by: VWC servicing " where: Q s = the daily service water requirement of the animal, m 3 {d; W = average live weight of the animal at the end of its lifespan, t;

Water Footprint of Industrial Processes
The crop and animal products as raw or auxiliary materials for industrial processes have not been taken into account in the calculation of water footprint of industrial sector, except for direct freshwater use, virtual water content, and water pollution [49].Thus, the water footprint of the industrial sector can be measured by: where V Iblue , V Igrey , V Ivirtual represent the blue, grey and virtual water footprint of the industrial processes.
The blue water footprint is regarded as the direct freshwater use in the industrial processes, calculated by [49]: where V draft indicates the water withdrawn by the stages related to production, transportation etc., m 3 {y; V effluent the waste water emissions, m 3 {y.
The grey water footprint is used to indicate the degree of water pollution, defined as the volume of freshwater which is required to dilute pollutants, calculated by [12]: where W l is the load of pollutant, tonne{y; C max the maximum acceptable concentration, tonne{y; C nat the natural concentration in the receiving water body, tonne{y.In its real application, the industrial waste water is treated by compliance with the National Standard of the third-class water quality standard, and then discharged into nearby rivers, no matter what treatment process is implemented [50].Thus, different cases can be only reflected by selecting the most typical pollutants of one industry, or the pollutants which have the most impact on the local water quality, to calculate the C max and C nat .
In the industrial sector, the most typical pollutants, e.g., COD, BOD 5 , ammonia nitrogen etc., are selected to reflect the degree of water pollution in the delineated region [49].Thus, Equation ( 10) is transformed as: where: V effl = the waste water emissions, m 3 {y; C effl = the concentration of the typical pollutants, tonne{y; C max = the maximum acceptable concentration, tonne{y; C nat = the natural concentration in the receiving water body, tonne{y.
The virtual water footprint of the industrial sector is incorporated into the water required by the raw and auxiliary materials, fuels consumption [49], which is measured as follows: where M i represents the consumption of the ith materials or fuels and k i the water footprint coefficient of the ith materials or fuels.

Water Footprint of Domestic Water
The urban domestic water is regarded as a fundamental source for provision of goods and services consumed by the inhabitants' daily life, including the water demand for household use, personal hygiene, drinking, washing clothes and dishes, flushing toilets, etc. [11,12].The water footprint of urban domestic water is measured as the direct water consumption in this study, by which the water withdrawal and purification for domestic use are incorporated.

Water Footprint of Eco-Environment
The eco-environmental water demand is defined as the amount of water used by the ecosystem to maintain the water balance of the living beings and improve the water environment as well as the environment in which human lives [51][52][53].In this paper, the water footprint of urban eco-environment comprises the ecological water use of urban green spaces, water use of rivers and lakes, and water use of urban sanitation, aiming at the improvement of water quality, ecological environment, and urban landscape etc., shown as follows [51]: where W gr indicates the water footprint of urban green spaces, m 3 {y; W rl the water footprint of urban rivers and lakes, m 3 {y; W sa the water footprint of the urban sanitation, m 3 {y.
The water footprint of urban green spaces W gr is calculated by the quota method, shown as follows [52]: W gr " q gr ˆcgr ( where: q gr = the water quota of urban green spaces, m 3 { `ym 2 ˘; c gr = the urban green coverage, m 2 .
The water footprint of urban rivers and lakes W rl is calculated by the water budget method, shown as follows [54]: where V e is the volume of evaporation from urban water surface, m 3 ; V s the volume of seepage from urban water body, m 3 ; V rl the volume of urban rivers and lakes, m 3 ; f c the period of water replacement.The evaporation from urban water surface V ep is calculated by the following equation [54]: where: A S = area of urban water surface, hm 2 ; E urb = urban evaporation potential, mm; The seepage from urban water body V sp is calculated by the following equation [54]: where: A S = area of urban water surface, hm 2 ; D sp = depth of seepage, mm; The water footprint of the urban sanitation W sa is calculated by the quota method, given as follows [54]: where: q sa = the water quota of urban sanitation, m 3 {hm 2 ; c sa = the urban area, hm 2 .
All of these coefficients, such as D sp , E urb , f c , are derived from the local statistical data.For example, the D sp represents the depth of seepage, which can be calculated by the permeability of sediment on the riverbed [54,55].According to the local standard of hydrogeology and engineering geology, the permeability of sediment in Southwestern China is approximately 0.6m/d.Other empirical values of the coefficients can be determined in a similar way.

Water Footprint of Virtual Water Trade
Virtual water trade is used to indicate the water embedded in traded water intensive products, aiming at improvement of water use efficiency and mitigation of water crisis [10,34].Especially in China, there is a closer connection between commodities trade among different regions or cities, which may result in virtual water trade.The water footprint of urban virtual water trade V trade is calculated by the difference between the gross virtual water export V exp and the gross virtual water import V imp , shown as follows [10,12]:

Case Background and Data Source
Leshan was selected as a typical case study city to assess its associated water footprint.It is a city with a history of more than 3,000 years, located in the southern part of Sichuan Province (102 ˝15'-104 ˝15' E, 28 ˝28'-29 ˝56' N).The administrative area's population is 12,827 including 4 districts, 7 country towns, and 211 village towns, with a total population of 3.56 million in 2013.Leshan's gross domestic production (GDP) has reached 113.479 billion Chinese Yuan, ranked as the 9th in Sichuan Province.As one of the most significant agricultural bases in Sichuan, Leshan City plays an important role in crop and livestock production, ranked the first based on the sales of farm products.In addition, the government of Leshan City pays greater attention to accelerating industrial development, by which a number of high-tech zones have been established in relation to electronics, poly-silicons, pharmaceutical etc [56,57].Leshan city is abundant with water resources with 11.52 billion cubic meters , and located at the confluence of Minjing, Dadu and Qhing Yijiang rivers [58].In this context, the dominant industries are water intensive [59].
In order to assist policy makers in the application of the above methodology for water footprint assessment, as well as to simplify the computational complexity, a calculation tool was developed using the C programming language.Figure 2 shows the user log-in interface of the developed tool.When a user has registered by filling out the request form and agreed on use for academic purposes, the log-on account and the corresponding password are provided by sending them to the user's email address.Once the tool is logged on successfully, calculation of the water footprint of all the above mentioned sectors can be implemented, e.g., on submitting the corresponding data to climate, soil, crops, respectively, the water footprint of crop production can be calculated, as shown in Figure 3.

Results
Table 1 shows the water footprint of each sector in the period 2001-2012, in which the total water footprint of 2012 is 7.79 billion m 3 , increased by 43.13% in contrast to 2001.Crop production accounts for the most at 41.53% on average, animal products 27.44%, virtual water trade 11.74%, Once the tool is logged on successfully, calculation of the water footprint of all the above mentioned sectors can be implemented, e.g., on submitting the corresponding data to climate, soil, crops, respectively, the water footprint of crop production can be calculated, as shown in Figure 3. Once the tool is logged on successfully, calculation of the water footprint of all the above mentioned sectors can be implemented, e.g., on submitting the corresponding data to climate, soil, crops, respectively, the water footprint of crop production can be calculated, as shown in Figure 3.

Results
Table 1 shows the water footprint of each sector in the period 2001-2012, in which the total water footprint of 2012 is 7.79 billion m 3 , increased by 43.13% in contrast to 2001.Crop production accounts for the most at 41.53% on average, animal products 27.44%, virtual water trade 11.74%,

Results
Table 1 shows the water footprint of each sector in the period 2001-2012, in which the total water footprint of 2012 is 7.79 billion m 3 , increased by 43.13% in contrast to 2001.Crop production accounts for the most at 41.53% on average, animal products 27.44%, virtual water trade 11.74%, industrial processes 7.71%, domestic water 6.75%, and eco-environment the least, at 4.82%.This indicates that crop production and animal products are the major water consumed sectors, which highlights Leshan City in its advantage of agricultural production within Sichuan Province [67].Thus, it is suggested that advanced technologies of water-saving irrigation should be developed to enhance the water productivity of crops and minimize loss of water delivery.In spite of the minor fluctuation, the water footprint of each sector presents an underlying trend of growth in the period 2001-2012, as shown in Figure 4.The water footprint of virtual water trade increases substantially, by approximately 55%, followed by industrial processes as 43%, domestic water as 31%, crop production and animal products as 24%.The water footprint of eco-environment is slightly changed, with its increase of 19% up to 2012.For Leshan city, the per capita water resource is 3326 m 3 , which is far greater than the water stress indicator, as the annual water supplies below 1700 m 3 per person [58,68].Since Leshan is rich in water resources, it is targeted to export water intensive products to other cities or regions, such as ceramic and pottery, stainless steel, phosphorus chemical products, etc. [59,69].
Sustainability 2015, 7, page-page 9 industrial processes 7.71%, domestic water 6.75%, and eco-environment the least, at 4.82%.This indicates that crop production and animal products are the major water consumed sectors, which highlights Leshan City in its advantage of agricultural production within Sichuan Province [67].Thus, it is suggested that advanced technologies of water-saving irrigation should be developed to enhance the water productivity of crops and minimize loss of water delivery.In spite of the minor fluctuation, the water footprint of each sector presents an underlying trend of growth in the period 2001-2012, as shown in Figure 4.The water footprint of virtual water trade increases substantially, by approximately 55%, followed by industrial processes as 43%, domestic water as 31%, crop production and animal products as 24%.The water footprint of ecoenvironment is slightly changed, with its increase of 19% up to 2012.For Leshan city, the per capita water resource is 3326 m 3 , which is far greater than the water stress indicator, as the annual water supplies below 1700 m 3 per person [58,68].Since Leshan is rich in water resources, it is targeted to export water intensive products to other cities or regions, such as ceramic and pottery, stainless steel, phosphorus chemical products, etc. [59,69].Dividing the total water footprint by the urban acreage, water footprint intensity (WFI) is thus calculated by using the ArcGIS software to create a water footprint map (WFM).The WFM is intended to visualize the water footprint intensity in order to reflect its variation in a determined spatial and temporal distribution, e.g., the spatial distribution is limited to the administrative area of Leshan City, whilst the temporal ranges from 2001 to 2012.The WFM is based on the "Three Elements of map design", which involves mathematical element, geographic elements and 0 0.5  Dividing the total water footprint by the urban acreage, water footprint intensity (WFI) is thus calculated by using the ArcGIS software to create a water footprint map (WFM).The WFM is intended to visualize the water footprint intensity in order to reflect its variation in a determined spatial and temporal distribution, e.g., the spatial distribution is limited to the administrative area of Leshan City, whilst the temporal ranges from 2001 to 2012.The WFM is based on the "Three Elements of map design", which involves mathematical element, geographic elements and decorative elements [70].The satellite imagery of Leshan City is rasterized to create the basic elements layers, such as administrative division, rivers, etc. Related data, e.g., regional area, water footprints of different sectors etc., are input into the database.By using the anti-gravity weighting method and interpolation, the preliminary WFM is generated through further colour rendering.The annotation is labeled on the map, such as legend, plotting scale, etc.
Figure 5 shows the WFI variation of Leshan City in the period 2001-2012.The ranges are highlighted using a background colour which gradually changes from blue to red, indicated that the water footprint intensity increased from 100 to 120,000.Accordingly, the features of water resources utilization in Leshan City can be summarized as follows from the figure : (1) Variation of spatial distribution: The spatial representation of water usage helpful in identifying the water scarce area.For instance, the WFI of the Northeastern area is apparently larger than that of the Southwestern area, which is consistent with the regional characteristics of demographic and geographic distribution, i.e., the Northeastern area is the comparatively economic developed area in Leshan City, with dense population and water intensive industries.

Discussion
In this study, the sectors division is similar to that has been proposed by Aldaya et al. (2010) and Rui et al. (2011), who place emphasis on evaluating regional water resources utilization from the perspective of factors of production [71,72].From the decomposition of water footprint, the sector of crop production accounts for the largest proportion, at 41.53% averagely.This has been validated by a number of similar studies using different accounting methods, to indicate that the agricultural sector consumes a large amount of freshwater, and contributes to the largest water

Discussion
In this study, the sectors division is similar to that has been proposed by Aldaya et al. (2010) and Rui et al. (2011), who place emphasis on evaluating regional water resources utilization from the perspective of factors of production [71,72].From the decomposition of water footprint, the sector of crop production accounts for the largest proportion, at 41.53% averagely.This has been validated by a number of similar studies using different accounting methods, to indicate that the agricultural sector consumes a large amount of freshwater, and contributes to the largest water footprint [7,37,73,74].At the same time, it has great potential to enhance agricultural water efficiency, by which water saving irrigation is thus strongly encouraged [75].For instance, drip irrigation is proposed by using networks of pipes to water the soil surface or root zone of plants directly, which is a management option to reduce water consumption and maintain agricultural production and productivity [76].Compared with the conventional surface flooding irrigation commonly used in China, drip irrigation is more efficient in reduction of energy consumption and carbon emissions [77].Significant efforts should be also made by policy makers regarding identification of advanced water resource management to promote water saving irrigation, e.g., remote sensing and satellite imagery are useful to identify losses in agricultural productivity by assessing agricultural water use, which supports better agricultural planning [78].Besides, a number of policy instruments, such as pricing strategy, economic sanctions, financial subsidies, etc., may give rise to internalization of the external cost of water, thus to drive development of less water intensive industries and services [7,79].
Water footprint accounting provides better implications to illustrate the potential links among human demand, water consumption and global trade, in order to improve water resource management [40,80].However, the accounting methods of water footprint are still in progress, which may result in uncertainties of the associated results, e.g., lack of a unified methodology for the accounting [81,82].The accounting methods can be mainly divided into the bottom-up and top-down methods in general, in which the former focuses on water footprint assessment in the process of manufacture, whilst the latter emphasizes the product consumption process [12,83].The incomparability of the final accounting result thus gives rise to a variety of data demands and division of systems boundary.In our study, water footprint is calculated by the sum of the water footprints of various sectors, which can be attributed to a typical application of the top-down method, to evaluate regional water resources utilization in a perspective of factors of production.It is expected to confirm the validity and sensitivity of this approach by applying it to a wider range of case scenarios.
Apart from that, it is noteworthy that no water footprints can completely enhance the understanding of water resource issues and contribute to policy analysis, except for bringing increased attention to discussions on water scarcity [84][85][86].A single indicator like a water footprint is not sufficient to provide policy guidance, but more information that involves variables regarding physical, social and economic dimensions should be taken into account in order to determine an optimal framework of policy making related to water resources management [86,87].

Conclusions
This study applies water footprint assessment to measuring the water consumption of Leshan City in the period 1992-2012.The water footprints of various sectors show an increasing trend to varying degrees, in which the virtual water trade increases substantially by approximately 55% up to 2012, followed by the industrial processes at 43%, domestic water at 31%, crop production and animal products at 24%, and eco-environment at 19%.In contrast to 2001, the total water footprint of 2012 increased by 43.13%.Crop production and animal products are identified as the major water consumption sectors, accounting for about 41.53% and 27.44% of the total water footprint.
From the perspective of the geographic distribution of Leshan City, the water footprint in the Northeastern area is greater than that of the Southwestern area in the period 1992-2012.In particular, the distribution of total water footprint has led to an expansion in Sha Wan and Wu Tongqiao Districts in the Northeastern area, with the development of urbanization.
It is expected that the water footprint assessment can provide insight into the improvement of urban water management.For example, crop production and animal products are identified as the major water consumed sectors in the study.Water saving irrigation is thus recommended to enhance the water productivity of crops and minimize loss of water delivery, in order to ultimately reduce the associated water footprint.
However, there are some limitations in this study.First, the specific water requirement from different sources in a designated sector, e.g., the origin of that water, has not been taken into account.Second, the internal structure of water footprint accounting, such as the water balance between blue and green water, has been omitted in this study.Further study will center on improvement of the assessment, in order to examine the internal structure and establish a more precise system boundary for the accounting, to promote sustainable water management.

Figure 1 .
Figure 1.System boundary of regional water footprint assessment.

Figure 2 .
Figure 2. Log-on interface of the calculation tool.

Figure 3 .
Figure 3. Water footprint calculation of crop production by using the calculation tool.

Figure 2 .
Figure 2. Log-on interface of the calculation tool.

Figure 3 .
Figure 3. Water footprint calculation of crop production by using the calculation tool.

Figure 3 .
Figure 3. Water footprint calculation of crop production by using the calculation tool.

Figure 4 .
Figure 4. Variation trend of water footprint for each sector in the period 2001-2012.

Figure 4 .
Figure 4. Variation trend of water footprint for each sector in the period 2001-2012.

Figure 5 .
Figure 5. Variation of water footprint intensity of Leshan City in the period 2001-2012.

Table 1 .
Water footprint of Leshan City in the period 2001-2012 (billion cubic meters).

Table 1 .
Water footprint of Leshan City in the period 2001-2012 (billion cubic meters).