Relationship between Sea Surface Drag Coefﬁcient and Wave State

: The sea surface drag coefﬁcient plays an important role in momentum transmission between the atmosphere and the ocean, which is affected by ocean waves. The total air–sea momentum ﬂux consists of effective momentum ﬂux and sea spray momentum ﬂux. Sea spray momentum ﬂux involves sea surface drag, which is largely affected by the ocean wave state. Under strong winds, the sea surface drag coefﬁcient ( C D ) does not increase linearly with the increasing wind speed, namely, the increase of C D is inhibited by strong winds. In this study, a sea surface drag coefﬁcient is constructed that can be applied to the calculation of the air–sea momentum ﬂux under high wind speed. The sea surface drag coefﬁcient also considers the inﬂuence of wave state and sea spray droplets generated by wave breaking. Specially, the wave-dependent sea spray generation function is employed to calculate sea spray momentum ﬂux. This facilitates the analysis not only on the sensitivity of the sea spray momentum ﬂux to wave age, but also on the effect of wave state on the effective C D ( C D , ef f ) under strong winds. Our results indicate that wave age plays an important role in determining C D . When the wave age is >0.4, C D decreases with the wave age. However, when the wave age is ≤ 0.4, C D increases with the wave age at low and moderate wind speeds but tends to decrease with the wave age at high wind speeds.


Introduction
Sea surface wind stress is the main factor that drives upper-ocean circulation and wind waves. An accurate estimate of sea surface wind stress is important for modelling and predicting the dynamic processes involving ocean and atmosphere interactions. Therefore, sea surface wind stress has been widely observed and studied in the past decades, which is usually parameterized by a sea surface drag coefficient (C D ) or sea surface aerodynamic roughness [1][2][3][4]. The study of sea surface wind stress is thus the study of C D or sea surface aerodynamic roughness.
The initial atmospheric boundary layer theory assumed that the interface was the fixed land. Therefore, under the condition of neutral stable stratification, air-sea flux only depended on the variation of wind speed, and the influence of wind waves was not taken into account. Researchers tended to parameterize C D only by wind speed [5][6][7][8]. However, these relationships between C D and wind speed obtained by different researchers differed greatly.
Ocean waves often exist on the sea surface, and their simulations and effects have been studied extensively [9,10]. In recent years, more and more evidence has indicated that the influence of wave state on wind stress cannot be ignored when investigating air-sea interaction [4,[11][12][13]. However, the results of laboratory and field measurements were not always consistent, and the relationships between C D and the wave age obtained by the two approaches were even contradictory [11,14]. considered the modified action of sea spray momentum flux on the total momentum flux at the air-sea interface, and the sea spray momentum flux representing the foam action was reflected in the high wind speed drag coefficient relation. In the expression of the high wind velocity drag coefficient of Makin [26], he only parametrized the drag coefficient using measurements from Powell [21]. Therefore, the high wind speed drag coefficient relation of Andreas [25] can clearly describe the decay effect on the drag coefficient at high wind speed. Therefore, the theory of Andreas [25] is worth further investigation.
It is worth noting that the generation of sea spray droplets is also related to the wave state [29][30][31][32]. As found by the studies mentioned above, when the wave age increases, the number of droplets produced by breaking waves increases as well. However, these studies have not considered the direct effect of sea spray droplets on the air-sea momentum flux. Therefore, under the condition of high wind speeds, the influence of wave state on the total air-sea momentum flux should be considered, as well as the influence of sea spray momentum flux.
The structure of this paper is as follows. Firstly, we provided the expression of the sea surface drag coefficient, which is able to describe the effect on the air-sea momentum flux and includes the sea spray generation function (SSGF) that calculates the sea spray momentum flux. Secondly, we carry out sensitivity experiments, analyzing the influence of wave state on the calculation of sea spray momentum flux. Thirdly, we analyze the influence of wave state on the sea surface drag coefficient and sea surface roughness length in different wave age conditions. Fourthly, we verified the calculated sea surface drag coefficient and sea surface roughness length by using observation data.

Calculations on Sea Surface Drag Coefficient (C D )
According to Andreas and DeCosmo [27], the momentum exchange between droplets and the surrounding atmosphere mainly occurs in the evaporative layer of sea spray droplets, which approximately has a thickness of one significant wave height. In the evaporative layer of sea spray droplets, the droplets draw momentum from the air, and the total air-sea momentum flux τ T can be divided into two parts: the effective air-sea interface momentum flux τ eff and the sea spray momentum flux τ sp . τ T = τ eff + τ sp (1) τ T can generally be obtained by the following formula: where u * is the friction velocity, and ρ a is the air density. The total air-sea momentum flux can also be expressed by the classical bulk formula: where C D is the drag coefficient that does not exclude the influence of droplets, and U 10 is the wind speed at 10 m above the sea surface. The effective air-sea interface momentum flux τ eff is expressed as follows: where τ e f f is the effective C D (C D, e f f ) without considering the influence of sea spray droplets. Substituting Equations (2) and (4) into Equation (1) yields Then: Equation (6) is the expression of C D,e f f considering the influence of the droplets. In Equation (6), it shows that the effective sea surface drag coefficient is affected by the sea spray droplets. The second term on the right is the influence of the sea spray droplets effect. With the increase of the sea spray droplets effect, sea surface drag coefficient is reduced. The sea spray droplet effect is mainly reflected in the ratio of the sea spray momentum flux and sea surface wind speed square. Therefore, accurate estimation of sea spray momentum flux is needed to make good use of Equation (6), to calculate the effective sea surface drag coefficient.

Calculations on Sea Spray Momentum Flux
Sea spray momentum flux τ sp can be calculated by using Andreas's [25] scheme, where dF dr 0 is called the SSGF, and ρ w is seawater density. r lo and r hi are the minimum and maximum radii of droplets when calculating τ sp , respectively. r 0 is the initial radius of the droplet, and u sp (r 0 ) is the horizontal velocity of the droplet before it falls back to the sea surface. The sea surface wind speed can be expressed as follows: Andreas [25] suggested that all droplets be subjected to local winds for a brief flight over the sea surface before falling back. Therefore, u sp (r 0 ) does not depend on the value of the initial droplet radius, but it depends on the height of the sea spray droplet, which can be expressed as follows: where z 0 is the sea surface roughness length, κ(= 0.4) is the von Karman constant, and u sp (z s ) is the wind speed at z s above the sea surface. As suggested by Iida et al. [33] and Zhao et al. [29], z s = 0.635H s , where H s is the significant wave height, which can be parameterized with the empirical relationship [14,34]: where β is the wave age parameter, which represents the state of wave growth. Sea surface roughness length z 0 can be obtained by using the parameterization scheme of Zhao and Li [4], which depends on the wave state and wind speed:

Sea Spray Generation Function
In Equation (7), the SSGF dF dr 0 is generally regarded as a function of the initial radius r 0 and U 10 . If the scale spectrum of sea spray droplets is independent of wind speed, the SSGF can be divided into two independent parts: one is only related to wind speed, and the other is only related to droplet size, i.e., where f 1 and f 2 are the independent functions of the SSGF. Andreas [25,35] simplified the sea spray momentum flux and assumed that the sea spray momentum flux was proportional to the fourth power of the friction velocity: Studies had shown that various functions were related to wind speed among different SSGFs, and there was a considerable difference in magnitude between these functions [28]. In fact, wind is not the only factor that affects the formation of sea spray droplets at the sea surface, as both wind and waves at the surface play important roles. Waves affect the exchanges of air-sea fluxes. Researchers have pointed out that the uncertainty of the calculated value of the sea spray momentum flux was caused by ignoring the influence of the wave state, and they believed that the SSGF should consider not only the influence of wind but also the influence of the wave state [29]. Assuming that the dependence of the SSGF on sea surface wind speed and wave state is independent of the sea spray droplet radius, the SSGF can be divided into two independent functions: Here, R B = u 2 * ω p ν is the windsea Reynolds number, in which ν is the coefficient of air viscosity, and ω p is the angular frequency corresponding to the peak of the wind wave spectrum. The windsea Reynolds number is widely used in the air-sea boundary layer to study the influence of the wave state on the gas and momentum transports on the sea surface. Together with Equations (2) and (3), R B can also be written as follows: Applying this theory on the measured data in the laboratory, Zhao et al. [29] proposed a new SSGF: The calculated results of Equation (16) can well cover the range of the calculated values of the mainstream droplet generation functions [29]. However, this formula is only applicable to sea spray droplets with a radius range of 30 ≤ r 0 ≤ 500 µm, which means that the droplets induced by bubbles (what are their sizes?) are not taken into account.
Based on the coverage rate of whitecaps, Monahan et al. [36] gave the SSGF applicable to the droplets generated by bubbles and the applicable radius of sea spray droplets with a relative humidity of 80% is 0.8-10 µm: where r 80 is the droplet radius when the reference relative humidity is 80%, τ d = 3.53 s is the typical decay time of whitecaps, and W is the whitecap coverage rate. Monahan et al. [36] suggested that W be only related to wind speed. dE dr 80 is the number of droplets generated per unit area of whitecap in its decay stage, per unit time and within the increment range of the unit radius corresponding to r 80 .
Monahan et al. [36]  The generation function established in Equation (18) by Monahan et al. [36] has been shown to be useful by several subsequent studies [37][38][39]. Taking advantage of the tank observation data for whitecaps, Woolf et al. [38] modified the droplet generation function by Monahan et al. [36], which is applicable to bubbles, and pointed out that the revised The fitting relationships in Equations (18) and (19) are complicated, and the droplet scale spectrum is different from the common spume droplets' spectrum in terms of math expression. In the present study, the least square method is adopted to fit the average values of Equations (18) and (19), yielding dE dr 80 with a simpler form. Let r 80 be in the range of 1-10 µm and substitute into Equations (18) and (19) to calculate the results, then average the results on each r 80 . Based on the average value, the least square method is used to obtain the relationship between r 80 and dE dr 80 where R is the correlation coefficient. The following relationship between r 80 and r 0 by Andreas [40] is considered: Then, the SSGF applicable to droplets with an initial radius of 2 ≤ r 0 ≤ 20 µm can be expressed as follows: dF According to Equation (22), the sea spray droplet scale spectrum is simplified as a function of r −3 0 , whereas f (R B ) represents a function related to the wave state. From Zhao et al. [29], we know that the function can be expressed by the wave state parameter R B .
Andreas [41] recommended the SSGF based on the data by Wu [5]: where C 1 , C 2 , and C 3 are functions of U 10 . Based on the calculated value at U 10 = 15 m/s of the SSGF SEA93 proposed by Smith et al. [42], Andreas [41] determined that C 1 = 1.955 × 10 3 , C 2 = 1.331 × 10 6 , and C 3 = 3.344 × 10 16 . Therefore, based on the sea spray droplet scale spectrum given by Equation (23), the formula by Zhao et al. [29] is revised here. The minimum applicable sea spray droplet radius is extended to 20 µm, giving: Considering the continuity of the SSGF, and using Equations (22) and (24), the following relation is obtained: According to Equation (25), dF/dr 0 is an expression composed of four sections, and its value is proportional to R B and inversely proportional to r 0 . The variation trend of spray droplet number with r 0 in the first section is different from that in the other three sections. In the first section, dF/dr 0 is proportional to r −3 0 , and the other three sections are mainly in a similar form to Wu [5], that is, dF/dr 0 is proportional to r −1 0 , r −3 0 , and r −8 0 as the range of radius increases. This is because the first section describes sea spray droplets produced by bubbles, and its generation mechanism is different from that described in the other three sections.
With Equations (7), (9)- (11), and (25), the momentum flux of sea spray droplets with a radius from 2 to 500 µm can be calculated, which incorporates the influence of the wave state.

Wave State Influence on Sea Spray Momentum Flux
From Equation (7), we can see that the core of sea spray momentum flux is the SSGF. Most of the existing sea spray generation functions do not consider the influence of waves. As shown in Figure 1, the sea spray generation function in Equation (25) has different values under different wave ages, which confirms that it is necessary to consider the effect of the wave state in the sea spray generation function. Note that the sea spray generation function increases with the increase of wave age, instead of showing a varied trend as a function of wave age, as indicated in Equation (11). This is due to the fact that the significant wave height in Equation (10) increases with the increase of wave age. Therefore, the sea spray generation function maintains an increasing trend when wave age increases, though with a faster growth when the wave age is smaller than 0.4 and a slower increase when wave age is larger than 0.4. According to Equation (25), / is an expression composed of four sections, and its value is proportional to and inversely proportional to . The variation trend of spray droplet number with in the first section is different from that in the other three sections. In the first section, / is proportional to , and the other three sections are mainly in a similar form to Wu [5], that is, / is proportional to , , and as the range of radius increases. This is because the first section describes sea spray droplets produced by bubbles, and its generation mechanism is different from that described in the other three sections.
With Equations (7), (9)- (11), and (25), the momentum flux of sea spray droplets with a radius from 2 to 500 µm can be calculated, which incorporates the influence of the wave state.

Wave State Influence on Sea Spray Momentum Flux
From Equation (7), we can see that the core of sea spray momentum flux is the SSGF. Most of the existing sea spray generation functions do not consider the influence of waves. As shown in Figure 1, the sea spray generation function in Equation (25) has different values under different wave ages, which confirms that it is necessary to consider the effect of the wave state in the sea spray generation function. Note that the sea spray generation function increases with the increase of wave age, instead of showing a varied trend as a function of wave age, as indicated in Equation (11). This is due to the fact that the significant wave height in Equation (10) increases with the increase of wave age. Therefore, the sea spray generation function maintains an increasing trend when wave age increases, though with a faster growth when the wave age is smaller than 0.4 and a slower increase when wave age is larger than 0.4. The droplets accelerated by wind extract momentum from the wind near the sea surface and slow down the wind speed. Figure 2 shows the variation of the ratio of sea spray momentum flux to the total momentum flux, with wind speed at different wave ages. We can see that when the wind speed is fixed, with the growth of the wave, the ratio increases, which further indicates that the wave state has a great impact on the generation of the sea spray momentum flux. As shown in Figure 2, the sea spray momentum flux can exceed the order of 10 1 in different wave growth states, which increases with the growth of waves when the wind speed is constant. The variation of the total momentum flux with wind The droplets accelerated by wind extract momentum from the wind near the sea surface and slow down the wind speed. Figure 2 shows the variation of the ratio of sea spray momentum flux to the total momentum flux, with wind speed at different wave ages. We can see that when the wind speed is fixed, with the growth of the wave, the ratio increases, which further indicates that the wave state has a great impact on the generation of the sea spray momentum flux. As shown in Figure 2, the sea spray momentum flux can exceed the order of 10 1 in different wave growth states, which increases with the growth of waves when the wind speed is constant. The variation of the total momentum flux with wind speed is displayed in Figure 2. When the wave age is varied, the variation of the total momentum flux is not as obvious as that of the sea spray momentum flux. When the wind speed reaches approximately 35-80 m/s, the sea spray momentum flux can reach the order of the total momentum flux under different wave age conditions. Therefore, it is known that the wave state has a great influence on the order of the sea spray momentum flux. The larger the wave age is, the earlier the momentum flux of the sea spray droplet reaches the same order as the total momentum flux. speed is displayed in Figure 2. When the wave age is varied, the variation of the total momentum flux is not as obvious as that of the sea spray momentum flux. When the wind speed reaches approximately 35-80 m/s, the sea spray momentum flux can reach the order of the total momentum flux under different wave age conditions. Therefore, it is known that the wave state has a great influence on the order of the sea spray momentum flux. The larger the wave age is, the earlier the momentum flux of the sea spray droplet reaches the same order as the total momentum flux. Since the sea spray momentum flux is related to the wave state, the wave state parameters can be selected for parameterization of the wave state [29]. is the wave breaking control parameter of the air-sea interface, which can express the joint impact of waves and wind on air-sea flux. Figure 3 shows the variation of the sea spray momentum flux as a function of under different wave ages. Different from Figure 2, the change of the sea spray droplet momentum flux under different wave states is small here. This means that, compared to the parameterization of the sea spray droplet momentum flux only by wind speed, can better fit the variation of the sea spray droplet momentum flux under various wave ages. Since the sea spray momentum flux is related to the wave state, the wave state parameters R B can be selected for parameterization of the wave state [29]. R B is the wave breaking control parameter of the air-sea interface, which can express the joint impact of waves and wind on air-sea flux. Figure Figure 4 also shows the variation of the droplet momentum flux to the total momentum flux ratio under different wave ages. When the value of reaches 2 × 10 6 , the ratio is close to one, which indicates that the two fluxes are of the same order of magnitude. When the value of is fixed, the flux ratio increases with the growth of the waves.  Figure 4 also shows the variation of the droplet momentum flux to the total momentum flux ratio under different wave ages. When the value of R B reaches 2 × 10 6 , the ratio is close to one, which indicates that the two fluxes are of the same order of magnitude. When the value of R B is fixed, the flux ratio increases with the growth of the waves.   The above analysis shows that the windsea Reynolds number can fit the sea spray droplet momentum flux well under various wave conditions. After transformation, the relationship between , friction velocity, and wave age becomes: When the friction velocity and wave age are fixed, the value of affects as well. The above analysis shows that the windsea Reynolds number R B can fit the sea spray droplet momentum flux well under various wave conditions. After transformation, the relationship between R B , friction velocity, and wave age becomes: When the friction velocity and wave age are fixed, the value of C D affects R B as well. Andreas [25,35] simplified the sea spray droplet momentum flux by considering that the sea surface droplet stress to be related to the fourth power of the friction velocity, i.e., τ sp = 6.2 × 10 −5 ρ w u 4 * The results from the formula above are compared with the results of Equation (25) proposed in this study (as shown in Figure 5). The average of R B is taken for estimation, and the results are shown in Figure 5. We can see that the sea spray droplets calculated by Equation (25) can cover those by Equation (27) under different wave ages. This not only demonstrates the correctness of Equation (25), but also shows that the wave state has an impact on the sea spray droplet momentum flux, which should be considered when calculating the droplet momentum flux.

Sea Surface Drag Coefficient
The relationship between sea surface roughness z 0 and C D can be obtained by combining Equations (2), (3), and (8): Therefore, C D is related to wave age β and surface wind speed at 10 m (U 10 ). Figure 6 displays the variation of C D with wind speed under different wave ages without considering the influence of sea spray droplets. When the wave age is less than 0.4, C D increases with the increase of wave age. When the wave age is greater than 0.4, C D decreases with the increase of wave age. The drag coefficient variation with wave age is identical to that of z 0 . When wind speed is less than 30 m/s, the calculated value of C D covers the result of Wu [5], which indicates that Equation (28) is correct at low or medium wind speeds. When the wind speed increases, the result becomes smaller than that of Wu [5]. The C D in Wu [5] increased linearly with wind speed, which overestimated C D at high wind speeds, as recent studies have shown that C D decays at high wind speeds. To clearly show the wind speed scenarios when the attenuation of C D occurs, Equation (6)  Andreas [25,35] simplified the sea spray droplet momentum flux by considering that the sea surface droplet stress to be related to the fourth power of the friction velocity, i.e., = 6.2×10 * (27) The results from the formula above are compared with the results of Equation (25) proposed in this study (as shown in Figure 5). The average of is taken for estimation, and the results are shown in Figure 5. We can see that the sea spray droplets calculated by Equation (25) can cover those by Equation (27) under different wave ages. This not only demonstrates the correctness of Equation (25), but also shows that the wave state has an impact on the sea spray droplet momentum flux, which should be considered when calculating the droplet momentum flux.

Sea Surface Drag Coefficient
The relationship between sea surface roughness and can be obtained by combining Equations (2), (3), and (8): Therefore, is related to wave age and surface wind speed at 10 m ( ). Figure  6 displays the variation of with wind speed under different wave ages without considering the influence of sea spray droplets. When the wave age is less than 0.4, increases with the increase of wave age. When the wave age is greater than 0.4, decreases with the increase of wave age. The drag coefficient variation with wave age is identical to that of 0 z . When wind speed is less than 30 m/s, the calculated value of covers the result of Wu [5], which indicates that Equation (28) is correct at low or medium wind speeds. When the wind speed increases, the result becomes smaller than that of Wu [5]. The in Wu [5] increased linearly with wind speed, which overestimated at high wind speeds, as recent studies have shown that decays at high wind speeds. To clearly show the wind speed scenarios when the attenuation of occurs, Equation (6) is analyzed, considering the influence of droplets on under high wind speeds. Based on Equation (6), the variation of with the presence of sea spray droplets can be analyzed. In other words, now becomes the effective , namely , . Figure   7 shows that when wind speed is greater than 20 m/s, , is somewhat attenuated, but the corresponding initial attenuation wind speeds differ remarkably for different wave ages. For example, when the wave age is 0.2, , starts its attenuation at a wind speed of about 50 m/s. When the wave age is 1.2, however, the attenuation of , occurs at a wind speed of 20 m/s, which indicates that the attenuation of , is affected by the wave state. When wind speed is higher than 45 m/s, the drag coefficient tends to be larger when the wave age is less than 0.4 than when the wave age is greater than 0.4. under different wave ages can basically reproduce the high-wind-speed offshore observation Based on Equation (6), the variation of C D with the presence of sea spray droplets can be analyzed. In other words, C D now becomes the effective C D , namely C D,e f f . Figure 7 shows that when wind speed is greater than 20 m/s, C D,e f f is somewhat attenuated, but the corresponding initial attenuation wind speeds differ remarkably for different wave ages. For example, when the wave age is 0.2, C D,e f f starts its attenuation at a wind speed of about 50 m/s. When the wave age is 1.2, however, the attenuation of C D,e f f occurs at a wind speed of 20 m/s, which indicates that the attenuation of C D,e f f is affected by the wave state. When wind speed is higher than 45 m/s, the drag coefficient tends to be larger when the wave age is less than 0.4 than when the wave age is greater than 0.4. C D under different wave ages can basically reproduce the high-wind-speed offshore observation data of Powell et al. [21] and Jarosz et al. [43] and the laboratory observation data of Donelan et al. [23] and Takagaki et al. [44].
Based on Equation (6), the variation of with the presence of sea spray droplets can be analyzed. In other words, now becomes the effective , namely , . Figure   7 shows that when wind speed is greater than 20 m/s, , is somewhat attenuated, but the corresponding initial attenuation wind speeds differ remarkably for different wave ages. For example, when the wave age is 0.2, , starts its attenuation at a wind speed of about 50 m/s. When the wave age is 1.2, however, the attenuation of , occurs at a wind speed of 20 m/s, which indicates that the attenuation of , is affected by the wave state. When wind speed is higher than 45 m/s, the drag coefficient tends to be larger when the wave age is less than 0.4 than when the wave age is greater than 0.4. under different wave ages can basically reproduce the high-wind-speed offshore observation data of Powell et al. [21] and Jarosz et al. [43] and the laboratory observation data of Donelan et al. [23] and Takagaki et al. [44].  In the laboratory, the fetch length is short, so the ocean waves have a small wave age. Meanwhile, the ocean waves have the bigger wave age in the offshore. Especially, in the typhoon condition, because of the change of the wind direction, the waves are difficult to fully grow. Therefore, we use the small wave age calculated results to analyze the errors with the laboratory observation data and use the medium wave age calculation results to analyze the errors with the offshore observation data, as shown in Table 1. As can be seen from the root mean square error, the error of Equation (6) in the laboratory is smaller than that in the offshore. However, it is much better than the result of the Large and Pond (1981) scheme. This analysis indicates that the influence of sea spray droplets and wave state on the sea surface drag coefficient should be considered under high wind speed conditions.

Sea Surface Roughness
There is a one-to-one relationship between C D and sea surface roughness. Equation (11) states that sea surface roughness increases with the increase of wind speed. As shown in Figure 8, when the wave age is less than 0.4, the sea surface roughness increases with the increase of wave age. However, when the wave age is greater than 0.4, sea surface roughness decreases with the increase of wave age.

Sea Surface Roughness
There is a one-to-one relationship between and sea surface roughness. Equation (11) states that sea surface roughness increases with the increase of wind speed. As shown in Figure 8, when the wave age is less than 0.4, the sea surface roughness increases with the increase of wave age. However, when the wave age is greater than 0.4, sea surface roughness decreases with the increase of wave age. Similarly, using Equation (26), the effective sea surface roughness ( , ) corresponding to , can be obtained as follows: Similarly, using Equation (26), the effective sea surface roughness (z 0,e f f ) corresponding to C D,e f f can be obtained as follows: According to Figure 9, sea surface roughness also decays with the increase of wind speed, and the wind speed, at which the roughness decays, varies at different wave ages. Compared with the high wind speed offshore observation data by Powell et al. [21] and the laboratory observation data by Takagaki et al. [44], Equation (29)  According to Figure 9, sea surface roughness also decays with the increase of wind speed, and the wind speed, at which the roughness decays, varies at different wave ages. Compared with the high wind speed offshore observation data by Powell et al. [21] and the laboratory observation data by Takagaki et al. [44], Equation (29) can well reproduce the measured values under different wave ages, with a variation trend similar to that of the measured values.

Discussion
To further analyze the influence of wave age on in detail, the range of wave age is selected as 0.01 to 1.20. When the wave age is smaller than 0.4 and wind speed is low, increases with the increase of wave age, as shown by the blue shading of low wind speed in Figure 10. With the increase of wind speed, this relationship is reversed, that is, the drag coefficient decreases with the increase of wave age. This is because the momentum flux of sea spray droplets keeps increasing with the increase of wave age (shown in

Discussion
To further analyze the influence of wave age on C D in detail, the range of wave age is selected as 0.01 to 1.20. When the wave age is smaller than 0.4 and wind speed is low, C D increases with the increase of wave age, as shown by the blue shading of low wind speed in Figure 10. With the increase of wind speed, this relationship is reversed, that is, the drag coefficient decreases with the increase of wave age. This is because the momentum flux of sea spray droplets keeps increasing with the increase of wave age (shown in red in Figure 2), i.e., the momentum flux acquired from the air increases with the increase of wind speed at high wind speeds. When the wave age is greater than 0.4 but the wind speed is changing from low-to-medium to high, C D,e f f decreases with the increase of wave age from low to high wind speeds (as shown by the blue shading of Figure 10). In this case, the sea spray droplet momentum flux increases with the increase of wave age, which highlights the relation that C D,e f f decreases with the increase of wave age. Through the comparison of Figures 7 and 9, it can be seen that there are some differences between the offshore measurement data and the laboratory measurement data, in terms of the variation trend and range as a function of wind speed. When the wind speed is low to medium, , and effective sea surface roughness in the open sea are lower than those measured in the laboratory. Holthuijsen et al. [22] stated that sea surface roughness was reduced due to the effect of a storm surge.
It is worth noting that in some wave states, may decrease to zero with the increase of wind speed, as shown in Figure 10. This phenomenon also appeared in Andreas [25]. However, when becomes zero, it means that the momentum of the atmosphere will not be directly fed into the ocean. Andreas [25] reported that when sea spray droplets returned to the sea surface, the short waves of the sea surface would disappear, thus inhibiting the growth of waves. In this scenario, it is difficult for the wave to grow sufficiently. In other words, it is not easy to see the drag coefficient corresponding to larger wave ages in Figure 10 when it is approaching zero.
When wind speed is greater than 50 m/s, the measured in the laboratory is significantly greater than that measured in the open sea. As shown in Figure 7, Donelan et al. [23] stated that the air flow would separate from the wave surface and jump directly from the crest of one breaking wave to another crest, without passing through the trough. The existence of these continuous breaking waves would restrict the growth of wind waves, reducing . In the open sea, the above-stated process holds as well. In addition, according to Holthuijsen et al. [22], when wind speed is high, a large number of sea spray droplets is generated in the open sea, forming a layer of sea spray droplets, which lowers sea surface roughness. As the wave height and wavelength in the tank are smaller than those in the open sea (see Donelan et al. [23], Figure 4), the droplets would return to the sea surface soon, making it difficult to form a spray droplet layer. Therefore, the reduction of Through the comparison of Figures 7 and 9, it can be seen that there are some differences between the offshore measurement data and the laboratory measurement data, in terms of the variation trend and range as a function of wind speed. When the wind speed is low to medium, C D,e f f and effective sea surface roughness in the open sea are lower than those measured in the laboratory. Holthuijsen et al. [22] stated that sea surface roughness was reduced due to the effect of a storm surge.
It is worth noting that in some wave states, C D may decrease to zero with the increase of wind speed, as shown in Figure 10. This phenomenon also appeared in Andreas [25]. However, when C D becomes zero, it means that the momentum of the atmosphere will not be directly fed into the ocean. Andreas [25] reported that when sea spray droplets returned to the sea surface, the short waves of the sea surface would disappear, thus inhibiting the growth of waves. In this scenario, it is difficult for the wave to grow sufficiently. In other words, it is not easy to see the drag coefficient corresponding to larger wave ages in Figure 10 when it is approaching zero.
When wind speed is greater than 50 m/s, the C D measured in the laboratory is significantly greater than that measured in the open sea. As shown in Figure 7, Donelan et al. [23] stated that the air flow would separate from the wave surface and jump directly from the crest of one breaking wave to another crest, without passing through the trough. The existence of these continuous breaking waves would restrict the growth of wind waves, reducing C D . In the open sea, the above-stated process holds as well. In addition, according to Holthuijsen et al. [22], when wind speed is high, a large number of sea spray droplets is generated in the open sea, forming a layer of sea spray droplets, which lowers sea surface roughness. As the wave height and wavelength in the tank are smaller than those in the open sea (see Donelan et al. [23], Figure 4), the droplets would return to the sea surface soon, making it difficult to form a spray droplet layer. Therefore, the reduction of the drag coefficient in the tank is smaller than that in the open sea.
According to the analysis results of open-sea and laboratory measurement data in Zhao et al. [4], the wave age in the open-sea measurement data is generally greater than 0.4, while the wave age in the laboratory measurement data is typically less than 0.4. Although these were measured under low and medium wind speeds, it also shows that there are differences in wave age between the laboratory tank and the open sea. Using the C D at high wind speeds discussed in this study (Figure 7), we see that when the wind speed is higher than 50 m/s, the laboratory measurement data are mainly distributed in the region with low wave ages. Therefore, it can be inferred that the wave age in the tank is still small at high wind speeds.

Conclusions
At high wind speeds, sea spray droplets are generated by wave breaking. The sea surface drag coefficient stops increasing and decays when the wind speed reaches a certain level. Although this phenomenon has been confirmed by observations in the open sea and laboratories, there still exist some controversies about the physical mechanism of this attenuation. In this paper, a new relationship between the sea surface drag coefficient and wave state is constructed from the influence of wave state on the drag coefficient and the influence of wave state on the generation of sea spray droplets.
Firstly, we establish an SSGF (Equation (25)), which can further calculate the momentum flux of sea spray droplets with a radius from 2 to 500 µm. We show that sea spray droplets absorb the wind energy above the sea surface, which produces the sea spray momentum flux. The sea spray momentum flux increases with the increase of wind speed and wave age, which attenuate the effective momentum actually transmitted at the air-sea interface.
Secondly, based on a new spray droplet generation function (Equation (25)), sea spray momentum flux can be calculated, and a new sea surface drag coefficient under high wind speeds is established (Equation (6)), which takes into account the influence of sea spray droplets and wave state. Compared with the sea surface drag coefficient at medium wind speeds, we find that Equation (6) also performs well in these cases. According to the calculation, the effective sea surface drag coefficient at high wind speeds is attenuated, and the attenuation is different under different wave age conditions. Comparing the observation data in the open sea and from the laboratory at high wind speeds, we find that the calculated drag coefficient under different wave age conditions can basically cover the measured values. Therefore, the sea surface drag coefficient given by Equation (6) under high wind speeds is recommended for calculating air-sea momentum fluxes in the open sea and laboratory tanks.
In the future, with the increase of observation data under high wind speeds, especially for sea spray droplets, the SSGF will be more reasonable, and the calculation accuracy of sea spray momentum flux can be further improved. At the same time, the influence of ocean current on the sea surface drag coefficient can be considered in the parameterization scheme, which can further improve the theory of air-sea interaction and provide an accurate airsea momentum flux for simulating marine environments in ocean-atmosphere coupling models under high wind speeds.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The reader can ask for all the related data from the first author (shijian@nudt.edu.cn) and the corresponding author (sunyuan17a@nudt.edu.cn).