Future droughts in northern Italy: high-resolution projections using EURO-CORDEX and MED-CORDEX ensembles

We analyse the expected characteristics of drought events in northern Italy for baseline (1971–2000), near (2021–2050), and far (2071–2100) future conditions, estimating the drought spatial extent and duration, the percentage of affected area, and the frequency of drought episodes. To this end, daily ensembles of precipitation and temperature records from Global Climate Models (GCMs) and Regional Climate Models (RCMs) pairs, extracted from EURO-CORDEX and MED-CORDEX for the RCP 4.5 and 8.5 scenarios, are collected at spatial resolution of 0.11 degrees. Before the analysis, model outputs are validated on daily weather station time series, and scaling factors for possible use in bias correction are identified. Annual temperature and precipitation anomalies for near and far future conditions are investigated; drought events are identified by the standardized precipitation evapotranspiration index and standardized precipitation index at the 12-, 24-, and 36-month timescales. This study highlights the importance of using multiple drought indicators in the detection of drought events, since the comparison reveals that evapotranspiration anomaly is the main triggering factor. For both scenarios, the results indicate an intensification of droughts in northern Italy for the period 2071–2100, with the Alpine chain being especially affected by an increase of drought severity. A North-to-South spatial gradient of drought duration is also observed.


Introduction
Droughts are natural hazards modulated by climate variability, and they occur on a variety of different spatial and temporal scales (Mishra and Singh, 2010). In contrast to aridity, that is a permanent climatic feature, droughts are generated by a temporary imbalance of water availability. Determining the spatial and temporal distribution of droughts is not immediate, as droughts develop slowly and are characterised by uncertain severity, duration, and frequency (Dubreuil, 1997;Vicente-Serrano et al. 2010). Cascade effects such as intensification of wildfires are also to be expected (Turco et al. 2018). The Mediterranean region is affected by strong drought episodes, owing also to the high interannual variability of precipitation that reduces freshwater stores (Drumond et al. 2017). Drought episodes during the wet season have significant impact on agriculture and economy, requiring massive use of surface and groundwater (Lamy and Dubreuil, 2013).
Prolonged Mediterranean and European droughts, associated to severe impacts, became more frequent since the 1970s. Manzano et al. (2019) found a strong link between the spatial and temporal distribution of droughts and the changes in atmospheric circulation patterns. The remarkable drought event recorded in summer 2003 affected most central and western Europe and it was characterised by persistent anticyclonic conditions, generated by a strong positive phase of the East Atlantic teleconnection pattern (Luterbacher et al. 2004). Another important factor in the insurgence of summer droughts is the presence of a dry soil moisture anomaly at the end of spring, associated for example with low winter precipitation (D'Andrea et al. 2006). An important drought event was recorded in 2011, mainly affecting south-western Europe, leading to a severe soil moisture deficit in northern Italy, in southern France and in the Iberian Peninsula (Spinoni et al. 2015). This event was characterised by a precipitation deficit that started in winter 2011, generated by the negative phase of the Western Mediterranean Oscillation circulation pattern. During summer 2017, most of southern and western Europe suffered from a major drought episode. The whole event was characterised by weakened zonal circulation and intensification of high-pressure systems. In Spain and Italy, crops were severely affected, with a resulting reduction in agricultural production (García-Herrera et al. 2019). In northern Italy, drought episodes have generally increased since 1983, with a longer duration of the lean period of the main Italian river, the Po (Baronetti et al. 2020).
Future climate projections can be obtained using the Representative Concentration Pathways (RCPs) for the global climate model (GCM) simulations developed in the framework of the Coupled Model Intercomparison Project Phase 5 (CMIP5, Taylor et al. 2012), with typical spatial resolution in the range 100-250 km. To represent climate variability on smaller scales, Regional Climate Models (RCMs) are nested into global models to provide a dynamical downscaling of the GCMs, with a spatial resolution in the range 10-50 km (Colmet-Daage et al. 2018). The CORDEX project, in particular, has provided high spatial resolution models for different European regions. The two sub-projects Euro-CORDEX (Jacob et al. 2014) and Med-CORDEX (Ruti et al. 2015) have produced present-day and future simulations at 12 km. A detailed analysis of future projections for the Mediterranean basin is provided by Spinoni et al (2018), who analysed a large number of simulations based on a combination of GCMs and RCMs. Their analysis highlighted a significant climatic change for the Mediterranean basin. The marked negative precipitation trend, linked with a significant temperature increase, could presumably lead to more severe drought episodes. For Italy, Desiato et al. (2015) have collected four RCMs models with spatial resolution of 50 km, belonging to the Med-CORDEX subproject. The authors found a temperature increase of 2.5 °C for RCP 4.5 and of 4.4 °C for RCP 8.5 in the period 2061-2090, compared to the reference period 1971-2010. Spatial distribution maps show that northern Italy, in particular the western sector, is most sensitive to such temperature increase, see e.g. Brussolo et al. (2022). For precipitation, a slight decrease of about 1.5% with respect to the reference period is expected for most of Italy by the end of the century. In the Alpine area, by contrast, weak positive precipitation anomalies are foreseen.
In past years, several studies have been devoted in the investigation of future precipitation and temperature changes in a continental-scale context, such as the whole Mediterranean Basin (Lopez-Bustins et al. 2013) or the whole of Europe (Teuling, 2018s). Future droughts and water scarcity were investigated for the European Alps, using RCM and GCM simulations (e.g. Terzi et al. 2021;Brunner et al. 2019). Owing to the significant precipitation decrease, future drought events for southern Italy were explored (Critto et al. 2016). On the contrary, fewer efforts were devoted to the study of past and future droughts in northern Italy, a borderline Mediterranean climate region. This area is historically rich in water resources and the availability of water for irrigated and rainfed farming is traditionally high. However, drought events have increased in the last twenty years, and such dry and warm conditions determine the prolongation of the Po River lean period. The present work addresses, for the first time, the high-resolution analysis of baseline (1971-2000), near (2021-2050), and far future (2071-2100) drought episodes for northern Italy, providing spatial maps of the estimated drought distribution. The study is performed at weekly scale by analysing GCMs-RCMs pairs belonging to the EURO-CORDEX and MED-COR-DEX sub-projects at 12 km resolution, considering both the historical runs and the future simulations nested into the RCP 4.5 and RCP 8.5 GCM scenarios. First, a regional model validation on ground data is performed, using an innovative application of two softwares (Co.Temp and Co.Rain). Model biasses are determined and the most reliable model configurations for northern Italy are identified and then used in the subsequent analysis. Near future and far future annual anomalies for precipitation and temperature are investigated and the SPEI and SPI drought indices at 12-, 24-, and 36-month timescales are used in the drought analysis. The expected frequency of drought events and their magnitude and the percentage of area affected for specific severity thresholds and the temporal trends are quantified. The spatial extent of baseline and future drought episodes and the possible mechanisms and triggering factors are finally discussed.

Study area
Northern Italy includes the Po Plain, surrounded by two mountain chains, the Alps at the northern and north-western boundaries and the Apennines at the southern boundary (Fig. 1). The northern Adriatic Sea constitutes the eastern boundary. The Po Plain represents the main Italian plain and it covers 46,000 km 2 , 71% of all the Italian plains. It is crossed by the 652km long river Po, the largest Italian river. The Po River flows eastward to the Adriatic Sea and it is fed by the water of all the tributaries streaming from the Apennines and the Alps. The orogenesis of the Alps and Apennines influenced the evolution of northern Italy, resulting in a complex topographical structure. The highest peaks are in the Alps, having mean elevation of 2500 m asl and maximum altitude at the Mont Blanc (4810 m asl). In the Apennines, the highest peak is Monte Cimone, with altitude 2165 m asl. This orographic configuration promotes the formation of a positive atmospheric pressure difference (4-8 hPa) between the upwind and the downwind sides of the Alps, with foehn airflows coming from North, North-West, and West especially during February, March, and April (Fratianni et al. 2009). Northern Italy hydrology has an important contribution from glacier melt and the total runoff from glacierized surfaces is expected to significantly increase in the next 30 years and subsequently to decrease at the end of the twenty-first century (Farinotti et al. 2016).
As indicated in the Koppen classification, the most frequent climate observed in northern Italy is hot summer temperate (Cfa), especially in the Po Plain Vallorani et al. 2018). At the foothills, climate is temperate oceanic (Cfb) and the presence of several lakes has a mitigating influence, allowing the cultivation of Mediterranean crops (wine, olives, citrus fruits). Along the mountain ranges a cold temperate climate (Dw) is present, and the mean temperature of the coldest month is > − 3 °C (Nigrelli et al. 2018). In northern Italy, winters are relatively dry and cold, and the mean temperature recorded in the Po Valley during the reference period 1971-2010 is between 1 and 4 °C. July temperatures range between 22 and 24 °C in the Po Valley (Bigi et al. 2012) and strong summer thunderstorms are common near the Alps. The foothill (Prealpine) zone is the rainiest sector, with annual precipitation of 1500-2000 mm, and in the Po Valley, the annual precipitation ranges between 600 and 800 mm. The highest rainfall amounts are reached in spring and autumn, while the lowest precipitation is observed in winter and summer. Snow episodes are common during the period from early December to early March .

Data set
RCMs simulations exhibit uncertainties that are especially marked in areas with complex topography. For this reason, it is recommended to analyse an ensemble of simulations (Teutschbein and Seibert, 2012). In the following, the RCM outputs of the two sub-projects EURO-CORDEX (https:// www. euro-cordex. net/) and MED-CORDEX (https:// www. medco rdex. eu/) were collected and processed. Here, we considered 18 daily maximum and minimum temperature and 18 daily precipitation ensembles at spatial resolution of 0.11 degree (Table 1)

Model validation and scaling factor estimate
We first considered a model validation procedure to detect which GCMs/RCMs pairs perform best in northern Italy. For this purpose, 10 daily precipitation and temperature series (reference series) were extracted from the GCMs/RCMs historical runs (1971-2000 period) and compared to the corresponding 30-year-long quality controlled and homogenised data series recorded at the ground (candidate series), obtained from the database produced by Baronetti et al. (2020) for northern Italy. The statistical comparison between reference and candidate series was performed with the innovative use of Co.Temp (Guenzi et al. 2019) for temperature and Co.Rain (Guenzi et al. 2017) software algorithms for precipitation. These algorithms were developed to detect discontinuities in precipitation and temperature series due to changes in measurement devices. Such methods compare a temperature (or precipitation) signal coming from an old instrument (taken as candidate) with that recorded by a more recent device (taken as reference), with at least 5 years of overlapping. The statistical approach includes the Root Mean Square Error (RMSE) to show differences or similarities between the series. This allows to estimate whether the two signals recorded (statistically) similar temperature (or precipitation) weak, mean, heavy, and extreme events. After that, the two algorithms estimate the number of candidate and reference events included in these ranges. Using the more reliable model pairs, we estimated the northern Italy bias scaling factor for precipitation and temperature, in case one wants to scale the RCM raw mean values using a linear scaling bias correction method (Shrestha et al. 2017). For this purpose, the precipitation scaling factor P Sf is estimated as: (1) Table1 GCM-RCM model pairs selected from the EURO-CORDEX and MED-CORDEX sub-projects. The models that provided the less reliable results for northern Italy during the validation phase are shaded in grey where P mean_obs is the mean annual value of observed precipitation at the ground and P mean_raw is the mean annual value of raw precipitation estimated from the historical runs of the GCM/RCM model pairs (Ghimire et al. 2019). The temperature scaling factor T Sf is given by: where T mean_obs is the mean annual value of observed temperature at the ground and T mean_raw is the mean annual value of raw temperature estimated from the historical runs of the GCM/RCM model pairs (Fang et al. 2015).
We also tested the significance of annual precipitation and temperature scaling factors for northern Italy. First, the ratio between the standard deviation (σ) of ground data and that of GCM/RCM simulations for the 1971-2000 period was considered. Afterwards, the ratio between the annual scaling factor and the model standard deviation for the baseline was obtained for each model pair. Note, however, that in the following drought analysis we do not use bias corrections as we compare the model (near or far) future with the model baseline.

Methods
To detect the main baseline , near (2021-2050), and far (2071-2100) future weekly drought episodes, the daily precipitation and temperature series of the best-performing models were averaged at weekly temporal scale. The estimate of drought indices requires homogenous periods; thus, it was not possible to use calendar weeks as reference periods because the first day of each year can fall on a different week day. To avoid the occurrence of lap years, each month was divided into 4 weekly periods: the first week goes from the 1 st to the 8 th day; the second from the 9 th to the 15 th day; the third from the 16 th to the 22 nd day; and the fourth week from the 23 rd until the end of the month (Vicente-Serrano et al. 2017).
Before the drought analysis, for each GCM/RCM pair, near and far future precipitation and temperature anomalies were investigated. For both concentration pathways, the mean annual maximum and minimum temperature anomalies were calculated as the difference between the mean annual average in the near (2021-2050) or far (2071-2100) future and the mean annual average of the baseline (1971-2000, Toreti and Desiato, 2008. The same approach was applied to the annual mean precipitation anomalies. The reference evapotranspiration (ET0) was estimated by means of the empirical Hargreaves and Samani (HS) method (Hargreaves and Samani, 1985), using a SPEI package written in R language . Vicente-Serrano et al. (2014) have shown that the HS method usually performed better than other empirical models. The HS method is based on the following expression: where H a and H e are empirical standard parameters; the H a value is 0.0023 and the H e value is 0.5. The solar radiation, R e , is expressed in mm·day −1 . T is the mean temperature ((Tmax + Tmin)/2 in °C) and ΔT is the difference between Tmax and Tmin.
The standardized precipitation index (SPI) and the standardized precipitation evapotranspiration index (SPEI) are the most common drought indicators used in the Mediterranean Basin. The SPI (McKee et al. 1993) is the most widely used index to quantify drought severity, and it is based on the assumption that rainfall variability is the main droughtdriving factor. It is calculated on the total cumulative weekly (or monthly) precipitation data. SPI values usually span the range ± 3, where negative values indicate dry periods and positive values stand for wet periods (Stagge et al., 2015). By contrast, the SPEI (Vicente-Serrano et al. 2010), based upon the same calculation procedure used for the SPI, considers both precipitation and temperature. It assumes that water balance, expressed as the difference between precipitation and evapotranspiration, is the best drought identifier (Tirivarombo et al. 2018). Both indices are able to detect droughts at several timescales. This allows to compare the SPEI and SPI results over space and time.
In the following, the weekly SPEI and SPI values were estimated at 12-, 24-, and 36-month timescales for the selected GCM/RCM pairs. The main drought episodes were identified based on three conditions (Baronetti et al. 2020): • Minimum duration of 3 consecutive weeks.
• Classification of drought episodes in heavy and extreme by means of thresholds. Heavy episodes were defined as those between the index value − 1.65 and − 1.28. Extreme episodes all the drought episodes characterised by an index value < − 1.65. • Drought episodes have to affect at least 25% of the study area (González-Hidalgo et al. 2018).
Baseline , near future (2021-2050), and far future (2071-2100) drought events for RCP 4.5 and 8.5 were obtained identifying their magnitude, temporal duration, and percentage of affected area. The spatial distribution of the expected drought trends was checked using the Mann Kendall test (Mann, 1945), and a statistical significance of 95% was adopted. Finally, the spatial behaviour of drought duration according to SPEI and SPI was compared, based on the comparison between baseline climatic conditions  and future states expected respectively at global mean temperatures of + 2° and + 3 °C above preindustrial levels (Donnelly et al. 2017;Turco et al. 2018).

Model validation and scaling factor estimate
The results of model validation indicated a precipitation overestimation by the models compared to the ground data, as reported in Table 2. Marked difference between model outputs and ground data are observed especially for the extreme event class. The pairs EC-EARTH-RACMO22E and MPI-ESM -CCLM4 provided the worse comparisons, with the highest Root Mean Square Error (RMSE) in the extreme event class (59.56 mm for the former and 56.11 mm for the latter), while the other GMC/RMC pairs gave RMSE that are no larger than 25.59 mm. In Table 2, the GCM/RCM pairs that are more biased in northern Italy are highlighted in grey. Table 3 reports the maximum and minimum temperature validation. In all cases, there was a clear bias of the model outputs to higher temperature values. This phenomenon is noticeable in the extreme cold and extreme warm classes, where the models significantly overestimated the temperature with RMSE values that range between 2.00 and 2.48 °C. The analysis of temperature differences did not identify GMC/ RCM pairs that are more biased than others. Figure 2a shows the spatial distribution of the scaling factors between ground data and GCM/RCM pairs for maximum and minimum temperature and for precipitation (further details on the performance of each model pair are given in Annex 1). For precipitation, values close to 0 indicate overestimation by the climate models, values close to 1 indicate agreement between simulations and ground data, and values larger than one represent underestimation by the GCM/RCM simulations. The precipitation scaling factors reported in Fig. 2a indicate a good agreement between simulations and ground data in the Po Plain and for the lowest parts of the Alpine chain. The standard deviations calculated for the ground data and for the models are also comparable in this area (Fig. 2b), and the scaling factors are no larger than the 0.3 times the GCM/RCMs standard deviation (Fig. 2c). However, for higher elevations in the Alps (Julian and Carnic Alps and western Alps), scaling  Fig. 2a, the spatial distribution of the scaling factors shows that the Po Valley and the foothills have the best agreement between simulated and measured temperatures. For maximum and minimum temperature, the scaling factors are between − 2 and + 2 °C. Figure 2b also shows that the standard deviations calculated on the simulated and ground data are comparable in this area; the scaling factors are up to 5 times the GCM/RCMs standard deviation (Fig. 2c). Model overestimation is observed for the southern part of the study area with scaling factors for maximum temperature of about − 2 °C and -4 °C. In a few areas along the western Ligurian coast, scaling factors between − 4 and − 6 °C were detected. The underestimation is mainly located in the eastern Alps (Julian and Carnic Alps) and western Alps, with scaling factors from about 4 to 6 °C for maximum and minimum temperature.  Figure 3 shows the ensemble average for the near (2021-2050) and far (2071-2100) future mean annual temperature and precipitation anomalies in northern Italy for the RCP 4.5 and 8.5 scenarios (further details on the performance of each model pair are provided in Annexes 2, 3, and 4). An increase of minimum and maximum temperature is evident for the whole study area, with a more complex behaviour of precipitation anomalies. For the near future, both scenarios indicate that the Alpine chain will be significantly affected by global warming, with annual maximum temperature anomalies between 3 and 4 °C and annual minimum temperature anomalies from 2 to 3 °C (Fig. 3a). In the same period, the Po Valley and coastal areas will experience no more than 2 °C of temperature increase. For precipitation, Fig. 3c shows that for RCP 4.5 most of the study area is expected to be interested by very small precipitation changes (anomalies close to 0 mm), while the Ligurian coast and the eastern Alps display positive precipitation anomalies of about 40 mm. A stronger precipitation increase is instead obtained for RCP 8.5 in the western and Julian Alps (anomalies between 80 and 120 mm).

Precipitation and temperature anomalies
For the last 30 years of the twenty-first century, the modelled temperature anomalies display the same spatial behaviour observed for the near future ( Fig. 3b) with annual maximum temperature anomalies between 4 and 6 °C and annual minimum temperature anomalies between 3 and 4 °C in the Alps. In the Po Plain, and especially along the Ligurian coast, the temperature anomalies will not exceed 3 °C. For precipitation, the eastern portion of the Alps (Carnic and Julian Alps) is expected to experience a precipitation increase of about 40 mm, while the western Alps will be affected by a decrease of precipitation (anomalies of − 40 mm). The RCP 8.5 scenario shows an intensification of this tendency, with negative annual precipitation anomalies in the entire Alpine chain, and anomalies of − 80 mm in the western Alps (Fig. 3d). Overall, these results indicate an unchanged, or slightly increasing, precipitation in the Alps for the near future, followed by a significant precipitation decrease at the end of the century. The RCP 4.5 and RCP 8.5 scenarios provide qualitatively analogous results, with the RCP 8.5 case showing enhanced changes. Figure 4 shows the statistics of extreme drought episodes for the baseline (1971-2000), near (2021-2050), and far (2071-2100) future 30-year periods as identified by SPEI and SPI at the 12-month timescale, for RCP 4.5 and RCP 8.5. The plain-coloured histograms represent the SPEI results, while the line pattern marks the SPI outcomes. For each combination of RCM and GMC, the figure reports the total number, the mean percentage of affected area, and the mean duration of the drought events for the considered 30-year period. Overall, Fig. 4 indicates that the SPEI and SPI detect the same number of drought events, with somehow different results in the percentage of affected area and drought duration. The drought analysis based on SPEI indicates that for the baseline period , 3 to 5 extreme drought events are identified. In the period 2021-2050, an increase of extreme drought events is present for the RCP 8.5 scenario, with 5 to 8 extreme drought events (Fig. 4b). For the far future (2071-2100), the number of severe drought episodes is expected to further increase for the majority of the reliable model pairs and for both scenarios.

Drought estimates
The percentage of area affected by droughts is reported in Fig. 4, showing that in the 1971-2000 period, the maximum percentage of area involved by an extreme drought event ranged between 32 and 47%. For the concentration pathway RCP 4.5 in the 2021-2050 period, the maximum percentage of area interested by an extreme drought event is expected to be comparable with that of the baseline period. On the other hand, an increase of the drought-affected area is expected for the last 30 years of the century, with a maximum of 55% of interested area during a single episode as simulated by the pair HadGEM2-RAC-MO22E for the RCP 4.5 scenario. This percentage is even larger in RCP 8.5, with 68% estimated by the same pair.
The maximum duration of the main drought events ranges from 13 to 25 weeks for the 1971-2000 period (Fig. 4). An increase of drought duration is detected for the near future in the RCP 4.5 scenario, with maximum duration ranging from 15 consecutive weeks for the CM5-ALADIN52 pair to 29 weeks for EC-EARTH-HIRHAM5. A further increase is expected for the last 30 years of the century; for the period 2071-2100, the maximum Fig. 4 Number of drought events, mean percentage of area affected by drought and mean number of consecutive drought weeks expected for the baseline , near (2021-2050), and far future . Colors identify different model pairs as indicated in the legend. Plain-coloured histogram represents 12-month SPEI results, and line patterns indicate 12-month SPI results. The less reliable models are underlined and values were not provided duration will be between 21 and 32 weeks, respectively estimated by CM5-ALADIN52 and EC-EARTH-HIRHAM5 in the RCP 4.5 scenario. Similar results, with larger changes, are obtained for RCP 8.5.
The comparison between the two drought indices calculated at 12-, 24-, and 36-month timescale, for the baseline, near, and far future periods is reported in Fig. 5. Positive values stand for SPEI-identified episodes that are more extended than those identified by SPI. In 1971-2000, for all GMC/RCM pairs, there is a tendency of SPEI to detect drought events that are more extended than for SPI; at 36-months SPEI detects events that are 25% and 30% more extended than those estimated by SPI (Fig. 5a). Consistently, in 2021-2050, for all GMC/RCM pairs, SPEI detects more severe drought events than SPI (Fig. 5b). This happens for all the timescales of integration; in particular at 36-months, the events detected by SPEI are 17% and 25% more extended that for SPI. For the 2071-2100 period, the difference between SPEI and SPI is more varied, and it was not possible to detect a consistent difference between the two indices (Fig. 5c).

Drought trends and spatial behaviour
The analysis of the most intense drought events identified an increase of severity in the 2071-2100 period, in terms of both the duration and percentage of drought-affected area. Figure 6 reports the 12-month SPEI trends and their statistical significance in the 2021-2100 period for RCP 4.5 and RCP 8.5, for the four most reliable GCM/RCM pairs. For RCP 4.5 (Fig. 6a), HadGEM2-CCLM4, HadGEM2-RACMO22, and CM5-ALADIN52 show a significant intensification of droughts along the Alpine chain. In this area, the drought index is expected to decrease between − 1.2 and − 1.8 ΔSPEI/30 years. In other parts of northern Italy, the three model pairs indicate no significant and close to zero trend Fig. 6 Spatial distribution and statistical significance of the 12-month SPEI trends for the 2021-2100 period in northern Italy for the following: a RCP 4.5 scenario, b RCP 8.5 scenario. The four most reliable model pairs are shown for the 2021-2100 period, while EC-EARTH-HIRHAM5 indicates a significant increase of droughts also for the central portion (North of the Po Plain, − 1.8 ΔSPEI/30 years). For RCP 8.5, the models estimate an increase of drought severity (Fig. 6b). For HadGEM2-CCLM4 and HadGEM2-RACMO22, the analysis showed significant trends for the whole study area, which becomes especially intense in the Alps (2.4 and -3 ΔSPEI/30 years).
The change in the spatial extent of drought episodes, compared to the baseline condition, was investigated considering also how droughts develop in case of a global warming of + 2 or + 3 °C above preindustrial levels, in keeping with the recent IPCC approach. Figure 7 shows the spatial extent of heavy and extreme drought events according to SPEI and SPI at 12 months. In northern Italy, heavy drought events are homogenously distributed in the 1971-2000 period with duration between 6 and 8 consecutive weeks. Longer extreme drought episodes are detected by SPI in the Alps (8 to 10 weeks, Fig. 7a). Both indicators reveal that, for a global temperature increase of + 2 °C, heavy drought events will be present in the study area, with duration from 5 to 15 consecutive weeks. In a few sectors along the Alps, extreme drought events with a duration of 20 weeks are also foreseen (Fig. 7b). Moving to a + 3 °C global temperature increase, no relevant enhancement of heavy drought events is observed. On the opposite, extreme drought episodes show a strong increase. In particular, the southern portion of the Po Plain is expected to experience extreme droughts with duration of 30 consecutive weeks (Fig. 7c). Once more, the difference in climate change impacts between + 2 and + 3 °C is striking.

Model validation and scaling factors estimate
Model validation was performed for the six historical GCM/RCM pairs in the period 1971-2000. The results provided by the innovative use of Co.Rain and Co.Temp software algorithms evidenced a well-defined model overestimation of the extreme event class. Moreover, the applied methodology indicated that, among the selected simulations, the pairs EC-EARTH-RACMO22E and MPI-ESM-CCLM4 provided the worse comparisons with the observations. The investigation of the precipitation and temperature annual scaling factors for northern Italy showed a good agreement between measured and simulated values in the Po Valley, while in mountain regions and coastal areas, the simulations are affected by considerable uncertainty. Precipitation overestimation and temperature underestimation were detected in proximity of the Carnian and Julian Alps (north-eastern Italy). This is one of the wettest sectors of northern Italy and its climate is strongly influenced by the Scirocco and Libeccio winds, which are two moist and mild flows coming from respectively South-East and South-West .
Underestimation of precipitation and overestimation of temperature by GCM/RCM pairs are observed in proximity of the Ligurian coast. The climate of this area is strongly affected by mesoscale convective patterns generated over the sea, because of the convergence between cold and dry flows coming from the North and warm and wet south-easterly flow. The interaction between these circulations and orography generates heavy precipitation events for prolonged periods (Fiori et al. 2017;Parodi et al. 2017).

Fig. 7
Spatial distribution of drought duration (in consecutive weeks) for "heavy" and "extreme" events (see text for definition), computed from SPI and SPEI: a baseline period ; b global warming of + 2 °C; c global warming of + 3 °C

Precipitation and temperature anomalies
The model results indicate a significant warming in the Alpine chain for the 2021-2050 period. This outcome is consistent with the study of Zimmermann et al. (2013) which found that, for the first 50 years of this century, the Alps will face higher warming trends with respect to the neighbouring sectors. Warming in the mountains is expected to continue in the last 30 years of the century, and this will be more marked in summer (Gobiet et al. 2014). This will be reflected in enhanced glacier volume loss, estimated to be about 80-90% by the end of the century for RCP 4.5 (Radić et al. 2014). Future precipitation anomalies are less clear than temperature changes. The models indicate that most of northern Italy is expected to experience only small precipitation changes, except for some regions (western and Julian Alps) with moister conditions in the near future, see also Brussolo et al. (2022) for a study of climate-induced hydrological changes in northwestern Italy.

Drought estimates
The comparison of future drought characteristics with the reference period  indicated that an increase in the percentage of affected area and drought duration will occur for the last 30 years of the twenty-first century, while the number of drought events remained approximately constant.
The application of SPEI and SPI has revealed a high complexity of the individual drought events: both indices were able to detect the same drought events, but marked differences were present in the estimated percentage of drought-affected area, suggesting that droughts may depend on different triggering factors (Trenberth et al. 2013). During the 1971-2000 period, most GCMs/RCMs pairs indicated that drought episodes appear to be related to an abovenormal evaporative demand. This is presumably linked with the significant positive temperature trends observed by Acquaotta et al. (2015) for 1961-2010 in Piedmont (northern Italy). Analogous results were obtained also for the near future (2021-2050), and Marcos-Garcia et al. (2017) suggested that the increase of temperature will play the most important role in future drought episodes in the first half of this century. On the contrary, for the far future (2071-2100), the relation between the two drought indices is not clear. Several studies (e.g. Merkenschlage andHertig, 2020 andBrogli et al. 2019) found that winter precipitation quantiles are expected to decrease at the end of the century in the whole Mediterranean area. More extended dry periods and extreme precipitation episodes are also expected, associated with drier conditions during summer.

Drought trends and spatial behaviour
This study provided, for the first time, drought spatial distribution maps for the whole of northern Italy. The trend analysis indicated that the entire Alpine Range, characterised by the rainiest and snowiest areas (Diodato et al. 2021 andTerzago et al. 2012), will also experience an increase in drought severity in the RCP 4.5 scenario.
For a + 2 °C increase of global temperature, extreme drought episodes will be characterised by very long dry spells in the Alps. This implies potentially severe impacts on permafrost degradation, influencing surface water quality, erosion and slope instability (Colombo et al. 2019). For a + 3 °C increase of global temperature, a North-to-South gradient of dry spells is identified, and the Po Plain will become more severely affected by extremely long drought events. Long dry periods have a marked impact on agriculture: as a consequence of a large number of consecutive dry days, there can be a significant reduction of crop yield (Raymond et al. 2019). In northern Italy, current drought episodes are mostly influenced by the North Atlantic Oscillation (NAO) and the Mediterranean Oscillation (MO), associated with the North-to-South propagation gradient (Baronetti et al. 2020). This drought gradient, estimated also for the twenty-first century, is probably the result of a stronger positive phase of NAO and MO, leading to drier conditions over northern Italy and most of the Mediterranean basin. Santos et al. (2007) estimated a future extension of the NAO's southern centre of action through the eastern part of the Mediterranean Basin, implying an intensification of NAO effects in southern Europe. This difference between the positive and negative phase of NAO can rationalize the expected, future North-to-South drought gradient and the decrease of precipitation for the last decades of the century.

Conclusions
This work analysed the expected properties of drought episodes in the near (2021-2050) and far (2071-2100) future compared to the baseline conditions  for northern Italy, a crucial region from the point of view of water resources and European economy. EURO-CORDEX and MED-CORDEX GCMs/RCMs pairs at spatial resolution of 0.11 degrees were analysed for the RCP 4.5 and RCP 8.5 concentration pathways. To increase the reliability of the results, model validation was addressed applying, for the first time, the Co.Rain and Co.temp software algorithms, revealing that all model pairs tended to overestimate precipitation and temperature extreme events, and the pairs EC-EARTH-HIR-HAM5, HadGEM2-CCLM4, HadGEM2-RACMO22, and CM5-ALADIN52 performed best for northern Italy. Such validation procedure provided estimates of the precipitation and temperature scaling factors which could be used in linear bias correction techniques, albeit they were not directly employed here as we compared model projections with model baseline (i.e., anomalies). The results indicated that the Po Valley and the foothills are the areas where the GCM/RCM pairs performed best, while in complex environments such as coastal areas and mountain regions, the simulations were affected by considerable uncertainty.
The results of the analysis indicated that, although northern Italy is historically rich in water resources, the trend towards dry conditions already observed in 1971-2000 will continue in the twenty-first century. Most GCM/RCM pairs indicated an increase of drought severity, in terms of duration and percentage of drought-affected area, especially for RCP 8.5 and for the later part of the century. In particular, the Alpine area (a water tower for the surrounding area) will be significantly affected by higher positive temperature anomalies and increasing drought conditions. This study also indicated the importance of adopting multiple drought indices, as the comparison revealed that only one index does not tell the full story. A positive evapotranspiration anomaly appeared to be the main drought-triggering factor for the baseline situation and for the first part of the century. Finally, spatial distribution maps of drought events were provided for northern Italy, showing that the Alpine chain is the most sensitive area in northern Italy.