An Overview of the Expected Shoreline Impact of the Marine Energy Farms Operating in Different Coastal Environments

The aim of the present work is to provide an overview of the possible implications involving the influence of a generic marine energy farm on the nearshore processes. Several case studies covering various European coastal areas are considered for illustration purposes. These include different nearshore areas, such as the Portuguese coast, Sardinia Island or a coastal sector close to the Danube Delta in the Black Sea. For the case studies related to the Portuguese coast, it is noted that a marine energy farm may reduce the velocity of the longshore currents, with a complete attenuation of the current velocity for some case studies in the coastal area from Leixoes region being observed. For the area located close to the Danube Delta, it is estimated that in the proposed configuration, a marine energy farm would provide an efficient protection against the wave action, but it will have a relatively negligible impact on the longshore currents. Summarizing the results, we can conclude that a marine energy farm seems to be beneficial for coastal protection, even in the case of the enclosed areas, such as the Mediterranean or Black seas, where the erosion generated by the wave action represents a real problem.


Introduction
A significant part of the global population lives near the coastal areas. It is estimated that almost 20% of the total population is located in a strip area of 25 km from the coast. The share increases to up to 40% of the total population when extending the considered strip area to 100 km. A particularity of these areas is that they are very dynamic environments, which present an annual urban growth of 2.6%. At the same time, the number of coastal cities has itself increased about 4.5 times since 1950 [1,2]. Although these areas are defined by numerous opportunities, they are also facing some threats coming from the surrounding environment, such as coastal erosion, a natural event during which the balance between accretion and erosion is continuously shifting [3][4][5][6]. These natural events occur on various spatial and temporal scales, being influenced by many factors, such as wind, waves and nearshore currents. Moreover, it is expected that climate change will have a negative impact on coastal areas in the future, with sea levels being supposed to rise and the marine conditions becoming more aggressive [7][8][9].
Although multiple parameters shape the coastline processes, the wave action has a significant impact on the coastal erosion [10][11][12]. This aspect was highlighted by Bacino et al. [13], when assessing a beach sector (Samborombón Bay) near the Argentinean coast. In this case, a direct correlation was established between the higher levels of the wave energy and the severe erosion. The wave run-up erosion was studied by Feagin et al. [14], who concluded that dune vegetation will significantly reduce the coastal erosion caused by the storm events. In the work of Serafim et al. [15], the importance of considering the wave action in the set-up of a coastal management plan was highlighted through a case study focused on the Santa Catarina region (southern Brazil). An extreme event is expected to have a negative impact on the beach stability, as in the case of the storm Xynthia that hit the French Atlantic coast in February 2010 [16]. During this event, most of the coastal dunes were affected by breaches, with a strong shoreline retreat being noticed for Saint-Hilaire-de-Riez (−5 m) and La Tranche-sur-Mer (−15 m). The correlation between the storm intensity and coastal erosion was also addressed in Montreuil et al. [17], where an erosion index was established.
At present, most of the erosion management strategies involve the use of hard solutions, such as seawalls, interlocking blocks, rubble mound, detached breakwaters, groins, dunes, tyres or jetties [18][19][20][21]. Nevertheless, there are studies indicating that this approach is not the best way to tackle coastal erosion. In this sense, either soft solutions on a standalone basis or an optimal combination of the soft and hard solutions are expected to provide better results. This aspect was highlighted by Liu et al. [22] when investigating China's coastal engineering system, as well as by Pranzini [23] when focusing on the Italian costal area and by Williams et al. [24] who provided a general overview of this topic.
As the industry evolves, some new technological solutions emerge, such as renewable systems, with wave energy converters (WECs) [25][26][27][28][29] being a prime example. Since the purpose of a WEC farm is to extract energy from the waves, the use of such projects for coastal protection was developed in recent years. Taking into account that at this moment there are no operational wave farms, most of the studies have been based on numerical simulations in which various "what-if" scenarios have been considered. This is the case of Bento et al. [30], Rusu and Onea [31] or Zanopol et al. [32], who considered for assessment the effects induced by a wave farm on different coastal environments. A common way to implement a wave farm is by using a line (or an area), defined as a single obstacle which is characterized by various absorption coefficients. This approach is commonly used for the Portuguese environment [33,34]. Furthermore, it is possible to simulate a wave farm that incorporates multiple WECs, this method being used in Diaconu and Rusu [35] or Rusu and Guedes Soares [36].
A full understanding of a wave farm's impact can be done by considering various transmission coefficients, that can go from 0% (compact farm) up to 90% (widely spaced WECs). This particular aspect was considered in Millar et al. [37] and in Zanopol et al. [38]. A more sophisticated approach involves the combination of a wave propagation model and a coastal circulation model in order to highlight the impact of a wave project on the beach stability. Following this approach, there are numerous studies indicating that the implementation of a marine energy farm is a viable solution for coastal protection [39][40][41]. Driven by the fact that the offshore wind sector is already a mature market, some other works are focused on the analysis of dual wind-wave farms considering various point of views. For example, in Astariz and Iglesias [42,43], the possibility to reduce the downtime period was evaluated by taking into account some specific wind projects, such as Alpha Ventus and Horns Rev I. Various wave farm configurations were also proposed in order to improve the accessibility for operation tasks. Since the wave energy sector is moving forward, the next step is to assemble various wave farms and test their operation. Therefore, it is important to predict in advance the expected coastal impact of such a project [44][45][46][47].
By looking at similar works, we can notice that most of the discussions have been focused on the nearshore impact related to the changes induced in the local wave conditions and little attention has been given to the coastal hydrodynamics and sediment transport processes. As a consequence, the present work aims to cover this gap by providing some insights regarding the littoral drift budget that may be encountered in the presence of a marine energy farm. From the knowledge of the authors, another element of novelty is that for the first time this work considers a generic marine energy farm defined by the same characteristics (length and distance to the shoreline) that are assumed to be implemented in different coastal environments (two located in the ocean and two in enclosed seas). Therefore, the objectives are: (a) to establish the spatial distribution of the significant wave heights in the presence of a generic wave farm. In addition to this, the shoreline variation of various wave parameters will be considered for investigation; (b) to identify the distribution of the longshore currents (especially as regards the maximum velocity); (c) to determine the implications for the coastal hydrodynamics by evaluating the connection between the significant wave heights, the longshore current velocities and the sediment transport.

Target Areas
Four different European coastal areas were considered in this work, two at the North Atlantic Ocean and other two in sea environments (Mediterranean and Black Sea, respectively). These four target areas are illustrated in Figure 1. The first two are related to the nearshore areas facing the ocean environment, more precisely on the Portuguese continental coast (zone A). In the south of this zone, a coastal sector located close to the Sines peninsula was evaluated, while, in the north, a sector located close to the Leixoes area (north of the city of Porto) was considered. Going east, we identified a second target area (zone B) that is located in the Mediterranean Sea. In this case we evaluated an island environment, Sardinia. The north-western part of the Black Sea was also taken into account (zone C), more precisely the Saint George sector (in the Romanian nearshore), which is part of the Danube Delta.
A first step in the evaluation of the coastal impact is related to the identification of some relevant environmental conditions. In this case, such conditions would be the most important wave parameters, namely: significant wave height-Hs (in meters); mean wave period-Tm (in seconds); mean wave direction-Dir (in degree). Table 1 summarizes the wave statistics. The main idea was to consider approximately the same wave conditions, which are identified as most relevant, in all areas in order to assess and compare the coastal response to the presence of the marine energy farms. From this perspective, it must be highlighted that, although all the four different wave conditions defined are realistic for all the coastal environments considered, they represent in general different categories for the ocean waves than for the sea waves. For example, what represents total time average for the ocean waves is very close to the winter time average in the case of the sea waves and the conditions corresponding to a regular storm in the European side of the North Atlantic is close to a high storm in the nearshores considered for the Mediterranean and Black seas. This assumption was based on various analyses that were performed on the general characteristics of the wave climates in the areas targeted. See, for example, the results presented in [31][32][33][34]38,48].

The ISSM Model System and the Case Studies Considered
The case studies presented in this work were processed by using the ISSM (Interface for SWAN and Surf Models) modelling tool that combines a wave model with a surf model [49][50][51]. As a first step, the SWAN spectral wave model was used to assess the wave transformation in the coastal areas. This model solves the wave action balance equation that describes the variation of the wave spectrum in geographical, time and spectral spaces. This equation can be summarized as [52]: where N represents the action density spectrum, σ is the relative frequency, θ is the wave direction, and → U is the velocity of the ambient current (assumed to be uniform). The propagation velocities of the wave energy are the group velocity → c g in physical space σ and c θ = θ in spectral space. S represents the sink and source terms. A more detailed assessment of the nearshore processes (including the longshore currents) was provided by the Navy Standard Surf Model (denoted as Surf) [50], that uses as input in the present case the results coming from the SWAN model, more precisely the parameters Hs, Tm and Dir. Tm and Dir are used as direct inputs, while for the wave height it is required to use the root mean square wave height (Hrms) and therefore this parameter is deduced from the significant wave height with the relationship Hrms = 0.707Hs. As for the longshore currents variations, the Surf model uses the following expression [53]: where τ r y represents the longshore directed radiation stress (induced by the incident waves), the second term is the horizontal mixing due to cross-shore gradients in the longshore current velocity V, τ b y is the wave averaged bottom stress, and τ w y is the long-shore wind stress. At this point, it is important to mention that this modelling system was implemented and validated for all the target areas considered for evaluation [54][55][56]. Some additional details regarding the main characteristics and physical processes are presented in Table 2. In this table, ∆x and ∆y are the resolutions in the geographical space,∆θ is the resolution in the directional space, n f is the number of frequencies in the spectral space, nθ is the number of directions in the spectral space, ngx and ngy are the number of grid points in x and y direction, and np is the total number of grid points. The input fields considered in the computational domain are: wave forcing (wave), tide forcing (tide), wind forcing (wind) and currents fields (crt). The physical processes activated are: generation by wind (gen), whitecapping process (wcap), quadruplet nonlinear interactions (quad), triad nonlinear interactions (triad), diffraction process (dif ), bottom friction (bfric), wave-induced setup (setup) and activation of the depth-induced wave breaking (br).  Figure 2 illustrates the case studies considered. For each target area, a line of 3 km in length was considered, aiming to replicate the influence of a generic marine energy farm. A 2 km distance between the marine energy farm and the shoreline was considered, while the orientation of the farm was made according to the particularity of each target area.
The variations of the wave conditions in the presence of the marine energy farm will be assessed in the geographical space through spatial maps, while a deeper analysis of these fluctuations will be highlighted along the L-reference lines, or by the analysis performed in some offshore and nearshore points (O-points or NP-points). In order to provide a complete picture of the influence of an energy farm, two case studies were evaluated, as can be noticed from Table 3. The first one involved a realistic scenario where the wave farm was defined by a moderate absorption (denoted with M-farm) considering an absorption percentage of only 20% of the incoming waves. The other case study was related to a high absorption scenario (denoted with H-farm) and involved absorption of almost 40% of the waves, this being the case of a wave farm defined by several lines of WECs.

Leixoes (North Atlantic)
The spatial variations of the waves in the Leixoes sector are presented in Figure 3, considering only the waves coming from the north-western sector. As we go from the no farm situation to the high absorption scenario, it is clear that the impact of the wave farm is more visible, especially in the case of a storm, where multiple wave fields are noticed. The offshore points from this area present the following Hs values, that can range from~1.4 m (total average) to almost 5.62 m (storm conditions). For the total time period, the waves behind the marine energy farm may be reduced to a minimum of 0.5 m, while the harbour area appears not to be influenced by the presence of the farm. Regarding higher energy case studies, it is possible to notice that the direction of the waves represents an important factor for the coastal areas, the shielding effect of the farm being more visible in the lower part of this region. Table 4 presents the evolution of the wave parameters corresponding to the nearshore point group NP. Besides the significant wave height, some other parameters were considered, namely: (a) wave forces (in N/m 2 ); (b) Vbot (orbital velocity at the bottom in m/s). These indicators are used to identify the expected impact on the local seabed and also to assess what will be the expected dynamics of the sediment transport induced by the waves. Regarding the wave heights, the presence of the marine energy farm is more visible in the case of the high absorption scenario, especially in the case of the nearshore points NP2 and NP3, where the values indicate a minimum of 0.76 m in the case of the total average/H-farm. As for the values of the wave forces, there is a significant difference between the values corresponding to the point NP1 and those corresponding to the other points, which is obvious even when there is no wave farm. A marine energy farm can definitely reduce the initial forces, the expected impact gradually reducing as we go from total average to storm conditions. For example, in the case of the point NP2, the forces can be in the range of 1.05 to 0.41 N/m 2 (total average), while for the same point the values can decrease from 9.93 to 7.97 N/m 2 (storm).
The values of the parameter Vbot increase as the wave conditions become more energetic, with a minimum value of 0.24 m/s (total average-NP5) and a maximum of 2.5 m/s (storm-NP1) being noticed in the absence of the wave farm. More significant variations are related to the total time data, with a maximum difference of 0.3 m/s being expected between the total time/no farm and total time/H-farm (points NP2 and NP3). As we go to higher wave energy conditions (ex. storm/no farm), the variations are smaller, with no change being noticed for the point NP2 and an attenuation of 0.1 m/s in the cases of the NP3 and NP5 points.  Figure 4 and Table 5 present the variation of the Sines wave conditions in a geographical space and near the NP points. From the spatial distribution of the total average conditions, we notice that a marine energy farm defined by a high absorption property will have a visible impact. By looking at the values corresponding to the nearshore points, we notice that the point NP2 reveals a decrease from 1.33 to 1.18 m in the case of the total average/M-farm.
Regarding the other wave conditions, the spatial variations are similar to those corresponding to the Leixoes area. As for the nearshore points, more important variations correspond to the non-storm and storm conditions, being noticed a maximum attenuation of 1.34 m in the case of the point NP2. The variations of the forces is more significant near the points NP2, NP3 and NP4, while in the case of the parameter Vbot the site NP2 indicates in general lower values, regardless of the scenarios considered for assessment.

Porto Ferro (Mediterranean Sea)
Going from the North Atlantic Ocean to the Mediterranean Sea, in Figure 5 is presented the western part of the island of Sardinia, where a generic marine energy farm was considered north of the Porto Ferro inlet. The offshore points located in front of the marine energy farm indicate Hs values of 1.35 m-winter average, 2.4 m-non-storm, 3.55 m-storm or 5 m-high-storm, respectively. For the scenario winter average, the presence of a high absorption wave farm is reflected by a shielding effect that will keep the wave heights below 1 m. A similar pattern is noticed for the rest of the case studies, with the mention that for the high absorption scenario, an additional wave field occurs near the shoreline area. From the analysis of the spatial maps, we see that the wave farm has no impact on the wave conditions encountered near the Porto Ferro inlet (point NP4). Nevertheless, by looking at the values of the significant wave heights illustrated in Table 6a, we notice that the significant wave height presents some changes. For this point, the Hs values decrease from 1.09 (winter average/no farm) to 1.01 m (winter average/H-farm) or from 3.66 (high-storm/no farm) to 3.57 m (high storm/H-farm), respectively. The wave forces (Table 6b)

Scenario (a) Hs values (m)
Winter average  Another area considered for evaluation is the Saint George sector (north-west of the Black Sea, close to the Danube Delta); the Hs spatial distribution is presented in Figure 6. In this case, the coastline is located on the left side, while the dominant direction of the incoming waves is from the north-east  Moving to the nearshore area, the significant wave heights are influenced by the presence of a marine energy farm, as can be observed in Table 7a. For the scenario winter average, the reference points NP2, NP3 and NP4 indicate attenuation that reaches a maximum value of 0.3 m in the case of a high absorbing farm, while a variation of 0.2 m is noticed in the case of the moderate scenario. In the case of non-storm, a maximum difference of 0.54 m (point NP2) is noticed between the no farm situation and non-storm/H-farm, while for the same point the differences reach maximum values of 0.72 (storm) or 0.47 m (high-storm), respectively. The wave forces (Table 7b)

Scenario (a) Hs values (m)
Winter average Non-storm

Assessment of the Longshore Currents
It is well known that coastal areas are also influenced by the presence of the longshore currents [55,56]. These currents are generated by the breaking waves that enter in the surf area, it being expected that in the case of storm events the current velocity will indicate higher values. Figure 7 is focused on such an analysis that illustrates the longshore current distribution in the Leixoes sector, in this case indicating the maximum value of the current velocity (Vcmax in m/s) along the L-lines. From the analysis of these case studies, we can notice that various patterns can be identified.     The line L2 indicates no significant variation, regardless of the case study considered, while as we go to the scenarios storm and high-storm it appears that a marine energy farm will not have a big impact on the nearshore currents (excepting the line L3). In addition, it is important to mention that for the case studies winter average and non-storm, the velocity will increase along the line L1, reaching maximum values of 0.99 (winter average) or 1.37 m/s (non-storm), respectively. Figure 10 illustrates the maximum values of the current velocity along the reference lines, by considering this time the Saint George area. By looking at these results, we notice that mixed patterns occur for each case study. Thus, in the case of the line L1, the velocity corresponding to the case studies winter average and high-storm increases in the presence of a high absorber farm and decreases for a moderate one. For the scenario winter average/no farm, the current velocity decreases near the lines L3, L4 and L5, reaching minimum values of 0.07 m/s for winter average/M-farm and 0.01 m/s for winter average/H-farm. Since the wave direction is a crucial parameter in the development of the nearshore currents, the fact that the marine energy farm may produce significant changes in terms of wave direction may also affect the longshore current velocity. That is why, although the waves lose energy in the presence of the marine energy farm, in certain situations, due to such changes in wave direction, the longshore current velocity can increase down-wave from a marine energy farm. A similar situation was identified in [33], for various wave conditions.
For the non-storm situation, the current velocity will increase near the line L1, while an opposite trend is noticed close to the lines L4 and L5. Smaller fluctuations are noticed near the lines L2 and L3, while a similar situation is reported by the scenario storm (L2, L3, L5) or by high-storm (L3, L4, L5).

Discussion of the Results
The current approaches for coastal protection mainly involve seawalls and breakwaters. Since the aim of a WEC farm is to extract energy from the waves, a wave farm represents a very suitable alternative. Nevertheless, such farms are not beneficial in every environment, as is the case of enclosed seas, which do not represent the best option for the development of a wave project. This is because the wave conditions are significantly reduced compared to the ocean environment. On the other hand, Europe is an active player in the development of the WEC systems [57][58][59][60], and since a large part of this region is surrounded by semi-enclosed seas, there are higher chances to see marine farms operating in these waters in the future.
In any coastal area, the sediment transport is divided between bedload and suspended load, it being estimated that the most important hydraulic parameter for littoral transport is represented by the wave conditions. The littoral drift is influenced by several parameters (significant wave height; wave direction; grain size, etc.) which can be used to determine the littoral transport rate, denoted with Q (m 3 /24 hrs) [61]. For example, in the case of an incident angle of 30 • (like the one considered in the present work -see Table 1), the transport rate of a beach sand is associated with the following values: Hs = 1 m; Q = 300 m 3 /24 hrs; Hs = 3 m; Q = 10,000 m 3 /24 hrs; Hs = 5 m; Q = 65,000 m 3 /24 hrs. This assumption can be applied for most of the target areas, taking into account that the Saint George sector is defined by quartz sands (medium-fine sands), with a similar situation being noticed for the Sines area, where the local rivers represent the main source of sediments [33,38].
By interpolating the results presented in reference [61], it is possible to assess the transport rate corresponding to the nearshore points (from NP1 to NP5), assuming that the wave conditions are considered for a time interval of 24 h. Table 8 present the transport rates for some representative nearshore points. According to these results, we can notice that even small changes in the wave heights may lead to a significant reduction in the littoral transport rate generated by the wave action.
In this work, the longshore current velocity represents another parameter considered for investigation, which was assessed in each target area (from Figures 7-10). The Hjulström curve [62,63] is frequently used by hydrologists to determine if a river will transport/deposit sediment or will erode by taking into account the water velocity and the sediment particle size. If the water velocity is below 3 m/s, the sediment will be transported or deposited based on their size (<0.01 mm-transportation), while after this threshold the erosion processes may occur. At this point, it is important to mention that the Hjulström diagram is able to provide only a first-order analysis of the interaction between flowing water and sediments. From the analysis of the results in the Leixoes area, we notice that the erosion process may occur in the sector located close to the line L1, where in fact the presence of the farm increases the current velocity. For the scenarios total average, winter average and non-storm, a marine energy farm significantly reduces the current velocity for the lines L4 and L5, contributing in this way to the protection of the coastline. Regarding the Sines area, a marine energy farm may increase the erosion processes in the upper part of the shielded region (lines L1 and L2) and will significantly decrease the current velocity in the lower part (lines L4 and L5). As for the Porto Ferro area, by looking at the results corresponding to the high-storm scenario, we notice that near the line L3, the current velocity reduces. This does not necessarily mean that this sector will be protected, taking into account that the wave farm has no impact on the current velocity reported in the adjacent sectors. Regarding the expected values from the Saint George sector, we notice that they cannot be considered a source of coastal erosion, being assumed to be involved more in the sediment flow. For this sector, the coastal erosion will probably be more directly related to the wave action that can be generated during a storm event (storm and high-storm), with orbital velocities in the range of 1.3 and 1.5 m/s.

Conclusions
The objective of the present work is to identify the expected impact of a generic marine energy farm that would be implemented in various coastal environments, such as the Portuguese coast, Sardinia or in the north-western part of the Black Sea (close to the Danube Delta). The idea is to use the layout of a single marine energy farm, namely 3 km in length and located at about 2 km from the shoreline, in order to determine what will be the shoreline impact for a moderate and high absorption scenario.
According to these results, it was found that although from the analysis of the spatial maps there was reported no significant variation in the Hs values in the case of the first two scenarios (e.g., Portuguese coast-total average and winter average), according to the nearshore reference points, the influence of a WEC line is visible, and thus significant variations close to the shoreline are expected. From the analysis of the Hs values corresponding to the Sines and Leixoes areas, it was noticed that the Sines sector indicates in general many resources close to the shoreline, although the considered case studies are identical. A more complete picture of the nearshore impact was provided by including some other relevant parameters, such as the wave forces and the orbital velocity that may occur near the NP-points. Regardless of the considered target area, it was noticed that these two parameters increase as the wave conditions become more energetic. By looking at these results, we can say that, in general, the Leixoes area indicate much higher values, being followed by the Saint George sector in the Black Sea.
In terms of the current velocity, with the exception of the Saint George sector, it is possible that during an extreme event (e.g., storm), the longshore currents generated by the breaking waves affect the coastal stability by eroding the sediment deposits and transporting in other beach sectors the sediment transported by the waves. For the Saint George sector, it is more likely that the erosion processes are associated with to the hydraulic action of the waves and to the fact that the large dams build on the Danube River significantly reduces the volume of alluvia.
Since at this moment there are no operational wave farms, it is difficult to say what the configuration of such a project will be, and as a consequence throughout various ''what-if" case studies it is possible to estimate the expected implications for the coastal protection. Finally, we can mention that the results from the present work are in the line with the current research that considers hybrid modelling systems, which combine the output of a wave model with other simulation tools (ex: sediment transport; longshore currents).