An assessment of the urban water footprint and blue water scarcity : A case study for Van ( Turkey )

This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Today, most of the world’s population faces water scarcity, while global warming, urbanization, industrialization and population increases continue to increase the severity of the pressure on water resources. Management of water resources plays a key role in the sustainability of agricultural production. The water footprint (WF) is different in comparison to other water statistics because it takes direct and indirect water consumption into account, and helps in the management of water resources. Within this context, the WF of Van province, which is Turkey’s most easterly located arid region, was calculated from 2004 to 2019. The study area covers lake Van, which is Turkey’s largest lake, and the Van basin with an area of 23.334 km2 and a population of 1.136.757 (2019). In the calculations, crop (WFcrop), livestock (WFlivestock), and domestic and industrial water footprints (WFdomestic+industrial) were evaluated separately, and blue and green water footprints (WFblue and WFgreen) were analyzed in detail. According to the results, the average WF of Van province was found to be 8.73 billion m3 year-1. Throughout the province, 87.6% of the WF is composed of WFcrop, 4.9% is WFlivestock and 7.5% is WFdomestic+industrial. Of the WFcrop, 62.5% depends on WFblue, i.e., freshwater. Most of the WFlivestock consisted of dairy cattle (49%) and sheep (38%). The average WFdomestic+industrial for 2004 to 2019 was 0.64 billion m3 year-1. The average per capita water footprint of Van province was found to be 889.9 m3 year-1 capita-1. In addition, the province is classified as severe water scarcity (257%). This study is one of the first province-based calculations of WF in Turkey and is the first study to bring a different aspect to published literature by including residual soil moisture from the winter months. As a result of this study, the WFblue of the WFcrop is above the worldwide average and should be reduced by changing the crop pattern or synchronizing the planting and harvest dates of the crops to a period that benefits from precipitation. In addition, this study is expected to contribute to new studies for calculating the provincial scale WF and will have positive effects on agricultural planning, water allocation and the sustainability of water resources.

2011a; Zhang and Anadon, 2014). However, due to the fact that the agricultural water footprint constitutes a large portion of the total water footprint, that climate and regional conditions play a significant role in agricultural production and water distribution and consumption vary based on regional characteristics Hoekstra, 2014: Zhuo et al., 2014), urban water footprint comes to the fore (Vanham et al. 2014;Paterson et al. 2015;Degefu et al. 2018). Furthermore, since the vast majority of the global population live in cities and the number of urban dwellers continuously increases, the importance of urban water footprint calculations is quite important (Aerts et al., 2009). In addition to, water footprint, which allows numerical analysis, will pioneer the distribution and comparative analysis of water use in these sectors, the green and blue water utilization rates, and the water management and planning stages (Yerli et al., 2019b).
The aim of the present study was to calculate and analyze the production water footprint in Van province covers the Van basin and delivering water resources to different regions, including all the districts of the province. In this context, crop and livestock production, and domestic and industrial water utilization were evaluated separately, and the blue and green water footprint figures were calculated. In adittion, per capita water footprint was also analyzed to compare to other different studies. The main research questions of the study were as follows: What are the water footprints for different production sectors and their share in the province? What is the water footprint and virtual water content of significant crops in the province? What is the effect of water footprint on the sustainability of blue water in arid regions and water consumption for the future? How is the agricultural water footprint affected in arid regions compared to other regions? What is the effect of the trends in the water footprint over time in an area that is constantly changing?
Although Van province provides indispensable production of crops and ranks first in Turkey in terms of small ruminants, there is no comprehensive research regarding water resources. This study is one of the first studies that analyzed the urban water footprint in detail not just in Van but also in Turkey, and it is also the first study contributing to published literature from a different point of view by calculating the water footprint of crop production with high and reliable data inputs that considered the residual soil moisture from winter months. The advantages of the present study compared to other water footprint studies included working in a specific area, the use of up-to-date data accepted by national and international authorities, and benefits from a wide time period (2004 to 2019). In addition, the wide ranging diversity of crops cultivated in the province, their properties and effective root depths, soil properties and water holding capacity of the soils, residual soil moisture from the winter months and the effective precipitation (the average precipitation

Introduction
The freshwater resources on earth have significantly reduced due to global warming caused by anthropogenic effects, such as uncontrolled industrialization and urbanization, excessive fossil fuel consumption and insensitivity to the environment (Cakmak et al., 2007;Song et al., 2018;Awange, 2021). Many studies have pointed out that water will be scarce in the future and accordingly, production will decrease and humanity will face economic problems (Calzadilla et al., 2011;Ercin & Hoekstra 2014;Hoekstra, 2014;Damania, 2020). In addition, the unconscious use of water-increasing water demand -as the population increases, adds to the severity of the pressure on water resources day by day (Sensoy et al., 2007;Sahin et al., 2016). Since the amount of water resources cannot be increased, protection and sustainability of water resources is a must. Urban, basin-oriented and global-scale integrated studies should be conducted to reduce the risk of water scarcity (Mekonnen and Hoekstra, 2016), and according to the results obtained in these studies, the development of strategies can only be achieved by managing water resources (Zeng et al., 2012). Management of water resources would only be possible through effective and planned water utilization, water quality controls (Bortolini et al., 2018), improvement of water distribution and operational efficiency, conservation of water resources, waste management and reduction of unnecessary irrigation (Yerli et al., 2019a) water productivity should be increased (Vanham and Mekonnen, 2021).
The concept of water footprint, which has been recently introduced to published literature is aimed at water resource management (Hoekstra, 2003;Ge et al., 2011;Galli et al., 2012;Morillo et al., 2015;Muratoğlu, 2019). The water footprint includes the measurement of the water volume required to produce certain goods or services, or to conduct the entire chain of production from raw material processing to consumption (Chico et al., 2013;Brindha, 2017). Water footprint is an approach that reveals the significance of water to the economy with the current consumption of society and allows environmental, economic and social analysis on various scales. Water consumption during the production of any product, any process, or any urban can be calculated with the water footprint approach.
Because water footprint is a new concept (Adetoro et al., 2021), the current unavailability of the required data (Ewaid et al., 2019), the high data requirement and the complexity of analyses, studies on water footprint are quite limited (Muratoglu, 2019). Generally, the studies have been conducted on countries and basins, and estimated data, general climatic inputs and major crop products have been used as the study data to conduct low resolution calculations (Chapagain and Hoekstra, 2004; Aldaya et al., 2010; Mekonnen and Hoekstra, e colheita das safras para um período que se beneficie da precipitação. Além disso, espera-se que este estudo contribua para novos estudos para o cálculo da escala provincial WF e terá efeitos positivos no planejamento agrícola, alocação de água e a sustentabilidade dos recursos hídricos. Palavras-chave: sustentabilidade agrícola, água azul, água verde, pegada hídrica, sustentabilidade hídrica. for many years) were used to calculate the water footprint for crop production. As a result, the study approach was more original and realistic. The study is expected to increase efficiency and sustainability of water resources with better planning in arid and semi arid provinces such as the Van province, and to pave the way for new studies.

Study area
Located in the Eastern Anatolia region of Turkey, Van province is one of the largest and most developed provinces ( Figure 1). The surface area of Van is 23.334 km 2 and 2.86% of the total land mass of Turkey. With this area, Van province constitutes a larger area than the Van Lake basin, which is one of Turkey's 26 basins. The altitude of the province, located at 38°29′39″ North latitude and 43°22′48″ East longitude, is approximately 1725 m. The total population of Van province was 942.771 in 2004 and 1.136.757 in 2019 according to data from the Turkish Statistical Institute (TSI, 2020), and has increased by approximately 21% over 15 years. The vast majority of the population is employed in agriculture and trade Industrial activities are almost non-existent in the province (Kanberoglu, 2016). Crop and livestock production are the main income sources. The Van province ranks first in Turkey in terms of small ruminant assets (TSI, 2020).
In Van province, where the terrestrial climate is dominant, winters are cold and summers are dry with little rainfall, although the number of frosty days is high throughout the year. However, climatic conditions are temperate due to the microclimate effect created by Lake Van, located in the center of the province. According to the long term meteorological reports from 1976 to 2019-from six different climate stations selected to represent the region-obtained from the Turkish State Meteorological Service, the annual average air temperature was 9.1°C, the average relative humidity was 54%, the total hours of sunshine was 94.8, the total evaporation was 1397.6 mm and the total precipitation was 387.5 mm (TSMS, 2020). According to the averages from the six climate stations, the monthly total and effective precipitation, and the average air temperature values are given in Figure 2. Accordingly, the total precipitation and average temperature were 27.2 mm and 20.5°C in the active crop production season (June-July-August). Inadequate precipitation and high temperatures led to the need for irrigation to achieve high yields.
Approximately 60% of the agricultural soil in the province consists of sandy clay loam, and the rest consists of sandy loam and sandy clay textures, i.e., medium-texture soils (Tufenkci et al., 2009). Approximately 24% of the agricultural area of the province is shallow and the rest has deep soil (Karaca et al., 2019). Van province has 3744 km 2 of agricultural land, and the amount of irrigated land is 1164 km 2 .

Method
The water footprint requires a high level of data input (Muratoglu, 2019). When calculating the water footprint for a certain area, the accuracy of the results is closely associated with the quality of the data entered (Godar et al., 2015). In the present study, current and realistic data over a long time period (2004 to 2019) were used instead of estimated data, and many national and international databases were used.
In this study, water footprint calculations were made according to the water footprint method developed by Hoekstra et al. (2011). The total water footprint of an area (WF total ) consists of the water footprint for crop production (WF crop ), the water footprint for livestock production (WF livestock ), and the domestic and industrial water footprint (WF domestic+industrial ) (Equation 1) (Muratoglu, 2019). The most important component of WF total is the WF crop due to intense water consumption (Lovarelli et al., 2016). WF livestock constitutes the amount of water consumed by the livestock in the region . WF domestic+industrial , which is usually the smallest component in the total water footprint (Klemeš, 2015), is the water volume consumed by individuals during domestic and industrial activities. The flow diagram used in the calculation of WF total for Van province is given in Figure 3.

Water footprint analysis of crop production
The amount of WF crop for a region depends on evapotranspiration (ET c ). ET c , which can be defined as the water requirement of the crops, changes depending on the region's climate conditions, soil properties and crop characteristics (Allen et al., 1998). The part of ET c for the crop that is met by effective precipitation, (P eff ), constitutes the green water amount for the crop (ET green ), and the remaining part is the irrigation requirement, i.e., the amount of blue water for the crop (ET blue ) (Hoekstra et al., 2011). When the effective precipitation (P eff ) is more than the crop evapotranspiration (P eff > ET c ), the irrigation requirement is eliminated, and by contrast, when P eff < ET c , the net irrigation requirement of the crops is determined as the evapotranspiration minus the effective precipitation (ET c -P eff ) (Equations 2 and 3) (Muratoglu, 2019). In order to determine the ET c values for all crops, the reference crop water consumption (ET o ) values calculated by the FAO-Penman Monteith approach for each ten days during growing periods. These were corrected with crop coefficients (kc) to calculate crop evapotranpirations (ET c = ET o × kc). ET c values calculated with this approach for the studied region were obtained directly from the last and current study (Turkey Irrigated Crop Water Consumption Guide) by the Turkish Ministry of Agriculture and Forestry (TAGEM, 2017). Planting-harvesting dates and vegetation periods of the crops were also obtained from the same study. In the daily ET o calculations in the guide mentioned, the data from the last 30 years, daily average, minimum and maximum temperatures, relative humidity, wind speed, hours of sunshine and sunshine intensity were considered. The scientific results of different field trials conducted in  Turkey and the study region, the data at the international level, graduate theses and other scientific publications and FAO-56 data were all used to obtain favorable kc values (Allen et al., 1998;TAGEM, 2017).
The benefit that crops can obtain from total precipitation can be defined as P eff (Aldaya & LIamas, 2008). Using the approach developed by the United States Department of Agriculture, Soil Conservation Service (USDA, 1993), monthly P eff values were calculated with CROPWAT 8.0 software, and all precipitation below 25 mm was accepted as P eff based on the same approach. In the water footprint studies in literature, P eff was included only in the the crop planting-harvest periods. However, at the time of planting, moisture in the soil (residual soil moisture from winter) (P winter ) could be used by the crops for germination and during the initial growing periods (Kanber, 2010). Especially in regions such as the Van province, where it is arid and the water resources are inadequate, irrigation is not possible at the time of planting, and P winter enables germination of the crops. P winter is the amount of moisture retained in soil from the winter months at the beginning of the growth period of the crops. P winter was calculated from the amount of moisture that the soil can hold. Retained rates of precipitation from winter months were obtained from the study conducted by Kanber (2010) from different regions of Turkey for every winter month. Every two weeks, a correlation was made using the assumption that there is a linear relationship between the amount of moisture retained in the soil over two months. P winter was determined by using this rate and finally the P winter value obtained was added to the total P eff in the vegetation period. In the present study, the water holding capacities were considered to relate to each crops's effective root depth. The water holding capacities of non-problematic soils were determined according to soil texture (Tufenkci et al., 2009), soil depths (Karaca et al., 2019) and the effective root depth of the crops obtained from the FAO database (FAO, 2018).
WF crop calculations were done for the twenty-one most cultivated crops (barley, green and dry beans, chickpeas, green lentils, potatoes, silage corn, sugar beet, alfalfa, sainfoin, cabbage, grapes, apples, apricots, strawberrys, melon, watermelon, tomatoes, eggplants and onions). Production volumes and the amount of production areas for these crops were obtained from the Turkish Statistical Institute (TSI, 2020). Average yield values (ton/da) were obtained by dividing the production volume of each crop by the amount of production area.
Finally, the ET green and ET blue values and the average yield values (Y) for each crop in the province were used to calculate the WF crop (Equations 4 and 5) (Muratoglu, 2019). The water consumption volume during crop production was calculated as the per unit volume (m 3 ton -1 ), and this process was applied to all crops used, while the blue and green water footprints (WF blue and WF green ) were found separately. By collecting WF blue and WF green for each crop, the water footprint of a crop was determined (Equation 6), and by collecting the water footprints of all crops, the total WF crop for Van province was obtained (Equation 7) (Muratoglu, 2019 .

Water footprint analysis of livestock production
The vast majority of water consumed in the WF livestock consists of water used in the production of feed crops, while the rest is the drinking water consumed by the livestock and the indirect water consumption. Water consumed for the production of feed crops is included in WF crop . However, there are also studies in literature where the production of feed crops is included in WF livestock (Muratoglu, 2019(Muratoglu, , 2020. In the current study, blue water footprints for animals were evaluated instead of animal products such as meat, milk and eggs. The calculation of WF livestock was conducted with global average data because regional data were not available. Average data on a global scale for beef cattle, dairy cattle, goat, sheep, poultry (including broiler and layer chickens, turkey, goose and duck) and horse (including donkey and camel) were obtained from published literature . Data on the number of livestock obtained from TSI (2020) was multiplied by the annual water consumption (m 3 year -1 ) of a animal to calculate WF livestock . Gray water is not included in WF livestock because there is no available data on the gray water footprint.

Domestic and industrial water footprint analysis
The domestic water footprint of an area includes blue water consumption by the domestic activities of individuals living in that area (Chapagain and Hoekstra, 2004). The industrial water footprint is the water consumed in industrial activities. Gray water is an important concept in regions where industrial activities are intense. However, in regions with low industrial capacity, gray water footprint does not need to be included in the total water footprint calculations (Chapagain and Hoekstra, 2004;Hoekstra et al., 2011;Muratoglu, 2019Muratoglu, , 2020. Therefore, since the industrial activities in Van province are low (Kanberoglu, 2016), the gray water footprint was not evaluated, and the industrial and domestic water footprint were analyzed together. For this purpose, the data obtained from TSI (2020) was utilized, and WF domestic+industrial was found by multiplying the amount of water consumption per capita in the districts by the number of people residing in the districts from 2004 to 2019.

Blue water scarcity
The importance of blue water for regions with insufficient water resources or for arid regions is clear. Therefore, many different methods are used and developed to identify and understand the water resources, especially in these regions. However, flow requirements and real run-off amounts are usually ignored in the studies related to the sustainability of fresh water (Mekonnen & Hoekstra, 2011b), but the approach of  provides a different perspective. Because, according to other water analyses, blue water is used instead of using water withdrawal, environmental flow requirements and employing natural undepleted run-off (Muratoglu, 2019). The sustainability of blue water in the study area was determined using this approach. Accordingly, blue water scarcity (WS blue ) was obtained by assessing WF blue in the region in proportion to blue water availability (WA blue ) (Equation 8). WA blue was defined as the natural runoff amount minus the environmental flow requirements. The natural runoff value was calculated by summing the real runoff and WF blue in the region. The real runoff value of the Van province was obtained from the study of Mekonnen and Hoekstra (2011b). The environmental flow requirement was taken to be 80% of natural runoff, as Hoekstra et al. (2011) proposed. The result obtained was classified as low (WS blue < 100%), moderate (100 < WS blue < 150%), significant (150 < WS blue < 200%) and severe (WS blue > 200%) water scarcities .

Results and Discussion
In this study, in which the WF total for the years 2004 to 2019 was calculated in the Van province, the average annual water footprint was found to be 8.73 billion m 3 . Of this, 87.6% (7.66 billion m 3 ), 4.9% (0.43 billion m 3 ) and 7.5% (0.64 billion m 3 ) of the WF total consisted of WF crop , WF livestock and WF domestic+industrial , respectively. 62.5% and 37.5% of this WF crop corresponded to WF blue and WF green (Figure 4). On a global scale the WF total calculations by Mekonnen and Hoekstra (2011c) showed that production was found to be based on WF green . However, the opposite situation was observed for Van province. As expected, WF crop created the largest water volume in the province, and the vast majority of WF crop was based on WF blue . It was clear that this situation was due to the increase in the need for irrigation water because of low annual precipitation in the province as well as the amount of precipitation falling during the crop development period, which was also very low ( Figure 2).

Per capita water footprint
The per capita water footprint values are the best approach (Muratoglu, 2020) for a comparison with the water footprints of different regions (Cheng et al., 2019). The per capita water footprint in developed regions are in the range of 1250 to 2850 m 3 . The average per capita water footprint of the Van province between 2004 and 2019 was found to be 1385 m 3 year -1 capita -1 and is lower than the world average, which was 1385 m 3 year -1 capita -1 , and the average in Turkey, which was 1642 m 3 year -1 capita -1 (Pegram et al., 2014) between 2000 and 2011. When the per capita water footprint of Van and other countries we are examined, it was seen that Van has a very low water footprint like other arid areas. This is related to the volume of consumption, consumption pattern (Liu et al., 2008), the impact of the climate and agricultural practices (Chapagain et al., 2006). It is also thought that high migration rate from different regions to Van creates a rapidly increasing population that has a high impact on it.
When the per capita water footprints in the study were analyzed, it was seen that the lowest value was in 2009 and the highest value was in 2012 ( Figure 5). This situation is closely related to the drought experienced in 2009 and the high precipitation experienced in 2012 (TSMS, 2020). Per capita water footprints in the province were not consistent after 2009 where the sudden decline of some products with high water consumption, such as chickpeas, lentils and grapes in the province (TSI, 2020), had an effect on the per capita water footprints. In addition, the increase in out-migration since 2010 and the decrease in the population (TSI, 2020) are considered to be important in the increase in the water footprint as of 2010.

Water footprint of crop production
In the agricultural sector, which has higher water consumption than other sectors (Mancosu et al., 2015; Karci and Ucar 2019), the majority of freshwater is used for irrigation of crops (Galan-Martin et al., 2017). Especially in arid and semi-arid regions, more irrigation water is needed, i.e., blue water. In Van province, 87.6% of the WF total consisted of WF crop , and this was supported by similar studies (Chapagain and Hoekstra, 2004;Aldaya and LIamas, 2008;Muratoglu, 2019Muratoglu, , 2020. The WF blue and WF green components of WF crop were found to be 4.81 billion m 3 and 2.85 billion m 3 , respectively. While 75% of the WF crop consisted of WF green on average in Turkey (Mekonnen and Hoekstra, 2011a), 78% of the WF crop consisted of WF green on average, worldwide (Mekonnen and Hoekstra, 2011c). The WF green of WF crop was found to be 37.5% in the Van province and this led to the conclusion that the demand for blue water for crop production was high. The demand for freshwater increases as a result of the high evapotranspiration combined with insufficient precipitation during the crop producing season (Siebert and Döll, 2010). High WF blue values restrict agricultural production, especially in regions with poor water resources (Mali et al., 2018). WF blue is a risk factor that predicts water scarcity in the future. For this reason, having WF green higher than WF blue is very important for the sustainability of water resources (Novoa et al., 2019).
The average virtual water content per crop in Van province was 478.3 m 3 ton -1 , and is given in Figure 6 as WF blue and WF green for each crop. Accordingly, the WF blue of all crops except for onions is more than WF green . This situation is a result of insufficient precipitation in the province during the vegetation period. The reason why the onion WF green is more than the WF blue is related to the crop's abundant use of precipitation as a result of the early planting date (15 May) and short vegetation period (75 days). The crop with the highest water footprint in the province is chickpeas, which have been found to have a high virtual water content in many studies (Chapagain and Hoekstra, 2004;Gupta, 2008;Mohammadi-Kanigolzar et al., 2014), while the lowest is cabbage. The WF green and WF blue of chickpeas, which have the highest amount of green water, are almost equal to each other. Lentils, beans (green), chickpeas, apples, grapes and walnuts are crops with the highest WF blue values in the province, respectively. As Fu et al. (2018) stated, the crops that have high yields according to other regions, also have a high virtual water content. In addition, Pegram et al. (2014) stated that the water footprint of fruits is generally found to be higher than vegetables.  The WF green and WF blue values of crop products, which are important in terms of production capacity in the Van province are compared with the Upper Tigris River Basin, the worldwide average and the Turkish average, and are graphically illustrated in Figure 7. The WF green and WF blue values for apricots, tomatoes and cabbage in Van province are similar to the worldwide average. It was seen that the Van and worldwide average WF green values for chickpeas, potatoes, sugar beet, watermelons, eggplants and onions were similar. In particular the WF green values for beans (green), lentils and barley in Van are well below the worldwide average, and the WF blue value for chickpeas, apples, beans (green) and lentils were found to be well above the worldwide average. In general, the WF green and WF blue values for the crops in Van province are similar to the worldwide average, but the total water footprints of some crops (chickpeas, apples and grapes) are found to be much higher. This situation may be related to crops with high water footprints that have low yields in the region (Mekonnen and Hoekstra, 2010). In the comparison of the WF green for Turkey and for Van province, it was seen that the WF green for barley is higher in Turkey, and the WF blue is similar, while WF green and WF blue for grapes and beans (green) are higher in Van province. The WF green and WF blue values for silage corn (except WF green ), potatoes, sugar beet and tomatoes in Van province are relatively similar to the Turkish average. The main reason for these similarities and differences arises from the differences due to soil distribution and water resources (Sun et al., 2013) and especially to climate (Mekonnen and Hoekstra, 2011a). In addition, the WF crop values for the Upper Tigris River Basin, located to the south of the study area, is compared to the Van province. Accordingly, WF green for crop production in Van province is found to be higher for almost all crops. The reason for this is due to the inclusion of residual soil moisture from winter when calculating the WF crop for Van, i.e., the difference in method. WF green and WF blue values for onions and sainfoin showed similarity between the two regions, and the WF blue value for barley, clover, apricots, melons, watermelons, tomatoes and eggplants were lower in Van province. This situation can be explained by differences in the irrigation regime of the regions (Boelens and Vos, 2012), the differences in the crop varieties and variability of water requirements depending on the climatic differences of the regions. To summarize the WF crop for the Van province, it is possible to say that there are some changes in evaluation compared to other different studies, but in general, these changes diverge according to the regional conditions (Fader et al., 2010) and are in good agreement with other studies.

Water footprint for livestock production
The volume of water consumed directly or indirectly by animals during the production of meat, milk, eggs and other products can be defined as WFl ivestock (Bosire et al., 2017;Ibidhi et al., 2017). In Van province, WF livestock has been found to be 0.43 billion m 3 and constitutes 4.9% of the WF total . Cattle, which have the highest water footprint in the province (Figure 8), increased its water footprint. Because the per calorie water footprint of beef corresponds to nearly 20 times that of cereals (Hoekstra, 2012), high meat consumption in any region may increase the water footprint of that region (Chapagain and Hoekstra, 2004). Similar results were obtained with Hoekstra (2012), who stated that more than half of the WF livestock value is made up of cattle on a global scale. The annual water consumption of dairy cattle is almost 30 times higher than that of sheep . Dairy cattle come to the fore in terms of water footprint because of this, even though the number of sheep is far greater than dairy cattle in Van province. The reason for the low water consumption of goats and poultry is associated with their lower annual water consumption (32 and 33 m 3 year -1 capita -1 ) compared to other livestock (Hoekstra, 2012).
When the WF livestock of Van province was analyzed by year, it was seen that the lowest value was in 2012 and the highest value was in 2017 (Figure 8). The reason for the increase in 2017 is that 2017 featured the highest number of animals compared to all other years (TSI, 2020). The reason for the decrease of the WF livestock value in 2012 was thought to be related to the major earthquake in September 2011, and accordingly to a decrease in the number of animals. Feed crops consumed by livestock animals correspond to approximately one third of the total WF livestock value (Mekonnen and Hoekstra, 2011c). However, feed crops were not included in the WF livestock value but were evaluated as part of WF crop . If the production of feed crops for Van province was included in WF livestock , 75.3%, 18.6% and 6.1% of the WF total value are WF crop , WF livestock and WF domestic + industrial , respectively ( Figure 9). This situation reveals that feed crops consume high amounts of water in Van province. Comparing the conditions in which feed crops were included and not included in WF crop , blue and green water distributions were almost unaffected (WF blue : 65.2% and WF green : 34.8%). Looking at WF livestock under conditions in which feed crops were included in WF livestock , it was seen that 64.4% and 35.6% of WF livestock consisted of WF blue and WF green , respectively (Figure 9). Mekonnen and Hoekstra (2012) stated that 87% of the global WF livestock value corresponded to WF green , but in the Van province, livestock production dependent on blue water was carried out like crop production. When the values of Van province were compared with the Upper Tigris River Basin (Muratoglu, 2019) and Diyarbakır province (Muratoglu, 2020), it was seen that the crop, livestock and domestic-industrial water distributions were similar but that the WF blue value within the WF livestock value was higher in the Van province. This may be related to the use of high levels of blue water to increase the yield in the feed crops due to the high number of animals-in addition to being an arid region of Van province. When calculating WF livestock for a region, considering that feed crops can be procured from other regions and also that ready-made feed may be used, the inclusion of feed crops in WF crop will reveal more precise and realistic results.

3.4.Domestic and industrial water footprint
Water used in domestic needs and industrial activities make up WF domestic+industrial . The average WF domestic+industrial of Van province between 2004 and 2019 was found to be 0.64 billion m 3 , and this value corresponded to 7.5% of the WF total . Figure 10 shows that domestic and industrial water use has been increasing steadily with increasing population, and this increase only decreased in 2012. It is thought that the reason for the decrease in 2012 was because of a population decrease due to political reasons. Because the Van province is one of the least industrialized regions in Turkey (Kanberoglu, 2016), it is estimated that domestic water consumption is higher than industrial water consumption. In addition, considering that Turkey's total gray water footprint is 17% (Pegram et al. 2014), the gray water footprint value for Van province, which has limited industrial activities, is likely to be a very low amount. WF domestic+industrial values obtained for Van province are similar to regions with low industrial capacity (Aldaya & LIamas, 2008;Zeng et al., 2012;Muratoglu, 2019). Van Oel et al. (2009), in their study, stated that WF domestic+industrial corresponded to 33% of the WF total . It was expected that the values for WF domestic+industrial in regions with high industrial activities and dense population, would be high.

Blue water scarcity
Although WS blue analysis can be applied on monthly and annual scales, its monthly determination provides a more detailed view as it will reveal seasonal variations    (Muratoglu, 2019). However, since the main output of this study was the annual WF blue average, WS blue analysis was evaluated annually. Accordingly, the average annual WS blue of Van province between 2004 and 2019 was found to be 257%, and it is a region with severe water scarcity in the classification stated by . This situation shows that 257% of blue water was used in Van province every year. If this situation continued, it is inevitable that it will experience serious water shortages in the near future.
It has been pointed out from WS blue analysis of many basins that four billion people face water scarcity (Mekonnen and Hoekstra 2016). Mekonnen and Hoekstra (2011b) stated that in a study involving WS blue analaysis on 223 basins, where 201 basins consumed more than twice the available blue water for at least one month a year. WS blue exceeds 200% in the Euphrates-Tigris basin, which is located to the south and west of the Van province, especially in the summer when crops are produced (Degefu et al., 2018). In the Van, WS blue was estimated to be higher in the summer months compared to locations like the Euphrates-Tigris basin, because the crop production season has inadequate precipitation and high temperatures ( Figure 2). The WS blue value for that part of the Tigris-Euphrates basin in Turkey was reported to have an annual average of 87% (Muratoglu, 2019). The reason for the WS blue difference between the two province is due to the differences in the actual runoff, in the amount of WF blue and in consumption habits and especially the climate. In addition, WS blue can change according to the feeding conditions of the basins from different streams and the conditions of the region they serve.

Conclusions
In this study, the water footprints in Van province, which is located in the most easterly part of Turkey and is Turkey's fifth largest province in terms of surface area, were calculated. Crop production, livestock production and domestic-industrial water footprints have been analyzed over a wide time period (2004 to 2019) using current data, and the blue and green water footprints have been discussed in detail. The average water footprint of Van province was found to be 8.73 billion m 3 . Crop production is the main water-consuming sector in Van province being responsible for 87.6% of the total water footprint. More than half of the water footprint for crop production consists of blue water, i.e., underground and aboveground freshwater sources. Livestock production and domestic-industrial water footprints accounted for 4.9% and 7.5% of the total water footprint, respectively.
The crops that require intense water consumption are chickpeas, lentils, apples, beans (green), walnuts, grapes and apricots, respectively ( Figure 6). These crops accounted for approximately 84% of the water footprint for crop production in Van. Therefore, the production of crops with low water consumption, such as onions, sugar beet, corn, potatoes, cabbage, watermelons, melons and eggplants, can be recommended for this province. When similar studies in literature were examined, it was seen that crop production generally consists of a green water footprint. However, it was concluded that the blue water footprints of crop products is higher than the green water footprints in Van province. Although this situation is caused by drought in the region, it may cause concerns in terms of the sustainability of the region's water resources in future years. Therefore, comparing the use of blue water with the economic value of the crops, and explaining the efficiency obtained from each unit of blue water, indicates the precautions to be taken in line with these results. In order to increase the efficiency of green water, changing the planting pattern of the crops or matching the planting-harvest dates of crops to the period that will benefit from precipitation is important for the protection of water resources. In arid and semi-arid regions, crop production using green water and the supply of other products from different regions may illustrate a different perspective on water footprint reduction. In addition, reducing evapotranspiration, rainwater harvesting and the use of modern irrigation methods can be considered as alternative approaches to reduce water footprint. In addition, biotechnological research is needed to develop crops that provide higher yields with less water consumption. In the Van, when feed crops are included in the water used by livestock, the water footprint of livestock production increased by 13.7% (Figure 9). This ratio revealed that feed crops have a high water consumption in this province. Therefore, it is important to use ready-made feed or use different regions to supply feed crops to help secure the future of water resources in arid and semi-arid regions. The study region, where migration is constantly experienced for a variety of reasons, exhibits a rapid population growth and low consumption, and the amount of freshwater per capita is rapidly decreasing-especially in the last four years (Figure 1)-while the domestic and industrial water footprints have increased over the years ( Figure 10). For these reasons, the authorities need to take measures to reduce migration and the resulting population increase, and develop strategies to distribute the population evenly over its area.
The blue water scarcity value for Van province was found to be 257%, which means 257% of available blue water is consumed every year. This value is an indication that Van province will experience a serious blue water crisis in the near future. For this, incentives should be provided for the consumption of food, goods and services with low water footprints, and steps should be taken in this direction by the authorities. The reduction of the water footprint for production in Van province is important for the sustainability of freshwater. Crops with high water footprints should be exported and detailed water footprint calculations should be renewed at certain periods. In addition, high resolution water footprint maps can be created for especially arid and semi-arid regions such as the Van province. Of course, data records should be established first for all these calculations. After these data are obtained, high temporal resolution studies on the blue water scarcity would provide better visibility of the sustainability of freshwater resources in the study region. Since there is a blue water crisis in the study region, reducing the use of blue water in the agricultural sector, carrying out studies on green water use and raising awareness of this is of great CAKMAK, B., KENDIRLI, B. and UCAR, Y., 2007.  importance for the sustainability of water resources. It is recommended to cultivate products with high economic value requiring a low amount of blue water instead of products with low economic value and high blue water usage such as forage crops. For future research, the water footprint for consumption and virtual water transfer in the Van is suggested for investigation to reveal the water savings or deficits.
With this urban-scale study, it was concluded that the water footprint is a very useful approach to quantify water resources. Therefore, increasing studies to calculate water footprints on an urban scale will make a great contribution to developing policies and strategies. Finally, urban and basin-oriented water footprints based on production and consumption should be investigated by using up to date and realistic data instead of estimated data, and so comprehensive information can be obtained by carrying out studies in this area.