Estimation of the Near Future Wind Power Potential in the Black Sea

The main objective of the present study is to quantify the recent past and explore the near future wind power potential in the Black Sea basin, evaluating the possible changes. Furthermore, an analysis of the wind climate in the target area was also performed. The wind resources have been assessed using the wind fields provided by various databases. Thus, the wind power potential from the recent past was assessed based two different sources covering each one the 30-year period (1981–2010). The first source is the ERA-Interim atmospheric reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF), while the second source represents the hindcast wind fields simulated by a Regional Climate Model (RCM) and provided by EURO-CORDEX databases. The estimation of the near future wind power potential was made based on wind fields simulated by the same RCM under future climate projections, considering two Representative Concentration Pathways (RCPs) scenarios (RCP4.5 and RCP8.5) and they cover also a 30-year time interval (2021–2050). Information in various reference points were analyzed in detail. Several conclusions resulted from the present work. Thus, as regards the mean wind power potential in winter season, in 51% of the locations a significant increase is projected in the near future (both scenarios). Besides providing a detailed description of the wind conditions from the recent past over the Black Sea basin considering two major sources, the novelty of the present work consists in the fact that it gives an estimation of the expected wind climate in the target area for the near future period and at the same time an evaluation of the climate change impacts on the wind speed and wind power potential.


Introduction
The global context generated by the rapid population growth and by the development of new technologies implies an increase in the energy needs.By combining this aspect with the actual energy producing capability, which uses mainly conventional methods (natural gas, lignite, petroleum, etc.) and has a negative impact to the environment, it results that it is urgently required to find effective methods of counteraction.From this perspective, the topic of energy remains an extremely important one [1].Globally, many policies aiming to put into practice, on medium and long-term, methods for the energy sector development were adopted.In this way, the energy sector will become more efficient and sustainable.Developing and implementing energy policies is also an important condition in order to achieve strategic goals, which will propagate further in the economic development.
Various studies show that the feed-in tariff policy mechanism is really effective in fostering the sustainability transition of the energy sector and also to promote the investments in renewable energy [2,3].Thus, the growth of renewable resource extraction (wind, wave, geothermal, photovoltaic, and hydro, etc.) for generating electricity is pursued in the detriment of the conventional methods.basin is a suitable environment for the green energy extraction (average wind power being estimated in the range 500-900 W/m 2 , with a deviation of ±50 W/m 2 during the future periods).
Onat et al. [25] conducted an analysis of the wind climate and of the wind energy potential for several regions in Turkey.The authors analyzed also a small region of the Black Sea, located in the west (Amasra) by using a five-layer Sugeno-type ANFIS model developed with MATLAB-Simulink software.The relationship between the wind speed and other climate variables was also determined and the resulted data confirm that the Amasra region is a location with a good potential for the wind energy extraction.According to this study, the average power density at 10 m height is 232 W/m 2 (nearly good), at 50 m height is 603 W/m 2 (good) and at 80 m height is 1300 W/m 2 (very good).
The wind pattern of the Black Sea was also evaluated by Onea and Rusu [17], considering 12 years of data from the U.S. National Centers for Environmental Prediction (NCEP).In that study, the authors analyzed the wind power distribution taking into account the diurnal versus nocturnal variations.According to the above mentioned study, the northwestern and northeastern sectors of the Black Sea are the most suitable for wind energy extraction.The northwestern sector of the Black Sea was also analyzed by Lin-Ye et al. [26].Their approach considered a hybrid methodology involving the Simulating WAve Nearshore (SWAN) spectral wave-model to produce wave-climate projections.The wave model was forced with wind-fields corresponding to the two climate change scenarios.
From this perspective, the present study aims to characterize the offshore wind power potential of the Black Sea during the present and the near future, by analyzing four databases.Following this objective, the structure of the proposed work includes first a presentation of the materials and methods considered, focused on the description of the target area and of the databases taken into account.The next section presents the results providing in some reference points the wind speed and also the wind power.These relate both the 30-year period considered from the past  and that estimated for the near future period analyzed (2021-2050).Finally, a discussion of the results was also carried out.Thus, the novelty of this paper can be summarized as follows: 1.
A detailed description of the wind conditions over the Black Sea basin from two major sources (Era-Interim and Euro-Cordex) covering the recent past (1981-2010).

3.
Evaluation of the climate change impacts on the wind speed and wind energy potential by performing comparisons between the past and the future projections.

Target Area
In this study, the target area is the Black Sea.The Black Sea basin is under the influence of the NAO (North Atlantic Oscillations) mechanism, which causes, due to the influence of the cold air arriving from the northern regions, significant storm events during the winter [27].
Figure 1 illustrates the Black Sea basin and the geographical locations of the reference points further considered.In this study, the wind speed at 100 m height was evaluated in twenty-four points.The reference points are divided into two different categories.The first category includes the shallow water locations.More precisely, these points are in the range 29-39 m water depth (average depth being 34.9 m) and the distance to shore in the range 1.5 to 56 km (average distance to shore being 11.8 km).The second category refers to the deep water locations.These points are relatively close to shore in the range 3.6 to 109 km (the average distance to shore being 35.7 km) and a depth in the range 114 to 140 m, with an average value for the depth of 125.8 m.
Taking into account some previous studies that showed the existence of various areas in the Black Sea presenting different wind conditions, the target area was divided into five geographical zones with similar characteristics/patterns, labeled from A to E. Zone A contains seven points related to the coastlines of Romania and Ukraine (shallow water: A.n.1 to A.n.  1 and 2).

ECMWF Dataset
One set of the data considered was obtained from the ERA-Interim database.This is a reanalysis project conducted by the European Centre for Medium-Range Weather Forecasts (ECMWF) [28].ECMWF uses forecast models and data assimilation techniques that include a 4D analysis that has 12 h analysis windows to describe the atmosphere and oceans [29,30].
ERA-Interim is an ongoing project that comprises datasets with a various number of marine and atmospheric parameters from 1979 to the present.Among many data, it contains information about the wind speed components (u-zonal velocity and v-meridional velocity) at 10 m height.The wind speed at 10 m will be denoted as U 10 .The data considered have a spatial resolution of 0.75 • × 0.75 • , cover a period of 30-years (1 January 1981 to 31 December 2010) and are dived in four hourly intervals (corresponding to 00:00:00, 06:00:00, 12:00:00, 18:00:00) for each day.

EURO-CORDEX Database
The Coordinated Regional Downscaling Experiment (CORDEX) initiated by the World Climate Research Program (WCRP), produces high-resolution 'downscaled' climate data through a global partnerships.EURO-CORDEX is the European branch of this project, where results from several RCMs were jointed to cover the European continent [31,32].The data available have the spatial resolutions of 0.11 • and 0.4 • , respectively.
The EURO-CORDEX wind fields used in this study have a maximum range of 27 • N ÷ 72 • N and 22 • W ÷ 45 • E, with a spatial resolution of 0.11 • .These data contain information about the wind speed components at 10 m height and are dived in four hourly intervals (corresponding to 00:00:00, 06:00:00, 12:00:00, 18:00:00) for each day.From this database, three types of wind fields were used in the present study.First, the evaluation (hindcast) wind fields that cover a period of 30-years ( The second and third datasets cover the 30-year interval 1 January 2021 to 31 December 2050, and there are wind fields simulated by the RCA4 model under future climate projections, considering two Representative Concentration Pathways (RCPs) emission scenarios, RCP4.5 and RCP8.5.More detailed information regarding the CORDEX scenarios for the European areas, as provided by the Rossby Centre regional climate model RCA4, are given in [33,34].These two scenarios were simulated with the same RCM model as in the case of evaluation data, but forced by a Global Climate Model (GCM), namely EC-EARTH.In the global climate models, the related temporal evolution of atmospheric greenhouse gas and aerosol concentrations are prescribed, which then simulate the response of the climate system to the forcing.Thus, a range of potential future climate evolutions can be projected.

Data Evaluation
First, the data used in this study were evaluated in terms of daily, seasonal and yearly average.Both ERA-Interim and EURO-CORDEX databases provide the wind speed at 10 m height in terms of its components.For all the sites available and for all time scales various statistical analyses of the wind speed are conducted.Also, for each time series of the wind speeds the 5th and 95th percentiles were computed.The statistical parameters considered are the root mean square error, bias and the Pearson correlation coefficient.The root mean square error and bias are computed with the following relationships: where n is the total number of data pairs, x is the wind speed provided by ERA-Interim, while y is the evaluation wind speed value.
Maintaining the same notation of the data, the Pearson correlation coefficient is computed with the next equation: where x and y are the mean values of the analyzed datasets.
The linear trend of the wind speeds was also estimated from a linear regression whose y dependent variable is represented by the wind speed time series, while the independent one is the time (x variable).The linear regression line has the following equation: where b is the slope of this line, and a is the intercept.The slope indicates the linear rate of wind speed change and its values can be positive or negative.A positive value corresponds to increasing trends while the negative value indicates a decreasing trend.
Taking into account that consistent wind speed measurements over the Black Sea basin are not available, a comparison between the present wind climate simulated by the RCM model with that resulted from a widely used reanalysis database is performed.In this way, the skill of the RCM wind fields to represent the present wind climate in the Black Sea basin can be determined.
Thus, the evaluation (hindcast) wind fields from EURO-CORDEX are compared with the wind fields from ERA-Interim reanalysis.The evaluation wind fields being provided by a RCM model, only statistical comparisons can be performed [35,36].Near the Romanian coast, the wind speeds are recorded at an offshore platform, and the evaluation wind speeds were compared with these measurements by Rusu et al. [37].

Wind Speed Analysis
This section provides a detailed analysis of the wind speed at 10 m height.First, comparisons between the present wind climate simulated by the RCM model with that resulted from a widely used reanalysis databases ERA-Interim were performed to evaluate the skill of the RCM model.Also, some comparisons between the wind fields simulated under RCP4.5 and RCP8.5 scenarios were performed in order to identify the differences induced by these two different scenarios.Actually, these are the most studied greenhouse gas concentration scenarios.Thus, RCP4.5 describes an intermediate concentration scenario with radiative forcing stabilized at around 4.5 W/m 2 , while RCP8.5 describes a high concentration scenario under which the radiative forcing is expected to be higher than 8.5 W/m 2 by the end of the year 2100 [38,39].On the other hand, the comparisons between the evaluation data and those obtained under various scenarios are indicating the possible further evolution of the wind conditions.

Hindcast Data
Figure 2 illustrates a comparative analysis in terms of maximum values, averages, 5th and 95th percentiles (lower than 5% and 95%) for the total time interval.In Figure 2a the values computed in each point for the present period are compared, considering the evaluation and ERA-Interim wind speeds.Some differences are observed between the magnitudes of the wind speeds.With the exception of the point B.n.3, located near the Bosporus Strait, the EURO-CORDEX Evaluation wind speeds are higher than those from ERA-Interim.The resulted values show that the difference is in average about 17% regarding the maximum computed values.At the same time, significant differences can be noticed regarding the total interval averages (12%) and for 95th (lower 95%) percentile analysis (14%).
Energies 2018, 11, x 7 of 21 Thus, the evaluation (hindcast) wind fields from EURO-CORDEX are compared with the wind fields from ERA-Interim reanalysis.The evaluation wind fields being provided by a RCM model, only statistical comparisons can be performed [35,36].Near the Romanian coast, the wind speeds are recorded at an offshore platform, and the evaluation wind speeds were compared with these measurements by Rusu et al. [37].

Wind Speed Analysis
This section provides a detailed analysis of the wind speed at 10 m height.First, comparisons between the present wind climate simulated by the RCM model with that resulted from a widely used reanalysis databases ERA-Interim were performed to evaluate the skill of the RCM model.Also, some comparisons between the wind fields simulated under RCP4.5 and RCP8.5 scenarios were performed in order to identify the differences induced by these two different scenarios.Actually, these are the most studied greenhouse gas concentration scenarios.Thus, RCP4.5 describes an intermediate concentration scenario with radiative forcing stabilized at around 4.5 W/m 2 , while RCP8.5 describes a high concentration scenario under which the radiative forcing is expected to be higher than 8.5 W/m 2 by the end of the year 2100 [38,39].On the other hand, the comparisons between the evaluation data and those obtained under various scenarios are indicating the possible further evolution of the wind conditions.

Hindcast Data
Figure 2 illustrates a comparative analysis in terms of maximum values, averages, 5th and 95th percentiles (lower than 5% and 95%) for the total time interval.In Figure 2a the values computed in each point for the present period are compared, considering the evaluation and ERA-Interim wind speeds.Some differences are observed between the magnitudes of the wind speeds.With the exception of the point B.n.3, located near the Bosporus Strait, the EURO-CORDEX Evaluation wind speeds are higher than those from ERA-Interim.The resulted values show that the difference is in average about 17% regarding the maximum computed values.At the same time, significant differences can be noticed regarding the total interval averages (12%) and for 95th (lower 95%) percentile analysis (14%).The analysis of the annual average values illustrated in Figures 3-6 (see also Figures A1 and A2) shows the fact that wind speed at 10 m height had and it is expected to have a different evolution in The analysis of the annual average values illustrated in Figures 3-6 (see also Figures A1 and A2) shows the fact that wind speed at 10 m height had and it is expected to have a different evolution in the areas considered.For instance, the data resulted in zone A (shallow water data, Figure 3) show that in the time interval 1981 to 2010, the wind speed has values varying in the range 4 to 7 m/s.It can be noticed that the evaluation wind speed average values from RCM are higher than those from ERA-Interim.The difference is in a range of about 9% (point A.n.1) to 26% (point A.n.4).
Figure 4 illustrates two of the most representative shallow water points, close to the coastline of Bulgaria.By analyzing the data resulted for this area, it can be observed that in two points, the evaluation data are higher than ERA-Interim (point B.n.1: ∼ = 26% and point B.n.2, not shown here ∼ =21%).In the case of the point B.n.3, the Evaluation data are lower with ∼ = 16%.Regarding the point B.n.4, not shown here, both models (ECMWF ERA-Interim and EURO-CORDEX Evaluation) present similar behavior of the wind speed annual pattern evolution.According to the linear tendency, the annual means of the wind speeds show a different evolution for the period 1981 to 2010.As in zone A, in this zone the near future scenarios assess differently the wind speed pattern.Comparing the data from 1981-2010 with the near future interval for point B.n.1 it can be noticed a substantial increase of the wind speed.It is expected that the mean wind will reach values in the range of 5.7-6.6 m/s.
The wind speed analysis of two of the most representative deep water points located in zones A and B (A.o.1) and B (B.o.1) is presented in Figure 5. Regarding the deep water locations for zone A, the Evaluation wind speeds for the interval 1981-2010 are higher with approximately 9-12% than those from ERA-Interim.If zone B is analyzed, the difference between both data grows at about 15% to 29%.The data corresponding to the points associated with the coastline of Turkey are presented in Figure 6.By analyzing the data for the interval 1981-2010, in comparison with the data previously presented, it can be noticed that the wind is less intense in this area.Figure 6           The seasonal distributions of the wind intensity (total, winter, spring, summer, and autumn) are presented in Figures 7 and 8.By looking at these data, it can be noticed that during the winter the magnitude of the wind speed average is higher.In almost 90% of the cases presented, the wind speeds during the summer interval have the lowest values.For example, the points C.n.1 and C.o.1, which are located north of Turkey, have during summer wind speed averages almost as high as in the winter interval.However, it can be noticed that in these regions the difference between the seasonal averages is small.High wind speed values are encountered in the nearshore and offshore of the Crimea Peninsula.In this area, the wind speed at 10 m height was assessed for three locations, but here only the most representative results are presented.Thus, Figure A2 illustrates the annual averages of the wind speeds at 10 m height of two points located in shallow water (E.n.2) and deep water (E.o.1) of the Crimea Peninsula.By analyzing the time interval 1981-2010, it can be noticed that the wind has an uptrend evolution.Small differences are observed between both data.In the case of the point E.n.1, the Evaluation wind is higher with about 2.8%, for point E.n.2 with 7.1%, while for point E.o.1 with 7.3%.
The seasonal distributions of the wind intensity (total, winter, spring, summer, and autumn) are presented in Figures 7 and 8.By looking at these data, it can be noticed that during the winter the magnitude of the wind speed average is higher.In almost 90% of the cases presented, the wind speeds during the summer interval have the lowest values.For example, the points C.n.1 and C.o.1, which are located north of Turkey, have during summer wind speed averages almost as high as in the winter interval.However, it can be noticed that in these regions the difference between the seasonal averages is small.The seasonal distributions of the wind intensity (total, winter, spring, summer, and autumn) are presented in Figures 7 and 8.By looking at these data, it can be noticed that during the winter the magnitude of the wind speed average is higher.In almost 90% of the cases presented, the wind speeds during the summer interval have the lowest values.For example, the points C.n.1 and C.o.1, which are located north of Turkey, have during summer wind speed averages almost as high as in the winter interval.However, it can be noticed that in these regions the difference between the seasonal averages is small.The seasonal distributions of the wind intensity (total, winter, spring, summer, and autumn) are presented in Figures 7 and 8.By looking at these data, it can be noticed that during the winter the magnitude of the wind speed average is higher.In almost 90% of the cases presented, the wind speeds during the summer interval have the lowest values.For example, the points C.n.1 and C.o.1, which are located north of Turkey, have during summer wind speed averages almost as high as in the winter interval.However, it can be noticed that in these regions the difference between the seasonal averages is small.The correlation level between the data provided by the ECWMF ERA-Interim databases and the EURO-CORDEX Evaluation climate model results is presented in Table 3 through the Pearson correlation coefficient.This analysis was performed considering the daily averages.By analyzing the data, it can be observed that the Pearson coefficient varies in the range 0.51-0.92.In order to evaluate the differences between the daily values, the root mean square error (RMSE) and the mean error (Bias) were also computed.The data presented in Table 3 show that RMSE is in the range 0.25-0.5 m/s.As regards the bias values, in all points they are negative (ranging from −0.25 to −0.06 m/s) indicating that the mean Evaluation wind speeds are higher than those computed for ERA-Interim.
According to the correlation levels defined in Table 4, it can be noticed that the correlation between both data ranges from a reasonable correlation to a very high correlation.More precisely, 71% of the data are in the interval of very high correlation, 25% are in the interval of high correlation and only 4% indicate a reasonable correlation.Thus, the lower RMSE and Bias values, together with the higher values of the correlation indexes, show that between the data provided by ERA-Interim and those from the RCM model there is a good agreement.Based on the annual averages of the wind speeds, the linear trend in each point was computed and presented in Table 5.The ERA-Interim data shows a slight increasing (ranging from 0 to 0.132 m/s per decade) or decreasing (ranging from −0.054 to −0.002 m/s per decade) trends.The decreasing trends are found only in the western part of the Black Sea basin.These values of the linear trend are in line with those computed by Torralba et al. [40].As regards the Evaluation data, in 23 points the linear trend has positive values from 0.004 to 0.119 m/s per decade.Another comparative analysis between the wind speeds projected by two of the future scenarios available in the framework of the EURO−CORDEX project, RCP4.5 and RCP8.5, is presented in Figure 2b.By analyzing Figure 2b it can be noticed that between both data no significant differences appear in terms of average values, 5th and 95th percentiles.On the other hand, major differences in the evaluation of the maximum wind speed values in the shallow water areas of Bulgaria and Romania are observed.
According to the near future data, both scenarios present approximately the same values in terms of annual averages, but with a lag of several years between the peaks.As previously presented, a slight difference occurs between the data from RCP4.5 and RCP8.5.In the case of zone A (deep water locations), the data for RCP4.5 scenarios are on average with 1% higher, while for zone B with almost 1% smaller than the 8.5 data.By conducting a comparative analysis of the shallow water and deep water locations for both zones A and B it can be concluded that there is no notable difference in terms of wind intensity.
The near future projections (time interval 2021-2050) near to the Turkish coast (see Figure 6) present similar average values of the wind speeds as those computed from the EURO-CORDEX Evaluation data for the present climate.The near future annual wind speed analysis does not highlight any trend.Moreover, the differences between the annual averages are lower than the previous cases.
Although the data for the interval 1981-2010 show an increase of the wind speed near to the Russian coastline (see Figure A1), both estimates an approximate constant pattern for the 30-year interval of the near future.Both RCP4.5 and RCP8.5 scenarios predict for the Crimea Peninsula region an increase of the wind at 10-m height (Figure A2).
Regarding the linear tendency, for the 2021-2050 time period, the RCP4.5 scenario shows that the wind speed will have a low decrease while the RCP8.5 a small increase.The values of the linear trends computed for near future data are presented in Table 6.The four datasets analyzed contain information about the zonal and meridional velocities needed for analyzing the wind main directions from which the winds are blowing.From the wind roses presented in Figure 9, it can be observed the dominant direction and also the predominant wind speed range for all five zones studied.At a first sight, the data show that over the Black Sea there are two main wind patterns.The first pattern is more spread and it is observed in the northwest and west of the Black Sea basin.According to all the models considered, the main directions for these zones are north, northeast and southwest.However, it can be noticed that all three data from EURO-CORDEX show also that the wind blows quite often from the west.Moreover, the most frequent values of the wind speeds are in the interval 0 to 12 m/s.

Wind Power Analysis
The near future projections of the wind power potential of the Black Sea for the time interval 2021-2050 were also computed and compared with those corresponding to the time interval 1981-2010.Nowadays, the typical hub heights for the offshore wind turbines are ranging about 100 m [41], ).More precisely, in these points, the averages computed for the 30-year period are higher than 500 W/m 2 .Regarding zone A, the higher increase of the average wind power in the near future is observed in point A.n.4 (Figure 10b), where this will pass the level of 400 W/m 2 , while in the present it is below this threshold.Regarding the zone B, two opposite cases are noticed.The averages for the points B.n.1 and B.o.1 present a significant increase.Thus, the wind power average for B.n.1 (Figure 10d) exceeds the value of 400 W/m 2 , while in the case of the point B.o.1 (Figure 10f) this will be higher than 500 W/m 2 .As above mentioned, an opposite behavior is found in the point B.n.3 (Figure 10e), which is located near Bosporus Strait.In the case of zone C, points C.n.1 (Figure 10g) and C.o.1 (Figure 10h), no major changes seem to be encountered in terms of the wind power averages.They will be maintained near the threshold of 200 W/m 2 .As regards the reference points from the zone D, the higher increase in the average wind power is observed at the point located in deep water, D.o.2 (Figure A3b).This point represents in fact one of the most energetic locations in the Black Sea basin.
Another zone that can have a significant impact on the wind energy exploitation is zone E, situated near the Crimea Peninsula.Here, the data in the points E.n.2 (Figure A3c) and E.o.1 (Figure A3d) were studied in more detail.By analyzing these data, it can be observed that even though these two points are relatively close, the growth of the average value is higher in the deep water point Regarding the zone B, two opposite cases are noticed.The averages for the points B.n.1 and B.o.1 present a significant increase.Thus, the wind power average for B.n.1 (Figure 10d) exceeds the value of 400 W/m 2 , while in the case of the point B.o.1 (Figure 10f) this will be higher than 500 W/m 2 .As above mentioned, an opposite behavior is found in the point B.n.3 (Figure 10e), which is located near Bosporus Strait.In the case of zone C, points C.n.1 (Figure 10g) and C.o.1 (Figure 10h), no major changes seem to be encountered in terms of the wind power averages.They will be maintained near the threshold of 200 W/m 2 .As regards the reference points from the zone D, the higher increase in the average wind power is observed at the point located in deep water, D.o.2 (Figure A3b).This point represents in fact one of the most energetic locations in the Black Sea basin.
Another zone that can have a significant impact on the wind energy exploitation is zone E, situated near the Crimea Peninsula.Here, the data in the points E.n.2 (Figure A3c) and E.o.1 (Figure A3d) were studied in more detail.By analyzing these data, it can be observed that even though these two points are relatively close, the growth of the average value is higher in the deep water point E.o.1 (about 100 W/m 2 ), than in the shallow water point E.n.2.
The wind power averages at 100 m computed for the entire 30-year period and corresponding to the 12 reference points are now compared with the seasonal values, the results being presented in Figure 11.In this analysis, the four datasets available are also considered.The results clearly show that the winter averages are the highest, while the least energetic season is the summer.The other two seasons present almost equal values of the wind power averages.

Conclusions
From the comparative analysis of the wind speed at 10 m over the sea level simulated by RCM from EURO-CORDEX, the RCA4 model respectively, against the data existent in ERA-Interim database, a good correlation was noticed between these two datasets in the Black Sea area.On the other hand, probably due to the fact that the spatial resolution of the RCM model is higher than that from ERA-Interim, the wind speeds provided by RCA4 in the reference points considered in the Black Sea basin are slightly higher than those indicated by ERA-Interim.
The analysis of the wind speed annual averages shows that the wind speeds at 10 m height have a different evolution in the areas considered in the Black Sea.In general, for the present period, it was observed that the Evaluation wind speed averages are higher than those from ERA-Interim with values ranging from 0.06 to 0.25 m/s.Through the linear trend, it was observed that in 75% of the cases both data have the same trend (upward trend-the maximum value is 0.144 m/s per decade).
The daily averages of the wind speeds provided by both databases present a very good correlation.The Pearson correlation coefficient varies in the range 0.51 to 0.92.More precisely 71% of the data are in the interval of very high correlation, 25% are in the interval of high correlation and only 4% have a reasonable correlation.The differences between the daily values were assessed by the analysis of the root mean square error and Bias.The results of the analysis show that RMSE is in the range 0.25 to 0.50 m/s, while the bias values, in all points they are negative (ranging from −0.25 to

Conclusions
From the comparative analysis of the wind speed at 10 m over the sea level simulated by RCM from EURO-CORDEX, the RCA4 model respectively, against the data existent in ERA-Interim database, a good correlation was noticed between these two datasets in the Black Sea area.On the other hand, probably due to the fact that the spatial resolution of the RCM model is higher than that from ERA-Interim, the wind speeds provided by RCA4 in the reference points considered in the Black Sea basin are slightly higher than those indicated by ERA-Interim.
The analysis of the wind speed annual averages shows that the wind speeds at 10 m height have a different evolution in the areas considered in the Black Sea.In general, for the present period, it was observed that the Evaluation wind speed averages are higher than those from ERA-Interim with values ranging from 0.06 to 0.25 m/s.Through the linear trend, it was observed that in 75% of the cases both data have the same trend (upward trend-the maximum value is 0.144 m/s per decade).
The daily averages of the wind speeds provided by both databases present a very good correlation.The Pearson correlation coefficient varies in the range 0.51 to 0.92.More precisely 71% of the data are in the interval of very high correlation, 25% are in the interval of high correlation and only 4% have a reasonable correlation.The differences between the daily values were assessed by the analysis of the root mean square error and Bias.The results of the analysis show that RMSE is in the range 0.25 to 0.50 m/s, while the bias values, in all points they are negative (ranging from −0.25 to −0.06 m/s).Thus, the low RMSE and Bias values together with the higher values of the correlation indexes show that there is a good agreement between the data provided by ERA-Interim and those from the RCM model (Evaluation data).
As the results of the present work show, in most of the cases there are no relevant differences in the average wind speeds (smaller than 0.4 m/s) simulated under both scenarios.This is probably due to the fact that until the mid-century the differences between the two RCPs considered are not very high.According to the near future data, both scenarios present approximately the same values in terms of annual averages.Regarding the linear tendency for the period 2021-2050, the RCP4.5 scenario shows that the wind speed will have a low decrease in 92% of the points (values ranging from −0.044 to −0.009 m/s per decade) while the RCP8.5 a small increase in 83% of the cases (values ranging from 0.008 to 0.044 m/s per decade).
As regards the average wind power potential in winter season, for 51% of the reference points a significant increase was observed for the near future (both scenarios) compared with the present values, while for 41% only a slight growth.For both scenarios, a small decrease was noticed only in the point B.n.1.The most energetic zones of the Black Sea are the western part of the basin (northwest and southwest areas) and also in the east and south of the Crimea Peninsula.More precisely, in these points, the total wind power averages for the 30-year time interval are higher than 500 W/m 2 .As expected, in the winter time the higher wind speeds and wind powers are encountered, while the lower values are found in the summer.
The fact that it was noticed an increase of the mean values of the wind resources in 95.6% of the reference points considered in the Black Sea, either under the RCP4.5 or RCP8.5 scenarios, can be considered beneficial from the perspective of the wind projects and this can give momentum to the installation of the wind farms in the areas already identified as having a good potential, as for example the western side of the basin.From this perspective, it is well known that some European areas are affected by a decrease of the wind resources (see for example [35,47,48]).
4 and deep water: A.o.1 to A.o.3), B.o.2) related to the coasts of Bulgaria's and the northwest of Turkey.Zone C contains two shallow water points C.n.1 and C.n.2 and only one deep water point C.o.1 located in the north and northeast of Turkey, zone D contains five points associated to the coastal environment of Georgia and Russia (shallow water: D.n.1 to D.n.3, deep water points: D.o.1 and D.o.2) and zone E with three points close to the Crimea Peninsula (shallow water: E.n.1 and E.n.2, deep water: E.o.1) (see Tables

Figure 1 .
Figure 1.The geographical locations of the reference points corresponding to the 5 zones considered: zone A (shallow water points: A.n.1 to A.n.4 and deep water points: A.o.1 to A.o.3), zone B (shallow water points: B.n.1 to B.n.4 and deep water points: B.o.1 and B.o.2), zone C (shallow water points C.n.1 and C.n.2 and deep water point C.o.1), zone D (shallow water points: D.n.1 to D.n.3, deep water points: D.o.1 and D.o.2) and zone E (shallow water points: E.n.1 and E.n.2, deep water point: E.o.1); (a) Overview of the target zones; (b).The location of the points corresponding to zones A, E and a section of zone D; (c) The location of the points corresponding to zones B, C and a section of zone D.

Figure 1 .
Figure 1.The geographical locations of the reference points corresponding to the 5 zones considered: zone A (shallow water points: A.n.1 to A.n.4 and deep water points: A.o.1 to A.o.3), zone B (shallow water points: B.n.1 to B.n.4 and deep water points: B.o.1 and B.o.2), zone C (shallow water points C.n.1 and C.n.2 and deep water point C.o.1), zone D (shallow water points: D.n.1 to D.n.3, deep water points: D.o.1 and D.o.2) and zone E (shallow water points: E.n.1 and E.n.2, deep water point: E.o.1); (a) Overview of the target zones; (b).The location of the points corresponding to zones A, E and a section of zone D; (c) The location of the points corresponding to zones B, C and a section of zone D.

21 Figure 3 .
Figure 3. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Romanian nearshore of.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 3 .
Figure 3. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Romanian nearshore of.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 3 .
Figure 3. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Romanian nearshore of.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 4 .
Figure 4. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Bulgarian nearshore.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 5 .
Figure 5. Annual averages of the wind speed at 10 m height annual averages for two of the most representative points located in deep water Bulgaria (B.o.1) and Romania (A.o.1).ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 4 .
Figure 4. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Bulgarian nearshore.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 3 .
Figure 3. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Romanian nearshore of.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 4 .
Figure 4. Annual averages of the wind speed at 10 m height for two of the most representative points located in the Bulgarian nearshore.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 5 .
Figure 5. Annual averages of the wind speed at 10 m height annual averages for two of the most representative points located in deep water Bulgaria (B.o.1) and Romania (A.o.1).ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 6 .
Figure 6.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (C.n.1) and in deep water (C.o.1) of Turkey.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 6 .
Figure 6.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (C.n.1) and in deep water (C.o.1) of Turkey.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 7 .Figure 8 .
Figure 7. Seasonal wind speed averages at 10 m height compared with the total time averages for some representative points located in shallow water and deep water of Bulgaria and Romania coasts; (a) ERA.Int and Eval data, time interval 1980-2010; (b) RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 7 .
Figure 7. Seasonal wind speed averages at 10 m height compared with the total time averages for some representative points located in shallow water and deep water of Bulgaria and Romania coasts; (a) ERA.Int and Eval data, time interval 1980-2010; (b) RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 6 .
Figure 6.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (C.n.1) and in deep water (C.o.1) of Turkey.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 7 .Figure 8 .
Figure 7. Seasonal wind speed averages at 10 m height compared with the total time averages for some representative points located in shallow water and deep water of Bulgaria and Romania coasts; (a) ERA.Int and Eval data, time interval 1980-2010; (b) RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure 8 .
Figure 8. Seasonal wind speed averages at 10 m height compared with the total time averages for some representative points located in shallow water and deep water of Crimea Peninsula, Russia and Turkey coasts; (a) ERA.Int and Eval data, time interval 1981-2010; (b) RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure A1 .
Figure A1.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (D.n.3) and in deep water (D.o.2) of Russia.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure A2 .
Figure A2.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (E.n.2) and in deep water (E.o.1) of Crimea Peninsula.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.

Figure A2 .Figure A3 .
Figure A2.Annual averages of the wind speed at 10 m height for two of the most representative points located in shallow water (E.n.2) and in deep water (E.o.1) of Crimea Peninsula.ERA.Int and Eval data, time interval 1981-2010, RCP4.5 and RCP8.5 data, time interval 2021-2050.Energies 2018, 11, x 19 of 21

Table 1 .
The geographical locations and the characteristics of the shallow water points.

Table 2 .
The geographical locations and the characteristics of the deep water points.
1 January December 2010) are analyzed.The wind fields are simulated by the Rossby Centre regional climate model-RCA4 model at the Swedish Meteorological and Hydrologic Institute (SMHI) and forced with initial and lateral boundary conditions provided by ECMWF ERA-Interim reanalysis data.

Table 3 .
Statistical evaluation in terms of the correlation coefficient (r) and RMSE for all 24 points.

Table 4 .
Levels of interpretation of the Pearson correlation coefficient (r).

Table 5 .
Linear trend (m/s per decade) values computed for all 24 points for recent past data.

Table 6 .
Linear trend (m/s per decade) values computed for all 24 points for near future data.