1. Introduction
The impact of the wind–wave interaction on wind stress has been widely documented in both theoretical and observational studies (e.g., Donelan et al. 1993; Rieder et al. 1996; Hara and Belcher 2002). The wave age, which describes the development state of local wind waves, has been preferentially used in the parameterization of α in several studies (e.g., Johnson et al. 1998; Jones and Toba 2001; Oost et al. 2002; Gao et al. 2006). However, different experiments have shown considerable changes in the wave-age dependence of α, and even contradictory results have been reported between field and laboratory measurements (Toba et al. 1990; Donelan 1990).
Drennan et al. (2005) performed a comprehensive assessment of the performance of the wave-age- (D03) and wave-steepness-dependent (TY01) scaling using a combination of eight datasets representing a wide range of sea states. They concluded that both scalings yield improved estimates compared to a sea-state-independent scaling under wind-wave dominant conditions, but their performance degrades significantly with an increasing contribution of the swell component.
Swell either is generated elsewhere or results from the slowing of local winds or a change in direction. In the open ocean, waves are commonly mixed with wind waves and swells. Swells are typically energetic and even dominant, significantly changing the wave age and wave steepness compared to pure wind-wave conditions (Hanley et al. 2010; Semedo et al. 2011; Li et al. 2018a), but swells do not contribute as much to the sea surface roughness. García-Nava et al. (2012b) argued in support of using only the wind sea part in the parameterization of sea surface roughness. Hwang et al. (2011a) proposed using a characteristic wave frequency defined as a weighted average of swells and wind waves to parameterize the drag coefficient in mixed waves. Smedman et al. (2003) derived from observation analysis that the drag coefficient is governed by the wave energy spectral ratio between swell and wind-wave, as well as the wave age. Carlsson et al. (2010) suggested that the impact of swell on the wind–wave interaction can explain the discrepancies among the different investigations. The impact of swell on the wind–wave interaction was discussed by Hwang et al. (2011b), who suggested that the development of wind waves tends to be enhanced in the presence of swell, and by García-Nava et al. (2009, 2012a), who suggested that swell can reduce the wave-supported stress by suppressing the intensity of short wind waves in the equilibrium range. Nevertheless, the wind–wave interaction in the presence of swell is much less understood than the interaction for pure wind waves, largely accounting for the slow development of sea state scaling under mixed wave conditions.
The goal of this study is to address the impact of swell on the wind–wave interaction and the wave state dependence of sea surface roughness. First, we analyze the wave age and wave steepness dependencies of the sea surface roughness at moderate wind speed under open ocean conditions, where mixed waves frequently occur (section 3), using data observed from a fixed platform (section 2). Then, a wind–wave coupling model is used to analyze the physics behind the wave state dependence, and the impact of swell is specifically addressed (section 4). Finally, the discussions and conclusions are presented in section 5.
2. Observations
a. Data
The measurements were taken from the Flux Observation Project in the South China Sea (FOPSCS), which was conducted in the coastal waters of the northern South China Sea (Fig. 1). An eddy covariance system was configured to measure the momentum, heat and gas flux across the air–sea interface, including an R.M. Young sonic anemometer mounted at a height of ~19 m above the sea surface, which was used to measure the wind velocity and sonic temperature at a sampling frequency of 10 Hz, and an open-path CO2/H2O analyzer operated synchronously to observe the water vapor concentration at the same sampling frequency. The wind stress can be directly determined from the measured wind fluctuations using the eddy covariance method, where the average time scale is determined hourly by the multiresolution decomposition method (Vickers and Mahrt 2003). A series of data processing and quality control procedures, including spike removal and tilt correction, were applied to determine the flux, as detailed in our previous studies (Li and Zhao 2016; Zou et al. 2017). Flux observations with winds blowing from the northwest were discarded as they may have been distorted by the platform. Four wind anemometers were mounted at heights of 31.3, 23.4, 20.0, and 16.4 m above mean sea level to measure the hourly averaged wind speed and direction.

Location of the platform.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Location of the platform.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Location of the platform.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
b. Processing method

Mean wind profile (stars) derived from the time-averaged normalized wind speed at different heights and the curve (blue line) fitted to the logarithmic profile of Eq. (9) using a regression method.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Mean wind profile (stars) derived from the time-averaged normalized wind speed at different heights and the curve (blue line) fitted to the logarithmic profile of Eq. (9) using a regression method.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Mean wind profile (stars) derived from the time-averaged normalized wind speed at different heights and the curve (blue line) fitted to the logarithmic profile of Eq. (9) using a regression method.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
c. Wind and wave conditions
The temporal variation in the wind and wave conditions are illustrated in Fig. 3, and the statistics were presented in Fig. 4. Figure 3a shows the variation in the wind speed (black stars and left-hand scale) and the significant wave height (blue circles and right-hand scale) during the observation time period. The winds generally lie in the range of 4–10 m s−1, representing moderate wind speeds that center at approximately 7 m s−1 (Fig. 4a). The significant wave height is in the range of 0.5–1.7 m. The swell index is estimated as the energy ratio between swell and wind wave, and ranges from 0 to 5 (Fig. 4b). The wave age, defined as the ratio of the phase velocity of the dominant wind wave to the friction velocity, covering a wide range from 9 to 35 (Fig. 4c), but the wave ages are dominated by relatively large values (>15) representing well-developed wind waves.

(a) Hourly mean wind speed U10 (black crosses) and significant wave height (blue circles) from January 2012 to July 2012; (b) wind direction (black crosses) and dominant wave direction (blue circles), where the directions signify where the winds and waves are from and the convention is clockwise relative to north; (c) phase speeds of dominant waves (black crosses) and wind waves (green crosses). Deep-water values are denoted by circles (blue for dominant waves and red for wind waves).
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

(a) Hourly mean wind speed U10 (black crosses) and significant wave height (blue circles) from January 2012 to July 2012; (b) wind direction (black crosses) and dominant wave direction (blue circles), where the directions signify where the winds and waves are from and the convention is clockwise relative to north; (c) phase speeds of dominant waves (black crosses) and wind waves (green crosses). Deep-water values are denoted by circles (blue for dominant waves and red for wind waves).
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
(a) Hourly mean wind speed U10 (black crosses) and significant wave height (blue circles) from January 2012 to July 2012; (b) wind direction (black crosses) and dominant wave direction (blue circles), where the directions signify where the winds and waves are from and the convention is clockwise relative to north; (c) phase speeds of dominant waves (black crosses) and wind waves (green crosses). Deep-water values are denoted by circles (blue for dominant waves and red for wind waves).
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Relevant wind and wave statistics in the dataset: Probability distribution of the (a) wind speed, (b) swell index, (c) wave age, and (d) deflection direction between wind waves and swells.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Relevant wind and wave statistics in the dataset: Probability distribution of the (a) wind speed, (b) swell index, (c) wave age, and (d) deflection direction between wind waves and swells.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Relevant wind and wave statistics in the dataset: Probability distribution of the (a) wind speed, (b) swell index, (c) wave age, and (d) deflection direction between wind waves and swells.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Figure 3b shows the directions of the winds (black stars) and dominant waves (spectral peaks, blue circles); the directions signify where the winds and waves from and the convention is clockwise relative to north. The winds generally come from the east with a mean azimuth of 107.5°, and the azimuths of the dominant waves are slightly larger with a mean value of 126.7°. It was found that ~63% of the spectral peaks are occupied by the swell components. The deflection angles between the wind waves and swells mostly lie in the range of 0°–30° (Fig. 4d), representing the predominance of following-swell conditions.
Figure 3c shows the phase speeds c0 of dominant waves and wind waves. The phase speed is calculated as c0 = ωp/kp, in which the subscript p denotes the spectral peak, and k is the wavenumber, which can be calculated interactively from the dispersion relation
3. The wave state dependence of the sea surface roughness
The wave state dependence of the normalized roughness length z0/Hs is given in Fig. 5. The wave age and wave steepness scalings were overlain for comparison. The bin-averaged data from Drennan et al. (2005) are also shown for reference. Note that the wave state dependence of sea surface roughness was investigated using either the data of the wave parameters corresponding to the wind-wave spectrum (e.g., D03) or those corresponding to the full spectrum (e.g., TY01). For the purpose of this test, these two choices were used separately. To check the wave-state-dependent scalings under different wave conditions, we partitioned the wave conditions according to the swell index, with values < 1 indicating wind-wave dominance and the remaining values indicating mixed wave condition.

Wave state dependence of the normalized z0: (a) wave steepness dependence using the full spectrum; (b) wave steepness dependence using the wind-wave spectrum, where the vertical dashed lines denote Hs/Lp = 0.02; (c) wave age dependence using the full spectrum; (d) wave age dependence using the wind-sea spectrum property. The bin-averaged data from Drennan et al. (2005) are overlain for reference (denoted by circles). The lines represent the proposed sea-state-dependent scalings. The error bars represent 60% of the standard deviation in a bin.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Wave state dependence of the normalized z0: (a) wave steepness dependence using the full spectrum; (b) wave steepness dependence using the wind-wave spectrum, where the vertical dashed lines denote Hs/Lp = 0.02; (c) wave age dependence using the full spectrum; (d) wave age dependence using the wind-sea spectrum property. The bin-averaged data from Drennan et al. (2005) are overlain for reference (denoted by circles). The lines represent the proposed sea-state-dependent scalings. The error bars represent 60% of the standard deviation in a bin.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Wave state dependence of the normalized z0: (a) wave steepness dependence using the full spectrum; (b) wave steepness dependence using the wind-wave spectrum, where the vertical dashed lines denote Hs/Lp = 0.02; (c) wave age dependence using the full spectrum; (d) wave age dependence using the wind-sea spectrum property. The bin-averaged data from Drennan et al. (2005) are overlain for reference (denoted by circles). The lines represent the proposed sea-state-dependent scalings. The error bars represent 60% of the standard deviation in a bin.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
For the wave steepness dependence, no clear trend can be identified when either the full spectrum or the wind-wave spectrum is considered (Figs. 5a,b), and the correlation coefficients are −0.03 and 0.09, respectively. The correlations are much higher when the wave age is considered (Figs. 5c,d): the correlation coefficient is 0.47 when the full spectrum is considered, and reaches 0.89 when the wind-wave spectrum is considered. Notable scatters are observed at a given wave age or wave steepness. Uncertainties may exist in discriminating wind waves and swells, Li and Zhao (2012) suggested that separated wind waves can be underestimated by up to ~6%, which would result in an overestimation of z0/Hs by up to ~12% in the wave steepness dependence and ~6% in the wave age dependence.
Drennan et al. (2005) proposed a wave steepness threshold of Hs/Lp = 0.02, and they suggested that the sea-state scalings (D03, TY01) hold for Hs/Lp ≥ 0.02. Since the wave steepness dependence is quite weak, this dependence cannot be readily described by either TY01 or T12. This may be because the observed Hs/Lp values are generally small with values generally less than 0.04, while TY01 and T12 are both developed for much steeper waves. The wave age dependence can be effectively described by D03 for Hs/Lp ≥ 0.02 when only the wind-wave spectrum is considered but is systematically underestimated by D03 for Hs/Lp < 0.02.
4. Physical analysis
a. Wind–wave interaction

Comparisons between the friction velocity observed on the platform and those estimated from the analytical model. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Comparisons between the friction velocity observed on the platform and those estimated from the analytical model. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Comparisons between the friction velocity observed on the platform and those estimated from the analytical model. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Long breaking wind-wave fronts generate as AFS region, in which the form stress supported by short wind waves and the viscous stress locally vanish (Reul et al. 2008; Buckley and Veron 2016). Considering the sheltering effect, the wave-induced stress is reduced to dτw-sheltering(k) = [(1 − q(k)] dτw(k), and the viscous stress is τυ-sheltering = [(1 − q(kb)]τυ, where kb ≈ 2π/0.15 rad m−1 is the wavenumber of the shortest breaking wind wave providing the AFS and q is the cumulative fraction of the sheltered surface describing the cumulative contribution of breaking sea waves to the sheltered zones and is determined from the local integration of k−1Λ(k)dk. Because q has an upper limit of 1, the threshold wavenumber k0 is introduced as q(k0) = 1. Assuming that breaking wind waves with the wavenumber k > k0 do not contribute to the total AFS stress as they are trapped in longer waves, the AFS stress is expressed as dτs-sheltering(k) = h(k − k0)dτs(k), where h(x) is the Heaviside step function.
The surface stress can be determined from the wind–wave coupling model via an iterative solution of Eqs. (12)–(17). The modeled surface stresses are compared with the measured momentum flux in Fig. 6 revealing a good agreement between them.
We plotted the normalized viscous stress and form drag separately against the wave steepness (Fig. 7). The viscous stress is generally larger than the form drag, particularly at low wave steepness. The viscous stress is clearly reduced by the AFS effect for Hs/Lp ≥ 0.02, as inferred from the difference between the viscous stress with and without considering the AFS effect. The AFS effect strongly depends on the wave steepness: it is stronger in cases with higher wave steepness. The AFS effect is rather weak for Hs/Lp < 0.02. This phenomenon is expected because the wave steepness is frequently used to characterize the probability of wave breaking: wave breaking is more likely to occur for a higher wave steepness. Banner et al. (2000) found that the occurrence of dominant wave breaking is closely related to the so-called significant spectral peak steepness δp, defined as δp = (Hpkp)/2, where

Normalized viscous stress (with and without considering the AFS effect for red and green crosses), the form drag (blue crosses) and the total surface stress (black crosses) of the model estimates vs Hs/Lp. VSA and VS denote the viscous stress with and without considering the AFS effect, respectively. FSA denotes the form drag considering the AFS effect, and TSS denotes the total surface stress. The data are averaged in bins of Hs/Lp.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Normalized viscous stress (with and without considering the AFS effect for red and green crosses), the form drag (blue crosses) and the total surface stress (black crosses) of the model estimates vs Hs/Lp. VSA and VS denote the viscous stress with and without considering the AFS effect, respectively. FSA denotes the form drag considering the AFS effect, and TSS denotes the total surface stress. The data are averaged in bins of Hs/Lp.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Normalized viscous stress (with and without considering the AFS effect for red and green crosses), the form drag (blue crosses) and the total surface stress (black crosses) of the model estimates vs Hs/Lp. VSA and VS denote the viscous stress with and without considering the AFS effect, respectively. FSA denotes the form drag considering the AFS effect, and TSS denotes the total surface stress. The data are averaged in bins of Hs/Lp.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
The form drag exhibits a generally increasing trend with the wave steepness. Wave-induced stress dominates the form drag and clearly increases with increasing wave steepness regardless of the increasing AFS effect. AFS stress generally increases with increasing wave steepness, and the trend appears to be more significant for large wave steepness (Hs/Lp > 0.02).
To summarize, for Hs/Lp < 0.02, the wave steepness dependence of the surface stress is mostly attributable to wave-induced stress, and the AFS effect is negligible, but the wave steepness dependence of the total drag is not significant because the form drag is much smaller than the viscous stress. For Hs/Lp > 0.02, a higher wave steepness corresponds to a higher form drag and a stronger AFS effect, the latter results in a lower viscous stress. As a result, the counteracting effect results in a relatively weaker dependence on the total surface stress than that on the form drag. This offers a possible explanation for the discrepancy between the observations and the two wave-steepness scalings (TY01; T12): these scalings were both developed for much steeper waves that typically occur under high winds or experience a strong shoaling effect, in which the form drag is dominant and the viscous stress is rather weak (Kudryavtsev and Makin 2007; Mueller and Veron 2009), which may contribute to a stronger wave steepness dependence.
b. The swell effect
To further interpret the wind–wave interaction in mixed waves, it was necessary to address the wind-wave properties in the presence of swell. In pure wind waves, there are similarities in the wind-wave growth (Zakharov et al. 2019): the dimensionless wave height and wave period are consistently combined with each other following H ~ TR, in which H = gHs/u2 and T = gTp/u, where u corresponds to either the wind speed or the frictional speed. Some well-known relations are those proposed by Hasselmann et al. (1976) with R = 5/3, Toba (1972) with R = 3/2, and Zakharov and Zaslavsky (1983) with R = 4/3. Theoretical analyses suggest that these relations are associated with different regimes of wind-wave coupling and represent different stages of wind-wave growth (Badulin and Grigorieva 2012; Badulin and Geogdzhaev 2019), as Hasselmann’s relation corresponds to a constant wave momentum flux and represents young waves, while the relation of Zakharov and Zaslavsky indicates the dominance of nonlinear wave–wave interactions corresponding to fully developed waves, and that of Toba corresponds to a constant wind flux into waves, representing an intermediate case between the abovementioned two cases.
Figure 8a shows a plot of the dimensionless wave height against the wave period using the wind-wave information and wind speed; the relations of Hasselmann et al. (1976) and Zakharov and Zaslavsky (1983) are overlain, Fig. 8b shows a similar dimensionless plot using the friction speed; here, the relation of Toba (1972) is overlain. In Fig. 8a, the best fit to the data points with Hs/Lp ≥ 0.02 results in an R value of 1.47, closer to that of Zakharov and Zaslavsky (1983). In Fig. 8b, the best fit to the data points with Hs/Lp ≥ 0.02 results in an R value of 1.50, which almost coincides with that of Toba (1972). The above results indicate that wind waves range from moderately developed to fully developed, which is consistent with the wave age conditions in Fig. 4c. The data points are consistently biased lower than those predicted by the wind-wave growth relations for Hs/Lp < 0.02, representing the suppression of wind-wave energy.

(a),(b) Normalized wave height H vs the normalized wave period T. The normalization uses the wind speed in (a) and friction speed in (b). The data are grouped with respect to the wave steepness threshold Hs/Lp = 0.02.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

(a),(b) Normalized wave height H vs the normalized wave period T. The normalization uses the wind speed in (a) and friction speed in (b). The data are grouped with respect to the wave steepness threshold Hs/Lp = 0.02.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
(a),(b) Normalized wave height H vs the normalized wave period T. The normalization uses the wind speed in (a) and friction speed in (b). The data are grouped with respect to the wave steepness threshold Hs/Lp = 0.02.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
There are three main impacts of swell on wind waves. First, swell can interact with the airflow and thus change the wind input (Chen and Belcher 2000), but data points from these cases have been excluded from analysis. The second impact is through the long wave modulation of the wave breaking process (e.g., Dulov et al. 2002); however, this modulation should not result in a reduction in the wind-wave energy because field data implied reduced breaking wave dissipation through whitecap coverage observation in the presence of swell (e.g., Sugihara et al. 2007; Goddijn-Murphy et al. 2011). Third, through the direct interaction between swell and wind wave, field data (Young 2006; Badulin and Grigorieva 2012) and theoretical models (Badulin 2008) have shown strong interactions, as swell can absorb energy from locally generated wind wave. The interaction can be related to the dominant wave period or height ratio between swell and wind wave (Masson 1993; Carlsson et al. 2010; Potter 2015).
To analyze the influence of swell on wind-wave growth, a height ratio HR is defined as the ratio of the measured wind-wave height Hs to that expected from a pure wind wave, which was estimated from the relation of Toba (1972). This ratio was plotted against the frequency ratio (FR = fp_swell/fp_windsea) and swell index (SI = Eswell/Ewindsea) between swell and wind wave, as shown in Figs. 9a and 9b, respectively. A negative correlation exists between HR and FR, with correlation coefficients of −0.38 for the scatter points and −0.95 for the bin-averaged values. No clear trend is observed between HR and SI, with correlation coefficients of −0.06 for the scatter points and 0.30 for the bin-averaged values. These results indicate that the interaction between wind wave and swell is stronger when their spectral peaks are closer. The height ratio HR is consistently below the value of 1 for FR > 0.6, which is in agreement with the modeling study of Masson (1993). At low FR values (FR < 0.4), for which the spectral peaks of swells and wind waves are well separated, the interactions between swells and the energy-containing wind waves are supposed to be rather weak.

Height ratio HR vs (a) the frequency ratio and (b) the swell index. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1

Height ratio HR vs (a) the frequency ratio and (b) the swell index. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
Height ratio HR vs (a) the frequency ratio and (b) the swell index. The error bars represent the standard deviations of the mean.
Citation: Journal of Physical Oceanography 50, 11; 10.1175/JPO-D-20-0102.1
5. Discussions and conclusions
The wave state dependence of sea surface roughness and its parameterization have been widely studied for many years, but they are mostly limited to wind-wave-dominant conditions. Directly extrapolating this relationship to mixed waves or swell-dominant conditions would introduce large uncertainty. Considering the frequent occurrence of mixed waves in the open ocean, it is not surprising that a wind-speed-dependent scaling is still popularly used (e.g., Foreman and Emeis 2010; Andreas et al. 2012; Edson et al. 2013).
In this study, we examined the wave age and wave steepness dependencies of the sea surface roughness using observations taken from a fixed platform. The observations were characterized by moderate wind speeds, moderately developed to fully developed wind waves and mixed wave conditions, which are quite common in the open ocean (e.g., Monahan 2006; Semedo et al. 2011). No obvious trend was found in the wave steepness dependence, using either the full wave spectrum or the wind-wave spectrum alone. The wave age scaling (D03), proposed for pure wind-wave conditions, was found to hold under mixed wave conditions when the wind-wave spectrum was used and the wave steepness was above a threshold (Hs/Lp > 0.02), in which the wind-wave properties conform to the wind-wave growth relation. A wind–wave coupling model analysis suggested that the viscous stress tends to decrease with increasing wave steepness due to the airflow separation induced by the breaking of the dominant wind-wave front, which counteracts the increasing trend of the form drag. The Hs/Lp threshold is assumed to be related to the onset of dominant wave breaking.
The development of wind wave is determined by the relative importance of the wind input and wave dissipation, and the nonlinear wave–wave interaction plays a critical role in the self-adjustment process during the evolution of wind waves. As a result, the wind-wave property follows the wind-wave growth relation. In the presence of swell, the measured wind-wave properties follow the wind-wave growth relation quite well for Hs/Lp > 0.02, but some deviations are indicated for Hs/Lp < 0.02, probably resulting from the nonlinear energy transfer from wind wave to swell. The dependency of the interaction between wind wave and swell on the dominant wave period ratio is illustrated, however, this is a preliminary result, and other possible dependencies on wave characteristic parameters, including the deflection angle and the directional spread and frequency bandwidth of each wave component, should be further diagnosed.
As a result of wave–wave interactions, wind waves are smoothed due to the “lost” energy, reducing the wave-induced stress, but the viscous stress is enhanced due to the reduced airflow separation effect. These combined effects enhance the total drag. Note that the swell effect is quite different from that proposed by García-Nava et al. (2009, 2012a), who suggested that the wind stress is reduced through the impact of swell on short waves. The different wave conditions between their studies and the FOPSCS may account for the different effects of swell. In their studies, the spectral peaks of swells and wind waves are well departed, and thus, the wave–wave interaction is expected to be weak, and the form drag dominates the total wind stress under strong winds and typically young wind waves (Kudryavtsev and Makin 2007; Mueller and Veron 2009). They further showed that the short waves that support the form drag are considerably suppressed by swells. We did not discover the overall attenuation or enhancement of the intensity of short waves in the presence of swells from the FOPSCS (not shown). This may be because swells act to enhance short waves by reducing the AFS effect, which may counteract the suppression of short waves by swell.
The generality of these results should be treated with caution when they were extrapolated to other mixed wave conditions. The observations were characterized by well-developed wind waves at moderate wind speeds, as a result, the wind waves were not as steep as those under strong wind forcing with limited wind fetch or duration (TY01; T12) or those that experience the shoaling effect (Anctil and Donelan 1996), in which the form drag can be dominant (e.g., Suzuki et al. 2014; Babanin et al. 2018). In addition, the swells in the study area are not as old as those propagating thousands of miles across the wide ocean because they are mostly generated within the South China Sea (Li et al. 2018b), and their steepness are enhanced because they experience the shoaling effect near the platform. Therefore, a more comprehensive dataset is required to achieve a better understanding of the wave state dependence of the surface roughness under mixed wave condition.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (41976010, U1706216, 41606024, 41876010, 41806028), the National Key Research and Development Program of China (Grants 2016YFC1402000, 2018YFC1401002), the CAS Strategic Priority Project (XDA19060202), the Public Science and Technology Research Funds Projects of Ocean (201505007), and the Open Fund Project of Key Laboratory of Marine Environmental Information Technology.
REFERENCES
Anctil, F., and M. A. Donelan, 1996: Air–water momentum flux observations over shoaling waves. J. Phys. Oceanogr., 26, 1344–1353, https://doi.org/10.1175/1520-0485(1996)026<1344:AMFOOS>2.0.CO;2.
Andreas, E. L, L. Mahrt, and D. Vickers, 2012: A new drag relation for aerodynamically rough flow over the ocean. J. Atmos. Sci., 69, 2520–2537, https://doi.org/10.1175/JAS-D-11-0312.1.
Babanin, A. V., and I. R. Young, 2005: Two-phase behavior of the spectral dissipation of wind waves. Proc. Fifth Int. Symp. on Ocean Wave Measurement and Analysis, Madrid, Spain, CEDEX and ASCI CORPI, 51.
Babanin, A. V., J. McConochie, and D. Chalikov, 2018: Winds near the surface of waves: Observations and modeling. J. Phys. Oceanogr., 48, 1079–1088, https://doi.org/10.1175/JPO-D-17-0009.1.
Badulin, S. I., A. O. Korotkevich, D. Resio, and V. E. Zakharov, 2008: Wave-wave interactions in wind-driven mixed seas. Proc. Rogue Waves 2008 Workshop, Brest, France, IFREMER, 77–85.
Badulin, S. I., and V. G. Grigorieva, 2012: On discriminating swell and wind-driven seas in Voluntary Observing Ship data. J. Geophys. Res., 117, C00J29, https://doi.org/10.1029/2012JC007937.
Badulin, S. I., and V. G. Geogdzhaev, 2019: Evaluation of wind wave growth parameters basing on spectral fluxes. Radiophys. Quantum Electron., 61, 545–552, https://doi.org/10.1007/s11141-019-09915-8.
Banner, M. L., 1990: Equilibrium spectra of wind waves. J. Phys. Oceanogr., 20, 966–984, https://doi.org/10.1175/1520-0485(1990)020<0966:ESOWW>2.0.CO;2.
Banner, M. L., A. V. Babanin, and I. R. Young, 2000: Breaking probability for dominant waves on the sea surface. J. Phys. Oceanogr., 30, 3145–3160, https://doi.org/10.1175/1520-0485(2000)030<3145:BPFDWO>2.0.CO;2.
Buckley, M. P., and F. Veron, 2016: Structure of the airflow above surface waves. J. Phys. Oceanogr., 46, 1377–1397, https://doi.org/10.1175/JPO-D-15-0135.1.
Carlsson, B., Y. Papadimitrakis, and A. Rutgersson, 2010: Evaluation of a roughness length model and sea surface properties with data from the Baltic Sea. J. Phys. Oceanogr., 40, 2007–2024, https://doi.org/10.1175/2010JPO4340.1.
Charnock, H., 1955: Wind stress on a water surface. Quart. J. Roy. Meteor. Soc., 81, 639–640, https://doi.org/10.1002/qj.49708135027.
Chen, G., and S. E. Belcher, 2000: Effects of long waves on wind-generated waves. J. Phys. Oceanogr., 30, 2246–2256, https://doi.org/10.1175/1520-0485(2000)030<2246:EOLWOW>2.0.CO;2.
Donelan, M. A., 1990: Air–sea interaction. Ocean Engineering Science, Parts A and B, B. LeMehaute and D. M. Hanes, Eds., The Sea—Ideas and Observations on Progress in the Study of the Seas, Vol. 9, John Wiley and Sons, 239–292.
Donelan, M. A., F. W. Dobson, S. D. Smith, and R. J. Anderson, 1993: On the dependence of sea surface roughness on wave development. J. Phys. Oceanogr., 23, 2143–2149, https://doi.org/10.1175/1520-0485(1993)023<2143:OTDOSS>2.0.CO;2.
Drennan, W. M., H. C. Graber, D. Hauser, and C. Quentin, 2003: On the wave age dependence of wind stress over pure wind seas. J. Geophys. Res., 108, 8062, https://doi.org/10.1029/2000JC000715.
Drennan, W. M., P. K. Taylor, and M. J. Yelland, 2005: Parameterizing the sea surface roughness. J. Phys. Oceanogr., 35, 835–848, https://doi.org/10.1175/JPO2704.1.
Dulov, V., V. Kudryavtsev, and A. Bol’Shakov, 2002: A field study of whitecap coverage and its modulations by energy containing surface waves. Gas Transfer at Water Surfaces, Geophys. Monogr., Vol. 127, Amer. Geophys. Union, 187–192.
Edson, J. B., and Coauthors, 2013: On the exchange of momentum over the open ocean. J. Phys. Oceanogr., 43, 1589–1610, https://doi.org/10.1175/JPO-D-12-0173.1.
Foreman, R. J., and S. Emeis, 2010: Revisiting the definition of the drag coefficient in the marine atmospheric boundary layer. J. Phys. Oceanogr., 40, 2325–2332, https://doi.org/10.1175/2010JPO4420.1.
Gao, Z., Q. Wang, and S. Wang, 2006: An alternative approach to sea surface aerodynamic roughness. J. Geophys. Res., 111, D22108, https://doi.org/10.1029/2006JD007323.
García-Nava, H., F. J. Ocampo-Torres, P. Osuna, and M. A. Donelan, 2009: Wind stress in the presence of swell under moderate to strong wind conditions. J. Geophys. Res., 114, C12008, https://doi.org/10.1029/2009JC005389.
García-Nava, H., F. J. Ocampo-Torres, P. A. Hwang, and P. Osuna, 2012a: Reduction of wind stress due to swell at high wind conditions. J. Geophys. Res., 117, C00J11, https://doi.org/10.1029/2011JC007833.
García-Nava, H., F. J. Ocampo-Torres, and P. A. Hwang, 2012b: On the parameterization of the drag coefficient in mixed seas. Sci. Mar., 76, 177–186, https://doi.org/10.3989/scimar.03615.19F.
Garratt, J. R. D., 1977: Review of drag coefficients over oceans and continents. Mon. Wea. Rev., 105, 915–929, https://doi.org/10.1175/1520-0493(1977)105<0915:RODCOO>2.0.CO;2.
Gemmrich, J., C. J. Zappa, M. L. Banner, and R. P. Morison, 2013: Wave breaking in developing and mature seas. J. Geophys. Res. Oceans, 118, 4542–4552, https://doi.org/10.1002/jgrc.20334.
Goddijn-Murphy, L., D. K. Woolf, and A. H. Callaghan, 2011: Parameterizations and algorithms for oceanic whitecap coverage. J. Phys. Oceanogr., 41, 742–756, https://doi.org/10.1175/2010JPO4533.1.
Grare, L., L. Lenain, and W. K. Melville, 2013: Wave-coherent airflow and critical layers over ocean waves. J. Phys. Oceanogr., 43, 2156–2172, https://doi.org/10.1175/JPO-D-13-056.1.
Guan, C., and L. Xie, 2004: On the linear parameterization of drag coefficient over sea surface. J. Phys. Oceanogr., 34, 2847–2851, https://doi.org/10.1175/JPO2664.1.
Hanley, K. E., S. E. Belcher, and P. P. Sullivan, 2010: A global climatology of wind–wave interaction. J. Phys. Oceanogr., 40, 1263–1282, https://doi.org/10.1175/2010JPO4377.1.
Hanson, J. L., and O. M. Phillips, 2001: Automated analysis of ocean surface directional wave spectra. J. Atmos. Oceanic Technol., 18, 277–293, https://doi.org/10.1175/1520-0426(2001)018<0277:AAOOSD>2.0.CO;2.
Hara, T., and S. E. Belcher, 2002: Wind forcing in the equilibrium range of wind-wave spectra. J. Fluid Mech., 470, 223–245, https://doi.org/10.1017/S0022112002001945.
Hasselmann, K., D. B. Ross, P. Müller, and W. Sell, 1976: A parametric wave prediction model. J. Phys. Oceanogr., 6, 200–228, https://doi.org/10.1175/1520-0485(1976)006<0200:APWPM>2.0.CO;2.
Hasselmann, S., C. Brüning, K. Hasselmann, and P. Heimbach, 1996: An improved algorithm for the retrieval of ocean wave spectra from synthetic aperture radar image spectra. J. Geophys. Res., 101, 16 615–16 629, https://doi.org/10.1029/96JC00798.
Högström, U., 1988: Non-dimensional wind and temperature profiles in the atmospheric surface layer: A re-evaluation. Bound.-Layer Meteor., 42, 55–78, https://doi.org/10.1007/BF00119875.
Högström, U., and Coauthors, 2008: Momentum fluxes and wind gradients in the marine boundary layer—A multi-platform study. Boreal Environ. Res., 13, 475–502.
Hsu, S. A., 1974: A dynamic roughness equation and its application to wind stress determination at the air–sea interface. J. Phys. Oceanogr., 4, 116–120, https://doi.org/10.1175/1520-0485(1974)004<0116:ADREAI>2.0.CO;2.
Hwang, P. A., H. García-Nava, and F. J. Ocampo-Torres, 2011a: Dimensionally consistent similarity relation of ocean surface friction coefficient in mixed seas. J. Phys. Oceanogr., 41, 1227–1238, https://doi.org/10.1175/2011JPO4566.1.
Hwang, P. A., H. García-Nava, and F. J. Ocampo-Torres, 2011b: Observations of wind wave development in mixed seas and unsteady wind forcing. J. Phys. Oceanogr., 41, 2343–2362, https://doi.org/10.1175/JPO-D-11-044.1.
Johnson, H. K., J. Hojstrup, H. J. Vested, and S. E. Larsen, 1998: On the dependence of sea roughness on wind waves. J. Phys. Oceanogr., 28, 1702–1716, https://doi.org/10.1175/1520-0485(1998)028<1702:OTDOSS>2.0.CO;2.
Jones, I. S. F., and Y. Toba, 2001: Wind Stress Over the Ocean. Cambridge University Press, 307 pp.
Kahma, K. K., M. A. Donelan, W. M. Drennan, and E. A. Terray, 2016: Evidence of energy and momentum flux from swell to wind. J. Phys. Oceanogr., 46, 2143–2156, https://doi.org/10.1175/JPO-D-15-0213.1.
Kudryavtsev, V. N., and V. K. Makin, 2001: The impact of air-flow separation on the drag of the sea surface. Bound.-Layer Meteor., 98, 155–171, https://doi.org/10.1023/A:1018719917275.
Kudryavtsev, V. N., and V. K. Makin, 2007: Aerodynamic roughness of the sea surface at high winds. Bound.-Layer Meteor., 125, 289–303, https://doi.org/10.1007/s10546-007-9184-7.
Kudryavtsev, V. N., B. Chapron, and V. Makin, 2014: Impact of wind waves on the air-sea fluxes: A coupled model. J. Geophys. Res. Oceans, 119, 1217–1236, https://doi.org/10.1002/2013JC009412.
Kukulka, T., and T. Hara, 2008a: The effect of breaking waves on a coupled model of wind and ocean surface waves. Part I: Mature seas. J. Phys. Oceanogr., 38, 2145–2163, https://doi.org/10.1175/2008JPO3961.1.
Kukulka, T., and T. Hara, 2008b: The effect of breaking waves on a coupled model of wind and ocean surface waves. Part II: Growing seas. J. Phys. Oceanogr., 38, 2164–2184, https://doi.org/10.1175/2008JPO3962.1.
Li, S., and D. Zhao, 2012: Comparison of spectral partitioning techniques for wind wave and swell. Mar. Sci. Bull., 14, 24–36.
Li, S., and D. Zhao, 2016: Gas transfer velocity in the presence of wave breaking. Tellus, 68B, 27034, https://doi.org/10.3402/tellusb.v68.27034.
Li, S., H. Shen, Y. Hou, Y. He, and F. Bi, 2018a: Sea surface wind speed and sea state retrievals from dual-frequency altimeter and its preliminary application in global view of wind-sea and swell distributions. Int. J. Remote Sens., 39, 3076–3093, https://doi.org/10.1080/01431161.2018.1433889.
Li, S., F. Bi, Y. Hou, and H. Yong, 2018b: Characterization of wind-sea and swell in the South China sea based on HY-2 satellite data. J. Coastal Res., 84, 58–62, https://doi.org/10.2112/SI84-008.1.
Masson, D., 1993: On the nonlinear coupling between swell and wind waves. J. Phys. Oceanogr., 23, 1249–1258, https://doi.org/10.1175/1520-0485(1993)023<1249:OTNCBS>2.0.CO;2.
Melville, W. K., and P. Matusov, 2002: Distribution of breaking waves at the ocean surface. Nature, 417, 58–63, https://doi.org/10.1038/417058a.
Monahan, A. H., 2006: The probability distribution of sea surface wind speeds. Part I: Theory and sea winds observations. J. Climate, 19, 497–520, https://doi.org/10.1175/JCLI3640.1.
Monin, A. S., and A. M. Yaglom, 1971: Statistical Fluid Mechanics. Vol. 1, MIT Press, 769 pp.
Mueller, J. A., and F. Veron, 2009: Nonlinear formulation of the bulk surface stress over breaking waves: Feedback mechanisms from air-flow separation. Bound.-Layer Meteor., 130, 117, https://doi.org/10.1007/s10546-008-9334-6.
Oost, W. A., G. J. Komen, C. M. J. Jacobs, and C. V. Oort, 2002: New evidence for a relation between wind stress and wave age from measurements during ASGAMAGE. Bound.-Layer Meteor., 103, 409–438, https://doi.org/10.1023/A:1014913624535.
Pedersen, T., and E. Siegel, 2008: Wave Measurements from Subsurface Buoys. 2008 IEEE/OES 9th Working Conf. on Current Measurement Technology, Charleston, SC, IEEE, 224–233, https://doi.org/10.1109/CCM.2008.4480872.
Pedersen, T., A. Lohrmann, and H. E. Krogstad, 2005: Wave measurements from a subsurface platform. Proc. WAVES 2005, Madrid, Spain, CEDEX, 1–10.
Phillips, O. M., 1985: Spectral and statistical properties of the equilibrium range in wind-generated gravity waves. J. Fluid Mech., 156, 505–531, https://doi.org/10.1017/S0022112085002221.
Potter, H., 2015: Swell and the drag coefficient. Ocean Dyn., 65, 375–384, https://doi.org/10.1007/s10236-015-0811-4.
Reul, N., H. Branger, and J. Giovanangeli, 2008: Air flow structure over short-gravity breaking water waves. Bound.-Layer Meteor., 126, 477–505, https://doi.org/10.1007/s10546-007-9240-3.
Rieder, K. F., J. A. Smith, and R. A. Weller, 1996: Some evidence of colinear wind stress and wave breaking. J. Phys. Oceanogr., 26, 2519–2524, https://doi.org/10.1175/1520-0485(1996)026<2519:SEOCWS>2.0.CO;2.
Semedo, A., K. Sušelj, A. Rutgersson, and A. Sterl, 2011: A global view on the wind sea and swell climate and variability from ERA-40. J. Climate, 24, 1461–1479, https://doi.org/10.1175/2010JCLI3718.1.
Smedman, A., U. Högström, H. Bergström, A. Rutgersson, K. K. Kahma, and H. Pettersson, 1999: A case study of air-sea interaction during swell conditions. J. Geophys. Res., 104, 25 833–25 851, https://doi.org/10.1029/1999JC900213.
Smedman, A., X. Guo-Larsén, U. Högström, K. K. Kahma, and H. Pettersson, 2003: Effect of sea state on the momentum exchange over the sea during neutral conditions. J. Geophys. Res., 108, 3367, https://doi.org/10.1029/2002JC001526.
Smith, S. D., and Coauthors, 1992: Sea surface wind stress and drag coefficients: The hexos results. Bound.-Layer Meteor., 60, 109–142, https://doi.org/10.1007/BF00122064.
Sugihara, Y., H. Tsumori, T. Ohga, H. Yoshioka, and S. Serizawa, 2007: Variation of whitecap coverage with wave-field conditions. J. Mar. Syst., 66, 47–60, https://doi.org/10.1016/j.jmarsys.2006.01.014.
Suzuki, N., T. Hara, and P. P. Sullivan, 2014: Impact of dominant breaking waves on air–sea momentum exchange and boundary layer turbulence at high winds. J. Phys. Oceanogr., 44, 1195–1212, https://doi.org/10.1175/JPO-D-13-0146.1.
Takagaki, N., S. Komori, N. Suzuki, K. Iwano, T. Kuramoto, S. Shimada, R. Kurose, and K. Takahashi, 2012: Strong correlation between the drag coefficient and the shape of the wind sea spectrum over a broad range of wind speeds. Geophys. Res. Lett., 39, L23604, https://doi.org/10.1029/2012GL053988.
Taylor, P. K., and M. J. Yelland, 2001: The dependence of sea surface roughness on the height and steepness of the waves. J. Phys. Oceanogr., 31, 572–590, https://doi.org/10.1175/1520-0485(2001)031<0572:TDOSSR>2.0.CO;2.
Thomson, J., J. R. Gemmrich, and A. T. Jessup, 2009: Energy dissipation and the spectral distribution of whitecaps. Geophys. Res. Lett., 36, L11601, https://doi.org/10.1029/2009GL038201.
Toba, Y., 1972: Local balance in the air-sea boundary processes III Part I. On the growth process of wind waves. J. Oceanogr., 29, 209–225.
Toba, Y., N. Iida, H. Kawamura, N. Ebuchi, and I. S. F. Jones, 1990: The wave dependence of sea-surface wind stress. J. Phys. Oceanogr., 20, 705–721, https://doi.org/10.1175/1520-0485(1990)020<0705:WDOSSW>2.0.CO;2.
Vickers, D., and L. Mahrt, 2003: The cospectral gap and turbulent flux calculations. J. Atmos. Oceanic Technol., 20, 660–672, https://doi.org/10.1175/1520-0426(2003)20<660:TCGATF>2.0.CO;2.
Voorrips, A. C., V. K. Makin, and S. Hasselmann, 1997: Assimilation of wave spectra from pitch-and-roll buoys in a North Sea wave model. J. Geophys. Res., 102, 5829–5849, https://doi.org/10.1029/96JC03242.
Wu, J., 1980: Wind-stress coefficients over sea surface near neutral conditions—A revisit. J. Phys. Oceanogr., 10, 727–740, https://doi.org/10.1175/1520-0485(1980)010<0727:WSCOSS>2.0.CO;2.
Young, I. R., 2006: Directional spectra of hurricane wind waves. J. Geophys. Res., 111, C08020, https://doi.org/10.1029/2006JC003540.
Young, I. R., and M. A. Donelan, 2018: On the determination of global ocean wind and wave climate from satellite observations. Remote Sens. Environ., 215, 228–241, https://doi.org/10.1016/j.rse.2018.06.006.
Zakharov, V. E., and M. M. Zaslavsky, 1983: Dependence of wave parameters on the wind velocity, duration of its action and fetch in the weak turbulence theory of water waves. Izv., Atmos. Ocean. Phys., 19, 300–306.
Zakharov, V. E., S. I. Badulin, V. V. Geogjaev, and A. N. Pushkarev, 2019: Weak-turbulent theory of wind-driven sea. Earth Space Sci., 6, 540–556, https://doi.org/10.1029/2018EA000471.
Zou, Z., D. Zhao, B. Liu, J. A. Zhang, and J. Huang, 2017: Observation-based parameterization of air-sea fluxes in terms of wind speed and atmospheric stability under low-to-moderate wind conditions. J. Geophys. Res. Oceans, 122, 4123–4142, https://doi.org/10.1002/2016JC012399.
Zou, Z., D. Zhao, J. A. Zhang, S. Li, Y. Cheng, H. Lv, and X. Ma, 2018: The influence of swell on the atmospheric boundary layer under nonneutral conditions. J. Phys. Oceanogr., 48, 925–936, https://doi.org/10.1175/JPO-D-17-0195.1.