Energetic Potential Assessment of Wind-Driven Waves on the South-Southeastern Brazilian Shelf

: Global electric energy demand is constantly growing, consequently leading towards the usage of renewable energy sources reducing pollution and increasing sustainability. The ocean is a poorly explored renewable energy source; thus, to evaluate the Brazilian wave energy budget, this study investigated the mean behaviour of the wave power rate on the south-southeastern Brazilian Shelf as well as analysed the temporal variability of the wave power rate at the most energetic locations near the coast. Three locations were examined, namely Laguna, Ilhabela and Farol Island, based on the criteria of high means and small standard deviations. The mean wave power rate was approximately 9.08 kW/m on Laguna, 10.01 kW/m on Ilhabela and 15.93 kW/m on Farol Island. The standard deviation identiﬁed in the three locations reached values of 6.47 kW/m on Laguna, 7.59 kW/m on Ilhabela and 13.51 kW/m on Farol Island. Temporal variability analysis was conducted through wavelet analysis. The results show a dominant yearly cycle with a background presence of synoptic cycles, with little deviation between the locations. The El Niño southern oscillation plays a minor role on the energy spectrum of Laguna and does not have a signiﬁcant inﬂuence on Ilhabela and Farol Island.


Introduction
Electrical energy has become essential for humanity and the demand of this resource grows rapidly. This fact has promoted research in this field and concerns about the environment. Several studies have focused on clean and renewable energy sources. Wave energy is one of the most environmentally friendly sources of energy because it is clean, inexhaustible and has no major impact on local hydrodynamics [1].
The first recorded attempt to use wave energy for human needs dates to the thirteenth century, but the efforts to convert it to electrical energy began in 1910, with a device similar to an Oscillating Water Column [2]. Driven by the oil crisis, both Europeans and Americans started studies and developments in wave energy technologies, but the funding of these studies declined after the decrease in oil prices a decade later. However, today, with the menace of oil shortage, wave energy is regaining its value in modern society [2].
Several studies have already been performed in the field of wave energy, from the development and enhancement of different conversion devices [3][4][5][6], to the theoretical assessment of wave energy availability, either by modelling [7][8][9] or by field data analysis [10][11][12].
To obtain an efficient farm to generate electricity from ocean waves, a thorough study has to be performed on different aspects of the implementation. In the global scope, this subject has already been widely studied by several authors. One of the first estimates was made by Isaacs and Seymour [13] who, based on the energy emitted by the Sun, guessed that the portion that corresponds to the waves is around 2.5 TW. More recently, Krogstad and Barstow [14] presented a map of the wave energy estimate around the world based on two years of data from Topex/Poseidon altimetry data. This result provided an early notion of the distribution of wave energy around the globe.
Later, with the increasing development of the wave models, wave energy resource studies became easier. Arinaga and Cheung [15] created an atlas of wave energy with detailed information on its temporal and spatial distribution. Their results, based on a 10-year dataset, present an interleaved winter-summer cycle between the northern and southern hemispheres, with higher energy in the winter. Reguero et al. [16] analysed a 60-year dataset and correlated the temporal variations of wave energy with the main climate anomaly indexes.
Gunn and Stock-Williams [17] used a directional approach for calculating the wave energy availability around the globe. Their study showed an even division of the wave energy budget for both hemispheres, with around 1 TW for each. On the other hand, a great disparity is noticed when comparing the energy availability between the continents, with 20 % of it reaching the North American coastline. Another relevant conclusion by Gunn and Stock-Williams [17] is that, except for North America and Asia, all of the continents have wave energy sufficient to make up for their consumption.
The present study focused on the theoretical wave energy availability along the south-southeastern Brazilian Shelf (SSBS). Brazil is among the countries with the longest coastlines in the world, approximately 7500 km long, in which the southernmost 2777 km of coastline were approached in this work.
Neill and Hashemi [18] made a variability assessment of the waves around the British Isles that complements Arinaga and Cheung's [15] study. These authors showed that, even though this location has an extremely high mean wave energy potential, the variability in this location is also elevated.
Several other studies on wave climate, focusing on their energetic potential, have been conducted over the last years in different parts of the world, e.g., by Stopa et al. [19] in Hawaii; Iglesias and Carballo [20] in El Hierro, part of the Canary Islands; Liang et al. [21] in east China; Hiles et al. [22] and Robertson et al. [23], both on the West Coast of Canada, which is more energetic than the East Coast; Hemer and Griffin [24] in south Australia; and Behrens et al. [25] on the entire Australian coast.
Wave power studies are extremely scarce on the Brazilian coast. Pianca et al. [26] assessed WAVEWATCH III data to evaluate the Brazilian wave climate. These authors classified the climate in terms of height, period and direction. They conducted the analysis for four seasons of the year and took measurements on six representative points of the Brazilian coast.
Parise and Farina [27] conducted a one-year simulation of the waves on the south Atlantic Ocean and collected data on the local wave climate. Losada et al. [28] conducted a 58-year-long analysis of the sea level properties of Latin America and quantified factors such as mean sea level rising trends, tidal variation, storm surges and the contribution of the El Niño southern oscillation.
The wave climate in Brazil was studied by Guimarães et al. [29,30] and Cuchiara et al. [31] with usage of the third generation model SWAN (Simulating Waves Nearshore) to simulate extreme events on the coast of the Rio Grande do Sul state. Parente et al. [32] approached dynamic patterns relating wave conditions to weather seasonality in the study region.
In terms of wave study along the Brazilian shelf, there is a visible dominance of studies on the waves of the Rio de Janeiro State due to its historical and economical importance. On the other hand, some studies on the Rio Grande do Sul State appear due to the importance of Rio Grande Harbour. Considering these factors, this study aimed to evaluate the wave climate on the south-southeastern Brazilian Shelf (Figure 1) to identify the energetic power of wind-driven waves and to provide information on variability and integrated values. Colour bars represent the bathymetry varying from 0 to 5000 m depth. Cube Helix color scheme developed by Green [33]. Satellite images from maps.google.com.
The following sections are organized as follows. Section 2 briefly describes the numerical model used in this work, the boundary conditions and the validation procedure. Section 3 is subdivided into three parts: a general analysis of the mean behaviour of the waves on the SSBS, a local mean analysis of the wave power at the most energetic point of the SSBS and the wavelet analysis to identify the temporal variability of these patterns at selected points. Section 4 contains the conclusions and final considerations of this study.

Materials and Methods
This study was based on the application of the third generation sea state model TOMAWAC to simulate wave conditions over the SSBS over a period of 18 years, between 1997 and 2014, using boundary conditions from WAVEWATCH III and wind data from NOAA's Reanalysis 1. The data used were the integrated parameters Significant Wave Height (H s ), Mean Period (T m ) and Mean Direction (D m ) on the oceanic boundaries and the wind velocity on the free surface.

Numerical Model
The model TOMAWAC (TELEMAC-Based Operational Model Addressing Wave Action Computation) was used to conduct the numerical simulations. This model is a part of the modelling system Open TELEMAC-Mascaret 1 . TOMAWAC is a third generation model based on wave action density conservation (Equation (1)) [34][35][36]: where N represents the directional spectrum of wave action density, x and y are the coordinate system, k x and k y are the components on x and y of the wave number vector, t is the time, and the dot above a symbol denotes its time derivative. Equation (1) represents that, in a general situation of waves propagating in a non-homogeneous and unsteady environment, the wave action density is preserved within the source and sink terms, defined by Q. TOMAWAC calculates wind-driven waves taking into account most of the main physical processes involved, such as shoaling, whitecapping, bottom friction-induced dissipation, nonlinear interactions between waves and depth-induced refraction, and all these effects were taken into account in the simulation. TOMAWAC, however, does not take diffraction and reflection into account [36].
To solve Equation (1), TOMAWAC splits the directional spectrum (N) into a finite number of wave frequencies ( f i ) and directions (θ i ) and solves Equation (1) for each component ( f i , θ i ). The directional spectrum of wave energy, denoted by E( f , θ), can be associated with the directional spectrum of wave action through Equation (2): where ρ is the specific mass of water, g is the gravity acceleration and σ is the angular frequency of the waves given by σ = 2π f . The integration of E( f , θ) along the discretized frequencies and directions yields the energy per unit area of the random multi-directional waves (Equation (3)):

Boundary Conditions and Computational Grid
To perform the numerical simulations, data from NOAA's (National Oceanic and Atmospheric Administration) wave forecasting model WAVEWATCH III 2 [37][38][39] and from the NCEP/NCAR Reanalysis 1 3 [40] were used to represent the liquid and superficial boundaries of the model, respectively.
The model started from the rest and the surface boundary was set to winds from the NCEP/NCAR Reanalysis 1, with a temporal resolution of 6 h and a spatial resolution of 1.875°(T62 Gaussian Grid). These data were bilinearly interpolated to the numerical domain.
The oceanic boundaries were set by the imposition of the integrated parameters H s , T p (peak period) and D m . The data were generated by WAVEWATCH III and have a spatial resolution of 30 arc minutes and a temporal resolution of 3 h. Bathymetric data were acquired from the Brazilian Navy Admiral Charts and from the GEBCO (General Bathymetric Chart of the Oceans) database.
The validation and the case study simulations used the same data sources. The validation simulation lasted four years, from 2011 to 2014, the period when data from PNBOIA (Programa Nacional de Boias) were available. The simulation for the case study lasted 18 years from 1997 to 2014. The computational grid used in the simulations extents from the south of Rio Grande do Sul, in the border between Brazil and Uruguay, to the north of the state of Espírito Santo. The mesh also extents 650 km towards the ocean (Figure 1). Figure 2 shows mesh used in the simulations. The mesh contains 405,833 triangular elements, totalling 205,617 nodes. The mesh has varying levels of refinement, from 8000 m between nodes in the oceanic region, 2500 m on the deep continental shelf, to 1000 m in regions shallower than 100 m. In regions of interest, the nodes are 300 m apart, approximately. With this configuration, using TOMAWAC version V7P2R2 in scalar mode, the case study simulation took approximately four weeks to run on a Intel ® Core ™ i7-4960X CPU @ 3.60 GHz × 12, using around 8 GiB of RAM. The results of the case study simulation were saved in a 75.7 GB file.

Validation
The validation of The PNBOIA is an ongoing project for gathering and providing oceanographic and meteorological data along the Brazilian coast and the adjacent Atlantic Ocean. This project aims to contribute to the characterization of the Brazilian marine environment as well as to attend to the necessities of navigation security, oceanic meteorological forecasting and human life protection. The project, which counts with the assistance of several governmental institutions, has more than 300 buoys in its database providing information about oceanic circulation, water properties and waves and the associated atmospheric conditions. The time series provide a visual representation of the comparison between TOMAWAC and the buoys. This comparison can be quantified for better understanding, using well known [41][42][43][44][45] error metrics in this field of study. The equations to compute these metrics are presented in Table 1. In the equations, T i are the data modeled by TOMAWAC, B i the data measured by the buoys and n the number of data entries in a time series. To complement the results, some statistics about the data from both TOMAWAC and the buoys are presented in Table 2 with the error metrics. Table 1. Mathematical equations to compute the error metrics. The comparison of the time series is shown in Figure 3, where their period range from 1 July 2011 to 31 April 2014. The time series shows that there is a close agreement between the model and the measured H s . To quantify the quality of the numerical results, some statistical parameters were calculated during periods in which there were data collected from the model and measurements (Table 2). Average values show that the model results approach the measured data. Nevertheless, on average, numerical results tend to underestimate empirical data.
These mean values are considerably robust for a wave prediction model. In their wave resource study, Edwards et al. [46] matched ECMWF's (European Centre for Medium-Range Weather Forecasts) Cy331 WAve Model (WAM) with a Datawell Waverider Buoy and obtained a mean difference of 0.17 m. Studying wave energy power, Santos [47] found mean differences of approximately 0.30 m both in the calibration and in the validation of the numerical model. Melo et al. [43] found differences of approximately 0.22 m while assessing the results of the WAVEWATCH III model. Similar pattern of comparison was found by Iuppa et al. [48] validating the performance of SWAM model against satellite data.
Another statistical parameter calculated is the Root Mean Square Error (RMSE). This parameter is used because it makes larger errors more important than minor errors in calculating final results. RMSE is a widely used statistical parameter [43,46,47] for this type of analysis because, as the models are never perfectly adequate, RMSE neglects small differences in a time series that are not considered significant errors. The RMSE of the H s is an average 0.5 m, a slightly lower value than values calculated in other calibration and validation studies [42,43,46,47].
The mean T p is underestimated by TOMAWAC with a difference of approximately 0.9 s. The results obtained by Lalbeharry [42], Melo [43], and Santos [47] show approximately the same mean difference. Although the calibration and validation of wave models is primarily based on the H s parameter, due to the amount of errors in the T p data measured by the buoy, some filters can be applied to clean the data, and this often results in a robust correlation between the model and measured T p . The RMSE is approximately 2.0 s, which is comparable with aforementioned studies.
The correlation coefficient for H s is between 0.88 and 0.91, the lowest value being found near the coast, indicating a less strong correlation at this region. For T p , the correlation coefficient is slightly higher than for H s (between 0.89 and 0.93) indicating the same level of correlation for both parameters.
To conduct the last verification of the model validation, the Scatter Index (SI) was calculated for the two parameters. The advantage of SI is that being dimensionless, contrary to the RMSE, it can be compared for the same or different variable and, among different magnitudes of data, allowing the comparison between ocean and coastal zones. Based on the SI values, H s calculated in deep regions is slightly more precise than H s calculated near the near shore. Additionally, the model output for T p approximates the measured data better than values of H s .
Finally, to visualise the correlation between the measured data and TOMAWAC's results, scatter plots were created correlating both datasets. Figure 4 shows the scatter plots of the three validation datasets used. The scatter plots were made using discrete bins to avoid the cluttering caused by the large amount of data points.
The plots in Figure 4 emphasize the conclusions drawn above. The data are distributed along the 45°line with some scattering, which lowers the correlation coefficient. The solid diagonal lines in each scatter plot shows the linear regression of the dataset comparison. The linear regression enhances the slight underestimation of the wave heights by TOMAWAC; Cavaleri [49] pointed out a few sources of underestimation from wave models as being related to whitecapping processes and non-linear interactions. Nevertheless, the correlation between measured and modelled wave heights is satisfactory, although the issue of the underestimation could, in a future study, be addressed with a different calibration.

Results and Discussion
To investigate the energetic power of wind-driven waves along the SSBS, a long-term simulation was conducted over 18 years from January 1997 to December 2014. The results were investigated using mean values to identify the most energetic points near the coastal region. Spatial and temporal variability of wave power and wind intensity were investigated using wavelet analysis that was conducted following the methods of Torrence and Compo [50] and the bias correction proposed by Liu [51].

Temporal Mean Analysis
To obtain initial insight into the available resource of wave energy on the SSBS, the temporal mean of the modelled result was extracted over the entire studied region to calculate the mean surfaces of the calculated parameters. The mean H s ( Figure 5) shows similar behaviour along all of the SSBS, with values approximately 1.6 m over the continental shelf and reaching up 2.5 m in the oceanic region.
The mean H s surface shows a discontinuity on the Santos Basin, with a gap within 1.5 m to 1.7 m, which begins at Farol Island (42°W, 23°S). This occurs because the island retains a part of the waves that propagate to the southwest, contributing to a decrease of H s in the Santos Basin.
Southwest of the domain, the region of the Rio Grande do Sul coast (30°S to 34°S) has small wave amplitudes that are approximately 0.6 m. These smaller wave amplitudes appear to be the result of the mildly sloping beaches in this area [52]. Mildly sloped beaches dissipate much of the wave energy and do not allow greater wave heights to occur in this region.
The first exception to this pattern occurs on the coast of Laguna (48.74°W, 28.58°S), located in the Santa Catarina State, which has a mean H s of 1. The next remarkable location is Farol Island, near Arraial do Cabo city, which has the highest mean H s of all of the studied coastlines (up to 1.6 m), which is also due to its projection towards the ocean. TOMAWAC determines the wave propagation direction based on the direction of the waves relative to true north in a clockwise direction [36]. The mean D m surface ( Figure 6) shows that the main direction of propagation of the waves is 235°(waves going to southwest), with a deflection eastward in the direction of the ocean (180°) and another deflection to the northwest (315°), perpendicular to the shoreline. This latter deflection appears to be due to the refraction of the waves. This behaviour of waves, going to the southwest, complements the discontinuity found in the mean H s surface. The waves arrive at Farol Island with a mean D m of 230°and are blocked by the island, thereby reducing the H s of the waves that arrive at the Santos Basin.
The mean surface of T m (Figure 7) shows waves with periods between 8 s and 9 s on the south and southeastern regions, respectively. These values are practically constant for the entire SSBS, except near the coastline, where it is possible to observe regions where the waves begin to be affected by the bottom. In these regions, the T m lowers to 5 s, and wave height and length decrease. The location of this variation is where dissipative processes are the least pronounced. These locations are Laguna, Ilhabela, and Farol Island. This finding leads to the prediction that these locations experience the most energetic waves on the SSBS.
Finally, the mean P w field (Figure 8) exhibits the same behaviour as H s because the wave energy is proportional to the square of the amplitude and the period [1]. The mean P w in this oceanic region is approximately 22 kW/m and reaches up to 30 kW/m in the southern part of the domain. These elevated mean values in the southern region, for both P w and H s , are due to the wind fetch needed for the waves to develop in the Antarctic Ocean [1,17].
The wave energy is close to zero at locations very near the coastline on the southernmost part of the SSBS and in the Santos Basin. On the other hand, the most elevated values were observed in Laguna and Ilhabela, with 10 kW/m and Farol Island, with a mean value that reach up to 15 kW/m. Generally, the mean values simulated by the TOMAWAC model were quite similar to the Brazilian wave energy described by Pianca et al. [26].

Local Wave Power Analysis
The installation of wave energy converters should be made as near the coast as possible to reduce costs during construction and operational periods. Therefore, installation could be conducted at the most energetic locations: Laguna, Ilhabela and Farol Island.
An analysis of the P w available for conversion was conducted in the area of study. The selected locations are highlighted in Figure 8. A colormap of each location is presented in Figure 9. The southernmost point in the coastal region is located near the city of Laguna. This point, located approximately 7 km from the coast has a mean H s of 1.2 m ( Figure 5), a relatively high value compared to other locations along the southern Brazilian coast.
The P w spatial mean around Laguna (Figure 9a) exhibits an approximately linear behaviour parallel to the shoreline. This characteristic, combined with the bent geometry of the beach, contributes to the elevated mean P w of this location.
Its mean surface shows that, of the entire coastline, the potential is higher perpendicular to the coast of the Santa Marta Cape. Examining the value of P w from the coastline towards the ocean, a substantial increase from 2.0 kW/m to 7.0 kW/m within the first 4 km, and a subtle rise up to 10.0 kW/m within 20 km.
An important aspect of P w is its temporal variability. Locations with smaller means and low variability are more efficient for energy conversion than locations with higher means and high variability. Neill and Hashemi [18], for instance, studied Great Britain, which has one of the highest P w in the world. P w in Great Britain is more than 60 kW/m, which is unevenly distributed between the winter (up to 110 kW/m in average) and summer (close to zero).
Laguna, on the other hand, does not provide P w as high as Great Britain but, in terms of variability, is far more stable. Temporal variability in Laguna is expressed in terms of standard deviation in Figure 10a. The general characteristics of the standard deviation are rather similar to those of the mean surface.
In the first few kilometres, there is a smaller increase in the standard deviation, from 1.5 kW/m to 4.5 kW/m, than in the mean power. Further away, there is a standard deviation of 8 kW/m, associated with a higher average of 10 kW/m. Thus, the suitability, based on mean values, is better in the first 5 km off the coastline. In contrast with the linear behaviour of the waves mentioned above, Ilhabela has completely different wave characteristics on each side (Figure 9b). The southwestern sector offers a starting mean P w of 4.5 kW/m on the first 2 km off the coastline and does not show considerable growth further away. The standard deviation (Figure 10b) of this sector exhibits the same behaviour as the mean power, starting at 3 kW/m but with slightly less growth advancing towards Santos Basin.
The northeastern side of Ilhabela has averages (Figure 9b) close to 1.75 kW/m near shore, which grow linearly to approximately 4 kW/m within 6 km from the coast. Despite being lower than the values found on the opposite side of the island, these are compensated by the lower standard deviation (Figure 10b) of 1.4 kW/m near the shore, which increase up to 2 kW/m further into the ocean.
The northeastern sector has steady mean values in the area that covers the inner continental shelf from 5 km to 10 km off the coast. Again, the standard deviation behaves similarly to the mean power but with a magnitude 1.5 to 1.8 times smaller.
The northernmost site selected was Farol Island, near Arraial do Cabo, in the Rio de Janeiro State. The entire region is located in a peninsula. Furthermore, Farol Island is projected towards the ocean, which makes dissipative processes much less effective on the waves that arrive to the island.
The mean P w (Figure 9c) immediately south of the island is 15 kW/m, with an almost insignificant growth up to 17 kW/m, away from the island going south. At the eastern part of the island, the mean P w is reduced to 12 kW/m at approximately 1 km off the coast. This value is constant further away into the sea. The behaviour is similar on the western side.
The standard deviation at these locations is fairly high (Figure 10c), but the high values are a consequence of high mean values. At the southern end of Farol Island, the standard deviation is 11 kW/m immediately beside the coastline, increasing up to 15 kW/m further towards the sea. The east side has a standard deviation of 7 kW/m and this value remains constant and similar to the mean values. The West sector, on the other hand, has a standard deviation of 11 kW/m.
Another notable location is the bay in between Farol Island and Armação dos Búzios (41.9 W, 22.72 S), located at approximately 10 km off the shoreline. This spot has the largest mean values compared to the standard deviation, twice as large in some locations. However, the mean power varies by approximately 7 kW/m and the standard deviation decreases to 3 kW/m.

Temporal Variability
To investigate patterns in temporal variability of the three selected sites, Cross-Wavelet analysis was conducted according to the method of Torrence and Compo [50] and Liu [51]. This analysis consists of the deconstruction of two time series in frequency ranges and the correlation of both decomposed series with each other to detect similarities in variability patterns of their respective periods.
The time series used for comparisons were extracted from points approximately 15 km away from the coastline of all the three studied sites ("X" in Figures 9 and 10). The wavelets were used to correlate the waves P w with the wind intensity. Figure 11a shows both extracted times series on Laguna. P w ranged from zero to 21.5 kW/m and the wind intensity ranged from 1.1 m/s to 8.1 m/s. Figure 11b illustrates the Morlet local wavelet spectrum for Laguna: the most energetic frequencies are shown in dark red. The black outline encloses the parts of the spectrum that have 95 % statistical confidence, while the regions below the dashed line are the regions where the border effects are present in the time series.
The Morlet local wavelet on Laguna (Figure 11b) shows a dominant annual cycle, which is existent for nearly the entire time series. Minor influences occur in 2004 and 2009, which indicates that there is variability in periods of 12-16 months. The local spectrum also shows synoptic variability with stronger effects during 2-7 days associated with the passage of meteorological systems. Figure 11c presents the time-averaged Morlet local spectrum (global spectrum) on the frequency domain (black continuous line) and the 95 % significance level (gray dashed line). The global spectrum has two dominant period ranges, indicating that synoptic variability is embedded in a more energetic annual cycle on the P w available on Laguna. The Mexican Hat local wavelet spectrum on Laguna (Figure 11d) shows predominance of the annual cycle and details the oscillatory pattern over time. Oscillations happen twice a year and indicate that a well-defined cycle lasts six months, thereby implying that winter-summer variation exists.
The mid-section of the local spectrum covers the annual cycle and extends itself to approximately 600 days (two years), which corresponds to variability in the El Niño Southern Oscillation (ENSO). The frequency of these oscillations ranges from occurrences every 1.5 years to less frequent oscillations. The first ENSO in the time series occurred in 1998, an extremely strong among the ENSO [53][54][55][56][57][58], followed by ENSO of 2002 [59,60], 2007 [61] and two events of smaller intensity in 2010 [62] and 2012, respectively. A time series with the Multivariate ENSO Index can be found at NOAA 4 . 4 esrl.noaa.gov/psd/enso/mei Figure 11e shows the global spectrum for the Mexican Hat wavelet and three distinct time scales are dominant. The shortest and less energetic is the 15 to 30 days oscillation, followed by the intra-seasonal time scales and the most energetic annual cycle. Again, this is the most important fraction of the spectrum.
The Mexican Hat wavelet also shows that the ENSO plays a significant role in the energy spectrum of this location.
Moving on to Ilhabela, Figure 12a shows P w in the interval between zero and 24.6 kW/m and wind intensity from 1.0 m/s to 7.7 m/s. The Morlet local wavelet spectrum (Figure 12b) shows a well-defined presence of an annual cycle along the entire time series associated with the largest part of the global spectrum (Figure 12c). The dominant cycle observed for Ilhabela on the Mexican Hat Wavelet (Figure 12d) is also the yearly cycle [63], which holds most of the energy spectrum. This is primarily due to the well-defined seasonality of this region [63,64]. Associated with the annual pattern, there are cycles with considerable influence from 6 to 18 months, the former being associated with the seasonal pattern and the latter with short-term ENSO influence.
Gan and Rao [64] analysed cyclogenesis frequency for several years on the South American coast and showed a strong tendency for cyclone generation during the winter and a weaker tendency during the summer, with the autumn and spring being transitional seasons. Because cyclones have a direct effect on wind-driven waves, this difference between seasons is certainly one source of the yearly variability patterns in waves.
Different from Laguna, the annual cycle in Ilhabela has four oscillations within one year, which implies that there is a distinct behaviour in each season. This seasonal behaviour can also be observed on the global spectrum (Figure 12e), which shows much stronger seasonal variability in Ilhabela than in Laguna. This pattern accounts for a significant part of the Ilhabela energy spectrum.
Finally, at Farol Island, there is a higher P w (Figure 13a), with common values up to 42.0 kW/m, caused by the intense average winds in this region, between 1.7 m/s and 9.0 m/s. The high intensity of winds appear to be due to the greater influence of the of south Atlantic Anticyclone at middle latitudes over this island [65,66], thereby making Farol Island have the highest values of P w in the study region.
The Morlet local wavelet spectrum on Farol Island (Figure 13b) shows a highly concentrated energy on the yearly cycles, with little spread to higher or lower frequency oscillations. The local spectrum also shows that, for this location, there are no cycles with periods between one and two weeks that limit the synoptic time scale oscillations to one week at most.
The global Morlet spectrum ( Figure 13c) shows a peak of energy on the annual cycle and a decrease in synoptic variability. In addition, the Mexican Hat local wavelet spectrum on Farol Island (Figure 13d) displays the same narrow interval for the yearly cycle found on the Morlet spectrum. A pattern similar to that found on Ilhabela can be observed here, with four high energy cells within a year and less energetic patterns from 6 to 18 months.
The Morlet spectrum (Figure 13e) also shows that the two-week pattern is not observable on Farol Island. The shortest time scale captured by the Mexican Hat spectrum is a monthly pattern, which still has low energy.
In addition, the seasonal cycle also has little influence on the spectrum. On Ilhabela, the long period ENSO has a weak presence; on Farol Island, ENSO has a single occurrence in the four-year time scale, in between 2006 and 2007.
The results show that, over 18 years, Farol Island is more stable during the annual cycle. On the other hand, the less stable location is Laguna, which exhibits no pattern at all over time. A comparison can be performed between global energy spectra.
The global energy spectra (Figures 11c,e-13c,e) show most of the energy located on the annual cycle for the three locations but, among the locations studied, the most energetic one is Farol Island, with 1.2 times more energy than Ilhabela and 1.9 times more than Laguna. Additionally, the Morlet global spectra (Figures 11c-13c) show that the yearly cycle has 23, 38 and 47 times more energy than the synoptic cycle on Laguna, Ilhabela and Farol Island, respectively.
Similarly, the Mexican Hat global spectra (Figures 11e-13e) determined that there is 2.7, 3.6 and 4.3 times more energy for the annual cycle than there is for the seasonal and monthly cycles on Laguna, Ilhabela and Farol Island, respectively. The mean Mexican Hat global energy spectra (Figures 11e-13e) show a similar behaviour for the three sites for cycles from 8 to 650 days. For longer cycles, the behaviour becomes quite different. Laguna has the most energetic effect due to the long period (more than two years) ENSO, with three significant appearances within 18 years.
The annual cycle is dominant on the three spectra and the strongest was verified in Farol Island. The synoptic and seasonal patterns are stronger on Ilhabela and have about the same energy contribution on the Farol Island and Laguna.
Finally, Table 3 shows basic statistical parameters calculated for the P w time series used on wavelet analysis. Corroborating previous results, all of the values calculated for the time series on Farol Island are substantially higher than those calculated for the other localities.
The mean P w on Laguna and Ilhabela are similar, but P w on Farol Island is approximately 65 % higher. The same behaviour was observed in the standard deviation, but with a difference of approximately 90 %.
The maximum values, on the other hand, showed the greatest difference between Laguna and Ilhabela, and less from Ilhabela to Farol Island. However, these peak values are single values in the time series and do not occur frequently enough to have a great impact on the final power output.
The most remarkable difference lies on the integrated values that show the similarities between Laguna and Ilhabela. The integrated value exhibits a 60 % increase on Farol Island, with 209.5 MW/m.

Conclusions
The third generation model TOMAWAC was validated using a four-year simulation to analyse wave energetic potential. Afterwards, a case study simulation of 18 years, from 1997 to 2014 was conducted on the SSBS.
The temporal mean results of the simulation show a generally constant behaviour of wave parameters H s , T m , D m and P w on the oceanic region. As the waves approach the coastline, the H s , T m and P w values decrease due to bottom friction, while the D m undergoes the refraction effect and deflects towards the coast.
Three study sites were chosen based on the H s and P w parameters. The locations are Laguna, with a mean H s of 1.3 m and P w , with a mean of 10 kW/m; Ilhabela, with mean H s and P w of 1.2 m and 10 kW/m, respectively; and Farol Island, with means of 1.6 m and 15 kW/m for H s and P w , respectively.
After selecting the study sites, Cross-Wavelet analysis was used on a time series of P w and wind intensity to define the variability cycles for each location. The main variability cycle observed at the three locations is the annual cycle, which was present over the 18 years simulated, with a strong seasonal presence on Ilhabela. The synoptic cycles are, as expected, present at all localities but with little influence on Farol Island.
ENSO has an important role on wavelet spectra, appearing as relatively short-term oscillations ranging from 1.5 to 2 years. Long-term ENSO have a smaller influence on Ilhabela and a major influence on Farol Island. These appear as 5 to 10-year oscillations on Laguna and as high energy events on Farol Island, with periods lasting approximately 10 years.
Thus, the most feasible location found in this study for a hypothetical installation of wave energy converters is the region adjacent to Farol Island. Higher mean values were obtained at this location, in addition to the largest integrated values. In terms of variability, Farol Island has the least energetic synoptic cycles, meaning that there is greater stability over time. Most of the energy oscillation is due to the annual cycle and the long-term variability caused by ENSO.

Acknowledgments:
The authors are grateful to the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under contract 304227/2016-1 and to the Agência Nacional do Petróleo (ANP) regarding the Programa de Recursos Humanos (PRH-27) for the fellowships. The authors are also grateful to the Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS) for sponsoring this research under the contracts 17/2551-001159-7, to the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for bursaries related to the Programa Nacional de Pós Doutorado. Additional acknowledgement goes to the Brazilian Navy for providing detailed bathymetric data for the coastal area, to NOAA for supplying the wind and wave datasets, and to Électricité de France (EDF) for the public licence for the Open TELEMAC-Mascaret system. A special thanks to the Supercomputing Center of the Federal University of Rio Grande do Sul (CESUP-UFRGS) and to the Sdumont Supercomputer from the Laboratório Nacional de Computação Científica (LNCC) (SDUMONT-2017-C01#166515) where most of the computational work was carried out.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: