Assessment of Offshore Wave Energy Resources in Taiwan Using Long-Term Dynamically Downscaled Winds from a Third-Generation Reanalysis Product

: In this study, long-term wind ﬁelds during 1991–2010 from the Climate Forecast System Reanalysis (CFSR) were dynamically downscaled over Taiwan and its offshore islands at a 5 km horizontal resolution using the Weather Research and Forecasting (WRF) model. Simulations of the 10 m (above sea level) dynamically downscaled winds served as the atmospheric forcing for driving a fully coupled wave-circulation model. The sea states of the waters surrounding Taiwan during 1991–2010 were hindcasted to evaluate the offshore wave energy resources and optimal wave energy hotspots. This study reveals that the southeastern offshore waters of Taiwan and the Central Taiwan Strait exhibited the highest mean wave power density (WPD), exceeding 20 kW/m. The annual mean WPD, incidence of the hourly WPD greater than or equal to 4 kW/m, monthly variability index and coefﬁcient of variation of the WPD indicated that the sea areas located between Green Island and Orchid Island (OH_1), southeast of Orchid Island (OH_2), south of the Hengchun Peninsula (OH_3), and north of the Penghu Islands (OH_4) were the optimal hotspots for deploying wave energy converters. The most energetic months were October for OH_1 and OH_2 and November for OH_3 and OH_4, while the wave power was weak from March to June for OH_1, OH_2 and OH_3 and in May for OH_4. The wave direction is prevailingly east-northeast for OH_1, OH_2 and OH_3 and nearly northeast for OH_4. These phenomena reveal that wave power in the waters offshore Taiwan is induced primarily by the northeast (winter) monsoon. The exploitable annual WPD was estimated to be 158.06, 182.89, 196.39 and 101.33 MWh/m for OH_1, OH_2, OH_3 and OH_4, respectively.


Introduction
Offshore winds and waves have the greatest energy potential among all types of renewable energy, as the oceans cover a majority of the Earth's surface [1]. Wave energy is unique because it is the most concentrated form of renewable energy on Earth, with its power density being much higher than that of either wind or solar energy [2]; hence, wave energy converters have been installed and operated in many countries for decades [3]. Waves are created when wind moves over the ocean surface. Even small ripples on the ocean surface allow the wind to transfer some of its energy to the waves, causing the waves to grow from wind interactions over long distances. Globally, wave power resources have been estimated to reach 2.11 ± 0.05 TW (TeraWatt) [4]. The most appropriate sea areas for exploitation, featuring power densities of 60-70 kW/m, are found at medium-high

Description and Configuration of SCHISM-WWM-III
The Semi-Implicit Cross-Scale Hydroscience Integrated System Model (SCHISM) was utilized to hindcast the long-term tides and currents in the waters surrounding Taiwan. SCHISM, which was originally implemented by [27], is derived from the Semi-Implicit Eulerian-Lagrangian Finite-Element (SELFE) model [28] with many improvements. The shallow water equations are solved by employing a highly efficient and accurate fi-

Description and Configuration of SCHISM-WWM-III
The Semi-Implicit Cross-Scale Hydroscience Integrated System Model (SCHISM) was utilized to hindcast the long-term tides and currents in the waters surrounding Taiwan. SCHISM, which was originally implemented by [27], is derived from the Semi-Implicit Eulerian-Lagrangian Finite-Element (SELFE) model [28] with many improvements. The shallow water equations are solved by employing a highly efficient and accurate finite-element/finite-volume method with the Eulerian-Lagrangian algorithm and the application of the hydrostatic and Boussinesq approximations. The stability restrictions avoid constraints associated with the Courant-Friedrichs-Lewy (CFL) condition by applying semi-implicit numerical schemes to the governing equations. Additionally, using the Eulerian-Lagrangian method for the nonlinear advection terms further relaxes the stability constraints of the numerical model. SCHISM introduces the no-mode-splitting technique, which eliminates the numerical errors derived from an internal-external mode splitting procedure [29]. Consequently, owing to their high efficiency and accuracy, the SELFE model and SCHISM have been widely used for investigating the transport and fate of fecal coliform contamination in a tidal estuarine system [30], predicting coastal inundation [31,32] and riverine flash floods [33], and assessing the tidal current energy in nearshore waters [34]. Compared to a three-dimensional model, a two-dimensional model requires fewer computing resources and shorter model execution times; in addition, the hydrodynamics in a large-scale ocean can be well simulated by a two-dimensional model. Thus, because a two-dimensional model is preferred over a three-dimensional model for long-term wind wave modeling, SCHISM-2D (a depth-integrated version of SCHISM) was employed in the present study. A Manning coefficient of 0.025 and a time step of 120 s were prescribed in SCHISM-2D based on the seafloor material in Taiwanese waters reported by previous studies [35][36][37][38].
A third-generation, state-of-the-art spectral surface wave model, namely, the Wind Wave Model-III (WWM-III), was adopted to hindcast the sea state in the nearshore and offshore waters of Taiwan. Phase-averaged spectral wave models have the ability to simulate wave parameters on the free surface of a large water body [39]. WWM-III is a derivative work of the original WWM-I [40] that solves the wave action equation (WAE) by using new numerical schemes and revised physics on an unstructured grid [41]. A fractional time step method [42] was applied in WWM-III to integrate the various terms of the WAE, while the integration in spectral space was accomplished through the ultimate quickest approach [43]. The number of frequency bins was 36, and the frequency limits for the discrete wave period were in the range of 0.03 Hz to 1.0 Hz. The directions for the simulation were divided into 36 directional bins with a minimum and maximum of 0 • and 360 • , respectively. A peak enhancement factor of 3.3 based on the report from the Joint North Sea Wave Project (JONSWAP, [44]), a wave breaking coefficient of 0.78 and a bottom friction coefficient of 0.067 are preset in WWM-III.
Utilizing the different time steps for ocean and wind wave models and sharing the same numerical mesh are better methods to enhance the computing efficiency of a wavecirculation coupled modeling system; for this reason, a time step of 600 s was chosen for WWM-III. In other words, SCHISM-2D and WWM-III exchanged the wind speed, depth-integrated tidal current, water surface elevation and wave radiation stress on each grid every 5 hydrodynamic time intervals. Details on the coupling procedure for the SCHISM-WWM-III coupled modeling system are available in [45]. Due to its high numerical efficiency and robust matrix solver, SCHISM-WWM-III has been successfully implemented to investigate wave power resources and produce potential storm wave hazard maps for Taiwan [35,36], evaluate the wind forcing effect on storm wave hindcasting [37], quantify storm tide-induced nonlinear interactions [38], simulate extreme waves caused by a super typhoon [46,47], and evaluate wave-current interactions [48].
A large computational domain is essential for long-term wave climate hindcasting, especially when hindcasting typhoon waves and surges. Therefore, the domain from 111 • E to 135 • E and 18 • N to 30 • N was created for SCHISM-WWM-III. This wide area has been shown to be sufficient for simulating typhoon-driven sea states in the offshore waters of Taiwan. The entire computational domain comprised 168,829 unstructured grids and 327,955 triangular elements (as shown in Figure 2a). A fine mesh with a resolution of approximately 200 m was distributed along the coast of Taiwan (as shown in Figure  2b) and its offshore islands (as shown in Figure 2c), while a coarse mesh with a resolution of approximately 40 km was configured in the open ocean far from the coastal zone. The fine mesh developed in the present study is able to characterize the surf zones and shallow waters [46,[48][49][50][51]. Notably, the performance of the JONSWAP spectrum in the simulation of long-term or extreme sea states in the offshore waters of Taiwan was impressive, and this spectrum has been widely used to simulate wave conditions worldwide [35][36][37][38][46][47][48]52], although the JONSWAP spectrum was originally applied in the North Sea.
where pe(x,y) is the value to be estimated (the bathymetric grid in SCHISM-WWM-III), pm(x,y) is the known value (the bathymetric grids in the gridded datasets), and di is the distance from the data point to the estimated point.
All the initial conditions, such as sea surface elevation and horizontal velocity, on each grid of SCHISM-WWM-III were set equal to zero, and the first 15 days were taken as the spin-up time for SCHISM-WWM-III. The tidal boundary conditions of SCHISM-WWM-III are M2, S2, N2, K2, K1, O1, P1 and Q1 tidal constituents, which were derived from a regional inverse tidal model (China Seas and Indonesia, [53]). The 3-h wave spectral from operational global WW3 of NOAA was applied to the wave boundary conditions for SCHISM-WWM-III.

Wave Buoy Data
Hourly measurements of the SWH, peak wave period and wave direction recorded at four wave buoys by the Central Weather Bureau (CWB) of Taiwan in 2010 were secured to evaluate the hindcasting performance of SCHISM-WWM-III using the wind fields from CFSR and CFSR_5km. The four wave buoys shown in Figure 3a, namely, the Longdong, Hualien, Hsinchi and Xiaoliuqiu buoys, are located in the northern, eastern, western and southern waters of Taiwan, respectively.  was interpolated by the inverse distance weighting (IDW) method from a combination of these two gridded datasets (as shown in Figure 3a). The IDW formula is: where p e (x,y) is the value to be estimated (the bathymetric grid in SCHISM-WWM-III), p m (x,y) is the known value (the bathymetric grids in the gridded datasets), and d i is the distance from the data point to the estimated point.

Computation of Offshore Wave Energy
The wave power per unit crest length in the direction of wave propagation, i.e., the wave power density (WDP), can be obtained through the spectral output of the wind wave  All the initial conditions, such as sea surface elevation and horizontal velocity, on each grid of SCHISM-WWM-III were set equal to zero, and the first 15 days were taken as the spin-up time for SCHISM-WWM-III. The tidal boundary conditions of SCHISM-WWM-III are M2, S2, N2, K2, K1, O1, P1 and Q1 tidal constituents, which were derived from a regional inverse tidal model (China Seas and Indonesia, [53]). The 3-h wave spectral from operational global WW3 of NOAA was applied to the wave boundary conditions for SCHISM-WWM-III.

Wave Buoy Data
Hourly measurements of the SWH, peak wave period and wave direction recorded at four wave buoys by the Central Weather Bureau (CWB) of Taiwan in 2010 were secured to evaluate the hindcasting performance of SCHISM-WWM-III using the wind fields from CFSR and CFSR_5km. The four wave buoys shown in Figure 3a, namely, the Longdong, Hualien, Hsinchi and Xiaoliuqiu buoys, are located in the northern, eastern, western and southern waters of Taiwan, respectively.

Computation of Offshore Wave Energy
The wave power per unit crest length in the direction of wave propagation, i.e., the wave power density (WDP), can be obtained through the spectral output of the wind wave model. The WPD (P) formula is given as: where ρ is the density of seawater, g is the gravitational constant of acceleration, σ is the wave relative angular frequency, θ is the wave propagation direction, S(σ,θ) is the directional wave variance density spectrum, and d is the depth of water. In addition, c g (σ,θ) is the wave group velocity, which can be expressed as: where κ = 2π/λ is the wavenumber and λ is the wavelength. In deep water, i.e., d > λ/2, Equation (3) can be simplified to c g (σ,d) = g/4πσ. Hence, Equation (2) is rewritten as: The nth order spectral moment (m) is defined as: The SWH H s and the energy period T e in terms of the spectral moment are expressed as: Energies 2021, 14, 653 8 of 25 Substituting Equations (5)- (7) into Equation (4), a simple formula for calculating the WPD (P) can be obtained: ρ is assumed to be 1025 kg/m 3 in the present study, and thus, Equation (8) becomes: T e is best evaluated through other variables because it is rarely specified by the measurements of wave parameters. Thus, the peak period T p was adopted to calculate the WPD through multiplication by the coefficient α: in which α = 0.9 was adopted from the standard JONSWAP spectrum with a peak enhancement factor of γ = 3.3 [35,[54][55][56][57][58]. Therefore, the theoretical WPD at an arbitrary point on the computational mesh was estimated by means of the following formula: The effects of wave refraction and shoaling on wave power can be ignored since the continental shelf of Taiwan is small. Accordingly, Equation (11) is sufficient to represent the magnitude of the offshore WPD in Taiwan [35,57]. Finally, the total annual wave energy per unit length can be produced for a given time interval: More details regarding Equations (2)-(8) can be found in [59].

Metrics for Determining the Optimal Wave Energy Hotspots
Four criteria, i.e., the long-term annual mean WPD, incidence of exploitable wave power, monthly variability index (MVI) and coefficient of variation (COV), were used to identify the optimal wave energy hotspots in the offshore waters of Taiwan. Offshore sea areas are considered to be optimal hotspots for the deployment of wave energy converters as long as there is an abundant WPD and a higher incidence of exploitable wave power but a lower MVI and COV. Exploitable wave power is defined as a WPD equal to or greater than 4 kW/m in the present study. The spatiotemporal variabilities were characterized by adopting the MVI and COV: where P max_mavg_lt and P mix_mavg_lt are the maximum and minimum long-term (20-year) monthly average WPDs, respectively, and P yavg_lt is the long-term (20-year) yearly mean WPD.
where σ wp is the standard deviation of the long-term (20-year) average WPD.

Hindcast of Wave Parameters Using the CFSR and CFSR_5km Wind Fields
To confirm that the regional dynamically downscaled wind fields are superior to the global reanalysis fields in hindcasting wave parameters, hourly winds from the original CFSR and from CFSR_5km were employed as atmospheric boundary conditions for  (Figure 4e-h). The largest differences for the SWH (hindcast minus measurement) range between 6 and 8 m are detected at the Longdong buoy when SCHISM-WWM-III was given the wind fields from the CFSR (Figure 4a). The SWH differences were always within ±2 m at the Hualien buoy due to the utilization of the CFSR_5km winds. The scatter plots of the hindcasted versus measured peak wave periods for the different wave buoys are shown in Figure 5a-h (Figure 5a,e is for the Longdong buoy, Figure 5b,f is for the Hualien buoy, Figure 5c,g is for the Hsinchu buoy, and Figure 5d,h is for the Xiaoliuqiu buoy). The differences in the peak wave period reached 4-6 s at the Longdong and Hualien buoys due to the presence of CFSR winds in SCHISM-WWM-III (Figure 5a,b), while introducing the CFSR_5km wind fields into SCHISM-WWM-III led to peak wave period differences of only ±2 s. The differences in the wave direction between the hindcasts and measurements using different wind boundary conditions are illustrated in Figure 6a-h ( Figure 6a,e is for the Longdong buoy, Figure 6b,f is for the Hualien buoy, Figure 6c,g is for the Hsinchu buoy, and Figure 6d,h is for the Xiaoliuqiu buoy). These scatter plots show that the wave directions hindcasted by the CFSR_5km winds are more reliable than those hindcasted by the CFSR winds. Tables 1-3 summarize the percentages of the model-data differences for different classes and buoys owing to the adoption of winds from CFSR and CFSR_5km. According to the statistical results from these tables, the combination of SCHISM-WWM-III and winds from CFSR_5km exhibits a superior overall hindcasting performance for the SWHs, peak wave periods and wave directions. Thus, the CFSR_5km winds served as the meteorological forcing of SCHISM-WWM-III for the long-term wave parameter hindcasting in the present study.

Spatial Distribution of the Annual Average Wave Power
SCHISM-WWM-III outputted the hourly SWH and peak wave period in each grid, and then the WPD was computed by means of Equation (11). The spatial distributions of the offshore annual average WPD (in kW/m) in the Taiwanese waters from 1991 to 2000 and from 2001 to 2010 are shown in Figures 7 and 8, respectively. The most energetic sea areas are located along the southern and southeastern waters of Taiwan, with the annual mean WPD reaching 28-30 kW/m in 1991 (Figure 7). The sea area in the Central Taiwan Strait (i.e., the Penghu Islands) is another sea area rich in wave power. The highest annual mean WPD was 26 kW/m in 2000 (Figure 7). The annual average WPD was more abundant in 1991,1992,1994,1996,1999,2000,2004,2005,2007 and 2008, while 2002 had the smallest annual mean WPD (Figure 8). The offshore WPD in each grid was averaged over 20 years from 1991 to 2010, and the corresponding spatial distribution is presented in Figure 9. The offshore areas lying between Green Island and Orchid Island (Figure 3b), southeast of Orchid Island (Figure 3c), south of the Hengchun Peninsula (Figure 3d), and north of the Penghu Islands (as shown in Figure 3e) were found to be more energetic, with long-term average WPDs of 20, 22, 24 and 22 kW/m, respectively.

Annual and Monthly Average Incidence for the Exploitable Wave Power
Due to the purchase and deployment costs of wave energy converters, the sea area is worthy of exploitation only when the WPD exceeds 2 kW/m [14,[60][61][62]. To guarantee this requirement, only WPDs equal to or greater than 4 kW/m were included in the calculation of the incidence of exploitable wave power. The spatial distribution of the incidence (in percentages) of the WPD exceeding 4 kW/m in the offshore waters of Taiwan averaged over 1991 to 2010 is depicted in Figure 11. The highest incidences appear to be 70%-90% southeast of Taiwan, while the lowest incidence is only 30% in the nearshore southwestern waters off Taiwan. The sea area around the northern Penghu Islands exhibited a 50% annual mean WPD incidence. Similarly, Figure 12 describes the monthly average incidence for the WPD exceeding 4 kW/m in the waters surrounding Taiwan from 1991 to 2010. The monthly mean incidences along the eastern waters of Taiwan are approximately 100% in November, and they are always greater than 50% along the southeastern offshore areas every month. In July, the minimum monthly average incidence of 10% occurred in the Central Taiwan Strait; however, the monthly average incidences reached 80% in October, November, December, and January in the same sea area. mean WPD was 26 kW/m in 2000 (Figure 7). The annual average WPD was more abundant in 1991, 1992, 1994, 1996, 1999, 2000, 2004, 2005, 2007 and 2008, while 2002 had the smallest annual mean WPD (Figure 8). The offshore WPD in each grid was averaged over 20 years from 1991 to 2010, and the corresponding spatial distribution is presented in Figure 9. The offshore areas lying between Green Island and Orchid Island (Figure 3b), southeast of Orchid Island (Figure 3c), south of the Hengchun Peninsula (Figure 3d), and north of the Penghu Islands (as shown in Figure 3e) were found to be more energetic, with long-term average WPDs of 20, 22, 24 and 22 kW/m, respectively.

Annual and Monthly Average Incidence for the Exploitable Wave Power
Due to the purchase and deployment costs of wave energy converters, the sea area is worthy of exploitation only when the WPD exceeds 2 kW/m [14,[60][61][62]. To guarantee this requirement, only WPDs equal to or greater than 4 kW/m were included in the calculation of the incidence of exploitable wave power. The spatial distribution of the incidence (in percentages) of the WPD exceeding 4 kW/m in the offshore waters of Taiwan averaged over 1991 to 2010 is depicted in Figure 11. The highest incidences appear to be 70%-90% southeast of Taiwan, while the lowest incidence is only 30% in the nearshore southwestern

Variability of Wave Power
Sea areas with a steady and constant WPD are more appropriate for the deployment of wave energy converters than WPD areas that are more energetic but unstable [63]. To evaluate the temporal variability of the WPD in the offshore waters of Taiwan, the MVI and COV were calculated via Equations (13) and (14), respectively. The offshore areas displayed a relatively low MVI, and the COV presented weaker temporal variability; therefore, these offshore areas are considered to be potential zones for installing wave energy converters [64,65]. The spatial distribution of the MVI for the WPD in the offshore waters of Taiwan averaged over 1991 to 2010 is depicted in Figure 13. The MVI ranges between 1.8 and 2.7; the maximum MVI is located near the northern Penghu Islands (in the Central Taiwan Strait), and the minimum MVI appears in the northwestern waters of Taiwan, while an MVI of 2.4 is found along the southern and southeastern Taiwanese waters. Figure 14 exhibits the spatial distribution of the COV for the WPD in the waters surrounding Taiwan averaged over 1991 to 2010. The highest and lowest COVs are 2.5 and 1.5 in the southwestern nearshore and northwestern offshore waters of Taiwan, respectively. Moreover, a COV of 2.0 is discovered in the Central Taiwan Strait and the southern and southwestern offshore areas of Taiwan.

Variability of Wave Power
Sea areas with a steady and constant WPD are more appropriate for the deployment of wave energy converters than WPD areas that are more energetic but unstable [63]. To evaluate the temporal variability of the WPD in the offshore waters of Taiwan, the MVI and COV were calculated via Equations (13) and (14), respectively. The offshore areas displayed a relatively low MVI, and the COV presented weaker temporal variability; therefore, these offshore areas are considered to be potential zones for installing wave energy converters [64,65]. The spatial distribution of the MVI for the WPD in the offshore waters of Taiwan averaged over 1991 to 2010 is depicted in Figure 13. The MVI ranges between 1.8 and 2.7; the maximum MVI is located near the northern Penghu Islands (in the Central Taiwan Strait), and the minimum MVI appears in the northwestern waters of Taiwan, while an MVI of 2.4 is found along the southern and southeastern Taiwanese waters. Figure 14 exhibits the spatial distribution of the COV for the WPD in the waters surrounding Taiwan averaged over 1991 to 2010. The highest and lowest COVs are 2.5 and 1.5 in the southwestern nearshore and northwestern offshore waters of Taiwan, respectively. Moreover, a COV of 2.0 is discovered in the Central Taiwan Strait and the southern and southwestern offshore areas of Taiwan.

Determination of Optimal Wave Energy Hotspots
The hotspots for deploying wave energy converters were determined according to the assessments of the annual mean WPD, incidence of the exploitable WPD, and MVI and COV of the WPD, as well as the distance between the coast and wave energy conversion devices. Four offshore regions located between Green Island and Orchid Island (hereafter OH_1, as shown in Figure 3b), southeast of Orchid Island (hereafter OH_2, as shown in Figure 3c), south of the Hengchun Peninsula (hereafter OH_3, as shown in Figure 3d), and north of the Penghu Islands (hereafter OH_4, as shown in Figure 3e) were considered the optimal hotspots for the deployment of wave energy converters. These areas are consistent with the spatial distributions of relatively high WPD and incidence but lower MVI and COV. The monthly mean, maximum and minimum WPDs averaged over 1991 to 2010 for OH_1, OH_2, OH_3 and OH_4 are illustrated in Figure 15. The monthly maximum WPD exceeded 90 kW/m for OH_1, OH_2 and OH_3 in October and November and was approximately 100 kW/m for OH_4 in November. The monthly average WPD reached 30 kW/m-50 kW/m for OH_1, OH_2, OH_3 and OH_4 in October to December, while the monthly minimum WPD approached only 20 kW/m for the four optimal hotspots in December. Furthermore, the long-term directional distributions of the WPD were analyzed because the wave direction is an important parameter for the placement of a wave energy converter [35,66]. Figure 16 presents rose diagrams of the WPD for each of the four optimal wave energy hotspots averaged over 1991 to 2010. The dominant WPD directions lie between the east and northeast for OH_1, OH_2 and OH_3; however, the northeast is the prevailing WPD direction within OH_4. The characteristics of the selected optimal wave energy hotspots, such as the coordinates and water depths, are listed in Table 4.     converter [35,66]. Figure 16 presents rose diagrams of the WPD for each of the four o mal wave energy hotspots averaged over 1991 to 2010. The dominant WPD directions between the east and northeast for OH_1, OH_2 and OH_3; however, the northeast is prevailing WPD direction within OH_4. The characteristics of the selected optimal w energy hotspots, such as the coordinates and water depths, are listed in Table 4. 14, x FOR PEER REVIEW 22 o      "-" is the symbol for water depth below mean sea level; WP m is the annual mean wave power; annual total wave energy outputs; WE t is the annual total wave energy outputs; WE e is the exploitable annual wave energy output.

Discussion
The wave power is weak from March to June for OH_1, OH_2 and OH_3 and from April to August for OH_4, while the wave power is stronger from October to February for all four optimal hotspots. Additionally, the prevailing WPD directions for all optimal hotspots are nearly northeast. This phenomenon might be attributable to the prevalence of the northeast monsoon during late September to early May in Taiwan [67,68]. Alternatively, the long-term statistical data (from 1991 to 2019) issued by the CWB reveal that July, August, and September constitute the typhoon season in Taiwan; however, higher WPDs occur outside of the typhoon season. Although large waves are indeed induced by typhoons, the wave power in the offshore waters of Taiwan is driven mainly by the strong northeast monsoon. Moreover, the deeper water bathymetry also contributes to an increase in the WPD at these hotspots.
By applying Equation (12) with a given ∆t i equal to 8760 h, the annual total wave energy outputs (WE t ) for the four optimal hotspots were calculated, and the results are summarized in Table 4. WE t values of 184. 22, 198.15, 218.21, and 195.61 MWh/m were estimated for OH_1, OH_2, OH_3, and OH_4, respectively, which are more than double the results of a previous study in the same offshore waters [35,57]. The reason for such underestimated values is that the previous studies used wind field data with a temporal resolution of 6 h and a spatial resolution of 0.125 • from the ERA-Interim dataset, and the ERA-Interim global reanalysis is inferior to a regional high-resolution dynamical scaling atmospheric model in the performance of wind speeds. This is particularly true under severe weather conditions, e.g., strong monsoons or typhoons. In other words, the wave energy resources and the optimal hotspots for wave energy evaluated in the present study are more reliable than the results from previous studies [35,57]. The exploitable annual wave energy outputs (WE e ) for the hotspots were also assessed, where only ∆t i was replaced by the total hours corresponding to the WPD equal to or greater than 4 kW/m, i.e., 8760 h multiplied by an annual mean incidence. The estimations of WE e for OH_1, OH_2, OH_3, and OH_4 are 158.06, 182.89, 196.39, and 101.33 MWh/m, respectively (as listed in Table 4).

Summary and Conclusions
A three-layer nested Weather Research and Forecast (WRF) model with 45 km, 15 km and 5 km spatial resolutions for Domain 1, Domain 2, and Domain 3, respectively, was implemented over Taiwan and its offshore islands. The high-resolution WRF model was used to downscale the Climate Forecast System Reanalysis (CFSR) product dynamically. The hourly winds 10 m above sea level derived from the WRF model (CFSR_5km) and extracted from CFSR during the period of 2010 served as the meteorological boundary conditions for SCHISM-WWM-III, a fully coupled wave-circulation model. The validation results revealed that the hindcasting performance of SCHISM-WWM-III for the wave parameters, such as significant wave height (SWH), peak wave period and wave direction, by employing the winds from CFSR_5km was superior to those adopting the CFSR wind field. Therefore, the long-term sea states covering the period from 1991 to 2010 were hindcasted through the combination of CFSR_5km winds with SCHISM-WWM-III for assessing the offshore wave energy resources in Taiwan. The spatial distribution of the annual wave power density (WPD) averaged over twenty years shows that the offshore areas between Green Island and Orchid Island (OH_1), southeast of Orchid Island (OH_2), south of the Hengchun Peninsula (OH_3), and north of the Penghu Islands (OH_4) are more energetic than other regions, with long-term annual average WPDs of 20, 22, 24 and 22 kW/m, respectively. An analysis of the monthly mean WPDs over twenty years showed that the greatest WPD is observed from October to January, which means that the wave power is driven mainly by the northeast monsoon in the surrounding waters of Taiwan, even though typhoons can induce extreme waves. The annual average WPD, incidence of the hourly WPD greater than or equal to 4 kW/m, monthly variability index (MVI) and coefficient of variation (COV) for the WPD indicate that sea areas OH_1, OH_2, OH_3 and OH_4 can be considered appropriate sites for deploying wave energy converters. The exploitable long-term annual average total wave energy outputs are estimated to be 158.06, 182.89, 196.39, and 101.33 MWh/m for hotspots OH_1, OH_2, OH_3, and OH_4, respectively. Additionally, northeast is the primary wave direction for each of the four hotspots. The approach proposed by the present study could be applied elsewhere in the world, especially for a regional-scale study of wave power resources. As for evaluating the global wave energy resources, the reanalysis product is sufficient to reflect the distribution of wave energy worldwide.