Journal topic
Hydrol. Earth Syst. Sci., 23, 5111–5131, 2019
https://doi.org/10.5194/hess-23-5111-2019
Hydrol. Earth Syst. Sci., 23, 5111–5131, 2019
https://doi.org/10.5194/hess-23-5111-2019

Research article 17 Dec 2019

Research article | 17 Dec 2019

# Evaluation of drought representation and propagation in regional climate model simulations across Spain

Evaluation of drought representation and propagation in regional climate model simulations across Spain
Anaïs Barella-Ortiz and Pere Quintana-Seguí Anaïs Barella-Ortiz and Pere Quintana-Seguí
• Observatori de l'Ebre, Universitat Ramon Llull – CSIC, Roquetes, Spain

Correspondence: Anaïs Barella-Ortiz (anais.barella@gmail.com)

Abstract

Drought is an important climatic risk that is expected to increase in frequency, duration, and severity as a result of a warmer climate. It is complex to model due to the interactions between atmospheric and continental processes. A better understanding of these processes and how the current modelling tools represent them and characterize drought is vital.

The aim of this study is to analyse how regional climate models (RCMs) represent meteorological, soil moisture, and hydrological drought as well as propagation from precipitation anomalies to soil moisture and streamflow anomalies. The analysis was carried out by means of standardized indices calculated using variables directly related to each type of drought: precipitation (SPI), soil moisture (SSMI), runoff (SRI), and streamflow (SSI).

The RCMs evaluated are the CNRM-RCSM4, COSMO-CLM, and PROMES. All of the simulations were obtained from the Med-CORDEX database and were forced with ERA-Interim. The following datasets were used as references: SAFRAN (meteorological drought), offline land surface model simulations from ISBA-3L and ORCHIDEE (soil moisture drought), a SIMPA hydrological model simulation, and observed streamflow (hydrological drought).

The results show that RCMs improve meteorological drought representation. However, uncertainties are identified in their characterization of soil moisture and hydrological drought, as well as in drought propagation. These are mainly explained by the model structure. For instance, model structure affects the temporal scale at which precipitation variability propagates to soil moisture and streamflow.

1 Introduction

Impacts from recent climate-related extremes show that ecosystems and human systems are significantly vulnerable and exposed to current climate variability (IPCC2014). In addition, most of these impacts have been observed to increase in climate change scenarios. For instance, drought frequency has increased in Mediterranean regions during recent decades .

Drought is a complex phenomenon with important impacts on the environment as well as on society and the economy . The latter case, corresponds to socioeconomic drought and is outside the scope of this paper. There are several types of drought depending on which part of the system suffers from water deficit:

• meteorological drought, which corresponds to a period during which precipitation is considerably lower than its average level;

• soil moisture drought (also known as agricultural drought), which is due to a deficit of water in the soil's unsaturated zone, making it impossible to cover the crop's needs;

• and hydrological drought, which occurs when runoff, streamflow, and water levels in rivers, lakes, and groundwater are low. A more precise definition refers to superficial and groundwater water availability decreases in managed systems, affecting water demand.

As precipitation deficits propagate through the hydrological system , drought types are interrelated. Soil moisture drought is slightly out of phase with meteorological drought in non-irrigated areas and depends on several factors, such as the soil type and capacity to retain water as well as actual evapotranspiration. This is also applicable to regions where irrigation is carried out under unusual circumstances, for example, to avoid crop loss. In irrigated areas, soil moisture drought is directly related to the availability of irrigation water and therefore depends on hydrological drought. Hydrological drought is also affected by meteorological drought but generally at a larger temporal scale than soil moisture drought. Therefore, each system component is characterized by its own propagation dynamics and memory.

A better understanding of the different types of drought and their propagation processes is key to improving current and future drought representation and, thus, drought prediction and management tools. Regional climate models (RCMs) can help in this task. One of their greatest advantages is their adaptation to the regional scale . This allows us to use regional observations and to increase the model's physical parametrization complexity. Such models are suitable for drought analyses performed at the scale of large river basins. Another important aspect of RCMs is that these models are used to develop regional climate change scenarios via downscaling processes . In addition, using high resolution permits RCMs to be used to perform studies of atmospheric phenomena at a small scale. However, drought is difficult to model, as complex interactions between atmospheric and continental surface processes must be combined with human action . Furthermore, the relationships among drought types add complexity to such modelling.

Drought studies have gained interest in recent years due to increasing concern regarding (i) complexity in hydrological resource management, (ii) the escalation of the frequency, severity, and duration of extreme events, (iii) the increase in both the vulnerability to such events and the probability of being affected by them, and (iv) climate change . Drought studies are normally performed using drought indices, among which the Palmer drought severity index (PDSI) (Palmer1965) and the standardized precipitation index (SPI) are two of the most widely used. According to , standardized indices can be computed using other variables besides precipitation. For instance, used data from the Soil Moisture Active Passive (SMAP) satellite and the North American Land Data Assimilation System (NLDAS) to compute a standardized soil moisture index for drought warning. calculated a standardized runoff index using a simulation from the Variable Infiltration Capacity (VIC) model. They compared its behaviour to that of the SPI and concluded that a runoff-based index has the potential to complement other climate indices. analysed a method to compute a standardized streamflow index to compare streamflow hydrological conditions both spatially and temporally in a precise way.

Several studies analyse meteorological drought by means of RCMs. For instance, used the SPI and a version of the PDSI to analyse drought representation by 14 RCMs from the ENSEMBLES project at a European scale. analysed drought characteristics using RCMs from the PRUDENCE project across Britain and at a European scale by means of the drought severity index (DSI). and analysed projected changes to drought aspects in Canadian prairies by means of RCMs. The first study used three severity drought indices, whereas the second one used the SPI and PDSI. Other studies employed RCMs to analyse drought trends or to study the capacity of downscaled data to represent drought spells . RCMs are also used in soil moisture and hydrological drought analyses. However, there have been fewer studies about these types of drought. and studied hydrometeorological drought by computing precipitation and runoff standardized indices using simulations of the Weather Research and Forecasting (WRF) RCM driven by general climate models (GCMs). studied the impact of climate change on drought. They analysed drought characteristics and propagation impacts on meteorological, soil moisture, and hydrological droughts in Illinois using GCM-RCM nested simulations. Soil moisture and hydrological drought are also analysed by setting a given accumulated period to compute the drought index. studied how the WRF model improved dry and wet period detection using the SPI and the standardized precipitation–evapotranspiration index (SPEI) computed at 3- and 12-month timescales to study episodes related to soil moisture and hydrological droughts. used the SPI to analyse drought spatial patterns as a function of the timescales used in the index computation. For this, the SPI was computed at several timescales to detect soil moisture and hydrological droughts. Other studies that analysed the effect of the accumulated period are and .

Taking into account the different types of drought analyses performed with RCMs and the need to better understand droughts at the current time, but also in the future, the regional modelling tools must be evaluated. The work presented here has two objectives:

• to analyse how RCMs characterize meteorological, soil moisture, and hydrological drought;

• and to analyse how RCMs represent the propagation from a precipitation anomaly to a soil moisture and streamflow anomaly.

While the first objective aims towards contributing to drought representation analysis by RCMs, the second one assesses how these models simulate drought propagation. The latter objective is a key issue for the understanding of drought that has not been fully analysed from a RCM perspective. In this study, each type of drought is characterized using standardized indices that define drought according to the variability in a given variable: precipitation (meteorological drought), soil moisture (soil moisture drought), and streamflow and total runoff (hydrological drought). The work carried out analyses three RCMs: (i) the “Centre National de Recherche Météorologique” regional climate system model (CNRM-RCSM4, ), (ii) the COnsortium for Small-scale MOdelling (COSMO) model in CLimate Mode (COSMO-CLM, ), and (iii) PROMES . The reference data used are (i) the “Système d'Analyse Fournissant des Renseignements Atmosphériques à la Neige” (SAFRAN) atmospheric analysis , (ii) LSM offline simulations from the three-layer version of the “Interaction Sol-Biosphère-Atmosphère” (ISBA-3L) and the Organising Carbon and Hydrology In Dynamic EcosystEms (ORCHIDEE) models, (iii) a simulation from the “Sistema Integrado de Modelización Precipitación-Aportación” (SIMPA) hydrological model , and (iv) observed streamflow.

2 Area and study period

This study was carried out across mainland Spain, which is known to experience frequent droughts (Olcina2001). The time period was limited to the availability of RCM simulations and observations: from 1989 to 2008 for the meteorological and soil moisture drought and from 1989 to 2005 for the hydrological drought.

Precipitation across mainland Spain is complex and highly influenced by relief. The main climatic regimes are oceanic and Mediterranean. However, a semiarid climate or even a desert-like climate can be identified in the southeast (Aemet2011). In general, Spain is a semiarid region and is not densely vegetated. As a result, soil moisture displays a large annual cycle. From a hydrological point of view, there is a strong dependence of the main rivers on the precipitation generated in the nearby relief and the resulting runoff. An example is in the Ebro Basin and the Pyrenees. In addition, the impact of the anthropic effect must be taken into account because there is a wide network of dams and river canals that, in some cases, operate between basins.

According to , drought in Spain has increased in severity and frequency. Although precipitation does not show significant annual trends, observations show a reduction in spring and summer as well as an increase in the number of consecutive dry days . Both aspects have an impact on soil moisture drought, which is also affected by an increase in annual and seasonal temperatures . This rise in temperature increases the atmospheric demand and thus evapotranspiration, reducing the soil's water content. For hydrological drought, we must also consider the advance in the thaw date and a thinning of the blanket of snow in mountainous areas, such as in the Pyrenees , which affect streamflow and increase this type of drought. However, it should be noted that snow melt can affect streamflow and thus hydrological drought in different ways depending on its timing .

3 Datasets

This section describes the RCMs (Sect. 3.1) and the reference products (Sect. 3.2) used in the study.

## 3.1 Regional climate models

In this study, drought representation and propagation in three RCMs are analysed. RCM simulations were downloaded from the Med-CORDEX database, which is a contribution to the Coordinated Regional Climate Downscaling Experiment (CORDEX, ) focusing on the Mediterranean region . The criterion used to select the models was that each one used a different surface scheme and therefore represented physical processes related to precipitation, soil moisture, and surface and sub-surface runoff in different ways. The three RCMs selected are listed in Table 1 and described below:

• The CNRM-RCSM4 is a RCM developed by the CNRM. It includes the regional climatic atmospheric model “Aire Limitée Adaptation dynamique Développement InterNational” (ALADIN-Climate, ), the three-layer version of the ISBA LSM , the Total Runoff Integrating Pathways (TRIP) routing scheme , and the regional ocean model NEMOMED8 . Hereafter, it will be referred to as RS4.

• The COSMO-CLM (CCLM) model is the climate version of the COSMO model developed by the Goethe Universität Frankfurt (GUF). The surface scheme is a multilayer version of the two-layer model. Hereafter, it will be referred to as CL4.

• The PROMES model was developed by the Universidad Complutense de Madrid (UCM) and the Universidad de Castilla-La Mancha (UCLM). It is coupled to the ORCHIDEE LSM . Hereafter, it will be referred to as PMS.

All of the RCM simulations are driven by the ECMWF Interim reanalysis (ERA-Interim) , which provides a global atmospheric reanalysis that starts in 1979 and is continuously updated in real time. In addition, it improves some important issues pertaining to ERA-40, such as the representation of the hydrological cycle. This reanalysis is performed by means of a data assimilation system based on a 2006 release of the ECMWF's Integrated Forecast System, IFS (Cy31r2), and uses a four-dimensional variational analysis (4D-Var) with a 12 h analysis window. The database has atmospheric and surface parameters with a temporal scale of 6 and 3 h, respectively. The spatial resolution is 80 km, with 60 vertical levels from the surface to 0.1 hPa.

ERA-Interim is a well-known atmospheric forcing used in a large number of studies. For instance, and have validated it across the Iberian Peninsula. However, biases in this type of forcing have a negative effect on LSM simulations, which can be corrected .

In this study, ERA-Interim is the driving data of the three RCMs analysed. In addition, it is also used to force LSM simulations used as a reference in the soil moisture drought analysis. Hereafter, it will be referred to as ERA.

Table 1RCMs analysed in this study.

## 3.2 Reference products

This section is divided into three subsections to describe the products used as reference data in the meteorological (Sect. 3.2.1), soil moisture (Sect. 3.2.2), and hydrological (Sect. 3.2.3) drought analyses. The products are listed in Table 2.

Table 2Reference products used in the meteorological, soil moisture, and hydrological drought analyses.

### 3.2.1 Meteorological drought reference

SAFRAN is a meteorological analysis system developed by Météo-France. It provides estimates of the following variables: precipitation, 2 m temperature, 10 m wind speed, 2 m relative humidity, and cloudiness, as well as modelled downward visible and infrared radiation following the radiation scheme of . For this, an optimal interpolation algorithm (Gandin1966) that combines observations and a first guess is used. The first guess employed is ERA for all variables, except for precipitation, which is obtained from observations, and this is the reason why it has been selected as the reference dataset for the meteorological drought analysis. The meteorological station data belong to the Spanish Meteorological State Agency (Agencia Estatal de METeorología, AEMET) network. Precipitation is analysed by means of daily observations, whereas the remaining variables are analysed every 6 h. The data are then interpolated hourly using different methods that depend on the variable. In Spain, SAFRAN was extended for a 35-year period (1979/1980–2013/2014) and was implemented and validated over the Ebro Basin (northeastern Spain) .

In this study, the standard version of SAFRAN, at a 5 km resolution, has been regridded to a 30 km resolution to compare it with the RCM simulations as the reference in meteorological drought analysis and to force LSM simulations which are the reference in the soil moisture drought analysis. Hereafter, it will be referred to as SLR (SAFRAN low resolution).

### 3.2.2 Soil moisture drought reference

Offline LSM simulations are used in this study as a reference to analyse soil moisture drought because there is a lack of soil moisture data across mainland Spain that are suitable for studies that require a large spatial coverage. Internally, in RCMs, LSMs are bidirectionally coupled to the atmospheric model to simulate surface processes. However, when LSMs are run offline (forced by gridded databases), biases due to the atmospheric model and the coupling between the atmospheric model and the LSM are avoided. This makes offline LSM simulations good reference datasets to study drought. We selected three offline LSM simulations that use the same surface schemes employed by two of the RCMs analysed in this study.

The ISBA LSM , developed by the CNRM, is composed of various modules that simulate heat and water transfer in the soil, vegetation, snow, and surface hydrology. This scheme has evolved using different approaches to model the soil. As a result, there are several versions that can be used. For this study, we selected the ISBA-3L , which considers a three-layered description of the soil. It should be stressed that the ISBA-3L is limited by certain aspects. For example, underground water is not represented, and there is no horizontal water transfer. Most of the ISBA's soil and vegetation parameters are derived from the ECOCLIMAP land cover database at a 1 km resolution , such as the cover types and soil texture. This database considers more than 550 land cover types from all around the world. Its vegetation variability depends on the location, climate, and phenology.

The ORCHIDEE LSM was developed by the Institut Pierre-Simon Laplace (IPSL). It can be run in a stand-alone mode or coupled to the Laboratoire de Météorologie Dynamique (LMD-Z) general circulation model (Li1999), which was developed by the LMD in Paris. Hydrology is approached by means of a diffusive equation with a multilayer scheme. For this, the Fokker–Planck equation is solved considering a soil depth of 2 m distributed across 11 layers. The fine resolution is key to better model the interaction between the root profile and the soil moisture distribution at different depths as well as infiltration processes. In addition, ORCHIDEE includes sub-grid variability in soil moisture. Each grid box is divided into three soil moisture profiles with different vegetation distributions, but the same soil texture and structure that are obtained from the Zobler map .

In this study, three offline LSM simulations are used in the soil moisture drought analysis. Two of them use the three-layer version of the ISBA model (one simulation was forced with SLR and another one with ERA), and the other simulation uses the ORCHIDEE LSM forced with ERA. Hereafter, the ISBA and ORCHIDEE LSMs will be referred to as the ISB and ORC, respectively.

### 3.2.3 Hydrological drought reference

For the hydrological drought analysis, modelled and observed streamflow data were used as references. Two issues should be stressed before explaining these datasets. The first is that the Med-CORDEX database does not provide simulated streamflow for any of the three RCMs, which would be the variable ideally suited for this study. In the absence of such data, it was decided to use modelled total runoff (hereafter referred to as runoff) corresponding to the subbasins defined by a selection of gauging stations. We believe that this approximation is valid because we use a coarse time step, with a larger time propagation than the flow propagation. In fact, other studies use this variable to analyse hydrological drought . The second issue is that Spanish basins are highly influenced by human management. However, RCMs do not simulate water management procedures but do simulate natural regime behaviour. Bearing these considerations in mind, only gauging stations satisfying the following criteria were considered:

1. A 95 % data completeness during the study period; this assures that the monthly series of observations has few gaps.

2. An area greater than 10 000 km2; the analysis was limited to large areas because streamflow is approximated by runoff, which is likely to perform poorly in small basins considering the coarse resolution of RCM simulations.

3. A KGE between SIMPA and the observations greater than 0.5 for the consideration of a near-natural regime. Stations had to be as natural as possible and their corresponding basins large enough to be compared to a low-resolution RCM. This is difficult in Spain due to the high degree of human influence. A high KGE value between SIMPA (naturalized flow) and the observations would indicate near-natural regime behaviour. Therefore, we performed a sensitivity analysis testing different values, and 0.5 was a reasonable compromise between the near-natural regime and a sufficient number of stations for the analysis. It is true that this threshold does not represent a completely natural regime, and thus some human influence can be present. However, it allows a fair comparison of the RCMs and observations.

The first criterion was satisfied by 87 stations, the second by 13 of the 87 stations, and the third by 8 of the 13 stations, which is the final selection. Figure 1 shows the locations of these stations as well as the relief and the river network. Tables 45, and 6 contain their code, their area, and the basin to which they belong.

SIMPA is a Spanish acronym meaning “Integrated System for Rainfall–Runoff Modelling” (“Sistema Integrado de Modelización Precipitación-Aportación”) . It is a conceptually distributed hydrological model for water management developed by the Spanish Centre for Studies and Experimentation on Public Works (“Centro de Estudios y Experimentación de Obras Públicas”, CEDEX). SIMPA provides estimates of the water cycle's main components, such as precipitation, evapotranspiration, and river discharge, at a monthly scale in a natural regime. This regime is characterized by the free flowing of water, with no aspect unrelated to the environment (such as dams for water resource management) affecting it. Data are provided on a 1 km2 grid. The model employs its own forcing, which uses an observational dataset similar to that of SAFRAN.

In this study SIMPA is used as a reference product for streamflow to analyse hydrological drought. Hereafter, it will be referred to as SMP.

Daily streamflow observations were also used as a reference in the hydrological drought analysis. These belong to the Spanish Ministry for the Ecological Transition, which provides data for basins comprising more than one region. Hereafter, they will be referred to as OBS.

Figure 1Map showing the relief, river network, and the locations of the gauging stations.

4 Methodology

The methodology employed in this study is the same as that used in . Therefore, the two studies complement each other, providing a wider analysis of drought representation by both LSMs and RCMs.

Drought can be characterized in several ways. Here, we used standardized indices that define drought according to the variability in a given variable. It should be noted that variability is a key aspect of drought analysis. For the meteorological and soil moisture drought analyses, we computed the SPI and a standardized soil moisture index (SSMI) using precipitation and soil moisture data, respectively. For the hydrological drought analysis we computed two indices: the standardized streamflow index (SSI) using streamflow and the standardized runoff index (SRI) using total runoff.

## 4.1 Drought index calculation

We followed the spirit of the SPI to compute the SSMI, SSI, and SRI. In these indices, the variable's time series is transformed from its original distribution to a normal distribution. The resulting values correspond to the number of standard deviations by which the anomaly deviates from the mean. Biases were not computed, as they are zero by construction. The computation was carried out using monthly data for all indices. In the case of the SPI, a time series of the accumulated precipitation from the previous n months (with n being the index scale) was also calculated to perform the drought propagation analyses. It should be noted that in the SPI methodology, data are fitted to the corresponding parametric distribution, which can be an issue in studies that standardize several variables, as in our case. To solve this problem, we used the nonparametric methodology described by .

On the one hand, meteorological drought representation by RCMs was addressed by comparing their SPI-12 time series with those of ERA and SLR. Using results from an accumulation period of 12 months is particularly robust for seasonal reasons. In addition, this scale is suitable for the study area defined from a hydrological resource management perspective. Special attention was paid to differences in duration, severity, and area. In addition, the temporal correlation of the RCMs' SPI-12 with that of ERA and SLR was computed to identify regions where drought representation improved, if there were any. On the other hand, the analysis of the RCMs' soil moisture and hydrological drought representation was performed via the calculation of root-mean-square difference (RMSD) and the Pearson correlation coefficient (r) at a 95 % significance level of the standardized indices computed using an accumulation period of 1 month. We would like to stress that an RMSD equal to or greater than 0.5 often implies a change in drought category (i.e. from moderate to severe, for example) according to the SPI drought classification scale . To facilitate the Pearson correlation analysis interpretation, the guideline proposed by was followed.

## 4.2 Meteorological drought propagation

The analysis of how meteorological drought propagates to soil moisture and hydrological drought is useful in hydrological resource management and allows for the detection of similarities and differences in the way models address the physical processes that drive this propagation. The methodology employed to analyse meteorological drought propagation is based on . The soil moisture memory of precedent precipitation will vary from one point to another depending on the location because soil moisture is controlled by different aspects regarding climate and soil as well as vegetation properties such as precipitation, actual evapotranspiration, soil texture, stomatal resistance, and root depth, among others. To determine the month scale at which the precipitation deficit propagates to a soil moisture deficit, we carried out the following analysis:

• computed the standardized soil moisture index with a time accumulation of 1 month (SSMI-1);

• computed the SPI with a time accumulation of 1 to 28 months (SPI-nx);

• identified the nx scale that maximizes the correlation between SPI-n and SSMI-1.

The same methodology was employed to analyse how meteorological drought propagates to hydrological drought. For this, the SSI-1 and the SRI-1 were computed instead of the SSMI-1, and the SPI computation was performed using the areal mean of the basin precipitation defined by each studied gauging station.

## 4.3 Streamflow validation

The streamflow validation was carried out using the Kling–Gupta efficiency (KGE) . The optimal KGE value is 1, whereas negative values are sign of a model's bad performance.

5 Results

First, meteorological drought representation by RCMs is evaluated (Sect. 5.1). Next, soil moisture drought representation by RCMs is analysed as well as how the models address the transition from meteorological to soil moisture drought (Sect. 5.2). Finally, hydrological drought representation by RCMs and the propagation of meteorological to hydrological drought are analysed (Sect. 5.3).

## 5.1 Meteorological drought

In this section, we will focus on precipitation because its decrease or absence are the main cause of meteorological drought. The SPI is computed for the RCMs, ERA (driving data), and SLR (reference dataset). On the one hand, the comparison of the RCMs' SPI with that of ERA allows us to identify and determine the extent to which these models reproduce the structures of their driving data. On the other hand, the comparison of the RCMs' and SLR's SPI shows the extent to which the RCMs improve these structures, as the SLR dataset is based on observations.

Figure 2Mean annual precipitation across the study area from 1989 to 2008: ERA (a), SLR (b), RS4 (d), CL4 (g), and PMS (j). Panel (c) shows the difference between ERA and SLR mean annual precipitation; panels (e)(h), and (k) show the difference between the RCMs and ERA mean annual precipitation; and panels (f)(i), and (l) show the difference between the RCMs and SLR mean annual precipitation.

### 5.1.1 Mean annual precipitation comparison

Figure 2a–c shows the mean annual precipitation of ERA and SLR as well as the difference in mean annual precipitation between them. Figure 2d–l shows the RCMs' mean annual precipitation and their difference with respect to ERA (Fig. 2e, h, k) and SLR (Fig. 2f, i, l). All products show greater precipitation in the northwestern and northern regions of the peninsula (exceeding 2000 mm yr−1) as well as over mountainous chains. The products also show that precipitation is lower along the main basin valleys (due to the orographic shadow effect) and minimal over the southeast, which is the driest region of the peninsula. For instance, precipitation across some areas of this region does not exceed 100 mm yr−1. The RCMs' mean precipitation spatial structures show similar behaviour to those from ERA and SLR. The fact that precipitation is high over mountainous chains indicates the strong influence of relief, which is key in the way water from precipitation is distributed. In fact, we would like to stress SLR's significant contrast in relief due to its use of data from AEMET's dense pluviometric network. This evidence indicates the complex spatial structure of precipitation in Spain.

Regarding the RCMs, RS4 and PMS show the greatest similarity and the highest contrast with CL4. When compared to ERA, both RCMs have higher precipitation, especially in mountainous areas. Modelled precipitation tends to overestimate precipitation compared with observations and ERA . This reflects the addition of water in the form of precipitation, improving the RCMs' spatial distribution of precipitation with respect to the driving data. However, when RS4 and PMS are compared to SLR, precipitation is underestimated over some areas (valleys and the coastline). It must be noted that SAFRAN is mainly based on rain gauge information. The CL4 model is a different matter, as it underestimates precipitation for almost all of Spain when compared with ERA and SLR.

Figure 3SPI-12 time series calculated with the spatially averaged time series of mainland Spain precipitation, as reproduced by ERA (a), SLR (b), RS4 (d), CL4 (g), and PMS (j). Panel (c) shows the difference between ERA and SLR SPI-12; panels (e)(h), and (k) show the difference between the RCMs and ERA SPI-12; and panels (f)(i), and (l) show the difference between the RCMs and SLR SPI-12.

### 5.1.2 SPI comparison

Once the spatial distribution of precipitation in the peninsula is described, we can study the variability in the different products and thus their capacity to reproduce drought spells.

Figure 4Differences in the proportion of area under drought (SPI-12 $<-\mathrm{1}$) estimated by (i) the RCMs and SLR with respect to ERA (a) and (ii) the RCMs and ERA with respect to SLR (b).

Figure 3 shows the time series of the SPI-12 calculated using mainland Spain's average precipitation as reproduced by ERA (panel a), SLR (panel b), and the RCMs (panels d, g, and j). The computation is performed for a time accumulation of 12 months. ERA and SLR show several drought spells which occurred during the 20 years that comprise the study period (the most severe occurred in 2005–2006). These spells coincide with those found by and also appear in the RCMs' SPI time series plots. Therefore, RCMs are capable of reproducing these spells. However, differences in duration and severity can be observed. For instance, the duration of the spell that occurred in 1992 and 1993 was 21 and 22 months according to CL4 and PMS, respectively. However, it lasted 19 and 17 months according to SLR and ERA, respectively. Another interesting result is found in the spell that took place in 2002. CL4's mean severity is similar to that of ERA. However, the mean severity of both RS4 and PMS is in better agreement with that of SLR than with their driving data. Other information that can be extracted from these results is the timing of drought between the RCMs and their driving data and reference data. For example, the RCMs, ERA, and SLR agree that the spell that occurred between 1994 and 1995 started in May 1994. RS4, CL4, ERA, and SLR agree that it ended in December 1995, but the PMS model indicates the end of the spell occurring a month before the rest (November 1995). The spell between 2004 and 2006 started in November 2004 according to SLR, ERA, and PMS, which show similar durations (ERA shows 23 months, and SLR and CL4 show 24 months). However, it started in January 2005 according to RS4 and CL4 and had a longer duration: 27 (RS4) and 25 (CL4) months. In a previous study, a spurious trend in ERA was detected . This can be observed in Fig. 3c, where the difference between the SPI-12 time series of ERA and SLR is represented. However, the differences between the RCMs and SLR SPI-12 time series show that these models do not drag this trend.

Figure 5Map of the temporal correlation (Pearson) between the SPI-12 time series produced by (i) the RCMs and SLR with respect to ERA (a) and (ii) the RCMs and ERA with respect to SLR (b). Portugal and regions from the study area (mainland Spain) whose values are not within the colour scale are represented in white.

To deepen the analysis of the spatial structure of meteorological drought, monthly SPI-12 maps were computed. The comparison of these structures, as shown by the RCMs, with those of their driving data and the reference dataset, provides more information about drought representation by RCMs. For instance, the temporal evolution of the spatial correlation of the SPI-12 maps of the RCMs with ERA and SLR (not shown) indicates the similarity between the RCMs and their driving data and reality (as approximated by SLR), respectively. In the first case, RS4 most resembles ERA. Despite not showing the highest correlations, it has less variability and is therefore more robust. The correlations of CL4 and PMS with ERA are more variable, reaching values close to zero in some months, in which the spatial structure of drought is not captured. In the second case, RCMs show worse correlation with SLR than with ERA, as expected. It should be stressed that RS4 also displays better correlations than the other RCMs when these are compared with SLR. Thus, out of the three RCMs studied, RS4 deviates less from the driving data and most resembles the reference. The correlation between ERA and SLR shows some variability in its temporal evolution, especially since 2000. This is likely due to the effect of the spurious trend identified in Fig. 3. To complement the spatial structure analysis, Fig. 4 provides information about the difference in the percentage of the area affected by drought (SPI-12 $<-\mathrm{1}$) between the RCMs and (i) ERA (panel a) and (ii) SLR (panel b). The differences are generally under 25 %. In general, the RCMs show similar behaviour during a drought spell, meaning that they all overestimate or underestimate the affected area. For example, all three underestimate the area in 1994, overestimate it in 1995, and underestimate it in 1996. The difference relates to the degree to which they deviate. On the one hand, in 1995, RS4 overestimates the percentage of area affected by drought by 20 %, CL4 by approximately 15 %, and PMS by less than 10 %. On the other hand, in 1996, PMS underestimates this percentage by more than 30 %, while RS4 and CL4 underestimate it by 10 %. Another example is how the RCMs overestimate the area in 2000 and how they underestimate it between 2001 and 2003. RS4 is the RCM that shows the lowest percentage difference, which is consistent with our previous results.

Finally, Fig. 5 provides an idea of the spatial structures of drought representation. For this, the temporal correlation (Pearson) of the SPI-12 between the RCMs and (i) ERA (maps from Fig. 5a) and (ii) SLR (maps from Fig. 5b) is computed for each grid point. The correlations of the RCMs with ERA are above 0.8 in the south and worsen towards the north and the eastern coast, reaching 0.2. The Ebro Basin is the region showing the poorest values and thus the largest differences in drought representation between the RCMs and ERA. Differences in this area are also identified in the map showing the correlation between ERA and SLR, where values are even negative. However, when the RCMs are compared with SLR (Fig. 5b), the correlations are higher across this basin. Consequently, RCMs improve drought representation across this region.

## 5.2 Soil moisture drought

This section analyses how RCMs reproduce soil moisture drought (Sect. 5.2.1) and its propagation from precipitation to soil moisture (Sect. 5.2.2). Offline LSM simulations from the ISBA-3L and ORCHIDEE are used as references. We should bear in mind that in these simulations, soil processes do not impact the atmosphere, whereas the RCM simulations are performed in coupled mode and thus interact with the atmosphere.

### 5.2.1 SSMI comparison

The RMSD and the Pearson correlation coefficient (r) are calculated by comparing the SSMI from the RCM simulations and the LSM simulations. All mesh points and time steps are included in the comparison. Therefore, the results provide information regarding spatial and temporal drought structures. It should be noted that biases are not calculated because they are zero by construction (the mean of the SSMI is zero). The results are shown in Table 3, the upper block of the table corresponds to the RMSD and the lower block corresponds to the r calculations. A colour scale consisting of a blue (largest similarity between models) to red (lowest similarity between models) gradient via white has been included to facilitate reading.

To put these results into context, we will consider the drought classification according to the SPI, which is divided into eight categories from “extremely wet” (SPI =2) to “extreme drought” (SPI $=-\mathrm{2}$). A RMSD equal to 1 is a standard deviation of the index studied (in this case, the SSMI). In the framework of drought analysis, a value higher than 0.5 would imply a change in category (for example, from “slightly wet” to “moderately dry”). Therefore, the upper block of Table 3 shows that there is a change in drought category when comparing the RCMs with the reference offline LSM simulations, as the RMSD is above 0.5 (fourth to sixth columns). In addition, the three RCMs compared among them also represent soil moisture droughts of different categories (second and third columns), as expected.

Going into detail, we can observe some similarity between the RS4 simulation, which uses the ISBA surface scheme, and the ISB simulations. Compared with ISB, RS4 reproduces drought better than the other RCMs used in this study. However, this does not occur when the PMS and ORC simulations are compared, despite using the same surface scheme. In fact, RS4 and CL4 are, in general, more similar to ORC than PMS, according to the statistics. The SSMI comparison of the RCMs with both ISB simulations (fifth column vs. sixth column) shows very similar statistics. Discrepancies could be explained by the RCMs' land–atmosphere coupling and forcing effects.

Table 3Comparison of the SSMI data from the RCM and LSM simulations. The upper block shows the RMSD, and the lower block shows the Pearson correlation (r). The colour scale is a gradient from blue (largest similarity between models) to red (lowest similarity between models) via white.

### 5.2.2 Propagation to soil moisture drought

nx maps from Fig. 6 indicate the scale in months at which the correlation between the SSMI-1 and SPI-n is maximal and thus the temporal scale at which meteorological drought propagates to soil moisture drought. Figure 6a–c shows RCM maps, whereas Fig. 6d–f shows ISB and ORC LSM (reference data) maps. The scale ranges from 0 to 28 months, with the dynamics of the model in regions with a yellowish tone being slower than those in regions with a bluish tone.

Figure 6a–c shows that the RCMs provide different results even though they use the same driving data, which indicates the predominance of model structure with respect to the driving data. This becomes more evident when a RCM is compared to the LSM that has the same surface scheme. For example, the RS4 and ISB (Fig. 6a, d, and e) maps show similar spatial patterns. These are very homogeneous, with scales that range from 1 to 4 months, implying that ISB reacts very quickly to precipitation. Another example is the comparison between PMS and ORC (Fig. 6c and f). Both models show greater heterogeneity than ISB, with scales from 1 to 20 (PMS) and 24 (ORC) months, highlighting the role of the continental surface. Finally, the CL4 (Fig. 6b) behaviour is quite homogeneous. The peninsula is divided into two areas, one over the northwest, where the nx scale ranges between 6 and 12 months, and a larger one with a fixed value of 20 months. It is interesting to note the similar spatial structures from the ERA-CL4 and ERA-ISB maps (Fig. 6b and e), indicating that soil moisture drought propagation by CL4 drags the same spatial structures as its driving data.

Figure 6The nx timescale maximizing the correlation between the SPI-nx and SSMI-1 for the RCM and LSM simulations. Portugal and regions from the study area (mainland Spain) whose values are not within the colour scale are represented in white.

## 5.3 Hydrological drought

In Sect. 5.3.1, streamflow and aggregated runoff will be compared. Next, the RCMs' capacity to simulate hydrological drought will be analysed by comparing the SSI and SRI (Sect. 5.3.2). Finally, we will analyse the way in which RCMs reproduce meteorological drought propagation to hydrological drought (Sect. 5.3.3).

Table 4KGE of the aggregated runoff simulated by the RCMs compared with the reference streamflow of the SMP model and OBS. Note that observations are affected by anthropic effects. Negative values (model's poor performance) are shown in red. The colour scale only displays positive KGE values, and consists of a blue (best performance between the RCMs and SMP or OBS) to red (worst performance between the RCMs and SMP or OBS) gradient via white. The bottom row shows the average KGE for each RCM.

### 5.3.1 Streamflow and aggregated runoff comparison

Table 4 shows the KGE comparing RCM aggregated runoff with SMP and OBS as well as the comparison between SMP and OBS. All of these comparisons are performed at a monthly scale. Negative values, indicating poor performance, are marked in red. For the positive KGE values, a colour scale consisting of a blue (best performance between the RCMs and SMP or OBS) to red (worst performance between the RCMs and SMP or OBS) gradient via white has been included to facilitate reading.

The RCMs show KGE values that are negative or close to zero at many stations. Focusing on the positive KGE values, CL4 has the worst performance, as the values are below 0.5. Although RS4 is the RCM showing the highest KGE value (0.7 at station number 9002 for the comparison between RS4 and OBS), PMS performs better across both basins according to the average value (bottom row of Table 4).

Table 5Comparison of the SRI from the RCMs with the SSI from SMP. The table shows the RMSD and the Pearson correlation (r). The colour scale is a gradient from blue (largest similarity between models) to red (lowest similarity between models) via white.

An analysis comparing the temporal series (not shown) of the RCMs' aggregated runoff with the reference streamflow (SMP and OBS) shows that CL4 sustains summer flows, which is a positive aspect. This is not the case for RS4 and PMS, which have steeper recession curves. In addition, RS4 is found to overestimate streamflow peaks at stations in the Duero Basin. Finally, we would like to note that RS4 and PMS show a different behaviour from that of the reference data at station number 9025 (which represents the Segre subbasin). Both models overestimate a large number of peaks. We believe that this may be explained by the fact that the Segre subbasin is nival.

### 5.3.2 SSI and SRI comparison

Table 5 shows the RMSD and the Pearson correlation coefficient, r, comparing the RCMs' SRI with SMP's SSI. A colour scale consisting of a blue (largest similarity between models) to red (lowest similarity between models) gradient via white has been included to facilitate reading.

The results are similar to those in Sect. 5.2.1, as the RMSD values are above 0.5 (which often indicates a change in drought category). Therefore, according to the SPI drought classification, the RCMs represent hydrological droughts of a different category from that of the reference dataset. Nevertheless, the RCMs show moderate positive correlations with SMP according to the guideline proposed by . We would like to note that the worst statistics are shown by station 9025, which is identified in the previous section as the station showing the worst performance. It is highly influenced by the effects of snow, which could explain the high RMSD and weak r values.

### 5.3.3 Propagation to hydrological drought

Table 6 shows the nx values that indicate the monthly scale at which the correlation between the SPI-nx and (i) SRI-1 (RCMs) and (ii) SSI-1 (SMP and OBS) is maximal. This can be interpreted as the temporal scale at which meteorological drought propagates to hydrological drought. To better understand these results and the extent to which they reflect the propagation of a precipitation anomaly to a streamflow anomaly, a colour scale is included to indicate the strength of the correlation. It follows the guideline given by : very strong (yellow), strong (blue), moderate (white), and weak (red).

In contrast to Sect. 5.2.2, RS4 and PMS show very similar scales: means of 3 months for stations in the Duero Basin and 2 months (RS4) and 1 month (PMS) for stations in the Ebro Basin. However, CL4 provides larger scales, from 9 to 13 months (Duero Basin) and from 1 to 8 months (Ebro Basin). The difference in scales shown by the RCMs is an indicator of the relevance of model structure in drought propagation. This result is also obtained in Sect. 5.2.2.

The RCMs and SMP show higher nx values and thus slower dynamics in the Duero Basin than in the Ebro Basin. This is in agreement with the nx values obtained using the SSI-1 computed with OBS, in which the mean nx values are 9 (Duero) and 4 (Ebro). Analysing these results, we can establish that the RS4 and PMS runoff responds quickly to precipitation anomalies. When compared to SMP and OBS, the monthly scales provided are in good agreement in the Ebro Basin but are too low in the Duero Basin. In contrast, CL4 shows larger scales and behaves inversely to RS4 and PMS, as it is in better agreement with SMP and OBS in the Duero Basin.

Focusing on the strength of the correlations, RS4 (coupled to ISB) shows strong positive correlations at six of the eight stations analysed, indicating that precipitation has a significant role in streamflow variability. However, the correlations of CL4 and PMS with the SPI-nx are moderate (even weak for station number 9025), which implies that there are also other factors driving this variability.

Table 6The nx timescale maximizing the correlation between the SPI-nx and (i) SRI-1 (RCMs) and (ii) SSI-1 (SMP and OBS). Scales longer than 12 months are marked in bold. A colour scale has been included to indicate the strength of the correlations between the SPI-nx and (i) SRI-1 and (ii) SSI-1 following the guide proposed by . The correlation ranges and the associated colours are as follows: (i) very strong, 0.80–1.0 (yellow); (ii) strong, 0.60–0.79 (blue); (iii) moderate, 0.40-0.59 (white); and (iv) weak, 0.20–0.39 (red).

6 Discussion

This study provides an assessment of how RCMs represent meteorological, soil moisture, and hydrological drought as well as the way in which precipitation anomalies propagate to soil moisture and streamflow anomalies. This assessment was performed using standardized indices, which we believe are a good option for performing drought analysis. The reason for this is that they describe drought based on the variability in a given variable; thus, each type of drought can be studied using the variable that best suits its characteristics. In this work, four indices were computed: SPI (precipitation), SSMI (soil moisture), SSI (streamflow), and SRI (total runoff).

RCMs provide a good representation of meteorological drought. The results show that they are capable of reproducing the same drought spells as those detected by their driving data and the reference data. However, they differ in terms of an event's duration, severity, and area, as expected. We have identified that RCMs improve drought representation with respect to the driving data in several aspects. For instance, the temporal evolution of the SPI-12 shows that the severity of some of the drought spells is closer to that of the reference data. In addition, they do not reproduce the spurious trend identified in the driving data, which could lead to a misrepresentation of the phenomenon. Finally, a temporal correlation analysis shows that drought representation is improved over the northeastern region of the Iberian Peninsula, which is a known limitation of global analysis across Spain. These results are consistent with previous studies showing that RCMs provide a suitable representation of drought using drought indices across Spain .

Unlike the previous analysis, the results regarding soil moisture and hydrological drought representation show differences when RCMs are compared among themselves and to the reference data. The analyses are carried out using the RMSD and Pearson correlations, and the observed uncertainty corresponds, in most cases, to a change in the drought category (according to the SPI drought classification) of RCMs with respect to the reference data. These differences are expected if we consider the following aspects:

• First, the reference data that are employed. LSM offline simulations are used as reference datasets for soil moisture drought analysis due to the lack of observations. In Spain, soil moisture data from the REMEDHUS network and the Valencia Anchor Station are available. However, these datasets are suitable for studies in which large spatial coverage is not an issue (such as model or satellite-derived data validation), which does not apply to our case. Remote sensing products could also be an option, but there are certain limitations that should be taken into account, for example, uncertainty sources, gaps in the data, and short time series . showed that different satellite products behave differently across a region representative of Mediterranean landscapes (Catalonia in the northeast Iberian Peninsula). Therefore, using these products would add more uncertainty to the study. In addition, the retrieved soil moisture data correspond to surface soil moisture, and, in this study, we considered root zone soil moisture. If the study was limited to the consideration of surface soil moisture, processes important to soil moisture drought would not be considered. Nevertheless, improvements in this discipline, as well as an increase in the length of time series, will convert these products into interesting alternatives for LSM simulations. Regarding the hydrological drought analysis, the observed streamflow in mainland Spain was used as a reference. However, the large number of dams and canals have a high anthropic impact on river systems, affecting the observed streamflow. As RCMs do not take these effects into account, a simulation from the SIMPA model that provides streamflow data considering a natural regime is necessary as a reference dataset. It would be interesting to include these effects in RCM modelling to perform drought analysis using the observed streamflow as reference data. This would provide an idea of the anthropic impact on hydrological drought.

• The second consideration is the use of LSMs. In this study offline LSM simulations are used as reference datasets to compare with the RCM simulated soil moisture, due to the lack of observed soil moisture data. This is natural as RCMs internally used a LSM coupled to the atmosphere. We used three simulations as references, instead of one, because LSMs vary from one to another in several important aspects, such as soil discretization, evapotranspiration formulation, runoff formulation, subgrid processes, and so on. All of these differences impact soil moisture and streamflow.

• Third, the use of simulated total runoff to compute the SRI is due to the lack of RCM modelled streamflow data to compute the SSI. Hydrological drought representation and propagation are likely to be affected by this approximation according to the KGE values shown in Sect. 5.3.1. When compared to SMP, the KGE using total runoff is under 0.25 in almost half of the cases, which indicates poor performance. However, we would like to stress that PMS is the RCM that best approximates streamflow according to the KGE values.

• Finally, the effects of RCM coupling with the atmosphere should be taken into account, as the soil moisture references are LSM simulations performed in an uncoupled mode and thus do not include atmospheric feedback.

We would also like to note the relevance of the timescale used to compute drought indices. According to and , the added value of the SPI is influenced by the accumulation period employed. This is interesting to consider because soil moisture and hydrological droughts have different scales of propagation, and differences are also found in the propagation analysis in Sect. 5.2.2 and 5.3.3.

The main objectives of this study are to evaluate drought properties and propagation in RCM simulations. At this point, we would like to mention uncertainty, which is an important issue in this kind of study. Although the analysis was performed using only three models, these models provide different results in terms of drought indices and, especially, drought propagation. Adding more models would show the spread with more detail, but we believe that the differences among the three models are large enough to show that RCM developers should look at these issues and that RCM users should take them into account.

A key result of the study is the relevance of the models' physics, which prevails over the driving data. This is shown in the soil moisture and hydrological drought representation evaluation as well as in the analyses of drought propagation. In the latter case, the model's structure influences the temporal scale at which the variability in precipitation affects that of soil moisture and streamflow. For instance, in the analysis in Sect. 5.2.2, the spatial patterns of RS4 (coupled to ISB) and ISB are very similar and show that ISB, and thus RS4, respond too quickly to precipitation. The spatial structures of PMS (coupled to ORC) and ORC differ to a greater extent. The temporal scale of PMS is shorter than that of ORC, which may be due to coupling effects and to the fact that the precipitation extremes in PMS are too strong . , in their analysis of climate change impacts on droughts, also highlight the fact that model structure is likely to contribute in an extensive way to different regional climate change projections.

The results obtained are coherent with those from , in which the soil moisture and hydrological drought analyses also showed differences among LSMs. Therefore, both studies provide clear proof that improvements in modelling concerning soil moisture and streamflow are needed.

7 Conclusions

In the context of a changing climate, it is necessary to evaluate the evolution of extremes, such as drought. Understanding the processes involved is therefore vital. To do so, the current modelling tools must first be evaluated. The work presented here analyses how RCMs represent meteorological, soil moisture, and hydrological drought and the propagation from a precipitation anomaly to soil moisture and streamflow anomalies.

Of the three RCMs analysed, RS4 shows the largest similarities with the reference data and deviates the least from its driving data in the meteorological and soil moisture drought representation analyses. However, none of the RCMs analysed can be clearly identified as showing the best performance with respect to hydrological drought representation. CL4 and RS4 show larger similarities with the reference data than the other RCMs across the Duero and Ebro basins, respectively. The propagation to soil moisture drought analysis indicates that RS4 reacts too quickly to precipitation anomalies, while PMS provides richer spatial patterns and a higher role of the continental surface. Finally, the propagation to hydrological drought analysis shows that RS4 and PMS react quickly to precipitation anomalies. These RCMs are in better agreement with the reference data than CL4 over the Ebro Basin, whereas CL4 is in better agreement with the reference data over the Duero Basin.

It is concluded that RCMs provide added value to meteorological drought representation, minimizing possible error sources from the driving data and ameliorating its characterization over areas that are known to pose certain problems to global driving data products. However, soil moisture and hydrological drought representation by RCMs show uncertainties. This is mainly due to the relevance of model physics and its prevalence to the driving data. Similar results were obtained for the propagation processes, in which model structure was found to influence the dynamics of drought propagation, showing different temporal scales depending on how precipitation variability is formulated within the model.

RCMs are a suitable tool for meteorological drought studies but should be used cautiously for soil moisture and hydrological drought analyses. Improvements regarding soil moisture modelling and streamflow-related processes (natural and anthropic) should be performed to better characterize drought events as well as their propagation.

Some prospective studies pertaining to this work could be to relate the results obtained to real drought impact data to study the relevance of the uncertainties found in the soil moisture and hydrological drought analyses. In addition, the analyses can be extended by including more RCMs and other drought indices. On the one hand, using more models would increase the information about drought simulation at a regional scale but also allow for the identification of further improvements in LSMs. On the other hand, if the analysis is performed with other indices, we can study the effects that other variables and processes, such as temperature and evapotranspiration, have on droughts. Studies of seasonal effects on droughts would also be interesting because these effects play an important role in drought propagation. Considering the soil moisture drought analysis, a comparison could be performed using in situ soil moisture data over a given point. Finally, it would be interesting to analyse the added value of drought indices as a function of the timescale used and how this may affect drought representation and its propagation.

Data availability
Data availability.

The forcing datasets and driving data used in this study can be accessed from their original source – ERA-Interim: https://www.ecmwf.int/en/forecasts/datasets/archive-datasets/reanalysis-datasets/era-interim (ECMWF2019); the SAFRAN dataset for Spain is available for research purposes from the MISTRALS HyMeX database . The RCM simulations were downloaded from the Med-CORDEX database: https://www.medcordex.eu/ (last access: 25 September 2017). The ORCHIDEE and SURFEX LSM simulations were produced for this study but can be reproduced using the corresponding release of the models: ORCHIDEE – https://forge.ipsl.jussieu.fr/orchidee (last access: 29 October 2019) (release no. 4676); SURFEX – https://www.umr-cnrm.fr/surfex/spip.php?rubrique (last access: 29 October 2019). The observed and modelled streamflow were provided by the Spanish Ministry for the Ecological Transition https://sig.mapama.gob.es/ (last access: 5 December 2019).

Author contributions
Author contributions.

ABO performed the computations and wrote the paper with support from PQS, who developed the code for the analytical calculations.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

The authors would like to thank Miguel Ángel Gaertner and Jan Polcher for their help. We are also very thankful to the reviewers whose comments have greatly improved the document. We thank the Spanish State Meteorological Agency (Agencia Estatal de Meteorología, AEMET) and the Spanish Ministry for the Ecological Transition for providing us with data. This work contributes to the HyMeX programme (Hydrological Cycle in the Mediterranean Experiment; http://www.hymex.org, last access: 29 October 2019).

Financial support
Financial support.

This research has been supported by the Ministerio de Ciencia, Innovación y Universidades (grant nos. CGL2013-47261-R and CGL2017-85687-R), and the European Regional Development Fund (grant no. EFA210/16/ PIRAGUA) in the framework of the INTERREG V-A España-Francia-Andorra programme.

Review statement
Review statement.

This paper was edited by Ryan Teuling and reviewed by two anonymous referees.

References

AEMET: Iberian Climate Atlas, Agencia Estatal de Meteorología, Madrid, Spain, 2011. a

AghaKouchak, A., Farahmand, A., Melton, F. S., Teixeira, J., Anderson, M. C., Wardlow, B. D., and Hain, C. R.: Remote sensing of drought: Progress, challenges and opportunities, Rev. Geophys., 53, 452–480, https://doi.org/10.1002/2014RG000456, 2015. a

Anagnostopoulou, C.: Drought episodes over Greece as simulated by dynamical and statistical downscaling approaches, Theor. Appl. Climatol., 129, 587–605, https://doi.org/10.1007/s00704-016-1799-5, 2017. a

Andreadis, K. and Lettenmaier, D.: Trends in 20th century drought over the continental United States, Geophys. Res. Lett., 33, L10403, https://doi.org/10.1029/2006GL025711, 2006. a

Balsamo, G., Albergel, C., Beljaars, A., Boussetta, S., Brun, E., Cloke, H., Dee, D., Dutra, E., Pappenberger, F., de Rosnay, P., Muñoz-Sabater, J., Stockdale, F., and Vitart, F.: ERA-Interim/Land: A global land-surface reanalysis based on ERA-Interim meteorological forcing, Tech. rep., ECMWF, 2012. a

Barker, L. J., Hannaford, J., Chiverton, A., and Svensson, C.: From meteorological to hydrological drought using standardised indicators, Hydrol. Earth Syst. Sci., 20, 2483–2505, https://doi.org/10.5194/hess-20-2483-2016, 2016. a

Barrera-Escoda, A., Goncalves, M., Guerreiro, D., Cunillera, J., and Baldasano, J.: Projections of temperature and precipitation extremes in the North Western Mediterranean Basin by dynamical downscaling of climate scenarios at high resolution (1971–2050), Climatic Change, 122, 567–582, https://doi.org/10.1007/s10584-013-1027-6, 2013. a

Belo-Pereira, M., Dutra, E., and Viterbo, P.: Evaluation of global precipitation data sets over the Iberian Peninsula, J. Geophys. Res., 116, D20101, https://doi.org/10.1029/2010JD015481, 2011. a, b

Beuvier, J., Sevault, F., Herrmann, M., Kontoyiannis, H., Ludwig, W., Rixen, M., Stanev, E., Béranger, K., and Somot, S.: Modeling the Mediterranean Sea interannual variability during 1961–2000: Focus on the Eastern Mediterranean Transient, J. Geophys. Res., 115, C08017, https://doi.org/10.1029/2009JC005950, 2010. a

Blenkinshop, S. and Fowler, H.: Changes in drought frequency, severity and duration for the British Isles projected by the PRUDENCE regional climate models, J. Hydrol., 342, 50–71, 2007a. a

Blenkinshop, S. and Fowler, H.: Changes in European drought characteristics projected by the PRUDENCE regional climate models, Int. J. Climatol., 27, 1595–1610, 2007b. a

Boone, A., Calvet, J. C., and Noilhan, J.: Inclusion of a Third Soil Layer in a Land Surface Scheme Using the Force–Restore Method, J. Appl. Meteorol., 38, 1611–1630, 1999. a

Bowden, J., Talgo, K., Spero, T., and Nolte, C.: Assessing the Added Value of Dynamical Downscaling Using the Standardized Precipitation Index, Adv. Meteorol., 2016, 8432064, https://doi.org/10.1155/2016/8432064, 2016. a

Castro, M., Fernández, C., and Gaertner, M. A.: Description of a Mesoscale Atmospheric Numerical Model, in: Mathematics, Climate and Environment, edited by: Diaz, J. I. and Lions, J. L., Rech. Math. Appl. Ser. Mason, 55, 230–253, 1993. a, b

Christensen, J., Carter, T., Rummukainen, M., and Amanatidis, G.: Evaluating the performance and utility of regional climate models: the PRUDENCE approach, Climatic Change, 81, 1–6, 2007. a

Colin, J., Déqué, M., Radu, R., and Somot, S.: Sensitivity study of heavy precipitations in Limited Area Model climate simulations: influence of the size of the domain and the use of the spectral nudging technique, Tellus A, 62, 591–604, https://doi.org/10.1111/j.1600-0870.2010.00467.x, 2010. a

Coll-Pajarón, M.: Distribución de la humedad del suelo mediante observaciones del satélite SMOS, modelización con SURFEX y medidas in situ sobre la Valencia Anchor Station, Phd, Universitat de València, available at: http://roderic.uv.es/handle/10550/60877 (last access: 15 December 2018), 2017. a

Decharme, B., Alkama, R., Douville, H., Becker, M., and Cazenave, A.: Global Evaluation of the ISBA-TRIP Continental Hydrological System. Part II: Uncertainties in River Routing Simulation Related to Flow Velocity and Groundwater Storage, J. Hydrometeor., 11, 601–617, https://doi.org/10.1175/2010JHM1212.1, 2010. a

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, https://doi.org/10.1002/qj.828, 2011. a

del Río, S., Herrero, L., Pinto-Gomes, C., and Penas, A.: Spatial analysis of mean temperature trends in Spain over the period 1961–2006, Global Planet. Change, 78, 65–75, https://doi.org/10.1016/j.gloplacha.2011.05.012, 2011. a

de Luis, M., Brunetti, M., González-Hidalgo, J. C., Longares, L. A., and Martín-Vide, J.: Changes in seasonal precipitation in the Iberian Peninsula during 1946–2005, Global Planet. Change, 74, 27–33, https://doi.org/10.1016/j.gloplacha.2010.06.006, 2010. a

Déqué, M. and Somot, S.: Analysis of heavy precipitation for France using ALADIN RCM simulations, Idöjaras Q. J. Hungarian Meteorological Service, 112, 179–190, 2008. a

de Rosnay, P. and Polcher, J.: Modelling root water uptake in a complex land surface scheme coupled to a GCM, Hydrol. Earth Syst. Sci., 2, 239–255, https://doi.org/10.5194/hess-2-239-1998, 1998. a, b

Domínguez, M., Gaertner, M. A., de Rosnay, P., and Losada, T.: A regional climate model simulation over West Africa: Parameterization tests and analysis of land-surface fields, Clim. Dynam., 35, 249–265, 2010. a, b

Domínguez, M., Romera, R., Sánchez, E., Fita, L., Fernández, J., Jiménez-Guerrero, P., Montávez, J., Cabos, W., Liguori, G., and Gaertner, M.: Present-climate precipitation and temperature extremes over Spain from a set of high resolution RCMs, Clim. Res., 58, 149–164, https://doi.org/10.3354/cr01186, 2013. a

Durand, Y., Brun, E., Merindol, L., Guyomarc'h, G., Lesaffre, B., and Martin, E.: A meteorological estimation of relevant parameters for snow models, Ann. Glaciol., 18, 65–71, 1993. a

Durand, Y., Giraud, G., Brun, E., Merindol, L., and Martin, E.: A computer-based system simulating snowpack structures as a tool for regional avalanche forecasting, J. Glaciol., 45, 469–484, 1999. a

ECMWF: ERA-Interim, available at: https://www.ecmwf.int/en/forecasts/datasets/archive-datasets/reanalysis-datasets/era-interim, last access: 29 October 2019. a

Edossa, D., Babel, M., and Das Gupta, A.: Drought Analysis in the Awash River Basin, Ethiopia, Water Resour. Manag., 24, 1441–1460, https://doi.org/10.1007/s11269-009-9508-0, 2009. a

Escorihuela, M. J. and Quintana-Seguí, P.: Comparison of remote sensing and simulated soil moisture datasets in Mediterranean landscapes, Remote Sens. Environ., 180, 99–114, https://doi.org/10.1016/j.rse.2016.02.046, 2016. a

Estrela, T. and Quintas, L.: El sistema integrado de modelización precipitación-aportación SIMPA, Revista de Ingeniería Civil, 104, 43–52, 1996. a, b

Evans, J.: Straightforward Statistics for the Behavioral Sciences, Pacific Grove, CA, Brooks/Cole Publishing, 1996. a, b, c, d

FAO: The State of Food Insecurity in the World, Rome, Italy, FAO, 2009. a

Farahmand, A. and AghaKouchak, A.: A generalized framework for deriving nonparametric standardized drought indicators, Adv. Water Resour., 76, 140–145, https://doi.org/10.1016/j.advwatres.2014.11.012, 2015. a, b

Farda, A., Déqué, M., Somot, S., Horányi, A., Spiridonov, V., and Tóth, H.: Model ALADIN as regional climate model for Central and Eastern Europe, Stud. Geophys. Geod., 54, 313–332, https://doi.org/10.1007/s11200-010-0017-7, 2010. a

Faroux, S., Kaptué Tchuenté, A. T., Roujean, J.-L., Masson, V., Martin, E., and Le Moigne, P.: ECOCLIMAP-II/Europe: a twofold database of ecosystems and surface parameters at 1 km resolution based on satellite information for use in land surface, meteorological and climate models, Geosci. Model Dev., 6, 563–582, https://doi.org/10.5194/gmd-6-563-2013, 2013. a

Feser, F., Rockel, B., von Storch, H., Winterfeldt, J., and Zahn, M.: Regional Climate Models Add Value to Global Model Data: A Review and Selected Examples, B. Am. Meteorol. Soc., 92, 1181–1192, https://doi.org/10.1175/2011BAMS3061.1, 2011. a

Gandin, L. S.: Objective analysis of meteorological fields. By L. S. Gandin. Translated from the Russian. Jerusalem (Israel Program for Scientific Translations), Q. J. Roy. Meteor. Soc., 92, 447–447, https://doi.org/10.1002/qj.49709239320, 1966. a

García-Valdecasas Ojeda, M., Gámiz-Fortis, S. R., Castro-Díez, Y., and Esteban-Parra, M.: Evaluation of WRF capability to detect dry and wet periods in Spain using drought indices, J. Geophys. Res.-Atmos., 122, 1569–1594, https://doi.org/10.1002/2016JD025683, 2017. a, b, c

Giorgi, F., Jones, C., and Asrar, G.: Addressing climate information needs at the regional level: The CORDEX framework, WMO Bull, 58, 175–183, 2009. a

Gupta, H. V., Kling, H., Yilmaz, K. K., and Martínez, G. F.: Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling, J. Hydrol., 377, 80–91, https://doi.org/10.1016/j.jhydrol.2009.08.003, 2009. a

Herrmann, M., Somot, S., Calmanti, S., Dubois, C., and Sevault, F.: Representation of spatial and temporal variability of daily wind speed and of intense wind events over the Mediterranean Sea using dynamical downscaling: impact of the regional climate model configuration, Nat. Hazards Earth Syst. Sci., 11, 1983–2001, https://doi.org/10.5194/nhess-11-1983-2011, 2011. a

Hoerling, M., Eischeid, J., Perlwitz, J., Quan, X., Zhang, T., and Pegion, P.: On the Increased Frequency of Mediterranean Drought, J. Climate, 25, 2146–2161, https://doi.org/10.1175/JCLI-D-11-00296.1, 2012. a

IPCC: Climate Change 2007: Impacts, Adaptation and Vulnerability, in: Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Parry, M. L., Canziani, O. F., Palutikof, J. P., van der Linden, P. J., and Hanson, C. E., Cambridge University Press, Cambridge, UK, 976 pp., 2007. a

IPCC: Climate change 2014: Synthesis Report, in: Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Core Writing Team, Pachauri, R. K., and Meyer, L. A., IPCC, Geneva, Switzerland, 151 pp., 2014. a

Jacobsen, I. and Heise, E.: A new economic method for the computation of the surface temperature in numerical models, Beitr. Phys. Atm., 55, 128–141, 1982. a

Jenkins, K. and Warren, R.: Quantifying the impact of climate change on drought regimes using the Standardised Precipitation Index, Theor. Appl. Climatol., 120, 41–54, https://doi.org/10.1007/s00704-014-1143-x, 2014. a

Jiménez-Guerrero, P., Montávez, J. P., Domínguez, M., Romera, R., Fita, L., Fernández, J., Cabos, W. D., Liguori, G., and Gaertner, M. A.: Description of mean fields and interannual variability in an ensemble of RCM evaluation simulations over Spain: results from the ESCENA project, Clim. Res., 57, 201–220, https://doi.org/10.3354/cr01165, 2013. a

Kaptue Tchuente, A. T., Roujean, J. L., and Faroux, S.: ECOCLIMAP-II: An ecosystem classification and land surface parameters database of Western Africa at 1 km resolution for the African Monsoon Multidisciplinary Analysis (AMMA) project, Remote Sens. Environ., 114, 961–976, https://doi.org/10.1016/j.rse.2009.12.008, 2010. a

Kenawy, A., López-Moreno, J. I., and Vicente-Serrano, S. M.: Summer temperature extremes in northeastern Spain: Spatial regionalization and links to atmospheric circulation (1960–2006), Theor. Appl. Climatol., 113, 387–405, https://doi.org/10.1007/s00704-012-0797-5, 2013. a

Koster, R. D., Guo, Z., Yang, R., Dirmeyer, P. A., Mitchell, K., and Puma, M. J.: On the Nature of Soil Moisture in Land Surface Models, J. Climate, 22, 4322–4335, https://doi.org/10.1175/2009JCLI2832.1, 2009. a

Krinner, G., Viovy, N., de Noblet-Ducoudré, N., Ogée, J., Polcher, J., Friedlingstein, P., Ciais, P., Sitch, S., and Prentice, I.: A dynamic global vegetation model for studies of the coupled atmosphere-biosphere system, Global Biogeochem. Cy., 19, GB1015, https://doi.org/10.1029/2003GB002199, 2005. a, b

Li, Z. X.: Ensemble Atmospheric GCM Simulation of Climate Interannual Variability from 1979 to 1994, J. Climate, 12, 986–1001, 1999. a

López-Bustins, J. A., Pascual, D., Pla, E., and Retana, J.: Future variability of droughts in three Mediterranean catchments, Nat. Hazards, 69, 1405–1421, https://doi.org/10.1007/s11069-013-0754-3, 2013. a

Mariotti, A.: Recent Changes in the Mediterranean Water Cycle: A Pathway toward Long-Term Regional Hydroclimatic Change?, J. Climate, 23, 1513–1525, https://doi.org/10.1175/2009JCLI3251.1, 2010. a

Martínez-Fernández, J., Sánchez, N., and Herrero-Jiménez, C. M.: Recent trends in rivers with near-natural flow regime: The case of the river headwaters in Spain, Prog. Phys. Geog., 37, 685–700, https://doi.org/10.1177/0309133313496834, 2013. a

Masson, V., Champeaux, J. L., Chauvin, F., Meriguet, C., and Lacaze, R.: A Global Database of Land Surface Parameters at 1-km Resolution in Meteorological and Climate Models, J. Climate, 16, 1261–1282, https://doi.org/10.1175/1520-0442-16.9.1261, 2003. a

Masud, M. B., Khaliq, M. N., and Wheater, H. S.: Future changes to drought characteristics over the Canadian Prairie Provinces based on NARCCAP multi-RCM ensemble, Clim. Dynam., 48, 2685–2705, https://doi.org/10.1007/s00382-016-3232-2, 2017. a

Maule, C. F., Thejll, P., Christensen, J. H., Svendsen, S. H., and Hannaford, J.: Improved confidence in regional climate model simulations of precipitation evaluated using drought statistics from the ENSEMBLES models, Clim. Dynam., 40, 155–173, https://doi.org/10.1007/s00382-012-1355-7, 2013. a, b

McKee, T. B., Doesken, N. J., and Kleist, J.: The Relationship of Drought Frequency and Duration Time Scales, Proceedings of the 8th Conference on Applied Climatology, 17–22 January 1993, American Meteorological Society Boston, MA, USA, vol. 17, 179–183, 1993. a, b

Meresa, H., Osuch, M., and Romanowicz, R.: Hydro-Meteorological Drought Projections into the 21-st Century for Selected Polish Catchments, Water, 8, 206, https://doi.org/10.3390/w8050206, 2016. a, b

Mishra, A. K. and Singh, V. P.: A review of drought concepts, J. Hydrol., 391, 202–216, https://doi.org/10.1016/j.jhydrol.2010.07.012, 2010. a, b

Mishra, A. and Singh, V. P.: Drought modeling – A review, J. Hydrol., 403, 157–175, https://doi.org/10.1016/j.jhydrol.2011.03.049, 2011. a

Morán-Tejeda, E., Herrera, S., López-Moreno, J. I., Revuelto, J., Lehmann, A., and Beniston, M.: Evolution and frequency (1970–2007) of combined temperature–precipitation modes in the Spanish mountains and sensitivity of snow cover, Reg. Environ. Change, 13, 873–885, https://doi.org/10.1007/s10113-012-0380-8, 2013. a

Nabat, P. S., Somot, S., Mallet, M., Sánchez-Lorenzo, A., and Wild, M.: Contribution of anthropogenic sulfate aerosols to the changing Euro-Mediterranean climate since 1980, Geophys. Res. Lett., 41, 5605–5611, https://doi.org/10.1002/2014GL060798, 2014. a, b

Ngo-Duc, T., Polcher, J., and Laval, K.: A 53-year forcing data set for land surface models, J. Geophys. Res., 110, D06116, https://doi.org/10.1029/2004JD005434, 2005. a

Noilhan, J. and Mahfouf, J. F.: The ISBA land surface parameterisation scheme, Global Planet. Change, 13, 145–159, https://doi.org/10.1016/0921-8181(95)00043-7, 1996. a, b

Noilhan, J. and Planton, S.: A Simple Parameterization of Land Surface Processes for Meteorological Models, Mon. Weather Rev., 117, 536–549, 1989. a, b

Olcina, J.: Tipología de sequías en España, Ería, 56, 201–227, 2001. a

Owens, T., Hoddinott, J., and Kinsey, B.: Ex-Ante Actions and Ex-Post Public Responses to Drought Shocks: Evidence and Simulations from Zimbabwe, World Dev., 31, 1239–1255, https://doi.org/10.1016/S0305-750X(03)00068-8, 2003. a

PaiMazumder, D., Sushama, L., Laprise, R., Khaliq, M., and Sauchyn, D.: Canadian RCM projected changes to short and long term drought characteristics over the Canadian Prairies, Int. J. Climatol., 33, 1409–1423, https://doi.org/10.1002/joc.3521, 2013. a

Palmer, W. C.: Meteorological Drought, vol. 30, US Department of Commerce, Weather Bureau Washington, DC, 1965. a

Post, W. and Zobler, L.: Global soil types, 0.5-degree grid (modified zobler), Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, USA, https://doi.org/10.3334/ORNLDAAC/540, 2000. a

Quintana-Seguí, P.: SAFRAN analysis over Spain, ESPRI/IPSL, https://doi.org/10.14768/MISTRALS-HYMEX.1388, 2015. a, b

Quintana-Seguí, P., Peral, M. C., Turco, M., Llasat, M. C., and Martin, E.: Meteorological Analysis Systems in North-East Spain: Validation of SAFRAN and SPAN, J. Environ. Inform., 27, 116–130, https://doi.org/10.3808/jei.201600335, 2016. a

Quintana-Seguí, P., Turco, M., Herrera, S., and Miguez-Macho, G.: Validation of a new SAFRAN-based gridded precipitation product for Spain and comparisons to Spain02 and ERA-Interim, Hydrol. Earth Syst. Sci., 21, 2187–2201, https://doi.org/10.5194/hess-21-2187-2017, 2017. a, b

Quintana-Seguí, P., Barella-Ortiz, A., Regueiro-Sanfiz, S., and Míguez-Macho, G.: The Utility of Land-Surface Model Simulations to Provide Drought Information in a Water Management Context Using Global and Local Forcing Datasets, Water Resour. Manag., https://doi.org/10.1007/s11269-018-2160-9, online first, 2019. a, b, c, d, e

Radu, R., Déqué, M., and Somot, S.: Spectral nudging in a spectral regional climate model, Tellus A, 60, 898–910, https://doi.org/10.1111/j.1600-0870.2008.00341.x, 2008. a

Reichstein, M., Bahn, M., Ciais, P., Frank, D., Mahecha, M., Seneviratne, S., Zscheischler, J., Beer, C., Buchmann, N., Frank, D., Papale, D., Rammig, A., Smith, P., Thonicke, K., van der Velde, M., Vicca, S., Walz, A., and Wattenbach, M.: Climate extremes and the carbon cycle, Nature, 500, 287–295, https://doi.org/10.1038/nature12350, 2013. a

Ritter, B. and Geleyn, J. F.: A comprehensive radiation scheme for numerical weather prediction models with potential applications in climate simulations, Mon. Weather Rev., 120, 303–325, https://doi.org/10.1175/1520-0493(1992)1200303:ACRSFN2.0.CO;2, 1992. a

Rockel, B., Will, A., and Hense, A.: The Regional Climate Model COSMO-CLM (CCLM), Meteorol. Z., 17, 347–348, https://doi.org/10.1127/0941-2948/2008/0309, 2008. a, b

Ruiz, J. M.: Modelo distribuido para la evaluación de recursos hídricos (modelo SIMPA), Madrid, CEDEX, Ministerio de Fomento Doctorado 180 pp., 1999. a, b

Ruti, P., Somot, S., Giorgi, F., Dubois, C., Flaounas, E., Obermann, A., Dell'Aquila, A., Pisacane, G., Harzallah, A., Lombardi, E., Ahrens, B., Akhtar, N., Alias, A., Arsouze, T., Aznar, R., Bastin, S., Bartholy, J., Béranger, K., Beuvier, J., Bouffies-Cloché, S., Brauch, J., Cabos, W., Calmanti, S., Calvet, J., Carillo, A., Conte, D., Coppola, E., Djurdjevic, V., Drobinski, P., Elizalde-Arellano, A., Gaertner, M., Galán, P., Gallardo, C., Gualdi, S., Goncalves, M., Jorba, O., Jordà, G., L’Heveder, B., Lebeaupin-Brossier, C., Li, L., Liguori, G., Lionello, P., Maciàs, D., Nabat, P., Onol, B., Raikovic, B., Ramage, K., Sevault, F., Sannino, G., Struglia, M., Sanna, A., Torma, C., and Vervatis, V.: Med-CORDEX Initiative for Mediterranean Climate Studies, B. Am. Meteorol. Soc., 97, 1187–1208, https://doi.org/10.1175/BAMS-D-14-00176.1, 2016. a

Sánchez, E., Gallardo, C., Gaertner, M., Arribas, A., and Castro, M.: Future climate extreme events in the Mediterranean simulated by a regional climate model: a first approach, Global Planet. Change, 44, 163–180, 2004. a, b

Serrano, A., García, J. A., Mateos, V. L., Cancillo, M. L., and Garrido, J.: Monthly modes of variation of precipitation over the Iberian peninsula, J. Climate, 12, 2894–2919, 1999. a

Sevault, F., Somot, S., Alias, A., Dubois, C., Lebeaupin-Brossier, C., Nabat, P., Adloff, F., Déqué, M., and Decharme, B.: A fully coupled Mediterranean regional climate system model: design and evaluation of the ocean component for the 1980–2012 period, Tellus A, 66, 23967, https://doi.org/10.3402/tellusa.v66.23967, 2014. a, b

Shukla, S. and Wood, A. W.: Use of a standardized runoff index for characterizing hydrologic drought, Geophys. Res. Lett., 35, 2405, https://doi.org/10.1029/2007GL032487, 2008 a

Sousa, P. M., Trigo, R. M., Aizpurua, P., Nieto, R., Gimeno, L., and Garcia-Herrera, R.: Trends and extremes of drought indices throughout the 20th century in the Mediterranean, Nat. Hazards Earth Syst. Sci., 11, 33–51, https://doi.org/10.5194/nhess-11-33-2011, 2011. a, b

Sylla, M. B., Coppola, E., Mariotti, L., Giorgi, F., Ruti, P. M., Dell'Aquila, A., and Bi, X.: Multiyear simulation of the African climate using a regional climate model (RegCM3) with the high resolution ERA-interim reanalysis, Clim. Dynam., 35, 231–247, https://doi.org/10.1007/s00382-009-0613-9, 2010. a

Turco, M. and Llasat, M. C.: Trends in indices of daily precipitation extremes in Catalonia (NE Spain), 1951–2003, Nat. Hazards Earth Syst. Sci., 11, 3213–3226, https://doi.org/10.5194/nhess-11-3213-2011, 2011. a

van der Linden, P. and Mitchell, J. F. B. (Eds.): ENSEMBLES: Climate Change and its Impacts: Summary of research and results from the ENSEMBLES project, MetOffice Hadley Centre, FitzRoy Road, Exeter EX1 3PB, UK, 160 pp., 2009. a

Van Loon, A.: Hydrological drought explained, WIREs Water, 2, 359–392, https://doi.org/10.1002/wat2.1085, 2015. a

Van Loon, A. F., Van Lanen, H. A. J., Hisdal, H., Tallaksen, L. M., Fendekova,M., Oosterwijk, J., Horvat, O., and Machlica, A.: Understanding hydrological winter drought in Europe, Global Change: Facing Risks and Threats to Water Resources, 189–197, Centre for Ecology and Hydrology, Wallingford, Oxfordshire, OX10 8BB, UK, IAHS, 2010. a

Van Loon, A. F. and Van Lanen, H. A. J.: A process-based typology of hydrological drought, Hydrol. Earth Syst. Sci., 16, 1915–1946, https://doi.org/10.5194/hess-16-1915-2012, 2012a. a

Van Loon, A. F., Van Huijgevoort, M. H. J., and Van Lanen, H. A. J.: Evaluation of drought propagation in an ensemble mean of large-scale hydrological models, Hydrol. Earth Syst. Sci., 16, 4057–4078, https://doi.org/10.5194/hess-16-4057-2012, 2012b. a

Vicente-Serrano, S.: Differences in spatial patterns of drought on different time scales: An analysis of the Iberian Peninsula, Water Resour. Manage., 20, 37–60, https://doi.org/10.1007/s11269-006-2974-8, 2006. a

Vicente-Serrano, S. M. and López-Moreno, J. I.: Hydrological response to different time scales of climatological drought: an evaluation of the Standardized Precipitation Index in a mountainous Mediterranean basin, Hydrol. Earth Syst. Sci., 9, 523–533, https://doi.org/10.5194/hess-9-523-2005, 2005. a

Vicente-Serrano, S., López-Moreno, J. I., Beguería, S., Lorenzo-Lacruz, J., Azorín-Molina, C., and Morán-Tejeda, E.: Accurate Computation of a Streamflow Drought Index, J. Hydrol. Eng., 17, 318–332, 2012.  a

Vicente-Serrano, S., Gouveia, C., Camarero, J. J., Beguería, S., Trigo, R., López-Moreno, J., Azorín-Molina, C., Pasho, E., Lorenzo-Lacruz, J., Revuelto, J., Morán-Tejeda, E., and Sánchez-Lorenzo, A.: Response of vegetation to drought timescales across global land biomes, P. Natl. Acad. Sci. USA, 110, 52–57, https://doi.org/10.1073/pnas.1207068110, 2013. a

Vicente-Serrano, S., López-Moreno, J. I., Beguería, S., Lorenzo-Lacruz, J., Sánchez-Lorenzo, A., García-Ruiz, J. M., Azorín-Molina, C., Morán-Tejeda, E., Revuelto, J., Trigo, R., Coelho, F., and Espejo, F.: Evidence of increasing drought severity caused by temperature rise in southern Europe, Environ. Res. Lett., 9, 44001, https://doi.org/10.1088/1748-9326/9/4/044001, 2014. a

Vu, M., Raghavan, V., and Liong, S.: Ensemble Climate Projection for Hydro-Meteorological Drought over a river basin in Central Highland, Vietnam, Journal of Civil Engineering, 19, 427–433, https://doi.org/10.1007/s12205-015-0506-x, 2015. a, b

Wang, D., Hejazi, M., Cai, X., and Valocchi, A. J.: Climate change impact on meteorological, agricultural, and hydrological drought in central Illinois, Water Resour. Res., 47, W09527, https://doi.org/10.1029/2010WR009845, 2011. a, b

Weedon, G., Gomes, S., Viterbo, P., Shuttleworth, W., Blyth, E., Osterle, H., Adam, J. C., Bellouin, N., Boucher, O., and Best, M.: Creation of the WATCH Forcing Data and Its Use to Assess Global and Regional Reference Crop Evaporation over Land during theTwentieth Century, J. Hydrometeorol., 12, 823–848, https://doi.org/10.1175/2011JHM1369.1, 2011. a

Wilhite, D. A.: Drought: a global assessment, Vols I and II, Routledge Hazards and Disasters Series, Routledge, London, UK, 2000. a

Wu, C., Xian, Z., and Huang, G.: Meteorological drought in the Beijiang River basin, South China: current observations and future projections, Stoch, Environ. Res. Risk Assess., 30, 1821–1834, https://doi.org/10.1007/s00477-015-1157-7, 2016. a

Xu, Y., Wang, L., Ross, K. W., Liu, C., and Berry, K.: Standardized soil moisture index for drought monitoring based on soil moisture active passive observations and 36 years of north American land data assimilation system data: A case study in the southeast United States, Remote. Sens., 10, 301, https://doi.org/10.3390/rs10020301, 2018. a