Characterisation of Extreme Precipitation Events in the Pyrenees: From the Local to the Synoptic Scale

: Mountain systems within the Mediterranean region, e.g., the Pyrenees, are very sensitive to climate change. In the present study, we quantiﬁed the magnitude of extreme precipitation events and the number of days with torrential precipitation (daily precipitation ≥ 100 mm) in all the rain gauges available in the Pyrenees for the 1981–2015 period, analyzing the contribution of the synoptic scale in this type of event. The easternmost (under Mediterranean inﬂuence) and north-westernmost (under Atlantic inﬂuence) areas of the Pyrenees registered the highest number of torrential events. The heaviest events are expected in the eastern part, i.e., 400 mm day − 1 for a return period of 200 years. Northerly advections over the Iberian Peninsula, which present a low zonal index, i.e., implying a stronger meridional component, give rise to torrential events over the western Pyrenees; and easterly advections favour extreme precipitation over the eastern Pyrenees. The air mass travels a long way, from the east coast of North America, bringing heavy rainfall to the western Pyrenees. In the case of the torrential events over the eastern Pyrenees, the trajectory of the air mass causing the events in these areas is very short and originates in the Mediterranean Basin. The North Atlantic Oscillation (NAO) index has no inﬂuence upon the occurrence of torrential events in the Pyrenees, but these events are closely related to certain Mediterranean teleconnections such as the Western Mediterranean Oscillation (WeMO).


Introduction
The climate of the Pyrenees, in southwestern Europe, is especially diverse, due to the orographic complexity of the sector. The Pyrenees are located in the geographical transition between the wet mild domain of the mid-latitudes and the arid area of the subtropical anticyclone belt. Both its marked topography and its location between two climatically differentiated basins, the Mediterranean and Atlantic basins, mean that this mountain area presents a highly variable climate. The precipitation variability in this area has been extensively studied [1][2][3], since it is clearly associated with different natural hazards in the Pyrenees, such as avalanches [4][5][6], heavy snowfall [7][8][9], or flood events [10][11][12], among many others.
These natural hazards, particularly floods, can be triggered by torrential rainfall [13,14]. This highlights the need to study atmospheric patterns involving this type of torrential events. However, most studies tending to characterize atmospheric patterns are based upon a "Circulation-to-Environment" approach; this means that the atmospheric circulation in a particular area is previously characterized and its implications for a given atmospheric/environmental variable at local scale is subsequently monitored [15]. Examples of this approach are widespread in the Pyrenees area. For example, Esteban et al. [16] modelled daily mean precipitation in Andorra for different circulation types; similarly, Lemus-Canovas [17] employed statistical regression techniques to perform a synoptic classification focused on SW Europe aimed at mapping the daily mean precipitation of each weather type over the whole Pyrenees. Other examples of this type of approach can be found in different parts of the Iberian Peninsula [18][19][20][21][22]. However, another approach for characterising atmospheric patterns related to torrential events is known as "Environment-to-Circulation". It involves selecting days of interest (e.g., days with precipitation ≥ 100 mm) and creating a synoptic classification with these torrential days [14] or simply averaging the synoptic variable of interest on the days in which a given threshold is surpassed at the surface in order to finally perform a synoptic composite [23]. This type of approach has been applied to NE Iberia [24], Iberia [25], Southern Europe [14,26] and to other areas of the world, such as Japan [27] and Northern Chile [28], among others. However, it had never been applied to the whole Pyrenees sector, a fact that constitutes a knowledge gap in relation to the synoptic characterization of torrential events in this region.
In this sense, the main objective of this study involves describing in detail the spatial and temporal behavior of torrential rainfall in the Pyrenees, in order to subsequently synoptically characterize such events and to establish the origin of the air masses driving such episodes. The specific objectives are: (i) to compute the annual maximum daily precipitation values from different return periods, and to quantify the number of days with torrential precipitation (daily precipitation ≥ 100 mm) and establish its temporality; (ii) to characterize the large-scale atmospheric patterns causing such torrential events (≥ 100 mm); and (iii) to employ a backward trajectory analysis in order to determine the main trajectories of the air masses giving rise to the occurrence of torrential events.

Data
Daily precipitation was extracted from the database of the CLIMPY Project (Characterisation of evolution of climate and provision of information for adaptation in the Pyrenees), a transboundary project intended to establish a detailed analysis of recent trends in temperature, precipitation, and snow cover in the Pyrenees [29]. Only the weather stations containing at least 90% of the data in the 1981-2015 period (158 in total) were considered on counting the torrential events (Fig. 1). For the computation of return periods, we only used observatories providing a complete record in their data series. All series were subjected to a previous quality control process described in Serrano-Notivoli et al. [30]. The spatial and altitudinal distribution of weather stations is generally homogeneous, except for the central Pyrenees, for which no historical records meet the criterion for data availability.
We used mean sea level pressure (mslp), temperature at 850 hPa (t850) and geopotential height at 500 hPa (z500) to compute the synoptic classification. These variables were supplied by the ERA-5 reanalysis [31] trough at daily scale, encompassing the area 30°N-60°N and 20°W-20°E at a horizontal resolution of 0.25° for the 1981-2015 period. Finally, backward trajectories were computed by means of the Hysplit model [32], fed by the NCEP/NCAR reanalysis database [33] at a horizontal spatial resolution of 2.5°.

Analysis of extreme events and definition of torrential events
To estimate the return level of precipitation given a specific return period (RP), we employed the generalized extreme value (GEV) distribution for extreme events. The cumulative distribution function of the GEV results from combining the Frechet, Gumbel and Weibull families of distributions in one single distribution function: where three parameters ( , and ) represent the shape, location and scale of the distribution function, respectively. The and 1 + (x − )/ must be greater than zero. The specification of will determine the behavior of the tail of the distribution; so that depending on the value of this parameter, any of the following distributions might be obtained: • ξ > 0 giving the heavy-tailed (Frechet) case • ξ = 0 giving the light-tailed (Gumbel) case. For this case: The GEV family of distributions enables the annual daily maximum precipitation to be modelled with the use of the block maxima approach. This procedure consists of grouping the data into blocks of equal size -one for each year-and then fitting the GEV distribution to the set of maxima corresponding to each of the blocks.
Apart from the modeling of extremes by means of GEV analysis, torrential events were defined as days when daily precipitation ≥ 100 mm. A precipitation equal to or greater than 100 mm day -1 was the threshold used to identify heavy rainfall events in Catalonia (NE Iberia) [24,38], associated with an increase in soil erosion processes [39].

Synoptic classification approach and circulation composites
We applied a PCA to a T-mode (temporal) matrix of mslp, t850 and z500, where the variables (columns) were the 220 days displaying torrential precipitation and the cases (rows) were the grid points of ERA-5. Once the PCA was applied to the standardized data matrix, new variables were obtained, the principal components (PCs), which are linear combinations of the original variables. Subsequently, the PCs explaining most of the variance of the original data must be retained, for which purpose we used the Scree Test [40]. Once the PCs were retained, the components were rotated with a Varimax rotation in order to readjust the orthogonal combination of each PC to obtain a greater variance explained by those PCs of lesser rank [41]. From the rotated PCs, we obtained the loading, i.e. the correlation matrix, which indicates the degree of correlation of each day with respect to each PC. In this sense, the assignment of each day to each of the PCs is based on the value of maximum positive correlation and minimum negative correlation. For example, day 1 is assigned to the highest absolute correlation, but subsequently retains the correlation symbol. For this reason, PC1 may be split into two groups, one for the days with the highest positive maximum correlation, and one for the days presenting the lowest negative correlation. This means that, if we retain 5 PCs, up to 10 weather types can be obtained. We used the R package synoptReg [15] to develop the synoptic classification.
To complement the explanation of the synoptic mechanisms that develop such torrential events, we computed the vertically integrated horizontal water vapor transport (IVT) between the 1000 and 300 hPa levels, following the methodology proposed by Lavers et al. [ IVT q u dp q v dp gg Where q is the specific humidity (kg kg -1 ); u and v are the layer-averaged zonal and meridional winds (m s -1 ), respectively; g is the acceleration due to gravity (m s -2 ) ; and p is the pressure difference between two adjacent pressure levels (Pa).
Precisely, heavy precipitation is the product of two significant factors: a large amount of water vapour in the atmosphere and the upward flow of air, which generates the condensation of water vapour and subsequent precipitation. The Clausius-Clapeyron equation [43] provides a strong thermodynamic relationship between temperature and water vapour. Future climate projections indicate a greater increase in extreme precipitation in the future than in the mean annual total precipitation [44]. Furthermore, the meridional and zonal components of the wind at 300 hPa were used to identify the positions of the polar jet in these torrential events.

Hysplit model
The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model enables trajectories of simple air particles to be computed, and their transport, dispersion, and deposition [32] to be simulated. It is one of the most widely used models in the atmospheric sciences to determine the origin of air masses and to define the source-receptor relationship [45,46]. The use of gridded climate data permits the HYSPLIT model to define the trajectory air parcel arriving at a specific location at a specific time on a specific date. We used the NOAA Hysplit model to determine the backward trajectories of episodes equal to or greater than 100 mm at Darnius (easternmost Pyrenees) and Artikutza (westernmost Pyrenees); both these locations possess the highest number of records above this threshold and in opposite geographical locations within the study area.

Maximum daily precipitation for different return periods
The GEV analysis provided the return levels for 6 RP: 5, 10, 20, 50, 100 and 200 years. For a 5-year RP, maximum daily precipitation was generally estimated to be below 100 mm, except in the extreme east, where values of 150 mm day -1 were reached. The eastern Pyrenees is one of the areas receiving most torrential episodes on the Iberian Peninsula [14,47]. For a 10-year RP, 100 mm day -1 was generally reached throughout the region, with over 150 mm per day in the northeastern Pyrenees. The same spatial pattern was repeated for a 20-year RP, with the northeastern region exceeding 200 mm day -1 . Interestingly, values of 150 mm day -1 were recorded in different observatories in the northwest of the study area. For the most extreme RPs (RP 100 and RP 200), 300 mm day -1 was easily exceeded on the eastern slope of the Pyrenees. Even in the extreme northeast, over 400 mm day -1 was recorded for a RP of 200 years. In the rest of the study area, several observatories in the west and south stand out, providing 200-300 mm day -1 for a 200-year RP. Intense NW wind circulation over the Iberian Peninsula can lead to heavy rainfall in non-typical Mediterranean climates such as that of the Western Pyrenees [35]. The areas presenting the lowest extremes were the central part of the northern slopes of the Pyrenees, where different observatories did not reach 150 mm day -1 for the highest RP (200 year).

Characterisation of torrential events: spatial distribution and seasonality
The frequency of torrential events exhibits certain well-defined spatial patterns. Fig.  3a clearly shows the dependence on distance to the sea of the days with amounts equal to or greater than 100 mm in the eastern part of the study area [48]. Moreover, in the extreme west, the maximum was located slightly southward of the axis of the mountain range (Artikutza), whereas in the extreme east the torrential episodes spread over a wider strip, both north and south of the axis of the mountain range, with a maximum in Darnius (Fig.  3a). The latter phenomenon might be caused by the N-S direction of the coast, which is at a right angle to the mountain range's W-E axis [49]. The remaining rain gauges depend to a certain degree upon altitude. The observatories situated at higher elevations also present a higher number of torrential days as a result of convective precipitation in summer. Fig. 3b shows a great variety of days with seasonally distributed torrential precipitation in in the Pyrenees, a fact that reinforces the pluviometric complexity of the study area. In addition, episodes are observed to be prevalent in winter in the westernmost extreme, when NW advections are frequent [17], and in autumn in the extreme east, when easterly flows are more intense [47]. On the Mediterranean coast, all seasons are perceived to present at least one episode of torrential rain. Moreover, the colour purple symbolises the season of the year with most torrential events (autumn). There are only a few cases in the north-eastern area where the torrential regime is led by winter (blue colour), followed by autumn. At the opposite end (west) of the study area, a large number of observatories indicate that the torrential event occurs mainly in winter, with a transition to the autumn months towards the centre and south of the study area. Spring and summer, however, have a lesser impact, only affecting the northernmost central part of the region and presenting no torrential episodes. Furthermore, the springtime precipitation regime has become less predominant over the Iberian Peninsula during the last few decades due to an increase in autumn rainfall [50]. The border thickness of the box indicates the number of seasons with torrential rainfall, and the boxes with the thickest line therefore correspond to weather stations that recorded at least one event in winter as well as in spring, summer and autumn.

Large-scale attribution of torrential events
Based on the methodology described in section 2.3, we obtained 9 weather types (WT), since 5 PCs were retained, accounting for 80 % of the variance of the original data (Fig.  S1). Fig. 4 presents the spatial representation of the mslp, t850 and z500 of the WTs. All the synoptic patterns presented resulted in days with 100 mm or greater in a rain gauge. The high number of synoptic pattern (9) is due to the geographical and pluviometric diversity of the study area, already seen in the previous maps (Fig. 2 & Fig. 3), and the synoptic variety which for this kind of events is already known in areas of the northeast of the Iberian Peninsula [24].

Figure 4.
Mean sea level pressure (black contours, hPa), 500hPa geopotential height anomalies (shaded colours, dam) and 850 hPa temperature anomalies (dashed contours) of the nine most frequent weather types associated with torrential precipitation in the western Mediterranean region. Above each map, the number of days assigned to each synoptic pattern and their percentages with respect to the total are shown in parentheses.
The most frequent synoptic pattern in the occurrence of torrential events is WT 1, which is associated with 19.1 % of the cases, and with an unequivocal seasonality in the cold half of the year (Fig. S2). This WT is characterized by an atmospheric pressure low in the Balearic Sea, which drives moist maritime winds to the easternmost area of the Pyrenees (Fig. 5). In addition, this evident surface low is also accompanied by a low at 500 hPa and a below-average temperature at 850 hPa over the Iberian Peninsula. This WT1 is very well characterized by the negative phase of the WeMO (Fig. 6) due to the combination of low surface pressure on the south of the Iberian Peninsula with high pressure in northern Italy in this type of situation [24]. Moreover, the MO presents clearly negative values when such a situation occurs. Fig. 5 shows that for WT1 there is evidence of moisture advection in the extreme east of the Pyrenees, which increases the potential of these torrential events [51], whilst a clear southern circulation of the jet stream can be observed. Precisely, this WT 1 is at its most effective in the eastern Pyrenees (Fig. 7b), as well as in other areas of the western façade of the Iberian Peninsula [52], causing a large amount of torrential events (Fig. 7a). Due to its winter seasonality and negative temperature anomalies, the WT1 causes a significant fraction of these episodes to occur as snowfall [53]. WT2 (Fig. 4) is characterized by a trough located at the northwest of the Iberian Peninsula and an increase in pressure towards eastern Europe, a situation which is very well characterized both by the WeMO and by the MO. The location of the low favours warm and moist advection from the Mediterranean over the study area, with an increase in IVT intensity (Fig.  5), giving rise to prefrontal precipitation [51]. This pattern can appear throughout the year, but its maximum frequency is in mid-autumn (Fig. S2). Torrential rainfall during this synoptic situation is concentrated on the southern slope, as shown in Fig. 7, as a result of this important moisture advection of subtropical origin. In fact, Insua-Costa et al. [54] recently showed that moisture advected via an atmospheric river from the tropical and subtropical Atlantic played a major role in a well-known case of this type, the November 1982 flood event. In addition, this type of episode has aroused interest in studying the physical mechanisms of high precipitation events related to this synoptic pattern in the eastern Pyrenees [10,11]. WT3 and WT4 exhibit similar surface characteristics, although the former occurs in late summer and during autumn, and the latter mainly during the winter months (Fig.  S2). In both cases north-northwesterly winds over the Pyrenees are involved. However, in WT4 the main action centres driving the northerly advection are more clearly defined: the Azores anticyclone in the west of the Iberian Peninsula and a depression in Italy, favouring positive values of the WeMO index (Fig. 6). This indicates a stronger advection in WT4 than in WT3. The aforementioned phenomenon can also be observed depending on the intensity and direction of the IVT (Fig. 5) and with a wind trajectory at 300 hPa (Fig.  S3), which is perpendicular to the Pyrenees, mainly in WT4. This brings abundant precipitation to most of the northern façade of the region (see Fig. 7). Precisely, WT3 displays less defined surface conditions, a fact that suggests that the torrential situations under this weather type might be not persistent but rather stormy, a situation also defined in Pineda et al. [55]. Figure 5. vertically integrated horizontal water vapor transport (IVT) intensity (shaded colours, kg m -1 s -1 ) and direction (vectors) of the 9 most frequent atmospheric patterns associated with torrential precipitation in the western Mediterranean region. Above each map, the number of days assigned to each synoptic pattern and their percentages with respect to the total are shown in parentheses.
As opposed to WT1, WT5 occurs more frequently in late spring and early autumn (Fig. S2). This pattern reveals a surface depression over the Balearic Sea driving an eastnortheasterly advection in the eastern area of the Pyrenees and a north-northeasterly advection in the westernmost area, favoring two poles of precipitation maxima in both areas (map 5, Fig. 7b). Such a situation is captured very well by the ULMO teleconnection pattern, which represents the MO dipole, but at 500 hPa (Fig. 5). WT6, with a depression to the northwest of the peninsula and an anticyclone in central Europe, favours a warm advection from the south-southeast with a significant moisture content (Fig. 5) that favors abundant precipitation in the eastern coastal area (Fig. 7b) and also in the southernmost mountainous sector of the study area. This is one of the situations that is captured well by the negative phase of the WeMO (Fig. 6), and it concisely explains the occurrence of torrential-type events in the extreme northeast of the Iberian Peninsula [47], and their pluviometric variability in general [49,56]. WT7 presents a high zonal gradient, which favours westerly flows over the Iberian Peninsula with a long ocean surface circulation, which causes significant moist advection over the Atlantic façade; this exhibits a temporal frequency centred in the months of October, November and December (Fig. S2). This situation promotes notable precipitation in the western half of the Pyrenees, whereas in the eastern half, this contribution is truncated because the advection loses all its moisture due to an adiabatic process of leeward compression. For this situation no teleconnection pattern clearly defines the aforementioned atmospheric pattern, a fact also reported in Merino et al. [51]. WT8 exemplifies a widespread depressional situation centred over France, both on the surface and at 500 hPa, driving a northerly advection over the whole Pyrenees range; this advection is reinforced by a southward shift of the jet stream (Fig. S3) and by default, a displacement of the polar vortex from the Arctic region. The MO and ULMO provide the best definition of such a situation: a large depression over the western Mediterranean and an anticyclonic centre over the eastern Mediterranean. This kind of configuration causes widespread precipitation along the northern slope of the Pyrenees (Fig.  7b), and in some rain gauges in this northernmost part of the Pyrenees this is the WT that causes most torrential episodes (Fig. 7a). Finally, WT9 presents a more residual frequency, and is characterised by a weak depression in the southeast of the peninsula that drives southeasterly surface winds transporting significant moisture and distributes the biggest amount of precipitation over the southeastern extreme of the Pyrenees (Fig. 7b). These results reveal hardly any influence of the NAO in any of the situations causing torrential precipitation events in the Pyrenees. Extreme precipitation events in northeastern Spain are not associated with the NAO index values [47]. The main influence of the NAO on rainfall over the Iberian Peninsula is observed in the central and southwestern parts [57,25]. Esteban et al. [58] found that winter rainfall over southern central Pyrenees comes under the influence of the negative phase of the NAO, which is associated with southerly and southwesterly synoptic circulation (WT2 in Fig. 7a). This can also be seen in the highest values of mean daily precipitation (25-50 mm) in WT2 (Fig. 7b), which are mostly concentrated over the central area of the southern slope of the Pyrenees. The amount reached for a 200-yr RP over the central-south slope of the Pyrenees is rather low (about 200 mm) in comparison with that amount over the extreme east (about 400 mm). This fact denotes a climate typically under Atlantic influence (NAO) in the central-south Pyrenees and a highly variable climate under Mediterranean influence (WeMO) in the easternmost Pyrenees. Esteban et al. [58] found no statistically significant correlation between precipitation over the north slope of the Pyrenees and the NAO index, a fact that is likely a result of the orography barrier.

Backward trajectory analysis
There is an unequivocal contrast between the backward trajectories of Artikutza and Darnius for the three altitudinal layers: 500 hPa (Fig. 8a), 850 hPa (Fig. 8b), and surface (Fig. 8c). The first rain gauge is located at the northwestern end of the Pyrenees (Fig. 3a), very close to the Atlantic coast which implies a clear influence of the westerlies and the rapid westerly circulation. The second rain gauge is at the opposite (eastern) end of the Pyrenees, very close to the coast (Fig. 3a); it is leeward of the westerly flows and is clearly characterized by the Mediterranean climate. Most of the trajectories at 500 hPa (Fig. 8a) for Artikutza originate in the North Atlantic Ocean, implying a long maritime trajectory of the air particles, which come to generate intense rainfall. Darnius is clearly different from the Artikutza, because at a height of 500 hPa a long path of air particles is to be expected, but results reveal that most of the trajectories at that height are short, mainly of Mediterranean and subtropical Atlantic origins. In general, for the case of Darnius, the trajectory is short, but quite variable with regard to its origin. The trajectories at 850 hPa ( Fig. 8b) do not reveal significant differences with respect to those at 500 hPa. They are shorter than those at 500 hPa as a result of atmospheric dynamics and surface friction, among others. In the case of Darnius, it is worth noting that almost all the trajectories at 850 hPa are of Mediterranean origin (see WT1, WT5, and WT9 in Fig. 4) or come from North Africa (see WT2 and WT6 in Fig. 4). Finally, at the surface level (Fig. 8c) the shortest trajectories are detected; these are clearly of North Atlantic origin in the case of Artikutza, although different situations exist which are characterized by trajectories of northern and even northeastern origin travelling over the Cantabrian Sea. In the case of Darnius, practically all the trajectories originate between the Balearic Sea, the Gulf of Lion and the Italic Peninsula; very occasionally, some trajectories originate in the North Atlantic. The surface trajectories para Darnius reveal a strong dependence on the local and mesoscale environments in the Mediterranean basin, as opposed to the rapid surface air-circulation over Artikutza coming from afar, such as western Canada. Another remarkable fact involves the greater homogeneity in Artikutza than in Darnius at all levels with regard to the direction of the backward trajectories; these are mostly zonal and present very diverse origins.

Conclusions
Estimation of return periods revealed that in the easternmost area of the Pyrenees there is a high risk of events greater than 150 mm day-1 for return periods of only 5 years. For extreme return periods (200 years), the amounts are extremely high in this eastern area, reaching 400 mm day -1 . It is specifically this eastern area and the extreme northwest area that provide the most records equal to or greater than 100 mm, i.e. torrential rainfall.
Together with a the climatic and geographic complexity of this region, a wide range of synoptic patterns accounts for these torrential events. Thus, in the northern and western areas of the Pyrenees, the main situations causing torrential events are characterized by advection from the north and northeast at surface level and a meridional circulation of the jet stream over the study area, usually accompanied by a significant moisture content. The westerlies can also cause torrential events when a deep depression is located over the north of the IP. This WT has a large air-circulation over the North Atlantic as shown by the analysis of the backward trajectories. In the central-southern sector of the study area, torrential rainfall is dominated by the second most recurrent WT, which is characterised by low pressure in the northwest of the Iberian Peninsula, and which drives very moist southerly winds. On the other hand, the Mediterranean is subjected to the highest number of torrential events with low pressures over the western Mediterranean, which drive easterly and southeasterly winds over the IP. Only one situation presented a significant amount of torrential events; this was weakly defined at surface level due to cold air at high altitude during the summer season. In the easternmost area of the Pyrenees, the air masses that generate torrential rainfall present a short trajectory, coming from subtropical latitudes, in contrast with what happens in the westernmost area, where the path of the air masses is much longer and usually originates on the east coasts of North America and Canada. The heaviest precipitation in the Pyrenees is mainly monitored by the WeMO index: the eastern and northwestern areas come under the influence of the negative and positive phases of the WeMO, respectively (WT1 and WT4 in Fig. 4 and 6). The NAO index has no effect upon the occurrence of torrential events in the Pyrenees, except for isolated areas on the southern façade of this mountain range.
Further research is required on these extreme precipitation events in Mediterranean mountain ranges such as the Pyrenees. This sector could be directly affected by climate change, with the risk of several associated natural hazards becoming more frequent.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Figure S1: Scree test for the first 30 principal components (PCs). The red dots show the explained variance (%) for each principal component, while the black bars show the cumulative explained variance., Figure  S2. Relative monthly frequencies in percentages of each weather type (WT) for the 1981-2015 period. Figure S3: Daily mean 300 hPa wind direction and speed (m s -1 ) of the 9 most frequent weather types associated with torrential precipitation in the western Mediterranean region. Above each map, the number of days assigned to each synoptic pattern and their percentages with respect to the total are shown in parentheses.