Long-Term Variability of Wind Speed and Direction in the Mediterranean Basin

: In this work, the ERA5 reanalysis dataset, with its ﬁne spatial and temporal resolution, is used to study the wind speed and direction characteristics in the Mediterranean basin from 1979 to 2020. Their variability, trend and mean values, as well as their association, are examined in the monthly, annual and interannual/decadal time scales. The long-term variability of the wind direction is assessed using the angular variance, while the mean annual and interannual variability are used for the wind speed. The most characteristic and constant ﬂow systems appear in the Gulf of Lion (Mistral) and the Aegean Sea during the summer (Etesian winds). The mean annual variability maximizes in the northern part of the basin, while the western part appears to have the most variable wind directions, while the fastest increasing wind speeds appear in the southwestern Levantine sub-basin. The long-term linear trend for the mean, the 95th and 99th quantiles of annual wind speed and the number of occurrences of extreme wind events have been also assessed. For particular areas, there is an increasing tendency in both intensity and frequency of occurrence of extreme wind events. The linear-circular correlation coefﬁcient has been implemented in order to quantify the linear association between the wind direction and the wind speed. This parameter demonstrates higher annual values for the Mistral wind in the Gulf of Lion and the Etesian winds in the Aegean Sea during the summer. Finally, the comparison of the results of ERA5 with the results of ERA-Interim highlighted signiﬁcant differences in the Mediterranean wind speed and direction characteristics.

Wind speed and direction change over different time scales, as is evident in longterm time series, as well as atmospheric circulation patterns [34]. In order to assess the long-term behavior, trends and variability, the long-term time series of wind speed and direction are necessary [35][36][37]. This necessity is obvious in works such as the study Wind 2022, 2 514 conducted by Pogumirskis [5], where an analysis of the wind direction is presented for the Baltic Sea, and it is consequently emphasized that wind direction variability is also important as it affects human activities, transportation and temperatures in coastal zones. Another important facet of the wind direction variability emerging from this study is its importance in retrieving offshore wind speeds from satellite data [38]. In addition, the assessment of wind direction variability during night-time low wind speeds demonstrated how low wind speed values do not hinder wind direction variability in any way, while large variability has been measured even at low wind speeds. The terrain over which the wind blows seems to play a significant role, as intricate patterns increase the variability for small time scales, as well as the stratification conditions in the atmosphere [39][40][41]. Moreover, the variability in wind direction also affects the fetch length and therefore, the wave propagation schemes and sea-state intensity [42], as well as the estimation of wind and/or wave design values [43,44]. The probabilistic modeling of wind direction seems to be increasingly important for wind farm design and micro-siting, both onshore and offshore [45][46][47], for structural stability and reliability [48][49][50] and long-term economic gains from such applications [51].
However, despite the importance of both wind speed and direction, researchers focus mostly on statistical assessment and the modeling of wind speed, while wind direction is examined implicitly by using wind roses [52,53]. Some exceptions to this are [54]: where wind direction variability has been assessed using angular variance [35,45], where circular statistics have been applied in order to study directional metocean variables, including wind direction and [55] where the variability of wind direction has been examined with respect to the mean and extreme wind events. In the same work [55], it is highlighted that wind direction variability per se is important in the analysis of extreme winds since the "deviations in the direction of the mean wind are not a good indicator for deviations in the direction of extreme winds".
Since wind variability is important for a variety of applications, several studies have focused on this subject. The evaluation of wind speed and direction from 216 stations in the United States from 1961 to 1990 indicated that monthly mean values are similar, while variance values exhibit a seasonal change [56]. McVicar [57] has identified a declining trend in wind speed of −0.014 m/s per year to a global extent. A different study resulted in a −0.011 m/s per year for the Northern Hemisphere [58], and Watson [59] concludes that there appears to be a "very slight" decline in the UK. Conversely, there has also been a study examining long-term time series for the whole of Europe that indicated no particular trend whatsoever, as the authors concluded that the measured variability falls in the area of natural variability, and that uncertainty of the results poses yet another factor for not concluding on the subject of long-term trends. However, the NE Atlantic exhibited the highest average wind speeds and daily variability [36]. A recent analysis regarding the assessment of wind variability, trend and extremes over Europe and the North Atlantic based on ERA5 data has been performed by Laurila [34], by introducing an extreme wind factor (EWF) for the assessment of extremeness of the wind speed. It has been concluded that EWF takes its higher values (on a monthly basis) in southern Europe. Finally, Yu [37] performed an analysis of ERA-Interim wind datasets for Antarctica and the Southern Ocean, focusing on the potential implications of wind energy.
Focusing on the Mediterranean Sea, it is considered to be a climate change hot-spot [60]. In addition, it is an area where offshore wind energy is currently in its development phase. In this respect, wind climate analysis and variability, estimation of trends, frequency and intensity of extreme weather events become important since it shapes to a large extent the economics of any offshore wind energy project during its lifetime. Therefore, in the last few years, what has mostly been the subject of interest for researchers is the study of wind speed variability in the span of 20-30 years, which is the normal lifetime of a wind farm. The existing wind climate and offshore wind power potential analyses are based on gridded datasets obtained from numerical models or remote sensing techniques that cover different periods with various horizontal and temporal resolutions. For example, wind data provided by the SeaWinds scatterometer for the period 2000-2004 with a resolution of 0.25 • over the entire Mediterranean and Black Seas were analyzed to produce the spatial distribution, the seasonal wind stress and the annual and seasonal vorticity [61]. In order to assess ocean surface thermal anomalies, the analysis of extreme wind fields in the western Mediterranean has been performed using a merged data product comprising remotely sensed observations and the ERA-Interim wind dataset for the period 1987-2015 with a spatial resolution of 0.25 • and a temporal resolution of 6 h [62]. Soukissian [35] analyzed wind speed variability and long-term offshore wind characteristics based on the ERA-Interim product for the period 1979-2014 with a spatial resolution of 79 km and a temporal resolution of 6 h. The surface wind climate features in the Mediterranean have been described using the ERA40 dataset for the period 1957 to 2002 with horizontal and temporal resolutions of 2.5 • and 6 h, respectively [54], while the interannual and multi-decadal variabilities in the eastern Mediterranean have been studied based on visual observations from ships [63]. Moreover, on the subject of the Mediterranean cyclones, an extended assessment of their dynamics in the Mediterranean basin using a variety of different datasets has been previously presented [64]; an important finding on cyclogenesis seasonality derived from a study by Flocas [65] where it is proven that after May the cyclone track frequency decreases to increase again after October. The same study indicates that the cyclone track frequency seems to decrease annually.
Long-term wind climate analyses have been also performed in the context of offshore wind energy assessment for the Mediterranean basin. For example, the effects of climate change on future offshore wind energy potential have been studied, based on hindcast wind data for the period 1979-2020 with a spatial resolution of 10 km and a temporal resolution of 1 h [66], and for the period 1961-1990 with a horizontal resolution of 25 km [16]. The offshore wind power potential in the Mediterranean Sea has been estimated using the satellite Blended SeaWinds product of NOAA (1995-2014), with 6 h and 0.25 • temporal and spatial resolutions, respectively [67]. The classification of the offshore wind power potential has been performed using 40 years of 1-hour wind data from MEERA-2 reanalysis [68], while the offshore wind power potential has been estimated using wind data from ERA5 reanalysis and a combination of remotely sensed data [69] and results from ERA-Interim and NCEP [70].
This paper aims to investigate the variability of wind speed and direction in the monthly, annual and interannual/decadal time scales, using the most recent and fine spatiotemporal resolution of the ERA5 reanalysis dataset, extending from 1979 up to 2020. ERA5 is one of the few existing reanalysis products with a 1 h temporal resolution and a 0.25 • × 0.25 • spatial resolution, which combined with the very sufficient length of the available time series provides a robust basis for long-term analysis of the annual, seasonal and decadal timeframes. The results provided in this work will be also compared with corresponding results obtained from the ERA-Interim dataset [35]. The latter dataset covered a slightly shorter period (1979-2014) with a coarser temporal resolution (6 h) and spatial resolution (79 km). See also Section 2.
The area of interest is the entire Mediterranean basin defined by a rectangle with 5 • W-36.5 • E and 46 • N-30 • S coordinates. Apart from the long-term variability of wind speed and its extremes, wind direction will also be studied due to its importance, as already mentioned above, along with their association. In this context, the following issues will be examined in detail: (1) the angular distance at the decadal time scale; (2) the slopes of the site-specific extreme wind speed events; (3) the linear-circular association of wind speed and direction. Moreover, the comparison in the long-term time scale of the results obtained by using the ERA-Interim and ERA5 datasets reveals important differences between the two datasets as regards variability characteristics of wind speed and direction and long-term trends and extremes of wind speed. The structure of this work is the following: Section 2 describes the ERA5 reanalysis dataset, while the theoretical background and mathematical tools used are introduced in Section 3. Section 4 presents the numerical results derived in the annual, monthly and interannual time scales along with some comparisons with the corresponding results obtained by the ERA-Interim dataset. Special emphasis is given to the variability characteristics of wind speed and direction, as well as the trends for the mean and extreme values of wind speed. Finally, in Section 5, a general discussion on the results is provided along with some concluding remarks.

The ERA5 Reanalysis Dataset
In this work, the ERA5 reanalysis dataset was utilized [71,72]. It has been implemented by the European Centre for Medium-Range Weather Forecasts (ECMWF), and it is a dynamical-physical coupled numerical model that uses historical records and 4D-Var data assimilation in order to produce estimations on the climate and its various parameters on a global scale. ERA5 is the product dataset of a long line of reanalysis produced by the ECMWF that was initiated in 1979 and is becoming increasingly sophisticated with higher resolution and more accurate forecasting. Its data are available and can be freely accessed from the Copernicus Climate Data Store (https://cds.climate.copernicus.eu/#! /home (accessed on 7 March 2022)). Compared to the previous ERA-Interim reanalysis, the ERA5 dataset is based on an improved numerical model and data assimilation schemes, and has a more sophisticated resolution with~31 km for the horizontal dimension, 137 levels in the vertical dimension and 1 h for the time scale resolution, making it an adequate source of data for this particular study. In this work, the wind speed time series at a 10 m height above sea level will be analyzed and some of the derived results will be also compared with the corresponding results obtained from the ERA-Interim [35].
The ERA5 dataset has been validated by a number of studies dedicated to this task. In the Mediterranean, a validation focusing on the Tyrrhenian and the Ligurian Seas by downscaling the ERA5 atmospheric reanalysis in order to produce a hindcast for the period 1990-2018 indicated, after comparison with the ERA5 data, that the downscaling improves the hindcast's reliability for wind speed and wind direction results [73]. The estimation of the global wind power potential has been performed using ERA5, whereas the suitability of the ERA5 wind dataset for long-term regional analyses has been highlighted [74]. ERA5 wind speeds have been also compared against lidar measurements at different heights [75]. It was concluded that the ERA5 wind speed data deviate up to 20% from lidar measurements. ERA5 wind and wave data have also been validated after comparison with the ERA-Interim data and advanced scatterometer and offshore platforms measurements obtained at the Caspian Sea. It was concluded that there is a good agreement for wind speed values of 2 m/s or more [76]. Wind speed results over the South Brazil coasts from ERA5 and two different reanalysis datasets were compared to in situ measurements, and ERA5 emerged as the most accurate reanalysis [77]. A global comparison conducted by Sharmar [74] between different reanalysis models' hindcasts (ERA5, MERRA2, ERA-I, CFSR and WAVEWATCH III) for 37 years, accompanied also by a comparison with satellite and buoy measurements, proved again that the ERA5 reanalysis is the most accurate model. Similar conclusions were also drawn by Olauson [78]. ERA5 wave and wind parameters have been evaluated against in situ buoy measurements in the Greek Seas and it was concluded that the ERA5 reanalysis dataset, in general, tends to underestimate the actual values of the examined parameters, but the (site-specific) deviations were found to be acceptable [79]. Let it be noted that site-dependent deviations of numerical model products when compared to measured or remotely sensed data is a very common feature [80]. The ERA5 wind dataset has been also evaluated by comparison against in situ measurements from tall towers in very different offshore and onshore areas of the world and results from other reanalysis products [81]. It has been concluded that "ERA5 scores proved to be finer than those achieved on comparable offshore sites when using higher-resolution regional reanalysis products or even NWP models" [81].

Theoretical Background and Methodology
For the spatial analysis of wind speed and its variability at the annual, monthly and decadal times scales (Sections 4.1 and 4.2), the terminology of Soukissian [35] will be adopted. Using the basic time index t, u(t i ) = u i and θ(t i ) = θ i , i = 1, 2, . . . will represent the time series of wind speed and wind direction, respectively. The scale of t i is hours, while months will be denoted with the letter m and years with the letter j. Accordingly, the annual mean wind speed for a particular year j is with N being the number of observations of the year (in our case N is the number of hours in a year). The mean annual wind speed m u,Y then is estimated as where J is the total number of years. Similarly, the monthly mean wind speed for a particular year j and month m is m u,j,m , m = 1, 2, · · · , 12.
The mean annual variability, denoted as CV (also referred to as MAV), provides a measure of wind speed variability within each year [35] and is estimated as the mean annual coefficient of variation of wind speed. Additionally, the wind speed variability from year to year is estimated by means of the interannual variability I AV, which is the ratio of the standard deviation of the annual mean wind speed to the overall mean wind speed.
For the linear trend identification of wind speed u and the estimation of the corresponding slope in the annual and monthly time scales (Sections 4.1.3 and 4.2.3), the Theil-Sen estimator was used. This estimator was introduced by Theil [82] and extended by Sen [83] and provides a value for the linear slope so that the corresponding Kendall's correlation tau (τ), between u i − λt i and t i equals zero [84], where t i is the time index corresponding to u. Specifically, for the estimation of the Theil-Sen slope, the slopes of all possible combinations of pairs of points y, t, should be estimated, where t i , t j are the time indices corresponding to y i , y j , respectively [84]. The total number of slopes is n(n − 1)/2, where n is the sample size. The median of these slopes λ ij , is the Theil-Sen estimator. This method was preferred in comparison to the ordinary least squares approach, for a variety of reasons as, e.g., the fact that it is resistant to outliers, while at the same time being very robust. The robustness of the estimation has been verified by calculating its breakdown point-the point after which the estimator becomes redundant-as being 0.293, indicating a robust estimator [85]. Moreover, the Theil-Sen estimator has high asymptotic efficiency and proven unbiasedness for continuous error distribution, while it is considered a good estimator for mean squares error, in comparison to other slope estimators.
Regarding the wind direction, the statistical analysis for vectorial data differs from linear data in the sense that one needs to utilize a directional statistics approach. This kind of data has a starting point and a direction (clockwise and counter-clockwise), as well as a value, making it intrinsically different from linear data, especially with respect to notions such as "larger", "increase", etc. The way to compare such data derives from deciding one direction as being positive. An introduction to the most important statistical parameters of directional variables (mean value, circular variance, etc.) is provided by Soukissian [45]. Moreover, the angular distance ∆θ between two directional variables θ 1 and θ 2 is defined [86] as: The angular distance takes values in [−π, π]. The variable ∆θ is set as positive for the clockwise direction and negative for the counter-clockwise one.
Finally, the linear-circular correlation coefficient r 2 uθ , indicates the degree of relationship between a linear (i.e., wind speed u) and a directional (i.e., wind direction θ) variable [35].

Numerical Results
This section includes the analysis of the wind speed and wind direction in the MS, on the annual, monthly and interannual time scales. Some of the results that are presented in this section will be compared with the corresponding results of the ERA-Interim [35] in order to reveal potential discrepancies between the two datasets.
Before proceeding to the presentation and discussion of the results, let us describe the main intense wind regimes that are encountered in the Mediterranean Sea. See also, Lionello [64] for a detailed description of the wind systems in the Mediterranean Sea. The most important wind systems are of local nature: for example, Sirocco is a wind that blows from the Sahara Desert towards the Mediterranean, reaching high speeds, especially during spring and autumn, Tramontana is a northern wind that appears in many Mediterranean countries, Mistral is a strong northwesterly wind that blows across the Rhône River valley into the Gulf of Lion in France, the Etesians are northerly seasonal winds that blow across the Aegean Sea during summer and Bora is a north-easterly wind blowing in the Adriatic Sea, especially from November to March. These wind systems shape the wind climate of the basin, and their appearance can be traced in the following results. Note also that the same wind systems may have different names in specific regions, e.g., the Sirocco is named Khamsin in Egypt.
In Figure 1, the map of the Mediterranean basin along with the areas that are often mentioned in the text are provided.

Annual and Interannual Time Scales
In the following Figure 2, the areas with the highest mean wind speed emerge, namely the Gulf of Lion, the Aegean Sea and the sea surrounding the western and eastern

Annual and Interannual Time Scales
In the following Figure 2, the areas with the highest mean wind speed emerge, namely the Gulf of Lion, the Aegean Sea and the sea surrounding the western and eastern coasts of Crete Isl., as well as the Alboran Sea. These areas are the most wind energetic areas in the Mediterranean and therefore, they are promising candidates for offshore wind energy development. The highest mean wind speed values occur in the Gulf of Lion (7.80 m/s) and the offshore area west of Crete Isl. (7.62 m/s). The spatial distribution patterns of the wind direction are similar to previous results [35,67,69]. For the areas where the wind speed values appear to be the largest, the wind direction is NW (Mistral wind) in the Gulf of Lion and NE in the N Aegean which then turns to NW in the south Aegean Sea. Another easily identified wind system is the NE Bora (blowing from the Dinaric Alps in the Adriatic Sea). The wind climate in the Aegean Sea is largely shaped by NE winds during the winter and Etesians during the summer, and in the Gulf of Lion by the Mistral and the Tramontana wind systems. Results from previous studies indicate that the areas with the highest wind speed, as well as the wind direction patterns, are the same.

Annual and Interannual Time Scales
In the following Figure 2, the areas with the highest mean wind speed emerge, namely the Gulf of Lion, the Aegean Sea and the sea surrounding the western and eastern coasts of Crete Isl., as well as the Alboran Sea. These areas are the most wind energetic areas in the Mediterranean and therefore, they are promising candidates for offshore wind energy development. The highest mean wind speed values occur in the Gulf of Lion (7.80 m/s) and the offshore area west of Crete Isl. (7.62 m/s). The spatial distribution patterns of the wind direction are similar to previous results [35,67,69]. For the areas where the wind speed values appear to be the largest, the wind direction is NW (Mistral wind) in the Gulf of Lion and NE in the N Aegean which then turns to NW in the south Aegean Sea. Another easily identified wind system is the NE Bora (blowing from the Dinaric Alps in the Adriatic Sea). The wind climate in the Aegean Sea is largely shaped by NE winds during the winter and Etesians during the summer, and in the Gulf of Lion by the Mistral and the Tramontana wind systems. Results from previous studies indicate that the areas with the highest wind speed, as well as the wind direction patterns, are the same. In order to reveal the important differences that occur in the mean annual wind speed estimation between ERA5 and ERA-Interim, in Figure 3, the relative difference Δ between the mean annual wind speed obtained from ERA-Interim and ERA5 is shown. In order to reveal the important differences that occur in the mean annual wind speed estimation between ERA5 and ERA-Interim, in Figure 3, the relative difference ∆u AN between the mean annual wind speed obtained from ERA-Interim and ERA5 is shown.
Wind 2022, 2, FOR PEER REVIEW 8 Figure 3. Relative difference of the mean annual wind speed obtained from ERA5 and ERA-Interim datasets.
Δ is defined as follows: where , and , denote the mean annual wind speed obtained from ERA5 and ERA-Interim, respectively. The common reference period for both datasets is 1979-2014. From this figure, it can be seen that for extended areas of the Mediterranean Sea, the absolute relative difference is greater than 25%. Important relative differences in the mean ∆u AN is defined as follows: Wind 2022, 2 520 where u AN,ERA5 and u AN,I NT denote the mean annual wind speed obtained from ERA5 and ERA-Interim, respectively. The common reference period for both datasets is 1979-2014. From this figure, it can be seen that for extended areas of the Mediterranean Sea, the absolute relative difference is greater than 25%. Important relative differences in the mean annual estimates occur in the Aegean and Adriatic Seas, as well as in several nearshore areas. Specifically, in the N Adriatic Sea, the relative difference is of the order of up to −50%, while in the W Aegean Sea the relative difference is of the order of up to 35%.

Variability Characteristics
The variability characteristics are provided through the mean annual variability, i.e., the variability within a year (MAV), and the interannual variability, i.e., the year-to-year variability, (I AV). The left plot of Figure 4 illustrates the spatial distribution of the MAV.
The highest values of the MAV appear in the Ligurian Sea and especially in the Gulf of Genoa, an area prone to regular cyclones, the N Levantine Basin and especially in the Gulf of Antalya, where tornadoes are frequent [87], the E Tyrrhenian Sea and the N Aegean Sea, while relatively high values of MAV are also encountered in the Adriatic Sea, especially in the Gulf of Venice. The southern part of the basin exhibits lower MAV values, and the general spatial pattern consists of declining values from the north to the south. The results presented here are in agreement with the results referring to wind power density, where the ERA5 dataset was used [69], and with the results of the Blended SeaWinds product for the period 1995-2011 [67]. On the other hand, results based on the ERA-Interim dataset [35] demonstrate that the highest MAV occurs in the Adriatic Sea with a value of 68%. and ERA-Interim, respectively. The common reference period for both datasets is 1979-2014. From this figure, it can be seen that for extended areas of the Mediterranean Sea, the absolute relative difference is greater than 25%. Important relative differences in the mean annual estimates occur in the Aegean and Adriatic Seas, as well as in several nearshore areas. Specifically, in the N Adriatic Sea, the relative difference is of the order of up to -50%, while in the W Aegean Sea the relative difference is of the order of up to 35%.

Variability Characteristics
The variability characteristics are provided through the mean annual variability, i.e., the variability within a year ( ), and the interannual variability, i.e., the year-to-year variability, ( ). The left plot of Figure 4 illustrates the spatial distribution of the . The highest values of the appear in the Ligurian Sea and especially in the Gulf of Genoa, an area prone to regular cyclones, the N Levantine Basin and especially in the Gulf of Antalya, where tornadoes are frequent [87], the E Tyrrhenian Sea and the N Aegean Sea, while relatively high values of are also encountered in the Adriatic Sea, especially in the Gulf of Venice. The southern part of the basin exhibits lower values, and the general spatial pattern consists of declining values from the north to the south. The results presented here are in agreement with the results referring to wind power density, where the ERA5 dataset was used [69], and with the results of the Blended SeaWinds product for the period 1995-2011 [67]. On the other hand, results based on the ERA-Interim dataset [35] demonstrate that the highest MAV occurs in the Adriatic Sea with a value of 68%.  The right plot of Figure 4 illustrates the spatial distribution of the I AV. Again, the interannual variability exhibits lower values in the southern part of the basin (except for the coasts of Algeria). The largest value (6%) appears in the S Tyrrhenian sub-basin and the N Levantine Basin. The regions with the overall largest I AV are the Balearic, the Ligurian and the Tyrrhenian Seas, as well as the N Aegean and the S coasts of Turkey. The lowest value of I AV appears off the W coasts of Crete Isl. Note that the results based on ERA-Interim [35] demonstrate higher maximum values of I AV (7.25%), appearing in the W Ligurian Sea. Figure 5 depicts the angular variance of the mean annual wind direction. The angular variance is a measure of the spread of the mean annual wind direction values around the overall annual mean. The smaller the angular variance, the more concentrated the circular distribution. The majority of the Mediterranean areas presents values close to 0 (0-0.02) with a few specific areas with values above 0.03. This result suggests that the mean annual wind direction remains, to a large extent, stable. The regions that exhibit the highest angular variance are found in the western part of the Mediterranean Sea and are very localized, e.g., the middle of the Alboran Sea, the W Algerian basin, the Balearic Sea (around and especially north of Mallorca) and the N Ligurian Sea (off the coasts of Monaco and offshore the Gulf of Genoa). Values of angular variance close to 1 suggest that the mean annual wind direction is highly variable. Lower values around 0.4 appear as well in specific spots in the Adriatic Sea and the Gulf of Mersin. These results are in general similar to the results by Soukissian [35], in terms of the areas that present the highest values, however, the more sophisticated the ERA5 dataset indicates that there are more areas with high values of Wind 2022, 2 521 angular variance (higher than 0.4) in the Alboran, Balearic and Algerian basins, that had not emerged previously. gular variance are found in the western part of the Mediterranean Sea and are very localized, e.g., the middle of the Alboran Sea, the W Algerian basin, the Balearic Sea (around and especially north of Mallorca) and the N Ligurian Sea (off the coasts of Monaco and offshore the Gulf of Genoa). Values of angular variance close to 1 suggest that the mean annual wind direction is highly variable. Lower values around 0.4 appear as well in specific spots in the Adriatic Sea and the Gulf of Mersin. These results are in general similar to the results by Soukissian [35], in terms of the areas that present the highest values, however, the more sophisticated the ERA5 dataset indicates that there are more areas with high values of angular variance (higher than 0.4) in the Alboran, Balearic and Algerian basins, that had not emerged previously. In Figure 6, the angular distance between the mean decadal wind direction for the 1st (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)) and 2nd (1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)   In Figure 6, the angular distance between the mean decadal wind direction for the 1st (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)) and 2nd (1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)  high values of angular variance (higher than 0.4) in the Alboran, Balearic and Algerian basins, that had not emerged previously. In Figure 6, the angular distance between the mean decadal wind direction for the 1st (1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)) and 2nd (1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)  The largest differences are observed for the decadal mean wind direction between the 3rd and 4th, as well as between the 1st and last decades. It is interesting to note that for some very particular regions positive and negative angular distances are simultaneously present in the form of dipoles. These areas are the entire Balearic Sea, the Ligurian Sea (offshore the Gulf of Genoa), across the Algerian coasts and the southern coasts of Turkey. Note also that the areas exhibiting high values of angular variance are largely the same as the areas exhibiting the highest decadal angular distances.

Association between Wind Speed and Direction
In Figure 7, the spatial distribution of the mean annual linear-circular correlation coefficient between the wind speed and direction is presented. The highest values (0.56) of appear in an extended offshore area in the Gulf of Lion, suggesting a positive as- The largest differences are observed for the decadal mean wind direction between the 3rd and 4th, as well as between the 1st and last decades. It is interesting to note that for some very particular regions positive and negative angular distances are simultaneously present in the form of dipoles. These areas are the entire Balearic Sea, the Ligurian Sea (offshore the Gulf of Genoa), across the Algerian coasts and the southern coasts of Turkey. Note also that the areas exhibiting high values of angular variance are largely the same as the areas exhibiting the highest decadal angular distances.

Association between Wind Speed and Direction
In Figure 7, the spatial distribution of the mean annual linear-circular correlation coefficient r 2 uθ between the wind speed and direction is presented. The highest values (0.56) of r 2 uθ appear in an extended offshore area in the Gulf of Lion, suggesting a positive association between wind speed and wind direction of the Tramontana and Mistral wind regimes. Other areas with relatively high circular correlation coefficients are the Alboran, the Balearic and S Tyrrhenian Seas, the Algerian basin, the S Levantine Basin, the N Adriatic and the Aegean Seas and the W coasts of Lebanon. The areas with the lowest values are the E Levantine Basin, the N Tyrrhenian, the Central Mediterranean sub-basin, the Ionian Sea and the S Adriatic. This result is in agreement with previous work [35] where the same regions emerged as having the highest correlation coefficient values, but the more detailed analysis provided by the ERA5 dataset gives higher values. Moreover, the NE Alboran Sea emerged as having one of the highest correlation coefficient values, which it had previously not. same as the areas exhibiting the highest decadal angular distances.

Association between Wind Speed and Direction
In Figure 7, the spatial distribution of the mean annual linear-circular correlation coefficient between the wind speed and direction is presented. The highest values (0.56) of appear in an extended offshore area in the Gulf of Lion, suggesting a positive association between wind speed and wind direction of the Tramontana and Mistral wind regimes. Other areas with relatively high circular correlation coefficients are the Alboran, the Balearic and S Tyrrhenian Seas, the Algerian basin, the S Levantine Basin, the N Adriatic and the Aegean Seas and the W coasts of Lebanon. The areas with the lowest values are the E Levantine Basin, the N Tyrrhenian, the Central Mediterranean sub-basin, the Ionian Sea and the S Adriatic. This result is in agreement with previous work [35] where the same regions emerged as having the highest correlation coefficient values, but the more detailed analysis provided by the ERA5 dataset gives higher values. Moreover, the NE Alboran Sea emerged as having one of the highest correlation coefficient values, which it had previously not.

Wind Speed Trend
In Figure 8, the linear slope for the mean wind speed on the annual scale is presented. The slopes are indicative of the wind speed tendency to increase, decrease or remain constant over the time span under analysis.

Wind Speed Trend
In Figure 8, the linear slope for the mean wind speed on the annual scale is presented. The slopes are indicative of the wind speed tendency to increase, decrease or remain constant over the time span under analysis. The highest positive slopes appear in the area between Crete Isl. and the African coasts, and then in the E Levantine Basin, the E Balearic Sea and offshore the W and E coasts of Sardinia Isl. The overall highest slope (0.087 m/s/year) appears offshore the SW coasts of Crete Isl. Previous results [35] showed the largest positive slope to appear in the Ionian Sea. On the other hand, the Ligurian Sea, the Central Mediterranean sub-basin and the rest of the Levantine Basin exhibit negative slopes. The minimum slope (−0.011 m/s/year) appears in the Aegean Sea, offshore the island of Rodos. The negative slopes of the Ligurian Sea, the central Mediterranean sub-basin and the Levantine Basin, as well as the increasing trends in the W Levantine Basin, the Ionian and the E Balearic-W Tyrrhenian Seas, generally match the results of Soukissian [35]. The central Aegean was expected to have decreasing slopes [35,88,89], yet the results presented here are in agreement with those provided by the recent work of Weiss [90].
As another example of the discrepancies between the ERA-Interim and ERA5 datasets, in Figure 9 the Theil-Sen linear slopes of the mean annual wind speeds for the The highest positive slopes appear in the area between Crete Isl. and the African coasts, and then in the E Levantine Basin, the E Balearic Sea and offshore the W and E coasts of Sardinia Isl. The overall highest slope (0.087 m/s/year) appears offshore the SW coasts of Crete Isl. Previous results [35] showed the largest positive slope to appear in the Ionian Sea. On the other hand, the Ligurian Sea, the Central Mediterranean sub-basin and the rest of the Levantine Basin exhibit negative slopes. The minimum slope (−0.011 m/s/year) appears in the Aegean Sea, offshore the island of Rodos. The negative slopes of the Ligurian Sea, the central Mediterranean sub-basin and the Levantine Basin, as well as the increasing trends in the W Levantine Basin, the Ionian and the E Balearic-W Tyrrhenian Seas, generally match the results of Soukissian [35]. The central Aegean was expected to have decreasing Wind 2022, 2 523 slopes [35,88,89], yet the results presented here are in agreement with those provided by the recent work of Weiss [90].
As another example of the discrepancies between the ERA-Interim and ERA5 datasets, in Figure 9 the Theil-Sen linear slopes of the mean annual wind speeds for the locations 43.25 • N, 7.75 • E (left) and 37.50 • N, 19.25 • E (right) are depicted. Let it be noted that the wind speed trend of these locations has been assessed by using the ERA-Interim dataset [35].
coasts of Sardinia Isl. The overall highest slope (0.087 m/s/year) appears offshore the SW coasts of Crete Isl. Previous results [35] showed the largest positive slope to appear in the Ionian Sea. On the other hand, the Ligurian Sea, the Central Mediterranean sub-basin and the rest of the Levantine Basin exhibit negative slopes. The minimum slope (−0.011 m/s/year) appears in the Aegean Sea, offshore the island of Rodos. The negative slopes of the Ligurian Sea, the central Mediterranean sub-basin and the Levantine Basin, as well as the increasing trends in the W Levantine Basin, the Ionian and the E Balearic-W Tyrrhenian Seas, generally match the results of Soukissian [35]. The central Aegean was expected to have decreasing slopes [35,88,89], yet the results presented here are in agreement with those provided by the recent work of Weiss [90].
As another example of the discrepancies between the ERA-Interim and ERA5 datasets, in Figure 9 the Theil-Sen linear slopes of the mean annual wind speeds for the locations 43.25° N, 7.75° E (left) and 37.50° N, 19.25° E (right) are depicted. Let it be noted that the wind speed trend of these locations has been assessed by using the ERA-Interim dataset [35]. Despite the fact that the indicated slopes are of the same sign, the ones that are provided by the ERA-Interim dataset are one order of magnitude different from the corresponding slopes provided by the ERA5 dataset. In addition, in Figure 10, the Theil-Sen linear slope for the mean annual wind speed obtained from the ERA-Interim and ERA5 datasets for the location 37.75° N, 25.00° E is presented. In this case, there is a difference in the sign of the slope (ERA5 provides a slightly positive slope, while ERA-Interim provides a negative slope), as well as important differences in the annual statistics. Specifically, the mean annual, standard deviation, minimum and maximum values from ERA5 Despite the fact that the indicated slopes are of the same sign, the ones that are provided by the ERA-Interim dataset are one order of magnitude different from the corresponding slopes provided by the ERA5 dataset. In addition, in Figure 10, the Theil-Sen linear slope for the mean annual wind speed obtained from the ERA-Interim and ERA5 datasets for the location 37.75 • N, 25.00 • E is presented. In this case, there is a difference in the sign of the slope (ERA5 provides a slightly positive slope, while ERA-Interim provides a negative slope), as well as important differences in the annual statistics.  Nevertheless, the qualitative differences between the ERA-Interim dataset that has previously been used, compared to the much more sophisticated ERA5 dataset, have to be taken into account when examining the differences in those sets of results.

Extreme Wind Events
In Figure 11, the linear slopes for the annual 95th (left) and 99th (right) quantiles of  Nevertheless, the qualitative differences between the ERA-Interim dataset that has previously been used, compared to the much more sophisticated ERA5 dataset, have to be taken into account when examining the differences in those sets of results.

Extreme Wind Events
In Figure 11, the linear slopes for the annual 95th (left) and 99th (right) quantiles of wind speed are presented. Figure 10. Linear slope of the mean annual wind speed for the location 37.75° N. 25.00° E, as obtained from the ERA-Interim and ERA5 datasets.
Nevertheless, the qualitative differences between the ERA-Interim dataset that has previously been used, compared to the much more sophisticated ERA5 dataset, have to be taken into account when examining the differences in those sets of results.

Extreme Wind Events
In Figure 11, the linear slopes for the annual 95th (left) and 99th (right) quantiles of wind speed are presented. By comparing Figure 11 with Figure 8, it can be seen that the rate of change (increase or decrease) of the extreme wind speeds is clearly greater (in the absolute sense) than the corresponding rate of the annual mean wind speeds. The areas where the 99th annual percentiles of wind speeds tend to increase at the fastest rates are the coasts of the Middle East, the central Aegean Sea, the sea of Alboran and the W Ionian Sea. Areas where extreme wind speeds tend to decrease are the S coasts of Turkey, the area between the S Ionian and the Libyan Sea, the NW offshore area of Corsica Isl., the N coasts of Sicily Isl. and the Ligurian Sea.
Extreme wind events have been further assessed. Since an "extreme event" is a local phenomenon that is clearly defined and related to the wind climate of a particular area, we introduce the following definition for extreme events: Let , , = 1,2, ⋯ denote the annual maxima series for a particular location. A local extreme event in this area will be considered any event with min , . Accordingly, in Figure 12, the linear slopes for the annual number of extreme events of wind speed are depicted. It is clear that the By comparing Figure 11 with Figure 8, it can be seen that the rate of change (increase or decrease) of the extreme wind speeds is clearly greater (in the absolute sense) than the corresponding rate of the annual mean wind speeds. The areas where the 99th annual percentiles of wind speeds tend to increase at the fastest rates are the coasts of the Middle East, the central Aegean Sea, the sea of Alboran and the W Ionian Sea. Areas where extreme wind speeds tend to decrease are the S coasts of Turkey, the area between the S Ionian and the Libyan Sea, the NW offshore area of Corsica Isl., the N coasts of Sicily Isl. and the Ligurian Sea.
Extreme wind events have been further assessed. Since an "extreme event" is a local phenomenon that is clearly defined and related to the wind climate of a particular area, we introduce the following definition for extreme events: Let u j,max , j = 1, 2, · · · J denote the annual maxima series for a particular location. A local extreme event in this area will be considered any event with u ≥ min j u j,max . Accordingly, in Figure 12, the linear slopes for the annual number of extreme events of wind speed are depicted. It is clear that the extended areas characterized by increasing the number of occurrences of extreme wind events per year are the central Aegean Sea, the Alboran Sea, the W Ionian Sea and the SE part of the Levantine Basin. Taking also into account the results of Figure 11, it can be concluded that these areas are characterized by both increased intensities and frequencies of occurrence of extreme wind speed events.
Wind 2022, 2, FOR PEER REVIEW 13 extended areas characterized by increasing the number of occurrences of extreme wind events per year are the central Aegean Sea, the Alboran Sea, the W Ionian Sea and the SE part of the Levantine Basin. Taking also into account the results of Figure 11, it can be concluded that these areas are characterized by both increased intensities and frequencies of occurrence of extreme wind speed events.

Monthly Time Scale
In Figure 13, the mean wind speed and wind direction for the monthly scale are presented for May, August, October and December, since these months are characterized by distinct features as explained below.

Monthly Time Scale
In Figure 13, the mean wind speed and wind direction for the monthly scale are presented for May, August, October and December, since these months are characterized by distinct features as explained below. The month with the highest wind speeds overall is February, and the overall highest mean wind speeds are encountered in the Gulf of Lion and the central Aegean Sea. The overall lowest wind speeds occur in August, notwithstanding the existence of the Etesians (a seasonal phenomenon that is discernible in Figure 13) in the Aegean Sea. The monthly spatial wind patterns that are presented by Laurila [34] are in agreement with the results presented here. We will focus now on five areas that exhibit some interesting seasonal patterns. In the Aegean Sea, during the summer (JJA), the wind direction is NE at the N Aegean turning to NW at the central and south Aegean Sea. This season is dominated by the Etesian winds which blow steadily during the summer [91][92][93]. In the Gulf of Lion, steady NW winds (Mistral) are present especially during the winter [62]. In the north part of the Ligurian Sea, there is a seasonal wind direction pattern. From September to March, there is a circular pattern over Corsica Isl. with winds that turn from SE off the E coasts to SW off the W coasts. From April to August, the main winds are SW and there appears no circular pattern. In the Levantine Basin and in particular, on the E coasts of Cyprus, there appears a circular pattern that advances towards the north and is not visible from May to September. In January, the pattern is located offshore the eastern coasts of Cyprus and moves northerly, until April when it is barely visible on the southern coasts of Turkey. There it appears again in October and it moves southerly until December when it is again located in the east of Cyprus. Finally, in the Adriatic Sea during the summer, the wind direction in the central part is NW while for the rest of the year, the direction changes from the east part to the west, from SE to NW, creating a semi-circular pattern. During the months of April, and October until December, there are also SW winds that complete the circular pattern. Figure 14 depicts the monthly mean angular variance. The highest values are observed during March, April, October and November. The months with the maximum angular variance values are the transitional ones, as they are representative of the transition from spring to summer and from autumn to winter. The general behavior also matches the one of previous results, with the exception of the four transitional months (March, April, October and November), where more extended areas of the basin have angular variance values that exceed 0.35. The period from October to April is characterized by higher mean angular variance values, while from May to September, the mean angular variance decreases. The reason for this pattern is that during higher temperature months (late spring-summer-early autumn) there are steadier wind pressure patterns [94,95]. During The month with the highest wind speeds overall is February, and the overall highest mean wind speeds are encountered in the Gulf of Lion and the central Aegean Sea. The overall lowest wind speeds occur in August, notwithstanding the existence of the Etesians (a seasonal phenomenon that is discernible in Figure 13) in the Aegean Sea. The monthly spatial wind patterns that are presented by Laurila [34] are in agreement with the results presented here.

Variability Characteristics
We will focus now on five areas that exhibit some interesting seasonal patterns. In the Aegean Sea, during the summer (JJA), the wind direction is NE at the N Aegean turning to NW at the central and south Aegean Sea. This season is dominated by the Etesian winds which blow steadily during the summer [91][92][93]. In the Gulf of Lion, steady NW winds (Mistral) are present especially during the winter [62]. In the north part of the Ligurian Sea, there is a seasonal wind direction pattern. From September to March, there is a circular pattern over Corsica Isl. with winds that turn from SE off the E coasts to SW off the W coasts. From April to August, the main winds are SW and there appears no circular pattern. In the Levantine Basin and in particular, on the E coasts of Cyprus, there appears a circular pattern that advances towards the north and is not visible from May to September. In January, the pattern is located offshore the eastern coasts of Cyprus and moves northerly, until April when it is barely visible on the southern coasts of Turkey. There it appears again in October and it moves southerly until December when it is again located in the east of Cyprus. Finally, in the Adriatic Sea during the summer, the wind direction in the central part is NW while for the rest of the year, the direction changes from the east part to the west, from SE to NW, creating a semi-circular pattern. During the months of April, and October until December, there are also SW winds that complete the circular pattern. Figure 14 depicts the monthly mean angular variance. The highest values are observed during March, April, October and November. The months with the maximum angular variance values are the transitional ones, as they are representative of the transition from spring to summer and from autumn to winter. The general behavior also matches the one of previous results, with the exception of the four transitional months (March, April, October and November), where more extended areas of the basin have angular variance values that exceed 0.35. The period from October to April is characterized by higher mean angular variance values, while from May to September, the mean angular variance decreases. The reason for this pattern is that during higher temperature months (late spring-summer-early autumn) there are steadier wind pressure patterns [94,95]. During the low angular variance period, only a few specified locations have values over 0.30. Moreover, the locations that were expected to exhibit high angular variance values (Gulf of Lion, Aegean, Adriatic and Levantine) confirm this expectation. Figure 14 depicts the monthly mean angular variance. The highest values are observed during March, April, October and November. The months with the maximum angular variance values are the transitional ones, as they are representative of the transition from spring to summer and from autumn to winter. The general behavior also matches the one of previous results, with the exception of the four transitional months (March, April, October and November), where more extended areas of the basin have angular variance values that exceed 0.35. The period from October to April is characterized by higher mean angular variance values, while from May to September, the mean angular variance decreases. The reason for this pattern is that during higher temperature months (late spring-summer-early autumn) there are steadier wind pressure patterns [94,95]. During the low angular variance period, only a few specified locations have values over 0.30. Moreover, the locations that were expected to exhibit high angular variance values (Gulf of Lion, Aegean, Adriatic and Levantine) confirm this expectation.

Association between Wind Speed and Direction
In Figure 15, the mean monthly linear-circular correlation coefficient between wind speed and direction for August is shown. This month was selected as it represents the highest monthly correlation coefficient value (0.78), which appears in the Aegean Sea and the W Levantine Basin.

Association between Wind Speed and Direction
In Figure 15, the mean monthly linear-circular correlation coefficient r 2 uθ between wind speed and direction for August is shown. This month was selected as it represents the highest monthly correlation coefficient value (0.78), which appears in the Aegean Sea and the W Levantine Basin.
For the month of August, a comparison between this and previously reported results [35] indicates that the basin presents similar behaviors. The areas with the highest linear-circular correlation coefficient values are the Aegean Sea, the Levantine Basin and the Gulf of Lion. The Aegean Sea has the highest correlation coefficient values during the summer months (JJA) due to the presence of steady Etesian winds over the area. The Etesians are intense persistent northern winds; therefore, the wind speed and direction are two highly associated variables in this case. On the other hand, the W Algerian basin, the Ionian and the Adriatic Seas, as well as the S Levantine Basin, have the lowest values.

Association between Wind Speed and Direction
In Figure 15, the mean monthly linear-circular correlation coefficient between wind speed and direction for August is shown. This month was selected as it represents the highest monthly correlation coefficient value (0.78), which appears in the Aegean Sea and the W Levantine Basin. For the month of August, a comparison between this and previously reported results [35] indicates that the basin presents similar behaviors. The areas with the highest linearcircular correlation coefficient values are the Aegean Sea, the Levantine Basin and the Gulf of Lion. The Aegean Sea has the highest correlation coefficient values during the summer months (JJA) due to the presence of steady Etesian winds over the area. The Etesians are intense persistent northern winds; therefore, the wind speed and direction are two highly associated variables in this case. On the other hand, the W Algerian basin, the Ionian and the Adriatic Seas, as well as the S Levantine Basin, have the lowest values.

Wind Speed Trend
In Figure 16, the linear slopes for the monthly mean wind speed are presented. The highest (0.035 m/s/year) and lowest (−0.045 m/s/year) slopes are encountered during May and April, respectively, both in the Gulf of Lion. In January, the highest slopes appear in the eastern part and the lowest in the western part of the basin, while in February it takes place the opposite. In March, high slopes are encountered in the Alboran Sea and across the African coasts. In April, the lowest slopes can be found in the extended offshore area between the Gulf of Lion and Sardinia Isl. and the highest in the Alboran and the Aegean Seas. In May, a large part of the basin is characterized by positive slopes. During June, Figure 15. Linear-circular correlation coefficient of wind speed and direction for August.

Wind Speed Trend
In Figure 16, the linear slopes for the monthly mean wind speed are presented. The highest (0.035 m/s/year) and lowest (−0.045 m/s/year) slopes are encountered during May and April, respectively, both in the Gulf of Lion. In January, the highest slopes appear in the eastern part and the lowest in the western part of the basin, while in February it takes place the opposite. In March, high slopes are encountered in the Alboran Sea and across the African coasts. In April, the lowest slopes can be found in the extended offshore area between the Gulf of Lion and Sardinia Isl. and the highest in the Alboran and the Aegean Seas. In May, a large part of the basin is characterized by positive slopes. During June, July and August, the slope spatial distribution is very similar with very low absolute slopes (close to 0) at the W Mediterranean and increasing positive slopes in the Aegean Sea (reaching values up to 0.034 m/s/month during August) due to the Etesian winds. October, November and December are characterized by rather mild slopes (either positive or negative), while October is a transitional month.
In Figure 17, the linear slopes for the monthly 99th percentiles of wind speed are depicted. The monthly 95th percentiles of wind speed have been also estimated but are not shown here. Clearly, the slopes are in this case greater than the corresponding slopes of the mean monthly values. In May, the sea area between Sardinia Isl. and the Gulf of Lion exhibits the overall maximum (positive) slopes. The slope patterns exhibit some very distinct regional features for the western, central and eastern Mediterranean. Specifically, for the western Mediterranean from January to March the 99% percentile slopes are over 0.04 for the majority of the area (except for January in the Tyrrhenian Sea). In April, the slope values abruptly decrease taking negative values, where the sea area between Sardinia Isl. and the Gulf of Lion and the Tyrrhenian Sea exhibits the overall minimum slopes, only to rise abruptly again to a maximum in May. This monthly pattern of abrupt changes from positive to negative slopes and vice versa last until December for the western part of the basin. In the central part, the months with the lowest slopes are January, March and December, while the maximum (positive) slopes appear in September and October. Finally, in the Eastern part of the basin, the months of January, May and September present increased slope values. On the contrary, July, August and November have very low values. Specifically, in the northern part of the Aegean Sea, the slopes take higher values than the rest of the east basin during August and December. For the majority of the months, the eastern part of the basin seems to present decreased slopes in comparison to the west part, with the only exception to this behavior being the month of January.
July and August, the slope spatial distribution is very similar with very low absolute slopes (close to 0) at the W Mediterranean and increasing positive slopes in the Aegean Sea (reaching values up to 0.034 m/s/month during August) due to the Etesian winds. October, November and December are characterized by rather mild slopes (either positive or negative), while October is a transitional month. In Figure 17, the linear slopes for the monthly 99th percentiles of wind speed are depicted. The monthly 95th percentiles of wind speed have been also estimated but are not shown here. Clearly, the slopes are in this case greater than the corresponding slopes of the mean monthly values. In May, the sea area between Sardinia Isl. and the Gulf of Lion exhibits the overall maximum (positive) slopes. The slope patterns exhibit some very distinct regional features for the western, central and eastern Mediterranean. Specifically, for the western Mediterranean from January to March the 99% percentile slopes are over 0.04 for the majority of the area (except for January in the Tyrrhenian Sea). In April, the slope values abruptly decrease taking negative values, where the sea area between Sardinia Isl. and the Gulf of Lion and the Tyrrhenian Sea exhibits the overall minimum slopes, December, while the maximum (positive) slopes appear in September and October. Finally, in the Eastern part of the basin, the months of January, May and September present increased slope values. On the contrary, July, August and November have very low values. Specifically, in the northern part of the Aegean Sea, the slopes take higher values than the rest of the east basin during August and December. For the majority of the months, the eastern part of the basin seems to present decreased slopes in comparison to the west part, with the only exception to this behavior being the month of January.

Discussion
Several studies have indicated that the Mediterranean is a region important in its cyclogenetic characteristics, as it is impacted largely by the North Atlantic Oscillation and as indicated by the above-mentioned areas is favorable for cyclogenesis [96]. Moreover, the E Mediterranean basin has been distinguished for its cyclogenesis [65], which is also

Discussion
Several studies have indicated that the Mediterranean is a region important in its cyclogenetic characteristics, as it is impacted largely by the North Atlantic Oscillation and as indicated by the above-mentioned areas is favorable for cyclogenesis [96]. Moreover, the E Mediterranean basin has been distinguished for its cyclogenesis [65], which is also the area where present results locate the maximum angular distances. A relationship between wind direction angular distance and cyclogenesis seems possible given the previous results and mentioned research, however, this finding is not solidly established here. Climate change, which affects wind direction as indicated by the decadal time scale differentiations, is bound to also affect the cyclogenesis of the region. However, it is still not clear in what way this might happen since there have been studies indicating that cyclones may occur more often [97] and others that predict a decrease in cyclone occurrence [98,99]. Nonetheless, it has been the result of previous research that cyclogenesis in the area is affected seasonally by the large-scale wind circulation there [100]. The potential identification of the areas with intense cyclogenesis using the angular distance, as well as any association between them, is the subject of other research by the authors. While the impact of wind direction is not yet fully understood, anthropogenic climate change is another important factor that will potentially affect both wind speed and direction in the coming decades. Studies have shown slight alterations in wind speed and we can expect this behavior to continue as climate change effects get increasingly intense.
The comparison analyzed above, of collocated ERA-Interim and ERA5 datasets results, revealed important differences in the wind speed and direction statistical parameters. In this work, only comparisons between the mean annual wind speed and the Theil-Sen linear trend line have been analytically shown, but relevant differences have been observed for all estimated quantities in all temporal scales. This was expected due to the improved numerical model and data assimilation schemes, as well as the better temporal and spatial resolution of the ERA5 product. Clearly, for a finer and more detailed analysis of the regional wind climate in the Mediterranean Sea, especially in nearshore areas, higher spatial resolution data are required. This can be achieved, e.g., by downscaling the ERA5 wind fields. Moreover, the definitive version of ERA5 hourly data on single levels from 1950 to 1978 will allow a more robust multi-decadal analysis of wind speed and direction and provide more insight into the long-term behavior of mean and extreme wind speeds, as well as the variability of wind direction in the Mediterranean basin.

Conclusions
The mean annual and the interannual variability of wind speed are reduced in the southern parts of the Mediterranean basin, compared to the corresponding northern parts. The angular variance of the mean annual wind direction exhibits an analogous behavior. The areas characterized by large variability as regards wind speed do not, in general, coincide with the areas of increased angular variance values. This is also confirmed by the spatial distribution of the mean annual linear-circular correlation coefficient between wind speed and direction, which assumes values below 0.3 and close to 0 for the majority of the basin areas. Decadal changes in wind direction, in the form of angular distance, are presented here for the first time as far as the authors are aware. The most extreme angular distances have been observed in very particular areas such as the Balearic and Alboran Seas, the coasts of Algeria, the Ligurian Sea (at the offshore area west of the Gulf of Genoa) and the northeastern part of the Levantine Sea.
The slopes of the annual mean wind speed are rather mild (either positive or negative) for the entire basin. The highest value of the slope (0.087 m/s/y) that appears on the southwestern coasts of Crete Isl. corresponds to a mean increase in wind speed of 3.65 m/s during the examined 42-year period, but this is a rather exceptional case. Additionally, the linear slopes of the extreme wind speed have been also examined: firstly, the slopes for the 95th and 99th percentiles of wind speed are clearly greater than the corresponding ones for the mean wind speed, suggesting that extreme wind states tend to increase in magnitude at the long-term time scale. Secondly, site-specific extreme wind events have also been defined and assessed. The areas with an increasing number of extreme wind events occurrences, except for the Alboran Sea, are located in the eastern Mediterranean, namely the north and central Aegean Sea, western Ionian Sea and the southeastern part of the Levantine Basin. The areas with increasing slopes for the 99th percentiles of wind speed coincide with the areas with increasing slopes for the number of extreme wind events occurrences.  Data Availability Statement: ERA5 reanalysis data on single levels from 1979 to the present can be freely accessed and downloaded at: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form (accessed on 13 July 2021).