A temporal analysis of the consequences of the drought regime on the water footprint of agriculture in the Guadalupe Valley, Mexico

Changes in water availability have a substantial impact on the sustainability and maintenance of agriculture, with water footprint (WF) being a robust methodology to assess these transformations. The Guadalupe Valley is one of the places with the highest agricultural production in Mexico. Despite its semi-arid climatic conditions, it provides high-quality crops that are well-positioned in the world. The historical trend of rainfall and temperatures between 1987 and 2017 was analyzed to identify climatic patterns in the territory. Through the calculations of the water footprint of Grapevine and Olive crops, the sensitivity of the crops to recurrent water deficit and their adaptation in their yields to drought episodes was identified. The reduction in precipitation and occurrence of extreme temperatures have contributed significantly towards augmenting crop evapotranspiration and, consequently, intensifying crop irrigation demands. As a result, there has been an apparent increase in the consumption of WFagricultural since 2007. Thus, the period of highest WFagricultural consumption was 2014 (Extremely dry), as opposed to 2011 (Very wet). In particular, the lowest WFgreen consumptions were observed in extremely dry years, that is, > 20% of the WFagricultural intensifying drought events. Therefore, these periods were compensated with higher uses of WFblue and WFgray, which are inversely correlated with precipitation, where vine crops consume 73% more WFagricultural compared to olive plantations, showing greater interannual variability. These results contribute to analyzing the temporal evolution of water consumption for agriculture, providing a basis for rational water use strategies.

Similarly, blue water serves a multitude of functions, including but not limited to supply, transport/chemical loading (nutrients and pollutants), regulation/control, and production 24,25 ; furthermore, irrigation accounts for about 70% of the world's blue water use 26 .
A complex interaction of these variables (climatic and productive) has led water-scarce regions to produce more water-intensive but higher-value crops.In regions where water resources for agricultural purposes are scarce, a potential reduction of 56% in water usage can be achieved by cultivating crops better adapted to the area's green and blue water availability 27 .
At the same time, it is essential to acknowledge that soils in arid regions are fragile (eroded or degraded).This is because they have less natural primary production, organic matter, nutrients, water retention capacity, and fertility.As a result, plant growth and biomass are negatively affected.Therefore, before adding excessive amounts of chemical fertilizers, it is essential to assess the volume of grey water needed to assimilate the pollutants generated during the cultivation process 28 .
Considering the previously mentioned, prioritizing the evaluation of water abstractions is crucial for understanding the relationship between productive activities, the water cycle, and the pressure on water resources to ensure sustainable development in agriculture.Comprehensive indicators such as estimating the water footprint can help achieve this goal 20 .
In addition, the agricultural water footprint (WF agriculture ) is a measure that indirectly evaluates the impact of aridity by considering the local water demands and supplies 23 .It considers atmospheric water requirements such as evaporation and evapotranspiration, as well as temperature, precipitation, and irrigation.The measure also considers soil moisture, soil characteristics, and crop yield 21,29 .
WF agriculture studies distinguish between the use of different water sources-blue (WF blue ), green (WF green ), and gray (WF gray )-where the amount of precipitation before the growing season can reduce the need for irrigation during crop growth.This can lead to a significant decrease in the use of blue water up to 96% for some crops 8,20,30 .
However, during drought events, scarcity affects water flows equally, becoming crucial to increase irrigation to compensate for crops' deficit or green water requirements (CWU green ).Nevertheless, more than blue water requirements (CWU blue ) are required, constituting deficit irrigation at the expense of environmental flows and/ or groundwater reserves [30][31][32] .
Thus, to contribute to land management programs for irrigated crops in water-stressed basins, it is crucial to have accurate information on how agricultural water footprints have changed over time.This requires refined estimates over extended periods, which can be used in decision-making to reduce the latency of input data and improve future projections.By incorporating the influence of intra-and inter-annual climate variability in drought mitigation or action plans, we can make more informed decisions based on water footprint results 33 .
Our research lies in understanding water use in agriculture, considering that in Mexico, approximately twothirds of its territory is located in arid or semi-arid areas facing natural water scarcity, where the pressure on water resources has experienced a severe drought condition, which has persisted in the last 27 years, associated with El Nino Southern Oscillation (ENSO) and by the increase in water demand (up to 35%), especially after the mid-1990s [34][35][36] .
The water deficit is particularly severe in northern Mexico, where semi-arid areas cover about 50% of the total surface 37 .The central forecasts correspond to a 10-20% decrease in mean annual precipitation, a 1.5 to 2.5 °C temperature increase, intensified incidence of heat waves, and a reduction in aquifer recharge 38 .
The Valle de Guadalupe is a semi-arid rural area in the Ensenada municipality in Baja California.Despite the critical water situation, it is one of Mexico's highest producers of olives and grapevines, with about 46.6% of the total cultivated area dedicated to high-quality, high-value grapes 39 .The area produces 75% of Mexico's wine and 30% of its olive oil, making export crops and by-products the primary economic and tourism livelihood source.However, the limited water resources and agricultural land make the continuity of these industries vulnerable.
Different crops require varying amounts of water, and it is crucial to assess their water footprint to understand the factors that influence these differences.This helps promote the efficient use of water in agricultural production.To this end, this study examines the water footprint for olive and vine cultivation, two of the most predominant crops.Additionally, the impact of climate variability on agricultural WF agricultural consumption is analyzed through temporal analysis.This analysis is based on the climatic variables recorded for very wet, wet, normal, dry, and arid years.

Determination of climatic variability
The interannual precipitation of the Guadalupe Valley for the 29 years analyzed shows a historical mean of 250 mm, with 62% of the analyzed years below this mean.The years 1989, 1999, 2002, 2007, 2009, and 2014  recorded amounts between 59-133 mm, which is significantly lower precipitation compared to the other years (F (30,341) = 1.03; p = 0.001) (Fig. 1a).That is, years with a deficit between 76 to 47% of precipitation regarding the 1987-2017 average, with a negative trend starting in 2005 (Fig. 1a).
In the annual cycle, the behavior of seasonal cycles and extreme values in daily and monthly observations indicate greater variability and rainfall amounts during the winter period (December-March), where about 74% of precipitation occurs (ANOVA F (11,360) = 16.8, p = 0.001), with a prolonged dry period in June-August (Fig. 1b).
According to the percentile analysis, extreme conditions of very dry years were observed in 1989, 1999, 2002,  2007, 2009, and 2014 with 0-20% of the precipitation percentile, with 1989 being the most extreme.Contrary to the years 1987, 1992, 1993, 1998, 1998, 2004, 2010, and 2011 which showed a very wet condition, between 80-100%.Only five years of the analyzed series were classified as normal (Fig. 2 3a).The evolution of the monthly mean temperature shows significantly higher values in the summer period (July, August, September) compared to other months of the year (ANOVA F (11,360) = 170.2,P = 0.001) (Fig. 3b).
Lastly, the precipitation and temperature data series are moderately inversely correlated (r = − 0.66; p < 0.05), i.e., precipitation is lower in the months where temperature increases, decreasing water availability in the Guadalupe Valley in the warmer periods.
In summary, high climatic variability of the Guadalupe Valley can be observed from 1996, increasing both 'very dry' and 'very wet' years, thus intensifying extreme weather conditions.As a result, the occurrence of extreme meteorological droughts has been evident every 3 to 5 years since 1999.
The annual precipitation trend has decreased over the last 12 years; however, 1989 stands out for its low rainfall and high temperatures.Meanwhile, the changes identified in the annual temperature result in an amplification of the thermal oscillation.Similarly, the IR of the crops is inversely correlated with precipitation and effective precipitation (r = − 0.84; p < 0.05), which were higher in 2010-2011.With the increase in precipitation, ETo and ETc values were drastically reduced, combined with the decrease in temperatures, as in the case of 2010-2011.In contrast, in 2014, although it was not the only extremely dry year with a precipitation deficit, the extreme temperatures recorded led to an increase in ETc and IR.
The WF agricultural of the crops analyzed in the Guadalupe Valley presented estimated consumptions between 645 and 1121 m 3 /ton.Although there was some interannual variability, a slight incremental trend was identified in the 11 years analyzed, and the year 2014, which had extremely dry rainfall, had the highest WF agricultural consumption.On the other hand, 2011, which was a very wet year, had the lowest consumption estimated (as shown in Fig. 4).
The highest consumption of WF green was observed in 2010, 2011, and 2015 (Extremely humid and normal), with 35%, 38%, and 27% of the total WF agricultural , being found in these periods the precipitation and soil moisture necessary for crop growth, a variable that is inversely correlated with IR (r = − 0.78; p < 0.05).
Grapevine crops show high interannual variability, that is, an average of 3% ± 22.5 in WF agricultural consumption, with a maximum of 922.1 m 3 t −1 in 2014 (Extremely dry), presenting a higher use of WF green in 2010, 2011, and 2015, which is above the 170 m 3 t −1 uses of WF blue in 2007, 2014, and 2017, exceeding 500 m 3 t −1 , and WF gray in 2007-2014, exceeding 150 m 3 t −1 .Maximum water uses of CWU green crops were evinced in the 2011 period (Extremely humid) (1883 m 3 ha −1 ) and for CWU blue in 2014 (Extremely dry) (4601 m 3 ha −1 ) (Fig. 6a).
On the other hand, olive crops presented a lower interannual variability of WF agricultural , that is -1% ± 5.3, with amounts up to 199 m 3 t −1 recorded in 2014 (Extremely dry).The highest consumptions of WF green were in 2010, 2011, and 2013, above 55 m 3 t −1 ; WF blue in 2014-2017 with values above 95 m 3 t −1 ; and WF gray in 2016 with more than 74 m 3 t −1 .In addition, in the 2010 period, the maximum uses of CWU green (1414 m 3 ha −1 ) and 2014-2017 of CWU blue (> 2000 m 3 ha −1 ) were estimated (Fig. 6b).
Crop conditions in the Guadalupe Valley reported yields below 10 t −1 ha −1 for grapevine 39,40 and 12 t −1 ha −1 for olive 41 .Therefore, the net crop production and annual WF agriculture consumptions, when compared to the amount of water granted (25 × 10 6 m 3 according to the Drinking Water Committee), it is possible to identify the annual use of 36.4 to 82.7% of the available water (Table 1), where the maximum Total WF grapes were evinced in 2012, 2013, 2014 and for Total WF olives in 2010, 2013, 2016 (Table 1).Both crops accounted for an average of 55% ± 12 of the water available for agriculture in the study period.

Discussion
Currently, drought management comprises measures aimed at restricting water supply, reducing demand, and minimizing impacts.These measures have multiple effects on the maintenance and conservation of ecosystems, optimization of water supply, availability, and recharge within watersheds 4 .
In order to use resources efficiently, it is important to understand the consequences of climate variability on hydrological responses.In terms of precipitation in the Guadalupe Valley, in the last 30 years analyzed, 62% of the period is below the historical average of 250 mm, i.e., a deficit of 76-47% of precipitation, with a clear negative annual/seasonal trend since 2005.
The yearly variation in precipitation levels can be attributed to the impact of El Niño-Southern Oscillation (ENSO) and Pacific Warm Decadal Oscillation (PDO) events 42 .Our analysis has identified extreme conditions such as very dry or wet years and their alternation as evidence of these events.
Therefore, it is important to note that the transition from El Niño during 1998 (the rainiest year in the hydro-climatic record, classified as very wet) to the La Niña event (the following year) (Fig. 1), which marked the beginning of prolonged drought periods, in which extreme drought became even more severe later, that is, in the years 1999, 2002, 2006, 2006, 2007, 2009, and  El Niño is linked to the increase in surface water temperature in the central Pacific Ocean, inducing severe droughts in Indonesia, Australia, and northeastern South America 5,43 .La Niña, on the other hand, occurs when the surface water temperature drops off the coast in the Pacific Ocean, contributing to drier than normal weather in the North and South American regions, which also influences drought episodes, with a more complex impact on weather patterns than El Niño 5,43,44 .
Although El Niño-Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO) and Sea Surface Temperatures (SST) of the tropical Atlantic have played an important role in recent drought in northern Mexico 44,45 .However, global warming has increased the demand for atmospheric humidity and has probably altered circulation patterns, contributing to sustained episodes of drought in the northwestern region of Mexico, observed since 1994 34 , where recurrent dry years have significantly diminished water resources, in which high interannual climate variability through changes in precipitation and temperature, aggravating the present situation in the Guadalupe Valley 10 .
Temperature and precipitation are the main climatic factors that produce changes in soil water content and, therefore, directly regulate crop growth.By including these variables, the water footprint analysis reflects the water requirements in agriculture and the sensitivity of the results, where the precipitation deficit and extreme increase in temperatures conditioned the increase in ETo, ETc and, IR 46 .
Therefore, of the 11 years evaluated and considered as the present conditions (2007-2017), a slight incremental trend of water consumption was observed, despite the interannual variability, where the highest consumption of WF agricultural (grapevine + olive) was in the year 2014, classified as Extremely dry.This trend was also identified in the high WF blue consumption in 2007, 2009, 2014, and 2017, since during drought years, there was an increase in crop ET related to climatic variables and irrigation requirements.
In contrast, in the years classified as extremely wet and normal, i.e., 2010, 2011, and 2015, more WF green was consumed; being essential to include the evapotranspiration flux link between the ecological and hydrological systems, understanding that biological water use is inexorably linked to ecosystem productivity (Fig. 7).
The reported results, as described in Marston and Konar 25 in the California Valley, differences in WF agricultural during periods of drought are also attributed to reasons for changes in crop area and yields, increases in ETc, irrigation demand to compensate for precipitation deficits, and increased groundwater consumption due to reduced surface water availability 47 .
In these wine and fruit growing regions, in addition to the seasonal climatic changes described above, it is also necessary to include the influence of climate on an intermediate geographical scale (< 50 km) since these local climatic differences can cause considerable changes in crop physiology and phenology.Therefore, it is suggested to grow deep-rooted plants in drylands, due to the lower probability of suffering stress due to effects or changes in soil moisture, compared to shallow-rooted crops with less root coverage 48,49 .
However, crops can also modify their water needs through physiological adjustments, that is, in the regulation of stomatal conductance and leaf area transformation 50 , being crucial to the adaptation of plants to changing hydrological regimes, where their growth depends on whether the altered water supply is sufficient.Thus, the highest water requirements per hectare of crops were recorded in 2011 (1883 m 3 ha −1 ) for CWU green and 2014 (4601 m 3 ha −1 ) for CWU blue , with olive crops being the best adapted to interannual variability according to WF agricultural .
As WF agricultural is associated with both yield and seasonal water consumption, grapevine growing conditions in the Guadalupe Valley are limited by low yields, i.e., less than 10 t -1 ha −1 , compared to the main wine producers, for example, Italy, France, and Spain (over 80 t -1 ha −1 ).For olives, yields are even more extreme, reaching maximums of 23.40 t −1 ha −1 and minimums of 0.10 t −1 ha −1 , respectively 35,37 .
Excessive use of nitrogen fertilizers is common in sandy soils such as those in the Guadalupe Valley.As a result, the soil contains nitrate levels of more than 15 mg kg −1 , while nearby aquifers show levels of 10 mg L −151 .This has led to severe contamination problems by leaching and accumulation of macronutrients, especially during extremely dry periods when more irrigation is needed.As a solution, significant uses of WF gray were established to prevent further contamination 52 .
For agriculture to be a sustainable activity over time, it is necessary to include all approaches for mitigation and adaptation to water-scarce environments, from broadly tolerant crops to the introduction of climate and seed varieties, diversification, and crop calendars.In addition to the above, better information www.nature.com/scientificreports/ on droughts and climate variables, planned and comprehensive agricultural relocation, support for land use changes, and marketing strategies 53 become necessary, with productive land and fertile soil being the most important natural capital asset 27 .Therefore, soil preparation is essential, such as zero tillage and composting, the introduction of cover crops, intensive grazing management, and reduction in the use of chemical and synthetic inputs.Banwell 54 insist that such measures should consider the participation of various stakeholders, especially with a bottom-up approach.Indeed, only by modifying the cropping pattern could water scarcity be significantly reduced and production increased.According to Davis 56 , it is possible to reduce green and blue water consumption by 13.6% and 12.1%, respectively, only by growing the crops recommended for the area.First, however, farmers and decision-makers must be aware of the problem.
Since drought or water scarcity are dynamic, multiple indicators for their monitoring have been established 57 .Although initially, the water footprint is a tool for the evaluation of water consumption in production chains, in this study, it is possible to establish the sensitivity of crop water consumption results to changing weather conditions [58][59][60] , contributing to be the starting point for timely monitoring of water deficit and agricultural development in the area, for adequate preparation of responses according to the effects on the environment, agriculture, and local economy, understanding that the beginning and end of dry periods are perceived in retrospect.However, food security and water availability need immediate action 48,61 .
Lastly, we must not forget that the water resource is a priority requirement on which we depend, fulfilling its function through its cycle, storage, transport, hydroclimatic and ecological regulation, in which water flows, in all its states (blue, green, and gray), intervene in energy balances and biogeochemical cycles, for example, critical ecosystem services require 90% of global evapotranspiration for their functioning 62,63 .

Conclusion
The increase in precipitation in the Guadalupe Valley has a significant impact on the drastic reduction in ETo and ETc values, also influenced by the decrease in temperatures.However, this situation was contrary to the precipitation deficit and extreme temperatures that caused an increase in ETc and IR in crops.
A clear incremental trend in the consumption of WF agricultural was identified during the 11 years analyzed.Moreover, 2014 (Extremely dry) was the period with the highest WF agricultural consumption, as opposed to 2011 (very wet), where the lowest consumption was estimated.
The quantification of green water flow must necessarily be included in the evaluation of water resource security since the lowest consumption of WF green was observed in extremely dry years, that is, > 20% of the WF agricultural , due to the limitation in the availability of green water, intensifying and aggravating drought events.Periods compensated with higher uses of WF blue , corresponding to 64-68% of the WF agricultural , periods where the need for irrigation increases.www.nature.com/scientificreports/ The dissolution of fertilizers (nitrogen) decreases in periods of water scarcity, thus increasing the probability of contamination of environmental matrices, where the highest WF gray was determined in Extremely dry years with a consumption of over 200 m 3 /ton, a variable correlated with the irrigation requirement and inversely correlated with precipitation.
Grapevine crops consume 73% more WF agricultural compared to olive plantations and also showed greater interannual variability, in other words, less adaptation to climate variability.Nevertheless, the commercialization and water productivity of these crops is responsible for the continuity of their production.
The results obtained are a significant step towards reducing water waste in agriculture in arid regions with limited resources and fragile ecological environments.Despite the efforts made so far, such as the National Program Against Drought, the Mexican Drought Monitor, and the Programs for Drought Prevention and Mitigation Measures, these measures are considered insufficient.We must formulate strategies involving the rational use of agricultural water to tackle this issue effectively.

Site of study
The Guadalupe Valley basin (32°03'N-116°37'W) is located northeast of the municipality of Ensenada, Baja California, Mexico (2380 km 2 ).Its population reaches 6,924 inhabitants, concentrated in 0.4% of the basin's territory, 38.7% of the inhabitants live in poverty, and 35.5% are vulnerable due to social deprivation.Economic activities are mainly related to livestock, tourism, and agriculture 64 .
Similarly, the little of the land in this area is used for agriculture (cultivated lands), accounting for 2.7% of the total land area (6272 hectares).Out of this, 60% is dedicated to grapevine and olive production.Other land uses include scrub (94.3%), forest (1.8%), grassland (1%) and water bodies (0.02%) 64,65 as seen in Fig. 8.
The Guadalupe Valley is classified as a semi-desert region with a semi-arid Mediterranean climate, directly influenced by the Pacific Ocean coastal system, which is only 25 km away, exposing Guadalupe Valley to dense fogs, breezes, and marine humidity.Climatic conditions resulted in mild temperatures year-round, with an average of 18 °C annually from 1986-2016.Monthly temperatures range from 5.4 to 31.6°C 67 .
Precipitation is a consequence of the weakening of the North Pacific system in the winter, which allows the entry of convective and frontal systems; average annual rainfall amounts reach 298 mm (1992-2016) 68 .Precipitation variations are mainly associated with El Niño/Southern Oscillation and the Pacific Decadal Oscillation (PDO).These hydroclimatic phenomena bring unusual precipitation events to the region, as under warm conditions during El Niño and the PDO, the precipitation rate tends to increase.In contrast, under cold conditions during El Niño and the PDO, the precipitation decreases 45,69 .The research methodology was developed in 3 stages.In the first stage, the climatic behavior of the Guadalupe Valley basin was analyzed.During the second www.nature.com/scientificreports/stage, the database was developed, and the classification of the drought cycles was conducted.Finally, the crops' water footprint was calculated.

Determination of climatic variability
The meteorological input data were selected from 2 meteorological and pluviometric stations maintained by the Information System for Irrigation Water Management of Baja California (Fomento Agropecuario de Baja California) and the meteorological network CICESE-CLICOM (Center for Scientific Research and Higher Education at Ensenada).These stations are located in the Guadalupe Valley, that is Agua Caliente (32°04'N-116°37'W) and Porvenir (32°06'N-116°27'W) stations, considering a 31-year data (1987-2017 period).Monthly, seasonal, and annual precipitation and temperature data were used.Exploratory and consistency analyses of the data were carried out using the double-mass curve method 70 .Boxplots were used to visualize the behavior, which provided information on the trend of seasonal cycles and extreme values in the observations.The classification of dry, normal, and wet years was made using the percentile method, which defines five categories: very dry (0-20 percentile), dry (20-40 percentile), normal (40-60 percentile), wet (60-80) and very wet (80-100) 71 .

Calculation of the water footprint of olive and grapevine crops (WF agricultural )
In the WF agricultural analysis, climate variables were included with georeferenced information from 2007-2017, with primary values of precipitation (mm/month), temperature (°C), wind (km/day), insolation (hrs), humidity (%), and secondary values of effective precipitation (mm/month), evapotranspiration (mm/day), solar radiation (MJ/m 2 /day), as well as standardized soil properties data.
In addition, grapevine and olive crop data were used, such as Kc values, stages, and critical depletion fraction by development period, along with estimates of crop evapotranspiration and irrigation requirements.
To estimate the WF agricultural , the methodology proposed by Hoekstra 72 was used, which considers the sum of its components: WF blue , WF green , and WF gray , of the assessed crops grapes (WF grapes ) and olives (WF olives ) (Eqs.(1) and ( 2).
The determination of WF blue and WF green was made from Eqs. ( 3) and ( 4), where Y corresponds to crop yield (t ha -1 ) and CWU blue , CWU green green and blue crop water use, according to its source (precipitation or irrigation), expressed in m 3 ha −1 .
For the calculation of blue and green crop water use (CWU blue , CWU green ), they were estimated based on the relationships of Eqs. ( 5) and ( 6).In the equations, Σ represents the crop growth cycle, that is, which starts from sowing (day 1) and ends at harvest.In addition, Igp represents the length in days of each cycle stage.ET blue is blue water evapotranspiration, and ET green is green water evapotranspiration.We use a conversion factor 10 to convert the water depths in millimeters into water volumes per land surface area in m 3 /ha.

Green and blue water evapotranspiration of crops
The water demand of grapevine and olive crops was determined considering the crop water requirement (CWR) using the CROPWAT 8.0 software 73,74 .This software allows the calculation of water supply for different crop patterns, amount of irrigation according to their typology, and crop yield, both under irrigated and non-irrigated conditions 75 .
The crop evapotranspiration (ETc) was calculated based on irrigation requirement (IR), assuming that water losses due to irrigation remain and return to the basin.
Where ETc (mm/dec) is adjusted to decadal values (dec), in relation to irrigation efficiency throughout the crop growing season, under ideal growing conditions, it is assumed that the water requirements of the plantations are met 76 , so that the sum of ET green and ET blue equals ETc.
ET blue was estimated according to Eq. ( 7), an equation that included the determination of I.R. (Eq.(8). the program calculated Effective Precipitation (Peff) according to the USDA S.C. method.This parameter refers to the fraction of the total precipitation used to meet the irrigation needs of the crop.If (Peff) is higher than the total evapotranspiration of the crop, ET blue is equal to zero indicating that no additional water is needed.
The value of ET green was determined using Eq. ( 9), which considers the ETc values calculated in Eqs. ( 10) and (11).The crop coefficient, Kc, considers the crop characteristics and the average effects of soil evaporation.The reference evapotranspiration, ETo (mm month −1 ), was calculated using the Penman-Monteith method with the CROPWAT 8.0 program.Climate data based on the latitude and analyzed period was also considered in the calculation 29 .

Calculation of the gray water footprint of crops (WF gray )
The WF gray is defined as the amount of water necessary to assimilate the residues and dilute them so that their concentration remains within the quality ranges, according to the regulations 38 .WF gray was estimated according to Eq. ( 12), where AR corresponds to the applied amount of fertilizer, α represents the leaching and runoff fraction of the product expressed as a percentage, C max is the maximum acceptable concentration defined by quality standards, C nat is the natural concentration of the pollutant, Y is the yield of agricultural production, and the agent used corresponds to nitrogen for each agricultural crop (AR) 72 .
The leaching fraction of this agent was 10%, and the C max was 0.03 kg m 3 , based on NOM-001-SEMARNAT Mexican Emission Standard.The C nat was set at 0.00001 kg m 3 through the best-case scenario assessment according to the methodology of Hoekstra 72 .

Statistical analyses
Correlation analyses were performed for the hydrological variables according to the Pearson coefficient using STATISTICA 8.0.In addition, for the WF agricultural ANOVA tests were computed.

Figure 1 .
Figure 1.Precipitation behavior from 1987-2017 of the Guadalupe Valley.(a) Percentage anomaly of total annual precipitation: blue bars represent positive values (surplus), and red bars represent negative values (deficit) and accumulated frequency of annual precipitation in both climatological stations Agua Caliente and El Porvenir, (b) Trends of monthly periods (boxplots): the horizontal line and square inside the box indicate the median and mean.The lower and upper ends of the box correspond to the quantiles 0.25 and 0.75.Lower and upper whiskers: quantile 0.05 and 0.95, respectively.

Figure 2 .
Figure 2. Climate classification according to the precipitation thresholds in the Guadalupe Valley.

Figure 3 .
Figure 3. Evolution of the average temperature of the Guadalupe Valley.(a) Absolute anomalies of the average annual temperature for the 1987-2017 period (red bars indicate warming, and blue bars indicate cooling) and annual average, (b) monthly average.

Figure 6 .
Figure 6.Annual variation of WF agricultural consumption of crops: a) grapes b) olive analyzed in 2007-2017 for the Valle de Guadalupe, Mexico.

Figure 7 .
Figure 7. Summary of the agricultural water footprint of the Guadalupe Valley (Reference: Author's elaboration).

Figure 8 .
Figure 8. Location, land use, and water bodies identified in the study area, Guadalupe Valley, Mexico (Reference: ArcGIS 9.3 software developed by ESRI (Environmental Systems Research Institute), https:// www.esri.com) 66 .

Table 1 .
WF grapes and WF olives consumption according to annual production and water availability in the Guadalupe Valley, Mexico.Source: SIAP, IMIP.