Bridging the Data Gap in the Water Scarcity Footprint by Using Crop-Specific AWARE Factors

The assessment of the water scarcity footprint of products emerged as an important step in supporting water management strategies. Among others, the AWARE methodology was published as a consensus-based indicator to perform such an assessment at a watershed level and monthly scale. The need to adopt such a detailed resolution, however, collides with the availability of data, so that general year and country-wide factors are commonly used. The objective of this study is to develop and verify the applicability of 26 crop-specific water scarcity characterization factors to help assess the water scarcity footprint when data and information availability is limited. To do so, a weighted average consumption approach was adopted, starting from local AWARE characterization factors and local crop-specific water consumption. The resulting factors, ranging from 0.19 m3/ton eq for “other perennial crop” in Brunei to 9997 m3/ton eq for “other annual crop” in Mauritania, illustrate the large variability of potential water scarcity impacts. Factors were applied to the water consumption of selected crops to assess their water scarcity footprint. The results of the study confirmed that the use of crop-specific factors is recommended as they are a better proxy of water scarcity in a region when compared to their national generic counterparts.


Introduction
In the framework of the Sustainable Development Goals [1], with the global population expected to be 9 billion people by 2050, agriculture is considered to play a fundamental role in human well-being and economic growth [2]. With a production that tripled between 1965 and 2015, the sector is a key contributor to global food security [2], renewable energy production (14% of global energy production) [3], and a variety of new emerging markets, such as bioplastics [4]. Its relevant position in the global economy, however, is paired with growing concerns about the environmental impacts related to agricultural production systems. With reference to climate change, for example, the sector is responsible for 11% of global direct greenhouse gas emissions, mainly resulting from land use practices and the use of fertilizers and pesticides. Moreover, with specific reference to water, agricultural products are responsible for the consumption of 1180 km 3 yr −1 [5], representing about 70% of the total water consumption [6], thus granting this sector a leading role in order to guarantee the sustainable use of this resource. In this context, the need for tools and methods to adequately quantify and address undesired environmental impacts has clearly emerged [7].
Among different tools and methods, Life Cycle Assessment (LCA) has been extensively applied for the quantification of potential environmental impacts of products and processes and to support decision making processes for the correct management and reduction of related negative impacts [8].
The strength of LCA resides in its capability to quantify a wide range of potential environmental impacts from a life cycle perspective, from raw material extraction to the end of life management of products, avoiding the issue of an environmental burden shift [9]. To do so, in LCA studies, the environmental inputs and outputs of each unit process that take place in the life cycle of a product are translated into potential environmental impacts using suitable environmental characterization factors (CF) [8]. With specific reference to water, also due to the introduction of ISO 14046 standard [10], the metric of the water footprint in LCA has emerged to quantify potential environmental impacts related to water. The assessment considers the potential effects of the release of a pollutant that can affect the quality of water resources, as well as issues related to the scarcity of water resources (also referred as the water scarcity footprint). The assessment of this last aspect has received significant attention from the scientific community in order to define a consensus method [11]. This is the case of the AWARE [12] methodology that was published as the consensus-based water consumption impact indicator to assess the LCA-based water scarcity footprint (WSF), as the outcome of the UNEP-SETAC Life Cycle Initiative flagship project Global Guidance for Life Cycle Impact Assessment Indicators and Methods (GLAM). The authors provided characterization factors for over 11,000 watersheds in the world (with the 34 largest watersheds subdivided into sub-watersheds). The method is based on the causality impact pathway approach, focusing on a generic scarcity-based proxy midpoint indicator (labeled proxy as it does not directly lead to endpoint impacts on human health or ecosystem quality specifically) applicable to water consumption (surface and groundwater) [12]. The development of the method confirmed the importance of adopting a resolved geographical and temporal scale when assessing the water scarcity issue. Whenever possible, performing the quantification starting from specific water use data at a watershed level and monthly scale is recommended. In addition, other studies have confirmed that the geographical scale significantly influences the accuracy of the results [13,14]. The need to adopt such a detailed geographical resolution, however, collides with the availability and completeness of water consumption data and information [12,15,16]. In water footprint and LCA studies, it is in fact usually necessary to use secondary and background data to determine the potential environmental impacts that take place during the life cycle of products [7].
To overcome this limit, the authors of AWARE also provided general factors for assessing the water scarcity footprint of water consumption at the country level, but did not recommend using them. Instead, for cases in which watershed-level factors cannot be used, they provided two sector-specific factors: agri or non-agri, based on whether the activity assessed is agricultural in nature (i.e., irrigation) or industrial/domestic. These factors, calculated using a consumption-weighted average of the watershed-scale characterization factor, better represent the temporal and geographical distribution of water consumption in the country for an agri or non-agri activity, and therefore provide a more representative weighted average than if weighted with all water consumption, as is done for the unknown factor [12].
On the same line, Northey et al. [17] presented specific characterization factors to address the water scarcity of 25 commodity-based products of mining industries, starting from AWARE factors. In this case, the characterization factors were developed by adopting a production-weighted-based approach. The study also confirmed that sector-specific characterization factors represent a better proxy of environmental impacts when compared to the more generic national ones; in some cases, results of the quantification can in fact be overestimated by up to 60%. The authors concluded that in the case where no watershed level data are available, commodity factors are therefore a better proxy than the national average factors.
Considering the interest of the agricultural sector in reducing their impacts related to water [18], the development of crop-specific AWARE characterization factors could improve the accuracy of assessment in the agricultural sector and facilitate the application of the water footprint and LCAs.
The aim of this research is to improve the accuracy of crop water scarcity footprint quantification when no watershed level data are available by

•
Calculating and providing crop-specific AWARE characterization factors for 26 crops' typology at a national level for 224 countries, and • Demonstrating their applicability in quantification of the water scarcity footprint of selected crops.

Materials and Methods
In this section, the AWARE method is summarized, as well as the mathematical approach for the calculation of the crop-weighted characterization factors (Section 2.1). Details on the data used (Section 2.2) and the ArcGIS calculations (Section 2.3) are provided below. Lastly, the calculation of the water scarcity footprint of different crops is shown (Section 2.4).

AWARE Method and Consumption-Weighted Average Approach
The AWARE method [12] is based on quantification of the relative available water remaining per area once the demand of humans and aquatic ecosystems has been met. This method builds on the assumption that the potential to deprive another user of water is directly proportional to the amount of water consumed and inversely proportional to the available water remaining per unit of surface and time in a region (watershed).
The method is therefore defined by starting from the difference between freshwater availability and demand (Availability Minus Demand-AMD) in a specific i-region (Equation (1)).
Availability was determined from WaterGAP2.2 [19], which models runoff based on the precipitation and evapotranspiration of >11,000 watershed in the world, for the period 1960-2010. Demand includes both human and ecosystem demands. The former (Human Water Requirements-HWR) was also determined by starting from the WaterGAP model and refers to data on the consumption of freshwater for domestic, industrial, irrigation, livestock, and energy uses (based on the year 2010). The latter is determined by Environmental Water Requirements (EWR) developed as a fraction of the available flow to maintain freshwater ecosystems in "fair" conditions, with respect to pristine flow (i.e., flow with no human influence), and where "fair" refers to an ecological state of the river between poor and good conditions [20].
The characterization factor relates the AMD i to the AMD worldavg , the consumption-weighted average of all AMD i in the world, as per Equation (2). This is valid for regions where the value of AMD i is not smaller than 1% of AMD worldavg and not larger than 1000% (10 times) of AMD worldavg . Outside of this range, the factor is set to the minimum (0.1) or maximum (100) value.
CF i values were calculated at a watershed level and by adopting a monthly temporal scale. They were than aggregated at a country and year level, as well as at a sector level (agri and non-agri per country), based on a consumption-weighted average approach.
In this paper, the same approach is adopted to develop crop-specific CF at a country level. To do so, Equation (3) is used: where CF j,z is the characterization factor of the j-th crop of the z-th country, CF i is the annual (agri) characterization factor of the i-th watershed, WC j,i (water consumption) is the water consumption required for the growth of the j-th crop within the i-th watershed, and WC j,z is the total water consumption required to grow the j-th crop in the z-th country.

Data
Data on crop-specific water consumption at a local scale covering the entire globe was taken from the Global Crop Water Model (GCWM) [21]. This model is a simulation of the yield and consumption (mm/yr) of green and blue water for 26 crops or different crop groups, for the years from 1998 to 2002, using climate data, the WaterGAP model. cultivated areas, cropping calendars, and evapotranspiration specific to each soil and crop (ETc). The crops include wheat, maize, rice, soybean, cotton, canola and rapeseed, barley, rye, managed meadows and pastures, sunflower, legumes, cassava, citrus, cocoa, coffee, date palm, grape and vine, peanuts, millet, coconut palm oil, potato, sorghum, sugar beet, sugar cane, and "other annuals", as well as "other perennials".
Blue water consumption data (i.e., water consumption for irrigation) refer to the irrigation water requirements for optimal crop growth and for the years from 1998 to 2002. Considering that climatic conditions vary over time, the blue water consumption data considered in this study are the averages of these years. Validation of the data was performed by comparing the total water consumption for agriculture obtained from WaterGap (and used in AWARE modeling), with the sum of the water consumption per country for the 26 crops above from GCWM. In this case, validation involves a comparison with a model output already published and validated [19,21]. In this paper, we are using these output data in AWARE [12].
Since the water consumption data per crop were only available at the yearly time scale, the annual AWARE agri CF was used, at the watershed scale, in Equation (3) above.

GIS Calculations
The calculation of crop-specific AWARE factors was done in ArcMap by cross-referencing blue water consumption data (five-degree raster scale, about 10 km at the equator) from GCWM and existing annual AWARE agri factors per watershed (see Figure 1), to create country and crop-specific AWARE factors. Data on crop-specific water consumption at a local scale covering the entire globe was taken from the Global Crop Water Model (GCWM) [21]. This model is a simulation of the yield and consumption (mm/yr) of green and blue water for 26 crops or different crop groups, for the years from 1998 to 2002, using climate data, the WaterGAP model. cultivated areas, cropping calendars, and evapotranspiration specific to each soil and crop (ETc). The crops include wheat, maize, rice, soybean, cotton, canola and rapeseed, barley, rye, managed meadows and pastures, sunflower, legumes, cassava, citrus, cocoa, coffee, date palm, grape and vine, peanuts, millet, coconut palm oil, potato, sorghum, sugar beet, sugar cane, and "other annuals", as well as "other perennials".
Blue water consumption data (i.e., water consumption for irrigation) refer to the irrigation water requirements for optimal crop growth and for the years from 1998 to 2002. Considering that climatic conditions vary over time, the blue water consumption data considered in this study are the averages of these years. Validation of the data was performed by comparing the total water consumption for agriculture obtained from WaterGap (and used in AWARE modeling), with the sum of the water consumption per country for the 26 crops above from GCWM. In this case, validation involves a comparison with a model output already published and validated [19,21]. In this paper, we are using these output data in AWARE [12].
Since the water consumption data per crop were only available at the yearly time scale, the annual AWARE agri CF was used, at the watershed scale, in Equation (3) above.

GIS Calculations
The calculation of crop-specific AWARE factors was done in ArcMap by cross-referencing blue water consumption data (five-degree raster scale, about 10 km at the equator) from GCWM and existing annual AWARE agri factors per watershed (see Figure 1), to create country and crop-specific AWARE factors. In order to allocate the water consumption of the raster to the different entities created within that raster, the simplest approach (option 1) would have been to assume an equal distribution of water consumption across the raster and allocate it proportionally to the area occupied by the entity. An entity refers to the region created from the overlay of different scales of data. However, upon analysis of the results with this approach, several "errors" were found. An error consists of the In order to allocate the water consumption of the raster to the different entities created within that raster, the simplest approach (option 1) would have been to assume an equal distribution of water consumption across the raster and allocate it proportionally to the area occupied by the entity. An entity refers to the region created from the overlay of different scales of data. However, upon analysis of the results with this approach, several "errors" were found. An error consists of the attribution of a fraction of a consumption raster to an entity (country or watershed) where no other water consumption for that crop is found. For this reason, another approach (option 2) was chosen: if the raster was divided in two, the consumption was allocated to the largest of the two entities. If it was divided into three or more parts and none of the parts were larger than 50% of the area, then the consumption of this raster was not considered in the calculation. Hence, because of the small surface areas neglected, this had little influence for the relevant crop-producing countries and allowed us to fix the errors created by option 1 and avoid wrongfully attributing crop irrigation to an entity where that crop is not produced or irrigated. For example, Andorra, which has no consumption raster entirely or even predominantly in its territory, was allocated an AWARE of 90 for wheat with option 1. This AWARE was corrected (and no wheat AWARE was attributed to Andorra) with option 2. On the other hand, for Albania, which does produce wheat with irrigation, the value varies insignificantly between the two options (49.39 with option 1 versus 49.02 with option 2). Subsequent manual correction was performed for country-overlapping rasters. The use of option 2 above did not avoid the creation of artefacts where the largest portion of a raster falls in a country which presents no other raster of irrigation of that crop, and that does not produce it (with irrigation) at all. Hence the correction was made based on the criteria that if country A shows no complete irrigation raster in its territory for a specific crop, and only fractions of rasters from a neighboring country B, which itself has at least one complete raster, then no water consumption for this specific crop is allocated to country A.

Water Scarcity Footprint of Selected Crop Calculations
The resulting factors were then used to calculate the water scarcity footprint of four specific crops for illustration purposes (wheat, maize, rice, and soybean) in different countries. This was done by multiplying the crop-specific AWARE CF at the country scale calculated above by the total water consumed for that crop's irrigation per ton, as obtained from the Water Footprint Network [22].

Data
While the AWARE model uses the total water consumption for irrigation (from WaterGap), the data obtained from GCWM provide irrigation data for each of the 26 crops. As can be expected, discrepancies can be found between the total irrigation volume per country versus the sum of irrigation volumes for the 26 crops. While, for most countries, WaterGap irrigation data are higher, for 39 countries (out of 187), the WaterGap value is lower. These differences are the result of the differences in the underlying data and hypotheses of the two models [19,21]. The difference between the two data sets is, on average, a factor of 6 (ratio of irrigation Watergap Data/GCWM data), and the mean difference coefficient [13] is 0.61. This is calculated as the arithmetic average of the difference divided by the average of each data pair. Therefore, it varies between 0 and 2. These differences are discussed in Section 4.
The total water consumption for irrigation per country is shown in Figure 2 for the 58 countries with the largest consumption volumes, for both WaterGap and GCWM data.

Crop-Specific Characterization Factors
The crop-specific AWARE factors for selected countries are shown in Figure 3 and reflect the spatial variation of crop production (CF values of other countries are available in the supporting information). The units of the CF are dimensionless, and expressed in m 3 world eq/m 3 according to the AWARE method [12]. The resulting factors illustrate the large variability of potential water scarcity impacts, depending on the crop assessed within a country, due to the associated location of the activity. The average difference across the 224 countries between the highest and lowest cropspecific AWARE factor is 12.6, with the highest values being 95.9, 77.0, 72.6, 65.9, and 58.7 for Italy, Australia, South Sudan, Argentina, and the USA, respectively (Figure 4). This means that in these regions, the information on the type of crop you are growing can influence your water scarcity results even more so than the country your activity is taking place in.
The characterization factors for rice, maize, soy, and wheat are shown on a map in Figure 5, and all CF values for all crops are available in SI. For example, as can be seen in Figure 5, in Argentina, the potential impacts in m 3 world-eq per m 3 of irrigation increase from rice (0.45), soy (3.9), and maize (9.5) to wheat (17.4). These differences represent the spatial variability of the location of the different cultures and indicate that wheat is grown in regions of higher scarcity than maize and so on.

Crop-Specific Characterization Factors
The crop-specific AWARE factors for selected countries are shown in Figure 3 and reflect the spatial variation of crop production (CF values of other countries are available in the supporting information). The units of the CF are dimensionless, and expressed in m 3 world eq/m 3 according to the AWARE method [12]. The resulting factors illustrate the large variability of potential water scarcity impacts, depending on the crop assessed within a country, due to the associated location of the activity. The average difference across the 224 countries between the highest and lowest crop-specific AWARE factor is 12.6, with the highest values being 95.9, 77.0, 72.6, 65.9, and 58.7 for Italy, Australia, South Sudan, Argentina, and the USA, respectively (Figure 4). This means that in these regions, the information on the type of crop you are growing can influence your water scarcity results even more so than the country your activity is taking place in.
The characterization factors for rice, maize, soy, and wheat are shown on a map in Figure 5, and all CF values for all crops are available in Supplementary Materials SI. For example, as can be seen in Figure 5, in Argentina, the potential impacts in m 3 world-eq per m 3 of irrigation increase from rice (0.45), soy (3.9), and maize (9.5) to wheat (17.4). These differences represent the spatial variability of the location of the different cultures and indicate that wheat is grown in regions of higher scarcity than maize and so on.

Crop-Specific Characterization Factors
The crop-specific AWARE factors for selected countries are shown in Figure 3 and reflect the spatial variation of crop production (CF values of other countries are available in the supporting information). The units of the CF are dimensionless, and expressed in m 3 world eq/m 3 according to the AWARE method [12]. The resulting factors illustrate the large variability of potential water scarcity impacts, depending on the crop assessed within a country, due to the associated location of the activity. The average difference across the 224 countries between the highest and lowest cropspecific AWARE factor is 12.6, with the highest values being 95.9, 77.0, 72.6, 65.9, and 58.7 for Italy, Australia, South Sudan, Argentina, and the USA, respectively (Figure 4). This means that in these regions, the information on the type of crop you are growing can influence your water scarcity results even more so than the country your activity is taking place in.
The characterization factors for rice, maize, soy, and wheat are shown on a map in Figure 5, and all CF values for all crops are available in SI. For example, as can be seen in Figure 5, in Argentina, the potential impacts in m 3 world-eq per m 3 of irrigation increase from rice (0.45), soy (3.9), and maize (9.5) to wheat (17.4). These differences represent the spatial variability of the location of the different cultures and indicate that wheat is grown in regions of higher scarcity than maize and so on.  Figure 3. Variability of crop-specific AWARE factors for selected countries (m 3 world eq/m 3 ).

Water Scarcity Footprint of Selected Crops
The characterization factors shown in Section 3.1 were used to calculate the water scarcity footprint of maize, rice, soybean, and wheat per producing country, in m 3 world-eq per ton of crop produced in country x. The results are shown in Figure 6. Within the same crop, the WSF can vary by more than three orders of magnitude across countries, following a similar trend as the CF themselves. Within one country, the results may also vary by as much, depending on the crop. This is the case for Bulgaria, Uruguay, Tajikistan, and France, which vary by three orders of magnitude or more. Uruguay, however, always presents values lower than 500 m 3 word-eq/ton in Figure 6, as the WSF in m 3 word-eq/ton of the four crops varies between 0.4 for wheat and 378 for rice, a value still in the lowest color-scale category. India, for example, varies by two orders of magnitude, as it produces soy with 348 m 3 world-eq per ton and wheat with 35,770 m 3 world-eq per ton. All values for these four crops are available in Supplementary Materials SI. When the CF for a specific region (e.g., Rice in Greenland) are not shown, it is either because of a lack of production or lack of specific data. The AWARE method is a relative one, so values are meant to be compared with each other and lower CF represent regions where the contribution to water scarcity from the cultivation of these crops is lower.
Water 2020, 12, x FOR PEER REVIEW 8 of 11

Water Scarcity Footprint of Selected Crops
The characterization factors shown in Section 3.1 were used to calculate the water scarcity footprint of maize, rice, soybean, and wheat per producing country, in m 3 world-eq per ton of crop produced in country x. The results are shown in Figure 6. Within the same crop, the WSF can vary by more than three orders of magnitude across countries, following a similar trend as the CF themselves. Within one country, the results may also vary by as much, depending on the crop. This is the case for Bulgaria, Uruguay, Tajikistan, and France, which vary by three orders of magnitude or more. Uruguay, however, always presents values lower than 500 m 3 word-eq/ton in Figure 6, as the WSF in m 3 word-eq/ton of the four crops varies between 0.4 for wheat and 378 for rice, a value still in the lowest color-scale category. India, for example, varies by two orders of magnitude, as it produces soy with 348 m 3 world-eq per ton and wheat with 35,770 m 3 world-eq per ton. All values for these four crops are available in SI. When the CF for a specific region (e.g., Rice in Greenland) are not shown, it is either because of a lack of production or lack of specific data. The AWARE method is a relative one, so values are meant to be compared with each other and lower CF represent regions where the contribution to water scarcity from the cultivation of these crops is lower.
. Figure 6. Water scarcity footprint per crop, shown for maize, rice, soy, and wheat, in m 3 world-eq/ton.

Discussion and Conclusions
The results of the quantification of the crop-specific CF proved that there is a high variability among crops in selected countries, suggesting that the use of these CF is a better proxy of water scarcity in a region when compared to the national generic CF counterparts ( Figure 5). Generic unknown CF will either overestimate or underestimate the water scarcity impacts of a specific activity. As a consequence, it is recommended that one uses sector-, commodity-, or crop-specific CF in the case where no watershed inventory-specific data are available. This provides a basis for improving understanding of the relative exposure of different crops to water scarcity-related risks. The tendency of over-or under-estimating water scarcity impacts was also proved with reference to generic agri-factors: for example, in the case of Italy, soybean CF values were 3.38 vs. 48.5 of the agri-CF and 43.3 of the general unknown-CF.
The relevance of using more specific CF rather that the general country-wide unknown factor has also been confirmed by other studies. Northey et al. [17] presented commodity-specific CF for the mining industries, proving that the use of sector-specific CF is more representative than the general country-wide unknown factors. It has to be noted that in this last work, a different weighting approach based on a production-weighting average was used to determine the specific CF. This is because the specific commodity water use data were not available and the authors assumed that water use was uniform across production.
The data used in this paper for specific crop water consumption are based on GCWM which model irrigation based on the years 1998-2002. Cropping patterns have certainly changed since then, as well as climatic data, and one could expect that updated data would influence the results. Already, when comparing GCWM with WaterGap data, which bases itself on production patterns of 2010, and climatic data over a 50-year period (1960-2010), discrepancies were observed. In addition to the models themselves, the climatic period, cropping patterns, and other input parameters, the data resolution and aggregation to country level can also introduce discrepancies, which could be responsible for these differences. The point to remember is that crop-specific water consumption is model-based, and that models are affected by their hypothesis, as well as input data validity and representativeness. Lastly, models quantify water consumption for optimal crop growth, whereas, in reality, more or less than optimal water may be used based on the water availability, technology used, regions, etc.
Crop-specific CF have a national geographical and annual temporal resolution. Considering that crop growth can be significantly influenced by local climatic conditions that change through space and time (this is particularly relevant, for example, for temperate climates), future research should focus on improving the temporal representativeness of the assessment. This could be achieved, for example, by weighting the monthly water use data of specific crops; however, no watershed crop-specific data with such a temporal resolution are available, according to the knowledge of the authors.
The crop-specific CF can also be interpreted as indicators of the contextual water risk of crop-specific productions and related industries. It can be stated that regions with higher crop-specific CF values can in fact encounter higher risks of water scarcity-related issues with that specific crop (e.g., more difficult access to freshwater supplies, increasing tension with competing users, and higher costs). This result is also in line with that of Northey et al. [17], who presented similar maps for the mining industries' commodities.
Focusing on the impact assessment presented in Section 3.1, the results showed that the trend of CF follows the overall values of WSF, and high variability can be identified among different locations and among different crops within the same country. This is in line with the fact that local water scarcity conditions influence the production of specific crops and thus the related water consumption (e.g., in more arid regions), as well as the fact that different water crops have different water needs.
Considering the result of the study, further improvement of the method can be outlined. In fact, the CF presented in this paper represent the first attempt to provide crop-specific factors. Future research, for example, should adopt other data validation references, such as specific field datasets and more recent agriculture management approaches, to determine local water use and climatic conditions [23,24].
The developed crop-specific CF could help crop producers increase the transparency and reporting of their water scarcity footprint results [18], as well as facilitate meaningful interpretations of data. To allow this, the developed crop-specific CF could be integrated into commercial LCA software. Practitioners should therefore be able to select crop-specific CF when performing impact assessment.

Funding:
We acknowledge the financial support of the industrial sponsors supporting WULCA during this work: Hydro-Québec, Veolia Environment, Danone, Unilever, Cottons Inc., Engie, and L'Oréal.