Supply and Demand Analysis of Water Resources. Case Study: Irrigation Water Demand in a Semi-Arid Zone in Mexico

To sustainably use water resources, it is important to quantify water availability in a certain region. Due to climate change, population increase, and economic development, water demand increases continuously. Consequently, the difference between supply and demand of water becomes a significant issue, especially in arid and semi-arid regions. In this research, the Soil and Water Assessment Tool (SWAT) model has been applied to the Guadalupe river basin, to assess supply and demand analysis of water resources in this area, specifically for the irrigation of agricultural crops and municipal uses. From the land use, soil type, and terrain slope maps, 763 Hydrostatic Release Units (HRU) were estimated, distributed in the diverse relief types making up the basin, featured by mountains, hills, plateaus, plains, and valleys. For the crop area, 159 HRU were found with the three slope classification types, where 57 HRU represent 91% of the cultivated area on slopes, from 0 to 15%, located in the Ojos Negros and Guadalupe Valleys. The Soil Conservation Service method (SCS) was used to estimate the average monthly runoff and soil moisture content. As a result, water resource parameters related to the supply were determined with this, e.g., runoff, aquifer recharge, flow, infiltration, and others. Crop coefficient values (Kc) were used to determine crop evapotranspiration (ETc), to estimate the water demand of these for each month, using the multi-year monthly average reference evapotranspiration (ETo) calculated with the SWAT model. Overall good performance was obtained considering average monthly discharges data from the Agua Caliente gauging station. The model was calibrated, modifying the parameters chosen according to sensitivity analysis: SCS curve number, base-flow factor, ground-flow delay, and the threshold for return-flow occurrence. The Soil and Water Assessment Tool–Calibration and Uncertainty Programs SWAT-CUP has different goodness-of-fit indicators for the model e.g., determination coefficient (R2), standard deviation of the measured data (RSR), Nash–Sutcliffe coefficient of efficiency (NSE), and others. Multiple iterations were performed, resulting in a ratio between the root mean square error and the standard deviation of the measured data (RSR) of 0.61, a coefficient of determination (R2) of 0.70, and a Nash–Sutcliffe efficiency coefficient (NSE) of 0.63. A supply–demand analysis of the volume generated by the runoff from the basin was performed using the method of estimating useful volume for a reservoir. It is observed in these results that only positive deviations were obtained, implying that runoff in this Agriculture 2020, 10, 333; doi:10.3390/agriculture10080333 www.mdpi.com/journal/agriculture Agriculture 2020, 10, 333 2 of 20 basin is not enough to meet monthly demand. Finally, the need to establish actions to ensure water management efficiency is highlighted, both for irrigation of agricultural crops and for supply to the region population.


Introduction
Water is a fundamental natural resource for socio-economic development, food production, and health and sanitation of communities. Similarly, population growth, pollution of water sources, and increasing demand are issues that threaten to exhaust the resource. Therefore, it is necessary to study hydrological cycle variables such as precipitation to be able to quantify liquid availability in a certain region and thus, seek strategies for the comprehensive management of water resources to guarantee socio-economic development, without affecting water availability to future generations [1][2][3].
These water scarcity issues are more evident in areas with arid and semi-arid climates such as those covering 49% of the Mexican territory since the distinctive limited rainfall in these regions aggravates the situation [4]. This occurs in the hydrographic basin of the Guadalupe River in Baja California, which has a precipitation regime with a dry season and another one of rains during the winter, recording an average of 254.6 mm per year [5]. Likewise, this said regime is affected by variability and climate change, contributing to the water resources restriction in the basin [6].
Not only climatic factors aggravate the scarcity issue occurring in the study area but also the demand increase for irrigation water as well, due to the rise in winemaking, the main regional economic activity representing 90% production of Mexican wines [7]. Since aquifers are the main source to meet water demand for irrigation and supply the community, there is currently a deficit in water volumes due to overexploitation and lack of planning. A situation that negatively affects the regional economy [8].
This shows the need to implement comprehensive management of water resources that promote coordinated development and administration of water, land, and related resources, to equitably maximize the economic results and social assistance without affecting the sustainability of essential ecosystems [9]. Therefore, it is first necessary to determine the resource availability, so that it can be quantified in all the phases and systems (soil, atmosphere, and surface water bodies) of the hydrological cycle, starting from precipitation as the input variable because it is measurable and observable [10,11].
Additionally, this set of phases obeys the water balance equation, which expresses the relationship between inputs, outputs, and internal variations in the soil moisture storage, occurring during a given period. The equation suggests a simple expression; however, its components quantification is usually complex due to lack of direct field measurements, as well as due to spatio-temporal variation thereof [12,13]. Therefore, various research has proposed water balance models through estimation of the hydrological cycle components [14][15][16]; among which those including transfer in the soil-vegetation-atmosphere system, known as Soil-Vegetation-Atmosphere Transfer (SVAT), stand out, which consider this system as a water continuum, made up of a main pathway transporting water between the soil and the atmosphere through the plant, and also consider possible effects of external factors on the physiognomy or physiology of plants [17,18].
In an SVAT system, the rise of the fluid through the plant tissues is due to the gradient of water potential (Ψ), the quantification of the water tendency to flow from one medium to another, passing from greater to lower potential areas (soil > Ψroots > Ψstem > Ψleaves > Ψatmosphere). In this journey, various obstacles must be overcome, starting with soil resistance (r soil ), since how the soil transfers water to the plants is not homogeneous as it depends on humidity and its properties, influencing the plant retention and water availability. Continuing the journey, the water faces root resistance (r root ), which opposes its passage from the roots to the stem. Likewise, stem resistance (r xilem ) to the water conduction through the xylem is presented through capillary action to reach the dermis of the leaves. Then, in the leaves. stomata resistance (r leaves ) must be overcome so that the plant can carry out the transpiration process. Finally, there is the aerodynamic resistance (r a ), corresponding to the impediment that the wind offers to the exit of water vapor from the epidermis of the leaves to the atmosphere [19][20][21][22][23].
One of the Soil-Vegetation-Atmosphere (SVAT) evaluation models used in various research is the SWAT (Soil Water Assessment Tool) model, developed by the United States Department of Agriculture to analyze the impact of soil management practices on the hydrological cycle, sediment production, and diffusion of chemical substances derived from agronomic practices; and to study parameters related to water quality, such as suspended algae biomass, biological oxygen demand (BOD), and dissolved oxygen concentration; thus allowing discovery of the current and future conditions of the water dynamics from a basin [24].
Recently, the applicability of this model has been verified in multiple scenarios worldwide, managing to develop algorithms for self-irrigation, improve crop yield, analyze precipitation behavior and evaluate water scarcity in a region, and other contributions [25][26][27][28][29]. Nonetheless, most SWAT applications in Mexico have been focused on the center and south of the country, in regions with different climatic conditions from those of the study area, estimating useful parameters in hydrological evaluations and the comprehensive management of the water resource [30][31][32][33]. Nonetheless, in the semi-arid zones, the study has been limited to the basins hydrological responses to the different climate change scenarios, and to the study of transboundary basins located between Arizona and Mexico [34,35], without proposing alternatives for decision-making regarding water resource management [36,37].
Because of the above, the performance of the SWAT model applied to the Guadalupe basin is evaluated, located in a semi-arid zone with a Mediterranean climate, determining the significant parameters (Curve Number II (CN2), Based Flow Alpha Factor (ALPHA_BF), Ground Water Dealy time (GW_DELAY) and Threshold Water Depth in the Shallow Aquifer for Flow (GWQMN)) in the calibration and validation of the model results. Subsequently, the model is calibrated from the data obtained from the Agua Caliente gauging station; to later analyze the scarcity problem of water resources by carrying out an evaluation of supply and demand for water in the area.
Consequently, in this research study, water balance of the semi-arid basin from the Guadalupe River was determined through the SWAT model, to assess the resource availability in the region and generate recommendations for the comprehensive management of water resources that could be replicated in regions with similar climatic conditions.

Theory
The SWAT model is based on the water balance equation (Equation (1) used to determine the input, output, and storage of water in a basin [38], where SW t is the final amount of soil water (mm); SW 0 is the initial soil water content (mm); t is the time (days); R day is the precipitation in the form of rain, snow, dew, or fog, received at the soil surface [24,38,39]; Q surf is the runoff that refers to the total water that reaches the main streams, composed of surface, sub-surface, and base flow [24,38,40]; E a is evapotranspiration that involves evaporation from the soil surface, a body of water or that intercepted by vegetation and transpiration from plants [10,41]; W seep is the quantity of water that trickles in the soil profile (mm); and Q gw is the quantity of underground flow [42,43].
From Equation (1), change in soil storage (∆SW) is estimated, defined as the sheet stored in a soil thickness, and corresponds to the difference between the initial and final content of soil moisture [10]. This way it is established that for any volume and for a certain time, the difference between the water inputs and outputs is conditioned by the stored volume variation. Likewise, SWAT is a semi-distributed model based on physical principles to quantify the amount of water inside a basin, since it allows consideration of the spatial distribution of precipitation and, in turn, the aspects that affect said distribution from the climatological information of the study area, the physical properties of the soil, topography, vegetation cover and soil management practices [15,38]. Therefore, as presented in Equation (1), precipitation is the main input variable for the SWAT model, data gathered from the meteorological stations, and with which the other parameters from the water balance equation, are estimated. From these input parameters, Hydrological Response Units (HRU) are generated, which are scattered areas of the same basin, disconnected from each other, but which maintain homogeneous features such as the slope (which determines the speed of runoff), elevation (precipitation volume varies with height), soil type, vegetation cover, and land use, so these units are expected to exhibit uniform hydrological behavior [30,44,45].

Materials and Methods
As previously mentioned, the study area corresponds to the Guadalupe river basin located in the northern state of Baja California, between latitudes 31 • 50 N, 32 • 16 N and longitudes 116 • 54 W, 115 • 52 W. It is a semi-arid zone with a Mediterranean climate, with an average annual temperature of approximately 16 • C [46] and average annual precipitation of approximately 254.6 mm [5]. It has an area of 2390 km 2 and is divided into three sub-basins, Ojos Negros located upstream, Guadalupe Valley in the central part, and La Misión sub-basin downstream that empties into the Pacific Ocean ( Figure 1). and, in turn, the aspects that affect said distribution from the climatological information of the study area, the physical properties of the soil, topography, vegetation cover and soil management practices [15,38]. Therefore, as presented in Equation (1), precipitation is the main input variable for the SWAT model, data gathered from the meteorological stations, and with which the other parameters from the water balance equation, are estimated. From these input parameters, Hydrological Response Units (HRU) are generated, which are scattered areas of the same basin, disconnected from each other, but which maintain homogeneous features such as the slope (which determines the speed of runoff), elevation (precipitation volume varies with height), soil type, vegetation cover, and land use, so these units are expected to exhibit uniform hydrological behavior [30,44,45].

Materials and Methods
As previously mentioned, the study area corresponds to the Guadalupe river basin located in the northern state of Baja California, between latitudes 31°50′ N, 32°16′ N and longitudes 116°54′ W, 115°52′ W. It is a semi-arid zone with a Mediterranean climate, with an average annual temperature of approximately 16 °C [46] and average annual precipitation of approximately 254.6 mm [5]. It has an area of 2390 km 2 and is divided into three sub-basins, Ojos Negros located upstream, Guadalupe Valley in the central part, and La Misión sub-basin downstream that empties into the Pacific Ocean ( Figure 1).  Figure 2 shows the structure of the SWAT model that allows use of a Digital Elevation Model (DEM), soil information and meteorological data to simulate the behavior of the components of the hydrological cycle in the study area and to calibrate the model results using the data from a gauging station.
For the modeling, the high-resolution digital elevation model-Light Detection and Ranging or Laser Imaging Detection and Ranging (LiDAR)-from the National Institute of Statistics and Geography (INEGI), corresponding to the Guadalupe river basin, was used as topographic information. With a 5 × 5 m resolution, where it is recorded that the study area elevations are in the  Figure 2 shows the structure of the SWAT model that allows use of a Digital Elevation Model (DEM), soil information and meteorological data to simulate the behavior of the components of the hydrological cycle in the study area and to calibrate the model results using the data from a gauging station.
For the modeling, the high-resolution digital elevation model-Light Detection and Ranging or Laser Imaging Detection and Ranging (LiDAR)-from the National Institute of Statistics and Geography (INEGI), corresponding to the Guadalupe river basin, was used as topographic information.
With a 5 × 5 m resolution, where it is recorded that the study area elevations are in the range of 0 to 1283 masl ( Figure 3). With this data, the physiographic characteristics of the basin were determined, including main currents, micro-catchment division, and slopes classification ( Figure 4).
Agriculture 2020, 10, x FOR PEER REVIEW 5 of 22 range of 0 to 1283 masl ( Figure 3). With this data, the physiographic characteristics of the basin were determined, including main currents, micro-catchment division, and slopes classification ( Figure 4).     The maps of use and type of land implemented were prepared from the INEGI Series V digital maps, the INEGI Edaphological chart, scale 1:250,000 [47], and were complemented with the land use data obtained from the State Inventory of Forest and Soils from Baja California [48]. According to Figure 5, it was found that the basin is made up of diverse types of vegetation, where desert scrub can be highlighted with 61%, pine forests with 15%, induced grassland with 11%, and agriculture The maps of use and type of land implemented were prepared from the INEGI Series V digital maps, the INEGI Edaphological chart, scale 1:250,000 [47], and were complemented with the land use data obtained from the State Inventory of Forest and Soils from Baja California [48]. According to Figure 5, it was found that the basin is made up of diverse types of vegetation, where desert scrub can be highlighted with 61%, pine forests with 15%, induced grassland with 11%, and agriculture with 11. Soil composition corresponds to 67% coarse texture, 32% medium texture, and 1% fine texture ( Figure 6).
Agriculture 2020, 10, x FOR PEER REVIEW 7 of 22 with 11. Soil composition corresponds to 67% coarse texture, 32% medium texture, and 1% fine texture ( Figure 6).   Hydrologic Response Units (HRU) were determined from the described data and use of map algebra, by superimposing the layers of land use, soil type, and slopes, estimating a total of 763 HRU in the basin area. Figure 7 shows the HRU that contemplate the most relevant land uses for this Hydrologic Response Units (HRU) were determined from the described data and use of map algebra, by superimposing the layers of land use, soil type, and slopes, estimating a total of 763 HRU in the basin area. Figure 7 shows the HRU that contemplate the most relevant land uses for this research, corresponding irrigation agriculture, to crops that require water supply from external sources (gravity, dripping, etc.) and temporary agriculture to that which depends on precipitation [49].
Agriculture 2020, 10, x FOR PEER REVIEW  8 of 22 research, corresponding irrigation agriculture, to crops that require water supply from external sources (gravity, dripping, etc.) and temporary agriculture to that which depends on precipitation [49].
On the other hand, the input variables required for modeling were: precipitation, maximum temperature, minimum temperature, solar radiation, relative humidity and wind speed; and for the climatological analysis, data from 17 stations located in the study area were used (Figure 1), which contain daily historical records for 32 years (1979-2011). It should be mentioned, also, that the records of the weather variables were obtained from the Quick Extraction of Weather Information tool (ERIC III), which contains historical databases from the national meteorological service (SMN); the Climate Computer Project (CLICOM), and from the national water commission (CONAGUA); in addition to the climatological database, The Modern-Era Retrospective analysis for Research and Applications (MERRA), attached to NASA.
From the land use, soil type, and terrain slope maps, (Figures 4-6), 763 HRU were obtained, distributed in the diverse relief types making up the basin, featured by mountains, hills, plateaus, plains, and valleys. For the crop area, 159 HRU were found with the three types of slope classification, where 57 HRU represent 91% of the cultivated area on slopes from 0 to 15%, located in the Ojos Negros valley and the Guadalupe Valley ( Figure 4). The soil type in the crop areas corresponding to planosol, regosol, and phaeozem soils of medium to coarse texture, show a humid hydrometric regime, very good for agricultural production [50]. Moreover, these valleys are settled on a geological basement of alluvial rocks, mainly made up of granular sediment deposits [51], allowing water retention and benefiting crops even more. This water retention is evident only in the rainy months (November to April) probably causing an annual runoff contribution to groundwater. The model results show a decrease in runoff in the months of low rainfall, mainly expected in the sandy soils where the crops, which represent 13% of the basin area, are located. Nonetheless, due to the rain characteristics in the area, in some years, rains often occur out of season benefiting the surface runoff found with humid soils and therefore, this should be exploited.  On the other hand, the input variables required for modeling were: precipitation, maximum temperature, minimum temperature, solar radiation, relative humidity and wind speed; and for the climatological analysis, data from 17 stations located in the study area were used (Figure 1), which contain daily historical records for 32 years . It should be mentioned, also, that the records of the weather variables were obtained from the Quick Extraction of Weather Information tool (ERIC III), which contains historical databases from the national meteorological service (SMN); the Climate Computer Project (CLICOM), and from the national water commission (CONAGUA); in addition to the climatological database, The Modern-Era Retrospective analysis for Research and Applications (MERRA), attached to NASA.
From the land use, soil type, and terrain slope maps, (Figures 4-6), 763 HRU were obtained, distributed in the diverse relief types making up the basin, featured by mountains, hills, plateaus, plains, and valleys. For the crop area, 159 HRU were found with the three types of slope classification, where 57 HRU represent 91% of the cultivated area on slopes from 0 to 15%, located in the Ojos Negros valley and the Guadalupe Valley ( Figure 4). The soil type in the crop areas corresponding to planosol, regosol, and phaeozem soils of medium to coarse texture, show a humid hydrometric regime, very good for agricultural production [50]. Moreover, these valleys are settled on a geological basement of alluvial rocks, mainly made up of granular sediment deposits [51], allowing water retention and benefiting crops even more. This water retention is evident only in the rainy months (November to April) probably causing an annual runoff contribution to groundwater. The model results show a decrease in runoff in the months of low rainfall, mainly expected in the sandy soils where the crops, which represent 13% of the basin area, are located. Nonetheless, due to the rain characteristics in the area, in some years, rains often occur out of season benefiting the surface runoff found with humid soils and therefore, this should be exploited.
Subsequently, once the input data had been organized, the selection of the Penman-Monteith methods for estimating potential evapotranspiration (PET) and the curve number method from the United States Soil Conservation Service [52] to determine runoff, were used; which allowed estimating the basin's monthly water balance. The SCS method is widely used to determine the average monthly runoff and soil moisture content, therefore, with this method, the water resource parameters related to the supply were determined, such as runoff, aquifer recharge, flow, infiltration, and others. Similarly, the Curve Number (CN) considers the relationship between soil cover, type of soil, and the response given in the runoff during a rain event, a situation affected by the previous soil humidity conditions (dry, humid, and normal).
The estimation of the basin water demand was made considering the requirement of the crops and population demand. In the case of crop demand, the crops present in the area were analyzed based on the National Agricultural Survey of the Insituto Nacional de Estadística y Geografía (INEGI by its initial in Spanish) and the general overview document from the Valle de Guadalupe, Baja California, carried out by the Ministry of Agricultural Development [50]. Table 1 lists the crops located in the study and planting areas.
Subsequently, crop coefficient values (K c ) shown in Table 1 for each crop established in the area, were used to determine crop evapotranspiration (ET c ), to estimate the water demand of these, for each month, by means of the multi-year monthly average reference evapotranspiration (ET o ) calculated with the SWAT model. The K c values were gathered from the document of the Information System for the management of irrigation water in Baja California, the entity that made the estimation based on the FAO methodology and by means of the BISms application, with which the irrigation schedules for these crops are gathered [54,55].
The crop evapotranspiration was estimated with Equation (2), considering the crop phases, which were taken from the Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP by its initial in Spanish) [56].
Where ET c is the crop evapotranspiration in mm; K c , the crop coefficient, either in the initial, middle, or final stage (with no dimensions); and ET o is the reference evapotranspiration in mm.
Regarding population demand, the demand established by the norm for population consumption in Baja California [57] was used. This was calculated with Equation (3).
Where Q med is the average daily flow in l/s, D wd is the domestic water demand in l per capita per day, P is number of inhabitants and 86,400 is a conversion factor.
Regarding the model calibration and validation, the SWAT-CUP (Calibration and Uncertainty Programs) computer program was used, which belongs to the SWAT family of programs and has 4 different algorithms: Sequential Uncertainty Fitting version 2 (SUFI-2), Generalized Likelihood Uncertainty Estimation (GLUE), applied Parameter Solution (ParaSol) and Particle Swarm Optimization (PSO). The SWAT-CUP was developed by the Swiss Federal Institute of Aquatic Science and Technology and is a public domain program. In this case, the SUFI-2 algorithm was used, based on a Bayesian framework and which determines uncertainties through a sequential and adjustment process. In this algorithm, all sources of uncertainty are considered, such as model input, model structure, model parameters, and measured data. Furthermore, SUFI-2 performs a combined optimization and uncertainty analysis using a global search method and handles many parameters through Latin Hypercube sampling [58,59].
The SWAT-CUP has different goodness-of-fit indicators for the model such as the determination coefficient (R 2 ), the standard deviation of the measured data (RSR), the Nash-Sutcliffe coefficient of efficiency (NSE), and others. In this case, R 2 was used as an objective function for calibration. In this way, following the methodology used by Weber and Ocampo [60], a univariate manual pre-adjustment was performed on a set of 9 initially selected parameters, in this way it was possible to detect those parameters that reflect a greater contribution to the improvement of fit statistics. Thus, the 4 parameters selected for model calibration were identified, which are: CN2, ALPHA_BF, GW_DELAY, and GWQMN. It was found that the most sensitive parameter corresponds to the SCS curve number because it is the method used by SWAT to determine surface runoff ( Table 2). Once the most sensitive parameters have been chosen, the SUFI-2 algorithm varies these parameters within a predetermined range until there is a good fit with the observed data [58].

Results and Discussion
First, the hydrological parameters of the SWAT model used for the calibration and simulation uncertainty analysis were selected. The initial and calibrated values of the adjusted parameters are given in Table 2. The data chosen for the calibration correspond to the average monthly flow data for the 1979-1986 period from gauging station 1023 AGUA CALIENTE (Figure 6), as it is the only station with available information on flows within the study area. Multiple iterations were performed, resulting in a relationship between the root mean square error and the standard deviation of the measured data (RSR) of 0.61, a coefficient of determination (R 2 ) of 0.70, and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.63. Therefore, according to the reviewed literature, the calibration result can be considered satisfactory [59,61]. In Figure 8, the comparison of the series of observed data (recorded in the gauging station) with the simulated data after the calibration of the model can be observed. considered satisfactory [59,61]. In Figure 8, the comparison of the series of observed data (recorded in the gauging station) with the simulated data after the calibration of the model can be observed.  Figure 9 shows the diagram of the water balance of the study area, where precipitation corresponds to 304.4 mm, evapotranspiration of the plants to 217.7 mm, and runoff to 3.6 mm. Similarly, it is seen that the number of the basin average curve is 35.1 and for this reason, an average percolation of 99.9 mm is shown, equivalent to 32% of the average precipitation.  Figure 9 shows the diagram of the water balance of the study area, where precipitation corresponds to 304.4 mm, evapotranspiration of the plants to 217.7 mm, and runoff to 3.6 mm. Similarly, it is seen that the number of the basin average curve is 35.1 and for this reason, an average percolation of 99.9 mm is shown, equivalent to 32% of the average precipitation.  Figure 9 shows the diagram of the water balance of the study area, where precipitation corresponds to 304.4 mm, evapotranspiration of the plants to 217.7 mm, and runoff to 3.6 mm. Similarly, it is seen that the number of the basin average curve is 35.1 and for this reason, an average percolation of 99.9 mm is shown, equivalent to 32% of the average precipitation. Similarly, according to the results estimated from the water balance shown in Table 3 and Figure  7, the maximum rainfall in the basin occurs from November to March, a period in which 68.35% of the rains of the year occurs, corresponding to 208.1 mm, with February being the rainiest month with 56.1 mm and July, the month of least rainfall with 8.7 mm. Similarly, the calculated total runoff, representing the sum of runoff, sub-surface flow, and base flow, minus transmission losses and abstractions, was 105.3 mm, which is 34.59% of the water precipitated annually. Likewise, it is demonstrated that in the period from November to March, precipitation was 48.9 mm, equivalent to 46.49% of the total annual runoff. Regarding the calculation of storage, it was found that there is a  Similarly, according to the results estimated from the water balance shown in Table 3 and Figure 7, the maximum rainfall in the basin occurs from November to March, a period in which 68.35% of the rains of the year occurs, corresponding to 208.1 mm, with February being the rainiest month with 56.1 mm and July, the month of least rainfall with 8.7 mm. Similarly, the calculated total runoff, representing the sum of runoff, sub-surface flow, and base flow, minus transmission losses and abstractions, was 105.3 mm, which is 34.59% of the water precipitated annually. Likewise, it is demonstrated that in the period from November to March, precipitation was 48.9 mm, equivalent to 46.49% of the total annual runoff. Regarding the calculation of storage, it was found that there is a deficit during the months with the least rainfall (March to July), showing the need to implement measures that allow efficient use of the water available in the study area.
Once the SWAT modeling was performed, and with the results from the ET o , the evapotranspiration of the crops, ET c (Table 4), and the demands of existing agricultural crops and population, were estimated (Table 5). Table 5 shows that the greatest demand is focused on pine forests, the vegetation from El Chaparral area, and pastures. Although the crops with the highest amount of water demand are vine and olive tree, in line with the economy of the area that is focused on these two agricultural activities, which is the reason the Valle de Guadalupe is recognized for its excellence in wine production.  Analogously, knowing the population of 11,871 inhabitants in the Guadalupe basin and its water supply of 200 L per inhabitant, per day [57]; an average monthly demand of 0.071 Mm 3 was estimated. Irrigation demand was also calculated from the real evapotranspiration estimated in the modeling and the cultivation coefficient (K c ) [62], obtaining a crop demand (Table 5). Based on these results, a supply-demand analysis of the volume generated by the runoff from the basin was performed using the method of estimating useful volume for a reservoir, shown in Table 6. In these results (Table 6), it is observed that only positive deviations were obtained, which implies that the runoff in this basin is not enough to meet the monthly demand. Also, the value of each deviation indicates the water volume necessary to address it, the reason it can be said that the accumulated runoff would not be enough to meet the demand for a month, proving the deficit occurring in the basin.
It was also identified that between April to September, evapotranspiration is greater than the total precipitation from each month, e.g., in June the actual estimated evapotranspiration was 23 mm while precipitation was 10.7 mm, a phenomenon that may occur since evapotranspiration is a continuous process that happens day after day, regardless whether precipitation events occur or not.
Subsequently, the supply and demand of the water resource were evaluated by considering the volumes provided by the aquifer. The multi-annual monthly contributions of the aquifer are considered from the water concessions granted by Comisión Nacional del Agua (CONAGUA by its initial in Spanish) in the area and were gathered from previous studies [63].
The results in Table 7 show that there is still a deficit in the resource, where demand exceeds supply. It should be noted, as can be seen in Table 8, that during this period the highest average temperatures and the highest records of solar radiation are shown. This is because these months correspond to the Summer season when the amount of water evaporated from the soil exceeds that transpired by plants, because in most Mediterranean climate species, there is a decrease in transpiration as a strategy to avoid or withstand the stresses that act on the plant due to water deficit, causing a significant variation in optimal living conditions (stress). According to the factor causing it, this can be classified as thermal stress due to extreme temperatures, stress due to high solar radiation, stress due to the scarcity of water available in the soil (water stress), or stress due to pollutants and biotic factors (pests) [18,[64][65][66][67]. Within the study area ( Figure 5), the species present in the chaparral, coastal Rosetophile scrub, and in the oak forests, have developed adaptive physiognomic and physiological strategies to survive. For instance, the Arctostaphylos pungens and Salvia apiana present fuzziness on its leaves that reflect solar radiation, while species such as Adenostoma fasciculatum, Ceanothus greggii, Quercus spp., and Ceanothus spp. can present the ability of curling, wax segregation, size reduction, senescence, and vertical orientation in their leaves [49].
When stress disturbs cultivated plants, production is the first thing affected, as can be seen in vine crops. Water stress decreases the quantity and quality of the harvested grapes, resulting in losses in wine production, proving that plants under stress conditions interrupt vital functions [65].
During the Summer the soil dries decreasing its matrix potential or the water cohesion force to the surface of the soil particles, preventing plants from continuing to absorb water and reaching the point of permanent wilting or cell death, i.e., the water level in the soil is less than the minimum necessary for the plant to absorb it (field capacity). Also, the vegetation generates a readjustment in the osmotic potential, the one that allows water movement due to differences in solute concentrations, without stomatal closure, causing unnecessary water losses, greater absorption by the roots, and increasing the volume necessary for irrigation. Furthermore, if the temperature and solar radiation are very high, the opposite effect occurs, stimulating stomata closure and reducing respiration and photosynthesis. A fact that leads the plant to use its carbon reserves to perform its vital functions, which could affect the sweet taste of fruit and vegetables.
Considering the above and the importance of agricultural activities for the study area, it is necessary to propose alternatives to mitigate the effect of the water deficit on crops. Therefore, the plant should maintain a low evaporation rate and a transpiration level compatible with the water available in the soil, i.e., that transpiration does not exceed the water amount that its roots can absorb, reducing its height, stem, roots, leaf area, and, in cases such as the study area where the water deficit is recurrent, not killing crops. Because of this and given the influence of solar radiation and ambient temperature on the physiology of plants, one of the ways to decrease evapotranspiration is to reduce the amount of net radiation that reaches the surface of the crop or soil using alternatives that increase albedos, such as the use of shade nets over the crop, crop residue mulches, or plastic sheeting on the soil surface.
Also, there are other alternatives such as using crops in protected environments, the collection of rainwater in such a way as to guarantee the availability of the resource during a critical period of greatest consumption for plants, i.e., the flowering stage of each species, and the implementation of soil moisture content controls to avoid excess irrigation and runoff leading to the leaching of nutrients such as potassium (K), magnesium (Mg) or Nitrogen (N), and the deficiency of aeration that causes root suffocation and increased disease in crops [68]. Furthermore, implementation of drip irrigation, micro sprinkling, sprinkling, furrow flooding, among other techniques might be sought, so that the available resource can be used most efficiently according to the needs of each crop.
Regarding the water used for irrigation, the water resource must be used during times of maximum rainfall, storing it properly, taking into account that in climates such as the one from the area, evaporation can consume a large part of the stored water if storage structures are kept open and/or exposed to the sun.
It is still worth noting that there are different methods of rainwater use (water harvests) that allow, in one way or another, to optimize irrigation work depending on the crop type and the features of the area, such as micro-basins [69,70], mini runoff catchments, macro catchments and stone lines [71,72], and other rainwater harvesting (RWH) techniques [73]. In the case of the study area, implementation of furrows or blind inlets to the contour of the crops is recommended, since these can reduce the speed of surface drainage, the sediments transport and, in the same way, facilitate the recharge of aquifers to maintain soil moisture and optimize water absorption by plant roots [74][75][76].
Finally, the need to establish actions to ensure water management efficiency is highlighted, both for irrigation and for supplying the community from the region. Therefore, training programs can be implemented, promoting strategies for administration, self-managing, and of equality regarding the water resource in the basin [77]; identifying the micro-catchments in which it is feasible to build catchment structures that guarantee access to the resource during the months of less rainfall; and preparing and implementing a basin management plan that incorporates the variables of climate change and climate vulnerability.

Conclusions
The water balance from the Guadalupe basin was estimated, finding evidence that there are conditions favoring deficit in the area, such as low rainfall and high evapotranspiration during different periods of the year. The above is added to the increase in the water demand necessary to meet the community needs and to guarantee crop irrigation. Therefore, promoting a comprehensive water resource management that guarantees agricultural sustainability in addition to the challenges that this activity faces under climate change conditions, require strengthening and development of efficient systems based on new technologies that guarantee crop irrigation with the minimum amount of water required by plants, paying attention to avoiding losses of water or soil moisture due to high temperatures, as well as developing projects for diversification, reconversion, and relocation of crops seeking to only use the existing water availability. Similarly, information systems must be integrated into the area to allow the establishment of accurate forecasts of hydrometeorological events and take actions to counteract them. Governance of water resources must also be reinforced with actions to strengthen agricultural information systems on the availability and allocation of the resource to allow decision-making by following the basin reality. This is done by promoting the organization of the basin users, conducting training, and raising awareness about the efficient use and saving of the water resources. Consequently, investing economic resources in irrigation infrastructure, and rainwater harvesting and use in agriculture, is required, incorporating also, social inclusion indicators that allow the optimization of the amount of water present in the study basin.
In general terms, it was found that the SWAT model is appropriate for the spatio-temporal analysis of the variables that make up the hydrological cycle, allowing establishment of basin water balance.
Nonetheless, it should be mentioned that the quality of the results will always be subject to information availability. In the case of the Guadalupe basin, it was found that the information related to runoff is scarce due to the few gauging stations installed, so it is expected that for future research a greater amount of hydrometric data will be available for the calibration of the model.