Impacts of Atmospheric Rivers in Extreme Precipitation on the European Macaronesian Islands

: The European Macaronesia Archipelagos (Azores, Madeira and Canary Islands) are struck frequently by extreme precipitation events. Here we present a comprehensive assessment on the relationship between atmospheric rivers and extreme precipitation events in these three Atlantic Archipelagos. The relationship between the daily precipitation from the various weather stations located in the different Macaronesia islands and the occurrence of atmospheric rivers (obtained from four different reanalyses datasets) are analysed. It is shown that the atmospheric rivers’ inﬂuence over extreme precipitation (above the 90th percentile) is higher in the Azores islands when compared to Madeira or Canary Islands. In Azores, for the most extreme precipitation days, the presence of atmospheric rivers is particularly signiﬁcant (up to 50%), while for Madeira, the importance of the atmospheric rivers is reduced (between 30% and 40%). For the Canary Islands, the occurrence of atmospheric rivers on extreme precipitation is even lower.


Introduction
Extreme precipitation events (EPEs) occurring in the western Atlantic coast of the Iberian Peninsula (IP) during the winter months have been historically linked to major socio-economic impacts such as flooding, landslides, extensive property damage and human casualties. These events are usually associated to the presence of low pressure systems of Atlantic origin [1][2][3].
The analysis of the contribution of atmospheric rivers (ARs) to EPEs has been restricted to a few areas of the world, with a strong focus on the eastern North Pacific and their associated impacts on the contiguous North American west coast [4][5][6]. ARs are characterized for being shallow (1-2.5 km in height) and narrow (~300-500 km in width) plumes with high water vapour content, stretching over distances of at least 2000 km, and are often associated with the pre-cold-frontal region of extratropical cyclones in the region of the warm conveyor belt [7]. In addition, for the North Pacific, it was shown the role of the East Asian jet stream which exhibits a strong linkage with the total frequency of ARs that make landfall on the western United States [8]. Over Europe, the large amount of water vapour (WV) that is usually transported by these ARs can also lead to EPEs and flooding as described for a few specific extreme events [3,9]. Likewise, ARs also play an important role in a climatological context [10,11], showing, in particular, that there is a strong relationship between ARs and the occurrence of annual maxima precipitation days in Western Europe [12]. This relationship is especially strong along the Western European seaboard, with some areas having eight of their top 10 annual maxima precipitation days related to the occurrence of ARs.
A close look at the work by [13], shows that the frequency of ARs in the North Atlantic basin is one of the highest in the globe (>12% percent of time steps) with a mean duration ranging from 14 h (in the mid to high latitudes of North Atlantic) to 20 h (in the tropical part of the North Atlantic Ocean). Therefore, it is expected that the ARs have a measurable impact in the European Macaronesian Archipelagos EPEs.
The European Macaronesian region consists of the Portuguese Archipelagos of Azores and Madeira and the Spanish Archipelago of Canary ( Figure 1). The three Archipelagos share regional features: a volcanic origin, a contrasting landscape and a gentle climate. The Azores Islands consist of nine islands located in the middle of the North Atlantic Ocean between 31.5 • W and 25 • W and 39.2 • N and 36.8 • N, while the Madeira Islands are located further south-east (near 17 • W and 32.5 • N). The Canary Archipelago, in the south-eastern European Macaronesian region, are located close to the African coast between 18 • W and 13.5 • W and 29.5 • N and 27.5 • N. There are very few studies looking simultaneously into several of these Macaronesian Archipelagos, although it is worth noticing the studies that do [14,15].
Atmosphere 2018, 9, x FOR PEER REVIEW 2 of 15 A close look at the work by [13], shows that the frequency of ARs in the North Atlantic basin is one of the highest in the globe (>12% percent of time steps) with a mean duration ranging from 14 h (in the mid to high latitudes of North Atlantic) to 20 h (in the tropical part of the North Atlantic Ocean). Therefore, it is expected that the ARs have a measurable impact in the European Macaronesian Archipelagos EPEs.
The European Macaronesian region consists of the Portuguese Archipelagos of Azores and Madeira and the Spanish Archipelago of Canary ( Figure 1). The three Archipelagos share regional features: a volcanic origin, a contrasting landscape and a gentle climate. The Azores Islands consist of nine islands located in the middle of the North Atlantic Ocean between 31.5° W and 25° W and 39.2° N and 36.8° N, while the Madeira Islands are located further south-east (near 17° W and 32.5° N). The Canary Archipelago, in the south-eastern European Macaronesian region, are located close to the African coast between 18° W and 13.5° W and 29.5° N and 27.5° N. There are very few studies looking simultaneously into several of these Macaronesian Archipelagos, although it is worth noticing the studies that do [14,15]. Despite the possible impact of the ARs in EPEs in the European Macaronesian Archipelagos, as far as we are aware, there is only one study that systematically analysed the association between EPEs and ARs in the Madeira Islands [16] for a 10 year period derived from satellite imagery. However, several studies have looked in depth into the atmospheric circulation conditions, including ARs, that led to the 20 February 2010 flash flood in the Madeira Islands, with a death toll of near 50 (e.g., [17]).
Therefore, it is of paramount interest to have a long-term systematic study of the influence of ARs in the EPEs in the three European Macaronesian Archipelagos. To do so, we use a precipitation dataset covering several stations located in various islands of each Archipelago. In addition, a global AR database using different reanalyses is used in order to access the days of occurrence of the ARs in each Archipelago. The main objectives of this work are to evaluate the occurrence and impact of ARs in EPEs on the European Macaronesian Archipelagos for a 30 year period and to evaluate its sensitivity to the AR tracking algorithm when using different reanalyses datasets with different spatial resolutions. Despite the possible impact of the ARs in EPEs in the European Macaronesian Archipelagos, as far as we are aware, there is only one study that systematically analysed the association between EPEs and ARs in the Madeira Islands [16] for a 10 year period derived from satellite imagery. However, several studies have looked in depth into the atmospheric circulation conditions, including ARs, that led to the 20 February 2010 flash flood in the Madeira Islands, with a death toll of near 50 (e.g., [17]).
Therefore, it is of paramount interest to have a long-term systematic study of the influence of ARs in the EPEs in the three European Macaronesian Archipelagos. To do so, we use a precipitation dataset covering several stations located in various islands of each Archipelago. In addition, a global AR database using different reanalyses is used in order to access the days of occurrence of the ARs in each Archipelago. The main objectives of this work are to evaluate the occurrence and impact of ARs in EPEs on the European Macaronesian Archipelagos for a 30 year period and to evaluate its sensitivity to the AR tracking algorithm when using different reanalyses datasets with different spatial resolutions.

Precipitation Dataset
The precipitation dataset includes daily accumulated data from 31 climatological weather stations representative of the three European Macaronesian Archipelagos (Azores, Madeira and Canary) covering the 1980 and 2010 period (for the extended winter months, October to March). The extended winter of a certain year n corresponds to October-November-December (OND) of year n − 1 and January-February-March (JFM) of year n. Overall, the chosen stations were based on a combination of tests for data length, completeness, quality and homogeneity, as in [18,19]. Concerning the Madeira and Azores Islands, the datasets were retrieved from the Portuguese Met Office (IPMA), while the Canary dataset was retrieved from the Spanish Met Office (AEMET).
For each weather station, we have calculated the deciles of the available data (disregarding values below 1 mm/day). These deciles are divided into ten equally spaced bins for later comparison with the ARs datasets.

Precipitation Dataset
The precipitation dataset includes daily accumulated data from 31 climatological weather stations representative of the three European Macaronesian Archipelagos (Azores, Madeira and Canary) covering the 1980 and 2010 period (for the extended winter months, October to March). The extended winter of a certain year n corresponds to October-November-December (OND) of year n − 1 and January-February-March (JFM) of year n. Overall, the chosen stations were based on a combination of tests for data length, completeness, quality and homogeneity, as in [18,19]. Concerning the Madeira and Azores Islands, the datasets were retrieved from the Portuguese Met Office (IPMA), while the Canary dataset was retrieved from the Spanish Met Office (AEMET).
For each weather station, we have calculated the deciles of the available data (disregarding values below 1 mm/day). These deciles are divided into ten equally spaced bins for later comparison with the ARs datasets.   From the western group of the Azores, with two islands (Flores and Corvo), the weather stations located at each airport were used (Figure 2b), while from the eastern group of the Azores, including the islands of São Miguel and Santa Maria, six weather stations were analysed (Figure 2c). Despite its important orographic features, none of the weather stations analysed in the different islands were located above 600 m, due to the lack of data with sufficient long-term series.

Madeira Archipelago
For the Madeira Archipelago ( Figure 3, Supplementary Table S2) we have used eight weather stations (seven in Madeira Island and one in Porto Santo). In the case of the Madeira Islands, three weather stations are placed above 600 m-Bico da Cana (1582 m), Poiso (1360 m) and Santo da Serra (668 m), which are located in the centre of Madeira Island and the highest peak at 1862 m. The island has a mean altitude of 1220 m with a complex orography with steep terrain, and therefore the inclusion of stations at higher altitudes is particularly important to better characterise the impact of ARs. From the western group of the Azores, with two islands (Flores and Corvo), the weather stations located at each airport were used (Figure 2b), while from the eastern group of the Azores, including the islands of São Miguel and Santa Maria, six weather stations were analysed (Figure 2c). Despite its important orographic features, none of the weather stations analysed in the different islands were located above 600 m, due to the lack of data with sufficient long-term series.

Madeira Archipelago
For the Madeira Archipelago ( Figure 3, Supplementary Table S2) we have used eight weather stations (seven in Madeira Island and one in Porto Santo). In the case of the Madeira Islands, three weather stations are placed above 600 m-Bico da Cana (1582 m), Poiso (1360 m) and Santo da Serra (668 m), which are located in the centre of Madeira Island and the highest peak at 1862 m. The island has a mean altitude of 1220 m with a complex orography with steep terrain, and therefore the inclusion of stations at higher altitudes is particularly important to better characterise the impact of ARs.

Atmospheric Rivers Dataset
The ARs detection scheme developed at the global scale by [13] is used here. For this method, there is no need for a reference starting point to search for the ARs; instead, the method isolates contiguous regions of the world with enhanced vertically integrated horizontal water vapour transport (IVT) exceeding a certain IVT threshold (>85th percentile or 100 kg m −1 s −1 , whichever is greater). Each of these regions are subsequently analysed to verify if they fulfil the geometry requirements of length >2000 km, length/width ratio >2 and other considerations indicative of AR conditions. For additional information regarding the methodology, please refer to [13]. The databases used here were kindly provided by these authors.
The global ARs database consists of four different sets, each one corresponding to the application of the method of [13] to four different reanalyses: (ERA-Interim with 1.5° × 1.5° spatial resolution; Climate Forecast System Reanalysis (CFSR) with 0.5° × 0.5° spatial resolution; Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) with 0.625° × 0.5° spatial resolution and NCEP/NCAR with 2.5° × 2.5° spatial resolution) between 1980 and 2010 for the winter extended months (October to March) at a six-hourly temporal resolution. The use of multiple reanalyses is advisable, instead of recurring to a single reanalysis, in order to ensure more robust results. In fact, unlike the large Mediterranean Islands, the European Macaronesian Islands are relatively small and therefore changes in the resolution of the reanalyses datasets, among others, could influence the results. It is also important to highlight that despite using two high resolution reanalyses (MERRA2 and CFSR), most of the islands in the different Archipelagos are not resolved by the four reanalyses, limiting the discussion regarding the impacts at higher altitudes [11,20] or even the AR axis orientation [21].

Relationship between Atmospheric Rivers and Precipitation
In order to extract the ARs that influence each island, we have located the nearest grid point of the different reanalyses datasets relative to the centre of each island. Afterwards, for each six-hourly time step between 1980 and 2010, we check if that grid point was classified as an AR, obtaining a list of days for each island that was under the influence of ARs.
Regarding the precipitation days, as mentioned in Section 2.1, we have divided the precipitation distribution for each weather station in deciles (ten equally spaced percentile bins, ignoring values below 1 mm/day). So the first decile bin would include the days that have precipitation between the 0 percentile and 9 percentile, while the second bin would include precipitation days falling between 10th percentile and 19th percentile, until the last decile bin which will include the most extreme

Atmospheric Rivers Dataset
The ARs detection scheme developed at the global scale by [13] is used here. For this method, there is no need for a reference starting point to search for the ARs; instead, the method isolates contiguous regions of the world with enhanced vertically integrated horizontal water vapour transport (IVT) exceeding a certain IVT threshold (>85th percentile or 100 kg m −1 s −1 , whichever is greater). Each of these regions are subsequently analysed to verify if they fulfil the geometry requirements of length >2000 km, length/width ratio >2 and other considerations indicative of AR conditions. For additional information regarding the methodology, please refer to [13]. The databases used here were kindly provided by these authors.
The global ARs database consists of four different sets, each one corresponding to the application of the method of [13] to four different reanalyses: (ERA-Interim with 1.5 • × 1.5 • spatial resolution; Climate Forecast System Reanalysis (CFSR) with 0.5 • × 0.5 • spatial resolution; Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) with 0.625 • × 0.5 • spatial resolution and NCEP/NCAR with 2.5 • × 2.5 • spatial resolution) between 1980 and 2010 for the winter extended months (October to March) at a six-hourly temporal resolution. The use of multiple reanalyses is advisable, instead of recurring to a single reanalysis, in order to ensure more robust results. In fact, unlike the large Mediterranean Islands, the European Macaronesian Islands are relatively small and therefore changes in the resolution of the reanalyses datasets, among others, could influence the results. It is also important to highlight that despite using two high resolution reanalyses (MERRA2 and CFSR), most of the islands in the different Archipelagos are not resolved by the four reanalyses, limiting the discussion regarding the impacts at higher altitudes [11,20] or even the AR axis orientation [21].

Relationship between Atmospheric Rivers and Precipitation
In order to extract the ARs that influence each island, we have located the nearest grid point of the different reanalyses datasets relative to the centre of each island. Afterwards, for each six-hourly time step between 1980 and 2010, we check if that grid point was classified as an AR, obtaining a list of days for each island that was under the influence of ARs.
Regarding the precipitation days, as mentioned in Section 2.1, we have divided the precipitation distribution for each weather station in deciles (ten equally spaced percentile bins, ignoring values below 1 mm/day). So the first decile bin would include the days that have precipitation between the 0 percentile and 9 percentile, while the second bin would include precipitation days falling between Atmosphere 2018, 9, 325 6 of 15 10th percentile and 19th percentile, until the last decile bin which will include the most extreme precipitation days, i.e., days with precipitation ≥90th percentile. For each weather station, and for each decile bin, we cross-check the dates between the precipitation dates within each bin, and if they are associated to an AR day for each of the four reanalyses.
Finally, and due to the large amount of results, we have also computed the mean and the standard deviation of the four reanalyses for the most extreme days in the ≥90th percentile in order to summarise the most important results.

Results and Discussion
This section is organized into three main sub-sections. Section 3.1 provides a brief overview regarding the EPEs in the Macaronesian Archipelagos, while Section 3.2 focuses on the extended winter ARs frequency between 1980 and 2010. Finally, Section 3.3 summarises the relationship between ARs and EPEs, especially for the most extreme precipitation days (≥90th percentile).

Precipitation Extremes
In order to give a sense of the magnitude of the EPEs in each of the weather stations, we have computed the deciles of the distribution, being divided into ten equally spaced bins. An example of the daily precipitation (mm) distribution of the 50th, 60th, 70th, 80th and 90th percentiles (respectively 5th, 6th, 7th, 8th and 9th decile) for three selected locations in each Archipelago (Ponta Delgada in Azores, Funchal in Madeira and Santa Cruz in Canary) are shown in Figure 5. One can see that for the extended winter (October to March), the 50th percentile (5th decile), which corresponds to the median of the distribution, is just short of 5 mm in Ponta Delgada, and slightly above 5 mm in Funchal and Santa Cruz.
Atmosphere 2018, 9, x FOR PEER REVIEW 6 of 15 precipitation days, i.e., days with precipitation ≥90th percentile. For each weather station, and for each decile bin, we cross-check the dates between the precipitation dates within each bin, and if they are associated to an AR day for each of the four reanalyses. Finally, and due to the large amount of results, we have also computed the mean and the standard deviation of the four reanalyses for the most extreme days in the ≥90th percentile in order to summarise the most important results.

Results and Discussion
This section is organized into three main sub-sections. Section 3.1 provides a brief overview regarding the EPEs in the Macaronesian Archipelagos, while Section 3.2 focuses on the extended winter ARs frequency between 1980 and 2010. Finally, Section 3.3 summarises the relationship between ARs and EPEs, especially for the most extreme precipitation days (≥90th percentile).

Precipitation Extremes
In order to give a sense of the magnitude of the EPEs in each of the weather stations, we have computed the deciles of the distribution, being divided into ten equally spaced bins. An example of the daily precipitation (mm) distribution of the 50th, 60th, 70th, 80th and 90th percentiles (respectively 5th, 6th, 7th, 8th and 9th decile) for three selected locations in each Archipelago (Ponta Delgada in Azores, Funchal in Madeira and Santa Cruz in Canary) are shown in Figure 5. One can see that for the extended winter (October to March), the 50th percentile (5th decile), which corresponds to the median of the distribution, is just short of 5 mm in Ponta Delgada, and slightly above 5 mm in Funchal and Santa Cruz. As expected, when higher deciles are evaluated, it is possible to observe a substantial increase in the value of the precipitation. When analysing the 90th percentile (9th decile) there is a difference of more than 10 mm between Ponta Delgada (near 19 mm) and Santa Cruz (near 30 mm), probably due to the different nature of the EPEs occurring in each region [15]. While the Azores is located in Figure 5. Distribution of daily precipitation (mm) for the 50th, 60th, 70th, 80th and 90th percentiles for three selected locations in each Archipelago: Ponta Delgada (Azores, blue), Funchal (Madeira, green) and Santa Cruz (Canary, yellow).
As expected, when higher deciles are evaluated, it is possible to observe a substantial increase in the value of the precipitation. When analysing the 90th percentile (9th decile) there is a difference Atmosphere 2018, 9, 325 7 of 15 of more than 10 mm between Ponta Delgada (near 19 mm) and Santa Cruz (near 30 mm), probably due to the different nature of the EPEs occurring in each region [15]. While the Azores is located in the middle of the North Atlantic Ocean, and therefore its precipitation is associated with the main storm track [22,23], the location of the Canary, on the tropical region, is more prone to convective precipitation even in winter [24].
In order to synthesize the EPE information, we have included in Figures 2-4 the information regarding the 90th percentile (in mm) of each weather station in brackets in its location. Regarding the Azores Islands (Figure 2), the 90th percentile of the precipitation distribution ranges from 16 mm in Graciosa airport (10 m) to 33 mm in Furnas, São Miguel (290 m). Regarding São Miguel and Terceira Island (where we have several stations), it is important to mention the orographic effect due to the volcano cones usually located in the centre of the island. These orographic features lead to generally higher values of the 90th percentile in locations with higher altitudes.
When analysing the Madeira Archipelago (Figure 3), a similar picture emerges with the orographic effect also playing a role on the 90th percentile of the precipitation distribution. While the higher altitude locations have values for this percentile that vary between 46 mm and 53 mm, the locations near the coast reveal lower values between 27 mm in the southern coast and 30 mm in the northern coast, with the lowest value being found in Machico (20 mm).
Once again, the vast majority of the weather stations in the Canary (Figure 4) are located near sea level, with the exception of Izaña and Rodeus (>600 m), located in the Tenerife Island. In this case, the 90th percentile of the precipitation distribution is relatively low in La Palma, Lanzarote and Fuerte Ventura with values below 20 mm. Regarding Tenerife, the direct effect of orography seems not to play a role when analysing the 90th percentile at higher altitude locations of Izaña (14 mm) and Rodeus (24 mm) when compared with Santa Cruz (30 mm), probably due to Izaña and Rodeus being located downstream of the Vulcano Teide, and therefore in the shadow and protected from the warm and humid winds from the west and north-west.

Atmospheric River Frequency
For each island, we have extracted the number of AR time steps relative to each reanalysis, taking into account the closest point of each gridded dataset to the corresponding island (see Section 2.3). In order to evaluate the number of days under the influence of ARs, the number of AR time steps was divided by four (corresponding to six-hourly time steps). For the sake of simplicity, the time evolution of the results from the winter months is summarised in Figure 6 for the islands of São Miguel, Madeira and Tenerife. On average, São Miguel presents a higher frequency of AR days compared to both Madeira and Tenerife.
It is possible to observe strong inter-annual variability on the occurrence of ARs for all three islands, with a maximum/minimum spanning between circa 10 and 40 AR days for the extended winter months (180 days) in São Miguel Island, corresponding approximately to 5.5% to 22% of the winter days. For the eastern Archipelagos (e.g., Madeira and Tenerife), the inter-annual variability is also noticeable, with a clear ARs maximum activity in the 2010 extended winter peaking to near 45 in Madeira (25% of the winter days) and near 40 AR days in Tenerife (40% of the winter days).
Moreover, when analysing the results of Figure 6 and Table S4, it is possible to state that the distinct spatial resolution of the four different reanalyses (ERA-Interim (blue line); CFSR (red line); MERRA2 (green line) and NCEP/NCAR (purple line)) may influence the frequency of the AR frequencies, not only in the example given in Figure 6, but also at the different islands of the Archipelagos (Table  S4). A closer look at Table S4 shows that for the Azores Archipelago, the ARs mean frequency using the NCEP/NCAR reanalysis is in the same range as for the other three reanalyses. However, in the Madeira and Canary Archipelagos, the ARs mean frequency is lower than the results obtained with the other three reanalyses. Therefore, the coarsest grid resolution of the NCEP/NCAR seems to have more impact in the ARs frequency at the southern Archipelagos (Madeira and Canary) than in the northern Archipelago (Azores). The influence of large scale teleconnections and other climate variability factors on the frequency of ARs over the Macaronesian Islands have not been studied in detail. However, over southern (northern) Europe during the negative (positive) phase of the North Atlantic Oscillation (NAO), there is an increase in the number of ARs [9,10]. In contrast, Ramos et al. [11] showed that the NAO pattern is not correlated with the number of ARs over the Iberian Peninsula but, instead, it is significantly associated with the East Atlantic (EA) pattern. Regarding the precipitation in the European Macaronesian Islands, Cropper and Hanna [15] show that the winter NAO correlates negatively with precipitation across the Azores, Madeira and western and central Canary Islands. In addition, Cropper [14] show that not only do the teleconnection patterns play a role in the climate of the European Macaronesian Islands, but also the semi-permanent Azores subtropical high-pressure system, trade winds and also (extra-/sub-) tropical storms. is not correlated with the number of ARs over the Iberian Peninsula but, instead, it is significantly associated with the East Atlantic (EA) pattern. Regarding the precipitation in the European Macaronesian Islands, Cropper and Hanna [15] show that the winter NAO correlates negatively with precipitation across the Azores, Madeira and western and central Canary Islands. In addition, Cropper [14] show that not only do the teleconnection patterns play a role in the climate of the European Macaronesian Islands, but also the semi-permanent Azores subtropical high-pressure system, trade winds and also (extra-/sub-) tropical storms.

Influence of Atmospheric Rivers on Extreme Precipitation
The relationship between ARs and EPEs is analysed for each individual weather station, taking into account the closest grid point of each reanalysis and applying a similar analysis to the previous works developed by some of us for Iberia [11] and South Africa [25]. To do so, we associate the ARs occurrence and the precipitation in the 31 weather stations of the three Archipelagos (see Section 2.3) for the top decile precipitation bins, starting on the 6th decile (50th-59th percentiles) and concluding in the most extreme precipitation days in the last decile (≥90th percentile). Due to the huge amount of results obtained using the many different datasets (31 weather stations times four reanalyses), we made an effort to synthesize without compromising the main results. Therefore, we focus on a representative example for three selected weather stations located in the three different Archipelagos (Figure 7): Ponta Delgada in São Miguel (Azores); Funchal in Madeira Island and Santa Cruz in Tenerife (Canary).

Influence of Atmospheric Rivers on Extreme Precipitation
The relationship between ARs and EPEs is analysed for each individual weather station, taking into account the closest grid point of each reanalysis and applying a similar analysis to the previous works developed by some of us for Iberia [11] and South Africa [25]. To do so, we associate the ARs occurrence and the precipitation in the 31 weather stations of the three Archipelagos (see Section 2.3) for the top decile precipitation bins, starting on the 6th decile (50th-59th percentiles) and concluding in the most extreme precipitation days in the last decile (≥90th percentile). Due to the huge amount of results obtained using the many different datasets (31 weather stations times four reanalyses), we made an effort to synthesize without compromising the main results. Therefore, we focus on a representative example for three selected weather stations located in the three different Archipelagos (Figure 7): Ponta Delgada in São Miguel (Azores); Funchal in Madeira Island and Santa Cruz in Tenerife (Canary).
Atmosphere 2018, 9, x FOR PEER REVIEW 9 of 15 For Ponta Delgada in São Miguel Island (Azores, Figure 7a), the analysis shows that in the 6th decile bin, the percentage of precipitation days associated with ARs occurrence falls between circa 30% (NCEP/NCAR) and 40% (ERA-Interim, CFRS and MERRA2). These values increase steadily for the higher decile (≥90th percentile) bin of precipitation, reaching values ranging from 45% (NCEP/NCAR) to 60% (MERRA2) of the most extreme precipitation days being associated with ARs. A similar analysis performed for the Funchal station (Madeira, Figure 7b) depicts a similar level of impact of ARs on extreme precipitation to the Ponta Delgada in São Miguel, with values up to 60% in the last decile (≥90th percentile). As for Ponta Delgada, it must be stressed that the results for Funchal are partially dependent on the use of different reanalyses, especially in the last decile bin, with values ranging from 45% in the NCEP/NCAR to 60% in the MERRA2.
In the case of Santa Cruz in Tenerife (Canary Archipelago, Figure 7c), the relationship between ARs and precipitation is lower than that found for both Ponta Delgada and Funchal. However, an overall similar increasing dependence of precipitation for higher deciles is also observed. A weak relationship between ARs and precipitation is found for the lower deciles, and a stronger relationship is found for the higher deciles. Values vary between 30% in the lower decile bins and can reach values over 50% in the bins with the highest precipitation days in the last decile (≥90th percentile).
From the analysis of Figure 7 and from the results of the remaining islands (not shown), it is apparent that the NCEP/NCAR reanalysis have, in general, a slightly lower correspondence between EPEs and AR days. This can be related to the coarser resolution of the NCEP/NCAR reanalysis when comparing to the other three reanalyses. On one hand, we have a lower frequency of AR days in the For Ponta Delgada in São Miguel Island (Azores, Figure 7a), the analysis shows that in the 6th decile bin, the percentage of precipitation days associated with ARs occurrence falls between circa 30% (NCEP/NCAR) and 40% (ERA-Interim, CFRS and MERRA2). These values increase steadily for the higher decile (≥90th percentile) bin of precipitation, reaching values ranging from 45% (NCEP/NCAR) to 60% (MERRA2) of the most extreme precipitation days being associated with ARs.
A similar analysis performed for the Funchal station (Madeira, Figure 7b) depicts a similar level of impact of ARs on extreme precipitation to the Ponta Delgada in São Miguel, with values up to 60% in the last decile (≥90th percentile). As for Ponta Delgada, it must be stressed that the results for Funchal are partially dependent on the use of different reanalyses, especially in the last decile bin, with values ranging from 45% in the NCEP/NCAR to 60% in the MERRA2.
In the case of Santa Cruz in Tenerife (Canary Archipelago, Figure 7c), the relationship between ARs and precipitation is lower than that found for both Ponta Delgada and Funchal. However, an overall similar increasing dependence of precipitation for higher deciles is also observed. A weak relationship between ARs and precipitation is found for the lower deciles, and a stronger relationship is found for the higher deciles. Values vary between 30% in the lower decile bins and can reach values over 50% in the bins with the highest precipitation days in the last decile (≥90th percentile).
From the analysis of Figure 7 and from the results of the remaining islands (not shown), it is apparent that the NCEP/NCAR reanalysis have, in general, a slightly lower correspondence between EPEs and AR days. This can be related to the coarser resolution of the NCEP/NCAR reanalysis when comparing to the other three reanalyses. On one hand, we have a lower frequency of AR days in the NCEP/NCAR reanalysis, especially on the Madeira and Canary Archipelagos as discussed in Section 3.2. On the other hand, the finer resolution reanalyses allow for grid points closer to the respective islands. The mean distance of the closest grid point to the island shows that the NCEP/NCAR presents the highest distance between them (~116 km), while the CFSR (~17.4 km) and the MERRA2 (~18.5 km) presents the lowest distance. Therefore, the use of a coarse resolution reanalysis like NCEP/NCAR may influence the results both in terms AR frequency (lower frequencies), as well as a worse representation of the ARs that strike each island, since the closest grid point analysed is located further from the island when compared with higher resolution reanalysis.
In order to summarise the main results, we have computed for each weather station the mean percentage of days (Figures 8-10) and the respective standard deviations (Supplementary Figures  S1-S3) associated with ARs using the all reanalyses for the most extreme precipitation days (≥90th percentile), allowing us to synthesize without compromising the main results.
Atmosphere 2018, 9, x FOR PEER REVIEW 10 of 15 NCEP/NCAR reanalysis, especially on the Madeira and Canary Archipelagos as discussed in Section 3.2. On the other hand, the finer resolution reanalyses allow for grid points closer to the respective islands. The mean distance of the closest grid point to the island shows that the NCEP/NCAR presents the highest distance between them (~116 km), while the CFSR (~17.4 km) and the MERRA2 (~18.5 km) presents the lowest distance. Therefore, the use of a coarse resolution reanalysis like NCEP/NCAR may influence the results both in terms AR frequency (lower frequencies), as well as a worse representation of the ARs that strike each island, since the closest grid point analysed is located further from the island when compared with higher resolution reanalysis. In order to summarise the main results, we have computed for each weather station the mean percentage of days (Figures 8-10) and the respective standard deviations (Supplementary Figures S1-S3) associated with ARs using the all reanalyses for the most extreme precipitation days (≥90th percentile), allowing us to synthesize without compromising the main results.  group (Figure 8a), for Faial, Pico and a couple of locations in Terceira, more than 50% of extreme precipitation days are associated with ARs. However, for all other locations, the impact of ARs remains below 40%. On the western group (Figure 8b), the values attained in Corvo and Flores are in line with the ones found for the central group. Regarding the eastern group (Figure 8c), the importance of the ARs in the EPEs seems to be higher in Santa Maria (>50%) but somewhat lower in the main island of São Miguel, with values always below <50% of the days with ARs. A similar analysis performed on the Madeira Archipelago ( Figure 9) reveals that the impact of the ARs for the extreme precipitation days (above the 90th percentile) is significantly lower when compared with the Azores. In this case, the percentage of extreme days associated with ARs always remains below 50%, with the highest value being found for Funchal (between 40% and 50%). The higher altitude stations of Bico da Cana and Poiso have their values ranging between 30% and 40%, while the northern locations present values below 30% on the days associated with ARs.
Finally, the results for the Canary Islands are shown in Figure 10, where the range of results is similar to those obtained for Madeira. For both islands of La Palma and Lanzarote, results show that the ARs importance reduces to less than 30% of extreme precipitation days, while for Gran Canaria and Fuerte Ventura, this value stays between 30% and 40% on those days. Focusing on Tenerife Island, where three datasets from three stations are available, results show that the higher altitude locations of Rodeus and Izana present the highest values in terms of importance of ARs in EPEs (between 40% to 50%) not only in Tenerife, but also when compared with the other locations in the Canary Archipelago. In this case, it seems that the effect of orography plays a relevant role since Rodeus (617 m) and Santa Cruz (36 m) are only a few km apart.

Conclusions
The main goal of this work was to evaluate the impacts of ARs in EPEs on the European Macaronesian Archipelagos between 1980 and 2010 for the extended winter months (October to March). To achieve this objective, we have used daily precipitation series covering most of the islands of the three Archipelagos in order to access the dates where EPEs occurred, taking into account the decile bins. Moreover, four AR global databases [13] corresponding to the different reanalyses were used: ERA-Interim; CFSR; MERRA2 and NCEP/NCAR between 1980 and 2010 at a six-hourly temporal resolution.
The number of ARs that strike each Archipelago is different, with higher frequencies in the Azores and lower frequencies in the Madeira and Canary Islands. The Azores are located in the middle of the North Atlantic basin, being often affected by the passage of cold fronts associated with the storm track in the North Atlantic [26]. Since the ARs are usually associated with the low-level jet stream ahead of the cold front of an extratropical cyclone [5], they can therefore be highly relevant when explaining most EPEs.
Regarding the Madeira and Canary Archipelagos, the frequency of ARs is relatively lower when compared with the Azores due to the fact that the storm track is weaker than the one found in the Azores sector whatever the storm track climatology used [27]. In addition, given its southern location, Regarding the Azores Archipelago (Figure 8), the association between ARs and EPEs (≥90th percentile) is higher, with values in the range of 30% and >50% of the days. Focusing on the central group (Figure 8a), for Faial, Pico and a couple of locations in Terceira, more than 50% of extreme precipitation days are associated with ARs. However, for all other locations, the impact of ARs remains below 40%. On the western group (Figure 8b), the values attained in Corvo and Flores are in line with the ones found for the central group. Regarding the eastern group (Figure 8c), the importance of the ARs in the EPEs seems to be higher in Santa Maria (>50%) but somewhat lower in the main island of São Miguel, with values always below <50% of the days with ARs.
A similar analysis performed on the Madeira Archipelago ( Figure 9) reveals that the impact of the ARs for the extreme precipitation days (above the 90th percentile) is significantly lower when compared with the Azores. In this case, the percentage of extreme days associated with ARs always remains below 50%, with the highest value being found for Funchal (between 40% and 50%). The higher altitude stations of Bico da Cana and Poiso have their values ranging between 30% and 40%, while the northern locations present values below 30% on the days associated with ARs.
Finally, the results for the Canary Islands are shown in Figure 10, where the range of results is similar to those obtained for Madeira. For both islands of La Palma and Lanzarote, results show that the ARs importance reduces to less than 30% of extreme precipitation days, while for Gran Canaria and Fuerte Ventura, this value stays between 30% and 40% on those days. Focusing on Tenerife Island, where three datasets from three stations are available, results show that the higher altitude locations of Rodeus and Izana present the highest values in terms of importance of ARs in EPEs (between 40% to 50%) not only in Tenerife, but also when compared with the other locations in the Canary Archipelago. In this case, it seems that the effect of orography plays a relevant role since Rodeus (617 m) and Santa Cruz (36 m) are only a few km apart.

Conclusions
The main goal of this work was to evaluate the impacts of ARs in EPEs on the European Macaronesian Archipelagos between 1980 and 2010 for the extended winter months (October to March). To achieve this objective, we have used daily precipitation series covering most of the islands of the three Archipelagos in order to access the dates where EPEs occurred, taking into account the decile bins. Moreover, four AR global databases [13] corresponding to the different reanalyses were used: ERA-Interim; CFSR; MERRA2 and NCEP/NCAR between 1980 and 2010 at a six-hourly temporal resolution.
The number of ARs that strike each Archipelago is different, with higher frequencies in the Azores and lower frequencies in the Madeira and Canary Islands. The Azores are located in the middle of the North Atlantic basin, being often affected by the passage of cold fronts associated with the storm track in the North Atlantic [26]. Since the ARs are usually associated with the low-level jet stream ahead of the cold front of an extratropical cyclone [5], they can therefore be highly relevant when explaining most EPEs.
Regarding the Madeira and Canary Archipelagos, the frequency of ARs is relatively lower when compared with the Azores due to the fact that the storm track is weaker than the one found in the Azores sector whatever the storm track climatology used [27]. In addition, given its southern location, both Archipelagos are also under the influence of other synoptic systems that can be responsible for intense precipitation, including mesoscale convective systems or cut off low systems [28]. Nevertheless, it was shown by [16] that for a 10 year period (September 2002 to November 2012) in the Madeira Islands, the ARs sometimes contribute to the occurrence of EPEs with high socio-economic impact (e.g., [29]).
We then investigated the relationship between ARs and EPEs between 1980 and 2010 on the European Macaronesian Islands. Regarding the Azores, our results are novel and clearly indicate that the presence of ARs plays a major role in explaining the occurrence of EPEs, particularly in the central and western groups.
Regarding Madeira, our results confirm those found by [16], namely that the island is sometimes struck by an AR capable of producing EPEs, especially for the capital city of the island (Funchal). The inter-annual variability of AR frequency also shows that the winter of 2010 was highly unusual for the Madeira, with many more ARs than the usual range of variability. In addition, the 2010 winter was also one of the wettest in Madeira, culminating with the ARs event in February 2010 that was responsible for numerous landslides and flash flood events in the urban area of Funchal, provoking a death toll of near 50 [16,29].
In the case of the Canary Archipelago, the results obtained are in line with those found for Madeira. However, when analysing the different locations on Tenerife Island, the importance of the ARs in explaining the EPEs is clearly stronger in stations located at higher altitudes (Rodeos and Izana) then in Santa Cruz, probably due to the orography effect associated with these stations.
To conclude, we have developed here the first systematic assessment of the influence of the ARs in the EPEs in the European Macaronesian Islands between 1980 and 2010. This study intends to shed light on quantifying to what extent EPEs observed in these Islands are associated with ARs. This is relevant when analysing past and present EPEs with wide socio-economic impacts in the European Macaronesian Islands, such as the case of February 2010 in Madeira. Likewise, it is also important to take into account when considering the foreseen increment of the occurrence of ARs in the North Atlantic basin under most likely representative climate change scenarios [32].
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4433/9/8/325/s1, Figure S1. Standard deviation of the mean percentage of days, using the four reanalyses, associated with atmospheric rivers for the most extreme precipitation days (≥90th percentile) in the Azores Archipelago: (a) central group, (b) western group and (c) eastern group; Figure S2. Standard deviation of the mean percentage of days, using the four reanalyses, associated with atmospheric rivers for the most extreme precipitation days (≥90th percentile) in the Madeira Archipelago; Figure S3. Standard deviation of the mean percentage of days, using the four reanalyses, associated with atmospheric rivers for the most extreme precipitation days (≥90th percentile) in the Canary Archipelago; Table S1. Information regarding the location and altitude of the weather stations in the Azores Archipelago; Table S2. Same as Table S1 for the Madeira Archipelago; Table S3. Same as Table S1 for the Canary Archipelago. Table S4. Atmospheric rivers mean frequency (days) considering the nearest grid point of each island taking into account the four reanalysis: ERA-Interim; Climate Forecast System Reanalysis (CFSR); Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and NCEP/NCAR.