Assessment of Hybrid Wind-Wave Energy Resource for the NW Coast of Iberian Peninsula in a Climate Change Context

Offshore renewable energy has a high potential for ensuring the successful implementation of the European decarbonization agenda planned for the near future. Hybrid wind-wave farms can reduce installation and maintenance costs, and increase the renewable energy availability of a location by compensating for the wind’s intermittent nature with good wave conditions. In addition, wave farms can provide protection to wind farms, and the combined wind/wave farm can provide coastal protection. This work aims to assess the future hybrid wind-wave energy resource for the northwest coast of Iberian Peninsula for the near future (2026–2045), under the RCP 8.5 greenhouse gas emission scenario. This assessment was accomplished by applying a Delphi classification method to define four categories, aiming to evaluate the richness (wind and wave energy availability, downtime), the variability (temporal variation), the environmental risk (extreme events), and cost parameters (water depth and distance to coast) of the wind and wave resources. The combined index (CI), which classifies the hybrid wind-wave resource, shows that most of the NW Iberian Peninsula presents good conditions (CI > 0.6) for exploiting energy from wind and wave resources simultaneously. Additionally, there are some particularly optimal areas (CI > 0.7), such as the region near Cape Roca, and the Galician coast.


Introduction
The development of marine energies has become essential for the green energy sector, in order to mitigate the human influence on climate change, and move towards 100% clean energy, as stated in the 2030 Agenda for Sustainable Development released by the UN in 2015 [1]. Pursuing this goal, both the European Parliament and the European Council have endorsed 60% emission reductions by 2030 [2], and long-term climate neutrality for the EU by 2050.
Europe's decarbonization plan to achieve the net-zero greenhouse gas emission target relies heavily on an offshore renewable energy strategy [3]. The strategic roadmap is to deliver up to 100 GW of capacity by 2050, which represents approximately 10% of Europe's electricity consumption [3,4]. The European Commission, through its member states, has made great efforts to achieve this goal by promoting projects that address ocean energy efficiency [3][4][5][6][7], and the competitiveness of offshore wind [8][9][10] and ocean energies [11,12]. In fact, the European Union is currently the world leader in The aim of this study was to assess the present and near future wind and wave energy resources along the Atlantic coast of the Iberian Peninsula. A Delphi classification method [48] that takes into account wind and wave resources, stability, risk, and installation and maintenance costs was applied for the analysis. The originality of the present study relies not only on the evaluation of the combined resources using a Delphi, as done in previous studies, but also in the consideration of other factors besides the wind/wave power resources, to assess the feasibility of a hybrid wind-wave energy farm at a particular location for future exploitation conditions. This methodology contributes to the knowledge and applicability for global-to-local scale studies, and predicts future wind-wave energy resource for researchers and ocean energy entrepreneurs.
The paper is organized as follows: in Section 2, the data is presented, followed by the methodology applied to assess energy resources. Results and discussion can be found in Section 3. Finally, conclusions are presented in the Section 4.

Wind Data
Daily wind speed at 10 m height was retrieved from regional climate simulations within the CORDEX project [36] (Table 1) under the RCP8.5 greenhouse emission scenario. RCMs at a spatial resolution of 0.11° were considered from 2026 to 2045.  The paper is organized as follows: in Section 2, the data is presented, followed by the methodology applied to assess energy resources. Results and discussion can be found in Section 3. Finally, conclusions are presented in the Section 4.

Wind Data
Daily wind speed at 10 m height was retrieved from regional climate simulations within the CORDEX project [36] (Table 1) under the RCP8.5 greenhouse emission scenario. RCMs at a spatial resolution of 0.11 • were considered from 2026 to 2045. Table 1. Regional climate models (RCMs) from the Coordinated Regional Climate Downscaling Experiment (CORDEX) project (http://www.cordex.org/) used in this study. GCM = global climate models.

GCM RCM INSTITUTE
CNRM-CM5 CCLM4- 8-17  CLMcom  CNRM-CM5  RCA4  SMHI  EC-EARTH  RACMO22E  KNMI  IPSL-CM5A-MR  RCA4  SMHI  IPSL-CM5A-MR  WRF331F  INERIS  MPI-ESM-LR  CCLM4-8-17  CLMcom  MPI-ESM-LR  REMO2009  MPI-CSC A multi-model ensemble approach was applied by using a consensus criterion [49,50] in order to reduce biases of each of the models that make up the multi-model. This consensus criterion determines the concordance of the multi-model ensemble wind speed difference (future minus past) with the different models, by imposing two conditions to evaluate the statistical significance of the pixels. The first condition states that at least 75% of all models present the same sign of the multi-model mean in every pixel. The second condition imposes that at least 80% of the models that met the first condition pass the Mann-Whitney-Wilcoxon test [51] (a 5% significant level was considered at each pixel). For further information on the consensus criterion, the reader is referred to [49,50]. The accuracy of this ensemble of RCMs was previously analyzed in the area under scope [52,53].

Wave Data
Wave data parameters (H s and T p ) for the future period, 2026-2045, were obtained through the simulating waves nearshore (SWAN) model [54]. The model was previously implemented and validated for the northwest Iberian Peninsula coast [30]. The simulations consider the recommendations of the International Electrotechnical Commission (IEC) for a reconnaissance-stage model [55], which are required for an energy resource assessment.
The SWAN simulations are the dynamical downscaling (0.11 • × 0.11 • ) of the MIROC5 GCM. The Commonwealth Scientific and Industrial Research Organisation (CSIRO) database for wave hindcast [56] and wind-wave climate projections [57] offers eight different GCMs simulations. These simulations were statistically analyzed through the overlap percentage between the probability distribution functions (PDFs) [58] of the waves hindcast and the GCMs projections. The MIROC5 was found to be the most accurate model [30]. The time required to perform each of the SWAN simulations limits the use of only one GCM out of the eight possible. Winds provided by MIROC5 CCLM4-8-17 RCM [36] were also used in the SWAN simulations.

Methods
The synergy between the wind and the waves can optimize the harvest of marine resources in the same area using hybrid farms. Thus, when assessing the energy resource, focusing on the harvest through hybrid systems, various factors must be taken into consideration, both specific to each resource and common to both. The assessment applied in this study is based on different indices representing the resource richness, variability, risk, and cost of wind energy resource, previously used by Costoya et al. [49,50], and applied in a similar way for the wave energy resource.

Wind Power Resource
The wind power density (WPD in watts m −2 ) was determined by: where ρ a is the air density (1.225 kgm −3 at 15 • C and 1000 hPa) and W H is the wind speed at the selected hub height (120 m). WPD is the power available to be converted by a wind turbine at a specific location. An important detail, when investigating wind energy changes due to global warming, is related to changes in the air density due to temperature increases. However, these are found not be relevant since a 5 k temperature increase represents only a 1% decrease of the WPD (for the same wind speed).
Since WPD is independent of the turbine power curve, it is a suitable wind energy metric for this study.

Wave Power Resource
Wave energy resource was evaluated by means of the wave power (WP in kWm −1 ) following the expression: where H s is the significant wave height, ρ is seawater density (1025 kgm −3 ), g is the gravitational acceleration, and T e is the energy period estimated on the peak period (T p ) [30,59]. Harvesting energy from waves is not straightforward, as each wave energy converter has a characteristic power matrix that varies with its design. In this sense, wave power can be applied to any wave energy converter since it is independent of the power matrix.

Assessment of Future Wind Energy Resource
Eight indices were considered to assess the future wind energy resource in the northwest Iberian Peninsula coast, following the methodology applied by [38,49,50] to assess the richness of the wind resource, the variability, the environmental risk, and cost parameters. These indices are: Annual average wind speed (W ann ), rich level occurrence (RLO), downtime (DWNT), coefficient of variation index (C v ), monthly variability index (M v ), risk of extreme wind speeds (Risk), water depth (WD), and distance to coast (DC).
The richness of the wind resource parameter accounts for the potential wind energy resource available, which is independent of the converter device. The W ann index is divided into seven partitions following the NREL classification scheme [60]. The rich level occurrence measured by means of the RLO index is defined as the frequency of WPD higher than 200 Wm −2 .
The effectiveness in the exploitation of the wind energy resource is also a key factor because wind turbines only operate in a particular range that allows the production of energy. The upper (cut-off speed) and lower (cut-in speed) thresholds of this range are usually around 4 and 25 ms −1 , respectively [49,50]. In this sense, the DWNT index provides a measure of the total time with useful wind speed occurrence (EWSO) in which a wind turbine is producing electricity.
The temporal variability of the wind energy resource is another important factor, since a stable energy supply throughout the year would optimize the energy harvest efficiency. In this context, the coefficient of variation (C v ) index is defined as: where S and x are the standard deviation and the WPD mean, respectively. Additionally, the monthly variability (M v ) index was also considered. This index is defined as: where P year is the annual mean of the WPD, and P M1 and P M12 represent the WPD for the most and least energetic months, respectively. The risk index was also considered, to take into account extreme wind speeds (EWS) that can highly affect the safety of ocean engineering. This index was calculated by means of the EWS over a 50-year period calculated based on the Gumbel curve method.
Finally, water depth and distance to coast indices were also considered for every grid pixel. These indices represent the cost factors related to grid connection and marine engineering. For this purpose, the ETOPO bathymetry and the Global Self-consistent, Hierarchical, High-resolution Geography (GSHHG) coastline databases were used to calculate WD and DC indices, respectively.
The assessment of the future wind resource based on the indices described above was carried out following a Delphi classification method similar to the one described in [37,49,61]. First, each index was normalized to encompass all of them in a single parameter (the C wind index) since, as mentioned above, each index has its magnitude and units. The normalization process followed is similar to the one applied by [49,50], in such a way that both positive and negative indices were turned into positive factors with the optimal value equal to 1 and the worst equal to 0. The partitions of the [0, 1] interval were calculated considering a 1/(n − 1) step, where n is the number of categories (10 or 5, depending on the index). Most of the indices were normalized using 10 partitions (Table 2) ,with the exception of the W ann index normalized using 7 (as mentioned before), and WD and DC which were normalized using 5 partitions (Table 3). Table 2. Normalization criterion used for rich level occurrence (RLO), downtime (DWNT), coefficient of variation index (C v ), monthly variability index (M v ), and extreme wind speeds (EWS) indices related to richness factors, variability factors, and environmental risk factors.

WD (m) DC (deg)
The weight coefficients for wind energy indices based on the Delphi classification method are described in Table 4. Table 4. Weight coefficient in the wind energy classification. The standard deviation was calculated /N, where A is the weight coefficient, A is the mean value of A and N is the number of experts (10) consulted in [37].

Assessment of Future Wave Energy Resource
The assessment of the future wave energy resource was similar to the one used for the future wind resource, since both wind and wave resources share the same characteristics and constraints. Previous studies assessed the wave energy resource [38,62,63], taking into account several parameters which evaluated the energy production (linked to a specific WEC), extreme events, availability of the optimal wave height, accessibility, and monthly variability of the resource. Future wave energy resource was classified following a method similar to previous studies [38,49], based on a Delphi technique. In particular, experts were selected from the steering committee of, and participants in, the COST Action WECANet: A pan-European Network for Marine Renewable Energy (https://www.wecanet.eu/). As a first step, experts were questioned about the factors that they considered to be influential for wave energy converters. Experts were consulted only about general factors that did not depend on the particular features of a particular type of device. These factors, which include drivers that depend on the resource (mainly the amount of wave energy and its variability), on the environmental risk, and on cost constraints, were proposed and agreed by the experts. In a second step, a normalization method was adopted to compare different factors. Finally, as the last step, experts were consulted again to determine the weight of each factor following the protocol described latter.
Similarly as previously done for wind, the C wave index was calculated considering seven indices that evaluated the wave resource richness, variability, environmental risk, and cost parameters: wave energy resource (WP), downtime (DWNT), coefficient of variation index (C v ), monthly variability index (M v ), risk of extreme wave heights (Risk), water depth (WD), and distance to coast (DC).
The richness of the wave resource parameter (WP) includes the potential wave energy resource available in the area under study, and is independent of the wave energy converter. Like for wind energy converters, wave converters also operate on a particular range of T p and H s . In this context, and following [64], the DWNT index was defined with upper and lower thresholds set to H s = 8 m and H s = 1 m, which represent the level of "power outages" due to extreme wave conditions, and the time of calms, respectively.
The temporal variability of the wave energy was also considered for the wave energy exploitation by means of the C v and M v indices that were calculated as for the wind resource, but changing the WPD for WP.
Similar to EWS, a risk factor that takes into account extreme wave heights (EWH) was also considered. This index was calculated by means of the maximum H s over a 50-year period, calculated based on the Gumbel curve method.
The calculation of the cost factors (WD and DC) was performed for every grid pixel, using the same method applied for the wind resource, and the same normalization process described in Table 3.
The other indices were also subjected to the same normalization process applied in the wind energy classification. The normalization criterion is described in Table 5. The weight coefficients for wave energy indices based on the Delphi classification method are described in Table 6. In particular, experts were selected from the steering committee and participants in the COST Action WECANet: A pan-European Network for Marine Renewable Energy (https://www.cost.eu/actions/CA17105/). As a first step, experts were questioned about the factors that they considered to be influential for wave energy converters. Experts were consulted only about general factors that did not depend on the particular features of a particular type of device. These factors, which include drivers that depend on the resource (mainly the amount of wave energy and its variability), on the environmental risk, and on cost constraints, were proposed and agreed by the experts. In a second step, a normalization method was adopted to compare different factors. Finally, in the last step, experts were consulted again to determine the weight of each factor following the protocol previously described.
where A is the weight coefficient, A is the mean value of A and N is the number of experts (15). In this case, the indices that account for the wave resource richness are only two, compared to the three for the wind resource. This is because the wave energy harvest is usually done in terms of the average WP [65,66], in contrast to the wind resource that is usually done in terms of the W ann and EWSO.

Richness
Similar to the wind energy classification index, the wave energy classification index (C wave ) was calculated for each grid point as: where I i and W i are the indices and weights described in Table 6. Finally, the classification of the hybrid wind-wave energy resource was determined by a combined index (CI) of wind and wave energy resources defined as:

Sensitivity Test of the Energy Resource Parameters
The classification indices used in Equations (5) and (6) were calculated for each grid pixel (k) by means of an expression that can be generalized as: where N is the number of indices, I k i is the normalized value of each index (i) calculated for every pixel (k), and W k i is the mean weight of each index obtained from Tables 2-6. Note that each weight in Tables 4 and 6 is described as a mean and a standard deviation. Although the mean weight value of the weight coefficients ( W k i ) will be considered throughout most of the study (Equations (5) and (6)), a sensitivity test was carried out following a Monte Carlo approach to analyze the dependence of the classification indices (C k ) on the variability of the weight coefficients. The procedure applied to each pixel can be summarized as follows: 1. Random variation. Each weight was allowed to vary randomly where R k is a random number with equal probability between ±1, and σ is the standard deviation provided by Tables 4 and 6 for wind and waves, respectively. A different random number was used for each pixel (k). 2. Renormalization. Before the random variation, weights were normalized in such a way that However, that normalization cannot be guaranteed after random variation. Consequently, weights must be renormalized by 3. The classification indices were recalculated using the new weights where the subscript (j) refers to the different realizations of the Monte Carlo approach. In particular, j = 100,000 realizations will be considered in the present analysis.
4. The RMSE k associated to each pixel was calculated between the reference index (C k ) and the different realizations, C k j .
5. In a similar way, the RMSE j associated to each realization was calculated as where N p is the number of pixels.
6. The mean value of the classification index was calculated as the average of all realizations by using the following expression

Results and Discussion
The assessment of the potential energy resource for a hybrid wind-wave farm was carried out along the northwest coast of the Iberian Peninsula for the near future (2026-2045). First, spatial maps of all the wind and wave indices that make up the categories of resource richness, variability, risk, and cost were calculated.
The method used in this study is similar to the one developed by [62], where the authors proposed a multi-criteria approach for wave energy converters. They used similar coefficients (resource, variability, downtime, risk, and cost) although their study was more focused on particular devices for three specific locations. In the case of this research, which is based on previous studies carried out for wind energy converters [37,49,50,61], the analysis is focused on a large area covering the Atlantic part of the Iberian Peninsula under near future conditions. Previous works also focused on particular locations [25,65], although as far as we know, none on them considered future wind and wave conditions. The richness category accounts for the available resource in terms of W ann , RLO, and DWNT for the wind resource (Figure 2), and WP and DWNT for the wave resource ( Figure 3). W ann index (Figure 2a) shows the maximum value of 1 in the ocean, and decreases to less than half near land. This is consistent with the WPD values shown for this region in [67]. The RLO index pattern is similar to the one of W ann , with the lowest values (around 0.1) being detected along the coast, and the maximum values, of around 0.7, in the northwest region. The DWNT index shows a smaller difference (0.2) between the ocean and the coastal regions. Wave richness resource shows some similarities to that of the wind, mostly due to the waves being generated as a result of the wind forcing on the water surface. WP index (Figure 3a   The indices that account for the wind power resource variability, Cv and Mv, are represented in Figure 4a,b, respectively. In general, the highest Cv values of 0.5 are in the ocean region south of Cape Wave richness resource shows some similarities to that of the wind, mostly due to the waves being generated as a result of the wind forcing on the water surface. WP index (Figure 3a   The indices that account for the wind power resource variability, Cv and Mv, are represented in Figure 4a,b, respectively. In general, the highest Cv values of 0.5 are in the ocean region south of Cape Finisterre, which is not optimal for energy production. However, the Mv index is very good (>0.8) for the same region. These Mv and Cv indices suggest the presence of short events linked to seasonal fronts. Wave richness resource shows some similarities to that of the wind, mostly due to the waves being generated as a result of the wind forcing on the water surface. WP index (Figure 3a) shows its maximum (0.7) on the northwest region, as the RLO index. The lowest WP values are found near the coast, in particular in the Gulf of Cadiz (<0.2). Regarding wave DWNT index (Figure 3b), it should be noted that it reaches the maximum value of 1 in most of the study area, except in the Gulf of Cadiz.
The indices that account for the wind power resource variability, C v and M v , are represented in Figure 4a,b, respectively. In general, the highest C v values of 0.5 are in the ocean region south of Cape Finisterre, which is not optimal for energy production. However, the M v index is very good (>0.8) for the same region. These M v and C v indices suggest the presence of short events linked to seasonal fronts.   Figure 6 represents the risk category that considers the EWS and EWH indices, which affect the  Figure 6 represents the risk category that considers the EWS and EWH indices, which affect the safety related to the maneuverability of energy devices. The wind risk index (Figure 6a The cost category, represented by DC ( Figure 7a) and WD (Figure 7b) is the same for both wind and wave resources, since they only depend on topographic constraints.
The behavior of the wave power resource variability shows a high resemblance with the one shown for the wind. A Cv index (Figure 5a) of around 0.4 is obtained for the northern part of the area, with higher values (~0.5) in the southern part. Mv (Figure 5b) shows values higher than 0.6 in the southwest region of the study area and in a small fringe along the north coast of Spain. For the rest of the area under study, Mv shows values of 0.5. Once again, the lowest values of Cv and Mv indices are located in the Gulf of Cadiz.     The cost category, represented by DC ( Figure 7a) and WD (Figure 7b) is the same for both wind and wave resources, since they only depend on topographic constraints. The classification indices calculated following Equation (5) for wind and Equation (6) for waves are shown in Figure 8. The Wind index (C wind , panel a) shows elevated values (>0.6 in most of the area). Only a very narrow fringe along the cost shows values lower than 0.5. High values can be observed in a wider fringe (~50 km wide) along the coast, the highest ones being (around 0.7) observed from Cape São Vicente to Cape Roca. The wave classification index (C wave , panel b) shows the same pattern as the one observed for the wind index but with lower values (between 0.4 and 0.7) in the whole region. In this case, the coastal fringe with the highest values covers from Cape Roca to approximately Cape Peña. The cost category, represented by DC (Figure 7a) and WD (Figure 7b) is the same for both wind and wave resources, since they only depend on topographic constraints. The patterns observed in Figure 8 are similar to the ones depicted in Figures 9 and 10 obtained through the sensitivity analysis described in Section 2.3.5. In particular, Figures 9a and 10a calculated with Equation (15) are similar to Figure 8a,b, respectively. The RMSE values (Figures 9b and 10b) obtained from Equation (13) are at least one order of magnitude smaller than the average classification indices. This fact is corroborated by the distribution of the RMSEj calculated from Equation (14) (Figures 9c and 10c).  Figure 8a,b, respectively. The RMSE values (Figures 9b and 10b) obtained from Equation (13) are at least one order of magnitude smaller than the average classification indices. This fact is corroborated by the distribution of the RMSE j calculated from Equation (14) (Figures 9c and 10c).
The patterns observed in Figure 8 are similar to the ones depicted in Figures 9 and 10 obtained through the sensitivity analysis described in Section 2.3.5. In particular, Figures 9a and 10a calculated with Equation (15) are similar to Figure 8a,b, respectively. The RMSE values (Figures 9b and 10b) obtained from Equation (13) are at least one order of magnitude smaller than the average classification indices. This fact is corroborated by the distribution of the RMSEj calculated from Equation (14) (Figures 9c and 10c).  The combined index (CI) to classify the hybrid wind-wave energy resource ( Figure 11) shows strong similarities to the patterns observed in Figure 8. In general, the CI is lower than Cwind ( Figure  8a) and higher than Cwave (Figure 8b). The narrow coastal fringe with low values, and the wider coastal fringe with high values (from Cape Peñas to Cape São Vicente) can also be observed in this case. Once again, the lowest values, around 0.5, are observed in the Gulf of Cadiz.  The combined index (CI) to classify the hybrid wind-wave energy resource ( Figure 11) shows strong similarities to the patterns observed in Figure 8. In general, the CI is lower than Cwind ( Figure  8a) and higher than Cwave (Figure 8b). The narrow coastal fringe with low values, and the wider coastal fringe with high values (from Cape Peñas to Cape São Vicente) can also be observed in this case. Once again, the lowest values, around 0.5, are observed in the Gulf of Cadiz. The combined index (CI) to classify the hybrid wind-wave energy resource ( Figure 11) shows strong similarities to the patterns observed in Figure 8. In general, the CI is lower than C wind (Figure 8a) and higher than C wave (Figure 8b). The narrow coastal fringe with low values, and the wider coastal fringe with high values (from Cape Peñas to Cape São Vicente) can also be observed in this case. Once again, the lowest values, around 0.5, are observed in the Gulf of Cadiz.

Conclusions
The viability of exploiting the energy to be produced in hybrid wind-wave farms in the NW Iberian Peninsula for the near future (2026-2045) was investigated using a Delphi method to classify wind and wave energy resources. This classification was based on different indices that represent the resource richness, variability, risk, and cost of the energy resources.
The richness category shows the same patterns for both energy resources, with the highest values in the NW corner of the study area. The variability of the resource also shows a similar pattern for both wind and waves, with the highest values located south of Cape Finisterre. Both resource and variability factors are higher (better) for wind than for wave resource, the downtime factor being better for waves.
The best risk conditions occur along the coastal region, both for wind and waves. In the oceanic part, the conditions improve from north to south. In general, the risk index is worst for wind than for waves.
Finally, a macroscopic evaluation of the hybrid wind-wave resource shows that most of the NW Iberian Peninsula presents good conditions (CI > 0.6) for harvesting energy from wind and wave resources simultaneously. In particular, there are some optimal areas (CI > 0.7), such as the region near Cape Roca and the Galician coast, which had already been highlighted in previous studies on wind [46] and wave [30] resources.
The assessment carried out in this study must be considered as a decision support toolkit in the identification of the best locations for hybrid wind-wave farms, or in the management, adaptation, and resilience of projected and ongoing plans. This assessment only considers the resource from a physical point of view and can be complemented with the technical features of different devices, legal, and environmental constraints, the costs linked to the support of the electric grid connection, and other indices that consider the relation between the maximum production values of both wind and wave resources and the grid energy demand.

Conclusions
The viability of exploiting the energy to be produced in hybrid wind-wave farms in the NW Iberian Peninsula for the near future (2026-2045) was investigated using a Delphi method to classify wind and wave energy resources. This classification was based on different indices that represent the resource richness, variability, risk, and cost of the energy resources.
The richness category shows the same patterns for both energy resources, with the highest values in the NW corner of the study area. The variability of the resource also shows a similar pattern for both wind and waves, with the highest values located south of Cape Finisterre. Both resource and variability factors are higher (better) for wind than for wave resource, the downtime factor being better for waves.
The best risk conditions occur along the coastal region, both for wind and waves. In the oceanic part, the conditions improve from north to south. In general, the risk index is worst for wind than for waves.
Finally, a macroscopic evaluation of the hybrid wind-wave resource shows that most of the NW Iberian Peninsula presents good conditions (CI > 0.6) for harvesting energy from wind and wave resources simultaneously. In particular, there are some optimal areas (CI > 0.7), such as the region near Cape Roca and the Galician coast, which had already been highlighted in previous studies on wind [46] and wave [30] resources.
The assessment carried out in this study must be considered as a decision support toolkit in the identification of the best locations for hybrid wind-wave farms, or in the management, adaptation, and resilience of projected and ongoing plans. This assessment only considers the resource from a physical point of view and can be complemented with the technical features of different devices, legal, and environmental constraints, the costs linked to the support of the electric grid connection, and other indices that consider the relation between the maximum production values of both wind and wave resources and the grid energy demand. 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 also 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 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.