Wind Characteristics in the Taiwan Strait: A Case Study of the First O ﬀ shore Wind Farm in Taiwan

: This study analyzed the wind speed data of the met mast in the ﬁrst commercial-scale o ﬀ shore wind farm of Taiwan from May 2017 to April 2018. The mean wind speed and standard deviation, wind rose, histogram, wind speed profile, and diurnal variation of wind speed with associated changes in wind direction revealed some noteworthy findings. First, the standard deviation of the corresponding mean wind speed is somewhat high. Second, the Hellmann exponent is as low as 0.05. Third, afternoons in winter and nights and early mornings in summer have the highest and lowest wind speed in a year, respectively. Regarding the histogram, the distribution probability of wind is bimodal, which can be depicted as a mixture of two gamma distributions. In addition, the corresponding change between the hourly mean wind speed and wind direction revealed that the land–sea breeze plays a signiﬁcant role in wind speed distribution, wind proﬁle, and wind energy production. The low Hellmann exponent is discussed in detail. To further clarify the e ﬀ ect of the land–sea breeze for facilitating future wind energy development in Taiwan, we propose some recommendations.


Introduction
Wind energy is a sustainable renewable energy source with lower impact on the environment than fossil fuel burning. Nowadays, wind power is becoming increasingly ubiquitous in countries. In 2018, wind power accounted for 4.8% of global electricity. In 2019, the total installed capacity of wind turbines increased by 10% to 620 gigawatts (GW), of which the total installed capacity of offshore wind power was 27 GW, with a net increase of 5 GW (19%) [1,2]. The Taiwan Strait accounts for most of the top 20 offshore wind farm sites in the world, making it an excellent location for developing offshore wind power. The primary goal of Taiwan is to raise its offshore wind power capacity to 5.7 GW by 2025. Furthermore, for the next 10 years from 2026 to 2035, the goal is to install 1 GW per year [3,4].
Wind energy is an intermittent energy source. Although the annual power generated is roughly the same, it can vary greatly on shorter time scales, such as hourly, daily, or quarterly. The power supply system must always maintain a balance between instantaneous power generation and consumption to protect grid stability; hence, this variability introduces substantial challenges for integrating a large amount of wind energy into the grid system. Therefore, in addition to forecasting wind energy after operating wind farms, assessing wind energy potential before developing wind farms is essential.
insightful technical details for further wind energy development in Taiwan. Specific objectives of this study include: (1) characterizing and modeling the probability distribution of the at-site wind speed, (2) characterizing the seasonal variation of wind speed and wind direction, (3) investigating the major climatological factors that affect the at-site wind speed characteristics and the performance of the Taiwan Strait wind farms.
Our study analyzed the wind speed data of the mast of the Formosa 1 wind farm to first determine the wind characteristics of the wind farm to obtain the appropriate probability distribution function for preparing evaluation indexes of wind farms, such as annual energy production (AEP) and operating probability. The wind profile formed through the combination of wind speeds at different heights collected from field observation can then be used to calculate the blade stress and can serve as the basis for blade design and selection [18][19][20]. For evaluating the strength of land-sea breeze in the coastal area of Taiwan, indicators have been recommended, which can also be used as factors for site selection and wind energy forecast. Overall, the results may contribute some useful insights regarding offshore wind farm development in the Taiwan Strait in the future.

Data Source and Analysis Methods
Data used in this study were 10-min average wind speed measured from May 2017 to April 2018 by the met mast of the first commercial-scale offshore wind farm in Taiwan-Formosa I (operating from December 2019 to date). The met mast is located approximately 5 km away from Zhunan Township, Miaoli County (Figure 1). Anemometers and wind vanes are installed on the met mast at 10, 30, 50, and 90 m above the mean sea level. Winds highly vary spatially and temporally. Thus, understanding wind characteristics is essential for developing wind energy. Wind characteristics include the wind speed and direction. Furthermore, wind energy is related to the wind speed cube. Therefore, only through correctly estimating the wind speed can the wind energy potential be appropriately evaluated. Consequently, the probability density function is generally used to interpret the variation of wind speed with time and space. Regarding the wind direction, wind rose is often adopted to express its associated changes with the wind speed. Nevertheless, the most necessary and straightforward procedure is to evaluate the mean and standard deviation. These statistical methods and preliminary analysis outcomes are briefly described in this section. In addition, the common land-sea breeze phenomenon in the coastal area is described. Winds highly vary spatially and temporally. Thus, understanding wind characteristics is essential for developing wind energy. Wind characteristics include the wind speed and direction. Furthermore, wind energy is related to the wind speed cube. Therefore, only through correctly estimating the wind speed can the wind energy potential be appropriately evaluated. Consequently, the probability density function is generally used to interpret the variation of wind speed with time and space. Regarding the wind direction, wind rose is often adopted to express its associated changes with the wind speed. Nevertheless, the most necessary and straightforward procedure is to evaluate the mean and standard Energies 2020, 13, 6492 4 of 21 deviation. These statistical methods and preliminary analysis outcomes are briefly described in this section. In addition, the common land-sea breeze phenomenon in the coastal area is described.

Annual Mean and Standard Deviation of Wind Speed
The annual mean and standard deviation of 10-min mean wind speed are shown in Figure 2. The average wind speed at each height is almost the same (the annual average wind speed difference between the heights of 90 and 10 m is <1 m/s). Moreover, the Hellmann exponent (a) of the vertical profile of the horizontal wind described by power law is 0.05, which is less than that recommended by design codes or studies [21][22][23][24][25][26], and the spreading (standard deviation) of wind speed is quite large. Power law can be expressed as follows: where z r is the reference height and a is the power law exponent.

Annual Mean and Standard Deviation of Wind Speed
The annual mean and standard deviation of 10-min mean wind speed are shown in Figure 2. The average wind speed at each height is almost the same (the annual average wind speed difference between the heights of 90 and 10 m is < 1 m/s). Moreover, the Hellmann exponent (a) of the vertical profile of the horizontal wind described by power law is 0.05, which is less than that recommended by design codes or studies [21][22][23][24][25][26], and the spreading (standard deviation) of wind speed is quite large. Power law can be expressed as follows: where zr is the reference height and a is the power law exponent.

Probability Density Function of Wind Speed
Wind speed data have been modeled using different probability distribution functions including the Weibull, log-normal, gamma, Burr, Gumbel, and Rayleigh distributions [27][28][29][30][31][32][33][34]. Among these distributions, the Weibull distribution appears to be most widely used. In addition to the univariate modeling of wind speed data, various mixture distributions have also been used when the wind speed data exhibited a bimodal pattern [28,32,35]. Thus, in this study, the Weibull, mixture Weibull, and mixture gamma distributions were considered for probability density modeling of the wind speed data in Taiwan. The properties of the Weibull and gamma distributions and the maximum likelihood estimators of their parameters are briefly described in the following subsections.

Weibull Distribution
The probability density function and cumulative distribution function of a random variable X with the Weibull distribution are respectively given by where and are the shape and scale parameters, respectively, and > 0, > 0 and > 0. The mean and variance of the Weibull distribution can also be expressed as the following functions of the scale and shape parameters,

Probability Density Function of Wind Speed
Wind speed data have been modeled using different probability distribution functions including the Weibull, log-normal, gamma, Burr, Gumbel, and Rayleigh distributions [27][28][29][30][31][32][33][34]. Among these distributions, the Weibull distribution appears to be most widely used. In addition to the univariate modeling of wind speed data, various mixture distributions have also been used when the wind speed data exhibited a bimodal pattern [28,32,35]. Thus, in this study, the Weibull, mixture Weibull, and mixture gamma distributions were considered for probability density modeling of the wind speed data in Taiwan. The properties of the Weibull and gamma distributions and the maximum likelihood estimators of their parameters are briefly described in the following subsections.

Weibull Distribution
The probability density function and cumulative distribution function of a random variable X with the Weibull distribution are respectively given by where α and β are the shape and scale parameters, respectively, and x > 0, α > 0 and β > 0. The mean and variance of the Weibull distribution can also be expressed as the following functions of the scale and shape parameters, where Γ represents the gamma function. The scale and shape parameters can only be estimated using an observed sample. The method of moments (MOM) and the method of maximum likelihood (MML) are widely applied for parameters estimation. The former estimates the parameters of the underlying distribution by using the relationship between the population moments and the parameters and substituting the sample moments for the population moments. The latter maximizes the likelihood of observing the data, given the model. Generally speaking, the MOM estimators can be easily calculated, although they are often biased. By contrast, the maximum likelihood estimators (MLE) are more computationally intensive and often require numeric iterations. However, the maximum likelihood estimators are asymptotically normal and unbiased and have lower variance than other methods.
Given a random sample of size n, i.e., x 1 , x 2 , · · · , x n , the MLE of the shape and scale parameters of the Weibull distribution are as follows:β The above MLE can only be estimated by numerical iteration.

Gamma Distribution
The probability density function and cumulative distribution function of the gamma distribution are respectively given by where γ is the lower incomplete gamma function. The mean and variance of the gamma distribution can also be expressed as the following functions of the scale and shape parameters, Given a random sample of size n, i.e., x 1 , x 2 , · · · , x n , the MLE of the shape and scale parameters of the gamma distribution are as follows [36,37]: where x n represents the sample mean and U = lnx n − ln (x 1 x 2 · · · x n ) 1 n .

Mixture Distribution
The necessity of mixture distribution modeling arises from the situation when one samples from a population that consists of several homogeneous subpopulations, which are called the components of the population. For a k-component mixture distribution, let f x; ξ j represents the probability density function of the j-th subpopulation, with a parameter vector ξ j , then the probability density function of the mixture population can be expressed by where π j is the proportion of the total population that is in the j-th component and is called the component weight. In the above equation, π = (π 1 , π 2 , · · · , π k ) and φ = (ξ 1 , ξ 2 , · · · , ξ k ) together form the set of parameters of the mixture distribution. The weights satisfy the constraints that π j ≥ 0 and Given a random sample of size n, i.e., x 1 , x 2 , · · · , x n , of a k-component mixture distribution, the maximum likelihood estimates,π andφ, are the parameter values that maximize the log-likelihood ( ) of the mixture distribution, i.e., The expectation-maximization (EM) algorithm has gained great popularity in calculating the maximum likelihood or maximum a posteriori (MAP) estimates of the mixture distribution parameters. The EM algorithm is an iterative numeric approach which involves two steps. The expectation (E) step assigns each observation to each subpopulation by using the current estimates of the parameters according to the relative density of the data points under each subpopulation. The maximization (M) step updates the estimates of the parameters by considering the expectation of a set of latent variables that represent the values of the parameter ξ corresponding to individually sampled components. Readers are referred to Lindsay [38] and Hastie et al. [39] for details of the EM algorithm. Furthermore, Figure 3a shows that the prevailing wind is obvious, mainly in the direction of 30 • , and the wind from the north-easterly direction (15-45 • ) accounts for most of the wind for approximately 40% of the year. The secondary direction is 215 • , and the south-westerly wind (210-240 • ) accounts for most of the wind for approximately 20% of the year. The black portion in the figure represents the wind speed <4 m/s. Usually, wind turbines do not operate within this range. The blue portion shows the wind speed at 4-12 m/s; wind turbine blades begin to rotate but do not reach the full load. The yellow portion indicates that the wind speed is 12-25 m/s. As a rule, wind turbines adjust the pitch angle of their blades to maintain the rated load (full load). The red portion indicates that the wind speed has exceeded 25 m/s, and it is the point at which wind turbines normally shut down. According to the wind speed recorded at 90-m height, 147 times in a year, the 10-min mean wind speed exceeds 25 m/s, most of which occurs in October and November, when the north-east monsoon is rather strong, and in summer, only when typhoons strike.

Land-Sea Breeze
Due to the varying land-sea thermal capacity of absorbing solar radiation energy, the mesoscale thermal circulation near the coast (as shown in Figure 4), which is called the land-sea breeze, occurs at many coastal locations throughout the world. Coupled with the different land-sea surface friction values and the effect of synoptic pressure gradients, the local wind circulation becomes more complex.
(a) (b)  Figure 3b shows the wind rose of December. The wind direction is stable. In this month, 60% and 15% of 30 • and 45 • wind is present, respectively. In addition to the northeast wind direction, there is hardly any wind in other directions. Figure 3c shows the wind rose in June. Most of the wind comes from the southwest (210-240 • ). Figure 3d shows the wind rose in September. The northeast monsoon begins to prevail, but the southwest monsoon is still strong.

Land-Sea Breeze
Due to the varying land-sea thermal capacity of absorbing solar radiation energy, the mesoscale thermal circulation near the coast (as shown in Figure 4), which is called the land-sea breeze, occurs at many coastal locations throughout the world. Coupled with the different land-sea surface friction values and the effect of synoptic pressure gradients, the local wind circulation becomes more complex. Due to the varying land-sea thermal capacity of absorbing solar radiation energy, the mesoscale thermal circulation near the coast (as shown in Figure 4), which is called the land-sea breeze, occurs at many coastal locations throughout the world. Coupled with the different land-sea surface friction values and the effect of synoptic pressure gradients, the local wind circulation becomes more complex.

Sea Breeze
Sea breeze occurs during daytime and blows from the sea to the land. This is mainly because the land surface and the air above are heated by the sun, raising the air due to its low density. Moreover, because the heat storage capacity of seawater is high, the air temperature above it increases to a lesser extent, resulting in a density difference between the air above the seawater and that above the land. Consequently, the wind blows from the sea to the land. The sea breeze has the following characteristics: • The sea breeze speed might reach 10-11 m/s at 100-m height [40].

•
The height of the return flow of sea breeze ranges from hundreds of meters to 1 or 2 km.

•
The wind direction of the sea breeze is constantly affected by the Coriolis force, which rotates clockwise and counterclockwise in the northern and southern hemispheres, respectively.

•
The ideal land-sea breeze hodograph can be represented by an ellipse.
According to different synoptic gradients, the sea breeze can be divided into the following three models [41]: 1.
Pure sea breeze The need for a synoptic gradient is less. Sea breeze formation relies mainly on the temperature difference between air and seawater.

2.
Corkscrew sea breeze This is the strongest sea breeze model. The synoptic gradient has an offshore component, favoring sea breeze formation. For example, according to Buys Ballot's law, the formation of corkscrew sea breeze in Taiwan is favorable when the northeast monsoon is prevailing in winter.

Backdoor sea breezes
The synoptic gradient has an onshore component, hindering the sea breeze formation. This is the case in Taiwan summer as the south-west monsoon is prevailing.

Land Breeze
During the nighttime, because the heat storage capacity of the land decreases, heat dissipation increases. The air above the land surface is denser than that above the sea surface. Therefore, the wind blows toward the sea, resulting in land breeze. In the past, sailors often sailed out of the port in the morning by using land breeze.

Results
The wind characteristics of Formosa I offshore wind farm located in the Taiwan Strait can be roughly understood from the annual mean and standard deviation of wind speed and the wind rose. The average wind speed at 90-m height is approximately 8 m/s, but the standard deviation is >5.5 m/s.

•
The exponent of the vertical profile of horizontal wind is approximately 0.05, which is relatively low.

•
The prevailing wind direction is obvious and concentrated and strongly correlated with the season. The prevailing wind directions in winter and summer are almost opposite.

•
Wind turbines may be fully loaded (wind speed > 12 m/s) > 30% of the year.
We clearly understand the phenomenon that the standard deviation is somewhat high to the corresponding mean wind speed, and there is a notable amount of full load (>12 m/s) time for the moderate mean wind speed (<8 m/s). The monthly mean and standard deviation of wind speed and the hourly mean wind speed and direction in a month were further investigated.

Monthly Mean and Standard Deviation of Wind Speed
The annual mean wind speed at 90-m height approximately equals the mean summer wind speed ( Figure 5). However, the mean wind speed in winter can increase to 12 m/s with an even higher standard deviation, yielding numerous high-speed winds and light winds.

Monthly Mean and Standard Deviation of Wind Speed
The annual mean wind speed at 90-m height approximately equals the mean summer wind speed ( Figure 5). However, the mean wind speed in winter can increase to 12 m/s with an even higher standard deviation, yielding numerous high-speed winds and light winds. The fluctuation of the monthly mean wind speed in May, June, July, and August is caused by the variation in the southwest monsoon intensity. Similarly, changes in the northeast monsoon intensity cause the monthly mean wind speed to oscillate in October, November, and December. However, to conclude whether these phenomena occur every year, more time series data are needed for investigation.

Hourly Mean Wind Speed and Direction in a Month
Although monthly means and standard deviations of wind speed help to understand the seasonal effect on the wind speed, hourly variation in the wind speed in each month provides diurnal changes. Due to space limitations, only the results of June and December are shown, representing summer and winter, respectively. December has the highest monthly mean wind speed in a year. June is halfway from December and has no special incidents (typhoons strike). Two typhoons (Typhoons Nesat and Haitang) passed by the Formosa I offshore wind farm as they hit Taiwan in July 2017.
The highest wind speed of a day occurs in the afternoon both in winter and summer, and the highest wind speed of a year occurs in the winter afternoon ( Figure 6). Moreover, the lowest wind speed in a year often occurs at night and early morning in the summer. The daily mean angular variation of wind direction in summer is approximately 40° compared with that in winter, which is only approximately 25°. The fluctuation of the monthly mean wind speed in May, June, July, and August is caused by the variation in the southwest monsoon intensity. Similarly, changes in the northeast monsoon intensity cause the monthly mean wind speed to oscillate in October, November, and December. However, to conclude whether these phenomena occur every year, more time series data are needed for investigation.

Hourly Mean Wind Speed and Direction in a Month
Although monthly means and standard deviations of wind speed help to understand the seasonal effect on the wind speed, hourly variation in the wind speed in each month provides diurnal changes. Due to space limitations, only the results of June and December are shown, representing summer and winter, respectively. December has the highest monthly mean wind speed in a year. June is halfway from December and has no special incidents (typhoons strike). Two typhoons (Typhoons Nesat and Haitang) passed by the Formosa I offshore wind farm as they hit Taiwan in July 2017.
The highest wind speed of a day occurs in the afternoon both in winter and summer, and the highest wind speed of a year occurs in the winter afternoon ( Figure 6). Moreover, the lowest wind speed in a year often occurs at night and early morning in the summer. The daily mean angular variation of wind direction in summer is approximately 40 • compared with that in winter, which is only approximately 25 • .

Probability Distribution Modeling of Wind Speed
Wind speed at different heights can be fit to certain probability density functions to fully describe its statistical characteristics and be used for wind power potential assessment. Figure 7 illustrates the relative frequency histograms and the fitted Weibull density functions of the wind speed at different heights of the met mast. Wind speed data were fit to a Weibull probability density function owing to its high popularity for depicting wind speed distribution.
For the wind speed data at the 90-m height (see Figure 7a), the maximum likelihood estimates of the scale and shape parameters are 8.75 and 1.44, respectively. Although the fitted distribution yields expected value (7.95) and variance (31.45) very close to the sample mean (7.94) and sample variance (30.77), it is also evident that the fitted Weibull density is lower than the sample density in the low wind speed range (1 m/s to 4 m/s) and high wind speed range (12 m/s to 16 m/s), and is higher than the sample density in the median wind speed range (5 m/s to 10 m/s). A further check on the coefficient of skewness found that the coefficient of skewness of the observed wind speed data (0.696) is much lower than that of the fitted Weibull density function (1.148). Similar patterns of the sample wind speed distribution and fitted Weibull density function can also be observed for wind speed data at other heights. Such mismatch suggests that the univariate Weibull density function is not appropriate for characterizing the wind speed data in our study site.

Probability Distribution Modeling of Wind Speed
Wind speed at different heights can be fit to certain probability density functions to fully describe its statistical characteristics and be used for wind power potential assessment. Figure 7 illustrates the relative frequency histograms and the fitted Weibull density functions of the wind speed at different heights of the met mast. Wind speed data were fit to a Weibull probability density function owing to its high popularity for depicting wind speed distribution.
For the wind speed data at the 90-m height (see Figure 7a), the maximum likelihood estimates of the scale and shape parameters are 8.75 and 1.44, respectively. Although the fitted distribution yields expected value (7.95) and variance (31.45) very close to the sample mean (7.94) and sample variance (30.77), it is also evident that the fitted Weibull density is lower than the sample density in the low wind speed range (1 m/s to 4 m/s) and high wind speed range (12 m/s to 16 m/s), and is higher than the sample density in the median wind speed range (5 m/s to 10 m/s). A further check on the coefficient of skewness found that the coefficient of skewness of the observed wind speed data (0.696) is much lower than that of the fitted Weibull density function (1.148). Similar patterns of the sample wind speed distribution and fitted Weibull density function can also be observed for wind speed data at   Table 1. A more detailed examination of the sample relative frequency histogram indicates that the wind speed data may have a bimodal distribution with one peak near 3 m/s and another near 13 m/s. Thus, the two-component mixture Weibull distribution and mixture gamma distribution were also fitted to wind speed data of different heights. Figure 8 demonstrates the fitted probability density function of the mixture gamma, mixture Weibull, and univariate Weibull distributions. Both the mixture gamma distribution and mixture Weibull distribution well represent the bimodal pattern of the wind speed data. Parameters of the fitted mixture gamma and mixture Weibull distributions are summarized in Table 2.  Table 1. A more detailed examination of the sample relative frequency histogram indicates that the wind speed data may have a bimodal distribution with one peak near 3 m/s and another near 13 m/s. Thus, the two-component mixture Weibull distribution and mixture gamma distribution were also fitted to wind speed data of different heights. Figure 8 demonstrates the fitted probability density function of the mixture gamma, mixture Weibull, and univariate Weibull distributions. Both the mixture gamma distribution and mixture Weibull distribution well represent the bimodal pattern of the wind speed data. Parameters of the fitted mixture gamma and mixture Weibull distributions are summarized in Table 2.

Model Selection Based on the Akaike Information Criterion
In order to select the best fit distribution among the univariate Weibull, the mixture Weibull, and the mixture gamma distribution, the Akaike information criterion (AIC) was adopted in this study.
A random sample of size n, i.e., , , ⋯ , , can be fitted by different probability distributions with estimated parameters . The AIC of the fitted distribution of the given sample is calculated by [42] AIC = 2 − 2ℓ (17)

Model Selection Based on the Akaike Information Criterion
In order to select the best fit distribution among the univariate Weibull, the mixture Weibull, and the mixture gamma distribution, the Akaike information criterion (AIC) was adopted in this study.
A random sample of size n, i.e., x 1 , x 2 , · · · , x n , can be fitted by different probability distributions with estimated parametersθ. The AIC of the fitted distribution of the given sample is calculated by [42]  where p is the number of parameters of the fitted distribution and θ is the log-likelihood of the random sample under the fitted distribution (See Equation (16)). The best fit model yields the lowest AIC values. Table 3 shows the AIC values of the fitted distributions (univariate Weibull, mixture Weibull, and mixture gamma) for wind speed data at different heights. The two mixture distribution models clearly outperform the univariate Weibull distribution and the mixture gamma distribution yields the lowest AIC value for wind speed data at all different heights.

The Interaction between Land-Sea Breezes and Monsoons
The analysis of daily variation in wind speed and direction indicates the following features: • Daily wind direction oscillation basically centers on the prevailing wind direction.

•
The stronger the prevailing wind speed is, the smaller the variation in wind direction.
In addition, considering the bimodal probability density function, the distance between the met mast and the coastline, and the land-sea breeze cyclically rotating in a day, the wind field, consisted of resultant wind vectors, appears to be highly affected by the monsoon and land-sea breeze.

In Winter
Taiwan has a noticeable northeast monsoon in winter. Generally, the wind speed on the ground could surpass Beaufort number 5 (or 8.0-10.7 m/s) and even reach Beaufort number 10 (or 24.5-28.4 m/s). The following are an idealized land-sea breeze hodograph and a northeast monsoon (with a wind orientation of 30 • ), with the illustrations of daily variations in the mean wind speed and wind direction for the month. We used December 2017 as an example (Figure 6b). The angle between the coastline nearby the met mast and the north is approximately 35-40 • (Figure 1). Schematics are shown with invariable monsoon and ellipse. The following should at least be considered in a real hodograph:

•
The wind speed and direction of the northeast monsoon are not constant and change with time.

•
The speed intensity of land-sea breeze wind varies due to the temperature difference between the sea surface and the ground. For example, in winter the sunshine is weaker and shorter; thus, the ground temperature is low in the early morning and evening, hence the sea breeze is usually weaker in winter than that in summer.

•
Changes in the wind direction caused by terrain.

•
The speed of the land breeze during nighttime is frequently lower than that of the sea breeze during daytime. Therefore, the minor axis of the night part of the ellipse is shorter than that of the ideal ellipse.
In Figure 9, the blue arrow represents the northeast monsoon in winter, and the land-sea breeze arrow rotates in the ellipse, with different colors for different times of the day. The gray arrow represents the wind speed resultant vector composed of the land-sea breeze and the winter monsoon. First, Figure 9a shows that the angle of the wind direction at 9 am is close to the maximum angle of the day (approximately 40 • ). Then, Figure 9b shows that the wind speed reaches the maximum value of the whole day at 3 pm, and the angle of the wind direction also gradually rises from the minimum value. Figure 9c shows that the wind speed is still strong at 9 pm, but the wind direction returns to the maximum value of the whole day.

In Summer
The southwest monsoon often occurs in summer in Taiwan, and the wind speed is usually < 8.0 m/s. When monsoon prevails, the wind speed can reach an approximate Beaufort number 5 (or 8.0-10.7 m/s) but rarely exceeds Beaufort number 6 (or 10.8-13.8 m/s). We used June 2017 as an example (Figure 6a).
In Figure 10, the red arrow represents the summer southwest monsoon, and the gray arrow represents the wind speed resultant vector composed of the land-sea breeze and the summer monsoon. First, Figure 10a shows that the angle of the wind direction at 7 am is close to the minimum angle of the day (approximately 210°). Next, Figure 10b shows that the wind speed reaches the maximum value of the whole day at 2 pm, when the angle of the wind speed reaches the maximum value of the whole day (approximately 260°). Finally, Figure 10c shows that when the wind direction has returned to nearly parallel to the coastline at 7 pm, when the wind speed drops and approaches the minimum value.

In Summer
The southwest monsoon often occurs in summer in Taiwan, and the wind speed is usually <8.0 m/s. When monsoon prevails, the wind speed can reach an approximate Beaufort number 5 (or 8.0-10.7 m/s) but rarely exceeds Beaufort number 6 (or 10.8-13.8 m/s). We used June 2017 as an example (Figure 6a).
In Figure 10, the red arrow represents the summer southwest monsoon, and the gray arrow represents the wind speed resultant vector composed of the land-sea breeze and the summer monsoon. First, Figure 10a shows that the angle of the wind direction at 7 am is close to the minimum angle of the day (approximately 210 • ). Next, Figure 10b shows that the wind speed reaches the maximum value of the whole day at 2 pm, when the angle of the wind speed reaches the maximum value of the whole day  Figure 10c shows that when the wind direction has returned to nearly parallel to the coastline at 7 pm, when the wind speed drops and approaches the minimum value. On the basis of the above illustrations, the bimodal probability distribution function is probably induced by the land-sea breeze, which speeds up or slows down the prevailing monsoon.

Low Hellmann Exponent
The composition of monsoon and land-sea wind vectors also alters the shape of the wind profile, and the exponent of the power law changes accordingly. Particularly, when subtracting two wind vectors of comparable magnitude and the smaller one having a larger exponent, the resultant wind profile will have a lower exponent. This situation is common in summer (Figure 10c). For example, when the southwest monsoon wind speed along the coastline is 2 m/s and the exponent is 0.1, whereas sea breeze speed along the coastline toward the southwest is 1.2 m/s with 0.14 exponents, the exponent of the resultant wind profile is then 0.046. This may partly be the reason for the low exponent of the wind field. A detailed explanation regarding the low Hellmann exponent is provided in the Discussion section. On the basis of the above illustrations, the bimodal probability distribution function is probably induced by the land-sea breeze, which speeds up or slows down the prevailing monsoon.

Low Hellmann Exponent
The composition of monsoon and land-sea wind vectors also alters the shape of the wind profile, and the exponent of the power law changes accordingly. Particularly, when subtracting two wind vectors of comparable magnitude and the smaller one having a larger exponent, the resultant wind profile will have a lower exponent. This situation is common in summer (Figure 10c). For example, when the southwest monsoon wind speed along the coastline is 2 m/s and the exponent is 0.1, whereas sea breeze speed along the coastline toward the southwest is 1.2 m/s with 0.14 exponents, the exponent of the resultant wind profile is then 0.046. This may partly be the reason for the low exponent of the wind field. A detailed explanation regarding the low Hellmann exponent is provided in the Discussion section.

Discussion
We can further explore impacts of the above results on wind energy production and present recommendations regarding wind energy assessment, wind farm siting, wind turbine selection, and wind energy forecasting.

Annual Energy Production (AEP)
The AEP of a wind turbine is the total amount of electrical energy it produces over a year, which is measured in megawatt hours or gigawatt hours (MWh or GWh). Similarly, AEP can be described as an index for assessing the economic feasibility of a commercial wind turbine for electricity generation. Accordingly, this index can be estimated by multiplying the power for each wind speed from the power curve with the wind speed frequency distribution experienced by the wind turbine, and the number of hours in a year. (18) where AEP is the annual energy production. Pt(v) is the power output curve of a wind turbine. v ci is the cut-in wind speed (minimum wind speed for the wind turbine to operate), and v co is the cut-out wind speed (maximum wind speed before being turned off the wind turbine for equipment protection). However, power curves are proprietary of manufactures and are therefore difficult to obtain. Hereby, we used the power curve of an 8-MW reference wind turbine from the EU FP7 project LEANWIND [43] with the identical cut-in and cut-out speeds as that of the wind turbine deployed in Formosa 1 offshore wind farm. The power curve is illustrated in Figure 11.

Discussion
We can further explore impacts of the above results on wind energy production and present recommendations regarding wind energy assessment, wind farm siting, wind turbine selection, and wind energy forecasting.

Annual Energy Production (AEP)
The AEP of a wind turbine is the total amount of electrical energy it produces over a year, which is measured in megawatt hours or gigawatt hours (MWh or GWh). Similarly, AEP can be described as an index for assessing the economic feasibility of a commercial wind turbine for electricity generation. Accordingly, this index can be estimated by multiplying the power for each wind speed from the power curve with the wind speed frequency distribution experienced by the wind turbine, and the number of hours in a year.
where AEP is the annual energy production. Pt(v) is the power output curve of a wind turbine. vci is the cut-in wind speed (minimum wind speed for the wind turbine to operate), and vco is the cut-out wind speed (maximum wind speed before being turned off the wind turbine for equipment protection). However, power curves are proprietary of manufactures and are therefore difficult to obtain. Hereby, we used the power curve of an 8-MW reference wind turbine from the EU FP7 project LEANWIND [43] with the identical cut-in and cut-out speeds as that of the wind turbine deployed in Formosa 1 offshore wind farm. The power curve is illustrated in Figure 11. As soon as the AEP is determined, the capacity factor Cf can be calculated using the ratio of the AEP to the rated power output Pr over a period of an entire year as follows: In our case, the rated AEP is 70.08 GWh when the wind turbine works at full capacity for the entire year. As shown in Table 4, the AEP and capacity factor of the mixture gamma function are 28.73% and 40.99%, respectively, which is the best fitting function because of the closest approximation to the values calculated using the histogram. Univariate Weibull distribution function is undervalued by approximately 6.8%. Contrarily, the mixture Weibull distribution function is overrated by approximately 7.1%. As soon as the AEP is determined, the capacity factor C f can be calculated using the ratio of the AEP to the rated power output Pr over a period of an entire year as follows: In our case, the rated AEP is 70.08 GWh when the wind turbine works at full capacity for the entire year. As shown in Table 4, the AEP and capacity factor of the mixture gamma function are 28.73% and 40.99%, respectively, which is the best fitting function because of the closest approximation to the values Energies 2020, 13, 6492 17 of 21 calculated using the histogram. Univariate Weibull distribution function is undervalued by approximately 6.8%. Contrarily, the mixture Weibull distribution function is overrated by approximately 7.1%.

Operating Probability of Wind Turbines
The operating probability is the difference in probabilities between the cut-in and cut-out wind speeds. The calculation introduces merely the features of the probability density function, making it another indicator for selecting wind turbines based on economic feasibility.
The Formosa I offshore wind farm deploys 20 Siemens Gamesa 6-MW SWT-6.0-154 wind turbines. Its cut-in and cut-out wind speeds are 4.0 and 25.0 m/s, respectively. Substitute mixture distributions into the aforementioned equation for calculating the probability between the cut-in and cut-out wind speeds. Table 5 shows that the operating probability of the mixture gamma function is 68.62%, which is relatively close to the value obtained by a histogram. Contrary to the AEP results, univariate Weibull distribution function is overrated by approximately 6.8%, and the mixture Weibull distribution function is undervalued by approximately 1.5%.

Strength, Direction, Depth and Extension of Land-Sea Breeze
The magnitude and direction of the land-sea breeze affect the resultant wind speed that actually drives wind turbines, which determines the capacity of wind power production. The peak electricity consumption period in Taiwan is summer afternoons. If the sea breeze speed is insufficient or even absent (Figure 10b), acquiring the power produced by other energy sources in advance is necessary to avoid grid instability.
Land-sea breeze influences the speed gradient (i.e., exponent) of the resultant wind profile, accordingly affecting the shear stress imposed on wind turbine blades. Additionally, the extent to which the land-sea breeze depth encompasses the swept area of blades is related to the estimation of the shear stress on the blades.
The distance and strength of the land-sea breeze are both important factors for the site selection of offshore wind farms. Therefore, the strength, direction, depth, and extension distance of the land-sea breeze are the topics of in-depth study.

Hellmann Exponent
From the aforementioned analysis, the wind profile exponent is low, which may be caused by the following:

Composition of Monsoons and Land-Sea Winds
Under certain conditions, the composition of wind speeds will cause the Hellmann exponent to decrease, refer to Section 3.6 for details.

Atmospheric Instability above the Sea Surface
The Hsinchu buoy, set up by the Central Weather Bureau approximately 6 km northeast to the met mast, provides seawater and air temperatures from 12:00 a.m. to 9:00 a.m. The temperature difference between the air (Ta) and the seawater (Tw) in summer is smaller than that in winter ( Figure 12). However, even in the early spring mornings, the temperature of the seawater is higher than that of the air, which means that the ocean surface keeps heating the air.
Energies 2020, 13, x FOR PEER REVIEW 18 of 21 12). However, even in the early spring mornings, the temperature of the seawater is higher than that of the air, which means that the ocean surface keeps heating the air. shows that the met mast in this study is located nearby the narrowest section of the Taiwan Strait, by which the wind accelerated, leading to the internal momentum exchange.
The latter two reasons indicate that the frequent momentum exchange among the upper and lower layers exists inside the wind field, combined with the first reason, hereby making the Hellmann exponent low.

Wind Energy Forecast
For improving wind energy forecasting, the following factors are recommended:


Temperature difference between the sea surface and the ground to assist the assessment of the formation and strength of land-sea winds.  Temperature difference between the air above the sea and the sea surface to assist the evaluation of the possible changes in the vertical wind profile by the atmospheric stability.

Conclusions
This study analyzed the wind speed data of the met mast located at the Formosa I offshore wind farm from May 2017 to April 2018 to understand the wind characteristics of the Taiwan Strait:


The annual mean wind speed at 90-m height is approximately 8 m/s, but the standard deviation is > 5.5 m/s. The wind speed has a somewhat large spreading.  Two prevailing wind directions exist, which are almost opposite, namely the north-east monsoon in winter and the south-west monsoon in summer. The intensity of the north-east monsoon is usually greater than that of the southwest monsoon.  Winter afternoon has the highest wind speed of the year. Conversely, summer night and early morning have the lowest wind speed of the year.  The exponent of the vertical profile of the horizontal wind is low at approximately 0.05.
Compared with Weibull distribution, the histogram has a higher distribution probability for wind speeds < 4 and > 11 m/s and a lower distribution probability for wind speeds from 4 to 11 m/s. It appears to be a bimodal distribution. Therefore, a mixture distribution consisting of two distributions is used to fit the histogram. The mixture gamma distribution has the best performance for depicting wind data, which is confirmed by Sections 4.1 and 4.2. In addition to the diurnal  Figure 1a shows that the met mast in this study is located nearby the narrowest section of the Taiwan Strait, by which the wind accelerated, leading to the internal momentum exchange.

Channel Effect Caused by the Contraction of the Taiwan Strait
The latter two reasons indicate that the frequent momentum exchange among the upper and lower layers exists inside the wind field, combined with the first reason, hereby making the Hellmann exponent low.

Wind Energy Forecast
For improving wind energy forecasting, the following factors are recommended: • Temperature difference between the sea surface and the ground to assist the assessment of the formation and strength of land-sea winds. • Temperature difference between the air above the sea and the sea surface to assist the evaluation of the possible changes in the vertical wind profile by the atmospheric stability.

Conclusions
This study analyzed the wind speed data of the met mast located at the Formosa I offshore wind farm from May 2017 to April 2018 to understand the wind characteristics of the Taiwan Strait: • The annual mean wind speed at 90-m height is approximately 8 m/s, but the standard deviation is >5.5 m/s. The wind speed has a somewhat large spreading. • Two prevailing wind directions exist, which are almost opposite, namely the north-east monsoon in winter and the south-west monsoon in summer. The intensity of the north-east monsoon is usually greater than that of the southwest monsoon.

•
Winter afternoon has the highest wind speed of the year. Conversely, summer night and early morning have the lowest wind speed of the year.

•
The exponent of the vertical profile of the horizontal wind is low at approximately 0.05.
Compared with Weibull distribution, the histogram has a higher distribution probability for wind speeds < 4 and > 11 m/s and a lower distribution probability for wind speeds from 4 to 11 m/s. It appears to be a bimodal distribution. Therefore, a mixture distribution consisting of two distributions is used to fit the histogram. The mixture gamma distribution has the best performance for depicting wind data, which is confirmed by Sections 4.1 and 4.2. In addition to the diurnal variation of wind speed with associated changes in wind direction (Figure 6), the met mast is only approximately 5 km away from the coastline. Thus, the wind field near this wind farm is greatly affected by the monsoon and land-sea breeze. Land-sea breeze affects not only the wind speed but also the vertical gradient of the wind.
On the basis of the aforementioned findings, we can conclude two more wind characteristics of the Taiwan Strait.

•
The bimodal probability distribution function depicts wind accurately. A correct probability distribution function is crucial. Only with the correct probability distribution function can we correctly assess the wind energy potential and select an economically feasible site. Accordingly, after selecting the wind turbine, the AEP can be calculated correctly.

•
The strength and direction of the land-sea breeze have significant effects on Taiwan's peak electricity consumption hours in the summer afternoons.
To improve the accuracy of wind energy forecasting, two factors should be considered: the temperature difference between the sea surface and the ground and the temperature difference between the air above the sea and the sea surface. Indubitably, to better understand the effect of land-sea breeze on wind power production, it is highly recommended to conduct an extensive land-sea breeze observation experiment and to integrate high-resolution numerical weather prediction models so as to thoroughly understand the strength and extension of land-sea breeze.
In the future, we hope to obtain high sampling frequency wind data to evaluate the effect of wind speed extremes on the structural safety of wind turbine blades. Additionally, because the vertical gradient of wind speed and the wave height will mutually interact with each other, the extreme value will help to estimate the wave height accurately [44], thus providing an accurate vertical gradient of wind speed. For the sustainability operation of wind farms, long-term observation data collected together with the model output (such as reanalysis data) will help to evaluate the wind climate transition, which is certainly a worthwhile research issue [19,[45][46][47].