Articles | Volume 25, issue 2
Hydrol. Earth Syst. Sci., 25, 755–768, 2021

Special issue: Data acquisition and modelling of hydrological, hydrogeological...

Hydrol. Earth Syst. Sci., 25, 755–768, 2021

Research article 18 Feb 2021

Research article | 18 Feb 2021

Long-term water stress and drought assessment of Mediterranean oak savanna vegetation using thermal remote sensing

Long-term water stress and drought assessment of Mediterranean oak savanna vegetation using thermal remote sensing
María P. González-Dugo1, Xuelong Chen2,3, Ana Andreu1, Elisabet Carpintero1, Pedro J. Gómez-Giraldez1, Arnaud Carrara4, and Zhongbo Su5 María P. González-Dugo et al.
  • 1IFAPA, Consejería de Agricultura, Ganadería, Pesca y Desarrollo Sostenible, Apdo. 3048, 14071, Córdoba, Spain
  • 2Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, China
  • 3CAS Center for Excellence in Tibetan Plateau Earth Sciences, Chinese Academy of Sciences, Beijing, China
  • 4Fundación CEAM. Parque Tecnológico, Calle Charles Darwin 14, 46980 Paterna, Valencia, Spain
  • 5Faculty of Geo-Information Science and Earth Observation, University of Twente, Enschede, the Netherlands

Correspondence: María P. González-Dugo (


Drought is a devastating natural hazard that is difficult to define, detect and quantify. The increased availability of both meteorological and remotely sensed data provides an opportunity to develop new methods to identify drought conditions and characterize how drought changes over space and time. In this paper, we applied the surface energy balance model, SEBS (Surface Energy Balance System), for the period 2001–2018, to estimate evapotranspiration and other energy fluxes over the dehesa area of the Iberian Peninsula, with a monthly temporal resolution and 0.05 pixel size. A satisfactory agreement was found between the fluxes modeled and the measurements obtained for 3 years by two flux towers located over representative sites (RMSD = 21 W m−2 and R2=0.76, on average, for all energy fluxes and both sites). The estimations of the convective fluxes (LE and H) showed higher deviations, with RMSD = 26 W m−2 on average, than Rn and G, with RMSD = 15 W m−2. At both sites, annual evapotranspiration (ET) was very close to total precipitation, with the exception of a few wet years in which intense precipitation events that produced high runoff were observed. The analysis of the anomalies of the ratio of ET to reference ET (ETo) was used as an indicator of agricultural drought on monthly and annual scales. The hydrological years 2004/2005 and 2011/2012 stood out for their negative values. The first one was the most severe of the series, with the highest impact observed on vegetation coverage and grain production. On a monthly scale, this event was also the longest and most intense, with peak negative values in January–February and April–May 2005, explaining its great impact on cereal production (up to 45 % reduction). During the drier events, the changes in the grasslands' and oak trees' ground cover allowed for a separate analysis of the strategies adopted by the two strata to cope with water stress. These results indicate that the drought events characterized for the period did not cause any permanent damage to the vegetation of dehesa systems. The approach tested has proven useful for providing insight into the characteristics of drought events over this ecosystem and will be helpful to identify areas of interest for future studies at finer resolutions.

1 Introduction

Drought, which is a devastating natural hazard and is globally widespread, has complex consequences across spatiotemporal scales and sectors. Unlike other disasters, it is still a challenge to define, detect and quantify droughts (Sheffield and Wood, 2011), impeding most prevention and mitigation actions. When droughts affect savannas, the two canopies of this ecosystem, grasslands and trees/shrubs, suffer from different stresses: (i) the pasture production is reduced or lost, with a direct economic consequence resulting from the need to supplement animal feeding and, in more severe situations, the death or premature sale of animals; (ii) the decline and dieback of trees affect the ecosystem structure, jeopardizing the long-term conservation of the system (Fenshan and Holman, 1999). Traditional agropastoral systems in arid and semi-arid areas have developed strategies to cope with drought, such as diversifying crops and livestock, adding different animal species and breeds or fluctuating herd sizes (Hazell et al., 2001). More recently, insurance services have started to offer insurance for damage to pasture production caused by water stress, providing farmers with a means to recover after a disaster. However, the slow onset of drought, the large extension of savanna areas and their complex canopy structure introduce additional difficulties to the challenge of monitoring drought and assessing its adverse effects.

The increasing availability of global meteorological data and new remote-sensing products, with advanced processing services and free and open data, offers an opportunity to characterize drought objectively and to extend its analysis in space and time. Many indicators of drought using remote-sensing inputs have been developed in the last decades (Wardlow et al., 2012). Surface energy balance models (SEBMs) provide a physically based rationale to combine the most often used remote-sensing retrievals for drought monitoring: vegetation indices (VIs) and land surface temperature (LST). The VIs provide information about the amount and condition of the vegetation (Jackson and Huete, 1991), while the LST describes the state of the surface and the partitioning of the available energy into sensible heat (H) and latent heat (LE) or evapotranspiration (ET) (Kustas and Norman, 1996). SEBMs have been used to provide ET estimations over agriculture (Anderson et al., 2015; Allen et al., 2011; Cammalleri et al., 2012; Andreu et al., 2015; Gonzalez-Dugo et al., 2009, 2012) and agroforestry systems (Andreu, 2018a, b; Guzinski et al., 2018; Carpintero et al., 2016). In particular, the SEBS (Surface Energy Balance System) model (Su, 2002) presents a good compromise between the detailed parametrization of the turbulent heat fluxes for different states of the land surface and the minimization of the input requirements of the model without the need of local calibration. The evapotranspiration of a canopy is a suitable indicator of its water status and a good measurement of the impact of water shortage on vegetation and the functioning of the ecosystem. Evapotranspiration and soil moisture anomalies have been widely used for the spatially distributed monitoring of agricultural drought (Anderson et al., 2016; Cammalleri el al., 2015; Sheffield et al., 2004). These anomalies underline the abnormally dry conditions when compared to the usual state of an ecosystem, derived from historical data. Evapotranspiration anomalies were used here to assess drought and vegetation water stress in the holm oak savanna area of the Iberian Peninsula over a period of 17 years.

The Mediterranean oak savanna, called dehesa in Spain and montado in Portugal, is the most extensive and representative agroforestry system in Europe, with an area of more than 3×106 ha in the Iberian Peninsula (Moreno and Pulido, 2009). It is a man-made ecosystem that maintains a fragile balance between its multiple uses (livestock, cereal crops, cork, hunting, etc.) and the conservation of its natural resources. The dehesa's diversity of habitats, giving refuge to a large number of species (Díaz et al., 1997), is especially recognized, and it is listed as having community-wide interest in the EU Habitats Directive (92/43/EEC). It is a water-controlled system, with its productivity directly dependent on water availability. Mediterranean oaks can minimize the effects of water scarcity through a combination of physiological mechanisms that occur over a range of timescales (Rambal, 1993). However, an additional problem to the recurrent water scarcity is the identification of low soil water content as an initiating factor involved in the severe oak decline affecting a large area of dehesa since the early 1980s (Sánchez et al., 2002). Drought events impede the growth of Quercus ilex seedlings and increase their susceptibility to Phytophthora cinnamomi (Corcobado et al., 2014), the main biotic factor responsible for this decline (Sánchez et al., 2002).

Similarly to other savanna ecosystems, the different components of dehesa structure (sparse tall vegetation, large areas of grasses, shrubs and bare soil) contribute differently to the turbulent exchange and radiative transfer, hindering its modeling, especially when compared with more homogeneous landscapes. In addition, these vegetation layers differ in phenology, physiology and function: while most trees are evergreen and have access to deep sources of water all year, the herbaceous layer only taps water from the first centimeters of soil and dries up during summer. The combination of the different functioning and characteristics of the system components affects the exchange of sensible and latent heat flux, resulting in a high spatial and temporal flux variability difficult to account for in model parametrization and algorithms. This structure appears to play an important role in savannas' resilience, making the system an efficient convector of sensible heat and keeping the canopy surface temperature inside the adequate range for survival (Baldocchi et al., 2004).

In this work, a surface energy balance model, SEBS (Surface Energy Balance System; Chen et al., 2013; Su, 2002), has been applied to estimate evapotranspiration and other energy fluxes from 2001 to 2018 over the dehesa areas of Spain and Portugal. The first objective was to validate the energy fluxes produced by this model over the dehesa landscape. The second was to analyze the anomalies of the ratio of ET to reference ET as an indicator of agricultural drought in this environment at monthly and annual scale and use it to characterize the main drought events occurring in this period in space and time.

Figure 1Distribution of oak savanna area in the Iberian Peninsula. Location of Sta.Clo (Santa Clotilde) and ES-LMa (Las Majadas) validation sites and pictures of both eddy covariance flux towers.

2 Data and methodology

The study was conducted over the oak savanna area of the Iberian Peninsula (Fig. 1) using data from January 2001 to August 2018. This ecosystem covered 3.12×106 ha in 2006 according to the European CORINE Land Cover inventory (CLC2006 100 m – version 12/2009; (last access: 5 February 2021) clc-2006-raster-4). The area has remained fairly stable during the study period, with changes of less than 1.5 % between CLC2006 and the previous and posterior inventories, in 2000 and 2012.

2.1 SEBS model description

A revised version of the surface energy balance system model known as SEBS (Su, 2002) was used to estimate land heat fluxes, integrating remote-sensing and meteorological forcing data. A brief description of the model is presented below (for further discussion, see Su, 2002, and Chen et al., 2013). The latent heat flux (LE) was computed as a residual of the surface energy balance equation:

(1) LE = R n - G - H ,

where Rn is the net radiation, G is the soil heat flux and H is the turbulent sensible heat flux. The net radiation is calculated using the following equation:

(2) R n = 1 - α SW d + ε LW d - ε σ LST 4 ,

where α is broadband albedo, SWd the downward shortwave radiation, LWd the downward longwave radiation, ε the land surface emissivity, σ the Stefan–Boltzmann constant and LST the land surface temperature.

The soil heat flux is derived from its ratio to the net radiation (Γ) using Eq. (3):

(3) G = R n Γ c + 1 - f c Γ S - Γ c .

This ratio is assumed to be equal to 0.05 (Monteith, 1973) for surfaces with fully covered vegetation (Γc) and 0.315 for bare soils (Γs) (Kustas and Daughtry, 1990). The green canopy cover, fc, is determined using the normalized difference vegetation index (NDVI) in Eq. (7).

Using Eqs. (1) to (3) and energy balance considerations for limiting cases, the following reductions can be applied: (i) under the dry limit (Eq. 4), the evapotranspiration, λEdry, is assumed to become zero due to the limitation of soil moisture and the sensible heat flux, Hdry, is at its maximum,

(4) λ E dry = R n - G - H dry 0 .

(ii) Under the wet limit (Eqs. 5 and 6), the evaporation takes place at a potential rate, λEwet, only limited by the available energy at the given surface and atmospheric conditions. The sensible heat takes its minimum value, Hwet, with the internal resistance of the Penman–Monteith combination equation in the form written by Menenti (1984), ri≡0, by definition.


where ρ is the density of air, Cp the specific heat at constant pressure, e and es the actual and saturation vapor pressure respectively, γ the psychrometric constant, Δ the rate of change of saturation vapor pressure with temperature and rew the external or aerodynamic resistance. The sensible heat is computed according to the Monin–Obukhov similarity theory and limited by the dry and wet conditions. A complete description of the model and the use of the dry and wet limits can be found in Su (2002).

2.2 Model parametrization and dataset preparation

For the application of SEBS over the dehesa area, two surface variables, fc and the height of the canopy (hc), have been adapted to the specific characteristics of this ecosystem. The green canopy cover (fc) and leaf area index (L) were calculated using the following equations (adapted from Choudhury et al., 1994):


where NDVImax and NDVImin represent a surface fully covered by vegetation ( 0.94) and completely bare ( 0.15), respectively. The parameter ξ represents the ratio of the canopy extinction coefficient (K) to a leaf angle distribution term (k). k was assumed to be equal to 0.5 for a random distribution of leaves, as the ecosystem contains erectophile grasses and planophile oak tree leaves (Andreu et al., 2019). K adopted a value of 0.8 obtained from experimental data and within the range proposed for NDVI by Baret and Guyot (1991). NDVI data were provided by the MODIS instrument, averaging the 16 d original product to a monthly scale.

The height of the canopy was computed to account for variations in the tree component. This variable is needed for calculating the momentum roughness length and, thus, important for the sensible heat calculation. The tree stratum of the dehesa is quite homogeneous in composition, dominated by mature Quercus ilex sp., and the grassland canopy has a very high variability of low-height herbaceous species. Considering these reasons, the ecosystem structure has been simplified to compute hc in the following way: a constant height of 8 m has been assigned to oak trees, which is multiplied by its ground coverage in each pixel. Oak fc is computed annually using summer NDVI in Eq. (7). During the summer, the grasslands are dry, and the only photosynthetically active vegetation contributing to the NDVI signal is the oak trees. The grassland height is low (< 1 m), affecting the effective canopy height of each pixel less than the trees, and it is also difficult to compute based on monthly vegetation indices given the high species variability. For this reason, the grassland height has been discarded, and only the contribution of trees was considered to compute hc. Thus, a single hc value was used for every month of a year. This simplification of a complex system certainly may contribute to the error of modeled fluxes. However, it was an operative solution considering the scale of this study.

Table 1Input datasets used to calculate the surface energy fluxes over the Iberian Peninsula from 2000 to 2018.

a (last access: 5 February 2021). b (last access: 5 February 2021). c (last access: 5 February 2021).

Download Print Version | Download XLSX

The SEBS model was originally designed for instantaneous applications. Monthly calculations using the same model were demonstrated by Chen et al. (2014). The structure of the model was not changed, and the implementation differed in the input datasets. The model was applied over the entire Iberian Peninsula with a spatial resolution of 0.05 and a monthly input dataset. Satellite and meteorological input datasets are described in Table 1. All datasets were spatially averaged or subdivided to a common resolution of 0.05.

The land surface temperature (LST) was provided by the MODIS instrument, using the monthly mean of the day and night LST product, which provides the most complete coverage. The accuracy of this product, a key variable in SEBMs, was evaluated by Chen et al. (2017), supporting its applicability for climate studies and numerical model evaluation.

Meteorological data were provided by the ERA-Interim, a global atmospheric reanalysis dataset from the European Centre for Medium-Range Weather Forecasts (ECMWF). Monthly means of daily means were produced by ECMWF as the average of the four main synoptic monthly means at 00:00, 06:00, 12:00 and 18:00 UTC. The forecast model, data assimilation method and input datasets used to produce ERA-Interim can be found in Dee et al. (2011) and a description of the product archive in Berrisford et al. (2011).

To analyze model results, the monthly rainfall gridded data of the Climatic Research Unit (CRU) Time-Series (TS) Version 3.21 (Harris et al., 2014), provided by the Global Climate Monitor System (Camarillo-Naranjo et al., 2019), have been averaged over the dehesa area of the Iberian Peninsula.

2.3 Validation sites and model evaluation

Two experimental sites (Fig. 1) with similar flux measurement instrumentation have been used to validate the evapotranspiration and other energy fluxes estimated using the SEBS model. Both eddy covariance towers, named Sta.Clo (Santa Clotilde, Andalusia; 3812 N, 417 W; 736 m a.s.l.) and ES-LMa (Boyal de Majadas del Tiétar, Extremadura; 3956 N, 546 W; 260 m a.s.l.) are located over dehesa-type ecosystems under similar management and a landscape of scattered oak trees with a fractional cover of around 20 %, in southern and southwestern Spain, respectively. The convective fluxes of the systems are measured above the tree height (at 17 m in Sta.Clo and 15 m in ES-LMa), with closure balance errors of 20 % and 14 %, both values being within the range found by other authors (Foken, 2008; Franssen et al., 2010). For ES_LMa the processing of the data corresponded to the procedure standardized by the FLUXNET network (, last access: 5 February 2021). For Sta.Clo, detailed information on the measurements and the processing of the data can be found in Andreu et al. (2018a, b). In this case, the comparison period was selected attending to the quality of the data, and some months (3 of 36) were discarded due to missing information. Soil moisture, precipitation and other complementary measurements of the vegetation (reflectance, L, green canopy cover) were used to characterize the dynamics of the vegetation and the soil water status throughout the year.

The area contributing most to the fluxes measured was estimated by using Schuepp et al. (1990) and varied between 1 and 2 km. These footprints are lower than the pixel size of 5 km used for the application of the SEBS model. However, the homogeneity of the system, with similar tree ground cover fraction and pasture management at several kilometers around the towers, supported the capacity of these sites to serve as a reference for the validation of modeled fluxes. In both cases, the good correspondence between the model input meteorological data at the tower's location and the ground measurements was verified (data not shown).

Monthly rainfall data for the 17 years of the study were provided by the closest weather station to each site, located 3 and 16 km from Sta.Clo and ES-LMa towers, respectively. Both of them are operated by the Spanish Meteorology Agency (AEMET).

Model performance was quantified via the root mean square difference (RMSD) and the coefficient of determination (R2) between the modeled and observed fluxes. In addition, the mean bias error (MBE), computed by taking the difference between predicted and observed fluxes, was used to assess model under- and overestimations.

2.4 Water stress calculations

The relative evapotranspiration is the ratio of actual to potential or reference ET (ET / ETo). It has been used as an indicator of crop water stress (Anderson et al., 2015, 2016), of drought (Anderson et al., 2011) and as a proxy for soil moisture (Su et al., 2003). The same approach is used worldwide in irrigation engineering to compute crop water requirements following FAO (24 and 56) guidelines (Doorenbos and Pruitt, 1977; Allen et al., 1998). The reason for normalizing ET by ETo is to separate the ET signal component responding to soil moisture from variations due to the available energy. Anderson et al. (2011) showed that anomalies in ET / ETo were more strongly correlated with other drought indices as were anomalies in ET for most US climatic divisions, showing strong agreements in the southwest of the country, with a similar climate to the study area. The comparison of both variables anomalies has also been performed here.

Anomalous water stress conditions indicating drought were assessed here with the standardized values of relative ET. FAO56 reference ET (Allen et al., 1998) was selected to estimate the atmospheric evaporative demand (AED), given the difficulties of reproducing the biological control of the transpiration, even at potential rates, of the different types of vegetation conforming this ecosystem.

The vegetation water stress caused by the long dry summers of the Mediterranean climate can be considered to be the “normal” state of the system for several months of the year. To identify unusually dry conditions indicating drought, standard (z) scores of this variable (ET / ETo) for a given month/year have been computed. This standardization procedure assumes that the data follow a normal distribution. Some authors (Sheffield et al., 2004; Cammalleri et al., 2015) have pointed out that soil moisture and the water deficit index derived from it are generally characterized by a skewed distribution and can be statistically better represented using the beta distribution. In this case, the analysis of ET and relative ET monthly histograms (shown in the Supplement) indicated that most months presented an approximately symmetric distribution, with skewness between −0.5 and 0.5 for both variables. Among the months studied, 3 months were moderately skewed, and only 1 month (for ET) and 2 months (for ET / ETo) were slightly above 1, backing up the use of z scores for the standardization of this variable. Annual drought analyses were performed by averaging monthly anomalies.

Drought intensity is defined here in terms of the maximum negative anomaly of relative ET values reached during an event (thus using the standard deviation as a measure of its departure from the mean) and the drought event duration as the successive number of months with negative anomalies. To classify the events occurred during the study period, the following thresholds have been used: severe drought (anomalies <=-1.5); moderate drought (anomalies between −1 and −1.5) and mild drought (anomalies between −1 and 0). These classes are used for both annual and monthly time steps.

Two variables, vegetation coverage (fc) and rain-fed wheat production, have been selected as drought impact indicators. The vegetation condition and the failure of crops are known consequences of a declining soil moisture, and both have been used previously as indicators of drought (Liu and Kogan, 1996; FAO, 1983). Winter cereals are the main cropping system of these areas, in which the low fertility of the soils does not allow for a more intense agricultural use. Its growth cycle is similar to that of the natural grasslands, with both of them escaping drought and coping with the long summer dry season by completing their life cycle before serious soil and plant water deficits develop. Given that no irrigation is provided, the impact of moisture deficits over its yield can be consider an indirect indicator of the impact of drought on dehesa herbaceous vegetation. Annual yield statistics (, last access: 5 February 2021) have been gathered and aggregated for the dehesa area (Fig. 1).

Figure 2Comparison of monthly energy fluxes of latent heat (LE), sensible heat (H), net radiation (Rn) and soil heat flux (G) estimated using the SEBS model at a monthly scale and observed fluxes at each oak savanna site: ES-LMa (LA) for the years 2009–2011 and Sta.Clo (SC) for the years 2015–2017.


Figure 3Evolution of annual rainfall, ET, ETo and ET / ETo at the ES-LMa site (a) and the Sta.Clo site (b) and annual runoff at Sta.Clo watershed from the hydrological years 2001/2002 to 2017/2018.


3 Results and discussion

3.1 Model validation

The comparison of SEBS model estimation of monthly energy fluxes with measurements at the two eddy covariance (EC) towers during a total of 6 years, 2009 to 2011 for ES-LMa and 2015 to 2017 for Sta.Clo, displayed in Fig. 2, generally showed good agreement, with an average root mean square difference (RMSD) of 21 W m−2 and R2 of 0.76, for all energy fluxes and both sites. The estimations of the convective fluxes (LE and H) show higher deviations, with RMSD = 26 W m−2 on average, than Rn and G, with RMSD = 15 W m−2. Model performance at ES-LMa site was, in general, superior to that at Sta.Clo, with all the statistics metrics computed for the comparison (RMSD, MBE and R2) presenting lesser dispersion and slightly lower errors. LE was slightly overestimated at both sites (MBE = 10.3 and 2.8 W m−2 at Sta.Clo and ES-LMa, respectively), which is in agreement with previous applications of the model (Michel et al., 2016). This overestimation was particularly significant for some springtime months at Sta.Clo, when the sensible heat was underestimated by the SEBS model (Chen et al., 2019). It is worth noting than the model forces the closure of the energy balance, and the error in LE can be attributed to the propagation of errors in all the other balance components. However, LE estimations presented a similar or lower RMSD than other applications of the SEBS model (Chen et al., 2014; Vinukollu et al., 2011). In particular, the work by Chen et al. (2014) estimated energy fluxes over China at the same temporal scale and with similar input databases. The comparison with measurements at 11 Chinese flux towers presented results that were very close to the ones obtained by this application. Mean RMSDs for all fluxes were alike (RMSD = 22 W m−2 was reported by Chen et al., 2014), with a marginally better performance for convective fluxes and a poorer one for Rn and G (RMSDs in China were 22 and 24 W m−2 for convective fluxes and, Rn and G, respectively).

Figure 4Annual anomalies of relative evapotranspiration at ES-LMa and Sta.Clo experimental sites estimated using the SEBS model from 2001/2002 to 2017/2018.


Figure 3 presents the evolution of modeled ET and ETo, ET / ETo and measured precipitation from 2001 to 2018, aggregating the hydrological year (between 1 October and 30 September) at the two experimental sites. It can be observed that annual ET variations for the period followed a similar pattern of precipitation at both sites, confirming the predominant control of water availability over the evaporation in these systems. This control is consequently extended to ecosystem productivity, and in most years the water consumption, coupled to biomass production, is close to the total rainfall. Tree density is similar at both sites, and the differences in water consumption between them are explained by variations in annual pasture production, due to differences in water availability and soil properties. Very wet years, and those with average rainfall but intense precipitation events producing an increase in runoff, did not follow this pattern. This can be observed by the runoff recorded at Sta.Clo watershed reservoir (Fig. 3a). The main land use of this small watershed (48.4 km2) is dehesa, but other uses can be found as well, such as olive orchards and field crops.

Annual runoff measurements followed a close relationship (data shown in the Supplement, Fig. S2) with the annual aridity index (Budyko, 1974) estimated at Sta.Clo following Arora (2002) as the ratio between potential evaporation and annual precipitation. On average, we found aridity indices of above 1 at both sites, indicating dry regions where the evaporative demand cannot be met by precipitation. In this case, AED was computed using Penman–Monteith for comparison purposes. Sta.Clo site is noticeably less arid than ES-LMa, with an aridity index equal to 2.9 and 3.75 on average for the 17 hydrological years at Sta.Clo and ES-LMa, respectively, with both of them falling under the category of a semi-arid climate regime (Ponce et al., 2000). The two sites presented similar annual ETo values for the period (Fig. 3), but annual precipitation was around 200 mm higher, on average, at Sta.Clo, with a higher and more variable ET / ETo throughout the years. What can also be observed in Fig. 3 is the complementary relationship between actual and reference evapotranspiration at this temporal scale, with the sum of annual ET and ETo approaching a constant value at both sites, confirming the complementary hypothesis (Bouchet, 1963; Morton, 1975; Brutsaert and Stricker, 1979).

3.2 Annual drought monitoring and impact assessment

Drought was characterized on an annual scale over the experimental sites and the whole area of the dehesa of the Iberian Peninsula using the relative evaporation anomalies. Figure 4 presents their evolution for the two sites throughout the study period. A clear similarity can be observed in the main negative anomalies, which identify the most severe droughts during the years 2004/2005 and 2011/2012 at both sites, despite the differences in aridity and the distance (Fig. 1) between them, indicating the extended area and intensity of both events. Differences are more evident in the case of the mild droughts, occurring at both sites but with different intensities during two periods, 2007 to 2009 and 2016 to 2018.

Figure 5Evolution from 2001/2002 to 2017/2018 of annual anomalies of relative evapotranspiration, energy balance components, air and surface temperature, vegetation ground fraction cover and rainfed wheat yield, aggregated for the whole oak savanna area of the Iberian Peninsula.


When the whole dehesa area is considered (Figs. 5 and 6), a more complete view of the general intensity, impact and spatial distribution of those dry periods can be obtained. Figure 5 aggregates, for the total dehesa area, the evolution of the relative ET anomalies, together with the exchanges of energy between the surface and the atmosphere, the green canopy cover and the production of rainfed wheat. The last two variables were selected as indicators of the impact of water scarcity on the system.

Figure 6Spatial distribution of annual anomalies of relative evapotranspiration for the oak savanna area of the Iberian Peninsula from 2001/2002 to 2017/2018, the average ET / ETo for the period and its standard deviation (SD).

The 2 severely dry years identified at the experimental sites were the driest ones for the entire dehesa area, with 2004/2005 standing out as the most severe event of the time series. None of them lasted more than 1 year. For these 2 dry years, a reduction in the latent heat can be observed when compared to the complete series, producing a swap with the sensible heat in the second position in magnitude of the energy balance components. A rise in the surface temperature, increasing the difference with the air temperature, is also observed for those dry years. The order of severity in dryness, established by the magnitude of negative values of ET / ETo anomalies, is also observed in their impacts over the system (Fig. 6). In 2004/2005, the wheat production in the area was reduced by almost half of the average (45 %) for the period analyzed, and the vegetation groundcover fraction fell by 20 % compared to the average of the same period. This severe drought affected the entire Iberian Peninsula, with Spanish and Portuguese cereal and hydroelectricity production decreasing by 40 % and 60 % with respect to the average (Garcia-Herrera et al., 2007) and a 10 % reduction in total EU cereal yields (UNEP, 2006). The event during 2011/2012 was among the largest and most severe ones in Europe for the 18-year simulation period analyzed by Cammalleri et al. (2015), contributing to a global decline in grain production.

Figure 7(a) Monthly evolution of evapotranspiration anomalies (blue line) of the oak savanna area of the Iberian Peninsula from January 2001 to August 2018, with negative values indicating drier than normal conditions (depicted in red) and green canopy cover (green line). The dashed green lines connect the annual maximum and minimum values of fc. (b) Monthly evolution of rainfall, ETo and ET in the same region and time interval.


Figure 6 shows maps of ET / ETo anomalies in Iberia for the 17 years of the study, highlighting the dehesa area of interest in this work. The spatial variability of these anomalies for most years is significant, although prevalently dry and wet years can be distinguished. In 2004/2005 and 2011/2012, the drought was severe and affected most of the area of interest, as the aggregated values of Fig. 5 also point out. In 2008/2009, the water stress was milder in the western area, as can be observed in Fig. 6, than at the experimental site of Sta.Clo (Fig. 4) located in this part of the region. The recovery of the vegetation water status, in most areas, was achieved the year following dry ones.

3.3 Monthly drought analysis

The monthly evolution of relative evapotranspiration anomalies is displayed in Fig. 7a, with negative values indicating water stress conditions highlighted in red. Absolute ET and ETo values, used to calculate these anomalies, are shown in Fig. 7b together with monthly rainfall for the period. One can observe the alternation of complementary and parallel characteristics of ET and ETo throughout the year. The longest complementary period indicating water-limited ET conditions, starting in May for most of the years, is confirmed by the decreasing trend in rainfall starting in that month. At the end of the summer when the first rains arrive, the trend of ET and ETo changes, producing a secondary peak in ET, much weaker than the one earlier in the year, that lasts until the energy-limited parallel phase starts in November. Both variables follow a concurrent rise from January until the soil water deficit limits ET again.

The annual fluctuations of the green canopy cover (thick green line in Fig. 7a) followed the expected seasonality of Mediterranean vegetation, corresponding to the dynamics of ET and ETo changes. The maximum coverage (March and April) corresponds to the peak of grassland production (and ET although with different shape), and the minimum appears during the dry summer, only endured by the oak trees. In some years, the growing season presents a bimodal shape, with an initial peak produced by autumn pastures, which is also reflected in ET values. It can be observed mostly in wet years (e.g., 2003, 2007, 2011), with the vegetation growth following a pattern that can be related to the soil water availability, represented here by the ET / ETo anomalies.

The duration and intensity of each drought event help to explain the response of the vegetation during these periods. In this sense, the two main drought events identified on an annual scale (2004/2005 and 2011/2012) presented drier than normal conditions during the whole or most of the year. The first event was longer (16 months in the first case, prolonging the drought to the beginning of the following year) and with higher negative values than the second one, of an 11-month duration, explaining the greater impacts detected on the vegetation and cereal yield. Other dry periods, in 2009, 2017 and 2018, presented consecutive negative anomalies for 10 to 11 months, but, in some cases, the non-homogeneous distribution of the drought, observed in Fig. 6, may have undermined the impact analysis on this aggregated spatial scale. In terms of impact assessment, the time of the year with peak negative anomalies is important, with springtime events producing greater impacts (e.g. in 2004/2005 the highest negative values corresponded to January, February, April and May 2005).

During the dry years, the annual vegetation growth pattern varies with respect to the typical one, depending on the duration and severity of drought events. The dynamics of the vegetation in this system allows for a separate analysis of the effect of water scarcity over trees and pastures. The dashed green lines (Fig. 7a) show the changes in annual maximum and minimum values of fc, with the maximum ones mostly expressing the impact on pasture, and the changes in the minimum ones representing only the impact over the tree canopy. The decreases in pasture fc are more pronounced than changes in oaks fc, as grasslands are more abundant, and their roots are mostly located in the first centimeters of soil. On the contrary, the rooting system of the oak tree is in fact adapted to the regular dry periods of the Mediterranean climate, exploring a large volume of soil that can reach maximum values of around 5 m in depth and 30 m in horizontal extension (Moreno et al., 2005). The small decreases, observed in oaks fc in Fig. 7a during dry years, generally recovered within 1 or 2 years. This response of the tree leaf area is associated with low-frequency oscillations, such as annual rainfall (Poole and Milles, 1981). This is also supported by the variance observed in fc that can be explained by the anomalies of relative evapotranspiration of previous months. During the spring, the highest correlation coefficients are obtained for the previous 2 or 3 months (e.g., average fc for the peak month, April, is correlated with average anomalies from February to April with an R2 equal to 0.76 and with anomalies of the previous year with an R2=0.52). However, during the summer, the coverage of the vegetation can be better explained by what has happened during the previous year (e.g., R2 is equal to 0.39 for average August fc and the anomalies of the two previous months and 0.64 for the anomalies of the year), suggesting that those values of fc might be linked to processes occurring at different timescales.

Figure 8Comparison of monthly negative anomalies of ET, ET / ETo and fc for the entire oak savanna area of the Iberian Peninsula from January 2001 to August 2018.


A more detailed analysis is required, but these results support the conclusion that the drought events characterized for this period did not cause any permanent damage to the vegetation, considering both the grasslands and the oak trees.

Similar results can be derived from the analysis of ET anomalies. Figure 8 presents a comparison of monthly anomalies of ET, ET / ETo and fc. The anomalies of ET and ET / ETo showed a high similarity for the conditions of the study, with correlations of R2=0.76 at monthly scale and R2=0.82 at seasonal scale (results presented in Figs. S3 and S4). It suggests that ET anomalies could be an option to monitor drought in dehesa areas. Nevertheless, the computation of ETo does not require additional variables than those already used by the energy balance models, with quite a straightforward computation. Once actual ET is estimated, the computation of ET / ETo takes very little effort and adds some confidence to the focus on the soil moisture signal. Regarding the evaluation of fc anomalies, it can be derived that the drought events identified using this variable would have been the same as using ET or ET / ETo but with different intensities and duration. The main differences can be found during the cold winter months when the vegetation is largely dormant. In these cases, the anomalies of fc, similar to the performance of other indices based on vegetation, such as the Vegetation Condition Index (VCI; Heim, 2002) have a limited utility. The results are more comparable and could be more useful during the growing season.

4 Conclusions

The SEBS model was used to estimate monthly energy fluxes over the dehesa area of the Iberian Peninsula from January 2001 to August 2018. There was a satisfactory agreement between modeled fluxes and measurements obtained for 3 years over two sites that are representative of the ecosystem.

At both sites annual ET was very close to total precipitation, with the exception of a few wet years and those in which intense precipitation events producing a high runoff were observed. Average aridity indices for the 17 hydrological years of 2.9 and 3.75 were computed at Sta.CLo and ES_LMa, respectively, indicating that their evaporative demand cannot be met by annual precipitation of these sites.

Drought has been characterized on an annual and monthly scale over the experimental sites and the whole area of dehesa of the Iberian Peninsula using relative evaporation anomalies (ET / ETo). At the annual scale, the negative anomalies of 2 years, 2004/2005 and 2011/2012, stood out during the study period at the experimental sites and the entire dehesa area. However, a recovery of average values is observed in the years following the dry ones, indicating the absence of prolonged droughts for the period. Maps of ET / ETo anomalies showed that most of the dehesa area was affected in those dry years. These maps complemented the averaged data, providing spatial information about regional impacts that could be useful for a more detailed analysis.

On the monthly scale, the drought event of 2004/05 is confirmed as being the longest and the most intense event, with 16 consecutive months of negative anomalies (from October 2004 to January 2006). Peak negative values in January–February and April–May 2005 explain the important impact on cereal production. The dynamics of the vegetation strata on a monthly scale allows for a separate assessment of water stress impacts on oaks and pastures. The different behavior observed in vegetation ground cover during the drier events in months with a preponderant presence of grasslands, compared with months in which only oaks were active, is consistent with the different strategies adopted by the two strata to cope with water stress. In addition, the correlation of monthly vegetation fractional coverage with previous short or medium-term anomalies (from 2 months to 1 year) suggest that those values might be linked to processes occurring on a different timescale, depending on whether the grassland or the tree is the predominant vegetation.

These results back up the conclusion that the drought events characterized for this period did not cause permanent damage to the vegetation of dehesa systems, considering both the grasslands and the oak trees. The approach proved useful for providing insights into the characteristics of drought events over this ecosystem and for defining and identifying areas of interest for future studies at finer resolutions.

Code and data availability

The SEBS code is available to download from the GitHub repository (, last access: 5 February 2021; Chen, 2020). Validation data of the ES-LMa site are available from the European Fluxes Database Cluster (, last access: 5 February 2021; Carrara, 2021), and data of the Sta.Clo site may be distributed on request to the principal investigator of the Sta. Clotilde experimental site (María P. González-Dugo, IFAPA,


The supplement related to this article is available online at:

Author contributions

MPGD conceived the original idea, analyzed the data and took the lead in writing the manuscript. XC and ZS designed the model and the computational framework and contributed to the interpretation of the data. MPGD and XC collected the input data and performed the numerical calculations. AA, EC, PJGG and AC collected and analyzed the validation data and reviewed the paper. All authors provided critical feedback and helped to shape the manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Data acquisition and modelling of hydrological, hydrogeological and ecohydrological processes in arid and semi-arid regions”. It is not associated with a conference.


We would like to thank the owners and workers of the Santa Clotilde experimental site, as well as the group managing the experimental site of Las Majadas for the eddy covariance measurements and the additional data. We also thank the anonymous reviewers, whose comments have improved the manuscript.

Financial support

This research has been supported by the OECD Cooperative Research Programme: Biological Resource Management for Sustainable Agricultural Systems (grant no. JA00084693) and the projects PP.PEI.IDF201601.16 and PP.PEI.IDF2019.004, 80 % cofunded by the European Regional Development Fund, AOP 2014–2020. Additional support was provided by RTA2014-00063 C04-02 INIA-FEDER and PID2019-107693RR-C22 projects (MCIU/AEI/FEDER, UE). XC was supported by the National Natural Science Foundation of China (41975009) and AA by EU Horizon 2020 Marie Skłodowska-Curie Action grant agreement no. 703978.

Review statement

This paper was edited by Harrie-Jan Hendricks Franssen and reviewed by three anonymous referees.


Allen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop evapotranspiration: Guidelines for computing crop requirements, Irrigation and Drainage Paper No. 56, 300 pp., FAO, Rome, Italy, 1998. 

Allen, R. G., Irmak, A., Trezza, R., Hendrickx, J. M. H., Bastiaanssen, W., and Kjaersgaard, J.: Satellite-based ET estimation in agriculture using SEBAL and METRIC, Hydrol. Proc., 25, 4011–4027,, 2011. 

Anderson, M. C., Hain, C. R., Wardlow, B., Pimstein, A., Mecikalski, J. R., and Kustas, W. P.: Evaluation of Drought Indices Based on Thermal Remote Sensing of Evapotranspiration over the Continental United States, J. Climate, 24, 2025–2044,, 2011. 

Anderson, M. C., Zolin, C. A., Hain, C. R., Semmens, K., Yilmaz, M. T., and Gao, F.: Comparison of satellite-derived LAI and precipitation anomalies over Brazil with a thermal infrared-based Evaporative Stress Index for 2003–2013, J. Hydrol., 526, 287–302,, 2015. 

Anderson, M. C., Zolin, C. A., Sentelhas, P. C., Hain, C. R., Semmens, K., Tugrul Yilmaz, M., Gao, F., Otkin, J. A., and Tetrault, R.: The Evaporative Stress Index as an indicator of agricultural drought in Brazil: An assessment based on crop yield impacts, Remote Sens. Environ., 174, 82–99,, 2016. 

Andreu, A., Timmermans, W. J., Skokovic, D., and González-Dugo, M. P.: Influence of component temperature derivation from dual angle thermal infrared observations on TSEB flux estimates over an irrigated vineyard, Acta Geophys., 63, 1540–1570,, 2015. 

Andreu, A., Kustas, W. P., Polo, M. J., Carrara, A., and González-Dugo, M. P.: Modeling Surface Energy Fluxes over a Dehesa (Oak Savanna) Ecosystem Using a Thermal Based Two-Source Energy Balance Model (TSEB) I., Remote Sens., 10, 567,, 2018a. 

Andreu, A., Kustas, W. P., Polo, M. J., Carrara, A., and González-Dugo, M. P.: Modeling Surface Energy Fluxes over a Dehesa (Oak Savanna) Ecosystem Using a Thermal Based Two Source Energy Balance Model (TSEB) II – Integration of Remote Sensing Medium and Low Spatial Resolution Satellite Images, Remote Sens., 10, 558,, 2018b. 

Andreu, A., Dube, T., Nieto, H., Mudau, A. E., González-Dugo, M. P., Guzinski, R., and Hülsmann, S.: Remote sensing of water use and water stress in the African savanna ecosystem at local scale – Development and validation of a monitoring tool, Phys. Chem. Earth, 112, 154–164,, 2019. 

Arora, V. K.: The use of the aridity index to assess climate change effect on annual runoff, J. Hidrol., 265, 164–177, 2002. 

Baldocchi, D. D., Xu, L., and Kiang, N.: How plant functional-type, weather, seasonal drought, and soil physical properties alter water and energy fluxes of an oak-grass savanna and an annual grassland, Agr. Forest Meteorol., 123, 13–39,, 2004. 

Baret, F. and Guyot, G.: Potentials and limits of vegetation indices for LAI and APAR assessments, Remote Sens. Environ., 35, 161–173, 1991. 

Berrisford, P., Dee, D., Brugge, R., Fielding, K., Fuentes, M., Kallberg, P., Kobayashi, S., Uppala, S., and Simmons, S.: The ERA-Interim Archive, ERA Report Series No. 1 Version 2.0, European Centre for Medium Range Weather Forecasts, Reading, U.K., availabe at: (last access: 10 December 2019), 2011. 

Bouchet, R. J.: Evapotranspiration reelle, evapotranspiration potentielle, et production agricole, Annales Agronomiques, 14, 743–824, 1963. 

Brutsaert, W. and Stricker, H.: An advection-aridity approach to estimate actual regional evapotranspiration, Water Resour. Res., 15, 443–450, 1979. 

Budyko, M. I.: Climate and life, Academic Press, Orlando, FL, 508 pp., 1974. 

Camarillo-Naranjo, J. M., Álvarez-Francoso, J. I., Limones-Rodríguez, N., Pita-López, M. F., and Aguilar-Alba, M.: The Global Climate Monitor System: From Climate Data-Handling to Knowledge Dissemination, Int. J. Digital Earth, 12, 394–414,, 2019. 

Cammalleri, C., Anderson, M. C., Ciraolo, G., D'Urso, G., Kustas, W. P., La Loggia, G., and Minacapilli, M.: Applications of a remote sensing-based two-source energy balance algorithm for mapping surface fluxes without in situ air temperature observations, Remote Sens. Environ., 124, 502–515,, 2012. 

Cammalleri, C., Micale, F., and Vogt, J.: A novel soil moisture-based drought severity index (DSI) combining water deficit magnitude and frequency, Hydrol. Proc., 30, 289–301,, 2015. 

Carrara, A.: Las Majadas del Tietar flux data, European Fluxes Database Cluster, available at:, last access: 5 February 2021. 

Carpintero, E., Gonzalez-Dugo, M. P., Hain, H., Nieto, H., Gao, F., Andreu, A., and Kustas, W. P.: Continuous evapotranspiration monitoring and water stress at watershed scale in a Mediterranean oak savanna, Proceeding SPIE 9998, Remote Sensing for Agriculture, Ecosystems, and Hydrology XVIII, 99980N, Edinburg, United Kingdom,, 2016. 

Chen, X.: SEBS_Spain, available at:, last access: 5 February 2021, 2020. 

Chen, X., Su, Z., Ma, Y., Yang, K., Wen, J., and Zhang, Y.: An Improvement of Roughness Height Parameterization of the Surface Energy Balance System (SEBS) over the Tibetan Plateau, J. Appl. Meteorol. Clim., 52, 607–622,, 2013. 

Chen, X., Su, Z., Ma, Y., Liu, S., Yu, Q., and Xu, Z.: Development of a 10-year (2001–2010) 0.1 data set of land-surface energy balance for mainland China, Atmos. Chem. Phys., 14, 13097–13117,, 2014. 

Chen, X., Su, Z., Ma, Y., Cleverly, J., and Liddell, M.: An accurate estimate of monthly mean land surface temperatures from MODIS clear-sky retrievals, J. Hydrometeorol., 18, 2827–2847, 2017. 

Chen, X., Su, Z., Ma, Y., and Middleton, E. M.: Optimization of a remote sensing energy balance method over different canopy applied at global scale, Agr. Forest Meteorol., 279, 107633,, 2019. 

Choudhury, B. J., Ahmed, N. U., Idso, S. B., Reginato, R. J., and Daughtry, C. R. S.: Relations between evaporation coefficients and vegetation indices studied by model simulations, Remote Sens. Environ., 50, 1–17, 1994. 

Corcobado, T., Cubera, E., Juárez, E., Moreno, G., and Solla, A.: Drought events determine performance of Quercus ilex seedlings and increase their susceptibility to Phytophthora cinnamomi, Agr. Forest Meteorol., 192, 1–8,, 2014. 

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597,, 2011. 

Díaz, M., Campos, P., and Pulido, F. J.: The Spanish dehesas: a diversity in land-use and wildlife, in: Farming and birds in Europe: The Common Agricultural Policy and its implications for bird conservation, edited by: Pain, D. and Pienkowski, M., Academic Press, London, 178–209, 1997. 

Doorenbos, J. and Pruitt, W. O.: Crop water requirement, FAO Irrigation and Drainage, Paper N. 24, FAO, Rome, 1977, Food and Agriculture Organization: Guidelines: Land evaluation for Rainfed Agriculture, FAO Soils Bulletin 52, Rome, 156 pp., 1983. 

Fenshan, R. J. and Holman, J. E.: Temporal and spatial patterns in drought-related tree dieback in Australian savanna, J. Appl. Ecol., 36, 1035–1050, 1999. 

Foken, T.: The energy balance closure problem: An overview, Ecol. Appl., 18, 1351–1367,, 2008. 

Franssen, H. J. H., Stöckli, R., Lehner, I., Rotenberg, E., and Seneviratne, S. I.: Energy balance closure of eddy-covariance data: A multisite analysis for European FLUXNET stations, Agr. Forest Meteorol., 150, 1553–1567,, 2010. 

García-Herrera, R., Paredes, D., Trigo, R. M., Franco-Trigo, I., Hernández, E., Barriopedro, D., and Mendes, M.: The outstanding 2004/05 drought in the Iberian Peninsula: associated atmospheric circulation, J. Hydrometeorol., 8, 483–498,, 2007. 

Gonzalez-Dugo, M. P., Neale, C. M. U., Mateos, L., Kustas, W. P., Prueger, J. H., Anderson, M. C., and Li, F.: A comparison of operational remote sensing-based models for estimating crop evapotranspiration, Agr. Forest Meteorol., 149, 1843–1853,, 2009. 

González-Dugo, M. P., González-Piqueras, J., Campos, I., Andreu, A., Balbotin, C., and Calera, A.: Evapotranspiration monitoring in vineyard using satellite-based thermal remote sensing, in: Remote Sensing for Agriculture, Ecosystems, and Hydrology XIV, edited by: Neale, C. M. U. and Maltese, A., Proc. of SPIE., 8531, 85310N-85310N-8,, 2012. 

Guzinski, R., Nieto, H., El-Madany, T., Migliavacca, M., and Carrara, A.: Validation of Fine Resolution Land-Surface Energy Fluxes Derived with Combined Sentinel-2 and Sentinel-3 Observations, IGARSS 2018 IEEE Intl. Geosci. Rem. Sens. Symposium, Valencia, 8711–8714,, 2018. 

Harris, I. P. D. J., Jones, P. D., Osborn, T. J., and Lister, D. H.: Updated high-resolution grids of monthly climatic observations–the CRU TS3. 10 Dataset, Int. J. Climatol., 34, 623–642,, 2014. 

Hazell, P. B. R., Oram, P., and Chaherli, N.: Managing droughts in the low-rainfall areas of the Middle East and North Africa, EPTD Discussion Paper, 78, International Food Policy Research Institute (IFPRI), 45 pp., available at:, last access: 5 February 2021, 2001. 

Heim, R. R.: A review of twentieth-century drought indices Used in the United States, B. Am. Meteorol. Soc., 83, 1149–1166,, 2002. 

Jackson, R. D. and Huete, A. R.: Interpreting vegetation indices, Pre. Vet. Med., 11, 185–200,, 1991. 

Kustas, W. P. and Daughtry, C. S. T.: Estimation of the soil heat flux/net radiation ratio from spectral data, Agr. Forest Meteorol., 49, 205–223, 1990. 

Kustas, W. P. and Norman, J. M.: Use of remote sensing for evapotranspiration monitoring over land surfaces, Hydrol. Sci. J., 41, 495–516,, 1996. 

Liu, W. T. and Kogan, F. N.: Monitoring regional drought using the vegetation condition index, Int. J. Remote Sens., 17, 2761–2782, 1996. 

Menenti, M.: Physical Aspects of and Determination of Evaporation in Deserts Applying Remote Sensing Techniques, Report 10 (special issue), Institute for Land and Water Management Research (ICW): The Netherlands, 202 pp., 1984. 

Michel, D., Jiménez, C., Miralles, D. G., Jung, M., Hirschi, M., Ershadi, A., Martens, B., McCabe, M. F., Fisher, J. B., Mu, Q., Seneviratne, S. I., Wood, E. F., and Fernández-Prieto, D.: The WACMOS-ET project – Part 1: Tower-scale evaluation of four remote-sensing-based evapotranspiration algorithms, Hydrol. Earth Syst. Sci., 20, 803–822,, 2016. 

Monteith, J. L.: Principles of Environmental Physics, Edward Arnold Press., London, 291 pp., 1973. 

Moreno, G. and Pulido, F. J.: The Functioning, Management and Persistence of Dehesas, in: Agroforesty in Europe, edited by: Rigueiro-Rodróguez, A., McAdam, J., and Mosquera-Losada, M. R., Springer Link, 127–160,, 2009. 

Moreno, G., Obrador, J. J., Cubera, E., and Dupra, C.: Fine root distribution in Dehesas of Central-Western Spain, Plant Soil, 277, 153–162,, 2005.  

Morton, F. I.: Estimating evaporation and transpiration from climatological observations, J. Appl. Meteorol., 14, 488–497, 1975. 

Ponce, V. M., Pandey, R. P., and Ercan, S.: Characterization of drought across the climate spectrum, J. Hydrol. Eng., 5, 222–224,, 2000. 

Poole, D. K. and Miller P. C.: The Distribution of Plant Water Stress and Vegetation Characteristics in Southern California Chaparral, The American Midland Naturalist, 105, 32–43, 1981. 

Rambal, S.: The differential role of mechanisms for drought resistance in a Mediterranean evergreen shrub: a simulation approach, Plant. Cell. Environ., 16, 35–44,, 1993. 

Sánchez, M. E., Caetano, P., Ferraz, J., and Trapero, A.: Phytophthora disease of Quercus ilex in south-western Spain, Forest Pathol., 32, 5–18,, 2002. 

Schuepp, P. H., Leclerc, M. Y., MacPherson, J. I., and Desjardins, R. L.: Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation, Bound.-Lay. Meteorol., 50, 355–373,, 1990. 

Sheffield, J. and Wood, E. F.: Drought: Past problems and future scenarios, Earthscan London, UK, 192 pp., 2011. 

Sheffield, J., Goteti, G., Wen, F., and Wood, E. F.: A Simulated Soil Moisture Based Drought Analysis for the United States, J. Geophys. Res., 109, D24108,, 2004. 

Su, Z.: The Surface Energy Balance System (SEBS) for estimation of turbulent heat fluxes, Hydrol. Earth Syst. Sci., 6, 85–100,, 2002. 

Su, Z., Yacob, A., Wen, J., Roerink, G., He, Y., Gao, B., Boogaard, H., and van Diepen, C.: Assessing relative soil moisture with remote sensing data: theory and experimental validation, Phys. Chem. Earth, 28, 89–101, 2003. 

United Nations Environment Programme, GEO Year Book 2006: An Overview of Our Changing Environment, Nairobi: UNEP Division of Early Warning and Assessment, 2006. 

Vinukollu, R. K., Wood, E. F., Ferguson, C. R., and Fisher, J. B.: Global estimates of evapotranspiration for climate studies using multi-sensor remote sensing data: Evaluation of three process-based approaches, Remote Sens. Environ., 115, 801–823,, 2011. 

Wardlow, B. D., Anderson, M. C., and Verdin, J. P.: Remote Sensing of Drought: Innovative Monitoring Approaches, CRC Press, Taylor and Francis Group, Boca Raton, FL, 484 pp., 2012. 

Short summary
Drought is a devastating natural hazard and difficult to define, detect and quantify. Global meteorological data and remote-sensing products present new opportunities to characterize drought in an objective way. In this paper, we applied the surface energy balance model SEBS to estimate monthly evapotranspiration (ET) from 2001 to 2018 over the dehesa area of the Iberian Peninsula. ET anomalies were used to identify the main drought events and analyze their impacts on dehesa vegetation.