Modeling of the Wind Potential in the Open Sea and Its Application to the Calculation of Energy

: A reliable estimate of the wind potential in the marine atmospheric boundary layer (MABL) is of great importance to justify the energetic viability of new offshore wind farms. The purpose of the study is to provide an additional tool for the prediction of the energy that a wind turbine would produces in the open sea from the usual way of measurements at sea, that is, when they are carried out with measuring masts, where the meteorological data are obtained at levels much lower than those of a wind turbine hub. For this, the variation in the wind speed with the height in the MABL is determined, based on the Monin–Obukhov similarity theory, according to the boundary conditions of the air–sea interface, where the input data for the Validation of the results are extracted from the German FINO 3 research platform during the years 2016, 2017, and 2018. It is applied to the production of electrical energy from a 6.0 MW commercial wind turbine, with the hub at 100 m above the sea surface. As a more prominent result, the deviations from the proposed method do not exceed 2.5% in the energy calculation. investigation, M.I.B.G. and S.P.L.; resources, M.I.B.G.; data curation, M.I.B.G.; writing—original draft preparation, M.I.B.G. and J.M.M.; writing—review and editing, M.I.B.G. and C.R.G.; supervision, M.I.B.G. and S.P.L. All authors have read agreed to the version of the


Introduction
Offshore wind energy, with a total installed capacity at the end of 2020 of 35,300 MW [1], has had a great expansion in recent years. In this field, worldwide, there are large unexploited marine extensions with moderate to strong winds, which require a reliable method of predicting wind potential, with the aim of studying the electrical energy that would be produced by the offshore wind farms to be installed in the future.
Although the wind of the atmospheric surface layer is usually simplified considering that it is a two-dimensional vector magnitude defined by the speed and the direction [2]. The modeling of the variation in the wind speed with the height in short periods of time presents serious difficulties as there is not always continuity in growth, sometimes it even decreases with height [3]. On the other hand, the same does not happen with the direction of the wind, whose value remains practically constant at different levels of the surface layer.
Currently, the installation and maintenance costs of offshore wind farms are higher than their onshore counterparts, and this difference can be compensated for by the better characteristics of the wind at sea. However, a deeper knowledge is necessary to determine the long-term wind potential that can be translated into new decision tools. Along this line, it must be considered that the surface layers of the land and the sea are different, both in terms of quantity of movement exchanges, as well as in heat and moisture fluxes [4], and, also because the sea surface presents a constant change having a continuous length of aerodynamic roughness due to swell [5].
On the other hand, the Monin-Obukhov (M-O) [9] similarity theory contemplates the MABL boundary conditions. Thus, under conditions of neutral atmospheric stratification, in relation to the momentum, Equation (2) represents a logarithmic profile of the variation in the wind speed with the height. U(z) = u * /k ln(z/z 0 ) (2) where k is the von Kármán constant, which usually takes a value of 0.4, z 0 is the aerodynamic roughness length and u * is the friction velocity, which is defined by the following expression: where τ is the surface tension, which acts parallel on the surface of the sea, ρ is the air density and u 2 * is the kinematic tension. Moreover, based on this theory, the diabatic profile of the wind is given by: U(z) = u * /k[ln(z/z 0 ) + Ψ M (z/L)] (4) where L is the Obukhov length [10], expressed as: where g/θ v is the buoyancy in m/s 2 K, and the heat kinematic flux, F Hs , in K m/s will be changing during the daily cycle and we will use the approximate expression [11]: where T a is the air temperature, T s is the water temperature on the sea surface, and C H is the coefficient of main transfer of sensible heat. The value of said dimensionless coefficient ranges between 0.001 and 0.005 for neutral conditions, takes a value of 0.006 for unstable atmosphere, and varies between 0.0009 and 0.000001 for stable conditions [12]. The atmospheric stratification processes in the surface layer are given by the empirical function, Ψ M (z/L) [3]; thus, for neutral stratification, z/L = 0, the function Ψ M (z/L) is zero in Equation (4). Instead, for unstable atmosphere, z/L < 0, the function Ψ M (z/L) is represented by the Paulson equation [13]: Energies 2022, 15 For stable stratification, z/L > 0, and the function Ψ M (z/L) is represented by the Businger-Dyer equation [14,15]: Graphically, the wind speed profile with neutral atmospheric stratification, Equation (2), appears as a straight line when represented with a semi-logarithmic scale. Whereas, according to Equation (4), for the non-neutral cases, the profile of the wind speed with the height deviates slightly from the logarithmic one, in such a way that under conditions of stable atmospheric stratification the wind profile is concave downward, while under conditions of unstable atmospheric stratification it is concave upward. That is, with respect to neutral stratification, the wind speed increases in the first case (stability) and decreases in the second case (instability).
Graphically, in Figure 1a, the wind speed profile with neutral atmospheric stratification, Equation (2), appears as a straight line when it is represented with a semi-logarithmic scale. Whereas, according to Equation (4), for the non-neutral cases, the profile of the wind speed with the height deviates slightly from the logarithmic one, in such a way that under conditions of stable atmospheric stratification the wind profile is concave towards below, while under conditions of unstable atmospheric stratification, it is concave upwards. That is, with respect to neutral stratification, the wind speed increases in the first case (stability) and decreases in the second case (instability). Likewise, Figure 1b represents the previous equations of the similarity theory and the power law (1), which presents a profile at higher elevations, from the measurement elevation, with less slope than the logarithmic equation under stratification conditions neutral, giving higher values of wind speeds for the same height.
For stable stratification, 0 L / z  , and the function ) L z ( M  is represented by the Businger-Dyer equation [14,15]: Graphically, the wind speed profile with neutral atmospheric stratification, Equation (2), appears as a straight line when represented with a semi-logarithmic scale. Whereas, according to Equation (4), for the non-neutral cases, the profile of the wind speed with the height deviates slightly from the logarithmic one, in such a way that under conditions of stable atmospheric stratification the wind profile is concave downward, while under conditions of unstable atmospheric stratification it is concave upward. That is, with respect to neutral stratification, the wind speed increases in the first case (stability) and decreases in the second case (instability).
Graphically, in Figure 1a, the wind speed profile with neutral atmospheric stratification, Equation (2), appears as a straight line when it is represented with a semi-logarithmic scale. Whereas, according to Equation (4), for the non-neutral cases, the profile of the wind speed with the height deviates slightly from the logarithmic one, in such a way that under conditions of stable atmospheric stratification the wind profile is concave towards below, while under conditions of unstable atmospheric stratification, it is concave upwards. That is, with respect to neutral stratification, the wind speed increases in the first case (stability) and decreases in the second case (instability). Likewise, Figure 1b represents the previous equations of the similarity theory and the power law (1), which presents a profile at higher elevations, from the measurement elevation, with less slope than the logarithmic equation under stratification conditions neutral, giving higher values of wind speeds for the same height.

Aerodynamic Roughness Length
The variation in wind speed with height through a logarithmic profile under neutral conditions and the general shape of the similarity functions under non-neutral conditions have been confirmed in observations on the water surface [3,16]. From these measure-

Aerodynamic Roughness Length
The variation in wind speed with height through a logarithmic profile under neutral conditions and the general shape of the similarity functions under non-neutral conditions have been confirmed in observations on the water surface [3,16]. From these measurements, the transfer coefficients and the aerodynamic roughness length can be determined, z 0 , which is defined as the height above a surface where the wind speed cancels out.
At sea, the continuous movement of the surface determines a variable roughness length, which influences the logarithmic wind profile of the similarity theory. Thus, with very light winds, less than 4 m/s, the sea surface approaches a smooth aerodynamic surface, independent of the geometry of the rough element [17] and is given by the equation: where υ = 1.461 · 10 −5 m 2 /s is the kinematic molecular viscosity, as parameter of air at sea level for a standard atmosphere [3]. With wind speeds between 4 and 50 m/s, Equation (10) is applicable, called Charnock's relationship [18], where the roughness length increases when the kinematic stress increases.
where α c , it is a constant, called the Charnock parameter, which is usually quoted between 0.01 and 0.04 [19], and which takes from this range low values for open sea and high values for sites near shore. Thus, due to said behavior and the data available at the marine sites, a simplification in the logarithmic profile of the wind speed in the MABL, considering a constant roughness length and a neutral atmospheric stratification in all marine sites, is not an approach in line with the experimental evidence [20].

Literature Review, Novelty, and Objectives of This Paper
The accumulated experience in the onshore wind industry was collected by the International Electrotechnical Commission (IEC) in the IEC 61400-1 standard [21] and extended for marine conditions in the IEC 61400-3 standard [7]. Likewise, certifying entities such as Germanischer Lloyd (GL) and Det Norske Veritas (DNV) applied their codes and standards of the offshore oil and gas industries to the offshore wind industry. Thus, in relation to the wind potential, IEC [7] and GL [22] apply the power law (1) in the model of the normal wind profile (NWP), where a neutral atmospheric stratification and a constant roughness of the sea surface are considered. On the other hand, DNV [23] indicates that the variation in the wind speed with height under normal conditions adopts a logarithmic profile given by Equation (2) with neutral stratification and with a variable value of the roughness length given by the Charnock equation [18]; it also expresses the use of the power law (1) as a possible alternative.
Based on these considerations, the use of the power law (1) and logarithmic Equation (2) do not fully describe the wind potential in a marine location at any latitude, so it can give deviations in the energy calculation for deciding on the energy viability of an offshore wind project. In this sense, Lange et al. [24] indicate that the estimates of the power output of wind turbines improve significantly if atmospheric stability is taken into account. Kalvig et al. [6] also indicate that there is a gap between the best knowledge of MABL and its practical applications, requiring a review of the influence on energy production in offshore wind turbines if the stability of the atmosphere is taken into account along with the interaction of the wind on the surface of the sea.
Therefore, the novelty of the present study of wind potential considers the use of the extended logarithmic Equation (4) for all cases of atmospheric stability, neutral and non-neutral, which occur in the MABL, with a continuous variation in ten-minute intervals. The logarithmic factor of this equation includes the aerodynamic roughness length, which is related to a continuous variation of the swell, which has a minimal effect on the wind speed profile in the open sea [25].
The objectives of the paper are, on the one hand, the study of the wind potential by means of the M-O extended equation based on data from a ten-minute series, and, on the other, its application to the calculation of the energy of a commercial marine wind turbine.
The characterization of the marine site in each year of the study includes: -A quantitative analysis of the atmospheric stratification, to ascertain the boundary conditions and assess their influence on the wind speed. - The calculation of the friction velocity and the aerodynamic roughness length as a function of the wind speed measured at the lower level. The deviations in the model with respect to the real data were obtained by means of the relation of the measured and predicted wind speeds to the height of the hub of a commercial wind turbine.
Finally, as an application of the proposed model, the validation of the results is summarized in the estimation of the energy that said wind turbine would produce in the open sea, through its operating parameters.

Materials and Methods
The study developed in this article is outlined in Figure 2. The input data were obtained from the German research platform FINO 3 [26] in a ten-minute series during the years 2016, 2017, and 2018, these are the temperature of the sea surface and meteorological data at 30 m above sea level, and the characteristic curves of the Siemens-Gamesa SWT-6.0 MW-154 wind turbine with the hub at 100 m. For data processing, several computer tools were developed in order to model the variation in the wind speed at the height of the hub of the wind turbine.
other, its application to the calculation of the energy of a commercial marine wind turbine.
The characterization of the marine site in each year of the study includes: -A quantitative analysis of the atmospheric stratification, to ascertain the boundary conditions and assess their influence on the wind speed. - The calculation of the friction velocity and the aerodynamic roughness length as a function of the wind speed measured at the lower level.
The deviations in the model with respect to the real data were obtained by means of the relation of the measured and predicted wind speeds to the height of the hub of a commercial wind turbine.
Finally, as an application of the proposed model, the validation of the results is summarized in the estimation of the energy that said wind turbine would produce in the open sea, through its operating parameters.

Materials and Methods
The study developed in this article is outlined in Figure 2. The input data were obtained from the German research platform FINO 3 [26] in a ten-minute series during the years 2016, 2017, and 2018, these are the temperature of the sea surface and meteorological data at 30 m above sea level, and the characteristic curves of the Siemens-Gamesa SWT-6.0 MW-154 wind turbine with the hub at 100 m. For data processing, several computer tools were developed in order to model the variation in the wind speed at the height of the hub of the wind turbine.

Initial Parameters
The dimensionless Richardson number, Ri , represents the balance between the thermal and mechanical effects in the generation of the turbulent state of the lower layers of the atmosphere. In particular, the atmospheric stability conditions were determined by the bulk Richardson number, b Ri , given by the Equation (11):

Initial Parameters
The dimensionless Richardson number, Ri, represents the balance between the thermal and mechanical effects in the generation of the turbulent state of the lower layers of the atmosphere. In particular, the atmospheric stability conditions were determined by the bulk Richardson number, Ri b , given by the Equation (11): where g is the acceleration due to gravity and U is the wind speed measured at a height z.
The choice of this method is due to the fact that it requires few measuring instruments to characterize the atmospheric stability in the high seas, with good results [27,28]. The grouping of data obtained by the bulk Richardson number has been performed in a range between −0.02 and +0.02 for neutral atmosphere, for values less than −0.02 unstable atmosphere, and for values greater than +0.02 atmosphere stable [29].
The calculation of the friction velocity, u * , was carried out by means of the numerical method of iteration of the fixed point, as follows: If U < 4 m/s, Equation (9) is substituted in Equation (4): If 4 ≤ U ≤ 50 m/s, Equation (10) is substituted in Equation (4): with the Charnock constant for open sea, α c = 0.011. Likewise, in Equations (12) and (13) the stability function Ψ M (z/L) is zero for a neutral atmosphere, or it takes the value given by Equation (7) for an unstable atmosphere or the value given by Equation (8) for a stable atmosphere. Known the friction velocity, the aerodynamic roughness length is obtained by Equations (9) and (10). The Obukhov length, L, is calculated below in each ten-minute interval using Equation (5), where the average virtual potential temperature, θ v , for unsaturated air, with mixing ratio r, is given by the following equation [3]: with r = 20 g/kg, and the potential temperature, θ, by the following expression: where P a is the air pressure and P o the reference pressure, which will be taken as 100 kPa.

Wind Speed at Hub Height
In the last stage of modeling, known u * , z 0 and L, the wind speed is determined, U h , at the new height, h, of the surface boundary layer, according to the general expression (4), for the following cases: -Neutral atmospheric stratification, h/L = 0, where Equation (4) reduces to Equation (2). -Unstable atmospheric stratification, h/L < 0, substituting Equation (7) in Equation (4), it results: where  (8) in Equation (4), we obtain:

Operating Parameters
The action of the wind on the rotor of a wind turbine during an annual period produces the available energy, E d,annual , expressed by discrete sums by the equation: being T i the time of occurrence of each interval of the wind speed and • W i the power of the incident wind on the turbine rotor.
Likewise, of the energy available in an annual period, only the energy used is effective, E a,annual , which is given by the equation: where C P,i is the rotor power coefficient of the wind turbine for each wind speed interval. A part of the energy used by the rotor is converted into drive energy, which transmits the motor torque to the electrical machine.
In the end, the conversion of the energy in the electric generator produces the useful energy, E u,annual , which is given by discrete sums by the equation: where P i is the value of the power, extracted from the power curve of the selected wind turbine for each incident wind speed interval in a year. The summary of the energy calculations of a wind turbine is determined by the following operating parameters [30]: -Capacity Factor (CF) is the ratio between the useful energy, E u,annual , produced by a wind turbine in a period of a year, T, and the energy that would have been produced if during this period it had been operating continuously at its rated power, P n , which can be expressed by the following equation: -Usage Time (UT) is the number of hours that a wind turbine would have to operate at full load to produce the same energy in a year that would be obtained during its actual operation. It is given by the equation: We will use these parameters to carry out a comparative study of the energy that the selected wind turbine would produce at the site where the data are measured. Finally, the results obtained with the data measured on the FINO 3 research platform at the same elevation as the wind turbine are validated, that is, at 100 m above sea level.

Atmospheric Stability
Since the wind energy depends on the type of atmospheric stratification in a site, the modeling of the wind potential in the open sea begins by determining it, through the bulk Richardson number, Ri b , given by Equation (11).
Thus, from ten-minute data, such as air and sea surface temperatures and the wind speed at a height of 30 m above sea level, the corresponding dimensionless values were obtained. The grouping of these in a band around zero, between −0.02 and +0.02, identifies the quasi-neutral stratification, which we will also call neutral. For dimensionless values lower than −0.02, it will generically be called unstable atmospheric stratification. For dimensionless values higher than +0.02 it will be called, in general, stable atmospheric stratification.
The results obtained according to the bulk Richardson number are grouped in Table 1 with the objective of knowing the atmospheric stratification in each year of study at the site of the FINO 3 platform. The efficacy data is the relationship between valid data and total data in a year. Most of the cases represent quasi-neutral and unstable atmospheric stratifications, with the lowest number of cases being stable stratification. Even though the valid data measured in each year are different, it can be observed that the types of atmospheric stability remain within a certain range of percentage values, which characterizes, together with the direction and speed of the wind, the wind potential of the marine site. The ten-minute values obtained are represented in Figures 4a, 5a and 6a as dispersion graphs of the bulk Richardson number, Ri b , as a function of the wind speed measured at 30 m in the years 2016, 2017, and 2018. In a narrow band between −0.02 and 0.02 there are cases of neutral stratification for the entire range of wind speeds up to 25 m/s. From the horizontal by 0.02 upwards there are all the stable cases that, depending on the proximity to said line, can be: slightly stable, stable, and very stable, and from the horizontal by −0.02 downwards all the unstable cases, here, also depending on the proximity to said line, they can be: slightly unstable, unstable, and very unstable [31,32]. are all the stable cases that, depending on the proximity to said line, can be: slightly stable, stable, and very stable, and from the horizontal by −0.02 downwards all the unstable cases, here, also depending on the proximity to said line, they can be: slightly unstable, unstable, and very unstable [31,32].
The reason for this behavior is due to the balance between the thermal effects (buoyancy force) and the mechanical effects (wind shear force) expressed in Equation (11), where the absolute value of depends mainly on the wind speed and its sign of the temperature difference [29]. If buoyancy is dominant, positive (stable atmosphere) or negative (unstable atmosphere), the Richardson bulk number is above 0.02 or below −0.02, respectively. On the other hand, if the shear force of the wind is dominant, the value of the Richardson bulk number is in the range between 0.02 and −0.02.
Likewise, Figures 4b, 5b and 6b represent the frequencies of the values obtained according to the bulk Richardson number in the years 2016, 2017, and 2018, respectively. The quasi-neutral data, drawn with a filled bar, add up to a total of 46.31%, 51.78%, and 43.37% in each of the related years, the rest of the data being drawn with unfilled bars, slightly unstable towards the negative part, and slightly stable towards the positive part.      Figure 6. Graphs relating to atmospheric stability in 2018. Table 2 shows the distribution of the atmospheric stratification data at the FINO 3 site in the operating range, R, of the selected wind turbine. Thus, discarding the values for wind speeds below 3 m/s, neutral stratification is predominant and intensifies with respect to the total data, it is closely followed by unstable stratification and stable stratification is very far behind.

Friction Velocity
Friction velocity, u  , scatter diagrams are shown in Figures 7-9, as a function of wind speed measured at 30 m above sea level. For light winds, less than 4 m/s, Equation (12) was applied and for winds between 4 and 25 m/s, Equation (13) was used. Typical friction velocity values range from 0 m/s for calm to 1 m/s for strong winds [3,16]. As can be seen in the first part of the figures, in the region of light winds, the friction velocity is less than 0.2 m/s and the atmospheric stratification zones are grouped. In the following region, the force of the wind on the sea surface, which originates from the waves, is translated into the friction velocity, which does not exceed 1 m/s for 25 m/s; in addition, the three regions are differentiated according to the atmospheric stratification, the values for neutral stratification appear highly concentrated and are given for all surface wind speeds over the sea. The reason for this behavior is due to the balance between the thermal effects (buoyancy force) and the mechanical effects (wind shear force) expressed in Equation (11), where the absolute value of depends mainly on the wind speed and its sign of the temperature difference [29]. If buoyancy is dominant, positive (stable atmosphere) or negative (unstable atmosphere), the Richardson bulk number is above 0.02 or below −0.02, respectively. On the other hand, if the shear force of the wind is dominant, the value of the Richardson bulk number is in the range between 0.02 and −0.02.
Likewise, Figures 4b, 5b and 6b represent the frequencies of the values obtained according to the bulk Richardson number in the years 2016, 2017, and 2018, respectively. The quasi-neutral data, drawn with a filled bar, add up to a total of 46.31%, 51.78%, and 43.37% in each of the related years, the rest of the data being drawn with unfilled bars, slightly unstable towards the negative part, and slightly stable towards the positive part. Table 2 shows the distribution of the atmospheric stratification data at the FINO 3 site in the operating range, R, of the selected wind turbine. Thus, discarding the values for wind speeds below 3 m/s, neutral stratification is predominant and intensifies with respect to the total data, it is closely followed by unstable stratification and stable stratification is very far behind.

Friction Velocity
Friction velocity, u * , scatter diagrams are shown in Figures 7-9, as a function of wind speed measured at 30 m above sea level. For light winds, less than 4 m/s, Equation (12) was applied and for winds between 4 and 25 m/s, Equation (13) was used. Typical friction velocity values range from 0 m/s for calm to 1 m/s for strong winds [3,16].   In general, the friction speed increases when the wind speed increases but at a different rate depending on the behavior of the sea surface in the limit of 4 m/s, as a smooth or rough surface, and atmospheric stratification. Thus, for wind speeds lower than 4 m/s, the friction speed growth slope is less than for wind speeds equal to, or greater than, 4 m/s;in addition, in this range, the function ) L z (  of Equation (4) also it has its effects,   In general, the friction speed increases when the wind speed increases but at a different rate depending on the behavior of the sea surface in the limit of 4 m/s, as a smooth or rough surface, and atmospheric stratification. Thus, for wind speeds lower than 4 m/s, the friction speed growth slope is less than for wind speeds equal to, or greater than, 4 m/s;in addition, in this range, the function (4) also it has its effects,   In general, the friction speed increases when the wind speed increases but at a different rate depending on the behavior of the sea surface in the limit of 4 m/s, as a smooth or rough surface, and atmospheric stratification. Thus, for wind speeds lower than 4 m/s, the friction speed growth slope is less than for wind speeds equal to, or greater than, 4 As can be seen in the first part of the figures, in the region of light winds, the friction velocity is less than 0.2 m/s and the atmospheric stratification zones are grouped. In the following region, the force of the wind on the sea surface, which originates from the waves, is translated into the friction velocity, which does not exceed 1 m/s for 25 m/s; in addition, the three regions are differentiated according to the atmospheric stratification, the values for neutral stratification appear highly concentrated and are given for all surface wind speeds over the sea.
In general, the friction speed increases when the wind speed increases but at a different rate depending on the behavior of the sea surface in the limit of 4 m/s, as a smooth or rough surface, and atmospheric stratification. Thus, for wind speeds lower than 4 m/s, the friction speed growth slope is less than for wind speeds equal to, or greater than, 4 m/s;in addition, in this range, the function Ψ M (z/L) of Equation (4) also it has its effects, distinguishing the three stratification zones.  (9) and (10), as a function of the wind speed measured at 30 m above sea level, distributed according to atmospheric stability. Here too, for winds of 25 m/s, the values of the roughness length of the sea surface are very low, that is, they slightly exceed 1 mm.

Aerodynamic Roughness Length
Energies 2022, 15, x FOR PEER REVIEW 12 of 18

Aerodynamic Roughness Length
Figures 10-12 show the scatter diagrams of the roughness length, by application of Equations (9) and (10), as a function of the wind speed measured at 30 m above sea level, distributed according to atmospheric stability. Here too, for winds of 25 m/s, the values of the roughness length of the sea surface are very low, that is, they slightly exceed 1 mm.    (9) and (10), as a function of the wind speed measured at 30 m above sea level, distributed according to atmospheric stability. Here too, for winds of 25 m/s, the values of the roughness length of the sea surface are very low, that is, they slightly exceed 1 mm.      In the graphs of Figures 10-12, the same trend is observed in the limit of 4 m/s; thus, for light winds, the sea surface behaves as a smooth surface, indicated by Equation (9), with a roughness length value close to zero and two regions atmospheric stratification: neutralstable and unstable. On the other hand, for winds between 4 and 25 m/s, Equation (10) has been applied and three clearly differentiated regions are represented, with a neutral atmosphere, the data obtained are concentrated in a very narrow band; on the other hand, for unstable and stable conditions the regions increase, being more compact with unstable stratification and more dispersed with stable stratification. In addition, in these figures, it can be observed that when the type of stability is neutral, the entire range of wind speeds occurs, that is, up to 25 m/s and higher than this. Instead, when the atmospheric stratification is unstable, the wind speed does not exceed 20 m/s, its frequent values being up to 15 m/s, and when it is stable, it does not exceed 15 m/s, its most frequent values being up to 10 m/s.

Ratio of Wind Speeds
Knowing the wind speeds at heights of 30 m and 100 m above sea level on the FINO 3 platform, Figures 13-15 represent the relationships of the averages of the wind speeds and their standard errors, measured (m) and predicted (p) at 100 m, between 3 and 25 m/s, where the modeling of the wind potential in the open sea is justified and compared. The data have been grouped for different types of atmospheric stability, in Figures 13a, 14a and 15a, neutral, and in Figures 13b, 14b and 15b, unstable and stable. In the years of study, for neutral stratification, the largest deviations occur between 3 and 9 m/s of the measured wind speed, for unstable stratification, the deviations are less than 5%, and for stable stratification, the largest deviations can be punctually up to the 15% for different measured wind speeds. In the graphs of Figures 10-12, the same trend is observed in the limit of 4 m/s; thus, for light winds, the sea surface behaves as a smooth surface, indicated by Equation (9), with a roughness length value close to zero and two regions atmospheric stratification: neutral-stable and unstable. On the other hand, for winds between 4 and 25 m/s, Equation (10) has been applied and three clearly differentiated regions are represented, with a neutral atmosphere, the data obtained are concentrated in a very narrow band; on the other hand, for unstable and stable conditions the regions increase, being more compact with unstable stratification and more dispersed with stable stratification. In addition, in these figures, it can be observed that when the type of stability is neutral, the entire range of wind speeds occurs, that is, up to 25 m/s and higher than this. Instead, when the atmospheric stratification is unstable, the wind speed does not exceed 20 m/s, its frequent values being up to 15 m/s, and when it is stable, it does not exceed 15 m/s, its most frequent values being up to 10 m/s.  Figures 13b, 14b and 15b, unstable and stable. In the years of study, for neutral stratification, the largest deviations occur between 3 and 9 m/s of the measured wind speed, for unstable stratification, the deviations are less than 5%, and for stable stratification, the largest deviations can be punctually up to the 15% for different measured wind speeds.

Ratio of Wind Speeds
(a) (b) Figure 13. Bin-averaged ratio of measured and predicted to 100 m in 2016. Figure 13. Bin-averaged ratio of measured and predicted to 100 m in 2016.  Thus, with neutral stratification, in homogeneous and stationary wind conditions, Equation (2) of the M-O similarity theory adequately models the variation in wind speed with height. The largest deviations between the predicted value and the measured values occur for low wind speeds. On the other hand, for unstable stratification conditions, in Equation (4), with the Paulson stability function, the predicted data of wind speeds are similar to the measured data. This type of stratification is frequent in this location in an annual period. While for stable stratification conditions, Equation (4), with the Businger-Dyer stability function, presents the largest deviations between the model and the measured values; however, this type of stratification is rare at the site of study.

Energy Production
The known wind speeds, measured and modeled at 100 m, to calculate the energy production, the power curves, and the power coefficient of the Siemens-Gamesa SWT-6.0 MW-154 marine wind turbine were selected, as represented in Figure 16 for an air density of 1.225 kg/m 3 . In this wind turbine, with the hub at 100 m above sea level, the starting speed is 3 m/s, the nominal speed is 14 m/s, and the cutting speed is 25 m/s. To approximate the integration of the energy calculation to discrete sums, the wind Thus, with neutral stratification, in homogeneous and stationary wind conditions, Equation (2) of the M-O similarity theory adequately models the variation in wind speed with height. The largest deviations between the predicted value and the measured values occur for low wind speeds. On the other hand, for unstable stratification conditions, in Equation (4), with the Paulson stability function, the predicted data of wind speeds are similar to the measured data. This type of stratification is frequent in this location in an annual period. While for stable stratification conditions, Equation (4), with the Businger-Dyer stability function, presents the largest deviations between the model and the measured values; however, this type of stratification is rare at the site of study.

Energy Production
The known wind speeds, measured and modeled at 100 m, to calculate the energy production, the power curves, and the power coefficient of the Siemens-Gamesa SWT-6.0 MW-154 marine wind turbine were selected, as represented in Figure 16 for an air density of 1.225 kg/m 3 . In this wind turbine, with the hub at 100 m above sea level, the starting speed is 3 m/s, the nominal speed is 14 m/s, and the cutting speed is 25 m/s.
To approximate the integration of the energy calculation to discrete sums, the wind speeds between two integers of the characteristic curves were discretized into twenty values. Table 3 shows the energy production of said wind turbine in relative terms according to the efficiency of the data cited in each year, which indicate that the real energy that a wind turbine would produce in said location would be higher.

Energy Production
The known wind speeds, measured and modeled at 100 m, to calculate the energy production, the power curves, and the power coefficient of the Siemens-Gamesa SWT-6.0 MW-154 marine wind turbine were selected, as represented in Figure 16 for an air density of 1.225 kg/m 3 . In this wind turbine, with the hub at 100 m above sea level, the starting speed is 3 m/s, the nominal speed is 14 m/s, and the cutting speed is 25 m/s. To approximate the integration of the energy calculation to discrete sums, the wind speeds between two integers of the characteristic curves were discretized into twenty values. Table 3 shows the energy production of said wind turbine in relative terms according to the efficiency of the data cited in each year, which indicate that the real energy that a wind turbine would produce in said location would be higher.
The wind speed at 100 m predicted by the method produces the input energy, or available, Ed,annual, in the wind turbine rotor, given by Equation (18). After passing through the rotor, the wind speed slows down and the rotor output energy, or harnessed energy, Ea,annual, is given by Equation (19). The last two columns show the useful energy, Eu,annual, produced by the wind turbine by applying Equation (20),with the values of the wind  The wind speed at 100 m predicted by the method produces the input energy, or available, E d,annual , in the wind turbine rotor, given by Equation (18). After passing through the rotor, the wind speed slows down and the rotor output energy, or harnessed energy, E a,annual , is given by Equation (19). The last two columns show the useful energy, E u,annual , produced by the wind turbine by applying Equation (20),with the values of the wind speed at 100 m predicted by the method and measured at the site. The results obtained show that the useful energies for each year are similar, verifying the validity of the proposed method. Table 4 shows the equivalent operating parameters, CF and UT, of said wind turbine, for the wind speeds at the hub height measured on the FINO 3 platform and obtained by the method, in the years of the study, where the relative deviations do not exceed 2.5%.

Discussion
Due to the difficulties in measuring the height of meteorological data in marine conditions, a model was proposed to determine the variation in wind speed with height from meteorological data measured at levels close to the sea surface, for its application to the energy production of the wind turbines to be installed in new offshore wind farms.
The information obtained from the German FINO 3 research platform was of great importance to carry out our study, since it has meteorological data at different heights and oceanographic data, with high efficiency, which allows comparing wind speeds measured and predicted by the method at the height of the hub of a commercial wind turbine, which in our case was 100 m above sea level.
The modeling of the exposed wind potential is based on the equations of the quantity of motion of the M-O similarity theory for neutral conditions, Equation (2), and non-neutral conditions, Equation (4), being said conditions determined by the bulk Richardson number. The greatest difficulty arose in the calculation of the friction speed as it was present both in the equations of the aerodynamic roughness length and in the Obukhov length, being solved by means of the numerical calculation of iteration of the fixed point.

Conclusions
Wind speed depends on the type of atmospheric stratification, Figures 5-10. In this way, when the stratification is neutral, all wind speeds were measured up to 25 m/s and higher, when it is unstable the most frequent values were up to 15 m/s, and when it is stable up to 10 m/s. The roughness length of the sea surface is a characteristic parameter of Equations (2) and (4), of a variable nature, and with reduced values; thus, for the most extreme conditions of the measured wind speeds, values were obtained slightly greater than 1 mm.
It is usually considered that in open sea conditions the atmospheric stratification is mainly neutral [33], but in the results presented it can be observed that the unstable stratification follows it very closely, so we can infer that the atmospheric stratification in an annual period can be a characteristic of the site that influences the wind potential.
Wind speeds at 100 m above sea level, measured on the platform and obtained by the method, were applied to the calculation of the energy of the SWT-6.0 MW-154 wind turbine presenting acceptable deviations, so they are considered satisfactory for an industrial application.
The proposed modeling can be extended to study the energy viability of new offshore wind farms based on in situ measurements during the research phase of the project, where the heights of the meteorological measurements are lower than the height of the hub of the wind turbine that will occupy the site.  Acknowledgments: To the Bundesministerium für Wirtschaft und Energie (BMWi) and Projektträger Jülich (PTJ), both of the Government of Germany, for facilitating the use of data of the FINO Project. To Siemens Gamesa Renewable Energy, for the use of the characteristic curves of the SWT-6.0-154 wind turbine.

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