3.1 Historical climate
For historical climate, both SPI and SPEI (Fig. 1a, top panel) show similar spatial patterns for decadal frequency of moderate droughts on all timescales, with SPEI revealing slightly higher frequencies for the 1- month timescale. In the case of severe droughts (Fig. 1a, middle panel) SPI and SPEI spatial patterns resemble a lot, showing slightly larger values for SPEI than SPI for most of the temporal scales (except for the 6-months). In the case of severe droughts (Fig. 1a, middle panel), SPI shows slightly higher decadal frequencies on all time scales. As expected, smaller timescales show larger decadal frequencies, independently of the analysed index. When analysing both SPI and SPEI for the same timescale, decadal frequencies are larger on the western coast and north-western areas. The frequencies of moderate drought peak at decadal values above 26 events for the one-month scale in the northwest, and at the 2-years scale values ranging mostly between 4 and 6. For extreme droughts, the decadal frequency of events broadly varies from 5 events at the 3-monthly scale to above 2 events at the yearly scale.
During the 1971–2000 period, droughts show higher values of duration when computed using SPEI (Fig. 1b), with a similar spatial characterization relative to SPI. Droughts with longer durations are found in the Iberian north-eastern and south-eastern regions, reaching 14 months long in some areas for moderate droughts and above 8 months for severe drought. Those durations refer to the 24-month accumulation drought time scales.
In terms of the spatial extent of drought for the 1971–2000 period, all accumulation time scales reveal similar distributions in each moderate, severe, or extreme drought type. For the SPEI, the distribution has slightly more variability reaching 80% of the Iberian Peninsula in drought for the 90th percentile (high whisker). In general, the drought spatial extent is lower than 20% and 10% of the Iberian territory for severe and extreme droughts, respectively.
3.2 Multiple scenarios
The SPI and SPEI indexes in this section are built from a 300-year period, thus featuring the historical 1971–2000 and the 2011–2100 from all RCP scenarios. Thus, differences arise relative to the results shown in the previous section. Figure S1 displays those differences for the 1971–2000 time slice, which stand for the historical period. Overall and as expected, the reference period with the 300-year indexes reveals a more humid panorama, underestimating droughts on the current climate. For instance, Figure S5a displays the differences of the mean value, with higher values among all timescales, namely for the higher ones. Figure S5b and S5c shows negative differences for the decadal frequency, and mean duration of droughts. Those differences imply less frequent and shorter events for the reference period from the 300-year indexes in comparison with the 30-year indexes. These differences are more evident for the shorter timescales for the decadal frequency, while for the mean drought duration the higher differences occur for the longer timescales.
Regarding future drought intensity (Fig. 2) both indices reveal a clear shift to larger drought intensities when considering the RCP2.6 to RCP4.5 and RCP8.5. This shift agrees with the known future projections of rainfall reduction and temperature increase, with the latter included on SPEI for Iberia (Soares et al. 2017a; Cardoso et al. 2019; Argueso et al. 2012; Cos et al. 2022). Moreover, responding to the projected rising temperatures, SPEI indicates larger anomalies of intensities than SPI, as expected. A clear shift with time throughout the century for higher negative values is not seen within all RCPs. In fact, this only occurs for RCP8.5.
Looking at the anomalies of drought intensities for the RCP2.6 (Fig. 2a), SPI and SPEI points to an increase until the mid-century and a reduction for the end-of-century of drought intensity in some areas of the Iberian Peninsula and for all drought timescales. These positive anomalies are considerably higher for SPEI than SPI. Moreover, in the case of SPI for 2071–2100 the projected positive anomalies correspond to an attenuation of drought intensities. SPEI anomalies show a small increase of drought intensity for all periods, but being milder for the end of the century, in general smaller than − 0.5.
For the intermediate scenario (RCP4.5, Fig. 2b), drought intensity time evolution shows a monotonous increase for SPEI, but with similar changes for mid and end century. For the shorter drought scales (1 to 6 months) those increases are up to -0.75 confined mostly to the southern areas. For the longer droughts scales (12 to 24 months), the anomalies reach a much larger value, up to -1 for SPEI, in a large extension of southern Iberia, for the end of the century. The differences between RCP2.6 and RCP4.5, for one- and two-year(s) droughts (12 and 24 months), at the end of the century are notorious. For SPEI, in the former case, drought intensities are projected to reach at most − 0.75 in small areas scattered in Iberia, and in the latter case anomalies above those values can reach − 1.25 over large portions of the Iberian mainland. Examining the RCP8.5 results, the monotonous increase in drought intensity is greatly enhanced with the increased drought timescale and throughout the century, according to both indices. SPEI for 12- and 24-months scales point to intensity anomalies rather severe, which attain values above − 1.5 in 2071–2100, corresponding to approximately 5x with respect to the 2011–2040 period. By far, the results for the RCP8.5 are much more severe when compared to the other scenarios, even for mid-century. In line with the major differences on the reported drought intensity increases among two of three RCPs (RCP4.5 and RCP8.5), the drought decadal intensity evolution for the 21st century is extremely worrying.
Figure 3 reveals the projections for the future anomalies of drought decadal frequencies for severe drought in Iberia. Similarly, Figures S1 and S2 show the drought decadal frequencies, respectively, for the projected anomalies of moderate and extreme drought. The examination of the SPI and SPEI results for all RCPs show an aggravation of the decadal occurrence of drought events at all time scales but particularly at the shorter ones. In the high-end case, for the RCP8.5 and at the end of the century, the increases surpass 15 (9) occurrences in a decade on the 1-month time scale as indicated by SPEI (SPI). These are followed closely by the other time scales, from 3- to 24-months, for which a troubling growth of more than 9 down to 3 drought events per decade are projected.
It is highly relevant to acknowledge the large discrepancy on the anomalies of drought decadal frequencies from RCP8.5 to RCP2.6, which strikingly point to completely different impacts and adaptation challenges. Moreover, SPI and SPEI show different increments, with the differences of severe drought decadal frequency depicted by SPEI being far larger than SPI, which are a result of SPEI trends being generally larger due to the regional warming effect (Páscoa et al. 2021). Based on the SPEI index which presents the mild-case scenario when it comes to the analysed indices (Figure 3), within the RCP2.6 an increase of severe drought is projected for the near- and mid-century followed by a reduction at the end of century, for all drought timescales. The mid-century increase is particularly large for the 1-month scale, with the difference being of around 6 more events per decade for most of Iberia. From 3 months to 2 years, increases between 1 and 3 events are projected for a large extent of the peninsula. For the 2071-2100 period a slight reduction in the increase of drought events is revealed, with 4 to 6 more events over large areas of the southern Peninsula at the 1-month timescale and an increase of at most 2 events for the other timescales. For the RCP4.5 scenario changes in severe drought frequency are projected to grow throughout the century reaching values above 9 for the 1-month scale. This increase is progressive throughout the 21st century. For the 3- to 12-months scales, the future anomalies reach above 4 events in large extensions of Iberia. These augments reflect a massive growth of drought occurrence in Iberia, which corresponds to multiplying by 2 and 3 the decadal frequency of severe droughts for the mid- and the end-of-century, respectively, when compared to the near future climate, for vast areas of Iberia.
For moderate drought the projected increases are larger for short drought scales, such as 1 and 3 months, reaching anomalies around 10 for the RCP8.5 and end-of-century (Figure S2). Under the RCP2.6, a slight variation in the decadal frequency of moderate droughts is observed throughout the 21st century, up to 3 (9) more events for SPI (SPEI). For the RCP4.5 and RCP8.5 a rising in decadal frequency is presented from mid-century, stronger for the non-mitigation scenario, from 3 (3) up to 9 (12) more moderate droughts considering the SPI (SPEI). In what concerns the projected evolution of extreme droughts, the results follow the ones shown for severe drought, but with maximum anomaly values in the range of 4 in RCP2.6, 6 in RCP4.5 and 9 in RCP8.5 considering the SPI (Figure S3).
Regarding the changes on severe drought mean event duration, for the three periods and the three RCPs, exceptional differences amongst emission scenarios and periods are depicted in Fig. 4. The changes are characterised generally by being greater for the longer timescales for both SPI and SPEI, with SPEI showing larger increments, independently of the emission scenario and for the three time periods. For RCP2.6, SPI reveals rather small increases or decreases for the severe drought durations, particularly at the end of the century. In the meanwhile, SPEI reveals at the most increases which can reach 5 months for the 24-month accumulation period. Moreover, for the RCP2.6, similar anomalies are displayed for the near- and mid-century cases, and even smaller amounts for the end of the century. For RCP4.5 the highest mean duration differences are seen for the mid-century, reaching values higher than 10 months in the centre of Iberian Peninsula for the SPEI 24-month timescale. For the lower drought timescales, the increase is not so intense, ranging from no change (1 to 3 months) to 3–5 months for the scales between 6-months and 1-year. or the beginning and end of the century, SPI shows heterogeneous anomalies, with some areas of the Peninsula revealing increases, while others display decreases. Conversely, SPEI mostly shows a consistent increase of mean duration for all the timescales both in the beginning and end of the century, although for the end of the century milder changes are expected, ranging from no change at the shorter timescales of 1 to 3 months and approximately an increase of 3–5 months for the 24-months accumulation, in comparison to the mid-century. Finally, regarding the RCP8.5, both SPI and SPEI generally show a consistent and important increase of mean duration for the longer timescale throughout the 21st century. These longer lasting drought increments change rather monotonously. For the 6-month drought timescale the mean durations are added up to 5 months, the one-year up to 15 and the 24-month above the 40 months over large areas of the domain.
Regarding the moderate and extreme drought mean event duration (Figures. S4 and S5, respectively), the change patterns follow the ones observed for severe drought events. For moderate drought, differences in mean duration are quite significant for the RCP4.5 and RCP8.5 from mid- to end-of-century. For the RCP4.5, the mean event duration can reach more than 10 (15) months in the mid-century 12-month timescale (24-month timescale), whilst for the RCP8.5 the moderate drought events can reach increases in average duration of 20 (40) months for the 24-month timescale considering the SPI (SPEI). Notice that the projected changes in average duration of moderate drought events is more pronounced in longer timescales than for shorter ones. For the RCP2.6 scenario, significant differences are found for the mid-century and for longer timescales, reaching more 5 (10) months in event duration for SPI (SPEI) considering the 24-month accumulation period. Similarly, it is expected an increase in average duration of extreme drought events, reaching more 5 months in RCP2.6, 15 months in both RCP4.5 and RCP8.5 considering the SPI and 1-month timescale.
Figure 5 displays the evolution of the spatial extent (in percentage of mainland Iberia) for moderate, severe, and extreme droughts, only for SPEI at the 12-month timescale. The projections for all the other accumulations and for the SPI index reveal similar results and similar conclusions may be inferred. As expected, the percentage of area in drought decreases from moderate to extreme droughts. Focusing on severe droughts (Figure 5b) the projected changes are in line with the previous results (Figure 2 to 4). The differences of the future changes amongst the three emission scenarios are striking. For the RCP2.6, the boxplots reveal higher percentage of area in drought at the beginning of the century and progressively reducing towards the end, reaching similar values relatively to the historical reference period reference period. As for the RCP4.5 the spatial variability peak occurs at the mid-century, slightly decreasing for the end of the century. As expected, it is for the RCP8.5 where the largest changes in drought spatial extension are projected, with a progressively increase towards the end of the 21st century, surpassing in mean 50 % of area in severe drought. For comparison, the high whisker for the historical reference period do not reach the 20% of area in drought, while for the end century RCP 8.5, the low (high) whisker has values of 10% (above 90%).
For the moderate and extreme droughts (Figure 5a and 5c), similar future spatial extensions can be found. Notice however, at the end-of-century, for the RCP8.5 scenario, the mean of the percentage of area in moderate drought can surpass the 70 % of the territory. Overall, the tendency of changes for moderate drought is similar as for severe drought. Considering the extreme drought, also for the end-of-century RCP8.5 scenario, the mean reveals more than 30% of land in extreme drought, contrasting with the other periods and scenarios. Indeed, this increase in area is noteworthy, even when comparing with the moderate and severe droughts.
Finally, Figures S6 to S12 show a measure of the uncertainty of the future drought changes, based on decadal frequencies and the mean event duration. For all cases, the models tend to agree with the projected changes, more for the SPEI rather than the SPI, particularly for the end-of-century and considering the RCP8.5 scenario. The inclusion of temperature through potential evapotranspiration in a drought index is thus paramount in improving the projections.