Analysis of Agronomic Drought in a Highly Anthropogenic Context Based on Satellite Monitoring of Vegetation and Soil Moisture

: This paper aims to analyze agronomic drought in a highly anthropogenic, semiarid region, the western Mediterranean region. The proposed study is based on Moderate-Resolution Imaging Spectroradiometer (MODIS) and Advanced SCATterometer (ASCAT) satellite data describing the dynamics of vegetation cover and soil water content through the Normalized Difference Vegetation Index (NDVI) and Soil Water Index (SWI). Two drought indices were analyzed: the Vegetation Anomaly Index (VAI) and the Moisture Anomaly Index (MAI). The dynamics of the VAI were analyzed as a function of land cover deduced from the Copernicus land cover map. The effect of land cover and anthropogenic agricultural activities such as irrigation on the estimation of the drought index VAI was analyzed. The VAI dynamics were very similar for the shrub and forest classes. The contribution of vegetation cover (VAI) was combined with the effect of soil water content (MAI) through a new drought index called the global drought index (GDI) to conduct a global analysis of drought conditions. The implementation of this combination on different test areas in the study region is discussed.


Introduction
In semiarid areas, drought is a very frequent phenomenon that generates very serious problems, especially for agriculture and food security [? ? ? ]. This scenario often worsens with climate change and an increase in the frequency of extreme events [? ? ? ]. It is essential to properly quantify the intensity and the beginning and end of a drought event to properly assess its effects. In this context, several scientific studies have developed drought indices, mainly based on the detection of climatic anomalies [? ? ? ], mainly related to a deficit in precipitation compared to the average precipitation [? ].
Based on meteorological data collected from weather stations, several works have proposed different types of drought indices, such as the Palmer drought severity index (PDSI; [? ? ]) or the standardized precipitation index (SPI; [? ]). These approaches are often considered operational by meteorological services. However, the information from these approaches, despite their considerable relevance, may be weakened by the low density of rain gauges in some parts of the world. They may also not be fully correlated with the reality of agronomic drought, which directly affects vegetation cover development.
Within this framework, several studies have been developed over the last twenty years based on a time series of remote sensing satellite data [? ? ? ]. In fact, remote sensing has shown high potential to retrieve land surface properties [? ? ? ? ]. These studies are mainly related to the existence of increasingly longer time series, sometimes more than thirty years that are able to statistically determine climatic anomalies. These data and the products proposed with them mainly describe the dynamics of vegetation cover [? ? ? ], soil water content [? ? ? ], and surface temperature [? ].

TERRA-MODIS Data
In this study, we used 250-m spatial resolution, 16-day composites of MODIS NDVI data (MOD13Q1, collection 5). This product was retrieved from atmosphere-corrected [? ], daily, bidirectional surface reflectance observations using a compositing technique based on product quality assurance metrics to remove low-quality pixels. The 250-m spatial resolution is the finest available from the MODIS NDVI dataset, and the 16-day composite was selected to ensure a high probability of having the best quality pixel representing the NDVI within each 16-day period and without the effect of clouds [? ]. Here, we used this product to analyze the influence of drought between September 2000 and the end of 2020.

ASCAT Data
The ASCAT scatterometer radar is one of the 12 instruments carried by the ESA's METOP-A, METOP-B, and METOP-C satellites and operates in the C-band (5.3 GHz), in the vertical polarization. Over land, the measured radar backscattering coefficient depends on the soil moisture, surface roughness, vegetation characteristics, and incidence angle of the transmitted radar beam. The soil moisture data are retrieved from the backscattering coefficient using a change detection method developed at the Institute of Photogrammetry and Remote Sensing (IPF), Vienna University of Technology (TU-Wien) [? ? ]. The derived soil moisture product, expressed in relative units and referred to as 'surface soil moisture' (ms), represents the water content in the first 5 cm of the soil and ranges between the extremes corresponding to totally dry conditions and a saturated water capacity. The spatial resolution is defined by cells of approximately 50 km, with a grid spacing of 12.5 km, and the temporal resolution of the data varies between approximately two and three measurements per week.
The Soil Water Index data (SWI) was derived from m using, Equation (1), and represents the root-zone soil moisture content in the first meter of the soil in relative units ranging between wilting point and field capacity [? ].
where ms is the surface soil moisture estimate from the ASCAT scatterometer at time t i . The parameter T, called the characteristic time length, represents the time scale of soil moisture variations in units of time. T equal to 20 days has shown the best fit to ground measurements. These products have already been used and validated in the studied site [? ].

Copernicus Land Cover Map
The land use map proposed by the Copernicus program [? ] provides spatial information for different land cover classes or natural cover of land surfaces, such as forests, agricultural areas, wetlands, and urban areas. To simplify our study, we considered only the land cover map proposed in 2020, assuming little change in land cover at the western Mediterranean scale. Figure 1 illustrates the land cover map for our study site, with four main land cover types: bare soil often related to desert areas, crop agricultural areas, shrub and herbaceous vegetation, and forests.

Rif
11% 17%  49%  470  Andalusia  69%  8%  14%  29 These sites were approximately 50 × 50 km with the different land uses (crops, shrub/herbaceous vegetation, and open forests), and the drought index VAI was calculated for the different land uses at a monthly scale. We assumed that the climate context was the same among the three land covers since the regions are limited in size. . This product was initially developed by Döll and Siebert [46] and then enhanced to obtain 5 arc minutes resolution maps [47]. The reference years for the statistics used to establish the map differ between countries but in most cases, the data refer to the period 2000-2008. Despite several shortcomings detailed by the authors, this map remains the FAO reference today. Figure 2 represents the percentage of area equipped with irrigation from 0 to 100%. The 0% corresponds to rainfed areas. The 100% corresponds to areas fully irrigated. Three types of land use were identified for the analysis of drought proposed in this study: agricultural areas, shrub and herbaceous vegetation areas, and open forest areas. A part of this analysis was based on an evaluation of the effects of land use at four sites with different characteristics in terms of land use, climatology, and type of agriculture. The four sites were the region of Kairouan in Tunisia, characterized by rainfed agriculture, a low rainfall amount (300 mm per year); the region of Ain Defla in Algeria (640 mm per year); the Rif, a mountainous region of northern Morocco, with a high precipitation amount (800 mm per year); and Andalusia in southern Spain (480 mm per year) with intensively irrigated agriculture. The climate of Kairouan is semiarid (BSh on the Köppen-Geiger classification), while the three other areas have a Mediterranean climate (CSa). Table 1 illustrates the percentages of crops, shrubs, herbaceous vegetation and open forests, and also the average altitude for each site. These sites were approximately 50 × 50 km with the different land uses (crops, shrub/herbaceous vegetation, and open forests), and the drought index VAI was calculated for the different land uses at a monthly scale. We assumed that the climate context was the same among the three land covers since the regions are limited in size. Figure 2 illustrates the Global Map of Irrigated Area (GMIA) version 5 for the study area [? ]. This product was initially developed by Döll and Siebert [? ] and then enhanced to obtain 5 arc minutes resolution maps [? ]. The reference years for the statistics used to establish the map differ between countries but in most cases, the data refer to the period 2000-2008. Despite several shortcomings detailed by the authors, this map remains the FAO reference today. Figure 2 represents the percentage of area equipped with irrigation from 0 to 100%. The 0% corresponds to rainfed areas. The 100% corresponds to areas fully irrigated.

Methodology
Three drought indices were used in this study to understand the drought phenomena, and the two VAI and MAI indices had already been tested in the study region [37]. The global drought index (GDI) was proposed in this study.

Vegetation Anomaly Index (VAI)
This index is based on statistics derived from the NDVI MODIS time series and is referred to as the 'vegetation anomaly index' (VAI) [31], written as: where NDVIi is the NDVI estimate for a given month i, (NDVIi)mean is the mean value of the NDVI during month i, derived from the previously described 20 years of NDVI time series, and σi corresponds to the standard deviation of the NDVI values estimated for month i over the same 20-year period. When the VAI is greater than zero, it means that NDVI is higher than normal, corresponding to a bountiful season. When the VAI is negative, NDVI is lower than normal, which in many cases may be attributed to a drought or a period with a lack of precipitation. We have computed this index for each processed pixel.

Moisture Anomaly Index (MAI)
With vegetation development anomalies in the context of drought conditions, the soil moisture content is complementary information essential for detecting agronomic drought.
For soil water content, we propose the use of ASCAT SWI products. With a 12.5 km spatial resolution, this product remains highly adapted to the monitoring of climatic effects related to precipitation. The MAI drought index [23] is based on statistics derived from the SWI time series and is referred to as the 'Moisture Anomaly Index' (MAI), written as: where SWIi is the SWI estimate for a given month i, (SWIi)mean is the mean value of the SWI during month i, derived from the previously described 13 years of SWI time series, and σi

Methodology
Three drought indices were used in this study to understand the drought phenomena, and the two VAI and MAI indices had already been tested in the study region [? ]. The global drought index (GDI) was proposed in this study.

Vegetation Anomaly Index (VAI)
This index is based on statistics derived from the NDVI MODIS time series and is referred to as the 'vegetation anomaly index' (VAI) [? ], written as: where NDVI i is the NDVI estimate for a given month i, (NDVI i ) mean is the mean value of the NDVI during month i, derived from the previously described 20 years of NDVI time series, and σ i corresponds to the standard deviation of the NDVI values estimated for month i over the same 20-year period. When the VAI is greater than zero, it means that NDVI is higher than normal, corresponding to a bountiful season. When the VAI is negative, NDVI is lower than normal, which in many cases may be attributed to a drought or a period with a lack of precipitation. We have computed this index for each processed pixel.

Moisture Anomaly Index (MAI)
With vegetation development anomalies in the context of drought conditions, the soil moisture content is complementary information essential for detecting agronomic drought.
For soil water content, we propose the use of ASCAT SWI products. With a 12.5 km spatial resolution, this product remains highly adapted to the monitoring of climatic effects related to precipitation. The MAI drought index [? ] is based on statistics derived from the SWI time series and is referred to as the 'Moisture Anomaly Index' (MAI), written as: where SWI i is the SWI estimate for a given month i, (SWI i ) mean is the mean value of the SWI during month i, derived from the previously described 13 years of SWI time series, and σ i corresponds to the standard deviation of the SWI values estimated for month i over the same 13-year period. When the MAI is greater than zero, a higher than normal SWI value is indicated, corresponding to a wet profile and the absence of drought. When the MAI is negative, Remote Sens. 2021, 13, 2698 6 of ?? a lower than normal SWI value, which is probably the result of drought or a period with a lack of precipitation, is indicated.

Global Drought Index (GDI)
To consolidate drought detection, we proposed a mixed index based on the two indices VAI and MAI [? ]. This type of combination has been proposed by different studies [? ?
]. However, as stated in the introduction, the weightings of the two indices are generally based on simple addition without discussing different weights.
The index is written as: where αi and βi are the weights given to the two drought indices during month i. Ind i is the Ind estimate for a given month i σ Ind,i corresponds to the standard deviation of the Ind values estimated for month i over the same 13-year period. The GDI i is a dimensionless index ranging theoretically from negative values (drier than normal) to unlimited positive values (wetter than normal). For the three drought indices, MAI, VAI, and GDI, we could consider a context of severe drought if the index was less than −1, a context close to normal for indices between −1 and 1, and finally in a wet situation for an index superior to 1.
Droughts are generally marked by a decrease in NDVI and SWI products. In a semiarid climate, we observe periods of the year where very low NDVI or SWI are present independently of drought. This is typical of the summer season. However, in the autumn and spring seasons, the averages of one of the two parameters are also equally low. In this context, weights αi and βi could be calculated according to the mean level of vegetation or moisture in month i. This means that a higher weight should be given to the vegetation anomaly when the vegetation is well developed. This is expressed by Equation (5), where i, min and max of NDVI correspond to the 12 months of the year. Thus, a low αi is given to poorly vegetated months; one weighs the VAI more than the MAI in the calculation of the global drought index. This could be introduced in the weights proposed for the two indices MAI and VAI in the calculation of GDI.
where NDVI i is the NDVI estimate for a given month i, NDVI min is the minimum value of the NDVI for the 12 months, derived from the previously described 20 years of NDVI time series, and NDVI max is the maximum value of the NDVI for the 12 months, derived from the previously described 20 years of the NDVI time series.
where SWI i is the SWI estimate for a given month i, SWI min is the minimum value of the SWI for the 12 months, derived from the previously described 13 years of SWI time series, and SWI max is the maximum value of the SWI for the 12 months, derived from the previously described 13 years of SWI time series. Interestingly, as these values are calculated independently, adding αi and βi is not equal to 1, so the normalization of Equation (4) is necessary to obtain the final GDI.

Statistical Parameters
In the present study, the comparison between indices is evaluated with the Pearson's correlation, and the RMSE (Root Mean Square Error).
where x i and y i are individual samples taken at points indexed with the variable i, N is the number of samples, x is the mean of the samples x i , and y is the mean of the samples y i . RMSE (Root Mean Square Error) where N is the number of data samples, P i is the first index value of sample i, and O i is the second index value of sample i. The data are checked to ensure they follow a normal distribution, and for their significance, prior to all analyses. The intercomparison between indices was assessed using the RMSE error.

Analysis of VAI Variation Function of Land Use
Drought indices using NDVI satellite time series are generally calculated in study areas or pixels, taking into account only the NDVI level. However, the temporal variations in the NDVI are strongly related to the type of land use. Thus, it is legitimate to ask what the effect of these land-use types is on the level of the drought index. This is particularly the case in agricultural areas, which can be affected by additional factors to climatic conditions, such as human activity, presence or absence of irrigation, and decisions made by farmers in the context of extreme conditions. These anthropogenic factors could affect the identified anomalies. Thus, the objective of this section is to evaluate the VAI index for the different types of land use and to assess the behaviors of each. Figure 3 shows a cross-comparison of the VAI between the different land uses at the four test sites identified in Section 2.1. All monthly VAI estimates are displayed, each point corresponds to the VAI calculated for a given month. In each figure, the statistical parameters R 2 , RMSE, and the linear relation between the VAI of the two considered land uses (Y = a.x, a is the slope) are illustrated. For the case of two land uses with very close VAI indices, we should find an RMSE close to zero, an R 2 close to 1, and a slope close to 1.
Apart from Kairouan, the slope between the cropland class and the two other crop classes was lower than 0.9. The slope was steeper for forests. The observed behaviors could be explained by a maximum anthropogenic effect in agricultural zones, which probably strengthens the extremes. In fact, in the event of drought, in several regions, including irrigated areas, farmers do not sow if the first rainfall is absent and may even abandon the fields to animals during a very dry season. This human action accentuates the decline in NDVI and thus the effects of drought. Human influence is obviously less important in natural areas.
On the other hand, the correlation when comparing the VAI values of crop zones with other cover types showed very distinctive patterns. The semiarid area with limited irrigation in Kairouan showed high correlations R 2 of 0.91 and 0.88 between crops and shrub/herbaceous vegetation and crops and forests, respectively, and little dispersion (RMSE of 0.2 and 0.23, respectively). In contrast, Andalusia had low correlations and high dispersion. The Moroccan and Algerian sites had moderate correlations. These correlations were weakened by anthropogenic effects, particularly due to irrigation, which could be important in some areas. be important in some areas.
For studied sites, we observed a slope generally close to 1 between the VAI estimated for shrub/herbaceous vegetation and that estimated for forests. In Andalusia, this slope was weaker. The limited % of these land uses could explain this result; however, it could also have occurred based on the dominance of olive groves in this area. A strong correlation was also observed between shrubs and forests for the first three sites, Kairouan, Ain Dafla, and Rif. This correlation R 2 ranged between 0.9 and 0.97 for the three test sites.  For studied sites, we observed a slope generally close to 1 between the VAI estimated for shrub/herbaceous vegetation and that estimated for forests. In Andalusia, this slope was weaker. The limited % of these land uses could explain this result; however, it could also have occurred based on the dominance of olive groves in this area. A strong correlation was also observed between shrubs and forests for the first three sites, Kairouan, Ain Dafla, and Rif. This correlation R 2 ranged between 0.9 and 0.97 for the three test sites.
Distributions of the VAI for the three types of land use at the four sites (Figure 4) confirmed the first results of the comparisons, with generally similar behavior for shrubs, herbaceous vegetation, and forests. This figure also shows a large contrast between the different regions. In particular, the Andalusia distribution was concentrated around an average situation, while the flattened Kairouan's distribution shows the frequent recurrence of extreme events.
(d) Distributions of the VAI for the three types of land use at the four sites (Figure 4) confirmed the first results of the comparisons, with generally similar behavior for shrubs, herbaceous vegetation, and forests. This figure also shows a large contrast between the different regions. In particular, the Andalusia distribution was concentrated around an average situation, while the flattened Kairouan's distribution shows the frequent recurrence of extreme events. In conclusion, we observe a strong resemblance between the two classes shrub, herbaceous vegetation and openforest, compared to cropland with rather different behaviors. So, two global classes were used to illustrate the VAI, one class of agricultural areas and a second class, natural vegetation merging the two classes of natural areas (shrub, herbaceous vegetation and forest). Figure 5 illustrates the two proposed classes, cropland and natural vegetation over the studied region. The effect of topography was not analyzed in this study. However, it could have an effect on drought-related behaviors. In fact, natural areas are more present in areas with high altitudes than in valleys. These areas also generally have lower temperatures, which could allow for better resistance to drought. a second class, natural vegetation merging the two classes of natural areas (shrub, herbaceous vegetation and forest). Figure 5 illustrates the two proposed classes, cropland and natural vegetation over the studied region. The effect of topography was not analyzed in this study. However, it could have an effect on drought-related behaviors. In fact, natura areas are more present in areas with high altitudes than in valleys. These areas also generally have lower temperatures, which could allow for better resistance to drought.  The VAI and irrigation mapping are magnified for two sites, Kairouan ( Figure 7) and Andalusia (Figure 8). The first site (Kairouan) had limited irrigated agriculture. The second site (Andalusia) was characterized by intense irrigation. Our objective was to observe the differences between VAI indices in the context of combined effects related to land use (crop or natural vegetation) and then the presence or not of irrigation in the crop context At the first site, for the considered dates, we observed small differences between the estimated VAI for cropland and natural surfaces for some dates. In the context of the two extremes, in 2001 with an exceptional drought and in 2018 with wet conditions, the two types of land uses showed approximately the same behavior.
At the second site, there was clearly a strong difference between the drought indices found in the natural and cropland areas. There were even reversed trends for some dates such as 2016. This effect was particularly accentuated due to the highly irrigated context The irrigation map over the study area confirmed this accentuated effect on the regions marked by the highest rate of irrigation.
In the context of this study, the considered irrigation map is the FAO reference map However, it is based on data collected more than fifteen years ago. This could have some lags compared to the current context, generally marked by an intensification of irrigation especially in North Africa. This could generate some errors in the interpretation of the behavior results of irrigated and non-irrigated areas.   The VAI and irrigation mapping are magnified for two sites, Kairouan ( Figure 7) and Andalusia (Figure 8). The first site (Kairouan) had limited irrigated agriculture. The second site (Andalusia) was characterized by intense irrigation. Our objective was to observe the differences between VAI indices in the context of combined effects related to land use (crop or natural vegetation) and then the presence or not of irrigation in the crop context.  Over the Andalusia site, Figure 9 illustrates the intercomparison between the estimated VAI in areas irrigated at +20% and other rainfed agricultural areas or areas with irrigation levels below 20%. Higher extremes were generally observed for areas with no or low irrigation levels. This result could be explained by a climate with more extreme conditions related to drought or heavy rainfall. The irrigated context was less variable (amplitude of 1.88 versus 2.65) because of the organized irrigation management. Both curves had a negative trend, which was much higher for non-irrigated areas (−0.0006 per month) than for irrigated areas (−0.0013). This result could be explained by the strong conversion of seasonal crops to olive groves in Andalusia. According to agricultural statistics, the share of olive groves increased by 30% from 1993 to 2018, while the share of seasonal crops decreased by 43% during the same period [49]. Marked differences between the two curves were observed, especially in the summer months. In 2005 and 2012, the VAI of the irrigated area was significantly higher (>0.5) than that of the non-irrigated area. In 2005, less rain occurred than had previously in the past 20 years, and according to agricultural statistics, the sunflower, a typical summer crop that represents an average of At the first site, for the considered dates, we observed small differences between the estimated VAI for cropland and natural surfaces for some dates. In the context of the two extremes, in 2001 with an exceptional drought and in 2018 with wet conditions, the two types of land uses showed approximately the same behavior.
At the second site, there was clearly a strong difference between the drought indices found in the natural and cropland areas. There were even reversed trends for some dates, such as 2016. This effect was particularly accentuated due to the highly irrigated context. The irrigation map over the study area confirmed this accentuated effect on the regions marked by the highest rate of irrigation.
In the context of this study, the considered irrigation map is the FAO reference map. However, it is based on data collected more than fifteen years ago. This could have some lags compared to the current context, generally marked by an intensification of irrigation, especially in North Africa. This could generate some errors in the interpretation of the behavior results of irrigated and non-irrigated areas.
Over the Andalusia site, Figure ?? illustrates the intercomparison between the estimated VAI in areas irrigated at +20% and other rainfed agricultural areas or areas with irrigation levels below 20%. Higher extremes were generally observed for areas with no or low irrigation levels. This result could be explained by a climate with more extreme Remote Sens. 2021, 13, 2698 13 of ?? conditions related to drought or heavy rainfall. The irrigated context was less variable (amplitude of 1.88 versus 2.65) because of the organized irrigation management. Both curves had a negative trend, which was much higher for non-irrigated areas (−0.0006 per month) than for irrigated areas (−0.0013). This result could be explained by the strong conversion of seasonal crops to olive groves in Andalusia. According to agricultural statistics, the share of olive groves increased by 30% from 1993 to 2018, while the share of seasonal crops decreased by 43% during the same period [? ]. Marked differences between the two curves were observed, especially in the summer months. In 2005 and 2012, the VAI of the irrigated area was significantly higher (>0.5) than that of the non-irrigated area. In 2005, less rain occurred than had previously in the past 20 years, and according to agricultural statistics, the sunflower, a typical summer crop that represents an average of 14% of the total cropped area of Andalusia, experienced a severe drop in the cropped area in those two years. Jimenez et al.  . VAI time series at the Andalusia site in pixels with more than 20% of irrigated area and in pixels mainly rainfed or less than 20% of irrigated area. Figure 10 illustrates the NDVI, SWI, VAI, and MAI estimated at two test sites, Kairouan and Andalusia. It is important to note that with a resolution of 12.5 km, the SWI cannot be identified in a completely separate way between the two types of land use. This scenario explains the strong similarities in the time series for the two types of land cover in fragmented landscapes, while there was a higher contrast in adjacent landscapes.

Application of GDI Index
At the Kairouan site (Figure 10a), the trends of the two indices seemed to be strongly correlated for both land cover types in cropland and natural vegetation. In fact, at this site, the type of rainfed agriculture explained this correlation between the water stock in the soil and the dynamics of the vegetation cover. On the other hand, at the irrigated site in Andalusia (Figure 10b), strong differences were first observed between the two types of land cover linked to different vegetation cycles. The MAI and VAI also showed different trends. This result could be explained by the effect of irrigation on the dynamics of the vegetation cover, independent of the actual rainfall context. In this context, the application of the GDI index aims to consider the effects related to vegetation and soil water content through the combination of VAI and MAI. Figure 9. VAI time series at the Andalusia site in pixels with more than 20% of irrigated area and in pixels mainly rainfed or less than 20% of irrigated area. Figure ?? illustrates the NDVI, SWI, VAI, and MAI estimated at two test sites, Kairouan and Andalusia. It is important to note that with a resolution of 12.5 km, the SWI cannot be identified in a completely separate way between the two types of land use. This scenario explains the strong similarities in the time series for the two types of land cover in fragmented landscapes, while there was a higher contrast in adjacent landscapes.  At the Kairouan site (Figure ??a), the trends of the two indices seemed to be strongly correlated for both land cover types in cropland and natural vegetation. In fact, at this site, the type of rainfed agriculture explained this correlation between the water stock in the soil and the dynamics of the vegetation cover. On the other hand, at the irrigated site in Andalusia (Figure ??b), strong differences were first observed between the two types of land cover linked to different vegetation cycles. The MAI and VAI also showed different trends. This result could be explained by the effect of irrigation on the dynamics of the vegetation cover, independent of the actual rainfall context. In this context, the application of the GDI index aims to consider the effects related to vegetation and soil water content through the combination of VAI and MAI. Figure ?? illustrates parameters α and β defined in Section 2.2.3 for α0 and β0, respectively, equal to 0.1 at the two study sites in Kairouan and Andalusia. The values of α0 and β0 were chosen empirically to define the minimum weights of the VAI and MAI indices. The larger these values, the more we move towards an identical weighting of the two indices as applied by other existing drought indices.

GDI Combined Index Application
For the former, slight discrepancies between the different weights proposed for the two types of cover were observed. This result was related to the rainfall context of this study area and the different spatial resolutions of the two remote sensing products used. In Kairouan, the alpha cycles were similar between cropland and natural vegetation; however, some differences were observed. In particular, 50% of the boxes were much smaller for the natural vegetation between October and January, which illustrates more variability during this period when the start of the agricultural season with the first cereal seeding was very variable between farmers and more variability from one year to another due to rainfall.
Minimum weights were observed in July and August and corresponded to the weight with the lowest vegetation cover and the lowest soil water storage. These weights were close to 0.1, the minimum, and illustrate the limitations of the calculation of a drought index in this context with too weak effects on the almost nonexistent vegetation cover and the very low water stock.
In Kairouan (Figure ??a), variations in α and β were rarely above 0.2 for cropland and natural vegetation; however, α showed wide variation for Andalusia croplands, probably due to the high percentage of irrigated areas and intense crop rotations, while α variation for the natural vegetation areas remained under 0.2. Slight differences in the weights existed in autumn and spring. The Andalusia site (Figure ??b) illustrates additional uncorrelated relationships between canopy dynamics and water content stock, particularly due to irrigation in agricultural areas. A higher dispersion of β was observed in the months of March and April, suggesting that the low-resolution remote sensing product of soil moisture could identify the impact of irrigation. In this framework, weights α and β are only correlated in the case of natural canopies. Figure ?? shows the values of VAI, MAI, and GDI at the two sites Kairouan and Andalusia, during the last 13 years. At the Kairouan site (Figure ??a), we see strong similarities between the drought indices estimated on the two classes, croplands and natural vegetation, in the context of mainly rainfed agriculture. At the Andalusia site (Figure ??b), the correlation between VAI and MAI is weaker. The effect of irrigation is probably very important. GDI more clearly combines the effects of drought on vegetation and soil water content.
ities between the drought indices estimated on the two classes, croplands and natural vegetation, in the context of mainly rainfed agriculture. At the Andalusia site (Figure 10b), the correlation between VAI and MAI is weaker. The effect of irrigation is probably very important. GDI more clearly combines the effects of drought on vegetation and soil water content.  Figure 12 illustrates the mapping of the GDI with the proposed weights described in Equations (5) and (6) in March 2016 (dry year) and March 2018 (wet year). The patterns were approximately similar to those in the VAI maps of Figure 6. However, compared to the VAI map, the new map accentuates the effect of drought in 2016 on the entire study region. It also shows a wetter index in 2018. As noted, with GDI, we have a drought map that combines the effects of drought on vegetation and on the water content in the soil. For example, on the Andalusia site, we observe a context with a positive GDI index in 2018 when it was negative with the VAI index. This is explained by a positive MAI index  Figure ?? illustrates the mapping of the GDI with the proposed weights described in Equations (5) and (6) in March 2016 (dry year) and March 2018 (wet year). The patterns were approximately similar to those in the VAI maps of Figure 6. However, compared to the VAI map, the new map accentuates the effect of drought in 2016 on the entire study region. It also shows a wetter index in 2018. As noted, with GDI, we have a drought map that combines the effects of drought on vegetation and on the water content in the soil. For example, on the Andalusia site, we observe a context with a positive GDI index in 2018 when it was negative with the VAI index. This is explained by a positive MAI index level as illustrated in Figure ??. When we look closely at the sites of Kairouan (Figure ??) and Andalusia (Figure ??), for example, during the entire hydrological year 2015-2016, which was particularly dry, we observed approximately the same spatial trends between the two indices VAI and GDI. However, the intensity of the indices was not the same. We observe similar behaviors during the summer dry season. On the other hand, during the other seasons, the intensity of the drought was not the same. The difference is more marked on the highly irrigated Andalusia site, sometimes with an intensification of the drought index as observed in autumn and inversely at the end of spring.

Conclusions
A drought analysis was conducted over regions in North Africa and sou using satellite time series data, with NDVI products from optical MODIS da water content products from ASCAT microwave measurements. Drought a conducted first using the VAI index based on the NDVI. The effect of land cov

Conclusions
A drought analysis was conducted over regions in North Africa and southern Spain using satellite time series data, with NDVI products from optical MODIS data and soil water content products from ASCAT microwave measurements. Drought analysis was conducted first using the VAI index based on the NDVI. The effect of land cover type was discussed. There was a general increase in the drought indices over agricultural areas compared to that over other natural land covers (shrub and forests). The drought indices calculated in these agricultural areas showed less correlation with the indices calculated for the other natural land cover types compared to comparisons made between the natural cover types. This result can be explained by the anthropogenic effect, which can be relatively important in agricultural areas, particularly due to irrigation. The mapping of the VAI index was estimated separately for the two land-use classes, croplands and natural areas, which were composed of shrubs and forests.
To complete the information related to the vegetation cover, a new index, the GDI that combined NDVI information with water content information through the SWI was also proposed. This index is the combination of the VAI and MAI indices, weighted with variable weights depending on the seasonal context of precipitation and vegetation development or human action such as irrigation. The purpose of this variable weighting was to take into account more or less important effects according to the monthly average of the variable in each study month. Mapping of the GDI parameters was proposed over the studied region.