Spatial and Temporal Analysis of Water Resources in the Olive-Growing Areas of Extremadura, Southwestern Spain

Abstract

The increasing variability of precipitation, higher temperatures, and recurring droughts in the semi-arid regions due to climate change are leading to increased aridity, resulting in scarcer water resources for crops. The present study aimed to analyse the spatial distribution of climate variables related to water resources in the olive-growing areas throughout Extremadura, southwestern Spain. To perform this task, three climate variables were used: the potential evapotranspiration of the crop, the FAO aridity index, and the annual water requirement. Considering data from 58 weather stations located throughout Extremadura and 17 along boundaries with at least a 30-year length (within the 1991–2021 period), each variable was computed at each station. After calculating some descriptive statistics, a multivariate geostatistical (regression-kriging) algorithm, incorporating secondary information on elevation and latitude, was used to accurately map each climate variable. Later, temporal trends and their magnitude were analysed using the Mann–Kendall test and the Sen’s estimator, respectively. The highest evapotranspiration and water requirements are located in the southern part of the region, which has large areas dedicated to olive cultivation. In the northern part of the region, there is greater spatial variability in evapotranspiration and, consequently, in water requirements for olive groves due to the more rugged topography. Similarly, the olive-growing areas with the highest aridity are also in the south of Extremadura. In most areas of Extremadura, olive cultivation requires appropriate irrigation for optimal productivity. According to evapotranspiration trends, the water requirements will become greater in the future. However, it is not guaranteed that the water supply will be sufficient in olive-growing areas where aridity is higher and water resources are scarce. The results of this study are very important for evaluating water deficit and water resources in vulnerable olive-growing areas throughout Extremadura.

1. Introduction

The reduction in water quantity and quality, along with changing distribution patterns, can have severe impacts on the economy everywhere. Both human- and natural-induced phenomena, whether temporary or permanent, can lead to low water availability and consequently water scarcity, which is associated with droughts, aridity, and desertification.

Annual precipitation is very likely to decrease in most areas of southern Europe [1]. Furthermore, the global mean surface temperature is rising. Europe’s average temperature has risen even faster than the global mean, and the southern part of Europe has dried by as much as 20% in the last century [1]. The water balance has been seriously affected by climate and land use changes.

The main consumers of water resources worldwide are the agricultural ecosystems [2]. In developed countries, water consumption for irrigation accounts for approximately 60% of available water resources, while in developing countries, this can reach up to 90% [3].

The olive tree (Olea europaea L.) is one of the most important perennial crops in the Mediterranean region, with Spain being the country with the largest area dedicated to its cultivation. In the southwest of Spain, in the region of Extremadura, olive groves occupy about 270,000 ha [4]. This crop plays a key role in the economy and landscape of a region with significant socioeconomic inequalities, where the olive tree is crucial for rural development and job creation in disadvantaged areas [5]. The cultivation of olive trees has been fundamental for the socioeconomic resilience of rural areas in Extremadura during recent economic crises [6]. Furthermore, traditional olive groves are considered a sustainable agricultural system, playing an essential role at the social, economic, and environmental levels [7]. However, water scarcity is one of the main limitations for the growth and productivity of this crop in Extremadura, a region characterised by a Mediterranean climate, which presents particular challenges for agriculture, particularly regarding the availability of water for irrigation.

The study of the water requirements of olive trees is crucial for optimising water use and improving irrigation efficiency in this crop. Evapotranspiration is a key parameter for estimating the water requirements of crops. Various studies have estimated olive evapotranspiration in different regions of Spain, with values ranging between 700 and 1100 mm annually, depending on local climatic conditions [8,9]. In the specific case of Extremadura, Orgaz et al. [10] reported olive potential evapotranspiration values ranging between 800 and 900 mm annually, with higher requirements during the summer months. These researchers also noted that under rainfed conditions, the average annual precipitation in the region is not sufficient to meet the water requirements, necessitating supplemental irrigation to maintain the productivity of olive groves.

Other climatic indices can be used to characterise water availability in a region and its suitability for olive cultivation. For instance, the FAO aridity index, which considers the ratio between precipitation and evapotranspiration, has been applied in studies throughout Extremadura [11]. The results of these studies have shown that Extremadura has semi-arid climate conditions, with aridity indices indicating a marked water scarcity for much of the year. Moral et al. [11] indicated that in Extremadura, aridity during spring and summer was high (semi-arid) and very high (arid) across almost the entire region. These seasons are critical for the olive tree’s vegetative activity, which limits the development and productivity of olive groves under rainfed conditions. Ullah et al. [12] recommended using aridity in indices at regional and local scales to obtain more robust results and formulate better planning, adaptation, and mitigation measures related to water resources.

Another widely used index quantifies the crop’s water requirements as the difference between evapotranspiration and precipitation [13,14]. This difference is crucial for determining the amount of supplementary irrigation needed to maintain the crop in optimal conditions.

In addition to climatic characterisation, the analysis of temporal trends can provide valuable information about the potential effects of climate change on water availability for olive cultivation. Several studies have reported decreasing precipitation trends and increasing evapotranspiration trends in various olive-growing regions of Spain, including Extremadura [15,16]. Net irrigation requirements are expected to increase by 18.5% (70 ± 28 mm per season), reaching up to 140 mm in southern Spain and some areas of Algeria and Morocco [17]. Although the olive tree is a species highly adapted to water scarcity, it performs better with high precipitation or irrigation [18,19]. These changes in climatic conditions could have a significant impact on the sustainability of olive cultivation in the future, as net water requirements could exceed potential evapotranspiration due to a decrease in precipitation and an even greater impact from droughts and dry periods [17].

Some studies have analysed how crop water requirements will be affected by the impacts of climate change, making it important to have an assessment of water resources based on recent climatic data. For instance, Yassen et al. [20] found that significant changes in the spatial distribution of reference evapotranspiration have begun in the last 35 years in Egypt, and Zhao et al. [21] studied the spatial and temporal variability of evapotranspiration considering different climatic scenarios in China. Similar studies have also been reported elsewhere [22,23].

The objective of this research is to determine the spatial distribution of climate variables related to water resources in the olive-growing areas throughout Extremadura using geostatistical techniques and to analyse their temporal trends. The goal is to develop future water management strategies to improve irrigation efficiency and the sustainability of this crop in the region.

2. Materials and Methods

2.1. Study Area

Extremadura is a region located in the southwest of Spain, bordering Portugal, covering an area of 41,634 km², one of the largest regions in Europe. It ranges from 37°57′ to 40°29′ N and from 4°39′ to 7°33′ W. This region has two provinces: Cáceres in the north and Badajoz in the south (Figure 1). Extremadura has vast agricultural and forest areas and is recognised as an important ecological enclave in Europe. The mean elevation in Extremadura is 425 m a.s.l., ranging from 116 m in the plain of the Guadiana river, near the border between Spain and Portugal, to 2404 m in the north of the region, at the highest elevation in the Sistema Central mountains.

Olive cultivation in Extremadura covers an area of 271,404 ha, representing 10.30% of the total olive-growing area in Spain. It is the third autonomous community in terms of the area dedicated to olive cultivation and olive oil production, and the second in terms of table olives [4]. The region is divided into 12 olive-growing subregions (OGSs), with 6 in each province (Figure 1). Olive cultivation is mainly located at elevations ranging between 200 and 600 m a.s.l.

The climate in Extremadura is Mediterranean, influenced by the nearby Atlantic ocean and its inland location. The mean annual temperature is 16.5 °C, but summers are extremely dry and hot, frequently reaching temperatures over 40 °C. Winters can be mild and rainy when low-pressure systems dominate or cold and dry when strong anticyclones from northern Europe prevail [24]. During spring and autumn, high climatic variability is usual. Although the mean annual precipitation is 668 mm, it averages less than 600 mm across many zones of the region. Moreover, its interannual variability is very important. Extremadura experiences a dry season, from June to September, and a wet season, from October to May, with 80% of the annual precipitation. The interannual variability of rainfall can be a constraint to olive cultivation, since 50% of the olive-growing area receives less rainfall than is considered adequate for optimum production [25]. The predominant olive cultivation system in the region is extensive and rainfed. However, in recent years, there has been a significant increase in the area dedicated to irrigated and super-intensive olive cultivation, resulting in a notable rise in average yield.

2.2. Database, Data Treatment, and Climate Variables

Daily meteorological observations were obtained at 75 georeferenced weather stations of the Meteorology Agency of the Spanish Government, AEMET [26], located in Extremadura and along boundaries (Figure 1). The initial database contained the UTM geographic coordinates, daily temperature (mean, maximum, and minimum), and rainfall over the 30-year study period (1991–2021). These data were analysed and the mean monthly, seasonal, and annual values of precipitation, air temperature, and evapotranspiration were computed, forming the database for subsequent calculations regarding other climate variables. All daily data sets were homogenised and verified according to the quality controls recommended by the World Meteorological Organisation [27] and the Royal Netherlands Meteorological Institute [28]. Missing data were filled in, and the homogenisation process was applied using the R package CLIMATOL [29]. The selected weather stations provide adequate spatial coverage and are representative of the climate’s behaviour and variability in the region. Additionally, using some weather stations along the boundaries of Extremadura can improve estimates in border areas.

Based on the climatic database, three indices were used in this research to calculate the annual averages for the period studied: potential evapotranspiration (PET), FAO aridity index (I_F), and annual water requirement (AWR).

As evapotranspiration is a simultaneous process of soil water evaporation and plant transpiration, Hargreaves and Samani [30] considered this principle to develop estimates of reference evapotranspiration with only air temperature data. The evapotranspiration that occurs in a disease-free crop, growing under optimal soil conditions, including water and fertility, is the potential crop evapotranspiration (ETc). ETc refers to the ideal climatic water demand, achieved by keeping the soil close to the field capacity, according to the prevailing climatic conditions. The ETc for olive cultivation was calculated according to the manual produced by Perez-Rodríguez et al. [31]. In this study, PET is calculated using the Hargreaves method [30]:

(1)$\mathrm{P}\mathrm{E}\mathrm{T}=0.0023{\mathrm{R}}_{\mathrm{a}}\left({\displaystyle \frac{{\mathrm{T}}_{\mathrm{m}\mathrm{a}\mathrm{x}}+{\mathrm{T}}_{\mathrm{m}\mathrm{i}\mathrm{n}}}{2}}+17.8\right){{(\mathrm{T}}_{\mathrm{m}\mathrm{a}\mathrm{x}}+{\mathrm{T}}_{\mathrm{m}\mathrm{i}\mathrm{n}})}^{0.5}$

where PET is the daily evapotranspiration, R_a is the extraterrestrial radiation, and T_max and T_min are the maximum and minimum daily temperature, respectively.

I_F is the ratio between the mean annual precipitation, P_a, and the mean potential evapotranspiration, PET_a [32,33]:

(2)${{\mathrm{I}}_{\mathrm{F}}={\displaystyle \frac{{\mathrm{P}}_{\mathrm{a}}}{\mathrm{P}\mathrm{E}\mathrm{T}}}}_{\mathrm{a}}$

The AWR of a crop is the difference between the PET of the crop and the effective precipitation (the portion of rainfall used by the soil), calculated according to the USDA method [34].

2.3. Interpolation Technique and Trend Analysis

While numerous interpolation methods exist for estimating values at unsampled locations, geostatistical algorithms generally outperform deterministic techniques [35,36]. Geostatistics offers a wide range of methods, collectively known as kriging, which are utilised to estimate values at unsampled locations. Various types of kriging involve estimating the value of the studied variable at a specific location based on neighbouring sampling points and a variogram model that considers the distance and degree of variation among all sampling locations. This approach aims to minimise the variance of the estimation error [37].

Univariate geostatistical algorithms, especially ordinary kriging, are commonly employed for estimating climatic variables at unsampled points [38,39]. However, incorporating auxiliary information along with primary data can enhance estimation accuracy through multivariate extensions of kriging methods [35,36]. In this study, the regression kriging algorithm was used. This method is suitable when comprehensive secondary information is available, meaning auxiliary data are accessible throughout the study area [36]. Regression kriging was chosen because it has been shown to produce fewer errors when used to estimate climate variables in Extremadura [35].

In regression kriging, predictions are conducted in two distinct stages: first for the trend component and then for the residuals. These two predictions are subsequently combined to derive the final estimated values. Therefore, the estimation of any climate variable, ${\mathrm{Z}}_{\mathrm{R}\mathrm{K}}^{\mathrm{*}}(\mathrm{x})$, at a new unsampled point, x, is performed as follows:

(3)${\mathrm{Z}}_{\mathrm{R}\mathrm{K}}^{\mathrm{*}}(\mathrm{x})=\mathrm{m}(\mathrm{x})+\mathrm{r}(\mathrm{x})$

Here, m(x), represents the trend fitted using linear regression analysis and r(x) represents the residuals estimated using ordinary kriging algorithm.

In this study, two predictors were used, elevation (h) and latitude (l). Thus,

m(x) = a + b h(x) + c l(x) (4)

and

(5)${\mathrm{Z}}_{\mathrm{R}\mathrm{K}}^{\mathrm{*}}\left(\mathrm{x}\right)=\mathrm{a}+\mathrm{b}\hspace{0.33em}\mathrm{h}\left(\mathrm{x}\right)+\mathrm{c}\mathrm{l}(\mathrm{x})+\sum _{\mathrm{i}=1}^{\mathrm{n}}{\mathrm{w}}_{\mathrm{i}}(\mathrm{x})\cdot \mathrm{r}({\mathrm{x}}_{\mathrm{i}})$

The difference between the value of the soil property, Z(x_i), and the estimate provided by the trend, m(x_i), is the residual at each sampling point, r(x_i). The weights, w_i(x), are computed by solving the ordinary kriging system of the regression residuals, r(x_i), for the n sample points.

A digital elevation model for Extremadura was used to extract elevation in raster format (1000 m × 1000 m resolution). The latitude for the entire region was also utilised in raster format, with the same resolution as the elevation map. Therefore, using point data from weather stations at sampling locations, estimates can be made for any other unsampled location. This approach enables the creation of a continuous surface that covers the entire region, allowing for the determination of values for each climate variable (PET, I_F and AWR) at every pixel of the smallest resolution unit. Digital models for each climate variable were generated in raster format with a resolution of 1000 m × 1000 m. All operations, including the spatial representation and visualisation of the climate variables, were performed using the GIS software ArcGIS v. 10.5. The geostatistical analyses were conducted utilising the Geostatistical Analyst extension within ArcGIS.

The Mann–Kendall test was used to detect temporal trends in the three time series of climate variables, as recommended by the World Meteorological Organisation [27]. After computing the statistic of the Mann–Kendall test, increasing trends are indicated by positive values, while decreasing trends are indicated by negative values. In this study, the 5% significance level was used, meaning the null hypothesis of no trend was rejected when the absolute value of the statistic was higher than 1.96.

With the aim of identifying the gradients and their directions, the Sen’s non-parametric method was used [40]. Thus, the MAKESENS 1.0 was utilised. It is a computer model introduced by Salmi et al. [41], which was built using Microsoft Excel 97 and macros coded with the Microsoft Visual Basic language.

3. Results and Discussion

3.1. Descriptive Analysis of Data and Spatial Distribution of the Climate Variables

Some descriptive statistics were used to characterise the data distribution (Table 1). The mean and median values for THE PET were similar. Moreover, the skewness value was low; in consequence, all these values are indicative of the PET data coming from a normal distribution. However, more differences between the median and the mean values were apparent for the AWR and I_F. Additionally, the skewness values were higher. These statistics suggest that the AWR and I_F data fitted a lognormal distribution.

The distributions of the PET and AWR were shifted to higher values, as indicated by the median values being higher than the mean values. Conversely, the distribution of I_F was shifted to lower values, with the mean value being higher than the median value. This is coherent because higher PET values lead to increased water requirements and, consequently, greater aridity, meaning lower I_F values should predominate. Furthermore, the wide difference between the minimum and maximum values and the high coefficients of variation for all climate variables denote significant climatic spatial variability in Extremadura.

Despite the previous description of the three climate variables in Extremadura, based on point values at the weather stations, it is evident that spatial variability can exist within all OGSs. Therefore, a more accurate analysis is necessary to properly understand the patterns of each variable considered in this study. Consequently, the spatial representation of each climate variable is an essential tool for precisely visualising the climatic characteristics in Extremadura.

The spatial patterns of the PET, AWR, and I_F in Extremadura are shown in Figure 2. Table 2, Table 3 and Table 4 show some statistics for each OGS, considering the three climatic variables.

If the mean PET is considered (Table 2), the highest values were found in the southern part of Extremadura, in six OGSs in the province of Badajoz (Alburquerque, Vegas del Guadiana, Tierra de Barros, La Siberia, La Serena, Jerez-Llerena) and in the province of Cáceres (Montánchez). These mean PET values range from 825.15 mm in Alburquerque to 854.14 mm in Vegas del Guadiana. Conversely, the lowest mean PET values were found in the northern part of the region, where elevation is higher due to the existence of some mountain systems. In these OGSs, the mean PET values range from 753.33 mm in LaVera-Jerte-Ambroz to 814.23 mm in the centre of Cáceres. Moreover, variability is higher in the northern OGS, as expected, due to significant differences in elevation.

When the AWR was considered, only one OGS (La Vera-Jerte-Ambroz) showed a positive value, indicating that irrigation is not necessary there. In all other OGSs, there is an annual water deficit, ranging from 60.64 mm to 392.16 mm, which is more pronounced in the southern OGS of Extremadura (Table 3). As with the PET, variability was also higher for the AWR in the northern OGS due to the more significant topographical contrast.

According to the IF, there are three categories in Extremadura. Dry sub-humid (0.50 ≤ IF < 0.65) characteristics, where the IF ranges from 0.55 to 0.64, are predominant (Table 4). In these OGSs, precipitation is well below the values of the PET and, consequently, the water requirements through irrigation are more pronounced. Logically, these OGSs correspond to those showing higher AWR values. Conversely, the OGSs in northern Extremadura, in the province of Cáceres, have the highest mean IF values, ranging between 0.85 and 1.23, falling into the more humid category. In these OGSs, there is a lower water requirement through irrigation, with the extreme case of La Vera-Jerte-Ambroz, where precipitation even exceeds the PET demand. Three other OGSs, located in the central part of the region, fall into a medium IF category, with values between 0.67 and 0.72, but they also require irrigation to meet the PET demand.

The spatial distribution of the PET, AWR, and I_F in the Extremaduran OGSs was properly provided by maps of kriged estimates, and, additionally, these surfaces were very accurate. When the ratio between the root mean square error and the mean value of each climate variable was computed, values lower than 3.5% were obtained for all variables. This fact denotes that all maps very closely characterised the real spatial patterns of each climate variable.

The climate variations observed in this research according to the three indices suggest that olive cultivation management, particularly regarding water management, needs to be approached differently in each OGS. This spatial variability could become even more significant in the near future, as a decrease in water resources of more than 50% is expected in approximately 90% of Extremadura [42], along with a significant increase in aridity conditions [11]. The productivity of olive cultivation will be significantly reduced under these new climatic conditions, to the point where it may no longer be viable to grow olives under rainfed conditions in many OGSs where it is currently feasible. Additionally, the optimal areas for olive cultivation could shift to OGSs with higher elevation and a greater availability of water resources.

3.2. Trend Analysis

After computing the statistic of the Mann–Kendall test and the magnitude of trends approximated by the Sen’s estimator, the results showed that only the PET at some locations exhibited significant increasing trends.

Figure 3 shows the weather stations in Extremadura and around the region where trends were significant at the 0.05 level. Most of these locations in Extremadura are in the western part of the region, with all trends being positive, indicating that higher PET values are expected. Only one location in the north and at a site with a more humid climate, according to the I_F, has a significant trend. In contrast, in the south locations, where more arid conditions prevail, the PET is the highest. Consequently, with an increase in the PET in these areas, the water availability decreases, leading to a drier climate, which impacts the water requirements of olive trees.

With respect to the increasing trend in the PET at those locations in Extremadura in which the trends were significant, it ranges from 0.96 mm to 2.47 mm. In particular, in the two northern locations, their values were 1.01 mm and 1.70 mm, and in the four located more to the south, their values were higher than 2 mm. This suggests that the southwestern area of Extremadura, and the existing OGSs there, can be more seriously affected by the increasing PET. Figure 3, where the spatial pattern of the magnitude of this trend is visualised, shows this fact more easily. Using the ordinary kriging algorithm, the kriged map shows how the north and east of Extremadura constitute the zones of the region in which the trends have a lower increasing rate, lower than 0.5 mm. Even in the easternmost area of the region, there is a large zone with negative trends, meaning that the PET decreases. In contrast, the western area has a higher increasing rate, with values higher than 1 mm.

As drier regions are more sensitive to the availability of water resources [43], changes in the PET can alter the intensity of desertification. This can more seriously affect the OGSs in the southwestern of Extremadura, where the PET is higher and the AWR is lower. In these OGSs, olive groves require a large amount of water, which will be larger in the near future, but may not be available. Moreover, aridity in Extremadura, particularly in the south of the region, will progressively increase during this century [44]. This makes the region more vulnerable to climate change, shifting the OGSs to drier conditions, more intensively in the south. These results are consistent with those obtained by Lozano-Parra and Sánchez-Martín [42], who indicated that water resources in Extremadura could decrease by up to 55% in agricultural areas.

Several previous studies have shown similar trends in the PET, indicating that the causes of the observed trends are the increases in both maximum and minimum temperatures [22], the increase in net radiation [45], and even the combined contribution of wind, temperatures, and deficits in vapour pressure and relative humidity [23].

4. Conclusions

The availability of water resources for olive cultivation, both under rainfed and irrigated conditions, is essential for achieving good productivity and, consequently, optimal profitability. Therefore, high-resolution knowledge of the spatial distribution of the main climatic variables that influence water availability is of great importance for evaluating this resource at both the regional and local scales.

The present study analyses these spatial variations in Extremadura using the PET, IF, and AWR as climatic indices, considering the OGSs throughout the region. Moreover, the temporal trends were also studied to analyse the expected evolution. Although significant increasing trends were observed only in the PET at some locations, the results of this study are crucial for evaluating water resources. This information can help predict practical measures to control more vulnerable areas, where olive cultivation could be seriously affected in the near future due to compromised water supply.

In most areas of Extremadura, the AWR is negative, indicating that olive cultivation requires appropriate irrigation for optimal productivity. According to PET trends, the AWR will become more negative in the future, necessitating even more water. However, it is not guaranteed that the water supply will be sufficient in areas where the climate is characterised by aridity and where water resources are highly scarce and under intense competition, primarily for non-agricultural uses.

The majority of olive cultivation in Extremadura is rainfed, and it is expected to face greater stress in the future, becoming unfeasible in many areas of the region. This situation will have significant consequences, including land abandonment and landscape deterioration, leading to a loss of soil and water conservation measures and increased erosion. Additionally, the social consequences of abandoning traditional olive groves could be enormous.

Olive evapotranspiration and irrigation requirements strongly depend on the type of cultivation (traditional, intensive, or super-intensive); that is, tree density and the fraction of ground cover. Consequently, crop water requirements could increase even more if super-intensive olive groves continue to expand. Due to potentially high water requirements and limited water availability, as it is currently developing, regulated deficit irrigation will likely become a common practise in most areas. However, achieving water savings and high yields requires careful irrigation scheduling.

Author Contributions

Conceptualization, F.J.M. and L.L.P.; methodology, F.J.M., A.G.-M. and L.L.P.; software, F.J.M. and F.J.R.; validation, F.J.M. and L.L.P.; formal analysis, F.J.M., L.L.P. and A.G.-M.; investigation, F.J.M., L.L.P., A.G.-M., F.J.R. and F.H.; resources, L.L.P. and F.H.; data curation, F.J.M., F.J.R. and L.L.P.; writing—original draft preparation, F.J.M.; writing—review and editing, F.J.M.; visualisation, F.J.R., A.G.-M. and L.L.P.; supervision, F.J.M.; project administration F.H.; funding acquisition, F.J.M., A.G.-M. and F.H. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Initial data used in this research can be found at https://www.aemet.es/es/serviciosclimaticos/datosclimatologicos, accessed on 16 February 2024.

Acknowledgments

This study has been possible thanks to the collaboration of the Agencia Estatal de Meteorología (AEMET). The authors are also very grateful to three anonymous reviewers for providing constructive comments that contributed to improving the final version of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Digital elevation model and olive growing subregions of Extremadura.

Figure 2. Spatial distribution of the average annual values for the potential evapotranspiration (PET), the FAO aridity index (IF), and the annual water requirement (AWR) in Extremadura considering the baseline climate, 1991–2021. Olive-growing subregions and olive areas (in green) are also indicated.

Figure 3. Spatial distribution of the magnitude of trend of the potential evapotranspiration in Extremadura. Black points show the locations in which significant trends at the 5% level were found.

Table 1

Descriptive statistics of the climate variables from 75 weather stations in Extremadura and along its boundaries.

Variable	Mean	Median	Standard Deviation	Minimum	Maximum	Coefficient of Variation (%)	Skewness
Potential evapotranspiration (PET)	821.62	828.89	56.32	653.75	902.06	6.85	−0.89
FAO aridity index (IF)	−0.76	0.66	0.30	0.47	1.88	39.47	2.01
Annual water requirement (AWR)	−228.01	−297.73	225.71	−475.02	519.06	98.99	1.71

Table 2

Descriptive statistics of the potential evapotranspiration in each Extremaduran olive-growing subregion.

Natural Region	Potential Evapotranspiration (PET)
	Minimum	Maximum	Mean	SD	CV (%)
Gata-Hurdes	624.69	826.17	768.84	39.01	5.07
La Vera-Jerte-Ambroz	497.62	856.37	753.33	73.32	9.73
Ibores	694.52	845.89	793.59	31.39	3.95
Logrosán-Guadalupe	660.92	868.87	809.80	34.06	4.21
Montánchez	754.91	865.73	829.08	19.68	2.37
Resto de Cáceres	731.44	862.48	814.23	13.19	1.62
Alburquerque	784.89	873.53	825.15	16.53	2.00
La Siberia	778.48	882.11	847.28	18.03	2.13
La Serena	784.34	871.36	840.90	12.61	1.49
Vegas del Guadiana	799.09	889.56	854.14	14.59	1.71
Tierra de Barros	775.24	870.14	834.57	15.33	1.84
Jerez-Llerena	733.83	886.29	832.74	21.29	2.56

Table 3

Descriptive statistics of the annual water requirements in each Extremaduran olive-growing subregion.

Natural Region	Annual Water Requirement (AWR)
	Minimum	Maximum	Mean	SD	CV (%)
Gata-Hurdes	−247.89	432.23	−60.64	112.71	185.86
La Vera-Jerte-Ambroz	−217.03	984.74	135.36	286.94	211.98
Ibores	−298.11	279.45	−145.49	101.39	69.68
Logrosán-Guadalupe	−364.63	396.48	−147.96	130.88	88.46
Montánchez	−405.37	−11.21	−258.79	79.27	30.63
Resto de Cáceres	−355.48	137.70	−221.89	53.82	24.26
Alburquerque	−389.61	−110.25	−276.81	58.31	21.06
La Siberia	−437.22	−86.48	−301.39	67.54	22.41
La Serena	−477.12	−216.27	−381.80	44.12	11.56
Vegas del Guadiana	−468.58	−219.28	−392.16	40.11	10.23
Tierra de Barros	−469.00	−147.15	−359.74	51.01	14.18
Jerez-Llerena	−479.71	17.64	−317.08	78.10	24.63

Table 4

Descriptive statistics of the FAO aridity index in each Extremaduran olive-growing subregion.

Natural Region	FAO Aridity Index (IF)
	Minimum	Maximum	Mean	SD	CV (%)
Gata-Hurdes	0.70	1.61	0.95	0.15	15.79
La Vera-Jerte-Ambroz	0.75	2.32	1.23	0.39	31.71
Ibores	0.64	1.42	0.85	0.13	15.29
Logrosán-Guadalupe	0.58	1.58	0.86	0.17	19.78
Montánchez	0.54	1.05	0.72	0.10	13.89
Resto de Cáceres	0.57	1.21	0.64	0.07	10.94
Alburquerque	0.55	0.91	0.69	0.07	10.14
La Siberia	0.49	0.94	0.67	0.09	13.43
La Serena	0.44	0.78	0.57	0.06	10.53
Vegas del Guadiana	0.46	0.77	0.55	0.05	9.09
Tierra de Barros	0.46	0.86	0.59	0.06	10.17
Jerez-Llerena	0.45	1.08	0.65	0.09	13.85

References

1. IPCC. Climate Change 2023: Synthesis Report. Contribution of Working Groups I, II and III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team Lee, H.; Romero, J. IPCC: Geneva, Switzerland, 2023.

2. Pedro-Monzonís, M.; Solera, A.; Ferrer, J.; Estrela, T.; Paredes-Arquiola, J. A review of water scarcity and drought indexes in water resources planning and management. J. Hydrol.; 2015; 527, pp. 482-493. [DOI: https://dx.doi.org/10.1016/j.jhydrol.2015.05.003]

3. Adeyemi, O.; Grove, I.; Peets, S.; Norton, T. Advanced monitoring and management systems for improving sustainability in precision irrigation. Sustainability; 2017; 9, 353. [DOI: https://dx.doi.org/10.3390/su9030353]

4. Ministerio de Agricultura, Pesca y Alimentación. Gobierno de España. Anuario de Estadística. Available online: https://www.mapa.gob.es/es/estadistica/temas/publicaciones/anuario-de-estadistica/2021/default.aspx?parte=3&capitulo=07&grupo=12 (accessed on 19 April 2024).

5. Beaufoy, G. The Environmental Impact of Olive Oil Production in the European Union: Practical Options for Improving the Environmental Impact: Final Report for the Environment Directorate-General of the European Commission. 2000; Available online: https://www.efncp.org/download/oliveoil.pdf (accessed on 5 May 2024).

6. Sánchez-Zamora, P.; Fernández-Portillo, L.A.; Gallardo-Cobos, R.; Céspedes-Lorente, A. Rural areas face the economic crisis: Analyzing the determinants of successful territorial dynamics. J. Rural Stud.; 2014; 35, pp. 11-25. [DOI: https://dx.doi.org/10.1016/j.jrurstud.2014.03.007]

7. Fernández-Escobar, R.; de la Rosa, R.; León, L.; Gómez, J.A.; Testi, F.; Orgaz, M.; Gil-Ribes, J.A.; Quesada-Moraga, E.; Trapero, A. Evolution and sustainability of the olive production systems. Present and Future of the Mediterranean Olive Sector; Arcas, N.; Arroyo López, F.N.; Caballero, J.; D’Andria, R.; Fernández, M.; Fernandez Escobar, R.; Garrido, A.; López-Miranda, J.; Msallem, M.; Parras, M. et al. Options Méditerranéennes: Série A. Séminaires Méditerranéens CIHEAM: Zaragoza, Spain, 2013; Volume 106, pp. 11-41.

8. Villalobos, F.J.; Orgaz, F.; Testi, L.; Fereres, E. Measurements and modeling of evapotranspiration of olive (Olea europaea L.) orchards. Eur. J. Agron.; 2000; 13, pp. 155-163. [DOI: https://dx.doi.org/10.1016/S1161-0301(00)00071-X]

9. Fernández, J.E.; Diaz-Espejo, A.; Romero, R.; Hernandez-Santana, V.; García, J.M.; Padilla-Díaz, C.M.; Cuevas, M.V. Precision irrigation in olive (Olea europaea L.) tree orchards. Water Scarcity and Sustainable Agriculture in Semiarid Environment; Academic Press: Cambridge, MA, USA, 2018; pp. 179-217.

10. Orgaz, F.; Testi, L.; Villalobos, F.J.; Fereres, E. Water requirements of olive orchards-II: Determination of crop coefficients for irrigation scheduling. Irrig. Sci.; 2006; 24, pp. 77-84. [DOI: https://dx.doi.org/10.1007/s00271-005-0012-x]

11. Moral, F.J.; Rebollo, F.J.; Paniagua, L.L.; García-Martín, A.; Honorio, F. Spatial distribution and comparison of aridity indices in Extremadura, southwestern Spain. Theor. Appl. Climatol.; 2016; 126, pp. 801-814. [DOI: https://dx.doi.org/10.1007/s00704-015-1615-7]

12. Ullah, S.; You, Q.; Sachindra, D.A.; Nowosad, M.; Ullah, W.; Bhatti, A.S.; Jin, Z.; Ali, A. Spatiotemporal Changes in Global Aridity in Terms of Multiple Aridity Indices: An Assessment Based on the CRU Data. Atmos. Res.; 2022; 268, 105998. [DOI: https://dx.doi.org/10.1016/j.atmosres.2021.105998]

13. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998.

14. Todorovic, M.; Albrizio, R.; Zivotic, L.; Abi Saab, M.T.; Stöckle, C.; Steduto, P. Assessment of AquaCrop, CropSyst, and WOFOST models in the simulation of sunflower growth under different water regimes. Agron. J.; 2009; 101, pp. 509-521. [DOI: https://dx.doi.org/10.2134/agronj2008.0166s]

15. Picón-Toro, J.; González-Dugo, V.; Uriarte, D.; Mancha, L.A.; Testi, L. Effects of canopy size and water stress over the crop coefficient of a “Tempranillo” vineyard in south-western Spain. Irrig. Sci.; 2012; 30, pp. 419-432. [DOI: https://dx.doi.org/10.1007/s00271-012-0351-3]

16. Moriondo, M.; Ferrise, R.; Trombi, G.; Brilli, L.; Dibari, C.; Bindi, M. Modelling olive trees and grapevines in a changing climate. Environ. Modell. Softw.; 2015; 72, pp. 387-401. [DOI: https://dx.doi.org/10.1016/j.envsoft.2014.12.016]

17. Tanasijevic, L.; Todorovic, M.; Pereira, L.S.; Pizzigalli, C.; Lionello, P. Impacts of climate change on olive crop evapotranspiration and irrigation requirements in the Mediterranean region. Agric. Water Manag.; 2014; 144, pp. 54-68. [DOI: https://dx.doi.org/10.1016/j.agwat.2014.05.019]

18. Iniesta, F.; Testi, L.; Orgaz, F.; Villalobos, F.J. The effects of regulated and continuous deficit irrigation on the water use, growth and yield of olive trees. Eur. J. Agron.; 2009; 30, pp. 258-265. [DOI: https://dx.doi.org/10.1016/j.eja.2008.12.004]

19. Palese, A.M.; Nuzzo, V.; Favati, F.; Pietrafesa, A.; Celano, G.; Xiloyannis, C. Effects of water deficit on the vegetative response, yield and oil quality of olive trees (Olea europaea L.; cv Coratina) grown under intensive cultivation. Sci. Hortic.; 2010; 125, pp. 222-229. [DOI: https://dx.doi.org/10.1016/j.scienta.2010.03.025]

20. Yassen, A.N.; Nama, W.H.; Honge, E.M. Impact of climate change on reference evapotranspiration in Egypt. Catena; 2020; 194, 104711. [DOI: https://dx.doi.org/10.1016/j.catena.2020.104711]

21. Zhao, J.; Xia, H.; Yue, Q.; Wang, Z. Spatiotemporal variation in reference evapotranspiration and its contributing climatic factors in China under future scenarios. Int. J. Climatol.; 2020; 40, pp. 3813-3831. [DOI: https://dx.doi.org/10.1002/joc.6429]

22. Moses, O.; Hambira, W.L. Effects of climate change on evapotranspiration over the Okavango Delta water resources. Phys. Chem. Earth.; 2018; 105, pp. 98-103. [DOI: https://dx.doi.org/10.1016/j.pce.2018.03.011]

23. Gaertner, B.A.; Zegre, N.; Warner, T.; Fernandez, R.; He, Y.; Merriam, E.R. Climate, forest growing season, and evapotranspiration changes in the central Appalachian Mountains, USA. Sci. Total Environ.; 2019; 650, pp. 1371-1381. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2018.09.129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30308824]

24. Hurrell, J.W.; Van Loon, H. Decadal variations in climate associated with the North Atlantic Oscillation. Clim. Chang.; 1997; 36, pp. 301-326. [DOI: https://dx.doi.org/10.1023/A:1005314315270]

25. Honorio, F.; Aguirado, C.; Paniagua, L.L.; García-Martín, A.; Rebollo, L.; Rebollo, F.J. Exploring the Climate and Topography of Olive Orchards in Extremadura, Southwestern Spain. Land; 2024; 13, 495. [DOI: https://dx.doi.org/10.3390/land13040495]

26. AEMET, Agencia Estatal de Meteorología. Available online: https://www.aemet.es/es/portada (accessed on 16 February 2024).

27. WMO (World Meteorological Organization). Guide to Hydrological Practices. Volume II—Management of Water Resources and Application of Hydrological Practices; 6th ed. WMO-No. 168 World Meteorological Organization: Geneva, Switzerland, 2009.

28. KNMI Royal Netherlands Meteorological Institute. ECAD European Climate Assessment & Dataset. Algorithm Theoretical Basis Document; Royal Netherlands Meteorological Institute KNMI: Utrecht, The Netherlands, 2013.

29. Guijarro, J. Homogenization of Climatic Series with Climatol; CLIMATOL, AEMET: Balearic Islands, Spain, 2018; Available online: https://climatol.eu/ (accessed on 8 April 2024).

30. Hargreaves, G.H.; Samani, Z.A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric.; 1985; 1, pp. 96-99. [DOI: https://dx.doi.org/10.13031/2013.26773]

31. Pérez-Rodriguez, J.M.; Parras, J.; Lara, E.; Prieto, M.H. Manual Práctico de Riego Olivar de Almazara; Centro de Investigaciones Científicas y Tecnológicas de Extremadura, CICYTEX: Badajoz, Spain, 2014.

32. Paniagua, L.L.; García-Martín, A.; Moral, F.J.; Rebollo, F.J. Aridity in the Iberian Peninsula (1960–2017): distribution, tendencies, and changes. Theor. Appl. Climatol.; 2019; 138, pp. 811-830. [DOI: https://dx.doi.org/10.1007/s00704-019-02866-0]

33. Nastos, P.T.; Politi, N.; Kapsomenakis, J. Spatial and temporal variability of the aridity index in Greece. Atmos. Res.; 2013; 119, pp. 140-152. [DOI: https://dx.doi.org/10.1016/j.atmosres.2011.06.017]

34. United States Department of Agriculture, USDA. National Engineering Handbook. Irrigation Guide Available online: https://www.nrcs.usda.gov/sites/default/files/2022-11/WA-Irrigation-Guide_4.pdf (accessed on 17 June 2024).

35. Moral, F.J. Comparison of different geostatistical approaches to map climate variables: Application to precipitation. Int. J. Climatol.; 2010; 30, pp. 620-631. [DOI: https://dx.doi.org/10.1002/joc.1913]

36. Goovaerts, P. Using elevation to aid the geostatistical mapping of rainfall erosivity. Catena; 1999; 34, pp. 227-242. [DOI: https://dx.doi.org/10.1016/S0341-8162(98)00116-7]

37. Goovaerts, P. Geostatistics for Natural Resources Evaluation; Oxford University Press: New York, NY, USA, 1997.

38. Hrnjak, I.; Lukić, T.; Gavrilov, M.B.; Marković, S.B.; Unkašević, M.; Tošić, I. Aridity in Vojvodina, Serbia. Theor. Appl. Climatol.; 2014; 115, pp. 323-332. [DOI: https://dx.doi.org/10.1007/s00704-013-0893-1]

39. Croitoru, А.Е.; Piticar, А.; Imbroane, А.М.; Burada, D.C. Spatiotemporal distribution of aridity indices based on temperature and precipitation in the extra-Carpathian regions of Romania. Theor. Appl. Climatol.; 2013; 112, pp. 597-607. [DOI: https://dx.doi.org/10.1007/s00704-012-0755-2]

40. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc.; 1968; 63, pp. 1379-1389. [DOI: https://dx.doi.org/10.1080/01621459.1968.10480934]

41. Salmi, T.; Määttä, A.; Anttila, P.; Ruoho-Airola, T.; Amnell, T. Publications on Air Quality no. 31. Detecting Trends of Annual Values of Atmospheric Pollutants by the Mann-Kendall Test and Sen’s Slope Estimates–The Excel Template Application MAKESENS; Finnish Meteorological Institute: Helsinki, Finland, 2002.

42. Lozano-Parra, J.; Sánchez-Martín, J.M. Exploring the Impact of Climate Change on Water Resources for Vegetation Covers in Extremadura (Spain). Water; 2024; 16, 1418. [DOI: https://dx.doi.org/10.3390/w16101418]

43. Zhang, Q.; Xu, C.-Y.; Tao, H.; Jiang, T.; Chen, Y.D. Climate changes and their impacts on water resources in the arid regions: A case study of the Tarim River basin, China. Stoch. Environ. Res. Risk Assess.; 2010; 24, pp. 349-358. [DOI: https://dx.doi.org/10.1007/s00477-009-0324-0]

44. Moral, F.J.; Aguirado, C.; Alberdi, V.; Paniagua, L.L.; García-Martín, A.; Rebollo, F.J. Future scenarios for aridity under conditions of global climate change in Extremadura, Southwestern Spain. Land; 2023; 12, 536. [DOI: https://dx.doi.org/10.3390/land12030536]

45. Lotfi, M.; Kamali, G.A.; Meshkatee, A.H.; Varshavian, V. Study on the impact of climate change on evapotranspiration in west of Iran. Arabian J. Geosci.; 2020; 13, 722. [DOI: https://dx.doi.org/10.1007/s12517-020-05715-x]