Aalborg Universitet Spatial Distribution of the Baltic Sea Near-Shore Wave Power Potential along the Coast of Klaipda

Wave power is an abundant source of energy that can be utilized to produce electricity. Therefore, assessments of wave power resources are being carried out worldwide. An overview of the recent assessments is presented in this paper, revealing the global distribution of these resources. Additionally, a study, which aims to assess the spatial distribution of the Baltic Sea near-shore wave power potential along the coast of Klaipėda (Lithuania), is introduced in this paper. The impacts of the wave propagation direction and decreasing depth on wave power resources were examined using the numerical wind-wave model MIKE 21 NSW. The wave height loss of the design waves propagating to shore was modelled, and the wave power fluxes in the studied depths were calculated using the JONSWAP wave spectrum modified for the Baltic Sea. The results revealed that all waves that propagate to the shore in the Baltic Sea near-shore area along the coast of Klaipėda from 30 m depth to 5 m depth lose at least 30% of their power. Still, most common waves in this area are low, and therefore, they start to lose their power while propagating to the shore at relatively low (10–14 m) depths. To turn this into an advantage the wave power converter would have to work efficiently under low power conditions.


Introduction
It is globally recognized that wave power has great potential to be a future source of electricity.More important is that this is a renewable and abundant resource.The attractiveness of wave power potential is reflected in the growth of the Ocean Energy Systems (IEA-OES) organization.This organization was founded in 2001 by Portugal, Denmark and UK representatives.By 2013, when Nigeria and Monaco joined IEA-OES, the organization had 21 state members [1].The increasing attention that the power of waves has attracted in the last century has resulted in a vast number of scientific studies.A fair amount of these studies are dedicated to the assessments of wave power potentials.
Global wave power resources are usually presented as within the interval from 1 to 10 TW [2].These values are quoted so frequently that it is difficult to trace the original source.Kinsman, the author of these values, has warned about the potential errors in his calculations [3].The latest attempts to revise the value of global theoretical wave power resources were published in 2010 and 2012, and the Mexico.Therefore, there is suggestion to supplement needed energy with electricity generated from the waves [14].
Uruguay and Brazil are two other countries in South America where attempts to assess the wave energy potential have been made.Studies were carried out in territories within close proximity to each other-the Rio la Plata and Uruguayan waters in the Atlantic Ocean [15] and the coast of the state of Santa Catarina in Brazil [16].Although different wave models were selected (WaveWatch III and MIKE 21 SW) for the near-shore wave energy assessment, the authors of both papers remarkably elected to analyse wave energy resources at 20 m depth.The results obtained were very similar: along the Uruguayan coast, the wave power flux ranges from 8 to 14 kW/m, while along the Brazilian coast, the flux ranges from 8 to 14.5 kW/m.Both studies highlight the low seasonal variability of the wave energy, which is what makes these coasts particularly attractive for electricity production from waves [15,16].

Asia
There are plenty of recently published studies concerning wave energy resource assessments in Asia.Firstly, we concentrate on the eastern part of the Mediterranean Sea.Here, Lebanon has joined a not inconsiderable number of European countries that are interested in electricity production from Mediterranean Sea waves.Although the average wave height along the coast of Beirut is 78.1 cm, the noted average annual wave power flux is 4.6 kW/m [17].Interestingly, the authors of this study are trying to evaluate the performance of Pelamis, WaveDragon and AquaBuoy wave energy converters (WECs) in these conditions.Unsurprisingly, the capacity factor for these converters on the eastern coast of the Mediterranean Sea is 4-6%, and the authors come to conclusion that these WECs are not optimum solutions in these waters [17].On the other hand, the scarcity of wave energy converter power matrixes is limiting the options for scientists in this type of study.Therefore, scientists in Iran, whose goal was to select the most suitable wave energy conversion technology for the Caspian Sea, chose a different approach.With knowledge of the wave power flux range in the Caspian Sea (5 to 14 kW/m) and specific features of the Caspian Sea (salinity, water level fluctuations, near-shore seabed composition, etc.), they selected the most suitable of over 20 state of the art WECs by assigning weights according to each WEC design parameter [18].Of course, this rating of WECs is quite subjective, but it is a methodically suitable solution.
Wave energy resources were assessed in detail for northern Iran, the southwestern region of the Caspian Sea [19].The study area was near the port of Anzali, and bathymetry data up to 700 m was used.With the help of SWAN for the wave climate model, it was found that an average annual wave power varies from 0.2 kW/m to 1.2 kW/m.Also, there is quite a high seasonal variation of the wave power in this region: wave power in winter may be up to 2.5 times higher than in summer [19].The findings of this study also imply that the authors of the previous reference may have overestimated the power potential of the Caspian Sea waves.
The interest in wave power in Iran does not end in the Caspian Sea: the wave power potential of the Persian Gulf was also analysed [20].Again, with the help of the SWAN model, not only the wave power potential but also the influence of the prevailing wind direction on the wave power in this shallow gulf (average depth 36 m) was examined.The findings were remarkable: a 6% change of the wind direction in this case can result in a 77% decrease of the wave power.The estimated average wave power in the Persian Gulf is not high: in the deepest areas, it reaches 2 kW/m [20].
Attention to wave energy is growing in China.This is proved by the number of recent original research [21,22] and review [23,24] papers.Unfortunately, the highest average value of the wave power flux along the Chinese shore in the northern part of the South China Sea is only 16 kW/m [21].The near-shore numerical model of the Shandong peninsula's coasts revealed that the average wave power flux of the Yellow Sea is only 5.1 kW/m [22].Wave power flux values are similar for the South Korea coasts, and although it is noted that in the winter the wave power flux here can reach Energies 2017, 10, 2170 4 of 18 25 kW/m [25], this value is more important for the wave energy converter load estimations than for actual electricity generation.
On the other hand, the peculiarity of the southeast Asian climate should not be forgotten.The winter monsoon season is noted not only by Korean scientists but also by Malaysian [26] and Taiwanese [27] scientists.It was estimated that up to 84% of the total wave power is generated during the winter monsoon season in Taiwan [27].The temporal distribution assessment of the wave energy resource in the Indian shelf seas was even performed by dividing the year into pre-monsoon, south-west monsoon and post-monsoon periods [28].Because wave power fluxes can exceed 100 kW/m during typhoons, the regular typhoons in Southeast Asia also should be taken into consideration.Additionally, adjustments to WEC designs should be made [27].

Australia and New Zealand
While the northern coasts of Australia are described as unsuitable places for electricity generation from waves because of their relatively low resources [29], the southern coasts can only be described as unique.Extremely long fetches provide the conditions for the highest wind waves in the world, and the length of this coast is 3000 km [30].With wave power fluxes reaching 30-50 kW/m, it is estimated that wave energy resources in Australia's southern coasts represent five times the energy requirements of the whole of Australia [31].Therefore, it is not surprising that the capacity factor for the same WECs assessed for their suitability to the east Mediterranean Sea's Lebanon coast exceeds 50%.Moreover, it is estimated that the cost of electricity produced from a wave energy converter situated 5 km from shore can be comparable to electricity generated from wind or sun power plants [31].Although New Zealand can appreciate similar wave power to that of the southern Australia coast, recent developments show that in this country there is no requirement for detailed wave energy assessments.New Zealand is applying a cautious approach to wave energy [1,32].

Africa
Two African countries are members of the IEA-OES: South Africa and Nigeria.An assessment of the wave energy resources in South Africa revealed that it has extensive wave energy resources, especially in the southwestern coast, where in offshore locations wave power fluxes vary from 33 kW/m to 41 kW/m.The near-shore average annual wave power flux is 26 kW/m [33].So far there is no information about such assessments in Nigeria [34].
An assessment of ocean renewable resources in the western Indian Ocean was published in 2012.Wave energy was discussed, along with tidal, ocean current and ocean thermal energy.The assessment was carried out for the coasts of Kenya, Tanzania and Mozambique; and the islands of Madagascar, Reunion and Mauritius [35].It is worth mentioning that this region has one of the lowest electrification rates in the world: the rates in Kenya, Tanzania, Mozambique and Madagascar are less than 20% [36].The average offshore wave power flux value here of 25 kW/m, and the future perspective is that wave power in this region can be used not only for electricity generation but also for desalination of sea water.Unfortunately, the strong tropical cyclones in this area will be a serious limiting factor for WECs installations [35].

Europe
Europe's wave energy resources were recently reviewed in reference [37], which focused more on the European semi-enclosed seas.A browse through the latest published studies reveals two main directions of research concerning wave energy potential in Europe.The first one is the North Atlantic coasts, where Europe has its highest wave energy potential.Here, an international scientific collaboration (the EnergyMare project) to create a map of wave energy resources along the European Atlantic coast was established [38].In the project, the results of WaveWatch III and SWAN models were compared, and wave energy resources in potential wave energy conversion or test sites in Scotland, France, Spain and Portugal were evaluated.According to the results, the four above mentioned countries in descending order of their wave energy resources are Scotland (up to 40 kW/m in Shetlands), France (up to 34 kW/m in Bretagne), Spain (up to 34 kW/m in Estaca de Bares) and Portugal (up to 24 kW/m in Nazaré) [38].Furthermore, the attractiveness of this region for electricity generation from sea waves is reflected in recently published wave energy resource assessments for Ría de Vigo (Spain) [39], the Iberian Peninsula [40], the Sea of Iroise (France) [41], Galway Bay (Ireland) [42] and Ireland [43].The last study not only focuses on wind and wave energy potential in Ireland but also evaluates the correlation between these two renewable power sources.This is important because synergy is possible when converting wind and wave energy into electricity via purposely designed wave power plants such as Wavestar [44].
The second main direction of the wave energy potential studies in Europe is the Mediterranean Sea.New studies are constantly being conducted to update knowledge of the whole Mediterranean Sea's wave energy potential via high resolution models [45].Additionally, high resolution models are used for the regional assessment of wave energy, e.g., the Greek offshore areas [46].A study concerning wave energy resources in the Balearic Sea revealed new "hot spots" in the Mediterranean Sea that are viable for WEC construction-the northern coast of Spain's Menorca Island and the northern, western, and eastern coasts of Spain's Mallorca Island [47].
"Hot spot" is usually understood as an area with the highest concentrated wave power potential, but in the Mediterranean Sea this perception is beginning to change.Because the seasonal variation of wave energy is very high here-in winter, the wave power flux can be up to 5-6 times that of the summer [47]-the actual "hot spot" for WEC in the Mediterranean Sea can be an area with the lowest coefficient of variation, which means a more consistent value of wave energy throughout the year [45].For this reason the studies of the lower wave energy potential regions-areas such as the Aegean Sea, where the average wave power flux is 2-2.5 kW/m [48]-can become more and more important in the future.

Islands
Islands in this overview are presented separately because they are the areas where the supply of energy is usually a very important and expensive issue.For example, in St. George Island (USA), located in the Bering Sea approximately 500 km from the Alaskan shores, energy is generated via four diesel generators.It has been evaluated that a WEC situated 3 km from the shore at 40 m depth, where the average wave power flux is 28 kW/m, would save the local community fuel costing approximately US$81,600 annually [49].There could be considerable support from the islanders (population 100), whose main economic activity is fishing.
Another island of the USA, whose remoteness is the reason for the pursuit of alternative energy sources, is Hawaii.Here, the average wave power flux at 60 m depth varies from 15 kW/m to 30 kW/m.The newly established Hawaii National Marine Renewable Energy Center has already selected an area designated for a wave energy test site [50].
The Azores Islands (Portugal) was one of the first areas where the renewed attention to the marine energy in Europe has materialized.Here, on Pico Island in 1999, an oscillating water column type WEC was built [51].Unfortunately, a lack of funding slowed down the development process in this case, although the wave energy resources in the Azores Islands are twice as large as in Hawaii [52].
The Canary Islands (Spain) are another area in the Northern Atlantic suitable for wave energy conversion.The average wave power flux here is similar to Hawaii at 25-30 kW/m [53].The assessments of the wave energy resources for these islands continue constantly [54].Additionally, the study of the effect of WECs on the incoming waves in the near-shore area of Tenerife was carried out in reference [55].Although the purpose of this study was different, the related idea to use WECs as a shore protection measure is becoming quite common [56].
A European island with one of the highest wave energy potentials is Iceland.The average wave energy flux on Iceland's coast at the 50 m depth can reach 45 kW/m [57].On the other hand, the severe Energies 2017, 10, 2170 6 of 18 climate conditions and the fact that Iceland already generates all of its electricity from renewable sources [58] raises doubt on whether Iceland will press ahead to realize the potential of wave energy.
Another group of islands with impressive wave energy resources is Fiji.The estimated wave energy potential of Fiji is 29 GW.If the country could realize 0.5% of this potential, it would still be sufficient to meet its electricity demand [56].
Tapping the wave energy potential will be important not only for the islands that are situated in the open ocean.In the Mediterranean Sea, the area with the highest wave energy is the western coast of Sardinia [59].Recently, a hot spot was discovered in Sicily's near-shore area [60].In the Baltic Sea the possibility to generate electricity from the waves was studied in the area close to the island of Aland [61].
The worldwide wave energy assessments are summarized in Table 1, where the coastal wave power fluxes are presented in ascending order.Efforts were made to select and present in the table only those wave power flux values that were presented as averaged in the different studies.Table 1 clearly shows that wave energy resources are distributed unevenly around the world and that these resources are highest in the coasts that have the longest fetches.1.7.Lithuania The theoretical wave power resource of the Baltic Sea (excluding areas where the wave power flux is lower than 5 kW/m and potentially ice-covered regions) is estimated to be 1 GW [4].Interestingly, the resource is confined to the southeastern region of the Baltic Sea, where the Lithuanian coast is located (Figure 1a).Swedish scientists estimated that the average wave power flux in the whole Baltic Sea can reach 4-5 kW/m [65], whereas at the near-shore area along the Lithuanian coast, the multi-year average wave power potential reaches 1-2 kW/m [37,66].The assessment of the temporal distribution of the Baltic Sea near-shore wave power resources along the Lithuanian coast was published in reference [37].The spatial distribution of the Baltic Sea near-shore wave power resources along the coast of Klaip ėda (Lithuania), presented in this paper, which is continuation of the work that was published in reference [37], was assessed for the first time.The impacts of wave propagation direction and decreasing depth on wave power resources in the area of low wave power potential were examined.This upgraded knowledge on the Baltic Sea near-shore wave power resources is also used in this paper to show where the Baltic Sea stands according to the wave power potential with reference to presented global distribution of these resources.

Numerical Model
The MIKE 21 Nearshore Spectral Wind-Wave (NSW) model [67] was used to model the dynamics of wave height loss.The results were consequently the initial data used to assess the decreasing amount of the energy of waves propagating to shore.MIKE 21 NSW is a stationary, directionally decoupled parametric model that takes into account the effects of refraction and shoaling due to varying depth, local wind generation and energy dissipation from bottom friction and wave breaking.
The main equation of MIKE 21 NSW is the conservation equation for the spectral wave action density, which is solved via parameterization of the conservation equation in the frequency domain and the use of the zeroth and the first moments of the action spectrum as dependent variables.

Numerical Model
The MIKE 21 Nearshore Spectral Wind-Wave (NSW) model [67] was used to model the dynamics of wave height loss.The results were consequently the initial data used to assess the decreasing amount of the energy of waves propagating to shore.MIKE 21 NSW is a stationary, directionally decoupled parametric model that takes into account the effects of refraction and shoaling due to varying depth, local wind generation and energy dissipation from bottom friction and wave breaking.
The main equation of MIKE 21 NSW is the conservation equation for the spectral wave action density, which is solved via parameterization of the conservation equation in the frequency domain and the use of the zeroth and the first moments of the action spectrum as dependent variables.
The wave energy dissipation due to decreasing depth is controlled in the MIKE 21 NSW via bottom friction, the process by which the wave loses some of its energy due to the effect of friction at the sea bottom.The model's equation for the energy dissipation due to the bottom friction is based on the quadratic friction law [68]: where E = (H rms ) 2 /8, ω is the frequency, H rms is the root mean square wave height, k is the wave number, h is the water depth, and C fw is the wave friction factor.The wave friction factor (C fw = f w /2) is specified by the Nikuradse roughness parameter [69]: k N is the Nikuradse roughness parameter, and a b is amplitude of the particle motion at the bottom.The Nikuradse roughness parameter can be estimated by [70]: where d 50 is the median grain size of the sediment.In the presence of ripples, the Nikuradse roughness parameter increases, and its estimate should include the ripple characteristics.Its value increases further if there is vegetation.It is difficult to assess this parameter, and therefore, it is used for model calibration.
The process of wave breaking due to large wave steepness and limiting depth in MIKE 21 NSW is based on the formulation for energy dissipation rate due to wave breaking [71]: where: where E is the total energy, ω is the frequency, H rms is the root mean square wave height, H m is the maximum allowable wave height, Q b is the fraction of breaking waves, and α is an adjustable constant.
The maximum wave height is estimated by [71]: where γ 1 and γ 2 are wave breaking parameters: γ 1 controls the steepness condition, and γ 2 controls the limiting water depth condition.In this case, α, γ 1 and γ 2 are specified as constants for the whole model area, and they have values of 1.0, 0.8 and 1.0, respectively.The location of the model area is presented in Figure 1b.This location was selected for its several advantages.Firstly, there is the possibility of using the test data of the bottom mounted wave recorder Aanderaa SeaGuard [72], operated by the Lithuanian Environmental Protection Agency's Department of Marine Research.Secondly, the Klaip ėda Seaport breakwaters could be a suitable site for the installation of a wave energy converter in Lithuania.

Model Calibration
The model calibration was performed using instrumental measurements of the wave heights at 10 m depth (55 • 43 55 N, 21 • 4 20 E), i.e., the test data of the bottom mounted wave recorder and the ECMWF ERA-Interim reanalysis data [73] from 50 m depth (55 • 50 N, 20 • 25 E; Figure 1b).Term reanalysis means that ECMWF forecast is assimilated with measurements.In such way updated results are kept as ERA-Interim data-a reliable information on atmosphere, land and oceans.
The model was run in quasi-stationary mode with the additional supplementation of the time series of the model coastal areas' wind speed and direction for the simulation period.ERA-Interim data were used as offshore wave height boundary conditions for the higher resolution near-shore model.Of all the parameters that can be used in MIKE 21 NSW for calibration, the bottom friction was chosen.The increase in the bottom friction coefficient of the near-shore waters is the cause of the higher energy dissipation rate of the waves and vice versa.In this study, the bottom friction coefficient was specified via the Nikuradse roughness parameter.
For calibration purposes, the modelled and measured significant wave height data sets at the bottom mounted wave recorder area point were compared for the period from the 1 to 5 September 2012 with a model time step of 6 h.The model accuracy is assessed by calculating the root mean square error, bias and the correlation coefficient [7,21].These values for various magnitudes of the Nikuradse roughness parameter are presented in the Table 2.The bias values indicate that the modelled significant wave heights are higher than the measured ones.Hence, the Nikuradse roughness parameter could be increased further, especially because this enhancement manifests in the constantly decreasing root mean square error and bias values.The decision to finalize the calibration process was reached with reference to the variations of the correlation coefficient values.Figure 2 at least partially explains these variations: it shows the comparison between the measured and modelled significant wave height values.
The model was run in quasi-stationary mode with the additional supplementation of the time series of the model coastal areas' wind speed and direction for the simulation period.ERA-Interim data were used as offshore wave height boundary conditions for the higher resolution near-shore model.
Of all the parameters that can be used in MIKE 21 NSW for calibration, the bottom friction was chosen.The increase in the bottom friction coefficient of the near-shore waters is the cause of the higher energy dissipation rate of the waves and vice versa.In this study, the bottom friction coefficient was specified via the Nikuradse roughness parameter.
For calibration purposes, the modelled and measured significant wave height data sets at the bottom mounted wave recorder area point were compared for the period from the 1 to 5 September 2012 with a model time step of 6 h.The model accuracy is assessed by calculating the root mean square error, bias and the correlation coefficient [7,21].These values for various magnitudes of the Nikuradse roughness parameter are presented in the Table 2.The bias values indicate that the modelled significant wave heights are higher than the measured ones.Hence, the Nikuradse roughness parameter could be increased further, especially because this enhancement manifests in the constantly decreasing root mean square error and bias values.The decision to finalize the calibration process was reached with reference to the variations of the correlation coefficient values.Figure 2 at least partially explains these variations: it shows the comparison between the measured and modelled significant wave height values.It is clear from Figure 2. that to obtain lower wave heights, it is not required to increase further the Nikuradse roughness parameter, since the lower wave heights are not necessary in all cases.Furthermore, in Figure 2 there is a visible shift between measured and modelled significant wave heights on 4 January and possible error in measured significant wave heights on 5 January.This is the reason for the variations of the correlation coefficient.Taking into account these factors, a Nikuradse roughness parameter value of 0.005 was selected (also indicated in bold in Table 2) for the further modelling of the dynamics of wave height loss in the Baltic Sea near-shore area.Applying this value to the calibration model runs not only the highest value of the correlation coefficient but also acceptable values of bias and root mean square error was reached.It is clear from Figure 2. that to obtain lower wave heights, it is not required to increase further the Nikuradse roughness parameter, since the lower wave heights are not necessary in all cases.Furthermore, in Figure 2 there is a visible shift between measured and modelled significant wave heights on 4 January and possible error in measured significant wave heights on 5 January.This is the reason for the variations of the correlation coefficient.Taking into account these factors, a Nikuradse roughness parameter value of 0.005 was selected (also indicated in bold in Table 2) for the further modelling of the dynamics of wave height loss in the Baltic Sea near-shore area.Applying this value to the calibration model runs not only the highest value of the correlation coefficient but also acceptable values of bias and root mean square error was reached.

Spatial Distribution of the Wave Energy Potential
Average seasonal and annual wave heights for the design years (Table 3) were selected as offshore conditions to model the dynamics of wave height loss.The methodology for the calculation of these wave heights was published in reference [37].
To improve the possibility of assessing the spatial distribution of the wave power potential in the Baltic Sea along the coast of Klaip ėda, the wave power of the design waves was estimated alongside the 5, 10, 15, 20, 25 and 30 m depth isobaths.Alongside these isobaths, calculation points (Figure 1c) were selected, and at each point, the wave power flux of the wave propagating to shore was calculated.Each wave was treated as a separate case, and the stationary mode of the MIKE 21 NSW model was adapted to assess the spatial distribution of the wave power potential along the coast of Klaip ėda, taking into account the impacts of depth, wave propagation direction and wind.The wave power flux in this study was using the equation: where ρ is the mass density of water (in this case it is brackish water of the Baltic Sea with density 1010 kg/m 3 [74], g is the gravitational constant (9.81 m/s 2 ), S(f ) is the spectral density, and c g is the group velocity.The spectral density was calculated using the JONSWAP wave spectrum modified for the Baltic Sea [75]: where K m is an empirically determined constant (0.1786), H s is a significant wave height, T p is the peak wave period, f p is the peak frequency, f is the wave frequency, and γ is the peak enhancement factor (4.0), σ = 0.07 for f ≤ f p , σ = 0.09 for f > f p , where σ is shape parameter.

The Impact of Wave Propagation Direction
To assess the impact of the wave propagation direction on the wave power potential, three wave propagation directions were selected: west, northwest and southwest.Examples of the change of design heights for the waves propagating to shore from different directions and reaching 5 m depth are presented in Figure 3.The first observation is that the highest waves reach the shore when they propagate from the west.This is because this direction and the shoreline are almost perpendicular.On the other hand, the considered length of the Lithuanian shore is not parallel to the north direction; therefore, the waves that propagate from the northwest and the southwest lose a different amount of energy due to the refraction process.A lesser amount of energy is lost when waves propagate from the southwest.Of course, for the low waves, the difference in wave height loss is only a few centimetres (for example, for a low intensity year's waves, the difference does not exceed 0.02 m), a quantity that does not correspond to high wave power losses.However, when higher waves are examined, the impact of wave direction becomes more evident.For example, the highest difference for the autumn average wave height during a high intensity wave year is 0.14 m.
Energies 2017, 10, 2170 11 of 18 the southwest.Of course, for the low waves, the difference in wave height loss is only a few centimetres (for example, for a low intensity year's waves, the difference does not exceed 0.02 m), a quantity that does not correspond to high wave power losses.However, when higher waves are examined, the impact of wave direction becomes more evident.For example, the highest difference for the autumn average wave height during a high intensity wave year is 0.14 m.Another important aspect of the finding that waves lose lesser amount of energy when they are propagating from the southwest than from the northwest becomes clear when these findings are compared with wave rose (Figure 4), created for ERA-Interim grid point (Figure 1b).The fact that not only there are more waves in total that are propagating from the southwest, but also higher amount of higher and more powerful waves propagating from the same direction can influence the selection of the site for wave energy converter installation.For example, the site in question could be the Klaipėda Seaport breakwaters (Figure 1c).There are two breakwaters currently constructed in the Klaipėda Seaport: northern breakwater which is open for waves that are propagating from the north and northwest and southern breakwater which is open for waves that propagate from the south and southwest.Available information from the wave rose suggest that firstly for WEC installation should be considered the southern breakwater.The fact that waves that are propagating from the southwest lose lesser amount of energy is additional great advantage in this case.The study clearly revealed that when lower wave power areas are under consideration, the assessment of the wave height loss for greater near-shore depths (30-50 m) is not that essential.Another important aspect of the finding that waves lose lesser amount of energy when they are propagating from the southwest than from the northwest becomes clear when these findings are compared with wave rose (Figure 4), created for ERA-Interim grid point (Figure 1b).The fact that not only there are more waves in total that are propagating from the southwest, but also higher amount of higher and more powerful waves propagating from the same direction can influence the selection of the site for wave energy converter installation.For example, the site in question could be the Klaip ėda Seaport breakwaters (Figure 1c).There are two breakwaters currently constructed in the Klaip ėda Seaport: northern breakwater which is open for waves that are propagating from the north and northwest and southern breakwater which is open for waves that propagate from the south and southwest.Available information from the wave rose suggest that firstly for WEC installation should be considered the southern breakwater.The fact that waves that are propagating from the southwest lose lesser amount of energy is additional great advantage in this case.the southwest.Of course, for the low waves, the difference in wave height loss is only a few centimetres (for example, for a low intensity year's waves, the difference does not exceed 0.02 m), a quantity that does not correspond to high wave power losses.However, when higher waves are examined, the impact of wave direction becomes more evident.For example, the highest difference for the autumn average wave height during a high intensity wave year is 0.14 m.Another important aspect of the finding that waves lose lesser amount of energy when they are propagating from the southwest than from the northwest becomes clear when these findings are compared with wave rose (Figure 4), created for ERA-Interim grid point (Figure 1b).The fact that not only there are more waves in total that are propagating from the southwest, but also higher amount of higher and more powerful waves propagating from the same direction can influence the selection of the site for wave energy converter installation.For example, the site in question could be the Klaipėda Seaport breakwaters (Figure 1c).There are two breakwaters currently constructed in the Klaipėda Seaport: northern breakwater which is open for waves that are propagating from the north and northwest and southern breakwater which is open for waves that propagate from the south and southwest.Available information from the wave rose suggest that firstly for WEC installation should be considered the southern breakwater.The fact that waves that are propagating from the southwest lose lesser amount of energy is additional great advantage in this case.The study clearly revealed that when lower wave power areas are under consideration, the assessment of the wave height loss for greater near-shore depths (30-50 m) is not that essential.The study clearly revealed that when lower wave power areas are under consideration, the assessment of the wave height loss for greater near-shore depths (30-50 m) is not that essential.During a low intensity wave year, all waves, with two exceptions (average spring and autumn waves, propagating from southwest), reached a depth of 10 m without losing any height.Still, this proves the significance of depth for installations of wave energy converters, especially when the areas of lower wave power potential are being studied.
Next, the wave power resources were estimated by summing up the wave power flux values along the isobaths (Figure 1c).It must be emphasized that these resources were calculated by assessing the wave height variations along the isobaths, that is, by summing up the values of the sections at which centre the calculation points were located.
The results showed that wave power resources in the Baltic Sea near-shore area along the coast of Klaip ėda decrease starting from 25 m depth.This is misleading since the isobaths of 25 m and 30 m in the model area are shorter in length (Figure 1c).According to the spectral wave energy equations, wave power flux must increase with the increasing depth.Therefore, wave power resources in this study were recalculated to demonstrate the amount of wave power per kilometre of isobath, which also corresponds to the averaged wave power fluxes of the design waves at the different depths in the Baltic Sea near-shore area.The recalculations are not shown in detail in this paper.The changes in wave power fluxes of the design waves corresponding to the wave height losses at the 5 m depth (presented in Figure 3) are shown in Figure 5.
Energies 2017, 10, 2170 12 of 18 During a low intensity wave year, all waves, with two exceptions (average spring and autumn waves, propagating from southwest), reached a depth of 10 m without losing any height.Still, this proves the significance of depth for installations of wave energy converters, especially when the areas of lower wave power potential are being studied.Next, the wave power resources were estimated by summing up the wave power flux values along the isobaths (Figure 1c).It must be emphasized that these resources were calculated by assessing the wave height variations along the isobaths, that is, by summing up the values of the sections at which centre the calculation points were located.
The results showed that wave power resources in the Baltic Sea near-shore area along the coast of Klaipėda decrease starting from 25 m depth.This is misleading since the isobaths of 25 m and 30 m in the model area are shorter in length (Figure 1c).According to the spectral wave energy equations, wave power flux must increase with the increasing depth.Therefore, wave power resources in this study were recalculated to demonstrate the amount of wave power potential per kilometre of isobath, which also corresponds to the averaged wave power fluxes of the design waves at the different depths in the Baltic Sea near-shore area.The recalculations are not shown in detail in this paper.The changes in wave power fluxes of the design waves corresponding to the wave height losses at the 5 m depth (presented in Figure 3) are shown in Figure 5.The impact of the wave propagation directions on the Baltic Sea near-shore wave power potential along the coast of Klaipėda is revealed best in Figure 4.When higher waves (~1.0 m) propagating to shore reach 5 m depth, they lose approximately 36% of their power if propagating from the west, 41% if propagating from the southwest, and even 50% if propagating from the northwest.Waves of average height (~0.7 m) lose 35%, 38% and 45%, respectively.Lower waves (~0.5 m) lose approximately 31% when propagating from the west and the southwest and 36% when propagating from the northwest.It is highly likely that due to the refraction process, waves propagating from the north and the south would lose even more of their power.

The Impact of Decreasing Depth
The impact of decreasing depth on the Baltic Sea near-shore wave power potential along the coast of Klaipėda is assessed for the 55°45′ N parallel starting from 20 m depth.Along this line, on the basis of modelled wave height loss, the wave power fluxes were calculated.The location of the line, north from of the Klaipėda Seaport, was selected because the territory south of the seaport is the Curonian Spit, which is not only a national designated area but also a UNESCO world heritage site.Therefore, the construction of WECs in this near-shore area will be strongly questionable.
To examine the impact of decreasing depth on wave energy resources, three wave heights were selected: the minimum seasonal (low intensity wave year summer average) wave height (0.48 m), the The impact of the wave propagation directions on the Baltic Sea near-shore wave power potential along the coast of Klaip ėda is revealed best in Figure 4.When higher waves (~1.0 m) propagating to shore reach 5 m depth, they lose approximately 36% of their power if propagating from the west, 41% if propagating from the southwest, and even 50% if propagating from the northwest.Waves of average height (~0.7 m) lose 35%, 38% and 45%, respectively.Lower waves (~0.5 m) lose approximately 31% when propagating from the west and the southwest and 36% when propagating from the northwest.It is highly likely that due to the refraction process, waves propagating from the north and the south would lose even more of their power.

The Impact of Decreasing Depth
The impact of decreasing depth on the Baltic Sea near-shore wave power potential along the coast of Klaip ėda is assessed for the 55 • 45 N parallel starting from 20 m depth.Along this line, on the basis of modelled wave height loss, the wave power fluxes were calculated.The location of the line, north from of the Klaip ėda Seaport, was selected because the territory south of the seaport is the Curonian Spit, which is not only a national designated area but also a UNESCO world heritage site.Therefore, the construction of WECs in this near-shore area will be strongly questionable.
To examine the impact of decreasing depth on wave energy resources, three wave heights were selected: the minimum seasonal (low intensity wave year summer average) wave height (0.48 m), the maximum seasonal (high intensity wave year autumn average) wave height (1.25 m), and the median intensity wave year annual average wave height (0.67 m).Into this range falls 69.9% of the average monthly wave heights from the 1970-2010 period, which were used to determine the distribution of wave heights during the design years.Hence, the obtained curves that describe the impact of the decreasing depth on the wave energy resources in this study reflect the situation during the majority of the year in the Baltic Sea near-shore area along the coast of Klaip ėda, when waves propagate from western directions.Figures 6-8 show the impact of decreasing depth for studied waves that propagate from western directions.
Energies 2017, 10, 2170 13 of 18 maximum seasonal (high intensity wave year autumn average) wave height (1.25 m), and the median intensity wave year annual average wave height (0.67 m).Into this range falls 69.9% of the average monthly wave heights from the 1970-2010 period, which were used to determine the distribution of wave heights during the design years.Hence, the obtained curves that describe the impact of the decreasing depth on the wave energy resources in this study reflect the situation during the majority of the year in the Baltic Sea near-shore area along the coast of Klaipėda, when waves propagate from western directions.Figures 6-8 show the impact of decreasing depth for studied waves that propagate from western directions.From Figures 6 and 7, it is easy to spot the depth at which a wave starts to get influenced by sea bottom and lose its power due to the decreasing depth.For a 0.67 m wave, this depth is 14 m, while for 0.48 m wave, it is 10 m.However, determining this depth for a 1.25 m wave is still difficult (Figure 8).The amount of wave power that will reach shallower depths depends on both the wave height and the wave propagating direction.From Figures 6 and 7, it is easy to spot the depth at which a wave starts to get influenced by sea bottom and lose its power due to the decreasing depth.For a 0.67 m wave, this depth is 14 m, while for 0.48 m wave, it is 10 m.However, determining this depth for a 1.25 m wave is still difficult (Figure 8).The amount of wave power that will reach shallower depths depends on both the wave height and the wave propagating direction.From Figures 6 and 7, it is easy to spot the depth at which a wave starts to get influenced by sea bottom and lose its power due to the decreasing depth.For a 0.67 m wave, this depth is 14 m, while for 0.48 m wave, it is 10 m.However, determining this depth for a 1.25 m wave is still difficult (Figure 8).The amount of wave power that will reach shallower depths depends on both the wave height and the wave propagating direction.The energy losses due to decreasing depth were calculated.These losses are much more substantial compared with the losses due to the propagation direction.If the 1.25 m wave propagates from the most energetically favourable direction (west) and reaches 5 m depth, it loses 0.91 kW/m or 34.0% of its power.If the same wave is propagating from the least favourable direction in this case (northwest), it loses an additional 0.58 kW/m or 21.6% of its power.It is a similar story for the 0.67 m and 0.48 m waves.The former loses 0.24 kW/m or 30.8% when propagating from the most favourable direction and an additional 0.08 kW/m or 10.2% when propagating from the least favourable direction.The latter loses 0.11 kW/m or 27.5% as well as an additional 0.01 kW/m or 2.5% of its power when propagating from the most favourable and least favourable directions, respectively.

Discussion and Conclusions
This paper is continuation of the study that was partially published in [37].This study is different from the usual assessments of near-shore wave power potentials.Not only it is the first spatiotemporal assessment of wave power potential for the Baltic Sea area that in theory has one of the highest wave power potentials, but it is also performed in the area that has almost no instrumental measurements, but where instead long term (over 50 years) visual near-shore wave observations are available.These were used as initial data for the study.Statistical method designed for calculating a distribution of annual hydrological variables using multi-year data was applied together with probability distribution analysis.In such way design waves for high intensity, median intensity and low intensity wave years were calculated and parametrized JONSWAP wave energy spectrum modified for the Baltic Sea was used to estimate wave power flux of design waves.Same design waves were used for assessment of spatial distribution of wave power potential, presented in this paper.From performed extensive literature review it can be stated that the assessment of spatial distribution of wave power potential taking into account wave propagation direction and decreasing depth is still a novelty, especially in low wave energy potential areas.
Since wave power resources depend on fetch, they are distributed unevenly around the globe.The Baltic Sea's wave power resources are low compared to those of oceans yet are similar to those of other semi-enclosed seas, such as the Black Sea or the Persian Gulf.Nevertheless, the amount of published studies on the wave power potential in the Baltic Sea is constantly growing.It is already known, that the average wave power flux in the Baltic Sea can reach 4-5 kW/m [64].In the coastal "hot spots" areas average wave power flux can reach 2.55 kW/m at non-sheltered condition such as island of Saaremaa [66] and approximately 1.6 kW/m in sheltered conditions such as Bay of Gdansk [76].In the Lithuanian coast of the Baltic Sea the average wave power flux can reach 1-2 kW/m [37,66], therefore it can be further assessed as possible location for electricity generation from waves.This study represents a first attempt to conduct a study of spatial distribution of the The energy losses due to decreasing depth were calculated.These losses are much more substantial compared with the losses due to the propagation direction.If the 1.25 m wave propagates from the most energetically favourable direction (west) and reaches 5 m depth, it loses 0.91 kW/m or 34.0% of its power.If the same wave is propagating from the least favourable direction in this case (northwest), it loses an additional 0.58 kW/m or 21.6% of its power.It is a similar story for the 0.67 m and 0.48 m waves.The former loses 0.24 kW/m or 30.8% when propagating from the most favourable direction and an additional 0.08 kW/m or 10.2% when propagating from the least favourable direction.The latter loses 0.11 kW/m or 27.5% as well as an additional 0.01 kW/m or 2.5% of its power when propagating from the most favourable and least favourable directions, respectively.

Discussion and Conclusions
This paper is continuation of the study that was partially published in [37].This study is different from the usual assessments of near-shore wave power potentials.Not only it is the first spatio-temporal assessment of wave power potential for the Baltic Sea area that in theory has one of the highest wave power potentials, but it is also performed in the area that has almost no instrumental measurements, but where instead long term (over 50 years) visual near-shore wave observations are available.These were used as initial data for the study.Statistical method designed for calculating a distribution of annual hydrological variables using multi-year data was applied together with probability distribution analysis.In such way design waves for high intensity, median intensity and low intensity wave years were calculated and parametrized JONSWAP wave energy spectrum modified for the Baltic Sea was used to estimate wave power flux of design waves.Same design waves were used for assessment of spatial distribution of wave power potential, presented in this paper.From performed extensive literature review it can be stated that the assessment of spatial distribution of wave power potential taking into account wave propagation direction and decreasing depth is still a novelty, especially in low wave energy potential areas.
Since wave power resources depend on fetch, they are distributed unevenly around the globe.The Baltic Sea's wave power resources are low compared to those of oceans yet are similar to those of other semi-enclosed seas, such as the Black Sea or the Persian Gulf.Nevertheless, the amount of published studies on the wave power potential in the Baltic Sea is constantly growing.It is already known, that the average wave power flux in the Baltic Sea can reach 4-5 kW/m [64].In the coastal "hot spots" areas average wave power flux can reach 2.55 kW/m at non-sheltered condition such as island of Saaremaa [66] and approximately 1.6 kW/m in sheltered conditions such as Bay of Gdansk [76].In the Lithuanian coast of the Baltic Sea the average wave power flux can reach 1-2 kW/m [37,66], therefore it can be further assessed as possible location for electricity generation from waves.This study represents a first attempt to conduct a study of spatial distribution of the near-shore wave power resources in the southeastern region of the Baltic Sea, taking into account wave propagation direction and decreasing depth.
The assessment of the impact of the wave propagation direction on wave power resources in the Baltic Sea near-shore area along the coast of Klaip ėda revealed that the most powerful waves reach the Lithuanian coast when they propagate from the west.The most powerful in this case means that these waves loses the least of their power while propagating to the shore.
When lower waves (~0.5 m) are propagating to the shore from the west starting from 30 m and reaching 5 m depth they lose approximately 31% of their power flux.Change of propagation direction decreases the wave height of approximately 0.02 m and increases the losses just to 36%.
When higher waves (e.g., 1.25 m) are studied, their initial percentage loss of power flux due to the decreasing depth is similar-34.0%,however the change of propagation direction can increase the loses to 50% or more.Of course, the numerical value of this power loss is significantly higher than for lower waves.
Still, typical waves in the Baltic Sea near-shore area along the coast of Klaip ėda are low, and therefore, they begin to be impacted by decreasing depth and lose their power at relatively low (10-14 m) depths.However, to turn this into an advantage and to avoid a rapid decrease in wave power flux due to the decreasing depth, the wave power converter would have to work efficiently under low power conditions.Having similar spatial distribution of the design waves' power fluxes in various depths and the calculations of their power flux losses due to the propagation to the shore can help to assess the most economically feasible site for the wave power converter.

Figure 1 .
Figure 1.(a) Location of the study area in the southeastern part of the Baltic Seal; (b) model area together with the available wave data points used for model calibration; (c) bathymetry of the model area and calculation points.

Figure 1 .
Figure 1.(a) Location of the study area in the southeastern part of the Baltic Seal; (b) model area together with the available wave data points used for model calibration; (c) bathymetry of the model area and calculation points.

Figure 2 .
Figure 2. Comparison between measured and modelled significant wave heights in the Baltic Sea near-shore area along the coast of Klaipėda for a Nikuradse roughness parameter value of 0.005.

Figure 2 .
Figure 2. Comparison between measured and modelled significant wave heights in the Baltic Sea near-shore area along the coast of Klaip ėda for a Nikuradse roughness parameter value of 0.005.

Figure 3 .
Figure 3.The change in height of waves propagating to the shore.

Figure 4 .
Figure 4. Wave rose for year 2012 (significant wave heights) at ERA-Interim grid point.

Figure 3 .
Figure 3.The change in height of waves propagating to the shore.

Figure 3 .
Figure 3.The change in height of waves propagating to the shore.

Figure 4 .
Figure 4. Wave rose for year 2012 (significant wave heights) at ERA-Interim grid point.

Figure 4 .
Figure 4. Wave rose for year 2012 (significant wave heights) at ERA-Interim grid point.

Figure 5 .
Figure 5.The change in power flux of waves propagating to the shore.

Figure 5 .
Figure 5.The change in power flux of waves propagating to the shore.

Figure 6 .
Figure 6.The impact of decreasing depth on wave power flux for a 1.25 m wave.

Figure 7 .
Figure 7.The impact of decreasing depth on wave power flux for a 0.67 m wave.

Figure 6 .
Figure 6.The impact of decreasing depth on wave power flux for a 1.25 m wave.

Figure 6 .
Figure 6.The impact of decreasing depth on wave power flux for a 1.25 m wave.

Figure 7 .
Figure 7.The impact of decreasing depth on wave power flux for a 0.67 m wave.

Figure 7 .
Figure 7.The impact of decreasing depth on wave power flux for a 0.67 m wave.

Figure 8 .
Figure 8.The impact of decreasing depth on wave power flux for a 0.48 m wave.

Figure 8 .
Figure 8.The impact of decreasing depth on wave power flux for a 0.48 m wave.

Table 1 .
The distribution of wave energy resources worldwide.

Table 2 .
The results of the Baltic Sea's near shore wave model calibration.

Table 2 .
The results of the Baltic Sea's near shore wave model calibration.

Table 3 .
Average seasonal and annual wave heights for the design years along the coast of Klaip ėda.