This study delivers a method and datasets for a global offshore atlas for turbulence intensity (TI) from 10 to 200 m. The method includes both surface-driven three-dimensional boundary-layer turbulence and large-scale two-dimensional turbulence. This systematically includes the effect of large-scale eddies, particularly at weak wind conditions, and hence significantly improves TI in weak to moderate wind conditions. This method describes water roughness length through a dependence on wave age and wind speed, which is suitable for moderate to strong wind conditions. The method also includes stability dependence through the Obukhov length. Based on theories and measurements in literature, algorithms for TI have been calibrated for heights up to 200 m. We use the ERA5 atmospheric and wave data to demonstrate the use of the method and create a global dataset. The results show satisfactory agreement with measurements and data from the literature.
HORIZON EUROPE Climate, Energy and Mobility101146689Energiteknologisk udviklings- og demonstrationsprogram65020-1043Introduction
Turbulence intensity (TI) is one of the most used parameters in wind energy. For example, it is needed to make decisions about the turbine class for a particular place, as required by the IEC-61400 standard . It is required to calculate loads and fatigue. It is also an input parameter for engineering wake modelling.
The most direct way to obtain TI is through wind speed measurements, from which the standard deviation and the mean are calculated during a given period T. For wind speed close to the ground, T is typically 10 min to 1 h. Cup anemometers, and lidar and sonic anemometers are often used to measure TI, with sonic anemometers usually providing the most reliable turbulence data due to their fast response to flow fluctuations, but the associated cost is also the highest.
Measurements are generally expensive to collect, particularly for offshore conditions due to costs related to e.g. operation, maintenance, and data transfer. Traditionally, met-mast or LIDAR are installed to measure the wind conditions for e.g. a year or more before the wind farm is being built. With the rapid and large-scale development of offshore wind energy, there is often limited time to take such measurements, and therefore the information on TI relies heavily on modelling.
When there is a lack of measurements, TI needs to be modelled. The Global Atlas for Siting Parameters (GASP) project provided a near-global modelled dataset of TI at a spatial grid spacing of 275 m from a height of 10 to 150 m. This dataset includes both land and water areas, with water areas reaching 200 km from coastlines, leaving most of open ocean excluded. The calculation of TI over water used the mean wind statistics from Global Wind Atlas data data as input to the Kaimal turbulence model , with the surface roughness length z0 described through the Charnock formulation. The dataset has been applied in the GASP project for defining turbine classes and thus only the moderate to strong wind ranges were used and validated.
This study aims to improve the methodology and datasets for TI for offshore conditions from the GASP project, covering all global offshore areas and addressing also low to moderate winds in addition to moderate to strong wind conditions, as in GASP.
The standard suggests a simple monotonic decrease of TI with wind speed hub height U:
TIIEC=Iref0.75+5.6U,
where 5.6 is in m s−1. Equation is obviously only valid over land, because over water it is expected that the surface becomes rougher as the winds increase, contributing to the higher intensity of turbulence in stronger winds, which is supported by measurements .
developed algorithms from boundary-layer surface scaling laws, where TI is a function of stability and boundary-layer height (zi), and depends on roughness length z0, surface momentum and heat fluxes. To take into account the wave effect, they used the Charnock formulation for the roughness length z0=αchu*2/g, where αch is the Charnock parameter, for which a typical value of 0.011 is taken. Thus, TI increases with wind speed over the water. The details of the equation for TI from are included in Appendix (Eq. ).
The standard recommends an increase of TI with wind speed and a decrease with height:
TIISO=0.061+0.43U10z10-0.22,
where U is the wind speed at height z.
Similar relationships between TI and U have been empirically derived through measurements. Before the ISO standard, used measurements from the Frøya site and proposed several similar expressions with U scaled with 10 m s−1 and z scaled with 10 m: (z/10)-0.2, primarily for the wind speed range of about 10 to 26 m s−1. The “linear model” from reads
TILin=0.0871+0.302U10-1z10-0.2.
The other two models, the “modified Vickery model” and the “drag coefficient model”, are included in Appendix (Eqs. and ).
The expressions from the ISO standard, from and from do not intend to include the decreasing dependence of TI on wind speed under lower wind speed conditions. added such a dependence at lower wind speed by imposing climatological unstable stratification, which corresponds to stronger vertical mixing and hence larger turbulence intensity.
This decreasing dependence of TI on U at low wind speeds was considered in when deriving an empirical relation based on measurements from the FINO masts:
TIC=c1U+c2+c3U.
The height dependence of TI is not explicitly included in the expression above but considered through the set of coefficients c1, c2 and c3. For the FINO cases, two sets of coefficients for each site were derived for 30 and 100 m, respectively.
The study of could be considered a summary for the above expressions, including the increase of TI with U for moderate to strong winds, a decrease for light winds, and a height dependence of TI similar to the ISO relationship and the Frøya expressions:
TIJeans=a1U(z)10+a2+a3U(z)10-1z10-0.22,
where a1=0.035, a2=0.0089 and a3=0.0402 are the default values. The coefficients a1, a2 and a3 can be tuned to represent site-specific data, and different values for the three coefficients are also provided for different stability conditions. This expression has been validated in with data up to 82 m.
In the above-mentioned TI expressions, only the scaling modelling approach from is based on physics arguments; the others are empirical, driven by fitting to measurement. This has caused a large variety of sets of coefficients to these expressions, depending on the sites and sometimes measurement heights.
With respect to the application of physics, high-fidelity models are sometimes used to calculate turbulence parameters from which TI can be calculated, such as the large-eddy simulation (LES). Due to its high computational cost, it can only be used for a very limited area, over a short period or in idealised settings. Mesoscale modelling outputs, although by design are not capable of resolving turbulence, are sometimes post-processed to match expected results e.g. in .
In this study, we create a global offshore atlas for TI using a new and cost-effective approach that is based on physical principles. As a first step, we model the turbulence from 1 h to 10 Hz. Usually, the calculation of turbulence and hence TI considers only surface-driven three-dimensional (3D) turbulence. The 3D turbulence is described through spectral models through, for instance, the Kaimal model , the Mann model and the various spectral models reviewed in . The dotted black curves in Fig. show examples of the wind speed power spectra driven by the 3D turbulence at the site Høvsøre at two heights (refer to Sect. for details of this figure). The Kaimal model requires input of surface information such as z0 or u*, and the Mann model requires turbulence spectra (either from measurements or from other models) as input to further calculate the three-dimensional flow. In most cases, these models are used under the assumption of neutral stability. In common, these models (represented by the black curves in Fig. ) do not include wind fluctuations from the larger-scales (blue curves for frequency f lower than approximately 10-3 Hz). With the current offshore wind energy development and deployment, wind farm clusters can be of the size of several hundreds of square kilometres, and the tallest wind turbine has a hub height of about 185 m, reaching a height of 310 m with rotors (26 MW turbine). At those scales, including the contribution of large-scale turbulence is necessary and is listed as one of the “grand challenges” for the application of wind energy in the latest review article on this subject by .
The study of the large-scale wind fluctuation dates back to the early 1950s , followed by several studies throughout the time, including , , , , and . The idea of the superposition of the large-scale and local turbulence was brought up by and argued in . demonstrated the validity of the theory in the spectral gap region through measurements. This strongly suggests that the large-scale and local turbulence are only weakly correlated, if correlated at all in the gap frequency range. A model of superposition of the two components is demonstrated to be able to reproduce the spectral behaviour across the microscale, spectral gap and mesoscale range. There have been several expressions in the literature for the large-scale wind fluctuation e.g. “two-dimensional turbulence” , “geostrophic turbulence” , “mesoscale turbulence” or simply “mesoscale fluctuation” . In this study, we follow and , and use the expression “two-dimensional (2D) turbulence”.
Thus, this study will take into account the contribution of the large-scale 2D turbulence, using the full-scale model provided by . In addition to the 3D and 2D turbulence, our calculation also includes the effect of stability, sea state and height dependence calibrated up to 200 m. The outcome of the method is a look-up table (LUT). We call the method developed here the “LUT” method.
The output is valuable for an initial evaluation of the cost related to the design of the wind turbine and the wind farm, particularly for measurement-sparse regions. The output variables include the mean characteristics of TI, the 90th percentile and the standard deviation of the TI-standard deviation σσ of TI.
The method is introduced in Sect. . Data used for demonstrate the method for an example global calculation are introduced in Sect. . Data and analysis from previous studies that are used for validation of the method are presented in Sect. . Results are shown in Sect. , followed by discussions and conclusions in Sects. and , respectively.
Method for calculation of turbulence intensity
To calculate the turbulence intensity (TI), we start with its definition here by
TI=σU/U‾,
where U‾ is the mean wind speed and σU is the standard deviation of the wind speed. Note that in practice, sometimes it is the standard deviation of the wind components rather than the standard deviation of the wind speed that is used to define the turbulence intensity.
In the absence of measured high-frequency time series of wind speed, we calculate σU by integrating the turbulence spectrum model for the wind speed from a lower cutoff frequency f1 to a higher cutoff frequency f2, so
σU=∫f1f2S(f)df1/2,
where S(f) is the turbulence power spectrum for wind speed, the modelling of which is introduced in Sect. . In boundary-layer turbulence, f1 is a value in the spectral gap, typically (10 min)−1 or (1 h)−1, while f2 is a value in the inertial dissipation range e.g. 10 Hz, corresponding to a time series sampled at 20 Hz, which is typical for sonic measurements.
The following subsections explain how both the effects of the typical boundary-layer turbulence (the 3D turbulence) and large-scale variability (the 2D turbulence) are included in the calculation of TI (Sect. ), how the stability effect is included in the algorithms (Sect. ), how the sea state dependence is introduced through parameterisation of z0 as a function of wave age cp/u* (Sect. ), how the 90th percentile of TI is calculated (Sect. ) and how the height dependence is added and calibrated (Sect. ).
These algorithms will then be run by using a set of corresponding parameters defined in specific ranges, thus constructing a final LUT with broad ranges of wind, wave and stability conditions combined. The set of parameters and their respective ranges are (1) wind speed at 10 m from 0.1 to 45 m s−1, (2) wave phase velocity at peak frequency cp from 0.1 to 30 m s−1, (3) stability parameter z/L at 10 m from -3 to 3, (4) 12 wind sectors and (5) height z=10, 50, 100, 150 and 200 m.
For application at a given site, with input of wind speed at a given height (Uz), wave phase velocity (cp) and stability (z/L at 10 m), we can find the corresponding TI from the LUT.
In this study, we use the fifth generation ECMWF reanalysis data (ERA5) , to demonstrate how these algorithms can be applied globally; the details and preparation of the ERA5 data are provided in Sect. .
Adding large-scale turbulence
As mentioned previously, usually when calculating variance and thereafter standard deviation and turbulence intensity, one integrates Eq. () from f1= (10 min)−1 or f1= (1 h)−1 to f2=10 Hz.
Figure illustrates the concept of the 3D and the 2D turbulence. The data in Fig. are exported from the study of (please refer to that study for the details of the data). The blue curves are the mean spectra from long-term sonic anemometer measured wind speed at 10 m (Fig. a) and 100 m (Fig. b). The Høvsøre site is located on the west coast of Denmark, where the winds from the west, namely the sea, dominate. The prevailing winds at 100 m represent the sea condition, while 10 m is sometimes affected by the underlying land. Over sea, the spectra will resemble more Fig. b, where the 3D turbulence is in general rather weak due to the smooth water surface, and therefore the relative contribution from 2D turbulence is bigger.
The dotted black lines show the typical boundary-layer 3D spectra for wind speed. For f<fp, with fp (the peak frequency of the boundary-layer spectrum of wind speed), the Kaimal model shape suggests a saturation of power spectrum at lower frequencies, shown as a 1:1 slope in the plot of fS(f) vs f in log-log coordinates:
fS1,u(f)=α102u*2fz/U(1+33fz/U)5/3,
where U is the mean wind speed, z is the height and α=1 by default. S1 is for the 3D turbulence, and here u is referring to the streamwise wind component, whose spectrum is very similar to that of the wind speed.
The dashed green curves are the spectral model for the large-scale variability of wind speed from :
fS2,U(f)=a1f-2/3+a2f-2,
where a1=3×10-4 m2 s−8∕3 and a2=3×10-11 m2 s−4 are derived from offshore climatological wind datasets. S2 is for the 2D turbulence.
The purple curve is a superposition of the Kaimal model and the large-scale variability from :
fSu(f)=f⋅S1,u(f)+S2,U(f).
The good agreement between the purple (model) and blue curves (measurements) suggests the success in bridging the 2D and 3D turbulence calculation using Eq. () – the full-spectrum model in the gap region.
The scales of 1 h and 10 min are marked in vertical grey lines in Fig. . If the integration of Eq. () uses the dotted black curves for the 3D turbulence only, using f1= (1 h)−1 gives slightly larger variance and hence slightly larger TI than using f1= (10 min)−1. However, the dotted black curves obviously miss out significant amounts of energy for f< (10 min)−1 compared to the blue curves, which are the measurements. Note that the curve from measurements is well described by the model (Eq. , the purple curve). In this study, we use the purple curve for the integration of the power spectrum, with f1=1 h−1 and f2=10 Hz, thus adding the large-scale turbulence.
Due to the similarity of the spectrum of the wind speed U and the along-wind component u, we refer to the wind speed in this study.
Illustration of the effect of large-scale wind variation on the variance from the gap region using data from Høvsøre (Fig. 8b, d from , with line colours edited). The power spectra of wind speed is plotted for (a) 10 m and (b) 100 m. The scale of 1 h and 10 min are marked with thin grey lines. The spectra from the measurements are shown in the blue curve, the 3D turbulence for f<fp is shown as the Kaimal model in black dots (the 3D turbulence for f>fp is from measurements), the dashed dark-green curve is the mesoscale spectral model from (here Eq. ) and the purple dash-dotted curve is the 2D + 3D turbulence model.
Adding stability effect
Atmospheric stability is introduced in the calculation through the wind speed profile:
U(z)=u*κlnz/z0-ψm,
where ψm is the stability function, dependent on the Obukhov length scale L=-u*2θ‾/(κgΘ*), with κ=0.41 (the Von Kármán constant), Θ*=w′θ′‾/u* and θ‾ the mean surface temperature in K.
For stable conditions z/L>0, we used
ψm=-5z/L,
and for unstable conditions z/L<0, we used
ψm=ln12(1+X2)(1+X2)2-2arctan(X)+π2,
where
X=(1-15z/L)1/4.
There are many recommendations in the literature on the coefficients that are used together with the stability parameter z/L in Eqs. () and (); see a summary in .
For neutral conditions ψm=0 and U=UN, where
UN(z)=u*κln(z/z0).
Thus the relationship between TI and TIN (turbulence intensity for the neutral condition) can be derived through Eqs. () and () accordingly:
TI=TIN⋅ln(z/z0)ln(z/z0)-ψm.
Adding ocean surface wave effect
The wave effect is added through the roughness length z0, which is used in the wind profile algorithm in Eq. (). The roughness of the water surface z0 is related to the sea state, water fetch and water depth.
In this study, we apply the simple sea state dependence, namely the wind sea, and focus on wind-generated waves in open sea conditions. There are many models to describe water roughness z0 through parameters of the sea state, such as significant wave height, wave length, wave steepness and wave age; see e.g. for a summary. Here we build the sea state dependence on the Charnock formulation, with consideration of the smooth water effect for z0:
z0=0.11ν/u*+αchu*2/g.
The first term on the right-hand side is the smooth flow effect due to viscosity, with ν the viscosity coefficient; it introduces a weak decrease of z0 with increasing wind speed for light wind conditions. The viscosity effect becomes negligible at moderate to strong winds. For the second term in Eq. (), we use the parameterisation scheme from through wave age cp/u* in connection with the Charnock parameter:
αch=a(cp/u*)b,
where a=0.023/1.0568U10 and b=0.012U10, with U10 the wind speed at 10 m. The reason for choosing the Fan scheme is based on the study of , in which they compared six parameterisation schemes with measurements at the Horns Rev site in the North Sea where the Fan scheme provides the calculation of z0 closest to the measurements.
The friction velocity u* can thus be obtained through iteration by using Eqs. (), () and ().
The wave phase velocity at the peak frequency cp is calculated through the wave dispersion relation through the peak wave number kp, considering the water depth h:
cp=gkptanh(kph).
Adding 90th percentile and standard deviation of TI
The 90th percentile of the distribution of σu, σu,90, at a given wind speed is a concept relevant for turbine design. The standard deviation of σu, σTI, at a given wind speed corresponds to the approximately 84.1th percentile.
In the IEC standard, σu,90,IEC=1.84⋅Iref is recommended to be the 90th percentile in relation to the mean value. Iref is the reference TI at a hub height wind speed of 15 m s−1, which is 0.18, 0.16, 0.14 and 0.12 for turbine type A+, A, B and C, respectively.
proposed the following expression to be the 90th percentile:
σu,90,W=0.0123⋅U+0.1221,
which suggests a narrower spread of TI at higher hub height wind speed because of the inverse dependence on U: σu,90/U=0.0123+0.1221/U.
also proposed an expression for σσu:
σσu,W=0.0108⋅U+0.1189,
whose magnitude is slightly smaller than σu,90.
Based on measurements from the FINO 1, 2 and 3 masts, derived the following expression for the standard deviation of TI varying with wind speed:
σTI,C=C1+C2⋅e-C3U.
They have provided detailed coefficients C1, C2 and C3, depending on the site (FINO 1, 2 or 3) and height. The difference between sites and height is minor, compared to the natural spread of TI from measurements. The estimate of σTI from is even larger than σu,90 from , suggesting the site-dependence characteristics of these empirical expressions.
We found with measurements from FINO1, 2 and 3 that the difference using the above algorithms for σu,90 and σTI is small compared with the scatter in the measurements. Both Eqs. () and () are suitable for adding variation to TI at a given wind speed.
Adding the height dependence of turbulence variance
Several equations mentioned in Sect. implemented the height dependence of TI. Measurements from the tall masts FINO 1, 2 and 3 suggest an average and approximate linear decrease of TI with a height of up to 100 m; see coloured dots and circles in Fig. a.
Over water, the height dependence of the power spectrum is rather simple for the 2D turbulence as the coefficients in Eq. () are invariant with height at normal measurement heights, following the analysis of the measurements at Horns Rev in . However, the height dependence of the 3D turbulence could be complicated. The Kaimal model provides a simple dependence of the power spectrum on z; however, its validity still needs to be verified with measurements in different ranges of wind speed and stability, as well as for heights above the surface layer. Thus, using the Kaimal model directly together with the 2D turbulence can sometimes generate a height dependence deviating significantly from the measurements and from theoretical relations that are based on measurements e.g. those in Fig. b. To solve this problem, we define a coefficient α in Eq. (), which is a function of both wind speed and height, for z≥50 m, to multiply Eq. (), so that our LUT results are consistent with theoretical relations based on measurements at FINO 1. Appendix provides a list of expressions of α for different wind speed ranges and different heights used for creating the LUT data.
The variances of wind speed and TI at heights 10, 50, 100, 150 and 200 m are obtained using the LUT. To compare with measurements from the three FINO masts, TI from LUT is organised for three wind speed ranges – 1–30, 1–33 and 1–37 m s−1 – which represent wind speed ranges in the three FINO sites. The TI values at z are normalised by those at 30 m (the approximate lowest measurement height at the FINO masts) and are shown in Fig. a (black curves). The height dependence from the ISO and its extended form (z/10)-0.22 is also shown in a similar manner in the same figure (blue curve). Due to different measurement periods, data coverage and meteorological conditions at the three sites, there are some differences between the results from the three locations. Nevertheless, the agreement is striking between them, and there is an overall good agreement between the measurements and the LUT data. The ISO and its extended expressions, with (z/10)-0.22, corresponds to a much stronger decrease of TI with height, particularly in the lowest 100 m, when compared with measurements. At higher elevations it is approaching a similar decreasing rate with height to the LUT results.
To verify our calculations for elevations higher than 100 m, data from a remote sensing technique (e.g. sounding or lidar measurements) can be used. Here we use derivations from two classical studies in and .
provided the following height dependence of wind variances:
σu2u*2=6(1-zzi)2+zziu′top2‾u*2,
where u′top2‾ is the variance at the top of the boundary layer. noted: “Although this ratio (u′top2‾/u*2) is expected to vary from situation to situation, during the KONTUR experiment it was found to equal 2.0 …”. To show Eq. () in the figure, we need to specify a number of variables. If we only show normalised variance σu(z) by the values at 10 m, then we only need to specify zi. Assumptions need to be made to estimate zi, and here we assume a neutral condition and apply the following expression for the atmospheric boundary layer zi=0.3u*/fc, with fc the Coriolis parameter. It should be noted that under neutral and stable conditions, the atmospheric boundary-layer height is often denoted by h, whereas zi commonly refers to the inversion height under convective conditions. Following , zi is often taken as an approximation of the planetary boundary-layer height h. For consistency with the cited formulations adopted here, we use zi as a unified notation for boundary-layer height across stability regimes. The coefficient 0.3 was used following , where no theoretical arguments were provided for why this coefficient was used; this expression nevertheless was suitable for interpreting derived measurement behaviours. In Fig. b, we plot σu(z)/σu(10m) for both u*=0.1 and 0.2 m s−1, to show the sensitivity. The use of a range of u* values brings the analysis to a qualitative level, thus also reducing the importance of the absolute value for the coefficient associated with zi.
provided the following theoretical expression for the height dependence of variance σu2 for the “very low wavenumbers”, specified as “range (iii)” in their study:
σu2u*2=2.5+ln0.46u*Aκzfc,
where A can be determined through A=fcUu*1nl and nl is the lower frequency limit of the surface eddy range in Hz. This “range (iii)” refers to eddies larger than the spectral gap and relevant for the entire boundary layer. With their measurements in the surface layer, showed the validity of the above expression in the surface layer down to a couple of metres above the ground. The validity of Eq. () at higher elevations requires more measurements to confirm. In the absence of information of stability, in Fig. b, we plot σu(z)/σu(10m) for both u*=0.1 and 0.2 m s−1, as well as for A=0.5 and 0.3, to show the sensitivity to these parameters.
Results from the four data sources – our LUT calculation, the ISO expression, Stull (1988) and Högström et al. (2002) – are shown in Fig. b, where the vertical distribution of the wind speed variance, normalised by the value at 10 m, σu(z)/σu(10m), is shown. As introduced before, both u*=0.1 and 0.2 m s−1 are used, and for Eq. (), A=0.3 and 0.5 are used, to show the variation in connection with the use of the two expressions. In general, deeper boundary layer (corresponding to larger u* and A) results in a slower decrease of wind variance with height. For comparison, to match the general wind conditions from and , we extracted the variance from the LUT data for wind speed at 10 m not exceeding 18 m s−1. The ISO expression of the height dependence does not distinguish between wind speed ranges; consistent with Fig. a, it provides the strongest decrease with height for z<100 m.
(a) TI normalised with TI at 30 m, with data from measurements at the three FINO sites (see data length from Table ), ISO expression and LUT method. (b) Comparison of the distribution of wind speed variance with height, normalised with value at 10 m, following the theoretical expressions of and .
Data
In this study, we use ERA5 data to demonstrate the calculation of global offshore turbulence intensity using the algorithms presented in Sect. . Relevant details of the ERA5 data are provided in Sect. .
At the same time, we use measurements to verify the calculations in the LUT method and use published results on the dependence of TI on U from several studies and sites. The corresponding details are provided in Sect. .
The base model data
The ERA5 data are chosen because of the global availability of both atmospheric and wave data. Other atmospheric and wave data can also be used.
In order to match the Global Wind Atlas data layers, including both resource and site parameters , the ERA5 reanalysis from the same period 2008–2017 (both 2008 and 2017 are included) was used here. ERA5 data are available with an hourly time resolution. The meteorological data (wind speed components, 2 m temperature, friction velocity) are available with horizontal spatial grid spacing of 0.25° × 0.25° (about 30 km in the mid-latitude), and the ocean wave data are only available on a 0.5° × 0.5° spatial grid. To avoid unnecessary interpolation of the wave data, our analysis is performed on the wave data grids only.
The following time series from the ERA5 data have been used: meridional and longitudinal wind components at 10 m (U10m, V10m), wave period at peak frequency (Tp), temperature at 2 m T2m, meridional and longitudinal stress components U′w′‾ and V′w′‾, and surface heat flux w′θ′‾. We also used static bathymetry data.
From the above data, the following variables are derived for 12 directional sectors for all water model grids: (1) mean wind speed at 10 m from the meridional and longitudinal components, (2) frequency of occurrence, (3) mean friction velocity u*=(U′w′‾2+V′w′‾2)1/4, (4) mean temperature scale (w′θ′‾/u*) and hence Obukhov length (L), and (5) mean wave phase velocity (cp) with water depth (h) effect considered through Eq. ().
Validation data
In this study, we use 10 min time series of wind speed and standard deviation from FINO 1, 2 and 3 at heights from 30 to about 100 m, to examine the calculation of TI, the 90th percentile and the standard deviation of TI, as well as the dependence of TI on height z. We chose the data periods before the wind farms around them began operations. The periods are specified in Table .
Some other measurements are not open source, like the FINO data, so we validate the LUT calculation using the published results on the dependence of the mean TI on the wind speed from , and .
Data details of these measurements are provided in Table .
Some details of the data from measurements used in the study. “TI–U distribution”: the distribution of mean values of TI in wind speed bins at a certain height, obtained from measurements in the reference studies. Please refer to the references for further details of the data used there. For FINO sites marked with *, we use the results directly from the corresponding reference study.
SiteHeightLatitude (° N),Data formReferencelongitude (° E)FINO 1100 m(54.0148, 6.5876)10 min time series (2004–2009)FINO 2100 m(55.0069, 13.1542)10 min time series (2009–2013)FINO 3100 m(55.1950, 7.1583)10 min time series (2010–2014)Høvsøre100 m(56.433, 8.15)TI–U distributionHorns Rev50 m(55.508, 7.875)TI–U distributionFrøya Sletringen46 m(63.6660, 8.2590)TI–U distributionFrøya Skipheia70 m(63.6680, 8.3270)TI–U distributionDogger Bank West83 m(55.0994, 2.7025)TI–U distributionDogger Bank East83 m(54.8670, 1.8200)TI–U distributionFINO3*81 m(55.1950, 7.1583)TI–U distributionFINO1*80 m(54.0148, 6.5876)TI–U distributionFINO2*82 m(55.0069, 13.1542)TI–U distributionIJmuiden91 m(52.8482, 3.4357)TI–U distributionKentish Flats80 m(52.6064, 4.3896)TI–U distributionLondon Array82 m(51.5850, 1.3940)TI–U distributionEgmond aan Zee70 m(51.4463, 1.0781)TI–U distributionResultsThe look-up table for TI
Using the algorithms presented in Sect. , a look-up table (LUT) is generated, where TI is a function of wind speed U, height z, wave phase velocity cp and stability z/L at 10 m. Figure shows an example of some distributions of TI from the LUT at 10 m. Neutral stability has been applied in Fig. a to c, and the stability effect is added in Fig. d.
The focus of Fig. a is to show the effect of waves on TI. Results from LUT using Eq. () in the black dots in Fig. a show the variation in TI at a given wind speed, caused by including the wave-age effect through cp/u*. At a given wind speed, for the growing waves, TI increases with cp/u*. The spread in TI is larger at stronger winds. At U10=15 m s−1, the difference caused by cp/u* in TI is on the order of 0.01. There is a weak decrease of TI with U in light winds due to the smooth flow effect, followed by an increase in TI with U up to approximately 31 m s−1. As can be seen in Fig. a, for winds stronger than 31 m s−1, for some groups, TI continues to increase with U but at a lower rate; most often, TI decreases with U again. We still need more measurements at very strong wind conditions to provide a more complete picture of how turbulence behaves there. In the Spectral WAve Nearshore (SWAN) model, the algorithm describing the relation between u* and U10 adopts the measurements collected in : CD=u*2/U102 is parameterised in terms of U10 following CD=(0.55+2.97(U10/31.5-1.49(U10/31.5)2)×10-3. In this expression at high winds U10>31.5 m s−1, CD saturates, followed by decrease in wind speed, as a simplified effect of the wave breaking process. This results in TI first increasing with U10 and shows similar dependence on the wind speed for U10>31.5 m s−1; see the thick blue curve in Fig. a. At the same time, recommended the following relationship between u* and U10: u*=0.239+0.0433(U10-8.271)+0.12(U10-8.271)2+0.181. Their study suggests that the above expression recognises the mechanisms by which heat and moisture cross the air–sea interface through molecular and micro-physical processes at the surface of sea spray droplets, and it has been validated with data for winds up to 25 m s−1 and extrapolated to hurricane wind strength. Compared to the SWAN expression, the Andreas expression suggests a weak increasing dependence of u* on U10 and, accordingly, when being used together with the Charnock formulation, an increasing dependence of TI with wind speed – shown as the dashed blue curve in Fig. a.
Figure a shows the calculation using only 3D turbulence (Eq. ). In Fig. b, 2D turbulence is added to the calculation of variance and hence TI. The effect of including the 2D turbulence is shown in Fig. b and c in pink dots. In Fig. c, the black dots are the same as in Fig. a; comparison of the black with the pink dots suggests that the effect of including 2D turbulence is most obvious at light to medium wind speeds, where TI decreases from much higher values with U up to about 7 m s−1. Measurements have consistently shown these high values of TI in low to moderate winds e.g. in , , and .
Based on the values in Fig. c, TI is calculated using Eq. (), with z/L varying from -3 (highest curve) to 3 (lowest curve). Figure d shows an example of the impact of stability at a given cp=10 m s−1.
Corresponding data and LUTs are also made available at 50, 100, 150 and 200 m. Values at heights in between can be interpolated.
An example of the distribution of the turbulence intensity (TI) with wind speed at 10 m. (a) Includes only micro-scale (3D) turbulence and shows difference of including wave-age dependence, using z0=0.011u*2/g together with the Andreas expression (dashed curve), and the SWAN algorithm (solid curve). (b) Similar to (a) but including the large-scale (2D) turbulence. (c) Shows the difference between including and excluding the 2D turbulence in the calculation of TI. (d) The green dots show the effect of stability from z/L=-3 to 3 for the group where cp=10 m s−1. The pink dots are the same as from (b) and (c), showing the effect of wave age and 2D turbulence but for neutral conditions.
The global offshore atlas of turbulence intensity
For a given offshore location, if we can get basic information on z/L and cp, we can use the LUT and identify TI as a function of the wind speed at five heights: z=10, 50, 100, 150 and 200 m. In the absence of certain information, one can always assume neutral stability, skip the wave-age dependence, and use e.g. the SWAN expression or the Andreas expression (e.g. the blue curves in Fig. b); the associated uncertainty needs to be assessed, however.
With additional information on the wind speed frequency distribution as a function of each wind condition, we can calculate the corresponding TI distribution at the site. The frequency of occurrence can also be used to weight the calculation in obtaining the mean values of TI.
The information about stability, waves and wind speed distribution can be obtained either from measurements or from modelling. We have prepared global data for these variables based on the ERA5 data. Here we prepared mean values of z/L, cp and wind speed in 12 sectors, along with their occurrence frequencies.
Figure a and b show an example of the global distribution of offshore mean values of TI at ERA5 grid points at 10 and 100 m, respectively. To better show the data variability, TI in the range of [0.04, 0.14] is plotted for 10 m (Fig. a), and TI in the range of [0.03, 0.1] is plotted for 100 m (Fig. b). Here, one can see systematically larger TI at 10 m than at 100 m. Another striking feature is the high TI in the tropical regions; here, based on the ERA5 data, the wind speed is low, and z/L is negative and of large magnitude, suggesting very convective conditions. Thus, these large values of TI are a combination of effect from z/L (see Fig. d) and the relatively large 2D turbulence contribution at such low winds (see Fig. c). Measurements from this region could be of great value to verify this.
Global atlases of the mean turbulence intensity at (a) 10 m (only values in the range [0.04, 0.14] shown) and (b) 100 m (only values in the range [0.03, 0.1] shown) for offshore grid points, including 2D turbulence, wave age and stability effect. ERA5 data during the period of 2008–2017 are used.
Validation at sitesTI–<inline-formula><mml:math id="M267" display="inline"><mml:mi>U</mml:mi></mml:math></inline-formula> distribution
We compare our calculation with the results from at 11 offshore stations, one from and one from .
First we extracted data from ; their Fig. 5 for the 11 sites and plotted here in Fig. a–k, including measurements (red dots) and their calculation using Eq. () (dashed blue curves). The mean TI distribution with wind speed based on measurements from Høvsøre and from Horns Rev are plotted in Fig. l and m, respectively. The site names and data heights are shown in the subplots' labels. With the given locations of these sites (Table ), we extracted the corresponding data from the ERA5-derived atlas as shown in Sect. , including z/L, cp, wind speed and occurrence frequency in 12 sectors. These data were used to identify the corresponding TI from the LUT. Three of the 11 sites from – Frøya Sletringen, Frøya Skipheia and Egmond aan Zee – correspond to ERA5 land grid points as there are small islands nearby. We used the neighbouring water grid point to represent each of the three sites.
The LUT provides data at five heights: 10, 50, 100, 150 and 200 m. To obtain the LUT results at the corresponding heights for the 13 sites, we used the data at two neighbouring heights and applied linear interpolation. The results are shown in solid black curves in Fig. .
Our calculation shows overall good agreement with both measurements and Jeans' calculation using Eq. () at these sites, and it does not depend on wind speed range. The LUT captures well the variation of TI with U for low to moderate winds, suggesting the success of including the 2D turbulence. On average, the LUT provides an improvement of TI of about 20 % for wind speed greater than 8 m s−1, which is on the order of 0.01, in comparison with Eq. (), with a larger TI value. In this wind speed range, with the mean TI magnitude less than 0.1, a difference on the order of 0.01 is in general not important, as it does not change the decision of the turbine class. The associated impact on fatigue load is also only a few percentage under normal conditions but could be non-negligible under special conditions such as yaw misalignment. However, such a difference is systematic and should not be ignored.
Among these sites, both the LUT and the Jeans algorithm underestimate TI at Høvsøre. Note that the LUT and the Jeans algorithm assume the site open sea condition. Høvsøre is a coastal site located on land, a few kilometres away from the shoreline. At a height of 100 m, it is assumed that the flow from the ocean represents offshore conditions. The measurement data, shown as red dots in Fig. l, were from the water fetch. However, it is not clear if and how much it is affected by the presence of land, which could have contributed to the discrepancies here.
Comparison of TI from the LUT method (black curves), with values from measurements (red dots) and Eq. () from (dashed blue curves). The measurements are extracted from (a–k), from (l) and from (m). The names of the sites are (a) Frøya Sletringen, (b) Frøya Skipheia, (c) Dogger Bank West, (d) Dogger Bank East, (e) FINO3, (f) IJmuiden, (g) FINO1, (h) FINO2, (i) Kentish Flats, (j) London Array, (k) Egmond aan Zee, (l) Høvsøre and (m) Horns Rev. The corresponding heights of measurements used in the plot are shown in the subplots.
90th percentile of TI
Using the mean values of wind, wave parameter and stability in 12 sectors derived from ERA5 data as input to the LUT method does not provide a representative scatter of TI. The LUT method here adopts the algorithms from the literature for describing the 90th percentile and the standard deviation of TI. Figure shows an example for the site FINO 1. The mean and 90th percentile of TI at wind speed bins are calculated from measurements, from the IEC standard and from the algorithm provided by . In addition, we also calculated σTI at these wind speed bins using the algorithm of through measurements from FINO sites. Both the algorithms from and can reasonably well describe the 90th percentile of TI from the measurements here. However, at 30 and 50 m, at wind speeds lower than 10 m, the measurements suggest an even larger variation of TI, which could be partly caused by the fact that we did not take the flow distortion of the tower base into consideration. This effect is clearly not so obvious at 100 m. In Fig. a, the mean TI–U curves from the LUT are presented for the FINO 1 site, from 10 to 200 m.
Presentation of TI from the LUT method at FINO 1. (a) Mean values of TI varying with wind speed at six heights from 10 to 200 m. (b) FINO 1 at 30 m, showing mean TI (measurements and LUT), 90th percentile of TI from measurements (orange curve), and derived from Eq. (dashed teal curve) and σTI (Eq. (dotted blue curve). (c) FINO 1, similar to (b) but at 50 m. (d) FINO 1, similar to (b) but at 100 m.
Discussion
This study provides a method as well as an example database for turbulence intensity for offshore conditions globally, from 10 to 200 m, which are heights relevant for existing turbine types. The data include sector-wise distributions of mean as well as standard deviation of TI as a function of wind speed from calm to gale wind conditions. Note that the current study considers TI in the absence of offshore wind farms. For this reason, direction measurements from the three FINO masts are only used during the period before their surrounding wind farms began operations.
From this method, an LUT is created, where TI can be identified with input of wind speed (Uz), stability (z/L at 10 m), wave phase velocity (cp) and height (z). The interface for input (Uz, z/L, cp and z) allows the use of data from any source for preparing these variables. It could be from measurements, modelled data or a combination.
Two major contributions from the LUT method are particularly worth highlighting in comparison with previous studies.
The first is the inclusion of the large-scale two-dimensional turbulence, which results in the decreasing dependence of TI with wind speed for lower wind speed ranges. Previously, this part was taken care of by imposing an empirical and artificial dependence e.g. (z-0.22) or an assumption of climatological unstable conditions at low wind speed . The large values of TI observed in light winds are often associated with convective conditions. Although we included the stability effect through ψ(z/L) as in Eq. (), this effect through the sector-wise mean ERA5 values is not obvious. Note that we did not include the stability effect in the derivation using the turbulence spectrum, which is used to calculate σU in TIN=σU/UN. Figure c shows that the most significant effect from the 2D turbulence on TI is for small to medium wind speeds. At strong winds, the surface-driven turbulence that is generated using e.g. the Kaimal spectral model under the neutral condition is significant in comparison to the 2D turbulence. At weak winds, the calculated 2D turbulence becomes relatively more important as the Kaimal model generates relatively little turbulence. The inclusion of large-scale wind variability through the 2D turbulence model for weak to medium wind speeds thus compensates for the lack of consideration of stability in the 3D spectrum.
The second is the systematic calibration of height dependence of the calculation of TI up to 200 m with historic-measurements-based studies and measurements from FINO 1. The outcome is satisfactory at the validation sites (Figs. a, and ). The height dependence such as in the ISO standard (z-0.22) seems to be too strong in the lowest 200 m; see Fig. a and b. Note that the validation of TI in is done mostly at one height.
The LUT method is based on physical principles; it implements algorithms that have been derived from theory and measurements, and validated with measurements. Nevertheless, to carry out the calculation for a global coverage, we used many assumptions and simplified solutions in the chain of algorithms and data. In the following, we discuss the limitations and associated uncertainties in several key assumptions.
First of all, the LUT method includes atmospheric stability effect explicitly through Monin–Obukhov similarity theory (MOST), through Eq. (). It is known that stability affects the power spectra of wind, and hence the standard deviation of wind speed. The current study did not include any stability effect explicitly in describing the power spectra of wind speed, because it is a very challenging task, particularly for stable conditions. It is relatively easier for unstable conditions due to more organised, larger eddy sizes. suggested an extended spectral model for the unstable condition using a scaling parameter – the boundary-layer height zi. However, this model has not been tested for scales larger than the gap range. It was shown through a case with the presence of convective open cells in . The model using zi is of very limited use in describing the large-scale wind fluctuation, whereas combining 2D and 3D turbulence, in its simple manner, better describes the power spectrum of wind speed for the convective condition. Over water, the stability is less influenced by the diurnal cycles compared to land conditions, except in the coastal zones where land effect is present. Over water, it is mostly influenced by the advection of air masses from another region or from land, associated with weather conditions that likely vary with season. We expect smaller uncertainties associated with stability calculation at strong winds, where it is often close to neutral conditions.
Second, the LUT method uses the Fan scheme to include the sea state effect through the roughness length z0, which is a function of wave age cp/u* and wind speed. The use of the Fan scheme successfully describes the increasing dependence of TI on wind speed for moderate to strong winds, here validated with measurements up to about 30 m s−1. Note that like most parameterisation schemes for z0, wave conditions are simplified using the Fan scheme, here Eqs. () to (). For instance, the fetch effect on the wave field is not considered. Only the most dominant wave with a peak frequency is considered when cp is used. The increasing dependence of z0 on the wave age is applicable for wind-generated winds while not swell conditions. For a swell-dominated sea state, the waves are detached from the local wind, the wave age cp/u*>25 or cp/U10>1.2, and the corresponding wind speed is often weak. Thus, our neglect of the impact of the swell in our algorithms mainly contributes to the uncertainties in weaker wind conditions. Through the dependence of z0 (or interchangeably u* and drag coefficient CD), our LUT provides the variation of TI with U, including the two schemes: the SWAN and the Andreas schemes. Both schemes are well established, and they use a simple dependence of u* (or CD) on U10. Even though the derivation of both addressed possible physics of air–sea interaction, they show contradicting results at the very strong wind conditions (Fig. a and b). These reflect the varying challenges in modelling the effect of sea state on atmospheric turbulence in different conditions.
Third, the LUT method calibrates the height dependence of the Kaimal model using measurements from the FINO 1 site and applies it for general use. At the same time, we apply Eq. () and MOST to calculate the height dependence of the mean wind speed. This is an oversimplification as Eq. () and MOST are expected to be useful in the atmospheric surface layer. Large uncertainty is therefore expected for stable conditions, where the height of the surface layer is very low. Even though the above height dependence descriptions contributed to quite good results at 13 stations (Fig. ), as well as good agreement with the theoretical expressions from and , we need more measurements from worldwide locations to provide information for further improvement to reduce uncertainties.
In this study, we used the modelled data for a demonstration of the use of the LUT and hence provided an example of a global dataset. We used the ERA5 data, from which Uz, z/L, cp and z can easily be prepared. One could also use e.g. the New European Wind Atlas (NEWA) data to prepare atmospheric parameters Uz,z/L and z, and use another source for wave parameter cp. This would of course bring additional uncertainty due to inconsistency in the data sources. Nevertheless, it is worth trying and being validated for individual cases, as a complete set of both atmospheric and wave data is not easy to obtain. If the wave data are absent or proven to be unreliable, one could also neglect the wave-age dependence in the calculation. This corresponds to the solid blue curves in Fig. b, instead of the pink dots.
The ERA5 data have a global coverage of both the atmospheric and wave data but are of relatively coarse spatial resolution. This leads to higher uncertainties in coastal zones where the detailed spatial features are not resolved. However, in the LUT method calculations, the ERA5 data can be straightforwardly replaced by other model data with a higher spatial resolution, when possible.
ERA5 surface fluxes are used to calculate the Obukhov length L, which is later used to convert TI from neutral conditions to the corresponding stability conditions. All 13 validation sites are in the North Sea region; the combined data show that there is a dependence of the mean values of z/L on the sectors, with the majority of the data in the range -0.1<z/L<0 in most sectors and 0<z/L<0.05 in southwest sectors. In general, the averaged stability effect is not significant at the sites shown here, at least from the ERA5 data, but it may be important at other locations.
At the same time, in connection with the use of the global ERA5 wave data, the wave parameters are only available approximately every 50 km for the sites we studied here. This implies that most of the ERA5 data are quite far from the coastlines and are likely more representative for open ocean conditions.
For a global calculation, using a look-up table simplifies the calculation significantly and saves computation. The mean statistics of TI–U relationships and its dependence with height is well captured by the LUT method. To handle the large global ERA5 dataset, only the mean values from 12 directional sectors are stored for wind speed, cp and L. The spread of TI, associated with wind, stability, wave, fetch and their sector-wise dependence is thus underestimated by using the sector-wise mean values of the ERA5 data. A simple solution is taken here preliminarily through the use of the existing algorithms for describing the standard deviation of TI or the 90th percentile.
In the future, other data with a higher resolution than the ERA5 data can be used to provide the input to apply the LUT method. One could also include more advanced statistics than the sector-wise mean values as input to the LUT method, thus improving the calculation of the spread of TI.
We acknowledge that here the validation of the LUT has only been done to 13 stations in the European seas, which are dominated by similar weather conditions. With the open source codes and data from this study, contributions of data validation from other regions in the world are expected, which will help in identifying issues and finding solutions. For instance, it would be of great interest to see how the LUT performs in areas affected by tropical cyclones, as well as places with special prevalent phenomenon such as water sprouts.
Conclusions
This study delivers the look-up table (LUT) method for the turbulence intensity (TI) of offshore conditions. The method describes TI through both the two- and three-dimensional turbulence, making it more suitable for heights above the atmospheric surface layer. At weak wind conditions, this systematically includes the effect of large eddies through the two-dimensional turbulence, which otherwise is neglected in using the three-dimensional turbulence only, under the assumption of neutral stability and weak wind. Accordingly, it captures the decreasing dependence of TI for light to moderate winds. This method describes the roughness length z0 through dependence on wind speed and wave age, which accurately describes TI for moderate to strong winds, including wave breaking effect. Stability effect is also included. The method for TI can be applied globally and was calibrated to a height of 200 m.
The global data prepared using this method together with the global ERA5 wind and wave data were shown to perform reasonably well with validations from 13 offshore sites and comparison with earlier studies. With a given location, we can find the corresponding ERA5 coordinates and extract the corresponding data. The LUT will then provide TI according to this information at the above-mentioned five heights. If the height of interest is in between those five heights, one can interpolate between neighbouring heights.
This method is flexible to use if one has special data as input. In the absence of other data sources, one can use the global ERA5 data prepared here by providing the coordinate for the site of interest.
Collection of expressions of TI
The detailed expressions for TI from and are provided here. Specifically, the equation of TI from reads
TIWang=κCNT2+CMX2(-zi/L)2/31/2ln(z/z0)-ψm(z/L),
where CNT and CMX are two constants as approximates for σ/u* and σ/w*, respectively. For neutral and mixed-layer conditions, u* is the friction velocity and w* is the convective velocity scale, and it is a function of temperature, surface heat flux and height of the boundary layer zi, with L as the Obukhov length.
There are several similar expressions in , primarily for the wind speed range from 10 to 26 m s−1. This includes the “modified Vickery model”
TIVic=0.085U100.421z10-0.2,
the “drag coefficient model”
TIDra=0.08571+0.758U10-10.5z10-0.2
and the “linear model”, which has already been introduced as Eq. () in Sect. .
The three expressions (Eqs. , and ) provide similar estimates up to a wind speed of 26 m s−1, which describe well the measurements at the Frøya site from about 10 to 46 m. showed the validity of Eq. () from 8 m s−1. Above 26 m s−1, there are no measurements, and the linear model gives considerably larger values, which was recommended by the authors for a slightly conservative design approach.
Calculating the height dependence of variance
In the process of calibration using measurements from FINO 1, we added a coefficient α to the Kaimal model expression for u for different z, starting from 50 m and above. The coefficient α is a function of wind speed at each height.
B1α50m=0.018U+0.54,U≤10ms-1B2=0.024U+0.48,10ms-1<U≤32ms-1B3=1.29,U>32ms-1B4α100m=0.035U+0.037,U≤35ms-1B5=1.26,U>35ms-1B6α150m=0.031U+0.033,U≤35ms-1B7=1.12,U>35ms-1B8α200m=0.029U+0.031,U≤35ms-1B9=1.05,U>35ms-1
Code and data availability
Data and codes are published at 10.11583/DTU.30575555.v1.
Author contributions
XL designed the study, developed the method, did the calculation of TI, and outlined and wrote the paper. MI contributed with the preparation of global ERA5 data for the wind and wave statistics, and editing the paper. RF contributed with the preparation of global ERA5 data for the stability parameters and editing the paper.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
Fino data were made available by the FINO (Forschungsplattformen in Nord- und Ostsee) initiative, which was funded by the German Federal Ministry of Economic Affairs and Climate Action (BMWK) on the basis of a decision by the German Bundestag, organised by the Projekttraeger Juelich (PTJ) and coordinated by the German Federal Maritime and Hydrographic Agency (BSH). The ERA5 data have been provided by the Copernicus Climate Change Service Climate Data Store .
Financial support
This research has been supported by the Horizon Europe Climate, Energy and Mobility (DTWO project, grant no. 101146689) and the Energiteknologisk udviklings- og demonstrationsprogram (GASPOC project, grant no. 65020-1043 and IDEA project 134-21029).
Review statement
This paper was edited by Etienne Cheynet and reviewed by two anonymous referees.
ReferencesAndersen, O. J. and Løvseth, J.: The Frøya database and matitime boundary layer wind description, Marine Structures, 10.1016/j.marstruc.2006.07.003, 2006.Andreas, E. L., Mahrt, L., and Vickers, D.: An improved bulk air-sea surface flux algorithm, including spray-mediated transfer, Q. J. Roy. Meteor. Soc., 141, 642–654, 10.1002/qj.2424, 2014. Arya, P. S.: Introduction to Micrometeorology, Academic Press, ISBN 0-12-064490-8, 2001.C3S: ERA5 hourly data on single levels from 1940 to present, Copernicus Climate Change Service (C3S) Climate Data Store, 10.24381/CDS.ADBB2D47, 2018. Charney, J. G.: Geostrophic turbulence, J. Atmos. Sci., 28, 1087–1095, 1971.Cheynet, E., Jakobsen, Jasna, B., and Reuder, J.: Velocity spectra and coherence estimates in the marine atmospheric boundary layer, Bound.-Lay. Meteorol., 169, 429–460, 10.1007/s10546-018-0382-2, 2018.Cheynet, E., Li, L., and Jiang, Z.: Metocean conditions at two Norwegian sites for development of offshore wind farms, Renewable Energy, 224, 120184, 10.1016/j.renene.2024.120184, 2024.Christakos, K., Mathiesen, M., Henrik, O., Holvik, S., and Meyer, A.: Turbulence Intensity Model for Offshore Wind Energy Applications, 10.13140/RG.2.2.29810.88000, 2016.Davis, N. N., Badger, J., Hahmann, A. N., Hansen, B. O., Mortensen, N. G., Kelly, M., Larsén, X. G., Olsen, B. T., Floors, R., Lizcano, G., Casso, P., Lacave, O., Bosch, A., Bauwens, I., Knight, O. J., van Loon, A. P., Fox, R., Parvanyan, T., Hansen, S. B. K., Heathfield, D., Onninen, M., and Drummond, R.: The Global Wind Atlas: A High-Resolution Dataset of Climatologies and Associated Web-Based Application, B. Am. Meteorol. Soc., 104, E1507–E1525, 10.1175/BAMS-D-21-0075.1, 2023.Dörenkämper, M., Olsen, B. T., Witha, B., Hahmann, A. N., Davis, N. N., Barcons, J., Ezber, Y., García-Bustamante, E., González-Rouco, J. F., Navarro, J., Sastre-Marugán, M., Sīle, T., Trei, W., Žagar, M., Badger, J., Gottschall, J., Sanz Rodrigo, J., and Mann, J.: The Making of the New European Wind Atlas – Part 2: Production and evaluation, Geosci. Model Dev., 13, 5079–5102, 10.5194/gmd-13-5079-2020, 2020. Fan, Y., Lin, S., Held, I., Yu, Z., and Tolman, H.: Global ocean surface wave simulation using a coupled atmosphere-wave model, J. Climate, 25, 6233–6252, 2012.Grant, A.: Observations of boundary layer structure made during the 1981 KONTUR experiment, Q. J. Roy. Meteor. Soc., 112, 10.1002/qj.49711247314, 1986.Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, 10.1002/qj.3803, 2020. Högström, U.: Review of some basic characteristics of the atmospheric surface layer, Bound.-Lay. Meteorol., 78, 215–246, 1996.Högström, U., Hunt, J. C., and Smedman, A.: Theory and measurements for turbulence spectra and variances in the atmospheric neutral surface layer, Bound.-Lay. Meteorol., 103, 101–124, 10.1023/A:1014579828712, 2002. Højstrup, J.: Velocity spectra in the unstable boundary layer, J. Atmos. Sci., 39, 2239–2248, 1982. IEC: IEC 61400-1 Ed4: Wind turbines - Part 1: Design requirements, standard, International Electrotechnical Commission, Geneva, Switzerland, 2019.ISO: Petroleum and natural gas industries – Specific requirements for offshore structures – Part 1: Metocean design and operating considerations, ISO 19901-1:2015, https://www.iso.org/standard/60183.html (last access: 11 May 2025), 2015.Jeans, G.: Converging profile relationships for offshore wind speed and turbulence intensity, Wind Energ. Sci., 9, 2001–2015, 10.5194/wes-9-2001-2024, 2024. Kaimal, J., Wyngaard, J., Izumi, Y., and Coté, O.: Spectral characteristics of surface-layer turbulence, Q. J. Roy. Meteor. Soc., 98, 563–589, 1972. Kim, K. and Adrian, R.: Very large-scale motion in the outer layer, Phys. Fluids, 11, 417–422, 1999.Kosović, B., Basu, S., Berg, J., Berg, L. K., Haupt, S. E., Larsén, X. G., Peinke, J., Stevens, R. J. A. M., Veers, P., and Watson, S.: Impact of atmospheric turbulence on performance and loads of wind turbines: knowledge gaps and research challenges, Wind Energ. Sci., 11, 509–555, 10.5194/wes-11-509-2026, 2026. Kraichnan, R.: Inertial ranges in two-dimensional turbulence, Phys. Fluids, 10, 1417–1423, 1967.Larsén, X. G., Vincent, C. L., and Larsen, S.: Spectral structure of the mesoscale winds over the water, Q. J. Roy. Meteor. Soc., 139, 685–700, 10.1002/qj.2003, 2013. Larsén, X. G., Larsen, S. E., and Petersen, E. L.: Full-scale spectrum of boundary-layer winds, Bound.-Lay. Meteorol., 159, 349–371, 2016.Larsén, X. G., Du, J., Bolanos, R., Imberger, M., Kelly, M., Badger, M., and Larsen, S. E.: Estimation of offshore extreme wind from wind-wave coupled modeling, Wind Energy, 22, 1043–1057, 10.1002/we.2339, 2019.Larsén, X. G., Larsen, S. E., Petersen, E. L., and Mikkelsen, T. K.: A model for the spectrum of the lateral velocity component from mesoscale to microscale and its application to wind-direction variation, Bound.-Lay. Meteorol., 178, 415–434, 10.1007/s10546-020-00575-0, 2021.Larsén, X., Davis, N., Hannesdóttir, Á., Kelly, M., Svenningsen, L., Slot, L., Imberger, M., Olsen, B., and Floors, R.: The Global Atlas for Siting Parameters project: Extreme wind, turbulence, and turbine classes, Wind Energy, 10.1002/we.2771, 2022.Larsén, X. G., Imberger, M., and Floors, R. R.: Global Offshore Turbulence Intensity, Technical University of Denmark [data set], 10.11583/DTU.30575555.v1, 2026. Lindborg, E.: Can the atmospheric kinetic energy spectrum be explained by two-dimensional turbulence?, J. Fluid Mech., 388, 259–288, 1999.Mann, J.: The spatial structure of neutral atmospheric surface-layer turbulence, J. Fluid Mech., 273, 141–168, 10.1017/S0022112094001886, 1994. Muñoz-Esparza, D., Kosovic, B., Mirocha, J., and van Beeck, J.: Bridging the transition from mesoscale to microscale turbulence in numerical waether prediction models, Bound.-Lay. Meteorol., 153, 409–440, 2014.Nastrom, G. and Gage, K.: A climatology of atmospheric wavenumber spectra of wind and temperature observed by commercial aircraft, J. Atmos. Sci, 42, 950–960, 1985. Nastrom, G., Gage, K., and Jasperson, W.: Kinetic energy spectrum of large- and mesoscale atmospheric processes, Nature, 310, 36–38, 1984. Panofsky, H. and der Hoven, I. V.: Spectra and cross-spectra of velocity components in the mesometeorological range, Q. J. Roy. Meteor. Soc., 81, 603–606, 1955.Peña, A., Floors, R., Sathe, A., Gryning, S.-E., Wagner, R., Courtney, M., Larsén, X., Hahmann, A., and Hasager, C.: Ten Years of Boundary-Layer and Wind-Power Meteorology at Høvsøre, Denmark, Bound.-Lay. Meteorol., 158, 1–26, 10.1007/s10546-015-0079-8, 2016.Stull, R.: An Introduction to Boundary Layer Meteorology, Springer Dordrecht, 10.1007/978-94-009-3027-8, 1988.Tai, S.-L., Berg, L. K., Krishnamurthy, R., Newsom, R., and Kirincich, A.: Validation of turbulence intensity as simulated by the Weather Research and Forecasting model off the US northeast coast, Wind Energ. Sci., 8, 433–448, 10.5194/wes-8-433-2023, 2023. Veers, P.: Three-Dimensional Wind Simulation, SANDIA REPORT, OSTI ID: 7102613, 1988.Wang, H., Barthelmie, R. J., Pryor, S. C., and Kim, H. G.: A new turbulence model for offshore wind turbine standards, Wind Energy, 10.1002/we.1654, 2014. Zijlema, M., Van Vledder, G., and Holthuijsen, L.: Bottom friction and wind drag for wave models, Coast Eng., 65, 19–26, 2012.