Adjusted spectral correction method for calculating extreme winds in tropicalcycloneaffected water areas
 Department of Wind and Energy Systems, Technical University of Denmark, Risø Campus, Roskilde, Denmark
 Department of Wind and Energy Systems, Technical University of Denmark, Risø Campus, Roskilde, Denmark
Correspondence: Xiaoli Guo Larsén (xgal@dtu.dk)
Hide author detailsCorrespondence: Xiaoli Guo Larsén (xgal@dtu.dk)
A method is developed to calculate the extreme wind for tropicalcycloneaffected water areas. The method is based on the spectral correction method by Larsén et al. (2012) in connection with the use of numerically modeled data, where an enhancement coefficient is derived as a function of wind speed to reflect the large wind fluctuation during tropical cyclones. This is done through calibration with the estimates from Ott (2005), who used the best track data and Holland model to estimate the extreme wind over the typhoonaffected area in the western North Pacific. The method is applied in the current study to three regions where the 50year winds with an effective temporal resolution of 10 min are obtained at 10, 50, 100 and 150 m. The results are in agreement with Ott (2005) over their study domain, though with much more spatial details of the extreme wind distribution.
The 50year wind at turbine hub height (U_{50}) is one of the most important siting parameters that needs to be estimated for the wind turbine design (e.g., IEC, 2019; Yu et al., 2011). Finding the correct type of wind turbines for a particular area is a necessary step for regional wind energy planning, as it is related to levelized cost of energy (LCoE), which involves wind resource, safety, and risk in connection with design, operation, and maintenance.
Wind energy development has been expeditious in the past decade, and wind turbines are being raised or planned in many regions globally, including tropicalcyclone affected offshore areas.
The most reliable method to estimate U_{50} should be to use goodquality measurements at the site that are sufficiently long, e.g., more than 10 years. This is however a condition that is almost never satisfied in these areas. Measurements under hurricane conditions are difficult and expensive to obtain, due mostly to technical difficulties at winds of such strength. A few meteorological stations exist in some tropicalcycloneaffected regions over coastal land areas. It is far from being enough for assessing the extreme conditions, particularly over water. Moreover, most of these measurements are not accessible for general research purposes.
Alternatively, people use modeling approaches to obtain extreme wind statistics related to tropical cyclones. The oftenused modeling approaches include the stochastic Monte Carlo simulation (MCS) (e.g., IEC, 2019), the simplified physical Holland hurricane model (Holland, 1980), and numerical weather models.
In this study we focus on physicsbased modeling approaches and therefore will not discuss MCS. In using the Holland model, Ott (2005) jointly used the best track data from 1977 to 2004 and calculated U_{50} at 10 m for the region of western North Pacific ocean; his method is referred to as the Ott method here. The values of U_{50} from Ott (2005) are of an equivalent temporal resolution of 10 min. The simple Holland model describes some of the most important parameters of a tropical cyclone, including the minimum center pressure, maximum wind speed, and the distance to the cyclone eye. It also assumes a spatially symmetric distribution of wind speed from the cyclone eye to outer region. The simplicity of the Ott method at the same time caused its limitation, e.g., the lack of detailed spatial information of the wind speed, horizontally and vertically. Moreover, the Ott method needs feedin information of the best track data.
Numerical modeling has been a powerful tool in simulating tropical cyclones, as being done in research and forecasting centers for hurricanes and typhoons. These weather models can be run at a spatial resolution from a few kilometers to a few tens of kilometers, with outputs mostoften saved hourly. These longterm simulation data are valuable for assessing the extreme wind in tropicalcycloneaffected areas. However, most of these simulations are forecasts with model schemes and setups that undergo regular and continuous updates and improvements. The updates are advantageous for forecasting purposes but not necessarily for calculating the extreme winds, because inconsistency is introduced into the longterm data, which will affect the extreme wind samples artificially. Moreover, the data from weather centers are not always accessible to wind energy users.
The reanalysis data thus have been a very attractive option, as they are generated using a consistent setup and they are globally available with open access. Some of them, e.g., Climate Forecast System Reanalysis (CFSR) and the fifthgeneration ECMWF reanalysis (ERA5), come with a spatial resolution of about 25–40 km, available on an hourly basis. For instance, Pryor and Bartelmie (2021) have used the ERA5 data and created a global atlas of extreme winds at 100 m. However, their validation suggests that the extreme winds in tropicalcycloneaffected areas are significantly underestimated. This is nevertheless expected, as Larsén et al. (2012) showed through spectral analysis of the modeled wind speed time series. Wind time series from mesoscale numerical models in general suffers from the numerical smoothing effect, introduced by a coarse grid that facilitates the convergence of the model. Larsén et al. (2012) showed that the modeled time series misses the highfrequency variability in comparison with measurements of similar resolution, causing systematic underestimation of the 50year wind. This is the case almost everpresent in data from mesoscale modeling, and the smoothing effect is expected to be significant for reanalysis data of tens of kilometers spatial resolution. To solve this problem, Larsén et al. (2012) developed a socalled spectral correction (SC) method to fill in the missing wind variability from modeled time series through a spectral model. Thus the corrected time series will follow the power spectrum down to the temporal resolution of the measurements. To meet the IEC standards, we can correct the spectrum to an equivalent resolution of 10 min. This SC method is briefly introduced in Sect. 2.2.1 here.
The SC method has been used in connection with different reanalysis data to predict extreme winds globally, including CFSR and Climate FourDimensional Data Assimilation (CFDDA) (Larsén and Kruger, 2014; Hansen et al., 2016; Larsén et al., 2022). The SC method has been shown to be reliable when validated with midlatitude measurements. However, applying SC in the same way as for midlatitude storms suggests a significant underestimation of U_{50} for the western North Pacific Ocean, when comparing with the results from Ott (2005).
The current study provides a simple approach for all tropicalcycloneaffected water areas, combining the strength of the method from Ott (2005) and Larsén et al. (2012). This method is denoted here as the spectral correction for tropical cyclone method, in short, the SCTC method.
This new method is described in Sect. 2. Section 3 shows results with discussions for the western Pacific Ocean as well as two other regions that are under the impact of tropical cyclones. Summary and conclusions are provided in Sect. 4.
2.1 The data
In this study, in connection with the development of the SCTC method, we have used the wind speed at 10 m from CFSR1 reanalysis data. The CFSR1 data are available from 1979–2010, and they are hourly values at a spatial resolution of about 40 km (Saha et al., 2010). There are version2 CFSR data, CFSv2, which are available from 2011, with a higher spatial resolution of about 25 km. However, the total data length of CFSv2 is too short for the calculation of the 50year wind. The CFSv2 data are, however, used to investigate the wind speed spectra during the case of Typhoon Megi in Sect. 2.2.2. Outputs for Megi from the mesoscale Weather Research and Forecasting (WRF) model, with a spatial resolution of 2 km, are also used.
When developing the SCTC method, we also used the estimation of the 50year wind at 10 m from Ott (2005), which is derived from best track data, and Holland model, for an area over the western North Pacific Ocean. These values, called U_{50,BT} here, where BT stands for best track, are available for a grid size of about 1^{∘}. They are from Fig. 13 in Ott (2005) and reproduced here as contour lines in Fig. 4a.
2.2 The spectral correction method for tropical cyclone conditions
2.2.1 The spectral correction method
As explained in Sect. 1, the smoothened time series from numerical models results in missing wind variability and, accordingly, translates into low spectral energy at higher frequencies when compared to measurements (Larsén et al., 2012). Figure 1 shows an example of the spectrum from modeled time series (black curve, from the 32year CFSR1 data) in comparison with expected spectral tail that has a slope of $\mathrm{5}/\mathrm{3}$ in a log–log coordination (solid red curve).
The effect of the smoothing of the time series on the estimation of the extreme wind is calculated in Larsén et al. (2012) by assuming the time series follows a Gaussian process. The onceperyear exceedance can thus be described using the Poisson distribution. Following this, in Larsén et al. (2012), the wind that occurs once a year ${\overline{U}}_{\mathrm{max}}$ was derived as a function of the zeroth and secondorder spectral moments m_{0} and m_{2}:
where $\overline{U}$ is the mean wind speed, T_{0} is the basis period of 1 year, and m_{i} is the ith spectral moment defined by
where f is the frequency in hertz and S(f) is the spectrum of the wind speed. From the above equations, it is clear that ${\overline{U}}_{\mathrm{max}}$ is significantly affected by the highfrequency part of the spectrum through m_{2}. Thus, we can calculate the smoothing effect on ${\overline{U}}_{\mathrm{max}}$ from the modeled time series by correcting its spectral tail.
Such a correction has been done by replacing the highfrequency part of the spectrum from the modeled time series with the following spectral model:
from a certain frequency f_{c} to the expected frequency f_{h}. In our study it is the 10 min temporal resolution we aim at so that f_{h}=72 d^{−1} (0.00083 Hz), which is the Nyquist frequency for 10 min. Equation (3) is the mesoscale part of the expression from Larsén et al. (2013), which includes both a synoptic and mesoscale range of spectrum: $S\left(f\right)=a\cdot {f}^{\mathrm{5}/\mathrm{3}}+b\cdot {f}^{\mathrm{3}}$. As shown by Fig. 1, a is determined by the time series of the model data where the black curve and the red solid curve meet. In Hansen et al. (2016), f_{c}=0.8 d^{−1} was used. In Larsén et al. (2022), f_{c} and S(f_{c}) were chosen from a regression line of ln S(f) with ln f for the range $\mathrm{0.6}<f<\mathrm{0.9}$ d^{−1}; thus the choice of f_{c} becomes less sensitive to the fluctuation in S(f).
Using the original spectrum and the corrected spectrum, together with the mean wind speed $\overline{U}$ to Eq. (3), we obtain ${\overline{U}}_{\mathrm{max},\mathrm{ori}}$ and ${\overline{U}}_{\mathrm{max},\mathrm{SC}}$, and their ratio $R={\overline{U}}_{\mathrm{max},\mathrm{SC}}/{\overline{U}}_{\mathrm{max},\mathrm{ori}}$. We use this ratio R to correct the annual maximum wind speed obtained from the time series directly. From a time series of Nyear data, we obtain N samples of the annual maximum wind speed. After applying the spectral correction, we use the annual maximum method with the Gumbel distribution to obtain the 50year wind (e.g., Abild et al., 1992; Larsén et al., 2019).
2.2.2 Example of the spectral behaviors during a tropical cyclone
The abovementioned SC method was shown to provide reasonable estimates for midlatitude storms (e.g., Hansen et al., 2016; Larsén and Kruger, 2014). However, using Eq. (3) results in underestimation of U_{50} for the tropicalcycloneaffected areas. This is because the energy level in the mesoscale range during tropical cyclones is significantly higher. This is demonstrated in Fig. 2 through an example of spectral analysis during the case of Typhoon Megi that passed Taiwan during the period 26–27 September 2016. It is not the focus of current study to investigate mesoscale modeling of the tropical cyclones; thus, details of the simulation of Typhoon Megi using the WRF model are given in Appendix A. As the period is too short for a Fourier analysis of the time series in time domain, we examine the power spectra in wavenumber domain, namely, E(k) vs. wavenumber (k). We compare the spectral behavior and energy level for the wind speed between this particular case and climatological conditions in Fig. 2. Spectra E(k) corresponding to Megi were calculated over an area of about 200 by 200 km over the Pacific Ocean as Megi approached Taiwan during the period 26 September 2016 at 10:00 UTC to 27 September 2016 at 09:00 UTC; both the reanalysis data CFSv2 (blue in Fig. 2) and the WRF data were used (red). We calculated and averaged the onedimensional spectra in the north–south (N–S) direction as well as west–east (W–E) direction. They are plotted as solid and dashed curves, respectively, in Fig. 2. Thus, the spectra represent an average energy level over the entire area, which is more representative for this particular case in comparison with a spectrum derived with time series at a certain grid point. Here the climatological condition is represented by two data sources: one with the CFSv2 data of the entire year 2016 (green curves) over the same domain and the other as the wavenumber spectrum from Gage and Nastrom (1986), the Gage–Näström spectrum (the black curve). The Gage–Näström spectrum was obtained from thousands of commercial airplane flight measurements, and it is often used to represent the climatological power spectrum in the troposphere. The Gage–Näström spectrum has also been verified by the theoretical work of Lindborg (1999) on general twodimensional turbulence behavior. For the climatological condition, the energy level and spectral slope are similar for the overlapping wave number range between the 1year CFSR1 data and the Gage–Näström spectrum, and both have significantly lower energy level than that for the Megi case. For this particular case, the higher resolution of the WRF model simulation, at 2 km, produced a spatial wind variability that is 3–5 times larger than that of the CFSv2 data at about 25 km resolution.
2.2.3 The SCTC method
To adjust the SC method for more general storm conditions, we revise Eq. (3) to
with n a coefficient reflecting weather types, and its value can be adjusted with additional data, e.g., measurements. Thus for midlatitude storms, n=1, and for tropical cyclones, we develop a systematic way to determine n, which is described in the following. Here we need to decide how to link the new spectral tail (red dashed curve in Fig. 1) to the longterm spectrum (the black curve). To keep it simple and to be on a relatively conservative side, we use the blue line to proceed with the calculation. When we are using Eq. (4), it is equivalent to generate a climatology with tropical cyclone everpresent during a year. Since we are only interested in the smoothing effect of the modeled time series on the onceayear event, with rest of the data not considered, such an assumption is acceptable. Here, to match the IEC requirement, as also applied in Ott (2005), we correct the values to the equivalent 10 min resolution; hence f_{h}=72 d^{−1} is used when calculating the spectral moments using Eq. (3).
With an input of n, the SCTC method would provide an output of the ratio of the corrected (${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$) and uncorrected (U_{50,uncorr}) value of the 50year wind, r. By varying n in Eq. (4) in connection with the use of the SC method, different ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$, and hence r, are obtained. Thus a relation between n and r can be derived. Such a relation depends on the initial wind speed time series, and it therefore depends on the location of the interest in relation to the tropical cyclone structure, with n increasing faster with r for stronger wind region. This is demonstrated through the calculation from five grid points which are marked in Fig. 4: P1 (17.5^{∘} N, 131^{∘} E); P2 (24^{∘} N, 130^{∘} E); P3 (26.5^{∘} N, 131^{∘} E); P4 (29^{∘} N, 133^{∘} E) and P5 (32.5^{∘} N, 136^{∘} E). The corresponding values of U_{50,BT} at the five grid points are 73, 67, 58, 53 and 45 m s^{−1}, respectively. P1 corresponds to the largest value over the domain. Here we used n from 1 to 30 to obtain 30 corrected ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ and hence r, and we fit the data with a secondorder polynomial function. The corresponding n−r relation is obtained as
For P1, α=28.28, $\mathit{\beta}=\mathrm{30.24}$ and $c=\mathrm{0.66}$. The 30 points and Eq. (5) are shown in Fig. 3 for P1. In the same figure, similar data for the other four grid points P2 to P5 are also plotted. It is clear that the n−r relation is not uniform over the space. To be on the conservative side, we use the n−r relation from P1. The uncertainty related to this choice will be discussed in Sect. 3.
In finding n for Eq. (4), here we derive r as a function of a wind parameter for which the data are globally available freely. Firstly, 32year annual maximum winds u_{max} are extracted from the CFSR1 data from 1979 to 2010. The 50year wind at its original resolution, U_{50,uncorr}, is calculated applying the Gumbel distribution to the 32 values with the annual maximum method. Secondly, we use the results of the 50year wind from Ott (2005), U_{50,BT} (see Sect. 2.1), to train the CFSR1 data. We regrid U_{50,uncorr} (spatial resolution of about 40 km) to the grid points corresponding to U_{50,BT} (spatial resolution of 1^{∘}) and calculate the ratio ${U}_{\mathrm{50},\mathrm{BT}}/{U}_{\mathrm{50},\mathrm{uncorr}}$. For the regridding, we simply find the value at the closest grid point without any interpolation. The values of the ratio cover a range from 1 to 2.6, as shown in Fig. 5a. Based on the scatter plot of the ratio with U_{50,BT} in Fig. 5a, we can describe r in a linear relation with U_{50,BT} for different wind speed ranges. A linear regression is made for $\mathrm{27.5}<{U}_{\mathrm{50},\mathrm{BT}}<\mathrm{65}$ m s^{−1} following
where a_{1}=0.0163 and b_{1}=0.569. This relation is shown in Fig. 5a as a red dashed line.
We use a similar dependence of r with uncorrected model value u to Eq. (6):
for $\mathrm{27.5}<u<\mathrm{60}$ m s^{−1}; u is the uncorrected 50year wind. Here we use U_{50,BT} to calibrate Eq. (7). We use the same slope as Eq. (6), a_{2}=a_{1}. By requiring the calculation of the 50year wind through Eq. (4), U_{50,corr}, to match U_{50,BT} at P1, we loop the calculation and obtain b_{2}=0.62. One constant of r, ${a}_{\mathrm{2}}\cdot \mathrm{27.5}+{b}_{\mathrm{2}}$, is then obtained for ${U}_{\mathrm{50},\mathrm{BT}}<\mathrm{27.5}$ m s^{−1} and another one, ${a}_{\mathrm{2}}\cdot \mathrm{60}+{b}_{\mathrm{2}}$, for an upper wind speed limit, e.g., u>60 m s^{−1}. Thus,
We stopped the loop at n=15, when ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}=\mathrm{73.24}$ m s^{−1}. Equation (8) is plotted in Fig. 5b together with Eq. (6). We then use r calculated from Eq. (8) as input to Eq. (5), to obtain n. For the initial, uncorrected 50year wind that is less than 27.5 m s^{−1}, we use n=1, which is the same as Eq. (3). Thus, the SCTC method is ready to be used and it coincides with the SC method for u<27.5 m s^{−1}.
2.3 Extreme wind at different heights
Usually, model wind outputs are provided at a few certain heights, e.g., 10 or 100 m, and most models provide outputs at 10 m. The calculation of U_{50} in Ott (2005) is at height z=10 m, and the CFSR1 data we used so far are also at 10 m. Modern offshore wind turbines can be as tall as 200 m. Therefore, we introduce a simple approach here to obtain winds at other heights U_{z} from 10 m winds for storm conditions.
In extreme wind conditions, the logarithmic wind law can be considered valid up to several hundreds of meters according to Sonde wind speed measurements during hurricane conditions, as shown in, e.g., Powell et al. (2003) and Giammanco et al. (2013). Thus we can calculate U_{z} using
where u_{*} is the frictional velocity and κ=0.4 is the von Kármán constant. With U_{10 m} at z=10 m known, U_{z} can thus be obtained analytically if we have one more relation for z_{0} and u_{*}.
There are several such relations established in the literature for describing z_{0} and u_{*}. Here we examined three of them, with the first one the Charnock formulation for welldeveloped wind sea waves (Charnock, 1955):
where α_{ch} is the Charnock coefficient, which depends on the sea state, and it is larger for rougher sea. In numerical models, it is often between 0.011 and 0.02, as a function of wind speed. Here for the strong winds, we use α_{ch}=0.02.
The second one is the Andreas algorithm that takes into account of spray effect at strong wind conditions (Andreas et al., 2015):
And the third are the algorithms that are used in wave model SWAN (e.g., Zijlema and van der Westhuysen, 2005). The SWAN algorithms are based on field measurements over a wide range of wind speed conditions, including hurricane conditions. It suggests a leveloff or decreasing of surface drag coefficient (C_{d}) or roughness length with wind speed at strong winds, interpreted as a result from wavebreaking processes:
In a test of calculating the wind speed at 100 m using wind speed at 10 m, it is found that the difference in using the three abovementioned formulations is negligible for winds at 100 m up to 25 m s^{−1}. While the Charnock formulation with α_{ch}=0.02 gives very similar results to the Andreas formulation up to 40 m s^{−1}, it provides larger values for stronger winds, with a difference of about 1 at 50 m s^{−1}. Both the Charnock and the Andreas formulations give larger estimations than the SWAN formulations for winds larger than 25 m s^{−1}, and the overestimation increases to about 3 at 50 m s^{−1}. In this study we proceed with our calculation using the SWAN formulation as it provides a surface wavebreaking effect at strong winds, which is calibrated with measurements (though very few) during hurricane conditions. We note that all of these formulations are an oversimplification of the sea surface at such strong wind conditions during tropical cyclones and uncertainty is embedded.
The SCTC method is applied to three areas in this study and they are marked in Fig. 6 in boxes, with box1 covering the same region as in Ott (2005) for typhoons and boxes 2 and 3 for hurricanes. These three areas are chosen as they host the most severe and frequent tropical cyclones, based on the International Best Track Archive from National Centers for Environmental Information (NOAA) (NOAA, 2017).
For the three areas, the 50year winds are calculated at 10, 50, 100 and 150 m. With values at the four heights available, one can also obtain values at any height in between using an interpolation or extrapolation method, e.g., polynomial. The data of these 50year winds at the four heights are provided in the database (Larsén and Ott, 2022).
In order to show the effect of the spectral correction in comparison with the result from Ott (2005) (reproduced here as Fig. 4a), in Fig. 4, we plotted the 50year wind at 10 m directly from the reanalysis data without spectral correction, U_{50,uncorr} (Fig. 4b), the 50year wind at 10 m using the spectral correction through Eq. (3) (Fig. 4c) and the 50year wind at 10 m using the spectral correction for tropical cyclones through Eq. (4) (Fig. 4d). Figure 4b, c and d share consistent wind distribution patterns. Compared to U_{50,BT}, the magnitude of the 50year wind U_{50,uncorr} is significantly underestimated in Fig. 4b (with the maximum value 50.8 m s^{−1}). The spectral correction method with Eq. (3) brings the maximum value of 50.8 to 58.8 m s^{−1} (Fig. 4c), which is still significantly smaller than the corresponding U_{50,BT}: 73.02 m s^{−1}. Comparing Fig. 4a and d, U_{50,BT} vs. ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$, one can see that, first of all, the wind speed range is comparable for the common area of the two plots, being from about 35 to 73 m s^{−1}. Secondly, the locations of the strongest winds are consistent in the two plots, being between 120–140^{∘} E and 15–25^{∘} N. The contour lines in Fig. 4a are very smooth due to the use of the Holland model and coarseresolution best track data (1^{∘}) in the Ott method, while Fig. 4b reveals richer spatial details.
When calculating ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$, the following three new equations have been used: Eqs. (4), (5) and (8). Among them, the coefficients in Eq. (5) are derived with input of the model time series in connection of the use of the spectral correction method; we have used the time series from P1 for this. The coefficients in Eq. (8) have been calibrated using U_{50,BT} so that ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ matches U_{50,BT} at P1. For both Eqs. (5) and (8), we have used data from P1. By doing so, we are on the conservative side of the calculation and at the same time have saved huge computational cost. Ideally such calibration should be done to all CFSR grid points. However, this can not be done due to the dismatch of grid size between U_{50,BT} and the reanalysis data. Moreover, we do not have U_{50,BT} outside the domain of Fig. 4a. This implies that it is inevitable that we will introduce assumptions for using Eqs. (5) and (8) for most of the reanalysis grid points. Using the calibration from one point, P1, thus ensures a consistent workflow.
It is nevertheless relevant to address the uncertainty related to the use of calibration at one point. In theory, if we calibrate at all grid points using the same procedure as for P1, we get similar values of U_{50,BT} and ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ at these grid points. Thus, the difference between Fig. 4a and d will be able to directly suggest the related uncertainty. In Sect. 2, we introduced five grid points (Fig. 4), whose U_{50,BT} varies from about 45 to 73 m s^{−1}; the exact numbers are listed in Table 1. Among the five locations, the agreement between ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ and U_{50,BT} is good for P1, P3 and P4, while at P2 ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ is still underestimated, and at P5 ${U}_{\mathrm{50},\mathrm{SC}\mathrm{TC}}$ is overestimated. We need to remember that the spatial distribution of U_{50,BT} is of very coarse resolution, which can miss spatial details and result in large uncertainties in some individual locations; P2 could be such a case. For the five points, we also calculated U_{50} using Eq. (8) with the same coefficients as for P1, but together with their respective n−r relations from the five points (Fig. 3). The values are also shown in Table 1. Here one can see that the estimates of the 50year wind are underestimated.
We point out that in deriving the spectral correction method (Sect. 2.2.1), we have assumed Gaussian process for the wind speed. As also pointed out and discussed in the original paper Larsén et al. (2012) that a Gaussian process is not ideal for describing the wind speed; it however makes it possible for us to derive analytically the effect on the onceinayear event in a time series that misses highfrequency fluctuations. This also adds to the overall uncertainty.
The agreement between Fig. 4a and d is encouragement to apply the SCTC method to other tropicalcycloneaffected areas.
Figure 7 shows 10 min U_{50} at 100 m, which is a more relevant height for modern offshore turbines, over the three areas as marked in Fig. 6. Note that the three boxes do not include all areas with the impact of tropical cyclones. Note also that the current study only addresses “water areas”; this is mainly for two reasons. Firstly, extrapolation of surface winds to other heights over land requires much more complicated modeling approaches due to spatially inhomogeneous surface conditions. Moreover, it is not a trivial task to obtain such data for the surface conditions, whereas over water we benefit from the dependence of roughness length on the wind speed, even though there is uncertainty introduced. In addition, the best track data used for calibration in connection with the SCTC method are mostly available over water.
When we extrapolate the wind speed from 10 m to higher altitudes, we have used the roughness length algorithms that are used in wave model SWAN, Eq. (12). These algorithms suggest that at hurricane strength, water surface roughness becomes smoother as a result of wave breaking and foaming processes. The corresponding roughness length is thus smaller than what the Charnock formulation Eq. (10) gives, and thus it will provide smaller wind speeds at higher altitudes. One may argue whether we should use a more conservative approach again such as the Charnock formulation. We choose not to use it for the following reasons. Firstly, it has been proven incorrect at very strong winds. Secondly, Eq. (12) has been supported by measurements, even though very few. Thirdly, the difference caused by Eqs. (10) and (12) is systematic; one can easily include this uncertainty in the assessment once there are further measurements for validation.
Even though we have used the CFSR1 data in this paper, the SCTC method can be applied similarly with other reanalysis data with sufficient data length. However, the quality of the reanalysis data, including the information of tropical cyclone parameters such as path, intensity and spatial structures, needs to be quality checked first. Imberger and Larsén (2022) show that the characteristics of tropical cyclones are quite different in the following reanalysis data: CFSR1, MERRA2 (ModernEra Retrospective analysis for Research and Applications), ERA5 and CFDDA. For instance, among the four, the patterns of tropical cyclones are very weak in the CFDDA data. This suggests an important source of uncertainty associated with the input reanalysis data.
It should be pointed out that using the dependence of n on the only one parameter, the wind speed, allows the use of the SCTC method in areas of strong winds that are not necessarily affected by tropical cyclones. It can be improved with further input of information when possible, to indicate whether it is an area with tropical cyclone impact or not. In addition, the use of Eq. (4) has now been calibrated for midlatitude storms and tropical cyclones; it has not yet been examined for other types of extreme wind events such as thunderstorms. Further studies are needed to extend the application of the method.
This study develops a method for calculating the extreme wind for tropicalcycloneaffected areas, here called the SCTC method. This is done by adjusting the spectral correction (SC) method from Larsén et al. (2012), through adding an enhancement coefficient to the spectral model, which is a function of the local extreme wind. Such a dependence is calibrated using the estimates from Ott (2005), who used the best track data and Holland model for the area of southern North Pacific Ocean.
Here we summarize the main steps of the recipe for this method, which can be used in connection with other data sources, for instance, a new atlas that is based on measurements, and/or new reanalysis data:

Calculate U_{50} from the original model time series over the domain: U_{50,uncorr}.

Find out the maximum value of U_{50} from U_{50,BT} and the corresponding location (p_{0a}). We call this maximum value ${U}_{\mathrm{50},\mathrm{BT},\mathrm{max}}$.

Find out the maximum value of U_{50} from U_{50,uncorr} and the corresponding location (p_{0b}). In our study, p_{0a} and p_{0b} are very close to each other and we used p_{0b}, which is P1.

Extract model time series from p_{0b} and apply the spectral correction method with n ranging from, e.g., 1 to 30 to obtain the corrected values (we call it U_{50,test} here). Derive a relation between n and ${U}_{\mathrm{50},\mathrm{test}}/{U}_{\mathrm{50},\mathrm{uncorr}}$; thus Eq. (5) can be defined.

Derive Eq. (8) from the scatter plot of ${U}_{\mathrm{50},\mathrm{BT}}/{U}_{\mathrm{50},\mathrm{uncorr}}$ vs. U_{50,BT}, and find the coefficients through looping the calculations until U_{50,test} becomes close enough to ${U}_{\mathrm{50},\mathrm{BT},\mathrm{max}}$.
The three equations Eqs. (4), (5) and (8) are thus established. With an input u as the uncorrected 50year wind at 10 m from a given location, r can be calculated using Eq. (8). With r as input to Eq. (5), we obtain n, which can now be used in the SCTC method using Eq. (4).
The results from the SCTC method provide consistent wind distribution of 10 min 50year wind at 10 m, and in addition much richer and more realistic spatial details, in comparison with the results from Ott (2005). The method can be applied for all water areas with tropical cyclones and can be used to obtain extreme winds from the surface to a few hundreds of meters.
We acknowledge that extension of the SC method to the SCTC method is quite empirical, and there might be other approaches for improving the calculation. In spite of the uncertainties, the method and the data produced in this study serve to fill in the void of extreme wind estimations in most tropicalcycloneaffected areas. We acknowledge that we do not have enough measurements for validation, other than the study of Ott (2005), and therefore encourage validation from measurement owners.
There are two purposes for using the mesoscale Weather Research and Forecasting (WRF) model. The first is to examine the spectral energy level during an example case of typhoon in comparison of climatological conditions. The second is to compare the spectral energy level from the CFSv2 data (with a spatial resolution of tens of kilometers) with that from the WRF data (with a spatial resolution of 2 km).
We used WRF version 4.0, configured with the moving nest function. Three nested domains are used, with spatial resolutions of 18, 6 and 2 km, respectively. The innermost domain has 339 by 342 grid points. We used 52 vertical model levels from the surface to a pressure level of 5000 Pa. We used the new Thompson microphysics scheme (Thompson et al., 2004), the RRTMG scheme for long and shortwave radiation physics scheme (Iacono et al., 2008), MYNN 3.0 PBL scheme (Nakanishi and Niino, 2009) and Noah Land Surface Model. The Kain–Fritsch cumulus scheme (Kain and Fritsch, 1993) was used for the outer domain (I) but not for domains II and III. We used ERA5 data as initial and boundary conditions for WRF. The daily 0.25^{∘} OISST data were used to define the sea surface temperature conditions. The simulation started at 26 September 2016 at 00:00 UTC and ended at 27 September 2016 at 12:00 UTC, with the first 12 h as a spinup period. The model outputs are recorded every 10 min.
Figure A1 shows the typhoon track from three data sources, the best track data, the CFSv2 data as the position of the lowest mean sea level pressure (MSLP) and the WRF data, also as the position of the lowest MSLP. They are in good agreement.
Both the wind speed data from CFRv2 and WRF from 12:00 UTC on 26 September to 12:00 UTC on 27 September over the innermost model domain are used to obtain the mean spectra of the longitudinal and meridional wind as shown in Fig. 2.
The CFSR data are publicly available at https://rda.ucar.edu/datasets/ds093.1 (NCAR/UCAR, 2022). The 50year wind values over the three areas at 10, 50, 100 and 150 m are available at https://doi.org/10.5281/zenodo.7089426 (Larsén and Ott, 2022).
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This study is supported by the EUDP projects GASP (EUDP J. no. 640180095) and GASPOC (EUDP J. no. 640201043). The CFSR products are downloaded from https://rda.ucar.edu/datasets/ds093.1 (last access: 28 February 2022).
This research has been supported by the EUDP project GASP (EUDP J. no. 640180095) and the GASPOC project (EUDP J. grant no. 640201043).
This paper was edited by Julia Gottschall and reviewed by Javier Sanz Rodrigo and two anonymous referees.
Abild, J., Mortensen, N., and Landberg, L.: Application of the wind atlas method to extreme wind speed data, J. Wind Eng. Ind. Aerodyn., 41, 473–484, 1992. a
Andreas, E. L., Mahrt, L., and Vickers, D.: An improved bulk airsea surface flux algorithm, including spraymediated transfer, Q. J. Roy. Meteorol. Soc., 141, 642–654, https://doi.org/10.1002/qj.2424, 2015. a
Charnock, H.: Wind stress on a water surface, Q. J. Roy. Meteorol. Soc., 81, 639–640, 1955. a
Gage, K. and Nastrom, G.: Theoretical interpretation of atmospheric wavenumber spectra of wind and temperature observed by commercial aircraft during GASP, J. Atmos. Sci., 43, 729–740, 1986. a, b
Giammanco, I. M., Schroeder, J. L., and Powell, M. D.: GPS Dropwindsonde and WSR88D Observations of Tropical Cyclone Vertical Wind Profiles and Their Characteristics, Weather Forecast., 28, 77–99, https://doi.org/10.1175/WAFD1100155.1, 2013. a
Hansen, B. O., Larsén, X. G., Kelly, M., Rathmann, O. S., Berg, J., Bechmann, A., Sempreviva, A. M., and Jørgensen, H. E.: Extreme wind calculation applying spectral correction method – test and validation, Tech. Rep., Wind Enery Department, Technical University of Denmark, DTU Wind Energy E0098, 2016. a, b, c
Holland, G.: An analytic model of the wind and pressure profiles in hurricanes, Mon. Weather Rev., 1212–1218 pp., 1980. a
Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., and Collins, W. D.: Radiative forcing by longlived greenhouse gases: Calculations with the AER radiative transfer models, J. Geophys. Res., 113, 13013, https://doi.org/10.1029/2008JD009944, 2008. a
IEC: IEC 614001 Ed4: Wind turbines – Part 1: Design requirements, standard, International Electrotechnical Commission, Geneva, Switzerland, ISBN 2831810252, 2019. a, b
Imberger, M. and Larsén, X. G.: Sensitivity and quality assessment of global 50year return winds using reanalysis products and measurements, WindEurope, https://windeurope.org/annual2022/conference/posters/PO254/ (last access: 4 December 2022), 2022. a
Kain, J. S. and Fritsch, J. M.: Convective parameterization for mesoscale models: The KainFritsch scheme. The representation of cumulus convection in numerical models, Meteor. Monogr., Ameri. Meteor. Soc., 24, 165–170, 1993. a
Larsén, X. G. and Kruger, A.: Application of the spectral correction method to reanalysis data in South Africa, J. Wind Eng. Ind. Aerodyn., 133, 110–122, 2014. a, b
Larsén, X. G. and Ott, S.: Dataset to adjusted spectral correction method for calculating extreme winds in tropical cyclone affected water areas, Zenodo [data set], https://doi.org/10.5281/zenodo.7089426, 2022. a, b
Larsén, X. G., Ott, S., Badger, J., Hahmann, A. H., and Mann, J.: Recipes for correcting the impact of effective mesoscale resolution on the estimation of extreme winds, J. Appl. Meteorol. Climat., 51, 521–533, https://doi.org/10.1175/JAMCD11090.1, 2012. a, b, c, d, e, f, g, h, i, j
Larsén, X. G., Vincent, C. L., and Larsen, S.: Spectral structure of the mesoscale winds over the water, Q. J. Roy. Meteorol. Soc., 139, 685–700, https://doi.org/10.1002/qj.2003, 2013. a
Larsén, X. G., Du, J., Bolaños, R., Imberger, M., Kelly, M. C., Badger, M., and Larsen, S.: Estimation of offshore extreme wind from windwave coupled modeling, Wind Energy, 22, 1043–1057, https://doi.org/10.1002/we.2339, 2019. a
Larsén, X. G., Davis, N., Hannesdottir, A., Kelly, M., Svenningsen, L., Slot, R., Imberger, M., Olsen, B., and Floors, R.: The Global Atlas for Siting Parameters (GASP) project: extreme wind, turbulence and turbine classes, Wind Energy, 25, 1841–1859, https://doi.org/10.1002/we.2771, 2022. a, b
Lindborg, E.: Can the atmospheric kinetic energy spectrum be explained by twodimensional turbulence?, J. Fluid Mech., 388, 259–288, 1999. a
Nakanishi, M. and Niino, H.: Development of an improved turbulence closure model for the atmospheric boundary layer, J. Meteorol. Soc. Jpn, 87, 895–912, 2009. a
NCAR/UCAR: CFSR data, CFSv2 [data set], https://rda.ucar.edu/datasets/ds094.1, last access: 28 February 2022. a
NOAA: An Inventory of Tropical Cyclone Tracks, NOAA [data set], https://doi.org/10.25921/82ty9e16, 2017. a
Ott, S.: Extreme winds in the western North Pacific, Tech. Rep. RisoeR1544(EN), Risø National Laboratory, Roskilde, Denmark, https://orbit.dtu.dk/en/publications/extremewindsinthewesternnorthpacific (last access: 4 December 2022), 2005. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s
Powell, M. D., Vickery, P. J., and Reinhold, T. A.: Reduced drag coefficient for high wind speeds in tropical cyclones, Nature, 422, 279–283, 2003. a
Pryor, S. and Bartelmie, R. J.: A global assessment of extreme wind speeds for wind energy applications, Nature Energy, 6, 268–276, https://doi.org/10.1038/s41560020007737, 2021. a
Saha, S., Moorthi, S., Pan, H.L., Wu, X., Wang, J., Nadiga, S., Tripp, P., Kistler, R., Woollen, J., Behringer, D., Liu, H., Stokes, D., Grumbine, R., Gayno, G., Wang, J., Hou, Y.T., ya Chuang, H., Juang, H.M. H., Sela, J., Iredell, M., Treadon, R., Kleist, D., Delst, P. V., Keyser, D., Derber, J., Ek, M., Meng, J., Wei, H., Yang, R., Lord, S., van den Dool, H., Kumar, A., Wang, W., Long, C., Chelliah, M., Xue, Y., Huang, B., Schemm, J.K., Ebisuzaki, W., Lin, R., Xie, P., Chen, M., Zhou, S., Higgins, W., Zou, C.Z., Liu, Q., Chen, Y., Han, Y., Cucurull, L., Reynolds, R. W., Rutledge, G., and Goldberg, M.: The NCEP Climate Forecast System Reanalysis, Bull. Am. Meteorol. Soc., 91, 1015–1058, https://doi.org/10.1175/2010BAMS3001.1, 2010. a
Thompson, G., Rasmussen, R. M., and Manning, K.: Explicit forecasts of winter precipitation using an improved bulk microphysics scheme, PartI: Description and sensitivity analysis, Mon. Weather Rev., 132, 519–542, 2004. a
Yu, Q., Kim, K., and Lo, T.: Design Standards for Offshore Wind Farms, Tech. Rep., American Bureau of Shipping, contract M10PC00105, https://www.bsee.gov/sites/bsee.gov/files/taptechnicalassessmentprogram/670aa.pdf (last access: 4 December 2022), 2011. a
Zijlema, M. and van der Westhuysen, A. J.: On convergence behaviour and numerical accuracy in stationary SWAN simulations of nearshore wind wave spectra, Coast. Eng., 52, 237–256, https://doi.org/10.1016/j.coastaleng.2004.12.006, 2005. a