Agricultural Drought Trends on the Iberian Peninsula: An Analysis Using Modeled and Reanalysis Soil Moisture Products

: Drought has a great impact on agriculture and food security, and climate change is increasing its frequency and exacerbating its intensity. Given the enormous interest in studying the recent drought evolution, in this work, agricultural drought trends over the past four decades on the Iberian Peninsula (IP) were analyzed. A trend analysis was performed with soil moisture based on the study of the evolution of anomalies and the Soil Water Deﬁcit Index (SWDI). Two soil moisture databases (Lisﬂood and ERA5-Land) were used and the analysis was performed at daily and weekly temporal scales. Climate characteristics and soil properties were also considered to detect whether a trend spatial pattern could be identiﬁed. The results have shown a clear predominance of negative trends. A marked temporal pattern with negative trends was obtained over a 10-month period that coincided with the growing season of most of the crops on the IP, while a positive trend was observed over 2 months. No differences were found based on the climatic zone or soil characteristics. However, negative trends were observed to decrease as the clay content increased. These results can provide useful information for better water management and agriculture of the IP and other Mediterranean areas.


Introduction
Drought is one of the natural hazards that has the most consequences and the greatest impact on human activities and population [1,2]. Under a scenario of uncertainty and climate change, it is expected that most water-related processes will be negatively affected, and the hydrological extremes will be exacerbated [3]. Global warming has accelerated hydrological processes, first by increasing the energy available for evapotranspiration and, second, by increasing the temperature and thus, the water holding capacity of the atmosphere [4]. Consequently, it results in more intense, widespread, and persistent extreme climatic events, such as droughts [5].
Among the three main categories of drought, i.e., meteorological, hydrological, and agricultural drought [6], agricultural drought has a more direct, immediate, and sensitive impact. An agricultural drought is considered to occur when the soil moisture availability to plants drops to such a level that it adversely affects the crop yield and, hence, agricultural production [7]. Consequently, drought is a major cause of limited productivity in rainfed agrosystems throughout the world, accounting for a large proportion of crop losses and annual yield variations in crops [8]. However, it is not just the direct impact on agriculture or food security [9]. Agricultural drought has been considered a direct factor of social and political conflicts in developing countries since it could have a catalytic effect, contributing to political unrest [10] and affecting the population in aspects that extend beyond feeding [11]. Southern Europe is expected to be particularly vulnerable to climate change, at least in all environmental and social aspects that depend on water resources [12]. Droughts can have strong impacts in water-limited areas, such as the Mediterranean region, particularly for countries relying on rainfed agricultural production and in areas in which irrigation plays an important role [13]. The Intergovernmental Panel on Climate Change (IPCC) warned about the possible intensification of droughts in this century, especially in some areas, such as the Mediterranean Basin [14]. Despite this proclivity, few studies have focused on studying agricultural drought trends in Mediterranean countries. Studies have focused on drought trends, particularly hydrological or meteorological drought [15]. Divrak et al. [16] studied the impact of drought on agricultural production in the Mediterranean Basin from a meteorological drought point of view. Jimenez-Donaire et al. [17] evaluated the impact of future changes in the temperature and precipitation regime of Spain on agricultural droughts. Peña-Gallardo et al. [18] studied the impact of drought on cereal productivity in Spain using several meteorological drought indices. Bouras et al. [19] evaluated the impact of agricultural drought on cereal yields in Morocco using remote sensing indices.
The Iberian Peninsula (IP) is among the especially vulnerable areas of southern Europe [14] due to its location as an ecotone between arid and humid regions, its climatic characteristics, and its agricultural predominance. Spain and Portugal account for 17.2% [20] of the utilized agricultural area of the European Union. Almost half (47.3%) of the IP territory has agricultural use, and only 4.8% is irrigated [21]. The rest of the agricultural land is rainfed, whether for permanent or annual crops or pastures, and exclusively depends on the water provided by the atmosphere and stored in the soil. Climate projections for the present century indicate overall increases in the area suitable for rainfed agriculture in temperate dryland regions, although the regional exception to this trend is Europe, especially the southern part, where suitability in temperate dryland areas will decline substantially [22].
According to its definition, agricultural drought is identified by a soil water shortage. Therefore, soil moisture observations are the appropriate approach for agricultural drought analysis, considering, as well, that plant water stress is more strongly related to the relative amount of plant-available water in the soil than to the absolute amount of soil moisture [23]. In the past, agricultural drought has usually been studied using methods based on longterm atmospheric data, such as rainfall and temperature, or on precipitation indices [24] that, in fact, were close to meteorological drought approaches. In recent years, new methods for agricultural drought monitoring and assessment using soil moisture have been proposed. These methods use soil moisture obtained from in situ measurements [25], modeling [26], reanalysis [27] or remote sensing [28] as the defining variable [29] or in a combined approach [30,31].
The main restriction for trend analysis is the availability of long enough data series. Therefore, agricultural drought trend analysis was difficult to address due to a lack of long series of appropriate soil moisture databases, in terms of spatial and temporal scales. However, long-term soil moisture databases are currently available, even globally, that allow for use in applications that were not possible until recently. These databases come from remote sensing [32], modeling or reanalysis [33,34], have been validated with very satisfactory results, and provide decades-long series of soil profile moisture data at everfiner temporal and spatial scales. In the case of remote sensing databases, the series still have a limited length since the launch of satellites dedicated specifically to soil moisture observations, such as Advanced Scatterometer (ASCAT, 2007), Soil Moisture and Ocean Salinity (SMOS, 2010) or Soil Moisture Active Pasive (SMAP, 2015), has been relatively recent. The Soil Moisture Climate Change Initiative soil moisture dataset [32], from the European Space Agency, has the longest series at more than 40 years and has been proven extremely useful for many applications [35]. However, the amount of data available before 2003 in most areas [36] is a limitation for their use in trend analysis. For example, on the IP, the percentage of data availability barely exceeded 20% in the first 20 years of the series.
Instead, modeling or reanalysis databases such as the Global Land Data Assimilation System (National Aeronautics and Space Administration) or ERA (European Centre for Medium-Range Weather Forecasts) provide several decades of series of global soil moisture data [37].
The main goal of this study has been to analyze agricultural drought trends over the past four decades on the IP, an area which is especially vulnerable and a period in which environmental conditions in this region are significantly changing due to climate change. Trend analysis has been performed with soil moisture based on the study of the evolution of anomalies, a statistic frequently used in this type of analysis, and the Soil Water Deficit Index (SWDI), an indicator that provides agronomic information, as it is based on the identification of plant available water [25]. For this work, two soil moisture databases that have proved to be suitable for the study of soil moisture trends (Lisflood and ERA5-Land) [38,39] were used and the analysis was performed at two temporal scales (daily and weekly). As a secondary objective, it was studied whether differences were observed in the results of trends depending on the climatic conditions or soil characteristics. Knowledge of agricultural drought trends over an environmentally vulnerable region, as the IP, can provide useful information for better water management and agriculture.

Databases and Methodology
The reanalysis dataset ERA5-Land (ERA5 hereafter) from the ECMWF provides temporal series of several land variables from 1981 to the present. These series are obtained by replaying the land component of ERA5 but using atmospheric variables as inputs. ERA5 variables, including soil moisture, are provided in a regular grid of 0.1 × 0.1 • , with an hourly temporal resolution [37]. This database has been validated with in situ soil moisture networks worldwide, including the IP [40]. For this study, soil moisture at 12 am and 12 pm and the three depth layers provided by ERA5 were averaged. Thus, the daily series of root zone soil moisture (0-100 cm) were obtained for the IP from 1981 to 2019.

Lisflood Soil Moisture
The Lisflood (LF) hydrological model, used by the European Commission for the prediction and monitoring of natural hazards in Europe, such as floods or droughts, provides outputs of hydrological variables throughout the entire European continent with a 5 × 5 km spatial resolution and a daily temporal resolution [41,42]. In this work, the soil moisture series corresponding to the IP was used from January 1991 to December 2019. This variable is provided in three different soil layers for each pixel. The first two layers were selected, containing the first 100 cm of soil depth. At each pixel, the daily soil moisture value of each of these layers was averaged, thus obtaining a time series of soil moisture for the root zone. This model has been widely used and validated in Europe [43].

Irrigation Mask
As soil moisture evolution in irrigated areas is controlled by irrigation management, a mask was constructed to exclude these regions from the study. The Food and Agriculture Organization (FAO) has developed the Digital Global Map of Irrigation Areas with a spatial resolution of 5 arc minutes [44]. This map provides the percentage of irrigated areas. In this study, version five of the irrigation map was used to create an irrigation area mask in the IP (Figure 1). First, this map was resampled into the two soil moisture database grids. Then, one mask for each soil moisture database was computed by excluding all the pixels with more than 10% irrigation area.

Climate Classification
The Köppen-Geiger climate classification was used to evaluate whether climatic conditions have an influence on agricultural drought trends. This climate classification assembles the different types of climate worldwide into five principal classes subdivided into up to 30 subcategories. The classification takes into account the temperature and precipitation characteristics of each kind of climate. This classification is widely used and different maps with different characteristics can be found in the literature [45][46][47]. In this study, the global map produced by [48] with a spatial resolution of 0.1° was used. This map was resampled into the two soil moisture database grids to obtain one classification for each soil moisture database ( Figure 1). In the resulting classification, only the principal classes were considered. Thus, pixels of the IP were classified in the three different types of climate that exist in the region: Arid, temperate, and cold.

Soil Data
The ISRIC-World Soil Information database includes the SoilGrids product, which provides information on different chemical and physical soil properties at a spatial resolution of 250 × 250 m [49]. The different soil properties are estimated for six depth intervals reaching down to 200 cm. In this study, the first five intervals (0-100 cm) of sand, clay, and organic matter content were averaged to obtain each depth profile and the resulting maps were resampled into ERA5 and LF grids. They were used to calculate the soil water parameters of the SWDI with the pedotransfer functions (PTFs) proposed by [50]. In addition, sand and clay maps ( Figure 2) were used to analyze whether a relationship between agricultural drought trends and the soil texture exists.

Climate Classification
The Köppen-Geiger climate classification was used to evaluate whether climatic conditions have an influence on agricultural drought trends. This climate classification assembles the different types of climate worldwide into five principal classes subdivided into up to 30 subcategories. The classification takes into account the temperature and precipitation characteristics of each kind of climate. This classification is widely used and different maps with different characteristics can be found in the literature [45][46][47]. In this study, the global map produced by [48] with a spatial resolution of 0.1 • was used. This map was resampled into the two soil moisture database grids to obtain one classification for each soil moisture database ( Figure 1). In the resulting classification, only the principal classes were considered. Thus, pixels of the IP were classified in the three different types of climate that exist in the region: Arid, temperate, and cold.

Soil Data
The ISRIC-World Soil Information database includes the SoilGrids product, which provides information on different chemical and physical soil properties at a spatial resolution of 250 × 250 m [49]. The different soil properties are estimated for six depth intervals reaching down to 200 cm. In this study, the first five intervals (0-100 cm) of sand, clay, and organic matter content were averaged to obtain each depth profile and the resulting maps were resampled into ERA5 and LF grids. They were used to calculate the soil water parameters of the SWDI with the pedotransfer functions (PTFs) proposed by [50]. In addition, sand and clay maps ( Figure 2) were used to analyze whether a relationship between agricultural drought trends and the soil texture exists.

Soil Moisture Climate Change Initiative (CCI) Dataset
The CCI soil moisture database [32] integrates satellite soil moisture products from different missions that have measured this variable since 1978. It was used to validate the ERA5 and LF series over the IP. The validation was carried out with the coincident period of the three datasets from 1991 to 2019 ( Figure 3). First, the reanalysis and modeled databases were resampled into the CCI soil moisture grid. The Pearson correlation coefficient was calculated for each pixel by comparing the CCI with the ERA5 series, the CCI with the LF series, and the ERA5 with the LF. In addition, the Köppen-Geiger climate classification was used to evaluate this correlation in each climatic type region ( Figure 3). Good concordance between the three databases was obtained, with correlation coefficients (R) above 0.8 most of the time. When ERA5 and LF were compared, R was approximately 0.9 for most of the series, which suggests that, despite the different methodologies behind those products, both are similarly suitable. Cold regions showed the worst correlation with CCI, which is Atmosphere 2021, 12, 236 5 of 17 due to the geographical location of this type of climate. In mountain regions, the retrieval of satellite soil moisture is difficult, which leads to a scarcity and lower accuracy of data. Nevertheless, cold areas are very small in the IP (Figure 1).

Soil Moisture Climate Change Initiative (CCI) Dataset
The CCI soil moisture database [32] integrates satellite soil moisture products from different missions that have measured this variable since 1978. It was used to validate the ERA5 and LF series over the IP. The validation was carried out with the coincident period of the three datasets from 1991 to 2019 ( Figure 3). First, the reanalysis and modeled databases were resampled into the CCI soil moisture grid. The Pearson correlation coefficient was calculated for each pixel by comparing the CCI with the ERA5 series, the CCI with the LF series, and the ERA5 with the LF. In addition, the Köppen-Geiger climate classification was used to evaluate this correlation in each climatic type region ( Figure 3). Good concordance between the three databases was obtained, with correlation coefficients (R) above 0.8 most of the time. When ERA5 and LF were compared, R was approximately 0.9 for most of the series, which suggests that, despite the different methodologies behind those products, both are similarly suitable. Cold regions showed the worst correlation with CCI, which is due to the geographical location of this type of climate. In mountain regions, the retrieval of satellite soil moisture is difficult, which leads to a scarcity and lower accuracy of data. Nevertheless, cold areas are very small in the IP (Figure 1). ure 3. Correlation of LF and ERA5 with climate change initiative (CCI) during the coinciding periods (1991-2019) and climate zones of the Iberian Peninsula.

Soil Moisture Anomalies
Anomalies of different hydrometeorological variables have been widely used as a drought indicator [51][52][53]. In this work, the root zone soil moisture anomalies were computed, and the trends of the resulting series were analyzed at daily and weekly resolutions. These two temporal scales were used considering that they are more suitable to analyze agricultural drought [54]. Anomalies were calculated by removing the mean annual cycle from the root zone soil moisture series. This annual cycle was calculated with a monthly average using the entire study period of each root zone soil moisture database,  The CCI soil moisture database [32] integrates satellite soil moisture products from different missions that have measured this variable since 1978. It was used to validate the ERA5 and LF series over the IP. The validation was carried out with the coincident period of the three datasets from 1991 to 2019 ( Figure 3). First, the reanalysis and modeled databases were resampled into the CCI soil moisture grid. The Pearson correlation coefficient was calculated for each pixel by comparing the CCI with the ERA5 series, the CCI with the LF series, and the ERA5 with the LF. In addition, the Köppen-Geiger climate classification was used to evaluate this correlation in each climatic type region ( Figure 3). Good concordance between the three databases was obtained, with correlation coefficients (R) above 0.8 most of the time. When ERA5 and LF were compared, R was approximately 0.9 for most of the series, which suggests that, despite the different methodologies behind those products, both are similarly suitable. Cold regions showed the worst correlation with CCI, which is due to the geographical location of this type of climate. In mountain regions, the retrieval of satellite soil moisture is difficult, which leads to a scarcity and lower accuracy of data. Nevertheless, cold areas are very small in the IP (Figure 1).

Soil Moisture Anomalies
Anomalies of different hydrometeorological variables have been widely used as a drought indicator [51][52][53]. In this work, the root zone soil moisture anomalies were computed, and the trends of the resulting series were analyzed at daily and weekly resolutions. These two temporal scales were used considering that they are more suitable to analyze agricultural drought [54]. Anomalies were calculated by removing the mean annual cycle from the root zone soil moisture series. This annual cycle was calculated with a monthly average using the entire study period of each root zone soil moisture database,  Anomalies of different hydrometeorological variables have been widely used as a drought indicator [51][52][53]. In this work, the root zone soil moisture anomalies were computed, and the trends of the resulting series were analyzed at daily and weekly resolutions. These two temporal scales were used considering that they are more suitable to analyze agricultural drought [54]. Anomalies were calculated by removing the mean annual cycle from the root zone soil moisture series. This annual cycle was calculated with a monthly average using the entire study period of each root zone soil moisture database, which is a common technique to calculate anomalies of seasonal variables [55]. For the daily temporal resolution, these monthly values were interpolated using a cubic function to obtain a value for each day of the year and the same method was used to calculate a weekly series, similar to what [56] did.

Soil Water Deficit Index
The Soil Water Deficit Index (SWDI) was developed for assessing and monitoring agricultural drought through soil moisture [25]. It is calculated as follows: where SM is the soil moisture content, FC is the field capacity value, and AWC (available water content) is the difference between FC and the wilting point (WP). The FC and WP values used to calculate the SWDI were obtained through the PTFs proposed by [50], and the sand, clay, and organic matter content values were obtained from the Soil Grids database. If the SWDI obtained is positive, this indicates that the soil has a water excess, while if the value obtained is negative, this indicates that there is agricultural drought. The water deficit is absolute (wilting point) when the SWDI reaches ≤−10, where the soil water content is below the lower limit of the available water for plants [57]. For this study, the SWDI was calculated with the soil moisture data from both the ERA5 and LF databases, both on a daily and weekly scale, for the study periods of each database.

Trend Analysis
The Mann-Kendall (MK) test [58,59] is a commonly used tool to detect monotonic trends [60][61][62][63] and in this research, it was used to study agricultural drought trends in the IP. MK is a nonparametric test able to detect whether a series has a significant monotonic trend and if it is positive or negative. For such detection, the S statistic is used and calculated as follows: where The sign of S determines whether the series has an upward (positive) or downward trend (negative), although the significance of this trend is given by the statistical parameter Z: For a Z value and a given significance level α, whether |Z| > Z α /2 is met, the trend is considered significant. In this study, the significance level was set to α = 0.05. Thus, a trend was considered significant for absolute values of Z greater than 1.96. The Z values were calculated for each pixel of the SWDI and anomaly series of each soil moisture database and each temporal resolution. The two agricultural drought index series considered have different characteristics. The SWDI is a seasonal variable, which means that the data are autocorrelated, as distinct from the anomaly series. This autocorrelation could increase the chance of obtaining a significant trend even if it does not exist [64]. Several methods have been proposed in the literature to avoid this problem [64][65][66]. However, [67] showed that for a large enough sample, the influence of this autocorrelation is not significant. Thus, Z values were first computed using the entire study period and the simple MK test and the number of significant, positive, and negative trends was analyzed. Then, each month series was individually analyzed and 12 Z values were obtained for each pixel. In this way, the annual dynamics of agricultural drought trends could be studied. Accordingly, the autocorrelation problem is avoided using this approach.

Results and Discussion
Analyzing the series as a whole and at the scale of the entire IP (Table 1), it is observed that the number of cases with a significant trend is very high, both for anomalies and for SWDI, ranging from 81.5% to 88.1% at the daily scale and from 48.9% to 70.8% at the weekly scale. For both anomalies and SWDI, and on the two scales of analysis, negative trends clearly dominate, ranging from 76.8% to 99.8% of the cases. These results clearly indicate that over the past four decades, agricultural drought has increased significantly in the main part of the IP. This is particularly interesting as both agricultural drought indices have a different nature. Soil moisture anomalies have a strictly statistical character, indicating whether the variable indicator significantly deviates from the usual condition for the specific site and period [68], while the SWDI has an agronomic and ecohydrological interpretation [25]. In addition, in view of these results, it was confirmed that the likely autocorrelation [67] that affects the SWDI series does not seem to influence the chance of obtaining significant trends. If the results are analyzed considering the climatic characteristics from the classification of Köppen-Geiger, the results are very similar (Table 2). In the three identified climate zones (arid, temperate, and cold), the predominance of negative trends is remarkable. The cold zone could be considered outside the analysis since it occupies a very small area of approximately 3% (Figure 2), mostly in high mountains. However, in both arid and temperate areas, which account for 37% and 60% of the IP surface area, respectively, the results are very similar. In all cases, and similar to the whole analysis, negative trends are prevalent and always above 76% of the cases. It should also be highlighted that in the case of the two indices calculated with the ERA5 database, negative trends are above 95%, both on a daily and weekly basis. Therefore, there do not seem to be noticeable differences between the climatic areas, depending on whether they are rainy or dry or whether the average temperature is high or low. Agricultural drought trends have also been analyzed for each month of the year (Tables 3 and 4 and Figures 4 and 5) using soil moisture anomalies and SWDI, with the aim of identifying whether any temporary patterns exist. There is a marked pattern that includes 10 months with negative trends and 2 months with positive trends, and it is similar for both agricultural drought indices calculated with ERA5 and LF. With the results obtained using ERA5 data, a clear predominance of negative trends was obtained from April to January and a positive trend was obtained from February to March. For the results obtained with LF, negative trends predominate from May to February and positive trends predominate in March-April. During the spring and summer months, negative trends usually exceed 70%, and in the case of ERA5, they are above 90%. The maximum number of significant cases was obtained from April-May to October, which coincides with the main part of the growing season.   A differentiation between the results obtained with both soil moisture products is well observed in the monthly maps (Figures 4 and 5). The clearest difference between the results of both products is observed in the summer months. While with ERA5, the spatial predominance of negative trends is absolute during those months in most IP, with LF, even if negative trends clearly predominate, there are certain areas with positive trends, especially from May to July. Most of these positively trend areas in those months coincide A differentiation between the results obtained with both soil moisture products is well observed in the monthly maps (Figures 4 and 5). The clearest difference between the results of both products is observed in the summer months. While with ERA5, the spatial predominance of negative trends is absolute during those months in most IP, with LF, even if negative trends clearly predominate, there are certain areas with positive trends, especially from May to July. Most of these positively trend areas in those months coincide with areas where soil texture is finer (Figure 2), which could indicate that this soil characteristic could be related to the behavior of agricultural drought indices calculated with LF. This differentiation may also be related to the intrinsic differences of both products. While LF soil moisture is a mere modeling product, ERA5 is a reanalysis product with satellite soil moisture assimilation. That is, in the first case, soil moisture is obtained from modeling from, fundamentally, climate data and the reliability of the results also depends on the suitability of the meteorological ground stations networks [68]. In the second, the modeling also implies the assimilation of the observed satellite soil moisture [69,70]. This different approach implies that with one approach, the existence of a precipitation period can generate a soil moisture behavior that could not be verified by another approach. This is often observed in water-limited environments, as is the case in much of the IP, where precipitation does not always cause significant changes in soil moisture evolution [71,72], especially when it is analyzed at the scale of the entire soil profile due to the existence of certain decoupling between the atmosphere and the soil system [73]. The increase in the frequency of low-precipitation events detected in Spain in recent decades [74] could also affect this dynamic.
Considering that soil plays an important role in relation to agricultural drought, as it is the reservoir from which plants are supplied, a trend analysis was also performed considering the soil characteristics. Thus, the evolution of the agricultural drought indices was analyzed according to soil textural fractions. The soil texture is largely responsible for porosity and plays an important role in hydraulic conductivity and soil water retention [75] and therefore in the availability of water for plants [76]. When agricultural drought trends were analyzed according to the sand content (Table 5), negative trends predominated for all textural percentage intervals. Therefore, no relationship was observed between the type of trend and the sand content of the soils ( Figure 6). The percentage of cases with a negative trend ranged from 69.1% to 100% for the two indices calculated with the two soil moisture products and for the two temporal scales studied.  However, in the case of soil clay content (Table 6), although there is a clear predominance of negative trends similarly, a pattern of relationship is observed. As the clay fraction increases, the percentage of significant results and cases with a negative trend decrease to the point that for LF-calculated indices, a positive trend predominates for soils with clay content above 40%. However, the percentage of territory with soils that have this high clay content is very low (0.2%). This inverse relationship was observed for both agricultural indices and soil moisture products (Figure 6), although it is highlighted in the case of LF. The coincidence between areas of a positive trend with LF in summer months However, in the case of soil clay content (Table 6), although there is a clear predominance of negative trends similarly, a pattern of relationship is observed. As the clay fraction increases, the percentage of significant results and cases with a negative trend decrease to the point that for LF-calculated indices, a positive trend predominates for soils with clay content above 40%. However, the percentage of territory with soils that have this high clay content is very low (0.2%). This inverse relationship was observed for both agricultural indices and soil moisture products (Figure 6), although it is highlighted in the case of LF. The coincidence between areas of a positive trend with LF in summer months was mentioned above, with areas where the clay content is higher. These results are in accordance with the role of the clay fraction in the hydrological functioning of soils. As the clay content increases, the water retention capacity increases, and the hydraulic conductivity decreases. Therefore, water flow and transfer to the atmosphere are also reduced [75]. Clay content has been shown to play a significant role in soil water dynamics, even if its percentage is not very high [71,77].   This predominance of negative trends in agricultural drought in the IP that have been observed over the past four decades has to be analyzed in view of the climate context that occurred in that period. In recent decades, there has been no clear trend in precipitation, both in Spain [78] or Portugal [79] or even during the last century [80], as was observed in all of Southwestern Europe [81]. This lack of precipitation trends was observed both in the amount and in the intensity of precipitation [82]. However, an increase in temperature that has been clearly detected globally has also been observed on the IP [83,84], even with very significant rates of up to 0.2 • C per decade [85]. This temperature increase is linked to a positive trend in evapotranspiration [86] and has a direct impact on soil water dynamics.
Therefore, if the amount of precipitation or its distribution has not changed but the temperature and thus evapotranspiration have experienced a noticeable increase, then an increase in aridity [87] and a negative evolution of the soil water balance are expected. This evolution appears to be reflected in the results of the analysis of the agricultural drought trends carried out over the IP. What has been observed in recent decades shows a troubling picture for agriculture in IP, but what is inferred from the different climate change scenarios makes the near future even more disturbing. In recent studies, the evolution of agricultural drought indicators showed how drought frequency, duration, and intensity will increase over the present century [17].
There is another process that is occurring in parallel with global warming and could contribute to the evolution of soil moisture and therefore agricultural drought in many areas of the IP, and in the same direction that the climate evolution and the results of the trend analysis point out. This is the abandonment of rural areas as a consequence of the depopulation that has been affecting large areas of Spain [88], Portugal [89], and all of Europe [90] in recent decades and will continue affecting in the next ones [91]. This abandonment has led to unprecedented vegetation regeneration and recovery of forested area [92]. It is well known that the type of land use and plant cover influences the dynamics and storage of soil moisture [93]. Under the same soil and bioclimatic conditions, a soil covered with herbaceous vegetation will retain and store more water than another covered with forest [94,95]. Therefore, this process of substitution of crops or pasture by permanent vegetation and, finally, forest, could also influence a gradual reduction in the soil water content in the areas affected by land abandonment.
This study has shown widespread predominance of negative trends in agricultural drought indices, which have primarily been detected for the period April-May to October. This period has great relevance considering that it coincides with the growing season of the main crops of the IP [96] and with the dry season in countries with a Mediterranean climate. This implies that an increase in soil water deficit occurs during the season of less water availability. This evolution could further restrict the conditions of rainfed agriculture, increase the need to increase the irrigation area to maintain productivity, and increase water consumption in the pre-existing irrigated areas [97]. This could mean a noticeable change in water resources management and a potential negative impact on natural water-linked systems [98]. Although some countries, such as Spain [99], have experienced a very positive evolution in the use of more efficient irrigation technologies, more effort and investment will still need to be made in this field. What is undoubted is that the impact of climate change is very evident and it is intensifying agricultural drought in regions such as the IP, where the rainfed area will remain the main part of the agricultural land.

Conclusions
Many studies in recent years, in addition to the findings of the IPCC analysis, prove that Southern Europe is expected to be particularly vulnerable to climate change. Therefore, studies on the evolution of processes directly affected by this major environmental problem are of great interest. Agricultural drought is at the heart of this interest. On the one hand, as it is a phenomenon in which water plays a prominent role and therefore is directly involved in many negative issues related to climate change. On the other hand, it is a natural hazard that affects aspects as sensitive as agricultural production and food security.
The Mediterranean region, particularly those countries relying on rainfed agricultural production, is especially vulnerable to these constraints triggered by agricultural drought and the observed evolution of the climate system.
The results of the study of agricultural drought trends on the IP over the past four decades have revealed a clear negative trend. Therefore, over this period, agricultural drought has become increasingly common. Whether using soil moisture anomalies, which are a mere statistical indicator, or using the SWDI, which is an index that expresses soil water availability for plants, have obtained similar results. In all cases, the number of significant trends was very high and more than 75% were negative. Substantial differences were also not observed when the analysis became specific to the climatic zones obtained by the classification of Köppen-Geiger. Regardless of the climatic characteristics of each zone, negative trends in agricultural drought predominate across the IP, which is consistent with the vulnerability to the climate drift of Mediterranean regions as a whole.
A temporal pattern was obtained from the results of the analysis throughout the months of the year. Despite some differences observed between the results of LF and ERA5, especially from May-July, and probably related to the specificity of the two soil moisture products used, the results were very similar. For both agricultural drought indicators, during 10 consecutive months, the negative trends were predominant and positive in 2 months. This finding is especially relevant since those months with a negative trend predominance coincide with the growing season of most of the crops of the IP. This indicates that for these crops, there is a clear decrease in the water available in the soil, just at the time when they need it most. The lower availability of water in this season is a defining feature of the Mediterranean climate. However, the results indicate that scarcity, water deficit and, therefore, water stress are being accentuated over time.
Soil is the reservoir from which crops are supplied with water. Therefore, it is interesting to determine whether there is any relationship between the evolution of agricultural drought and soil characteristics. Soil texture was the attribute that was chosen due to its enormous influence on the dynamics of soil moisture. A predominance of negative trends was obtained for all textural fractions. In the case of sand content, no related pattern was observed, and the results were similar for all textural intervals. However, a pattern for clay was obtained, so that as its content increased, the percentage of negative trends decreased. This is in line with the behavior of clay soils in relation to water flow and retention capacity.
The results of the analysis of agricultural drought trends on the IP must be interpreted in view of the evolution of climate variables and the evolution of land use that have been observed in this region in recent decades. All these changes directly and negatively affect the soil water balance. Therefore, and if this continues, as climate change projections point out, the agriculture of the IP will be seriously affected in the coming years. Not only rainfed agriculture, which directly depends on the water that supplies the atmosphere and is stored in the soil, but also irrigation since the water demand will increase. Currently, approximately three-quarters of the water demand on the IP is used in agriculture. Therefore, the results of this study can be useful for properly planning water management in agriculture and implementing appropriate climate change adaptation measures.