The Modelling of the Evapotranspiration Portion of the Water Footprint: A Global Sensitivity Analysis in the Brazilian Serra Gaúcha

: Water footprints have been widely used to illustrate the consumption of water in many situations, for instance, in products, processes, or regions of interest. In this work, we analyzed—using a sensitivity analysis approach—the effect of some variables in the calculation of the water footprint in the viticulture in the Brazilian Serra Gaúcha (the major producing region of Brazilian wine). The classical Penman–Monteith model for evapotransporation was considered, with uncertainties in some parameters (dead mulch covering a fraction of the vineyard, maximum temperatures for some months, the altitudes and latitudes of the site). A sensitivity analysis was conducted using the SAFE toolbox under Octave framework. The results indicated that the the portion of the water footprint corresponding to evapotranspiration is more sensitive to the values of the mulch-covered fraction and the altitude of the site in comparison with the latitude and the maximum temperatures.


Introduction
Water consumption in the viticulture is gaining more attention every year [1,2].Thus, it is necessary to understand the factors that affect it and in which ones we can interfere to reduce water use.In recent decades, Brazil has been receiving attention for its capacity to produce wine [3].Furthermore, the most traditional and main wine hub in Brazil is the Brazilian Serra Gaúcha [4][5][6].
A well-established measure of water consumption is called water footprint.The water footprint concept was introduced by Hoekstra and Hung [7] and represents the consumption of water by a crop, product, region, or nation (for instance).The water footprint is usually classified as green, blue, and gray water.In essence, green water indicates the portion of water arising from rain and snow, whereas blue water corresponds to the water coming from irrigation.Finally, the gray portion refers to the necessary quantity of water employed in dilution processes (in order to achieve the concentration of each pollutant in a river, for instance).In the wine industry, the study of the water footprint is fairly recent, with calculations presented for New Zealand, Italy, South Africa, and Portugal, for instance [8][9][10][11].Water consumption in wineries and in vineyards must be considered.Water is consumed through evapotranspiration, irrigation, and the use of fertilizers and agrochemicals, among others [9].
An important part of the water footprint calculation of a crop is usually conducted considering the evapotranspiration concept.Evapotranspiration is a combined process of evaporation and transpiration by a crop.The initial step in the calculation of evapotranspiration is the determination of the reference evapotranspiration (ET 0 ).The reference evapotranspiration is commonly calculated using the Penman-Monteith model [12], which is a well-established methodology.Using a small set of parameters (latitude, altitude, wind speed, temperatures of the site), this model permits the evaluation of the reference evapotranspiration ET 0 .The crop evapotranspiration (ET c ) is calculated using the value of ET 0 and a crop coefficient, which considers the different growth stages of a particular crop and the existence of dead mulch (for instance) [12].Clearly, the different quantities employed in the calculation of the ET 0 and ET c parameters are subject to uncertainty (for instance, temperatures and wind speeds), and the main objective of a sensitivity analysis (in this context) is to evaluate the effect of these uncertainties on the calculated values of evapotranspiration [13][14][15].Furthermore, the productivity of the crop is also subject to uncertainty, which can have a direct impact on the calculation of the water footprint.

Previous Studies
Many studies deal with sensitivity analysis using the Penman-Monteith model for the reference evapotranspiration, considering different parameters and locations.Beven [16] studied the sources of errors and presented a sensitivity analysis study for the Penman-Monteith equation for three meteorological sites in England and Wales.The author observed a higher dependence of the vegetation type on the evapotranspiration, in comparison to climatic differences between the three sites analyzed.Ndiaye et al. [14] analyzed the effects of maximum/minimum temperatures, solar radiation, wind speed, and maximum/minimum relative humidity on evapotranspiration in different localities in Burkina Faso.Irmak et al. [17] conducted a sensitivity analysis study of several parameters (maximum/minimum wind speed, maximum/minimum air temperatures, vapor pressure deficit, and solar radiation) in some regions of the United States of America, concluding that ET 0 is more sensitive to variations in the vapor pressure deficit.Debnath et al. [18] conducted a sensitivity analysis study for the reference evapotranspiration value in India, indicating that ET 0 is more sensitive to solar radiation.Furthermore, Debnath et al. [18] noted a linear relationship between the value of ET 0 and the variables under consideration (maximum and minimum temperatures, solar radiation, relative humidity, and wind speed).Biazar et al. [19] conducted a sensitivity analysis study for the reference evapotranspiration ET 0 , using different models, in a humid region of Iran.The findings of Biazar et al. [19] indicate a high sensitivity of ET 0 to the maximum temperatures and number of sunshine hours, and some lower sensitivities with respect to the minimum temperature and wind speed.Arunrat et al. [20] analyzed the impact of climate change on the water footprint for rice, cassava, maize, soybean, and mung bean in Thailand.The authors considered changes in demography and economics, as well as the greenhouse effect, in the studied scenarios (Shared Socioeconomic Pathways-SSPs).Arunrat et al. [20] concluded that the water footprint for the studied crops in the future will be closely related to the yields, considering the greenhouse effect.Furthermore, some crops (maize, soybean, mung bean, and cassava) appear as alternatives to rice during the dry seasons.Rossi et al. [21] studied the water footprints for olive production in Italy (considering the green, blue, and gray portions of the water footprint).The authors observed a predominance of the green water footprint in comparison to the blue and gray portions.Furthermore, it was observed by Rossi et al. [21] that the water footprint was severely influenced by the productivity, even in some years with low precipitation, which is in agreement with the results of Arunrat et al. [20] (i.e., the effects of some parameters in the reference evapotranspiration can be compensated by an increase in the productivity, considering the calculation of the water footprint).Arunrat et al. [22] presented a study analyzing the nitrogen, carbon, and water footprints of organic and conventional rice in Thailand in the 2018-2021 period.The authors highlighted that the water footprint for the conventional crop is higher than the value for the organic rice (mainly with respect to the gray water footprint).Recently, Yong et al. [23] presented a review study concerning the different sensitivity analysis techniques applied in the calculation of the reference evapotransporation.The results obtained by Yong et al. [23] indicated that the global sensitivity analysis studies are recent and restricted to a few countries (Australia, Malaysia, the United States, and China).A recent work by Sabino and de Souza [24] presented a study of global sensitivity analysis for the Penman-Monteith model to calculate the reference evapotranspiration in the state of Mato Grosso, Brazil, using the Sobol' method.
Specifically considering the calculation of the water footprint, Zhuo et al. [25] presented a study of uncertainty and sensitivity analysis in the water footprint for the production of maize, soybean, rice, and wheat in China (Yellow River Basin).Recently, Li et al. [26] conducted a sensitivity analysis study in the production of wheat, rice, maize, and soybean in China.In this study, the authors considered the value of the reference evapotranspiration ET 0 as one of the parameters.

Proposed Current Study
In this work, we propose a detailed sensitivity analysis procedure in the water footprint calculation for the viticulture in the Serra Gaúcha, Brazil, located in the Rio Grande do Sul state.As mentioned previously, the Serra Gaúcha is Brazil's main wine-producing region.We consider, as parameters, the altitude and the latitude of the site (which impacts the reference evapotranspiration) and the existence of dead mulch in the vineyard (affecting the crop evapotranspiration).In addition to these, the maximum temperatures for the months of October, November, and December are also considered as parameters.Obviously, there is a huge scope for choosing other parameters for the sensitivity analysis (for instance, wind speed and humidity).On the other hand, the main focus of this work is to detail the influence of some parameters which can be objects of decision (for instance, the place to set up a vineyard and the use of a soil management technique in order to minimize the water consumption).Thus, the results of the water footprint sensitivity analysis provide supplementary data for decision-making purposes.Moreover, considering the rising temperatures in the Rio Grande do Sul state (Brazil) [27,28], the effects of maximum temperatures on the water footprint are also taken into account.Moreover, Biazar et al. [19] indicated that the most sensitive parameter for ET 0 in a humid region in Iran is the maximum temperature, which justifies this choice.It should be noted that this study only covers the portion of the water footprint that refers to crop evapotranspiration in the vineyard.One must also bear in mind that the aim of this study is not to present comparisons between evapotranspiration models, despite the importance of such studies (see, for instance, [29][30][31][32]).
Although the analysis is specific to a particular region, the procedures used are applicable to any location worldwide.Here, we employ the SAFE toolbox (under Octave [33]) for sensitivity analysis developed by Pianosi et al. [34].Global sensitivity analysis tools are also available on other platforms; see, for instance, Borges et al. [35].Finally, it should be emphasized that the use of global sensitivity analysis techniques in water footprint calculations has not been frequently found in the literature.
The main innovation of this work is to combine global sensitivity analysis techniques with evapotranspiration models in order to verify the effect of geographical and climatic parameters on the water footprint, with the aim of informing decision-making.In addition to pointing out the direction of the effects of changes in each of the variables under study, our sensitivity analysis study will be employed to indicate possible offsetting effects between the input parameters of the model.

Modelling
In this section, the mathematical models under consideration are detailed, namely the reference evapotransporation model (Penman-Monteith model) (ET 0 ), the crop evapotranspiration model (ET c ), and the water footprint model (WF).

The Reference Evapotranspiration Model
The reference evapotranspiration (ET 0 , in mm/day or mm/month, depending on the time step employed) was represented by the Penman-Monteith model [12]: where R n is the net flux in solar radiation (MJ/m 2 /day) , G is the heat flux in the soil (MJ/m 2 /day), T is the mean air temperature at 2 m height ( • C), e 0 s is the saturation vapor pressure (kPa), e a is the actual vapor pressure (kPa), ∆ is the slope of the vapor pressure curve, γ p is the psychrometric constant (kPa/K), and u 2 is the wind speed at 2 m (m/s).
In order to clarify the objectives of the sensitivity analysis study, the main terms of Equation ( 1) must be detailed.Even though this equation is well known, a detailed analysis of the model permits a clear understanding of the nonlinearities of the equations involved in the problem.Therefore, the calculations of e 0 s , ∆, γ p , R n , and G will be explained in detail, as well as the correlated quantities.A detailed description of these parameters is presented in the Appendix A.
The Penman-Monteith equation was chosen due to the large number of scientific studies involving this model.

Crop Evapotranspiration Model
The crop evapotranspiration can be obtained using the value of the reference evapotranspiration, considering the different development stages of the crop, according to [12]: where Kc is the crop coefficient.Usually, three different stages are considered (initial, midseason, and late season stages), with three values for Kc.The FAO [12] recommendation for grape wines are Kc init = 0.30, Kc mid = 0.70, and Kc late = 0.45, where the subscripts init, mid, and late refer, respectively, to the initial, mid, and late stages of the development of the crop.On the other hand, some parameters can affect the values of Kc, such as the existence of dead mulch [36].According to Fonseca [37], the effect of dead mulch in the value of Kc init can be represented by the following expression: where f dm is the dead mulch coverage fraction.Kc init,dc represents the crop coefficient considering the existence of dead mulch.Thus, an obvious effect of the use of dead mulch is the reduction in the values of the crop evapotranspiration (ET c ).It must be emphasized that the use of mulch is a common practice in the Brazilian Serra Gaúcha, as pointed out by Rosa et al. [38], which justifies the study of the impact of this parameter on the water footprint.The use of organic mulching in vineyards to reduce ET c was also reported by López-Urrea et al. [39].Furthermore, Fonseca [37] pointed out that although the use of mulch is very common for wine production in the tropical regions of Brazil, this practice has also been observed in temperate regions.
In this work, it was assumed that the water demand by the crop was fully met.Thus, the water footprint can be calculated as [7] where Y is the grape productivity by hectare (ton/ha) and ET c is in mm/month.The factor 10 converts the evapotranspiration into m 3 /ha.The sum refers to the months of the year.We are not differing between the types of water footprint (green and blue).We considered the latitudes between the cities of Nova Araçá (−28.654549• ) and Gramado (−29.433653• ). Figure 1 presents an altimetric map of the Rio Grande do Sul state, in Brazil, as well as the localization of the Serra Gaúcha, the region under study.The main objective of these maps is to detail the ranges of latitudes and longitudes for the sensitivity analysis.In the Serra Gaúcha, the elevations are in the range from 206 m to 926 m (considering the existence of valleys in the region).In recent years, there has been an increase in the production of wine at altitudes above 700 m (for example, the so-called Altos de Pinto Bandeira wine).Furthermore, the soil management techniques in vineyards usually employ mulching practices [40].With respect to the dead-mulch-covered fraction, f dm , this parameter is in the interval [0, 1].Calculations were conducted using information from the municipality of Bento Gonçalves, Brazil, the main wine-producing city in the Serra Gaúcha.Data for temperatures, relative humidities, and number of hours of sunlight per day were obtained from Climate Data (Available online: https://pt.climate-data.org/america-do-sul/brasil/rio-grandedo-sul/bento-goncalves-1386/,accessed on 9 April 2024).Wind velocities were obtained from Weather Spark (Available online: https://pt.weatherspark.com/y/29709/Clima-caracter%C3%ADstico-em-Bento-Gon%C3%A7alves-Brasil-durante-o-ano,accessed on 9 April 2024).We considered variations of 3 • C in the maximum temperatures of the months of October, November, and December (an arbitrary choice in order to evaluate the effect of the maximum temperatures on the water footprint).These variations are far greater than the differences of the maximum temperatures in the region under study (even considering the differences in altitude and latitude).In order to subsidize the choice of this temperature range, Figure 2 presents the maximum temperatures in the Serra Gaúcha for October, November, and December in the last decade (2013-2023).We can observe that the maximum temperatures for these months are in the interval T base ± 3 • C, where T base refers to the temperature in the city of Bento Gonçalves.
A medium productivity value Y for grape production was assumed, Y = 10 ton/ha.Usually, three scenarios may be considered for grape production in the Serra Gaúcha: Y max = 12 ton/ha, Y avg = 10 ton/ha, and Y min = 8 ton/ha [41] (obviously, the effect of the productivity exhibits a direct effect on the water footprint, and a sensitivity analysis is not necessary for this variable).Furthermore, we used a typical conversion factor of 70% to calculate the number of liters of wine by ton of grape, i.e., a ton of grapes is capable of producing 700 L of wine.

Global Sensitivity Analysis Techniques
As mentioned previously, we employed the SAFE package [34,42], under Octave, to perform the sensitivity analysis.Here, some of the steps proposed by Noacco et al. [42] in conducting the sensitivity analysis were taken.

Sampling Strategy
The most popular sampling strategy is Latin Hypercube Sampling (LHS) [43].Some recent works employed LHS in sensitivity analysis studies dealing with evapotranspiration calculations [44][45][46].Considering that each coordinate of a vector is in the interval [0, 1], a Latin Hypercube Sample with N intervals shows a random distribution of N values, with one value in each of the intervals 0, 1  N , 1 N , 2 N , . . ., 1 − 1 N , 1 .In practice, the Latin Hypercube strategy has been optimized, for instance, with random permutations of the values (coordinate exchange) and with a maximization of the minimum distance between sample points (called the maximin criterion), as implemented in the SAFE toolbox [34].The sampling strategy employed here was Latin Hypercube Sampling [47,48], with N = 3000.
A simple substitution of the Latin Hypercube Sample in the model (in our work, the calculation of the water footprint in the wine viticulture) permitted us to obtain the responses of the model using, for instance, scatter plots.

Analysis of Elementary Effects (EEs)
The first global sensitivity analysis technique applied in our computational experiments was the analysis of the elementary effects (EEs).As pointed out by Morris [49], the main objective of the analysis of elementary effects (also called the Morris method [49]) is to determine whether an effect of such a variable on the response of the model is (i) negligible, (ii) linear and additive, (iii) nonlinear, or (iv) involved in interactions with other variables.
An elementary effect EE i for an input i is defined by the following [49]: where X is the vector representing the variables of the problem and ∆ X is a perturbation in the vector X (only in the position i).Each value of X i may vary in some levels.We aim to compute r EEs, using r trajectories [50].The mean of the effects for r trajectories is then calculated by [50] The standard deviation is represented by [50] In our computations, we computed 500 elementary effects (r = 500).A bootstraping approach (i.e., a resampling strategy) was also applied [51], with N boot = 100.Furthermore, the SAFE package employed the "revised version" of µ, represented by µ * (using the absolute values of the elementary effects) [52].Considering r elementary effects, the number of evaluations of the model is n = r(M + 1) = 500(6 + 1) = 3500, where M is the number of variables considered in the analysis [53].

Fourier Amplitude Sensitivity Testing-FAST
Fourier Amplitude Sensitivity Testing (FAST) and the Variance-Based Sensitivity Analysis (VBSA, or Sobol' method [54]) are variance-based global sensitivity tests [55].The so-called first-order FAST sensitivity index is calculated by where V i is the first-order conditional variance in Y considering the input X i with a value ) and V(y) is the total variance in the model output [56].

Assumptions of This Study
The following assumptions were made in the present study: 1.
The analysis considers only latitude, altitude, the fraction of mulch covering the soil, and the temperatures during three months (October, November, and December) in the water footprint for the wine production; 2.
The water footprint only considers the evapotranspiration portion of the viticulture of wine production; 3.
Temperatures, relative humidities, and wind speeds are considered to be the same for the different latitudes and altitudes (this assumption may be reasonable considering the small size of the region under consideration; on the other hand, new studies can be conducted considering the uncertainties in temperatures and wind speeds, for instance).Additionally, as already evidenced, the range of variation in the maximum temperatures is higher than the real differences in the regions under study.
As mentioned previously, a large number of studies analyzed the effects of temperatures, relative humidities, and wind speeds in the reference evapotranspiration calculated by the Penman-Monteith model.In our study, we focused on the sensitivity of the model with respect to the parameters that can be effectively chosen by the decision-makers, i.e., the site (latitude and altitude) and the use of a mulch covering.

Results
In this section, the results of the computational experiments using the SAFE toolbox applied in the calculation of the water footprint for wine production in the Serra Gaúcha are presented.
The first analysis was conducted with only the sampling strategy and the evaluation of the model, producing a set of scatter plots.Figure 3 presents the scatter plots for the water footprint (WF, in L/bottle of wine) with respect to the covered fraction, the altitude and the latitude of the site, and to the temperatures over three months.As pointed out by Wang and Solomatine [57], the randomness of the distribution in a scatter plot indicates a low sensitivity, while a prominent pattern points to high sensitivity.In this sense, a clear effect of the covered fraction in the water footprint can be observed.We also noticed that the water footprint was directly affected by the altitude (even considering the nonlinear relationship between the altitude and the reference evapotranspiration).On the other hand, the effects of latitude and of the maximum temperatures for October, November, and December seemed to be less important, according to Figure 3. On the other hand, a direct comparison of the effects is not possible only using the scatter plots.On account of that, the global sensitivity analysis tools must be employed.First, we considered the calculation of the EE (Morris indexes).Figure 4 illustrates the convergence of the Morris indexes as a function of the number of model evaluations.The convergence with respect to the screening, to the ranking, and to the sensitivity indexes [53] was obtained with a low number of model evaluations.
Figure 5 presents the means and the standard deviations in the EEs for the water footprint calculation.It can be seen that the water footprint is most sensitive to the covered fraction (highest Morris index value, red circle in the figure), followed by the altitude of the site (yellow circle).In addition, it can be concluded that the water footprint is less sensitive (under the assumptions of this study) to the latitude of the location and to the maximum temperatures of the months under consideration.
An important feature of sensitivity analysis is the identification of interactions between parameters.In this sense, Figure 6 presents the interaction effects between the parameters under study, in order to highlight pairwise interactions.As mentioned by Noacco et al. [42] and Puy et al. [58], the existence of clear patterns (a separation of the colors used in the figure colorbars) in the two-dimensional scatter plots indicates an interaction between the corresponding parameters.We observed more expressive interactions between the covered fraction and the other parameters (indicated , for instance, by the yellow points at the left-upper corner of the first line in the figure).The first diagram in Figure 6 indicates that a combination of a low covered fraction and a high altitude (a second-order interaction) implies a maximum response in the water footprint; on the other hand, a high value for the covered fraction in conjunction with a low altitude value (the right-lower corner of this first diagram) minimizes the water footprint.The other parameters did not show an important correlation, i.e., the points/colors are essentially dispersed in the diagrams.In order to corroborate the results obtained by the EEs, we also employed the FAST sensitivity technique.Figure 7 illustrates the convergence pattern for the FAST indexes for each variable considered in our computations.Clearly, a convergence of the FAST indexes was achieved.In accordance with the results obtained, for instance, in the elementary effects, one can note that the main effects for the covered fraction and for the altitude were more important when compared to the maximum temperatures and to the latitude of the site, since a higher FAST index implies in a higher sensitivity.The accordance of the results obtained by the Morris method and by a Variance-Based Sensitivity Analysis technique (such as FAST) is expected, as pointed out by Campolongo and Cariboni [59].Furthermore, Paleari et al. [60] identified a concordance of the results (in terms of parameter ranking) obtained by the Morris method and e-FAST, although there was no linearity between the parameters obtained by the two methods.Table 1 presents the numerical values for the Morris EE (using µ * ) and the FAST indices.

Discussion
The results obtained indicate that the water footprint in the viticulture in the Serra Gaúcha (under all the hypotheses considered) is more affected by two parameters-the covered fraction and the altitude-when compared to the other quantities evaluated in the study (latitude and maximum temperatures during the months of October, November, and December).This situation is clearly supported by the results of Figures 3, 5 and 7.These findings are compatible, for instance, with the conclusions of Debnath et al. [18], which obtained an important influence of solar radiation in the reference evapotranspiration (in our computational experiments, this can be noted by the high sensitivity to the altitude of the site).This situation indicates that this kind of approach can be extended to other crops/regions, with the necessary adjustments, in order to predict the particular response (in terms of water footprint), for instance, to severe temperature variations.
It must be emphasized that the effect of the temperatures in the reference evapotranspiration model is absolutely nonlinear, as indicated, for instance, by Equations (A2) and (A8).In the same way, the altitude and latitude of the site also appear in a nonlinear form in the reference evapotranspiration model.Considering this fact, it is clear that some parameters which can be manipulated by the decision-maker are able, for instance, to compensate the effect of temperature variations with respect to the water footprint.To some degree, this situation is similar to that described by Arunrat et al. [20], with respect to some crops in Thailand.As discussed previously, Arunrat et al. [20] argued that the greenhouse effect on the water footprint can be compensated for by an increase in crop productivity.In our numerical experiments, it was shown that this effect (a reduction in the water footprint) can also be produced by an increase in the fraction of vegetation cover, even considering an increase in the temperatures.
With regard to the variables related to a possible decision-making process (covered fraction, altitude and latitude of the site), one can observe that the effect of altitude can be mitigated by the extent of coverage, but only up to a certain threshold (see Figure 6).Secondly, as anticipated, the fraction of mulch covering plays a decisive role across the investigated latitude range.Lower altitudes can decrease the water footprint, even at different latitudes.
It should be clear that the evaluation of the water footprint is extremely important from an environmental point of view, mainly considering climate change.For instance, Saraiva et al. [2] recently discussed, in the form of a review in an initial section of their study, the water footprint of wine production for some Mediterranean countries, highlighting the issue of water scarcity in the region.With respect to the Brazilian case, even considering that we are dealing with a tropical region, the Rio Grande do Sul state is subject to irregular rainfall events [61].In this sense, a correct understanding of the main variables that are capable of affecting the evapotranspiration and, by consequence, the water footprint, is a fundamental tool both from an environmental point of view and when making decisions to set up a vineyard.Once again, the main finding of this type of analysis is not the direction of the effect (which may be obvious in some cases; for instance, unsurprisingly, an increase in the maximum temperature for a specific month provokes an increase in the water footprint), but the different magnitudes (which may indicate offsetting effects).These magnitudes are indicated, in a pictorial way, by Figure 3 and are proved by the Morris EE and by FAST indexes, according to Table 1.

Conclusions
In this work, we conducted a global sensitivity analysis study regarding the water footprint in the viticulture in the Brazilian Serra Gaúcha (only considering the contribution of evapotranspiration to the water footprint).The results indicated that the mulch fraction and the altitude of the site are the two most important parameters in the calculation of the water footprint (under the assumptions of this study).Furthermore, the numerical values obtained using different global sensitivity methods (Morris method and FAST) are consistent and convergent.Our results clearly indicate possible compensation between temperature increase (caused by the greenhouse effect) and the use of vegetation cover.In other words, an increase in the dead mulch covering fraction (a land management technique, and, therefore, a decision) is capable of compensating (in this specific situation) for the impact of the elevation in the temperatures with respect to the water consumption.Finally, in Equation (1), G is the heat flux in the soil, calculated by [12] G = 0.07(T month,i+1 − T month,i−1 ), (A15) where T month,i refers to the mean air temperature of the month i.It can be noted that the calculation of the soil heat flux G demands information regarding the mean air temperature of every month of the year.On the other hand, as mentioned by Martins [63], the value of G can be approximated using the simple arithmetic mean between the maximum and minimum temperatures for a month as a good approach.In fact, in many cases, is reasonable to accept that G ≈ 0 [12].

Figure 1 .
Figure 1.Altimetric map of Rio Grande do Sul, Brazil, and localization map of Serra Gaúcha.

Figure 2 .
Figure 2. Maximum temperatures for the months of October, November, and December in the Serra Gaúcha (2013-2023).

Figure 5 .
Figure 5. Elementary effects (EEs) for the water footprint (L/bottle of wine) with r = 500 and N boot = 100.

Figure 6 .
Figure 6.Scatter plots (interaction effects) for the water footprint (L/bottle of wine).

Figure 7 .
Figure 7. Convergence pattern of main effects or first-order effects using FAST.

Table 1 .
Morris EE effects and FAST indices for the water footprint.