Evaluating the Future E ﬃ ciency of Wave Energy Converters along the NW Coast of the Iberian Peninsula

: The e ﬃ ciency of wave energy converters (WECs) is generally evaluated in terms of historical wave conditions that do not necessarily represent the conditions that those devices will encounter when put into operation. The main objective of the study is to assess the historical and near future e ﬃ ciency and energy cost of two WECs (Aqua Buoy and Pelamis). A SWAN model was used to downscale the wave parameters along the NW coast of the Iberian Peninsula both for a historical period (1979–2005) and the near future (2026–2045) under the RCP 8.5 greenhouse scenario. The past and future e ﬃ ciency of both WECs were computed in terms of two parameters that capture the relationship between sea states and the WEC power matrices: the load factor and the capture width. The wave power resource and the electric power capacity of both the WECs will decrease in the near future. The load factor for Aqua Buoy will decrease in the entire area, while it will remain unchanged for Pelamis in most of the area, except north of 43.5 ◦ N. The capture width and cost of energy will increase for both devices. The methodology here applied can be easily applied to any device and coastal domain under di ﬀ erent climate change scenarios.


Introduction
The oceans are the largest source of unexploited renewable energy in the world. The development of ocean energy is one of the five strategic priorities of Europe to meet growing energy demand; to significantly improve energy security, reducing the dependence on energy imported from outside Europe; to support sustainable jobs and growth in the "blue economy"; to create a new internal market capable of exporting marine renewable energy products and services; and to reduce carbon emission from burning fossil fuels, improving also the air quality. Renewable energy located in marine environments has the advantage of larger areas for marine energy farms, with a higher and more stable power density. The limitations are deep waters for technological reasons, and that marine energy farms must coexist with other sea uses, like fishing and shipping, and with restricted or protected areas. However, ocean renewable resources may be exploited without harming the marine environment or negatively impacting the other sea uses if farms are sited and scaled properly following environmental guidelines [1]. In addition, future marine wave energy farms could play a key role in coastal protection [2][3][4][5][6][7] by extracting a considerable portion of wave energy, diminishing therefore the amount of energy in the surf zone, which attenuates the intensity of coastal processes, such as local erosion and coastline retreat.
In the marine environment, there are various energy resources that can be tapped using different technologies depending on the energy source. The main energy sources exploited worldwide are offshore wind, tidal-stream, and ocean waves. From these sources, waves are considered the most energetic one, since they are constantly formed and transport and retain energy for thousands of miles with minimum energy loss. Over the last decades, a huge effort has been devoted both to assess the main characteristics of the wave energy resource (a complete summary can be found in [8,9]) and to develop technologies allowing the sustainable extraction and exploitation of the wave energy [10][11][12][13][14][15][16]. Several wave energy converters (WECs) were developed during recent decades within almost 150 WEC projects reported at a global scale, 50% of them developed in Europe [11]. One of the main concerns is to determine the more suitable WEC for each coastal location, which depends on the wave characteristics usually found in that area. Attending to the directional characteristics, WECs can be classified as point absorber (pilot plants are set up in Portugal, Scotland, and Japan), attenuator (Portugal), and terminator (prototypes have already been installed in Denmark and Norway) [17].
Despite the enormous effort dedicated to analyze the wave energy resource both at global [18] and regional scales [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34] and to evaluate the adequacy of different WECs located in the worldwide regions with the highest wave power [10,11], little has been done to analyze the effect of climate change on wave power resources. The accurate prediction of the future energy field of a wave farm needs to take into account both the natural and anthropogenic variability of climate, since waves are quite vulnerable to changes in wind forcing [35,36]. In this sense, the impact of climate change on wave energy in the Black Sea basin was analyzed by [37] for the near future (2021-2050).
Furthermore, the wave energy harvest viability relies on the rate between the cost of the energy and the capital required for the projects, which may vary depending on the targeted market and the revenue [38][39][40][41][42] predicted at the present currency. Several authors have assessed the calculation of the life-cycle cost of the floating offshore wave energy farm taking into account different phases, such as development, manufacturing, installation, exploitation, and dismantling [39,40,43,44].
The northwest Iberian Peninsula coast (NWIP) is the target area of this study ( Figure 1a) and is one of the worldwide regions with the highest wave power resource, since waves break after traveling thousands of miles, gaining energy as they move through the Atlantic Ocean [11,16,45,46]. In fact, the first place worldwide to install a commercial-scale wave energy device in the water was the "Pilot Zone" located off the Portuguese coast just north of Lisbon. The "Pilot Zone" serves as a testing ground for ocean energy companies. In this region, the current [47,48] and future [49] wave energy resources were analyzed, and different evaluations of several WECs were performed to determine which technology is the most suitable for commercial exploitation in this area [10,11,14,47]. The area is affected by important bathymetric constraints (Figure 1a), since 80% of the points have a depth greater than 50 m and only 5.2% have depths between 25 and 50 m, which limits the type of WECs that can be used. Thence, the efficiency of the WECs will be analyzed for two types of devices whose technical specificities allow them to operate in intermediate water depths: a point absorber (Aqua Buoy) and an attenuator (Pelamis). A picture of both devices can be found in [50] (see Figure 2 in that publication). The aim of this study is to assess the efficiency and cost of energy over the life-cycle [41,42,51] of two WEC types along the northwest coast of the Iberian Peninsula for the near future. The efficiency of WECs in transforming wave energy into electricity is evaluated using two parameters, the load factor and the capture width, following previous studies [10,11,47]. The lifetime cost energy of the WECs can be estimated by the sum of the capital, operational, and maintenance costs. Other aspects, such as the size of the device or the design of the wave farm, were not considered in the present study. These parameters are calculated from the wave power resource and the wave power of each wave device. The wave power of each WEC is calculated, taking into account the historical and future percentages of sea states (structured into bins of Hs and Tp) and the power matrix and nominal power provided by the manufacturer. The wave states are obtained by means of a dynamical downscaling (~0.037° × 0.037°) of the wave height (Hs) and the peak period (Tp). Wave propagation models were driven by historical  wind data and near future (2026-2045) wind data produced following the Coupled Model Intercomparison Project 5 (CMIP5) protocol [52] under the Representative Concentration Pathway (RCP) 8.5 greenhouse scenario [53,54].

Methodology
Comparing the efficiency of different WECs is not straightforward due to the wide variety of technologies (point absorber, terminator or attenuator) for harvesting wave energy. Each WEC has a characteristic nominal power and a power matrix that represents the maximum energy that the WEC could harvest. This confers a range of power that each WEC can efficiently capture from the wave energy flux for different sea states. These sea states can be characterized in terms of the Hs and Tp parameters obtained from in situ measurements or numerical modeling for the past and only through numerical modeling for future projections.  assessment of the projected changes in terms of mean wave power was made by computing the difference ΔP in the wave power resource between both periods ( Figure 2). For the historical period, P increases northward, with the highest value (>65 kWm -1 ) observed at the northwestern corner of the IP (Figure 2a). In general, P gradually decreases landward. The pattern is similar to that previously obtained by other authors [69]. A similar trend is shown for the near future, although the maximum values do not exceed 60 kWm -1 (Figure 2b). ΔP is negative for the entire region, which suggests that P will decrease for the near future, although at different rates ( Figure 2c). In general, |ΔP| increases with the distance to the coast and from south to north. The decrease in P for the near future is consistent with previous research by [70], who found that despite the overall increase for the Atlantic basin since 1948, P has decreased over recent years for the Northeast Atlantic.  The P decrease can compromise the overall WECs efficiency and consequently lower the total wave power harvested in a wave farm located or planned for these locations. In addition, the average lifespan of WECs is assumed to be around 25 years [42,71] with the current technology, leading to a weighted investment due to the constraints for an optimized wave energy harvest. However, despite the P decrease obtained for the near future, the potential of this region for the wave harvest is still viable due to WEC technical specificities such as Pmax, which represents the maximum electric power that can be harvested from different sea states. Several authors have assessed the efficiency of WEC prototypes for commercial purposes in different coastal environments along the NWIP [47,66,69]. Some of the WECs assessed by those authors are Pelamis and Aqua Buoy, due to technical specificities that allow them to operate between intermediate water depths (50-70 m) [47], and the available The aim of this study is to assess the efficiency and cost of energy over the life-cycle [41,42,51] of two WEC types along the northwest coast of the Iberian Peninsula for the near future. The efficiency of WECs in transforming wave energy into electricity is evaluated using two parameters, the load factor and the capture width, following previous studies [10,11,47]. The lifetime cost energy of the WECs can be estimated by the sum of the capital, operational, and maintenance costs. Other aspects, such as the size of the device or the design of the wave farm, were not considered in the present study. These parameters are calculated from the wave power resource and the wave power of each wave device. The wave power of each WEC is calculated, taking into account the historical and future percentages of sea states (structured into bins of H s and T p ) and the power matrix and nominal power provided by the manufacturer. The wave states are obtained by means of a dynamical downscaling (~0.037 • × 0.037 • ) of the wave height (H s ) and the peak period (T p ). Wave propagation models were driven by historical  wind data and near future (2026-2045) wind data produced following the Coupled Model Intercomparison Project 5 (CMIP5) protocol [52] under the Representative Concentration Pathway (RCP) 8.5 greenhouse scenario [53,54].

Methodology
Comparing the efficiency of different WECs is not straightforward due to the wide variety of technologies (point absorber, terminator or attenuator) for harvesting wave energy. Each WEC has a characteristic nominal power and a power matrix that represents the maximum energy that the WEC could harvest. This confers a range of power that each WEC can efficiently capture from the wave energy flux for different sea states. These sea states can be characterized in terms of the H s and T p parameters obtained from in situ measurements or numerical modeling for the past and only through numerical modeling for future projections.

SWAN Wave Model
The wave parameters (H s and T p ) for historical and near-future climate data were obtained through the third-generation spectral wave model SWAN [55] embedded in the Delft3D WAVE module [56]. SWAN computes the evolution of random waves in coastal regions with deep, intermediate, and shallow water, which includes wave propagation, refraction due to depth, generation by wind, dissipation due to whitecapping, bottom friction, and depth-induced wave breaking and nonlinear wave-wave interactions [55]. The model was implemented for the northwest Iberian Peninsula coast, taking into account the recommendations of the International Electrotechnical Commission (IEC) for the development of a model to assess the wave energy resource [57]. The IEC set standards related to the modelling and analysis of wave resources depend on the objectives and accuracy of the study. These standards are divided into three categories: reconnaissance, feasibility, and design. The physical processes for a reconnaissance-stage model [58] were taken into account, including triads, bottom friction, depth-induced breaking, and quadruplets. The wave spectrum was discretized to 25 frequencies (0.0418-1 Hz) and 36 directional bands.
The dynamical downscaling considered in the present work will be briefly described here. In a first step, WAVEWATCH III (WWIII) coarse resolution simulations (1 • × 1 • ) provided by Commonwealth Scientific and Industrial Research Organisation (CSIRO) were analyzed. For hindcast purposes [59], CSIRO offers 8 different historical realizations forced with 8 different Global Circulation Models (GCMs) and an additional simulation forced with reanalysis winds obtained from the Climate Forecast System Reanalysis (CFSR). Both H s and T p were considered to assess the accuracy of the models to reproduce the wavefield in the area under study. The statistical analysis was carried out through the overlap percentage (OP) between the probability distribution functions (PDFs) of the wavefield generated with hindcast (CFSR) and GCMs winds, following [60]. The WWIII simulations driven with wind from the MIROC5 GCM showed to be the most accurate for the area under study, providing the highest OPs both for H s (87%) and T p (90%). Then, coarse (1 • × 1 • ) resolution wavefields were downscaled by means of the SWAN model at three different domains (1/3 • , 1/9 • , and 1/27 • , respectively) ( Figure 1b). SWAN was forced at the open boundaries with the wave field provided by WWIII forced with winds from the most accurate model (MIROC5 GCM). In addition, SWAN was also forced in the inner domain with regional winds provided by MIROC5 CCLM4-8-17 RCM [61]. The highest spatial resolution (~0.037 • × 0.037 • ) was used for a fringe~100 km width along the western coast of the Iberian Peninsula. The simulations covered the historical  and the near-future (2026-2045) periods with an hourly data output. The accuracy of the SWAN model in reproducing wave parameters near the coast was also statistically analyzed for the historical period in terms of the OP between in situ buoys ( Figure 1a) and numerical PDFs both for H s and T p , following [60], as mentioned above. Overlap values (Table 1) ranging from 78% to 83% were obtained for H s and from 74% to 84% for T p . The OPs between winds at coastal buoys moored around the Iberian Peninsula were analyzed in previous research [62], showing OP values of 69±10% for GCMs and 83±3% for RCMs.

Wave Power Resource
The Wave Power resource [63] (in kWm −1 ) can be obtained from the wave spectral parameters H s and T e , and is depicted by: where H s is the significant wave height, ρ is the seawater density (1025 kgm −3 ), g is the gravitational acceleration, and T e is the energy period, which can be estimated based on the peak period (T p ) as: where α depends on the shape of the wave spectrum. In the present study, α = 0.9 was adopted because a JONSWAP spectrum with a peak enhancement factor of γ = 3.3 was considered [64,65].

Parameters to Assess the Efficiency of WECs
Two parameters have been widely used in the literature to assess and compare the efficiency of different devices for specific sites, namely the load factor and the capture width [11,14].
The power load factor (ε) for a specific WEC can be expressed in % by: where P max represents the maximum electric power, provided by the manufacturer [47], that can be harvested with a specific WEC, and P E the average electric power that can be extracted with a certain device.
where p ij is the percentage of occurrence of a given sea state bined in terms of H s and T p , and P ij is the electric power for the same bin obtained from the power matrix provided by the manufacturer for that specific WEC [11,14]. The capture width (Cw, in meters) represents the width of the wavefront from which the energy is extracted, and it can be defined as the ratio of the WEC electric power (P E ) to the wave power resource (P) defined above: Energies 2020, 13, 3563 6 of 15 According to [66] these indices are the most appropriate to compare the efficiency of different WECs. Actually, they have been used to evaluate the performance of the different technologies along the Portuguese continental coastal environment [11].
The amount of electricity generated (E p , in MWh) [41] during the WEC lifetime can be expressed by: where Lifetime represents the total service life of the WEC (25 years). The cost of energy (C E , in EUR/MWh) [42] is the total cost involved during the lifespan use of the WEC, including the capital, operation, and maintenance cost, divided by E p : where C A is the capital cost and C OM is the operation and maintenance cost of a specific WEC. As mentioned before, two types of WECs are studied, Pelamis [42] and Aqua [67], whose sum of capital, operation, and maintenance costs can be estimated to be around 8.9 and 2.2 million euros, respectively, for a 25-year lifetime (see [42,68]). The difference between the future and the historical period for all the variables (P, ε, P E , Cw, and C E ) will be calculated as ∆X = X F − X H , where X is the variable under study and the subscripts F and H refer to the future and historical data, respectively.

Results and Discussion
The wave power energy was used to identify the wave power resource available along the NWIP coast both for the historical period  and the near future (2026-2045), and also the assessment of the projected changes in terms of mean wave power was made by computing the difference ∆P in the wave power resource between both periods (Figure 2). For the historical period, P increases northward, with the highest value (>65 kWm −1 ) observed at the northwestern corner of the IP (Figure 2a). In general, P gradually decreases landward. The pattern is similar to that previously obtained by other authors [69]. A similar trend is shown for the near future, although the maximum values do not exceed 60 kWm −1 (Figure 2b). ∆P is negative for the entire region, which suggests that P will decrease for the near future, although at different rates ( Figure 2c). In general, |∆P| increases with the distance to the coast and from south to north. The decrease in P for the near future is consistent with previous research by [70], who found that despite the overall increase for the Atlantic basin since 1948, P has decreased over recent years for the Northeast Atlantic.
The P decrease can compromise the overall WECs efficiency and consequently lower the total wave power harvested in a wave farm located or planned for these locations. In addition, the average lifespan of WECs is assumed to be around 25 years [42,71] with the current technology, leading to a weighted investment due to the constraints for an optimized wave energy harvest. However, despite the P decrease obtained for the near future, the potential of this region for the wave harvest is still viable due to WEC technical specificities such as P max , which represents the maximum electric power that can be harvested from different sea states. Several authors have assessed the efficiency of WEC prototypes for commercial purposes in different coastal environments along the NWIP [47,66,69]. Some of the WECs assessed by those authors are Pelamis and Aqua Buoy, due to technical specificities that allow them to operate between intermediate water depths (50-70 m) [47], and the available information of the power matrix [14] and the corresponding P max are 750 kW for Pelamis and 250 kW for Aqua Buoy [47].
Thus, the Pelamis and Aqua Buoy electric power capacity (P E ) was calculated based on their power matrix and the sea states previously obtained from the SWAN simulations for the historical and near-future periods (Figures 3 and 4). Both figures show the historical values, the near-future projections, and the change in P E between both periods for Pelamis and Aqua Buoy, respectively. The results show that the P E patterns for Pelamis and Aqua Buoy are similar to those previously found for P. In addition, the decreasing pattern previously observed for ∆P is also detected for ∆P E . Pelamis shows a maximum ∆P E reduction (>10 kW in absolute value) at ocean locations in the northern part of the study area (Figure 3c). The ∆P E values are almost negligible from Roca Cape to 43 • N. Only a narrow fringe north of 43.5 • N shows a decrease higher than 4 kW in absolute value. Aqua Buoy shows a maximum ∆P E reduction (~−7 kW) also in the northern part (Figure 4c). ∆P E will decrease between −2.5 and −4 kW south of 43.5 • N, reaching the lowest differences close to the shore between Roca Cape and São Vicente Cape. Overall, both devices show a different ∆P E pattern south of 43.5 • N, with the pattern similar north of 43.5 • N. These patterns are explained by the differences in the power matrices of the WECs, whose diagrams of scattered T p -H s joint distributions have different ranges, between 1-5.5 s and 5-17 m for Aqua Buoy and 0.5-8 s and 5-13 m for Pelamis [14]. These varieties of ranges confer each device a diverse variety of sea states under which they can operate retrieving the maximum electric power energy. The performance of WECs may vary depending on the topographic constraints and sea states, The performance of WECs may vary depending on the topographic constraints and sea states, in such a way that just a portion of the total energy is converted into electric power. The power load factor (ε) was calculated for both devices to assess the total energy (in %) harvested under the different sea states. Figures 5 and 6 show the power load factor for Pelamis and Aqua Buoy for the historical period (a) and the near future (b), and (c) the difference between both periods (∆ε). The power load factor for Pelamis and Aqua Buoy presents maximum values (25% and 30%, respectively) at ocean locations north of 43.5 • N and low values nearshore both for historical and near-future periods (Figure 5a,b and Figure 6a,b). Although the ε pattern is similar for both periods, differences can be highlighted in terms of the ∆ε parameter (Figures 5c and 6c). The Pelamis power load factor will remain practically unchanged (<0.5% in absolute value) (Figure 5c), with the exception of the region north of 43.5 • N, where reductions between −1% and −1.5% are expected. The Aqua Buoy power load factor will decrease between −1% and −2% across the region, as shown in Figure 5c.
The second indicator of the efficiency of a WEC is the capture width (Cw), which represents the width of wave front where the energy is extracted by the WEC [29,43]. Low values correspond to a low efficiency of transforming the wave energy into electricity. The Cw of both wave devices shows higher values near the coast than at ocean locations both for the historical and future time periods. In addition, they also tend to increase southward. Cw values between 2.4 and 2.8 m were obtained for Pelamis (Figure 7a) and between 1.2 and 1.4 m for Aqua buoy (Figure 8a) during the historical period. The highest Cw values, 3 m for Pelamis and 1.5 m for Aqua Buoy, were obtained for a small region located on the north coast of Spain. Overall, the capture width will increase for the near future ( Figures  7b and 8b), showing a pattern similar to the one observed for the historical period but with higher Cw values, especially near the shore. The ∆Cw ranges from 0.2 to 0.25 m between Roca Cape and Peña Cape for Pelamis (Figure 7c). The highest ∆Cw is detected in Roca Cape, with values higher than 0.3 m. Aqua Buoy shows the same ∆Cw pattern as Pelamis, although with a more modest increment (between 0.04 and 0.06 m). Aqua Buoy also shows a maximum ∆Cw (~0.07 m) near Cape Roca (Figure 8c).  The performance of WECs may vary depending on the topographic constraints and sea states, in such a way that just a portion of the total energy is converted into electric power. The power load factor (ε) was calculated for both devices to assess the total energy (in %) harvested under the different sea states. Figures 5 and 6 show the power load factor for Pelamis and Aqua Buoy for the historical period (a) and the near future (b), and (c) the difference between both periods (Δε). The power load factor for Pelamis and Aqua Buoy presents maximum values (25% and 30%, respectively) at ocean locations north of 43.5 ° N and low values nearshore both for historical and near-future periods (Figures 5a,b and 6a,b). Although the ε pattern is similar for both periods, differences can be highlighted in terms of the Δε parameter (Figures 5c and 6c). The Pelamis power load factor will remain practically unchanged (<0.5% in absolute value) (Figure 5c), with the exception of the region north of 43.5 ° N, where reductions between −1% and −1.5% are expected. The Aqua Buoy power load factor will decrease between −1% and −2% across the region, as shown in Figure 5c.   The second indicator of the efficiency of a WEC is the capture width (Cw), which represents the width of wave front where the energy is extracted by the WEC [29,43]. Low values correspond to a low efficiency of transforming the wave energy into electricity. The Cw of both wave devices shows higher values near the coast than at ocean locations both for the historical and future time periods. In addition, they also tend to increase southward. Cw values between 2.4 and 2.8 m were obtained for Pelamis ( Figure 7a) and between 1.2 and 1.4 m for Aqua buoy (Figure 8a) during the historical period. The highest Cw values, 3 m for Pelamis and 1.5 m for Aqua Buoy, were obtained for a small region located on the north coast of Spain. Overall, the capture width will increase for the near future (Figures 7b and 8b), showing a pattern similar to the one observed for the historical period but with higher Cw values, especially near the shore. The ΔCw ranges from 0.2 to 0.25 m between Roca Cape and Peña Cape for Pelamis (Figure 7c). The highest ΔCw is detected in Roca Cape, with values higher than 0.3 m. Aqua Buoy shows the same ΔCw pattern as Pelamis, although with a more modest increment (between 0.04 and 0.06 m). Aqua Buoy also shows a maximum ΔCw (~0.07 m) near Cape Roca (Figure 8c).   As far as we know, little is known about the energy cost of WECs under future climate conditions. Here, the cost of energy for the near future will be calculated by considering the current capital (C A ), operation, and maintenance (C OM ) costs and the generated energy during a life-cycle of 25 years, as described above (Equations (6) and (7)). Figures 9 and 10 show the estimated cost of energy (C E ) for Pelamis and Aqua Buoy, respectively. Costs are calculated for the historical period (a) and the near future (b). In addition, the cost increase between both periods is depicted in frame (c). The C E patterns are similar for both WECs, with higher values confined to the coast and lower values offshore. In general, the C E is higher for Pelamis (400-1000 EUR/MWh) than for Aqua Buoy (120-300 EUR/MWh) due to the lower production, operation, and maintenance costs. Pelamis shows a maximum ∆C E increase (>20 EUR/MWh) south of Cape Roca, between 42 • N and Cape Finisterre, and in the northeast area (Figure 9c). Only in the ocean section between 39 • N and 43 • N is the cost of energy the same as at present. The Aqua Buoy ∆C E also increases in the same areas, with the exception of the northeast region. As far as we know, little is known about the energy cost of WECs under future climate conditions. Here, the cost of energy for the near future will be calculated by considering the current capital (CA), operation, and maintenance (COM) costs and the generated energy during a life-cycle of 25 years, as described above (Equations (6) and (7)). Figures 9 and 10 show the estimated cost of energy (CE) for Pelamis and Aqua Buoy, respectively. Costs are calculated for the historical period (a) and the near future (b). In addition, the cost increase between both periods is depicted in frame (c). The CE patterns are similar for both WECs, with higher values confined to the coast and lower values offshore. In general, the CE is higher for Pelamis (400−1000 EUR/MWh) than for Aqua Buoy (120−300 EUR/MWh) due to the lower production, operation, and maintenance costs. Pelamis shows a maximum ΔCE increase (>20 EUR/MWh) south of Cape Roca, between 42° N and Cape Finisterre, and in the northeast area (Figure 9c). Only in the ocean section between 39° N and 43° N is the cost of energy the same as at present. The Aqua Buoy ΔCE also increases in the same areas, with the exception of the northeast region.

Conclusions
The main novelty of this study is that the electric power capacity, the efficiency, and the cost of energy of WECs were calculated both for the present and near future. The calculations were performed for the entire domain (from 37.01° N to 44.43° N and from 7°W to 9.88°W) and not only at certain points, as carried out in previous studies in the same area. The SWAN model was used to provide high spatial resolution (~0.037° × 0.037°) in a fringe ~100 km-wide along the west coast of Iberia, resulting in around 7000 points that constitute possible locations to install WEC farms. The study provides then a very detailed picture on the most suitable locations to moor a particular device both at present and in the near future. The main conclusions on the impact of projected changes on the efficiency of WECs along the west coast of Iberia for the near future under the RCP 8.5 can be summarized as follows: First, the study highlighted a decrease in the wave power resource for the near future throughout the domain, which is more intense in the northern corner of the Iberian Peninsula, where the highest wave power resource was detected both for the historical and the future period.
The electric power capacity follows the same spatial pattern as the wave power resource across the region both for the Pelamis and Aqua Buoy devices. The maximum decrease in electric power for the near future was detected at ocean locations in the northern part both for Pelamis (decrease > 10kW) and Aqua Buoy (decrease ~7 kW). South of 43° N, Pelamis shows a negligible decrease, while Aqua Buoy shows a decrease ranging from 2.4 kW to 4.5 kW, which is the lowest decrease observed at southern locations.
The power load factor for Pelamis will remain practically unchanged (< 0.5% in absolute value) throughout the area, except for the region north of 43.5° N, where reductions between 1% and 1.5% are expected. The power load factor for the Aqua Buoy will decrease between 1% and 2% across the region.
The capture width will increase from 0.2 to 0.25 m between Cape Roca and Cape Peña for Pelamis. The greatest increase is detected in Cape Roca, with values greater than 0.3 m. Aqua Buoy shows the same spatial distribution for the capture width as Pelamis, although with a more modest increase (between 0.04 and 0.06 m). Aqua Buoy also will show a maximum increase (~0.07 m) near Cape Roca.
The cost of energy was projected to increase in the near future for the entire region for both devices mainly due to the decrease in the electric power (PE). A maximum increase nearshore south

Conclusions
The main novelty of this study is that the electric power capacity, the efficiency, and the cost of energy of WECs were calculated both for the present and near future. The calculations were performed for the entire domain (from 37.01 • N to 44.43 • N and from 7 • W to 9.88 • W) and not only at certain points, as carried out in previous studies in the same area. The SWAN model was used to provide high spatial resolution (~0.037 • × 0.037 • ) in a fringe~100 km-wide along the west coast of Iberia, resulting in around 7000 points that constitute possible locations to install WEC farms. The study provides then a very detailed picture on the most suitable locations to moor a particular device both at present and in the near future. The main conclusions on the impact of projected changes on the efficiency of WECs along the west coast of Iberia for the near future under the RCP 8.5 can be summarized as follows: First, the study highlighted a decrease in the wave power resource for the near future throughout the domain, which is more intense in the northern corner of the Iberian Peninsula, where the highest wave power resource was detected both for the historical and the future period.
The electric power capacity follows the same spatial pattern as the wave power resource across the region both for the Pelamis and Aqua Buoy devices. The maximum decrease in electric power for the near future was detected at ocean locations in the northern part both for Pelamis (decrease > 10 kW) and Aqua Buoy (decrease~7 kW). South of 43 • N, Pelamis shows a negligible decrease, while Aqua Buoy shows a decrease ranging from 2.4 kW to 4.5 kW, which is the lowest decrease observed at southern locations.
The power load factor for Pelamis will remain practically unchanged (<0.5% in absolute value) throughout the area, except for the region north of 43.5 • N, where reductions between 1% and 1.5% are expected. The power load factor for the Aqua Buoy will decrease between 1% and 2% across the region.
The capture width will increase from 0.2 to 0.25 m between Cape Roca and Cape Peña for Pelamis. The greatest increase is detected in Cape Roca, with values greater than 0.3 m. Aqua Buoy shows the same spatial distribution for the capture width as Pelamis, although with a more modest increase (between 0.04 and 0.06 m). Aqua Buoy also will show a maximum increase (~0.07 m) near Cape Roca.
The cost of energy was projected to increase in the near future for the entire region for both devices mainly due to the decrease in the electric power (P E ). A maximum increase nearshore south of Cape Roca and Cape Finisterre was projected for both devices. In addition, a noticeable future increase in C E is also observed in the northeast area for Pelamis.
The results have shown that the efficiency of the different devices is not only dependent on the wave resource at a certain study area but also on the relationship between the sea states and the power matrix of the WECs that can operate in that area. The projected changes in wave parameters that affect wave resources can be magnified or softened depending on the particular features of the device. The electric power capacity and its derived parameters (the power load factor and the capture width) have shown to be the pertinent variables to be considered to analyze the current and near future efficiency of the different converters.
In summary, this work provides a valuable methodology to analyze the projected changes in efficiency for different WECs along any coastal domain under different climate change scenarios. In addition, the estimation of energy cost in the future for different devices constitutes an added value of the study that allows assessing the economic viability of wave farms under future climate conditions. Funding: The first author of this work has been supported by the Portuguese Science Foundation (FCT) through a doctoral grant (SFRH/BD/114919/2016). Thanks are due to FCT/MCTES for the financial support to CESAM (UIDB/50017/2020+UIDP/50017/2020), through national funds. This work was partially supported by Xunta de Galicia under project ED431C 2017/64-GRC (Grupos de Referencia Competitiva) and by Ministry of Economy and Competitiveness of the Government of Spain under the project "WELCOME ENE2016-75074-C2-1-R" funded by European Regional Development Fund (ERDF). This work was also partially supported by the project CLIMENA-CLimate change Impacts on the Marine Environment of the North Atlantic-funded by the Portuguese Foundation for Science and Technology (FCT-Fundação Portuguesa para a Ciência e a Tecnologia) under contract n. • : PTDC/EAM_OCE/28561/2017. The present study is also part of the project "WECAnet: A pan-European network for Marine Renewable Energy" (CA17105), which received funding from the HORIZON2020 Framework Programme by COST (European Cooperation in Science and Technology), a funding agency for research and innovation networks.