Wet Spells and Associated Moisture Sources Anomalies across Danube River Basin

: The Danube River Basin is the second longest catchment basin in Europe and exhibits intense climatological diversity. In recent decades, the frequency and intensity of daily precipitation extremes have suffered from an increment in many parts of the world, including Central and Eastern Europe. Wet spells are deﬁned by the number of consecutive rainy days with different thresholds. The identiﬁcation of wet spells and their trends in the rainfall time is very important for many sectors, such as agriculture, ecology, hydrology and water resources. Wet spells can lead to extreme events and cause ﬂoods and other disasters. In this study, we will attempt to characterise global precipitation in the context of wet spells and associated precipitation depth of wet spells in the Danube River Basin area using daily precipitation data, as well as analysing different approaches to identifying wet spells. The ten most intense wet spells were detected, and the most intense, which occurred on 23 September 1996, was studied in depth in terms of precipitation and associated anomalies, the synoptic situation and the anomalous transport of moisture using a Lagrangian approach. The existence of a marked west-east dipole in the ﬁeld of sea level pressure between the Atlantic Ocean and the eastern Mediterranean leads to the anomalous moisture transport from the Northern Atlantic Ocean to the Mediterranean Sea, where a higher available amount of moisture existed, and subsequently penetrated within the low positioned over the Danube River Basin. In addition, an Atmospheric River was also responsible for the wet conditions in the Danube River Basin. The combination of all these factors was responsible for the extreme precipitation linked with the wet spell.


Introduction
Extreme precipitation events in Central Europe, including the area inside of the Danube River Basin (DRB), have become very common in the last few decades and have usually been associated with related phenomena, such as flooding, landslides, storms, significant material damage and human sacrifices. One of the primary causes of extreme weather events may be a consequence of global climate changes [1]. It is known that climate change is one of the major causes of increasing temperatures, precipitation amounts and variability of precipitation events. The Danube River Basin has a very diverse climate and notably variable precipitation characteristics because of the proximity of the Atlantic Ocean, Mediterranean Sea and Alps mountain range [2]. Thus, the significant amount of precipitation in the basin area throughout the whole year may be due to the impact of these factors. When the annual precipitation scale is considered, it is clear that those months with the maximum amount of precipitation usually occur during the summer. This phenomenon is especially highlighted in the low-lying part of the Danube River Basin, where convective precipitation makes To understand the mechanism of increase in mean and extreme precipitation, it is important to analyse wet spell events on different time scales and their associated precipitation magnitude [5]. Europe is similar to global trends, in that increases in average and extreme rainfall and their variability are expected for areas in the middle northern portion, signifying an increasing risk of flooding; alternatively, in southern regions, less precipitation and frequent dry spell periods may lead to increasing drought trends [1].
The primary and most common reason for river flooding in Central Europe, including the area of the Danube River Basin, is heavy precipitation events, with the exception being snow melting during the winter period. It is important to stress that those heavy rain events of different duration (in the range of one to several days) usually result in a large amount of precipitation per square metre. The most significant conditions are certainly intensity, magnitude and spatial-temporal distribution of precipitation, but also, significant roles have specific conditions within the river basin at the moment of heavy precipitation [6]. In recent years, the largest floods in the area of the Danube River Basin occurred in 2002, 2006, 2013 and 2014. As a consequence, all of the associated areas experienced economic losses in the range of billions of euros and, what is much worse, suffered human causalities [7]. Many authors have studied these flood events, their causes and impacts [8,9]. Generally, two main reasons are cited as leading to the extreme flood events. The floods which occurred in 2002 and 2006 were induced by a huge amount of precipitation falling over a short time period, which produced a massive single flood event in the territory of the Danube River Basin. In contrast, the 2010 flood occurred as a consequence of a high number of rainfall events during the whole year, which caused a large number of flood events throughout the area of the Danube, economic losses around two billion euros, and 35 causalities [7].
As previously indicated, when extreme precipitation events are analysed, two terms (or factors) should be taken into consideration: (i) if the event occurs during a short period (hours or less than one day in extreme situations) which happens as a result of a strong convergence of atmospheric water vapour with local dynamic processes, or (ii) if it occurs during a prolonged time period, when the extreme precipitation is related to huge amounts of precipitation occurring over several weeks, To understand the mechanism of increase in mean and extreme precipitation, it is important to analyse wet spell events on different time scales and their associated precipitation magnitude [5]. Europe is similar to global trends, in that increases in average and extreme rainfall and their variability are expected for areas in the middle northern portion, signifying an increasing risk of flooding; alternatively, in southern regions, less precipitation and frequent dry spell periods may lead to increasing drought trends [1].
The primary and most common reason for river flooding in Central Europe, including the area of the Danube River Basin, is heavy precipitation events, with the exception being snow melting during the winter period. It is important to stress that those heavy rain events of different duration (in the range of one to several days) usually result in a large amount of precipitation per square metre. The most significant conditions are certainly intensity, magnitude and spatial-temporal distribution of precipitation, but also, significant roles have specific conditions within the river basin at the moment of heavy precipitation [6]. In recent years, the largest floods in the area of the Danube River Basin occurred in 2002, 2006, 2013 and 2014. As a consequence, all of the associated areas experienced economic losses in the range of billions of euros and, what is much worse, suffered human causalities [7]. Many authors have studied these flood events, their causes and impacts [8,9]. Generally, two main reasons are cited as leading to the extreme flood events. The floods which occurred in 2002 and 2006 were induced by a huge amount of precipitation falling over a short time period, which produced a massive single flood event in the territory of the Danube River Basin. In contrast, the 2010 flood occurred as a consequence of a high number of rainfall events during the whole year, which caused a large number of flood events throughout the area of the Danube, economic losses around two billion euros, and 35 causalities [7].
As previously indicated, when extreme precipitation events are analysed, two terms (or factors) should be taken into consideration: (i) if the event occurs during a short period (hours or less than one day in extreme situations) which happens as a result of a strong convergence of atmospheric water vapour with local dynamic processes, or (ii) if it occurs during a prolonged time period, when the extreme precipitation is related to huge amounts of precipitation occurring over several weeks, months or seasons; in this last case, the duration of rainfall is the primary reason for the accumulation of large amounts of precipitation, which impacts nature and society through flood events [10].
There is no precise and generally accepted method for the calculation and identification of extreme precipitation events. Some authors [11,12] have used the highest recorded precipitation amount at some rain gauge stations as the selection criteria, whereas others have utilised the socioeconomic impacts of extreme precipitation events as their criterion [13]. Typically, all methodologies for ranking and identifying extreme precipitation events at least require daily precipitation values [1]. However, ranking and analysing extreme precipitation events (wet spell events) does not depend solely on available daily precipitation data, but additionally requires dataset resolution to characterise the precipitation [14] and the chosen criteria for the methodology. There have been many studies in recent years that have given attention to extreme precipitation events in the context of wet spells and associated accumulation precipitation in wet spells events. In these studies, the authors have used different methodologies for the analysis of wet spells. One method uses IDF curves (Intensity-Duration-Frequency) for characterisation and the study of wet spells behaviour [15]. In this study, they used the term "duration" to refer to a number of consecutive rainy days in the context of one day or more, not in reference to the actual duration of rainfall events. Another approach has used precipitation anomalies to identify extreme wet or dry spells and defines wet spell events as events with a minimum of three consecutive days with precipitation anomalies more than one standard deviation (std) from daily precipitation [16].
Understanding the atmospheric moisture transport is fundamental for explaining the nature of the precipitation during extreme events [17]. The Lagrangian approach has been broadly and satisfactorily used during the last several years to compute changes in moisture along trajectories and to identify sources of moisture or sinks, all around the globe [18]. Although other approaches (such as box models and isotopes) could be used with a similar purpose, the Lagrangian model supports an important benefit: it is able to compute the track of the moisture in time and permits the identification of the main moisture sources. More information about the comparison between certain methodologies can be found in Gimeno et al. [19].
In this study, for the identification and selection of wet spell events, we used the methodology developed by Ramos et al. [20,21] for the ranking of high-resolution daily precipitation extremes. The long time period analysed was from 1981 to 2015. The main objectives of this work are (i) to rank wet spell events in regards to different time scales in duration from 1 to 10 days for the whole area of the Danube River Basin using a daily precipitation dataset in high resolution and (ii) to analyse moisture source anomalies for the most intense wet spell event using a Lagrangian approach.

Precipitation Dataset
In this work, we used the daily precipitation dataset from the Climate Hazards Group Infra-Red Precipitation with Station data (CHIRPS) project [22]. CHIRPS is a relatively new precipitation database, which has been accessible since the beginning of 2015. This database was developed by the USAID Famine Early Warning System Network with the support of scientists at the University of California Santa Barbara (Santa Barbara, CA, USA). To calculate the wet spell events with major precision, the resolution of the daily precipitation database used is critical, which is the main reason to choose the CHIRPS dataset instead of other available daily databases. It is considered the new environmental record for analysing and monitoring extreme events [22]. The CHIRPS database is in a 0.05 • × 0.05 • degree spatial resolution in latitude and longitude from 6-hourly to 3-monthly aggregates, which makes it a unique daily database.
This database presents a combination of three types of precipitation information: global climatologies, satellite-based measurements and in situ rain gauge data. The database used for calculation wet spells events covers a temporal period from January 1981 until December 2015.
The CHIRPS dataset has been successfully used in recent publications to validate other common datasets in areas that present extreme climate or complex topography. For instance, CHIRPS has been Water 2017, 9, 615 4 of 17 used to quantify the impact of decreasing precipitation trends and increasing temperatures trends in the Greater Horn in Africa [22], and it validated using the data from 21 ground stations in Northeast Brazil [23]. CHIRPS dataset has also been used it to obtain a comprehensive evaluation of eight high spatial resolution precipitation products in an Alpine catchment, the Adige Basin in Italy [24], and to validate a hydrological model to simulate stream flow in a complex topography [25]. Additionally, the same database was successfully used for the analysis of precipitation extremes over Cyprus [26] and Bhutan [27].

Ranking of Wet Spell Extreme Precipitation Events
This ranking of wet spell events in the area of the Danube River Basin is based on the method developed by Ramos et al. [20,21], which used daily-normalised precipitation anomalies for the ranking of multi-day extreme precipitation events according to the accumulated amount of precipitation and the spatial distribution in the Iberian Peninsula.
The ranking is based on the magnitude of an event (R), which is obtained after considering the area affected as well as its intensity in every grid point, and taking into account the daily-normalised departure from climatology. This method [20] for the Iberian Peninsula was partially adapted from another approach [28], that has suggested some criteria to classify each day in terms of extremeness using different meteorological variables. The use of normalised precipitation departures from the seasonal climatology allows us to measure the rarity of an event given by the standardised precipitation anomalies. With this standardisation, we are assured of the different statistical distribution of daily areal precipitation among different areas that are being studied. Therefore, it [20,21] can be applied easily to other regions of the world using different gridded precipitation datasets (model data or observations), and the use of different time scales [21] is directly applicable to assess persistent precipitation episodes over a certain region.
To obtain a final R index to calculate the wet spell events in different durations, several steps should be applied prior to the final ranking (as in [20,21]). In the first step, we calculated daily normalised precipitation anomalies (N) for each grid point as the difference between the precipitation values for the day under analysis and the daily mean climatological value and later divided by the daily standard deviation climatological value. With this we ensure the different statistical distribution of daily areal precipitation among different areas of the Danube basin are taken into account and can easily be compared. The use of standardised precipitation anomalies to evaluate heavy large-scale rainfall events was already used with success by other authors in central Europe [6]. Therefore, we are confident that this methodology [21] reflects wet spell events in the Danube River Basin particularly well.
Furthermore, for this computation, only grid points with precipitation amounts above 1 mm were considered. The reference period that we took into account is the complete period of the CHIRPS precipitation data from 1981 until 2015. The noise in both time series has been smoothed by applying a 7-day running mean to the climatological series. Thus, climatological normalised precipitation anomalies are computed, taking into account each day and each grid point. The final daily index, according to which wet spell events were ranked, was: where A denotes the area in percentage which has precipitation anomalies above two standard deviations, and M is the mean value of these precipitation anomalies over A.
In a second step, we have also computed the accumulated precipitation anomalies for a certain period (NCC): that is, N added during different time periods. NCC represents in our case the wet spell events on different time scales. Finally, in the last step, we performed the ranking of wet spell events according to the final index R, which computed the magnitude of precipitation for each wet spell event in different durations. The main idea is to sum the daily normalised anomalies (see previous paragraph) over different time scales (2 to 10 days) to allow ranking the different anomalous precipitation on multi-day periods. For each time scale's (2 to 10) accumulated precipitation standardised anomalies, the magnitude of the wet spell event is R obtained after multiplying: (1) the area (A), expressed as the percentage that has accumulated precipitation anomalies (computed over different time scales) higher than two standard deviations and (2) the mean value of these accumulated precipitation anomalies (M), considering only grid points with precipitation anomalies of more than two standard deviations.
Wet spell events on the time scales from 1 day until a maximum of 10 days have been calculated, but for the sake of simplicity only 1, 3, 5, 7, and 10 days of length are presented in this study (time periods between 7 and 10 days are the typically synoptic time scale). These different multi-day extreme rankings will show that specific extreme events at shorter time scales (e.g., 2-3 days) may be absent from the top ranks at longer time scales lists (e.g., 10 days) and others that appear at long time scales are not present at shorter time scales.
Furthermore, we should make clear that the ranking of wet spell events illustrated in this paper is not related to economic impacts and/or human impacts, such as causalities, injuries or homelessness. The rankings of wet spell events are constituted by the daily high precipitation amount and the associated spatial extension.

Moisture Sources Anomalies
The second part of this work is to compute precipitation anomalies of moisture sources for the identified wet spell events in different time scales. The moisture source anomalies are computed using the Lagrangian FLEXPART V9.0 model. To track the changes in atmospheric moisture along trajectories the approach uses ERA-Interim reanalysis data from ECMWF (European Centre for Medium-Range Weather Forecast), which are accessible at approximately 80 km (T255 spectral) on 60 vertical levels from the surface up to 0.1 hPa, available at each 6-hour time interval [29]. The model developed by Stohl and James [30,31] consists of dividing the global atmosphere into approximately 2 million particles (air masses), which are transported by the three dimensional wind field along their trajectories. Transport time of the trajectories were limited to 10 days, as this is the mean water vapour lifetime in the atmosphere [32]. Changes in specific humidity (q) and locations of air particles were recorded every 6 h. The increases (e) and decreases (p) in moisture along each trajectory may be expressed by changes in specific humidity (q) by the equation: where m is the mass of the particle. By adding (e − p) values of all air particles residing at each time step over a specific area (in this case, over an area of 1.0 × 1.0 degrees in latitude and longitude), it is possible to obtain the instantaneous values of the (E − P) balance, where (E) denotes evaporation and (P) the precipitation rate per unit area. According to Stohl and James [30,31], this approach has two main disadvantages: (1) it is not possible to make separate calculations of E and P, and (2) the results are highly dependent on the input data quality. In addition, the fluctuations in q along individual trajectories may also occur for numerical reasons (e.g., because of the interpolation of q). However, such numerical noise may be partly mitigated by the large numbers of particles contained in an atmospheric column. A more detailed description of the use of this Lagrangian approach and backward/forward analysis for tracking moisture can be found in many research studies which have used this method for tracking variation and/or identification of moisture sources in many different worldwide regions, such as the Mediterranean region [33], Central America [34], Iberian Peninsula [35], Iceland [36], Greenland [37], China [38], and Niger River basin [39].
To calculate the anomalies for the wet spells detected, the first step is to compute the 35-years climatology (1981-2015) for the moisture sources over the Danube River Basin for those specific days of interest (daily climatological value). Precipitation anomalies of moisture sources were calculated as the difference between the E − P > 0 value (backward analysis) for the day/days of the wet spell events and the daily climatological value.  Table 1 presents the top ten extreme wet spells for the area of interest, the Danube River basin, for five accumulated periods of the lengths of 1, 3, 5, 7 and 10 days. The date that appears for each wet spell indicates the final day of each event. For instance, the wet spell event with a length of 3 days on 24 September 1996 (top 1) signifies that this event includes accumulated precipitation anomalies for that day (24 September 1996) and the two previous days: 23 and 22 September 1996, that is, 3 days in total. The same is true for longer periods. Table 1. The top ten wet spell events in the Danube River Basin (DRB) according to the different length of the events (1, 3, 5, 7, and 10 days). Column denoted with A corresponds to area of the DRB in percentages that had precipitation anomalies above 2 std (standard deviation). Column M corresponds to mean magnitude of precipitation anomalies of A, and the last column denoted with R represents the final index of ranking of wet spell events. The final column (R) is the magnitude of the events. In bold are represented wet spell events in different durations that we took into consideration and all events include the most anomalous one for the 1-day ranking identified event with an initial day on 23 September 1996. Focusing on the results (Table 1), it can be observed that the most significant wet spell event for the length of 1 day occurred on 23 September 1996 (in bold). This event is present in all ranking time scales, although in different positions. For the 3-day ranking, it appears in the second position. The top event is, in this case, 24 September 1996, but it is worth noting that it includes the two previous days (23 and 22 September 1996 in the calculation). Thus, it is ultimately the same event. For the 5-day ranking, the same event occurs for the 1, 2, 3, 6 and 7 positions; and for the 7-day ranking, it appears in the 6, 7, 9 and 10 position. Accordingly, we focused our attention on this extreme event, 23 September 1996, to analyse the anomalies in the accumulated precipitation field at different time scales, as well as the moisture source anomalies.

A(%) = Area M = Mean Magnitude R = index of Ranking
As previously mentioned, the R index corresponds to R = A × M, where A is the percentage of the area with precipitation anomalies higher than 2 std, and M is the mean magnitude for the area A. For instance, the top ranked wet spell event of one day duration (23 September 1996) had a R index of R = 198.65, which corresponds to the 44.87% of area (A) of the Danube River Basin with precipitation anomalies above 2 std and 4.43 mean magnitude (M) of the area marked with A. All magnitude values for the top 10 positions of wet spell events in duration from 1 to 7 days are shown in Table 1.
Moreover, we also need to emphasise that a specific wet spell events at a shorter time scale could not appear in the top ranking of wet spell events with longer lengths due to accumulated amount of precipitation over multi-day extreme precipitation events. That finding means that any of individual precipitation days which are included in the multi-day wet spell events are not equally represented on the highest position of the ranking at the individual daily scale.
From the results in Table 1, we can compare the domain of the affected area that shows anomalous precipitation for the wet spell events at shorter and longer time scales. The top ten events for 1-day length ranking have affected less than 50% of the Danube River Basin area. On this 1-day time scale, the top wet spell event in the ranking does not affect the largest percentage of affected area A = 44.87%; instead, the second ranked event affects the largest percentage of area (A = 50.92%) in comparison with the top ten events in the 1-day time scale. On the other side, for wet spell events calculated for longer time scales, for instance, for the 7-day length, the first ranking event exhibited the biggest affected area (A = 70.41%). Additionally, from Table 1, it is possible to say that the area affected by the wet spells is higher with longer lengths of calculation. Wet spell events of shorter or longer lengths of time affect smaller or bigger percentages of the DRB area respectively. Also, from Table 1, we can notice that the most anomalous event for the 1-day ranking (23 September 1996) does not appear in the top 10 ranking wet spell events on time scale for the duration of 10 days. In this last time scale the event of 23 September 1996 has position 18 and 19. So, we can conclude that the most anomalous wet spell for the 1-day ranking was also anomalous at a 7-days period of ranking, but at the 10-days period of ranking other events appear more extreme.

Precipitation
As previously discussed, this paper focuses from this point forward on the 23 September 1996 wet spell event. Figure 2 shows the accumulated precipitation during this wet spell event for 1-day's duration, as well as 3-day, 5-day and 7-day (Figure 2a-d, respectively). The accumulated precipitation maximum was recorded on the western part of the Danube River Basin area with an amount of precipitation over 100 mm/day, reaching 150 mm/day in the all-time scales and reaches the southern part of the DRB for the wet spell event with duration of 7 days. Intense precipitation with values over 100 mm/day also occurred over the eastern part of the domain. affected by the wet spells is higher with longer lengths of calculation. Wet spell events of shorter or longer lengths of time affect smaller or bigger percentages of the DRB area respectively. Also, from Table 1, we can notice that the most anomalous event for the 1-day ranking (23 September 1996) does not appear in the top 10 ranking wet spell events on time scale for the duration of 10 days. In this last time scale the event of 23 September 1996 has position 18 and 19. So, we can conclude that the most anomalous wet spell for the 1-day ranking was also anomalous at a 7-days period of ranking, but at the 10-days period of ranking other events appear more extreme.

Precipitation
As previously discussed, this paper focuses from this point forward on the 23 September 1996 wet spell event. Figure 2 shows the accumulated precipitation during this wet spell event for 1-day's duration, as well as 3-day, 5-day and 7-day (Figure 2a-d, respectively). The accumulated precipitation maximum was recorded on the western part of the Danube River Basin area with an amount of precipitation over 100 mm/day, reaching 150 mm/day in the all-time scales and reaches the southern part of the DRB for the wet spell event with duration of 7 days. Intense precipitation with values over 100 mm/day also occurred over the eastern part of the domain.  Figure 3 shows the mean daily precipitation anomalies for the Danube River Basin area for the whole period of the CHIRPS precipitation dataset for the 23 September 1996 wet spell event at durations of one, three, five and seven days. The largest positive precipitation anomalies are stressed for the wet spell events with shorter durations, one and three days (Figure 3a,b respectively), where the magnitude of the anomalies are in a range of 30 until more than 60 mm/day over the main area of the Danube River Basin. Wet spell events on the longer time scales (Figure 3c,d) also showed positive values of mean daily precipitation anomalies but with considerably less value, lower than 20 mm/day.  Figure 3 shows the mean daily precipitation anomalies for the Danube River Basin area for the whole period of the CHIRPS precipitation dataset for the 23 September 1996 wet spell event at durations of one, three, five and seven days. The largest positive precipitation anomalies are stressed for the wet spell events with shorter durations, one and three days (Figure 3a,b respectively), where the magnitude of the anomalies are in a range of 30 until more than 60 mm/day over the main area of the Danube River Basin. Wet spell events on the longer time scales (Figure 3c,d) also showed positive values of mean daily precipitation anomalies but with considerably less value, lower than 20 mm/day.

Meteorological Configuration
The synoptic situation for the day 23 September 1996 is shown in Figure 4. The field of daily composite mean sea level pressure (SLP) and total 850 hPa geopotential height, and their anomalies related to the wet spell event analysed are plotted. Those panels to the left (Figure 4a,c) show the climatology for 23 September 1981-2015, and the right panels (Figure 4b,d) show the anomalies. The climatological SLP and geopotential at 850 hPa (Figure 4a,c) denote that over the Atlantic region an anticyclone dominates the general pattern, extending the situation over Europe and the Mediterranean region. However, the anomalies (Figure 4b,d) show an intense low-pressure system over the Italian Peninsula, the Adriatic Sea, Croatia and Slovenia. The central low peaks at 988.3 hPa in SLP (1245.18 hPa in 850 hPa geopotential high level). This confirms that the most anomalous wet spell event for the 1-day ranking was characterised by a cyclone over Southern Europe, crossing the area of the Danube River and leading to heavy precipitation in this region. However, one strong anticyclone is positioned over the Northern Atlantic Ocean.

Meteorological Configuration
The synoptic situation for the day 23 September 1996 is shown in Figure 4. The field of daily composite mean sea level pressure (SLP) and total 850 hPa geopotential height, and their anomalies related to the wet spell event analysed are plotted. Those panels to the left (Figure 4a,c) show the climatology for 23 September 1981-2015, and the right panels (Figure 4b,d) show the anomalies. The climatological SLP and geopotential at 850 hPa (Figure 4a,c) denote that over the Atlantic region an anticyclone dominates the general pattern, extending the situation over Europe and the Mediterranean region. However, the anomalies (Figure 4b,d) show an intense low-pressure system over the Italian Peninsula, the Adriatic Sea, Croatia and Slovenia. The central low peaks at 988.3 hPa in SLP (1245.18 hPa in 850 hPa geopotential high level). This confirms that the most anomalous wet spell event for the 1-day ranking was characterised by a cyclone over Southern Europe, crossing the area of the Danube River and leading to heavy precipitation in this region. However, one strong anticyclone is positioned over the Northern Atlantic Ocean.

Meteorological Configuration
The synoptic situation for the day 23 September 1996 is shown in Figure 4. The field of daily composite mean sea level pressure (SLP) and total 850 hPa geopotential height, and their anomalies related to the wet spell event analysed are plotted. Those panels to the left (Figure 4a,c) show the climatology for 23 September 1981-2015, and the right panels (Figure 4b,d) show the anomalies. The climatological SLP and geopotential at 850 hPa (Figure 4a,c) denote that over the Atlantic region an anticyclone dominates the general pattern, extending the situation over Europe and the Mediterranean region. However, the anomalies (Figure 4b,d) show an intense low-pressure system over the Italian Peninsula, the Adriatic Sea, Croatia and Slovenia. The central low peaks at 988.3 hPa in SLP (1245.18 hPa in 850 hPa geopotential high level). This confirms that the most anomalous wet spell event for the 1-day ranking was characterised by a cyclone over Southern Europe, crossing the area of the Danube River and leading to heavy precipitation in this region. However, one strong anticyclone is positioned over the Northern Atlantic Ocean.  To show the moisture flux and vertical motion, Figure 5 shows plots for the VIMF (Vertically Integrated Moisture Flux) and its divergence for the climatology (Figure 5a) for all those 23 September dates in the period (1981-2015) and the anomaly (Figure 5b) for our studied case (23 September 1996). Between both figures, a larger difference in the distribution of the divergence-converge pattern and its values is notable. The anomalies show (Figure 5a) that over the Danube River Basin area, there is a significant region of convergence (bluish colours) that is concordant with the highest amounts of precipitation (Figures 2 and 3) experienced over the area. Conversely, Figure 5 shows two main areas of anomalous divergence (in red), one over the Mediterranean Sea and another over the Northern Atlantic Ocean, near the NW coast of the Iberian Peninsula. The convergent area over the Danube River Basin and the divergence over the Mediterranean are positioned around the low level pressure shown in Figure 4 (marked in Figure 5  To show the moisture flux and vertical motion, Figure 5 shows plots for the VIMF (Vertically Integrated Moisture Flux) and its divergence for the climatology (Figure 5a) for all those 23 September dates in the period  and the anomaly (Figure 5b) for our studied case (23 September 1996). Between both figures, a larger difference in the distribution of the divergence-converge pattern and its values is notable. The anomalies show (Figure 5a) that over the Danube River Basin area, there is a significant region of convergence (bluish colours) that is concordant with the highest amounts of precipitation (Figures 2 and 3) experienced over the area. Conversely, Figure 5 shows two main areas of anomalous divergence (in red), one over the Mediterranean Sea and another over the Northern Atlantic Ocean, near the NW coast of the Iberian Peninsula. The convergent area over the Danube River Basin and the divergence over the Mediterranean are positioned around the low level pressure shown in Figure 4 (marked in Figure 5 with a black cross). The anomalous anticlockwise circulation is clear in the VIMF plot, showing an NW-SE direction over the area of convergence over the Mediterranean, and an S-N flow over the divergent zone in the Danube area. Over the convergence area in the Atlantic, a NW dominant flux is evident.

Anomalous Moisture Uptake during the 23 September 1996 Wet Spell Event
Once the top-ranked wet spell event was identified for the DRB and the meteorological situation, the next step was to investigate changes in the moisture transport during the lifetime of the extreme event using a Lagrangian approach. Figure 6 shows the mean climatological sources of moisture for the 1-day, 3-day, 5-day and 7-day lengths of the wet spell. The climatology is calculated using a 35year period (1981-2015) obtained through backward trajectories from the Danube River Basin. The backward analysis allows us to track where the particles gain humidity during their trajectories towards the area of the Danube River Basin. The colours with positive values represent areas where evaporation is greater than precipitation (E − P) > 0; thus, these areas are moisture sources for the DRB. On the other hand, areas where precipitation is greater than evaporation are moisture sinks (E -P < 0). In Figure 6, these areas are marked with bluish colours that represent negative values. For the 1-day length, for instance, to compute the moisture climatological field, all of the values for 23 September along the 35 years are taken into account, and for the 3-day length, the three days involved (22)(23)(24) for the 35 years are taken into account. The moisture source patterns ( Figure 6) for the wet spell event at the different time scales exhibit similar behaviour. For all plots of the regions that show the major positive values of (E − P), the moisture sources are the northern-central and

Anomalous Moisture Uptake during the 23 September 1996 Wet Spell Event
Once the top-ranked wet spell event was identified for the DRB and the meteorological situation, the next step was to investigate changes in the moisture transport during the lifetime of the extreme event using a Lagrangian approach. Figure 6 shows the mean climatological sources of moisture for the 1-day, 3-day, 5-day and 7-day lengths of the wet spell. The climatology is calculated using a 35-year period (1981-2015) obtained through backward trajectories from the Danube River Basin. The backward analysis allows us to track where the particles gain humidity during their trajectories towards the area of the Danube River Basin. The colours with positive values represent areas where evaporation is greater than precipitation (E − P) > 0; thus, these areas are moisture sources for the DRB. On the other hand, areas where precipitation is greater than evaporation are moisture sinks (E -P < 0). In Figure 6, these areas are marked with bluish colours that represent negative values. For the 1-day length, for instance, to compute the moisture climatological field, all of the values for 23 September along the 35 years are taken into account, and for the 3-day length, the three days involved (22)(23)(24) for the 35 years are taken into account. The moisture source patterns ( Figure 6) for the wet spell event at the different time scales exhibit similar behaviour. For all plots of the regions that show the major positive values of (E − P), the moisture sources are the northern-central and western Mediterranean Sea, the Danube River Basin itself, and the northern and western Black Sea. This result is in concordance with a recently published paper [2] in which the sources of moisture for the DRB were analysed in depth.
Water 2017, 9, 615 11 of 17 western Mediterranean Sea, the Danube River Basin itself, and the northern and western Black Sea. This result is in concordance with a recently published paper [2] in which the sources of moisture for the DRB were analysed in depth. In general, for studies about extreme events, it is critical to analyse the differences between the event studied and the climatological conditions. Therefore, the anomalies in the moisture sources field for the 23 September 1996 event on the different time scales were investigated. Figure 7 shows that during this particular wet spell event, some areas reinforced their role as sources of moisture. This phenomenon is observed for the western Black Sea and the southern and most western areas of the Mediterranean Sea. The positive anomalies values reached values higher than 2.5 mm/day over these regions. In addition, it is important to note that a region that climatologically acted as a source can now be a sink, as is the case for the Danube River Basin and the climatological source over the northern Mediterranean Sea (Liguria Sea); and that the areas around central Italy exhibit negative anomalies. Other areas appear as effective sources of moisture, as is the case for the band over the Northern Atlantic Ocean that was not a primary climatological moisture source in the period considered (see Figure 6), but for the event analysed the anomaly pattern showed a positive signal. It is highlighted that for the shorter time scales (1 day) the anomalies are more intense than those for longer lengths (7 days). In general, for studies about extreme events, it is critical to analyse the differences between the event studied and the climatological conditions. Therefore, the anomalies in the moisture sources field for the 23 September 1996 event on the different time scales were investigated. Figure 7 shows that during this particular wet spell event, some areas reinforced their role as sources of moisture. This phenomenon is observed for the western Black Sea and the southern and most western areas of the Mediterranean Sea. The positive anomalies values reached values higher than 2.5 mm/day over these regions. In addition, it is important to note that a region that climatologically acted as a source can now be a sink, as is the case for the Danube River Basin and the climatological source over the northern Mediterranean Sea (Liguria Sea); and that the areas around central Italy exhibit negative anomalies. Other areas appear as effective sources of moisture, as is the case for the band over the Northern Atlantic Ocean that was not a primary climatological moisture source in the period considered (see Figure 6), but for the event analysed the anomaly pattern showed a positive signal. It is highlighted that for the shorter time scales (1 day) the anomalies are more intense than those for longer lengths (7 days). The negative values in the field E − P > 0 anomalies over the Danube River Basin (Figure 7 top right hand) are in concordance with the convergence in Figure 5 (right panel) and the highest amounts of precipitation.

Discussion
The configuration with an anticyclone on the left and the intense low level pressure system on the right acts as a belt of transport for the moisture over the Atlantic flowing to the Mediterranean Sea and it is available to penetrate within the storm and activates the processes for intense precipitation over the Danube River Basin (schematic process in Figure 8). Despite the Mediterranean being the main moisture source for the event (about 55% of the total), the anomalous moisture from the North Atlantic Ocean deserves a little more attention. During the days prior to the low over the Danube River Basin a hurricane occurred in the North Atlantic Ocean, the Hortense hurricane (3-16 September 1996). It has been shown that the hurricane events over the Atlantic may directly impact not only western Europe (e.g., [40]), but also the Mediterranean region

Discussion
The configuration with an anticyclone on the left and the intense low level pressure system on the right acts as a belt of transport for the moisture over the Atlantic flowing to the Mediterranean Sea and it is available to penetrate within the storm and activates the processes for intense precipitation over the Danube River Basin (schematic process in Figure 8).

Discussion
The configuration with an anticyclone on the left and the intense low level pressure system on the right acts as a belt of transport for the moisture over the Atlantic flowing to the Mediterranean Sea and it is available to penetrate within the storm and activates the processes for intense precipitation over the Danube River Basin (schematic process in Figure 8). Despite the Mediterranean being the main moisture source for the event (about 55% of the total), the anomalous moisture from the North Atlantic Ocean deserves a little more attention. During the days prior to the low over the Danube River Basin a hurricane occurred in the North Atlantic Ocean, the Hortense hurricane (3-16 September 1996). It has been shown that the hurricane events over the Atlantic may directly impact not only western Europe (e.g., [40]), but also the Mediterranean region Despite the Mediterranean being the main moisture source for the event (about 55% of the total), the anomalous moisture from the North Atlantic Ocean deserves a little more attention. During the days prior to the low over the Danube River Basin a hurricane occurred in the North Atlantic Ocean, the Hortense hurricane (3-16 September 1996). It has been shown that the hurricane events over the Atlantic may directly impact not only western Europe (e.g., [40]), but also the Mediterranean region (e.g., [41]). Other studies [42] and references therein] state the important role of intense transports of moist air from the tropical and subtropical Atlantic in the occurrence of cold season extreme Water 2017, 9,615 13 of 17 precipitation events in the Mediterranean region. These extreme precipitation events may originate from atmospheric processes associated with the formation of hurricanes or intense cyclones over the Atlantic Ocean. The development of these events is characterised by intense convergence of moist air from the tropics [43] that is followed by a fast intrusion of moist air into the Mediterranean region without significant mixing with the surrounding air. A significant number of these events in the Mediterranean region appear to take place during (or immediately after) Atlantic hurricanes or storms [42].
Hortense was a wet hurricane, as the National Hurricane Center reports (http://www.nhc. noaa.gov/data/tcr/AL081996_Hortense.pdf). It started as low-pressure near Africa on 30 August, it moved westward and it became in a tropical storm on 7 September near the Guadeloupe Islands. Hortense became a hurricane on 9 September over Puerto Rico, then it moved northward intensifying, and it became an extratropical low on 15 September. Figure 9 shows the synoptic configuration for SLP, wind and specific humidity at 900hPa using data from ERA-Interim during 14-16 September 1996, the last days during the hurricane situation and the first two days during the transition to an extratropical cyclone. It is clear that the core of the hurricane transported a higher amount of humidity to extratropical latitudes and it was available for its transport during the following 10 days, the period used for computing the E − P anomaly.
Water 2017, 9,615 13 of 17 (e.g., [41]). Other studies [42] and references therein] state the important role of intense transports of moist air from the tropical and subtropical Atlantic in the occurrence of cold season extreme precipitation events in the Mediterranean region. These extreme precipitation events may originate from atmospheric processes associated with the formation of hurricanes or intense cyclones over the Atlantic Ocean. The development of these events is characterised by intense convergence of moist air from the tropics [43] that is followed by a fast intrusion of moist air into the Mediterranean region without significant mixing with the surrounding air. A significant number of these events in the Mediterranean region appear to take place during (or immediately after) Atlantic hurricanes or storms [42].
Hortense was a wet hurricane, as the National Hurricane Center reports (http://www.nhc.noaa.gov/data/tcr/AL081996_Hortense.pdf). It started as low-pressure near Africa on 30 August, it moved westward and it became in a tropical storm on 7 September near the Guadeloupe Islands. Hortense became a hurricane on 9 September over Puerto Rico, then it moved northward intensifying, and it became an extratropical low on 15 September. Figure 9 shows the synoptic configuration for SLP, wind and specific humidity at 900hPa using data from ERA-Interim during 14-16 September 1996, the last days during the hurricane situation and the first two days during the transition to an extratropical cyclone. It is clear that the core of the hurricane transported a higher amount of humidity to extratropical latitudes and it was available for its transport during the following 10 days, the period used for computing the E − P anomaly.  The tropical cyclones also impact downstream wave breaking, and this remote impact was demonstrated in some episodes of intense rainfall over the Mediterranean in autumn [41]. This suggests that the interaction of tropical cyclones with a midlatitude flow over the western North Atlantic may be considered a perturbation to, rather than a source of, downstream wave breaking [41].
The extratropical cyclone derived from Hortense disappeared during 19 September 1996 over the northern Iberian Peninsula (see Supplementary Material Video S1 green box), but the low over the Danube River Basin was due to a second cyclone within the wave train (Supplementary Material Video S1 blue box) that started on 18 September 1996 over the middle North Atlantic, reached the Mediterranean Sea during the next day, and it was reinforced when it was situated on the Danube region.
The extra tropical cyclone that occurs immediately after the Hortense (highlighted with the blue box on Supplementary Material Video S1) struck the Northwest Iberian Peninsula on 20 September 1996. Figure 10a shows that associated with it there was a long corridor of vertically integrated horizontal water vapour transport (IVT) with the characteristics of an Atmospheric River (AR) around 38.5ºN. ARs are relatively long, narrow regions in the atmosphere that transport most of the water vapour outside of the tropics, they contain high amounts of water vapour and they could be associated with extreme events in terms of rainfall and floods where they impact [17,43,44]. In fact, the ERA-Interim ARs database developed by Ramos et al. [45] for the Iberian Peninsula identifies this particular AR. Then, this fast corridor of moisture (but with lower intensity) penetrated the Mediterranean basin during the following day (Figure 10b), and on the 23 September 1996 ( Figure 10d) the anomalous transport of moisture reached the Danube River Basin, associated with the low situated over our area of study. A close look at the global ARs database of Guan and Waliser [46] indicates that this particular pattern of IVT along the Mediterranean and North Africa towards the Danube is also an AR event.
The tropical cyclones also impact downstream wave breaking, and this remote impact was demonstrated in some episodes of intense rainfall over the Mediterranean in autumn [41]. This suggests that the interaction of tropical cyclones with a midlatitude flow over the western North Atlantic may be considered a perturbation to, rather than a source of, downstream wave breaking [41].
The extratropical cyclone derived from Hortense disappeared during 19 September 1996 over the northern Iberian Peninsula (see Supplementary Material Video S1 green box), but the low over the Danube River Basin was due to a second cyclone within the wave train (Supplementary Material Video S1 blue box) that started on 18 September 1996 over the middle North Atlantic, reached the Mediterranean Sea during the next day, and it was reinforced when it was situated on the Danube region.
The extra tropical cyclone that occurs immediately after the Hortense (highlighted with the blue box on Supplementary Material Video S1) struck the Northwest Iberian Peninsula on 20th September 1996. Figure 10a shows that associated with it there was a long corridor of vertically integrated horizontal water vapour transport (IVT) with the characteristics of an Atmospheric River (AR) around 38.5ºN. ARs are relatively long, narrow regions in the atmosphere that transport most of the water vapour outside of the tropics, they contain high amounts of water vapour and they could be associated with extreme events in terms of rainfall and floods where they impact [17,43,44]. In fact, the ERA-Interim ARs database developed by Ramos et al. [45] for the Iberian Peninsula identifies this particular AR. Then, this fast corridor of moisture (but with lower intensity) penetrated the Mediterranean basin during the following day (Figure 10b), and on the 23 September 1996 ( Figure  10d) the anomalous transport of moisture reached the Danube River Basin, associated with the low situated over our area of study. A close look at the global ARs database of Guan and Waliser [46] indicates that this particular pattern of IVT along the Mediterranean and North Africa towards the Danube is also an AR event.

Conclusions
In this work, we calculated the top 10 ranking wet spell events over the Danube River Basin region using the CHIRPS daily precipitation dataset at 0.05 degrees spatial resolution and for the temporal period from 1981 until 2015, considering different time scales from 1 to 7 days using a method developed by Ramos et al. [20,21]. The most intense event detected occurred on 23 September 1996. This top ranked wet spell events is present in the first position of the ranking for the length of 1 day but also between the top 10 events on 3-, 5-and 7-day time scales. During this event, the existence of a low level surface pressure (988.3 hPa in SLP) caused a situation of extreme precipitation over the DRB, reaching values up to 100 mm/day with peaks of 150 mm/day. The magnitude of the anomalies fluctuated from 30 to 60 mm/day over the main area of the Danube River Basin.
The analysis of the moisture source anomalies using a Lagrangian model of particles (FLEXPART) showed that the system was fed with anomalous moisture from the western and south-central Mediterranean Sea, the western Black Sea, and the northern Atlantic Ocean. This pattern is a consequence of several synoptic conditions: the occurrence of hurricane Hortense during the days prior to the event, which lets available moisture over the North Atlantic, an anticyclone positioned over the Atlantic during at least 10 days, and a low-level pressure system that occurred immediately after this hurricane with an associated Atmospheric River that directly struck the Iberian Peninsula on 20 September and then reached the DRB through the Mediterranean Sea on 23 September 1996. The dipole in the SLP and geopotential in low tropospheric level fields due the low-level pressure and the anticyclone made this transport possible, and the occurrence of an AR was responsible for the anomalous moisture availability in the area. This anomalous moisture transport affected the Danube Basin on 23 September 1996 and produced the extreme and anomalous precipitation event.