Key Role of Subdaily Wind Variability for Tropical Surface Wind Stress

Yunwei Yan aKey Laboratory of Marine Hazards Forecasting, Ministry of Natural Resources, Hohai University, Nanjing, China
bCollege of Oceanography, Hohai University, Nanjing, China

Search for other papers by Yunwei Yan in
Current site
Google Scholar
PubMed
Close
,
Xiangzhou Song aKey Laboratory of Marine Hazards Forecasting, Ministry of Natural Resources, Hohai University, Nanjing, China
bCollege of Oceanography, Hohai University, Nanjing, China

Search for other papers by Xiangzhou Song in
Current site
Google Scholar
PubMed
Close
, and
Marilena Oltmanns cNational Oceanography Centre, Southampton, United Kingdom

Search for other papers by Marilena Oltmanns in
Current site
Google Scholar
PubMed
Close
Free access

Abstract

High-frequency observations of surface winds over the open ocean are available only at limited locations. However, these observations are essential for assessing atmospheric influences on the ocean, validating reanalysis products, and building parameterization schemes. By analyzing high-frequency measurements from the Global Tropical Moored Buoy Array, the effects of subdaily winds on the mean surface wind stress magnitude are systematically examined. Subdaily winds account for 12.4% of the total stress magnitude on average. The contribution is enhanced over the intertropical convergence zone and reaches a maximum (28.5%) in the equatorial western Pacific. The magnitude of the contribution is primarily determined by the kinetic energy of subdaily winds. Compared to the buoy observations, the ERA5 and MERRA2 subdaily winds underestimate this contribution by 51% and 63% due to underestimations of subdaily kinetic energy, leading to 7% and 8% underestimations in the total stress magnitude, respectively. Two new gustiness parameterization schemes related to precipitation are developed to account for the effect of subdaily winds, explaining ∼80% of the contribution from subdaily winds. Considering the importance of wind stress for ocean–atmosphere interactions, the inclusion of these parameterization schemes in climate models is expected to substantially improve simulations of large-scale climate variability.

Significance Statement

Surface wind stress drives upper-ocean circulation, which is critical for the redistribution of mass, momentum, and energy in the ocean. Moreover, it is one of the key factors controlling oceanic turbulent mixing and therefore has significant impacts on the distribution of temperature, salinity, and associated ocean variability. Using high-resolution buoy observations, this study highlights the importance of subdaily winds for integrated wind stress estimates. In addition, it finds that current state-of-the-art and widely used reanalysis products largely underestimate the effect of subdaily winds. Two new parameterization schemes are developed, leading to a better representation of the effect of subdaily winds. Including the proposed parameterization schemes in climate models is expected to substantially improve their simulations of large-scale climate variability.

© 2023 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Yunwei Yan, yunwei.yan@hhu.edu.cn

Abstract

High-frequency observations of surface winds over the open ocean are available only at limited locations. However, these observations are essential for assessing atmospheric influences on the ocean, validating reanalysis products, and building parameterization schemes. By analyzing high-frequency measurements from the Global Tropical Moored Buoy Array, the effects of subdaily winds on the mean surface wind stress magnitude are systematically examined. Subdaily winds account for 12.4% of the total stress magnitude on average. The contribution is enhanced over the intertropical convergence zone and reaches a maximum (28.5%) in the equatorial western Pacific. The magnitude of the contribution is primarily determined by the kinetic energy of subdaily winds. Compared to the buoy observations, the ERA5 and MERRA2 subdaily winds underestimate this contribution by 51% and 63% due to underestimations of subdaily kinetic energy, leading to 7% and 8% underestimations in the total stress magnitude, respectively. Two new gustiness parameterization schemes related to precipitation are developed to account for the effect of subdaily winds, explaining ∼80% of the contribution from subdaily winds. Considering the importance of wind stress for ocean–atmosphere interactions, the inclusion of these parameterization schemes in climate models is expected to substantially improve simulations of large-scale climate variability.

Significance Statement

Surface wind stress drives upper-ocean circulation, which is critical for the redistribution of mass, momentum, and energy in the ocean. Moreover, it is one of the key factors controlling oceanic turbulent mixing and therefore has significant impacts on the distribution of temperature, salinity, and associated ocean variability. Using high-resolution buoy observations, this study highlights the importance of subdaily winds for integrated wind stress estimates. In addition, it finds that current state-of-the-art and widely used reanalysis products largely underestimate the effect of subdaily winds. Two new parameterization schemes are developed, leading to a better representation of the effect of subdaily winds. Including the proposed parameterization schemes in climate models is expected to substantially improve their simulations of large-scale climate variability.

© 2023 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Yunwei Yan, yunwei.yan@hhu.edu.cn

1. Introduction

Wind stress at the sea surface is the main driver of global ocean circulation, redistributing mass, momentum, and energy in the ocean (e.g., Gill 1982; Wunsch and Ferrari 2004). Moreover, representing the momentum flux between the atmosphere and the ocean, surface wind stress is one of the key factors controlling oceanic turbulent mixing and therefore exerts significant impacts on ocean temperature and salinity (e.g., Ouni et al. 2021; Zhou et al. 2018). Due to the nonlinear dependence of wind stress on wind, the high-frequency variability of winds has a large influence on low-frequency wind stress (e.g., Ponte and Rosen 2004; Zhai et al. 2012). Many studies have shown that the synoptic variability in surface winds makes a significant contribution to the mean wind stress, with far-reaching implications for ocean circulation and turbulent mixing (e.g., Duteil 2019; Lin et al. 2018, 2020; Ponte and Rosen 2004; Wu et al. 2016, 2020; Zhai 2013; Zhai and Wunsch 2013; Zhai et al. 2012). For example, Duteil (2019) found that removing synoptic winds (from 2 days to 1 month) in the tropical Pacific decreased the mean wind stress by up to 20%, weakening the wind-driven ocean circulation by 20%.

Likewise, neglecting subdaily wind variability may have severe effects on mean wind stress estimates (e.g., Eymard et al. 2003; Ledvina et al. 1993; Ouni et al. 2021; Zhai et al. 2012; Zhou et al. 2018). Based on ship/buoy measurements during March and April 1998 in the northwestern Mediterranean Sea, Eymard et al. (2003) evaluated the impact of the subdaily variability of air–sea parameters on air–sea flux estimates. They found a strong underestimation of the mean momentum flux when the vector-averaged winds were employed; the mean friction velocity decreased by 22% (Table 3 in their paper). Due to the nonlinear effect, subdaily winds exert significant impacts on wind power input to the ocean and ocean temperature and salinity at longer time scales (e.g., Lee and Liu 2005; Ouni et al. 2021; Zhai et al. 2012; Zhou et al. 2018). Zhai et al. (2012) showed that including subdaily winds in analyses led to a one-third increase in the estimated total wind power input to the ocean through an increased mean wind stress (Table 1 in their paper). Zhou et al. (2018) conducted sensitivity experiments using a one-dimensional ocean model at 15.5°N, 61.5°E and found that including subdaily signals in the wind field lowered the daily mean sea surface temperature (SST) by 0.8°C on average, primarily due to intensified wind stress magnitude and turbulent heat fluxes.

To date, however, the contribution of subdaily winds to the mean wind stress has not been quantified on a global scale. The Global Tropical Moored Buoy Array (GTMBA) provides long-term high-frequency (10 min and hourly) measurements of air–sea variables across the entirety of the tropical oceans (McPhaden et al. 2010), allowing us to address this problem. Buoy observations have been used to investigate some of the characteristics of subdaily winds, such as diurnal and semidiurnal cycles and wind gusts (e.g., Christophersen et al. 2020; Deser and Smith 1998; Giglio et al. 2017, 2022; Joseph et al. 2021; Serra et al. 2007; Turk et al. 2021; Ueyama and Deser 2008). Here, these data are used to investigate the relation between subdaily winds and the mean wind stress. Thus, the first objective of this paper is to systematically quantify the effects of subdaily winds on the mean wind stress magnitude across the entirety of the tropical oceans using GTMBA observations.

Given the importance of subdaily winds, their effect on the ocean has been taken into account in numerical simulations by employing wind forcing with subdaily variability (e.g., Ouni et al. 2021; Thushara and Vinayachandran 2014; Yu et al. 2020). Wind forcing with subdaily variability is largely available from the most recent atmospheric reanalysis products with hourly temporal resolutions, such as the fifth-generation ECMWF reanalysis product (ERA5; Hersbach et al. 2020) and the Modern-Era Retrospective analysis for Research and Applications product, Version 2 (MERRA2; Gelaro et al. 2017). However, the accuracy of subdaily winds in these reanalysis products is still unknown. Thus, the second objective of this paper is to evaluate the effect of subdaily winds on the mean wind stress magnitude in the ERA5 and MERRA2 reanalysis products.

A more efficient and less expensive approach to account for subdaily winds is to parameterize them by using a gustiness correction for the daily vector-averaged wind speed (e.g., Cronin et al. 2006; Praveen Kumar et al. 2012, 2013). For example, Praveen Kumar et al. (2012) proposed a climatological gustiness parameterization scheme related to SST and applied the gustiness correction to account for unresolved subdaily winds in the TropFLUX reanalysis product. An evaluation showed that including the gustiness correction increased the mean wind stress magnitude by ∼5% (Fig. 2d in Praveen Kumar et al. 2013). In this study, we relate gustiness to the precipitation rate, which has a close relation to the momentum flux associated with subdaily winds (see section 3). Two new gustiness parameterization schemes are proposed: one is time independent, similar to that of Praveen Kumar et al. (2012), and the other is time dependent. This is the third objective of this study.

This paper is structured as follows: after a short description of the data and methods in section 2, section 3 will present the major results of this study. First, the contribution of subdaily winds to the mean wind stress magnitude is estimated from the GTMBA observations, and then the effects of subdaily winds in the ERA5 and MERRA2 products are evaluated. Finally, two new parameterization schemes to represent unresolved subdaily winds are developed. Conclusions are drawn in section 4.

2. Data and methods

To estimate air–sea fluxes, the GTMBA observations are used as inputs for the Coupled Ocean–Atmosphere Response Experiment (COARE) 3.5 bulk flux algorithm (Edson et al. 2013). The GTMBA observations provide long-term, high-frequency measurements (10 min and hourly) of air–sea variables, including surface vector winds, surface air temperature and relative humidity, and water temperature near the surface (McPhaden et al. 2010). The 10-min measurements were averaged into hourly means before the flux estimate. For surface winds, we used vector averaging. More information about the GTMBA observations was provided in Yan et al. (2021, 2022).

The wind stress magnitude (τ) is estimated using the COARE 3.5 bulk flux algorithm (Edson et al. 2013):
τ=ρCd(|V|)|V|S=ρCd(|V|)|V|2fG,
where ρ is the near-surface air density, Cd(|V|) is the drag coefficient, V is the hourly vector-averaged surface winds, |V| is the vector-averaged wind speed ( |V|=u2+υ2, where u and υ are the zonal and meridional components of the hourly vector-averaged surface winds), S is the hourly scalar-averaged wind speed, and fG is the convective gustiness factor (fG = S/|V|, representing the ratio of the hourly scalar-averaged wind speed to the vector-averaged wind speed). The factor attempts to eliminate the estimation error when the vector-averaged winds approach zero, but the scalar-averaged wind speed is nonzero (Fairall et al. 1996). The drag coefficient is dependent on wind speed and atmospheric stability. Here, for convenience, the drag coefficient is expressed as a function of wind speed only. The true wind speed, rather than the relative wind speed, is used due to the lack of credible data on tropical surface currents with a diurnal jet (Masich et al. 2021).

Next, the wind stress magnitude induced by the mean daily wind vector is calculated using the COARE 3.5 algorithm [ τ_MDW=ρCd(|V¯|)|V¯|2fG, where V¯ is the mean daily wind vector]. The difference between τ and τ_MDW is defined as the wind stress magnitude induced by subdaily winds (τ_SDW). The annual climatology of τ_SDW is obtained by averaging all available samples to represent the effects of subdaily winds on the mean wind stress magnitude. Only buoys with more than two years of available data are selected in this study. Note that convective gustiness is also considered in the estimation of τ_MDW. Therefore, the nominal subdaily winds in this study represent winds at time scales between 1 day and 1 h.

To examine how subdaily winds influence the mean wind stress magnitude, the variables in Eq. (1) are divided into two components:
Cd(|V|)=Cd(|V|)¯+Cd(|V|)
|V|=|V|¯+|V|,
where an overbar denotes the daily mean and a prime symbol denotes a subdaily fluctuation. The contribution of subdaily winds to the daily mean wind stress magnitude ( τ_SDW¯) is written as
τ_SDW¯=τ¯τ_MDW¯=ρfGCd(|V|)¯(|V|¯2|V¯|2)meanwindspeedchangeterm+ρfG[Cd(|V|)¯Cd(|V¯|)¯]|V¯|2meandragcoefficientchangeterm+ρfGCd(|V|)|V|2¯windspeedvarianceterm+2ρfG|V|¯Cd(|V|)|V|¯covarianceterm
Here, the convective gustiness factor for estimating wind stress associated with the mean daily wind vector (τ_MDW) is assumed to be equal to that in the estimation of the hourly wind stress (τ). This assumption leads to a 3% increase in the total contribution of subdaily winds. The terms on the right-hand side of Eq. (3) represent the contributions from changes in the mean wind speed and drag coefficient, wind speed variance, and covariance between wind speed and drag coefficient due to subdaily winds. Here, these terms are defined as the mean wind speed change (MWSC) term, mean drag coefficient change (MDCC) term, wind speed variance (WSV) term, and covariance (CV) term. The annual climatologies of all the terms are calculated using the GTMBA observations to assess their corresponding effects on the mean wind stress magnitude.

Likewise, all the terms in Eq. (3) are estimated using the hourly ERA5 and MERRA2 reanalysis wind values to evaluate their performance. Note that the wind stresses estimated from ERA5 and MERRA2 are based on the COARE algorithm for the convenience of comparison with those derived from the GTMBA observations. However, wind stresses directly downloaded from these two reanalyses are estimated using their own algorithms. For example, the ERA5 air–sea fluxes are calculated using the ECMWF algorithm (ECMWF 2021). Therefore, the wind stress magnitudes of the reanalysis products in this study refer to our calculated wind stresses via the COARE algorithm.

Finally, we develop two new parameterization schemes to account for the effect of subdaily winds: one is time dependent, and the other is not. We use the gustiness approach, which has been widely applied to parameterize unresolved wind variability at different spatial and temporal scales (e.g., Bessac et al. 2019, 2021; Blein et al. 2020, 2022; Godfrey and Beljaars 1991; Hourdin et al. 2020; Redelsperger et al. 2000; Zeng et al. 2002). To establish the time-independent scheme, the climatological gustiness value (Gclim) is first calculated from the GTMBA observations per the methods of Cronin et al. (2006): Gclim=|V|¯2|V¯|2, where an angle bracket denotes a record-length average at each buoy site. The gustiness value (Gclim) was defined as the mesoscale gustiness, representing the impact of subdaily winds (at time scales between 1 day and 1 h) on wind speed (Cronin et al. 2006). It is nearly twice as large as the convective gustiness in the equatorial eastern Pacific Ocean (Cronin et al. 2006). Next, a simple linear function is built between the square of the gustiness value ( Gclim2) and the corresponding precipitation rate (Pclim). Here, we select Gclim2, rather than Gclim, to build the regression model due to the quadratic relation between wind stress magnitude and wind speed [Eq. (1)]. Likewise, the time-dependent scheme using monthly values is created based on the expression of the Michaels–Menten equation. We also apply the Michaels–Menten equation to build the time-independent scheme, and it performs as well as the linear model. For simplicity, only the linear regression model is given for the time-independent scheme in this study. The precipitation rate used here is the research-level product of the Integrated Multi-satellitE Retrievals for GPM (IMERG) Final Run.

3. Results

We start by examining the mean surface wind stress magnitude obtained from the GTMBA observations (Fig. 1a). Spatially, the wind stress magnitude is relatively small in the calm equatorial belt, especially in the equatorial Indian Ocean, western Pacific, and equatorial cold tongue region of the Pacific and Atlantic, reaching amplitudes of less than 4 × 10−2 N m−2. In contrast, the wind stress magnitude is larger in areas farther from the equator, such as the trade winds region, where it reaches amplitudes of over 8 × 10−2 N m−2. Subdaily winds make a positive contribution to the mean wind stress magnitude (Figs. 1b,c), with an average contribution of 6.5 × 10−3 N m−2 (12.4%) in the tropics. The contribution is particularly large over the intertropical convergence zone (ITCZ), a region characterized by strong precipitation, and it increases to 8.4 × 10−3 N m−2 (17.4%) in regions with a surface precipitation rate greater than 5 mm day−1. The maximum percentage contribution from subdaily winds reaches 28.5% in the equatorial western Pacific, where the mean wind stress magnitude is small. The estimates from buoy observations suggest that subdaily winds have a nonnegligible effect on the mean surface wind stress magnitude over the tropics, especially over the ITCZ.

Fig. 1.
Fig. 1.

(a) Mean surface wind stress magnitude estimated from GTMBA observations. (b),(c) Contribution from subdaily winds to the mean wind stress magnitude and the corresponding percentage. Gray arrows and contours denote the satellite-derived mean surface wind stress in (a) and precipitation rate (mm day−1) in (b) and (c), respectively. The satellite-derived mean wind stress was calculated from the Scatterometer Climatology of Ocean Winds (SCOW; Risien and Chelton 2008).

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

The contribution from subdaily winds to the mean wind stress magnitude arises from four aspects, which are associated with the mean wind speed change, mean drag coefficient change, wind speed variance, and covariance between wind speed and the drag coefficient [see Eq. (3)]. The MWSC, MDCC, WSV, and CV terms are all estimated from the GTMBA observations and amount to 3.5 × 10−3, 0.8 × 10−3, 2.1 × 10−3, and 0.4 × 10−3 N m−2 (51.9%, 12.3%, 31.2%, and 7.8%), respectively, on average, using all buoy sites (Fig. 2). The two largest contributors are associated with the MWSC and WSV terms, which account for approximately 80% of the overall contribution. The other two terms, which are related to the variable drag coefficient, make an approximately 20% total contribution. Moreover, the MWSC and WSV terms exhibit significant spatial variations, with large values located over the ITCZ, while the spatial distributions of the other two terms are relatively uniform (see also Fig. S1 in the online supplemental material). The enhanced MWSC and WSV terms over the ITCZ are thus the main contributors to the increased wind stress magnitude in these regions. The spatial correlation coefficient between the MWSC (WSV) term and the total contribution is 0.98 (0.93). This result implies that the contribution from subdaily winds to the mean wind stress magnitude largely results from the MWSC and WSV terms, which are related to the change in the mean wind speed due to subdaily winds and the variance in subdaily wind speed variability.

Fig. 2.
Fig. 2.

All the terms in Eq. (3) estimated from the GTMBA observations: The total contribution (black curve, representing the values at all buoy sites as shown in Fig. 1b but ordered from west to the east and from south to north; the order is listed in Table S1; black dots on the x axis denote the buoys, and the corresponding labels are their longitudes), mean wind speed change term (red curve), mean drag coefficient change term (cyan curve), wind speed variance term (blue curve), and covariance term (green curve).

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

Having calculated the contribution of subdaily winds to the mean wind stress, we next investigate the total change in wind speed, averaged over the full record, that results from the subdaily winds ( |V|¯|V¯|, Fig. 3a). The mean wind speed increases over the entire tropical region, with the increase being enhanced over the ITCZ. Specifically, the average value of the increase is 0.29 m s−1 in the tropics and rises to 0.41 m s−1 over the ITCZ (surface precipitation rate > 5 mm day−1). The variance in the subdaily wind speed variability ( |V|2) exhibits a similar spatial structure, with large values located over the ITCZ (1.35 m2 s−2, Fig. 3b). The enhanced mean wind speed change (subdaily wind speed variance) over the ITCZ leads to the increase in the MWSC (WSV) term in these regions, indicated by a high correlation coefficient of 0.95 (0.99) between them.

Fig. 3.
Fig. 3.

(a) Change in the mean wind speed due to subdaily winds ( |V|¯|V¯|). (b) Variance in subdaily wind speed variability ( |V|2). (c) Kinetic energy of subdaily winds [ ρ(u2+υ2)/2]. Gray contours in (a)–(c) denote the satellite-derived mean surface precipitation rate (mm day−1).

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

Considering that the mean wind speed change and subdaily wind speed variance are both dependent on the strength of subdaily winds, the kinetic energy of subdaily winds is examined next. The kinetic energy is calculated as ρ(u2+υ2)/2, where u′ and υ′ are the zonal and meridional components of subdaily winds. It is evident that the subdaily winds are strong over the ITCZ, which is characterized by high kinetic energy (Fig. 3c). The strong subdaily winds over the ITCZ induce large mean wind speed changes and subdaily wind speed variances, resulting in large subdaily wind contributions. Additionally, there is a linear relationship between the subdaily wind contribution (〈τ_SDW〉, units: 10−3 N m−2) and the kinetic energy of subdaily winds (〈KE_SDW〉, units: kg m−1 s−2):
τ_SDW=3.554KE_SDW+0.58,
where the explained variance (determination coefficient) is 0.93 and the root-mean-square error (RMSE) is 0.57 (Fig. 4). The linear relationship can be derived from the formula for the mean wind stress magnitude from subdaily winds, assuming constant convective gustiness and drag coefficient:
τ_SDW¯=ρCd(u2+υ2¯)fGρCd(u¯2+υ¯2)fG=2fGCdKE_SDW¯.
The linear expression can also be obtained from Eq. (3) (not shown). Thus, the contribution of subdaily winds to the mean wind stress magnitude primarily depends on the kinetic energy (strength) of the subdaily winds. The high kinetic energy over the ITCZ may be associated with the frequently occurring mesoscale convective system there (e.g., Yuan and Houze 2010; Feng et al. 2021), which can enhance the subdaily winds by generating cold pools (density currents) at the sea surface (e.g., Garg et al. 2021; Joseph et al. 2021). Future analyses quantifying the physical processes responsible for subdaily winds are desirable.
Fig. 4.
Fig. 4.

Linear relationship between the mean wind stress magnitude from subdaily winds and the kinetic energy of subdaily winds.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

To evaluate the role of the subdaily winds in the most recent atmospheric reanalysis products, the contribution from the subdaily winds to the mean wind stress magnitude is calculated from ERA5 and MERRA2 data. Notably, both products underestimate the contribution of subdaily winds (Fig. 5a). The mean absolute percentage error (MAE) in the tropics is 51% and 63% for ERA5 and MERRA2, respectively, where MAE = |(reanalysis − observation)/observation| × 100%. The value of the error in the contribution corresponds to the errors for the MWSC and WSV terms (the two largest contributors), which are underestimated by 50% and 49% (61% and 64%) in the ERA5 (MERRA2) product (Figs. 5b,c). The underestimated MWSC and WSV terms in these two products arise from the underestimations of the mean wind speed change due to subdaily winds (the MAEs are 47% and 59%) and the variance in subdaily wind speed variability (the MAEs are 56% and 71%), which are associated with the underestimation of the (strength) kinetic energy of subdaily winds (the MAEs are 59% and 72%, Fig. 5d).

Fig. 5.
Fig. 5.

(a) Contribution from subdaily winds to the mean wind stress magnitude at all buoy sites. (b) Change in the mean wind speed due to subdaily winds (dashed curve) and the associated contribution (solid curve). (c) Variance in subdaily wind speed variability (dashed curve) and the associated contribution (solid curve). (d) Kinetic energy of subdaily winds. Black, red, and blue curves denote values estimated from the GTMBA observations, ERA5 product, and MERRA2 product, respectively. The x axes in (a)–(d) are the same as those in Fig. 2.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

To assess in more detail the extent to which ERA5 and MERRA2 products underestimate the effects of subdaily winds relative to the GTMBA observations, we next examine the kinetic energy spectrum of winds. The kinetic energy spectrum is calculated by averaging the zonal and meridional wind power spectra, which are obtained using Welch’s overlapped segment averaging estimator with a 60-day Hamming window and 0% overlap. The results include peaks at diurnal (24-h) and semidiurnal (12-h) time scales (Fig. 6), which have been reported in previous studies (Deser and Smith 1998; Ueyama and Deser 2008). The corresponding diurnal and semidiurnal winds are induced by atmospheric tides and/or land–sea breezes (Christophersen et al. 2020; Takahashi 2012; Ueyama and Deser 2008). ERA5 and MERRA2 are capable of reproducing diurnal (84% and 77%) and semidiurnal (111% and 91%) winds. However, at other subdaily time scales, these two reanalysis products underestimate wind variability, and the estimation performance decreases as the time scale becomes shorter. For example, the MAE of kinetic energy is 48% (55%) between the 24- and 23-h time scales and increases to approximately 94% (100%) between the 3- and 2-h time scales in the ERA5 (MERRA2) product. The underestimation of the kinetic energy at these time scales results in underestimated subdaily winds in the ERA5 and MERRA2 products. A reason for the underestimated kinetic energy in atmospheric reanalyses is that, in numerical models, the energy in the highest resolved wavenumbers is damped due to the nature of numerical solutions and parameterizations (Abdalla et al. 2013; Bolgiani et al. 2022; Klaver et al. 2020; Skamarock 2004). The model’s effective resolution is typically 7 times the model grid spacing (Skamarock 2004). Recently, Bolgiani et al. (2022) showed that ERA5 has an effective resolution from 1300 to 1100 km in the tropical band and hence does not adequately capture mesoscale variability. In contrast, the buoy observations at fixed points reflect the contributions due to mesoscale spatial variability.

Fig. 6.
Fig. 6.

Kinetic energy spectrum of surface winds (variance-preserving plots, solid curve) and the ratio between the reanalysis estimates and observations (percentage, dashed curve). Black, red, and blue curves denote values estimated from the GTMBA observations, ERA5 product, and MERRA2 product, respectively.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

Importantly, the underestimated contribution from the ERA5 (MERRA2) subdaily winds results in a 7% (8%) error in the mean wind stress magnitude, on average, in the tropics (black curve in Fig. 7). The error increases to 12% (14%) in the tropical western Pacific. The total underestimation of the mean wind stress magnitude is 11% (18%) in the ERA5 (MERRA2) product based on GTMBA observations (red curve in Fig. 7). The underestimation of the subdaily wind contribution accounts for 60% (44%) of the total error in the mean ERA5 (MERRA2) wind stress magnitude. This percentage ratio rises to 73% (88%) in the tropical western Pacific, implying that accurate simulations of subdaily winds in ERA5 and MERRA2 will greatly improve estimates of the mean wind stress magnitude in the tropics, especially in the tropical western Pacific.

Fig. 7.
Fig. 7.

Percentage error (red curve) of the mean wind stress magnitude in the (a) ERA5 and (b) MERRA2 products based on the GTMBA observations. The black curve denotes the percentage relative error associated with inaccurate subdaily winds. The x axes are the same as those in Fig. 2.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

The preceding analysis shows that two state-of-the-art reanalysis products largely underestimate the impact of subdaily winds, suggesting potential biases in ocean numerical simulations. We therefore develop a gustiness parameterization scheme to properly represent the effect of subdaily winds. Two new schemes are proposed in this study: one is time dependent, and the other is not. Due to the close relation between flux enhancement and precipitation (Fig. 1b), these two schemes are parameterized with the precipitation rate. The detailed procedures can be seen in section 2. The linear regression model for the time-independent scheme is expressed as
Gclim2=0.37Pclim+1.78,
where the explained variance is 0.93 and the RMSE is 0.37 (Fig. 8). Applying the new gustiness correction to account for subdaily winds [ τ_Gclim=ρCd(|V¯|)Gclim2] increases the mean wind stress magnitude by 9.6% at the GTMBA sites, which explains 81.0% of the contribution of subdaily winds. Compared with the climatological parameterization scheme based on SST by Praveen Kumar et al. (2012, 2013), the new scheme performs better. An ∼20% underestimation of the subdaily wind contribution remains due to the use of invariable gustiness in the new scheme.
Fig. 8.
Fig. 8.

Regression model between the square of the gustiness value (G2) and precipitation rate (P) for the (left) time-independent and (right) time-dependent schemes.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

To overcome the limitation of the time-independent scheme, a time-dependent scheme using monthly gustiness and precipitation (Gmon and Pmon) is built based on the following regression model:
Gmon2=aPmonPmon+b+c,
where a = 6.88, b = 8.05, and c = 1.64. As shown in Fig. 8, the constant c represents the gustiness correction when there is no precipitation, a represents the maximum correction associated with precipitation, and b represents the precipitation rate corresponding to the half-maximum value of the precipitation-induced gustiness (a/2). The explained variance of the model is 0.68, and the RMSE is 1.12. Applying the time-dependent scheme explains 78.9% of the contribution of subdaily winds to the mean wind stress magnitude. The scheme has almost the same performance as the time-independent scheme. Including temporal variability does not improve the scheme performance because the relation between gustiness and precipitation becomes weaker at this time scale, as indicated by the explained variance of Eq. (7). We also tested a parameterization scheme based on wind speed, but the performance was lower than that of the scheme with precipitation. Thus, we have not included it here.

4. Conclusions and discussion

Previous studies suggested that subdaily winds may contribute significantly to the mean surface wind stress, but direct estimates over the open ocean are lacking, and to date, there is no quantification of this contribution on a global scale. Here, the effects of subdaily winds on the mean surface wind stress magnitude were systematically quantified in the entirety of the tropical oceans using high-frequency measurements of air–sea variables from the GTMBA. We found that subdaily winds contribute 12.4% to the mean wind stress magnitude on average, and the contribution is large over the ITCZ. The magnitude of the contribution is primarily determined by the kinetic energy of subdaily winds. We also estimated the effect of subdaily winds on the mean wind stress vector (Fig. 9). The inclusion of subdaily winds clearly significantly enhances the mean stress vector. The zonal and meridional components increase by 3.2% and 9.4%, respectively. This suggests that including subdaily winds in climate models would improve their simulations of the equatorial current regime and the resultant simulated ocean variability.

Fig. 9.
Fig. 9.

Effect of subdaily winds on wind stress vector (red arrows). Black arrows denote the long-term mean wind stress vector. Red arrows are amplified by 10.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

Thus, our study highlights the importance of subdaily winds for the total air–sea momentum flux in the tropics. Two state-of-the-art and widely used reanalyses, ERA5 and MERRA2, largely underestimate this effect by 51% and 63%, leading to 7% and 8% underestimations in the mean wind stress magnitude. Considering the importance of accurate wind stress estimates for simulating air–sea interactions, the identified biases are expected to have far-reaching implications for the simulated ocean and climate variability in models. Notably, the wind stresses derived from ERA5 and MERRA2 use the COARE algorithm, rather than their own bulk algorithms. To investigate the sensitivity to different algorithms, we also employ the ERA5 own algorithm (ECMWF algorithm) to estimate the effect of ERA5 subdaily winds. The estimated contribution of subdaily winds using the ECMWF algorithm is almost the same as that using the COARE algorithm (Fig. 10a). This finding suggests that our main conclusions are insensitive to the bulk algorithm. It is worth mentioning that the ECMWF algorithm slightly overestimates the contribution of the ERA5 subdaily winds with respect to the COARE algorithm (Fig. 10). The bulk formula [Eq. (1)] shows that the discrepancy between different bulk algorithms for estimating wind stress mainly relies on the parameterizations for drag coefficient and convective gustiness. The parameterizations for convective gustiness are the same in the COARE and ECWMF algorithms (Brodeau et al. 2017; also see black curve in Fig. 10b). Therefore, the discrepancy results from the different parameterizations for drag coefficient typically differ by approximately 15% between the bulk algorithms (blue curve in Fig. 10b; Bonino et al. 2022; Brodeau et al. 2017).

Fig. 10.
Fig. 10.

(a) Contribution from subdaily winds to the mean wind stress magnitude estimated from the GTMBA observation using the COARE algorithm (black curve) and from the ERA5 reanalysis using the COARE (red curve) and ECMWF (blue curve) algorithms. (b) Mean absolute percentage errors of the subdaily wind contribution (red curve), convective gustiness (black curve) and drag coefficient (blue curve) calculated using the ECMWF algorithm with respect to those calculated using the COARE algorithm.

Citation: Journal of Physical Oceanography 53, 3; 10.1175/JPO-D-22-0156.1

To improve the representation of subdaily winds, two parameterization schemes related to precipitation were developed in this study: one is time dependent, and the other is not. These two new schemes have almost the same performance regarding the effect of subdaily winds, explaining ∼80% of the mean wind stress magnitude associated with subdaily winds. The inclusion of these parameterization schemes in climate models is expected to substantially improve simulations of large-scale climate variability. Note that there still exists an ∼20% underestimation of the subdaily wind contribution when using these two schemes. Further improvements may potentially be achieved by future analyses that develop schemes based on different variables and techniques, for instance, using machine learning approaches. However, these analyses are beyond the scope of this study.

Acknowledgments.

This work is supported by the National Natural Science Foundation of China (41976028, 42122040, and 42076016).

Data availability statement.

The GTMBA measurements are available at https://www.pmel.noaa.gov/tao/drupal/disdel/, the ERA5 reanalysis at https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5, the MERRA2 reanalysis at https://disc.gsfc.nasa.gov/datasets?project=MERRA-2, and the IMERG products at https://gpm.nasa.gov/data/directory.

REFERENCES

  • Abdalla, S., L. Isaksen, P. Janssen, and N. Wedi, 2013: Effective spectral resolution of ECMWF atmospheric forecast models. ECMWF Newsletter, No. 137, ECMWF, Reading, United Kingdom, 19–22, https://www.ecmwf.int/node/17358.

  • Bessac, J., A. H. Monahan, H. M. Christensen, and N. Weitzel, 2019: Stochastic parameterization of subgrid-scale velocity enhancement of sea surface fluxes. Mon. Wea. Rev., 147, 14471469, https://doi.org/10.1175/MWR-D-18-0384.1.

    • Search Google Scholar
    • Export Citation
  • Bessac, J., H. M. Christensen, K. Endo, A. H. Monahan, and N. Weitzel, 2021: Scale-aware space-time stochastic parameterization of subgrid-scale velocity enhancement of sea surface fluxes. J. Adv. Model. Earth Syst., 13, e2020MS002367, https://doi.org/10.1029/2020MS002367.

    • Search Google Scholar
    • Export Citation
  • Blein, S., R. Roehrig, A. Voldoire, and G. Faure, 2020: Meso-scale contribution to air-sea turbulent fluxes at GCM scale. Quart. J. Roy. Meteor. Soc., 146, 24662495, https://doi.org/10.1002/qj.3804.

    • Search Google Scholar
    • Export Citation
  • Blein, S., R. Roehrig, and A. Voldoire, 2022: Parameterizing the meso-scale enhancement of oceanic surface turbulent fluxes: A physical-statistical approach. Quart. J. Roy. Meteor. Soc., 148, 16831708, https://doi.org/10.1002/qj.4273.

    • Search Google Scholar
    • Export Citation
  • Bolgiani, P., and Coauthors, 2022: Wind kinetic energy climatology and effective resolution for the ERA5 reanalysis. Climate Dyn., 59, 737752, https://doi.org/10.1007/s00382-022-06154-y.

    • Search Google Scholar
    • Export Citation
  • Bonino, G., D. Iovino, L. Brodeau, and S. Masina, 2022: The bulk parameterizations of turbulent air–sea fluxes in NEMO4: The origin of sea surface temperature differences in a global model study. Geosci. Model Dev., 15, 68736889, https://doi.org/10.5194/gmd-15-6873-2022.

    • Search Google Scholar
    • Export Citation
  • Brodeau, L., B. Barnier, S. K. Gulev, and C. Woods, 2017: Climatologically significant effects of some approximations in the bulk parameterizations of turbulent air–sea fluxes. J. Phys. Oceanogr., 47, 528, https://doi.org/10.1175/JPO-D-16-0169.1.

    • Search Google Scholar
    • Export Citation
  • Christophersen, J. A., G. R. Foltz, and R. C. Perez, 2020: Surface expressions of atmospheric thermal tides in the tropical Atlantic and their impact on open-ocean precipitation. J. Geophys. Res. Atmos., 125, e2019JD031997, https://doi.org/10.1029/2019JD031997.

    • Search Google Scholar
    • Export Citation
  • Cronin, M. F., C. W. Fairall, and M. J. McPhaden, 2006: An assessment of buoy-derived and numerical weather prediction surface heat fluxes in the tropical Pacific. J. Geophys. Res., 111, C06038, https://doi.org/10.1029/2005JC003324.

    • Search Google Scholar
    • Export Citation
  • Deser, C., and C. A. Smith, 1998: Diurnal and semidiurnal variations of the surface wind field over the tropical Pacific Ocean. J. Climate, 11, 17301748, https://doi.org/10.1175/1520-0442(1998)011<1730:DASVOT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Duteil, O., 2019: Wind synoptic activity increases oxygen levels in the tropical Pacific Ocean. Geophys. Res. Lett., 46, 27152725, https://doi.org/10.1029/2018GL081041.

    • Search Google Scholar
    • Export Citation
  • ECMWF, 2021: IFS Documentation Cy47r3. Part IV: Physical processes. ECMWF, 232 pp., https://www.ecmwf.int/sites/default/files/elibrary/2021/20198-ifs-documentation-cy47r3-part-vi-physical-processes.pdf.

  • Edson, J. B., and Coauthors, 2013: On the exchange of momentum over the open ocean. J. Phys. Oceanogr., 43, 15891610, https://doi.org/10.1175/JPO-D-12-0173.1.

    • Search Google Scholar
    • Export Citation
  • Eymard, L., A. Weill, D. Bourras, C. Guérin, P. Le Borgne, and J. M. Lefévre, 2003: Use of ship mean data for validating model and satellite flux fields during the FETCH experiment. J. Geophys. Res., 108, 8060, https://doi.org/10.1029/2001JC001207.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., E. F. Bradley, D. P. Rogers, J. B. Edson, and G. S. Young, 1996: Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere coupled-ocean atmosphere response experiment. J. Geophys. Res., 101, 37473764, https://doi.org/10.1029/95JC03205.

    • Search Google Scholar
    • Export Citation
  • Feng, Z., and Coauthors, 2021: A global high-resolution mesoscale convective system database using satellite-derived cloud tops, surface precipitation, and tracking. J. Geophys. Res. Atmos., 126, e2020JD034202, https://doi.org/10.1029/2020JD034202.

    • Search Google Scholar
    • Export Citation
  • Garg, P., S. W. Nesbitt, T. J. Lang, and G. Priftis, 2021: Diurnal cycle of tropical oceanic mesoscale cold pools. J. Climate, 34, 93059326, https://doi.org/10.1175/JCLI-D-20-0909.1.

    • Search Google Scholar
    • Export Citation
  • Gelaro, R., and Coauthors, 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). J. Climate, 30, 54195454, https://doi.org/10.1175/JCLI-D-16-0758.1.

    • Search Google Scholar
    • Export Citation
  • Giglio, D., S. T. Gille, A. C. Subramanian, and S. Nguyen, 2017: The role of wind gusts in upper ocean diurnal variability. J. Geophys. Res. Oceans, 122, 77517764, https://doi.org/10.1002/2017JC012794.

    • Search Google Scholar
    • Export Citation
  • Giglio, D., S. T. Gille, B. D. Cornuelle, A. C. Subramanian, F. J. Turk, S. Hristova-Veleva, and D. Northcott, 2022: Annual modulation of diurnal winds in the tropical oceans. Remote Sens., 14, 459, https://doi.org/10.3390/rs14030459.

    • Search Google Scholar
    • Export Citation
  • Gill, A. E., 1982: Atmosphere-Ocean Dynamics. Academic Press, 662 pp.

  • Godfrey, J. S., and A. C. M. Beljaars, 1991: On the turbulent fluxes of buoyancy, heat and moisture at the air-sea interface at low wind speeds. J. Geophys. Res., 96, 22 04322 048, https://doi.org/10.1029/91JC02015.

    • Search Google Scholar
    • Export Citation
  • Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 19992049, https://doi.org/10.1002/qj.3803.

    • Search Google Scholar
    • Export Citation
  • Hourdin, F., and Coauthors, 2020: LMDZ6A: The atmospheric component of the IPSL climate model with improved and better tuned physics. J. Adv. Model. Earth Syst., 12, e2019MS001892, https://doi.org/10.1029/2019MS001892.

    • Search Google Scholar
    • Export Citation
  • Joseph, J., M. S. Girishkumar, M. J. McPhaden, and E. P. R. Rao, 2021: Diurnal variability of atmospheric cold pool events and associated air-sea interactions in the Bay of Bengal during the summer monsoon. Climate Dyn., 56, 837853, https://doi.org/10.1007/s00382-020-05506-w.

    • Search Google Scholar
    • Export Citation
  • Klaver, R., R. Haarsma, P. L. Vidale, and W. Hazeleger, 2020: Effective resolution in high resolution global atmospheric models for climate studies. Atmos. Sci. Lett., 21, e952, https://doi.org/10.1002/asl.952.

    • Search Google Scholar
    • Export Citation
  • Ledvina, D. V., G. S. Young, R. A. Miller, and C. W. Fairall, 1993: The effect of averaging on bulk estimates of heat and momentum fluxes for the tropical western Pacific Ocean. J. Geophys. Res., 98, 20 21120 217, https://doi.org/10.1029/93JC01856.

    • Search Google Scholar
    • Export Citation
  • Lee, T., and W. T. Liu, 2005: Effects of high-frequency wind sampling on simulated mixed layer depth and upper ocean temperature. J. Geophys. Res., 110, C05002, https://doi.org/10.1029/2004JC002746.

    • Search Google Scholar
    • Export Citation
  • Lin, X., X. Zhai, Z. Wang, and D. R. Munday, 2018: Mean, variability and trend of Southern Ocean wind stress: Role of wind fluctuations. J. Climate, 31, 35573573, https://doi.org/10.1175/JCLI-D-17-0481.1.

    • Search Google Scholar
    • Export Citation
  • Lin, X., X. Zhai, Z. Wang, and D. R. Munday, 2020: Southern Ocean wind stress in CMIP5 models: Role of wind fluctuations. J. Climate, 33, 12091226, https://doi.org/10.1175/JCLI-D-19-0466.1.

    • Search Google Scholar
    • Export Citation
  • Masich, J., W. S. Kessler, M. F. Cronin, and K. R. Grissom, 2021: Diurnal cycles of near-surface currents across the tropical Pacific. J. Geophys. Res. Oceans, 126, e2020JC016982, https://doi.org/10.1029/2020JC016982.

    • Search Google Scholar
    • Export Citation
  • McPhaden, M. J., and Coauthors, 2010: The global tropical moored buoy array. Proc. OceanObs’09: Sustained Ocean Observations and Information for Society, Venice, Italy, ESA Publication, WPP-306, https://doi.org/10.5270/OceanObs09.cwp.61.

  • Ouni, M., A. Mehra, and A. Harzallah, 2021: The effect of the temporal resolution of the wind forcing on a central Mediterranean Sea model. Dyn. Atmos. Oceans, 96, 101262, https://doi.org/10.1016/j.dynatmoce.2021.101262.

    • Search Google Scholar
    • Export Citation
  • Ponte, R. M., and R. D. Rosen, 2004: Nonlinear effects of variable winds on ocean stress climatologies. J. Climate, 17, 12831293, https://doi.org/10.1175/1520-0442(2004)017<1283:NEOVWO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Praveen Kumar, B., J. Vialard, M. Lengaigne, V. S. N. Murty, and M. J. McPhaden, 2012: Tropflux: Air-sea fluxes for the global tropical oceans-description and evaluation. Climate Dyn., 38, 15211543, https://doi.org/10.1007/s00382-011-1115-0.

    • Search Google Scholar
    • Export Citation
  • Praveen Kumar, B., J. Vialard, M. Lengaigne, V. Murty, M. J. McPhaden, M. F. Cronin, F. Pinsard, and K. Gopala Reddy, 2013: Tropflux wind stresses over the tropical oceans: Evaluation and comparison with other products. Climate Dyn., 40, 20492071, https://doi.org/10.1007/s00382-012-1455-4.

    • Search Google Scholar
    • Export Citation
  • Redelsperger, J.-L., F. Guichard, and S. Mondon, 2000: A parameterization of mesoscale enhancement of surface fluxes for large-scale models. J. Climate, 13, 402421, https://doi.org/10.1175/1520-0442(2000)013<0402:APOMEO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Risien, C. M., and D. B. Chelton, 2008: A global climatology of surface wind and wind stress fields from eight years of QuikSCAT scatterometer data. J. Phys. Oceanogr., 38, 23792413, https://doi.org/10.1175/2008JPO3881.1.

    • Search Google Scholar
    • Export Citation
  • Serra, Y. L., M. F. Cronin, and G. N. Kiladis, 2007: Sub-seasonal variance of surface meteorological parameters in buoy observations and reanalyses. Geophys. Res. Lett., 34, L12708, https://doi.org/10.1029/2007GL029506.

    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., 2004: Evaluating mesoscale NWP models using kinetic energy spectra. Mon. Wea. Rev., 132, 30193032, https://doi.org/10.1175/MWR2830.1.

    • Search Google Scholar
    • Export Citation
  • Takahashi, K., 2012: Thermotidal and land-heating forcing of the diurnal cycle of oceanic surface winds in the eastern tropical Pacific. Geophys. Res. Lett., 39, L04805, https://doi.org/10.1029/2011GL050692.

    • Search Google Scholar
    • Export Citation
  • Thushara, V., and P. N. Vinayachandran, 2014: Impact of diurnal forcing on intraseasonal sea surface temperature oscillations in the Bay of Bengal. J. Geophys. Res. Oceans, 119, 82218241, https://doi.org/10.1002/2013JC009746.

    • Search Google Scholar
    • Export Citation
  • Turk, F. J., S. Hristovaveleva, and D. Giglio, 2021: Examination of the daily cycle wind vector modes of variability from the constellation of microwave scatterometers and radiometers. Remote Sens., 13, 141, https://doi.org/10.3390/rs13010141.

    • Search Google Scholar
    • Export Citation
  • Ueyama, R., and C. Deser, 2008: A climatology of diurnal and semidiurnal surface wind variations over the tropical Pacific Ocean based on the tropical atmosphere ocean moored buoy array. J. Climate, 21, 593607, https://doi.org/10.1175/JCLI1666.1.

    • Search Google Scholar
    • Export Citation
  • Wu, Y., X. Zhai, and Z. Wang, 2016: Impact of synoptic atmospheric forcing on the mean ocean circulation. J. Climate, 29, 57095724, https://doi.org/10.1175/JCLI-D-15-0819.1.

    • Search Google Scholar
    • Export Citation
  • Wu, Y., Z. Wang, C. Liu, and X. Lin, 2020: Impacts of high-frequency atmospheric forcing on Southern Ocean circulation and Antarctic Sea ice. Adv. Atmos. Sci., 37, 515531, https://doi.org/10.1007/s00376-020-9203-x.

    • Search Google Scholar
    • Export Citation
  • Wunsch, C., and R. Ferrari, 2004: Vertical mixing, energy, and the general circulation of the oceans. Annu. Rev. Fluid Mech., 36, 281314, https://doi.org/10.1146/annurev.fluid.36.050802.122121.

    • Search Google Scholar
    • Export Citation
  • Yan, Y., L. Zhang, X. Song, G. Wang, and C. Chen, 2021: Diurnal variation in surface latent heat flux and the effect of diurnal variability on the climatological latent heat flux over the tropical oceans. J. Phys. Oceanogr., 51, 34013415, https://doi.org/10.1175/JPO-D-21-0128.1.

    • Search Google Scholar
    • Export Citation
  • Yan, Y., X. Song, G. Wang, and C. Chen, 2022: Importance of high-frequency (≤30-day) wind variability to the annual climatology of the surface latent heat flux inferred from the global tropical moored buoy array. J. Geophys. Res. Oceans, 127, e2021JC018094, https://doi.org/10.1029/2021JC018094.

    • Search Google Scholar
    • Export Citation
  • Yu, Y., and Coauthors, 2020: Importance of diurnal forcing on the summer salinity variability in the East China Sea. J. Phys. Oceanogr., 50, 633653, https://doi.org/10.1175/JPO-D-19-0200.1.

    • Search Google Scholar
    • Export Citation
  • Yuan, J., and R. A. Houze Jr., 2010: Global variability of mesoscale convective system anvil structure from A-Train satellite data. J. Climate, 23, 58645888, https://doi.org/10.1175/2010JCLI3671.1.

    • Search Google Scholar
    • Export Citation
  • Zeng, X., Q. Zhang, D. Johnson, and W.-K. Tao, 2002: Parameterization of wind gustiness for the computation of ocean surface fluxes at different spatial scales. Mon. Wea. Rev., 130, 21252133, https://doi.org/10.1175/1520-0493(2002)130<2125:POWGFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., 2013: On the wind mechanical forcing of the ocean general circulation. J. Geophys. Res. Oceans, 118, 65616577, https://doi.org/10.1002/2013JC009086.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., and C. Wunsch, 2013: On the variability of wind power input to the oceans with a focus on the subpolar North Atlantic. J. Climate, 26, 38923903, https://doi.org/10.1175/JCLI-D-12-00472.1.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., H. L. Johnson, D. P. Marshall, and C. Wunsch, 2012: On the wind power input to the ocean general circulation. J. Phys. Oceanogr., 42, 13571365, https://doi.org/10.1175/JPO-D-12-09.1.

    • Search Google Scholar
    • Export Citation
  • Zhou, S., X. Zhai, and I. A. Renfrew, 2018: The impact of high-frequency weather systems on SST and surface mixed layer in the central Arabian Sea. J. Geophys. Res. Oceans, 123, 10911104, https://doi.org/10.1002/2017JC013609.

    • Search Google Scholar
    • Export Citation

Supplementary Materials

Save
  • Abdalla, S., L. Isaksen, P. Janssen, and N. Wedi, 2013: Effective spectral resolution of ECMWF atmospheric forecast models. ECMWF Newsletter, No. 137, ECMWF, Reading, United Kingdom, 19–22, https://www.ecmwf.int/node/17358.

  • Bessac, J., A. H. Monahan, H. M. Christensen, and N. Weitzel, 2019: Stochastic parameterization of subgrid-scale velocity enhancement of sea surface fluxes. Mon. Wea. Rev., 147, 14471469, https://doi.org/10.1175/MWR-D-18-0384.1.

    • Search Google Scholar
    • Export Citation
  • Bessac, J., H. M. Christensen, K. Endo, A. H. Monahan, and N. Weitzel, 2021: Scale-aware space-time stochastic parameterization of subgrid-scale velocity enhancement of sea surface fluxes. J. Adv. Model. Earth Syst., 13, e2020MS002367, https://doi.org/10.1029/2020MS002367.

    • Search Google Scholar
    • Export Citation
  • Blein, S., R. Roehrig, A. Voldoire, and G. Faure, 2020: Meso-scale contribution to air-sea turbulent fluxes at GCM scale. Quart. J. Roy. Meteor. Soc., 146, 24662495, https://doi.org/10.1002/qj.3804.

    • Search Google Scholar
    • Export Citation
  • Blein, S., R. Roehrig, and A. Voldoire, 2022: Parameterizing the meso-scale enhancement of oceanic surface turbulent fluxes: A physical-statistical approach. Quart. J. Roy. Meteor. Soc., 148, 16831708, https://doi.org/10.1002/qj.4273.

    • Search Google Scholar
    • Export Citation
  • Bolgiani, P., and Coauthors, 2022: Wind kinetic energy climatology and effective resolution for the ERA5 reanalysis. Climate Dyn., 59, 737752, https://doi.org/10.1007/s00382-022-06154-y.

    • Search Google Scholar
    • Export Citation
  • Bonino, G., D. Iovino, L. Brodeau, and S. Masina, 2022: The bulk parameterizations of turbulent air–sea fluxes in NEMO4: The origin of sea surface temperature differences in a global model study. Geosci. Model Dev., 15, 68736889, https://doi.org/10.5194/gmd-15-6873-2022.

    • Search Google Scholar
    • Export Citation
  • Brodeau, L., B. Barnier, S. K. Gulev, and C. Woods, 2017: Climatologically significant effects of some approximations in the bulk parameterizations of turbulent air–sea fluxes. J. Phys. Oceanogr., 47, 528, https://doi.org/10.1175/JPO-D-16-0169.1.

    • Search Google Scholar
    • Export Citation
  • Christophersen, J. A., G. R. Foltz, and R. C. Perez, 2020: Surface expressions of atmospheric thermal tides in the tropical Atlantic and their impact on open-ocean precipitation. J. Geophys. Res. Atmos., 125, e2019JD031997, https://doi.org/10.1029/2019JD031997.

    • Search Google Scholar
    • Export Citation
  • Cronin, M. F., C. W. Fairall, and M. J. McPhaden, 2006: An assessment of buoy-derived and numerical weather prediction surface heat fluxes in the tropical Pacific. J. Geophys. Res., 111, C06038, https://doi.org/10.1029/2005JC003324.

    • Search Google Scholar
    • Export Citation
  • Deser, C., and C. A. Smith, 1998: Diurnal and semidiurnal variations of the surface wind field over the tropical Pacific Ocean. J. Climate, 11, 17301748, https://doi.org/10.1175/1520-0442(1998)011<1730:DASVOT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Duteil, O., 2019: Wind synoptic activity increases oxygen levels in the tropical Pacific Ocean. Geophys. Res. Lett., 46, 27152725, https://doi.org/10.1029/2018GL081041.

    • Search Google Scholar
    • Export Citation
  • ECMWF, 2021: IFS Documentation Cy47r3. Part IV: Physical processes. ECMWF, 232 pp., https://www.ecmwf.int/sites/default/files/elibrary/2021/20198-ifs-documentation-cy47r3-part-vi-physical-processes.pdf.

  • Edson, J. B., and Coauthors, 2013: On the exchange of momentum over the open ocean. J. Phys. Oceanogr., 43, 15891610, https://doi.org/10.1175/JPO-D-12-0173.1.

    • Search Google Scholar
    • Export Citation
  • Eymard, L., A. Weill, D. Bourras, C. Guérin, P. Le Borgne, and J. M. Lefévre, 2003: Use of ship mean data for validating model and satellite flux fields during the FETCH experiment. J. Geophys. Res., 108, 8060, https://doi.org/10.1029/2001JC001207.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., E. F. Bradley, D. P. Rogers, J. B. Edson, and G. S. Young, 1996: Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere coupled-ocean atmosphere response experiment. J. Geophys. Res., 101, 37473764, https://doi.org/10.1029/95JC03205.

    • Search Google Scholar
    • Export Citation
  • Feng, Z., and Coauthors, 2021: A global high-resolution mesoscale convective system database using satellite-derived cloud tops, surface precipitation, and tracking. J. Geophys. Res. Atmos., 126, e2020JD034202, https://doi.org/10.1029/2020JD034202.

    • Search Google Scholar
    • Export Citation
  • Garg, P., S. W. Nesbitt, T. J. Lang, and G. Priftis, 2021: Diurnal cycle of tropical oceanic mesoscale cold pools. J. Climate, 34, 93059326, https://doi.org/10.1175/JCLI-D-20-0909.1.

    • Search Google Scholar
    • Export Citation
  • Gelaro, R., and Coauthors, 2017: The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). J. Climate, 30, 54195454, https://doi.org/10.1175/JCLI-D-16-0758.1.

    • Search Google Scholar
    • Export Citation
  • Giglio, D., S. T. Gille, A. C. Subramanian, and S. Nguyen, 2017: The role of wind gusts in upper ocean diurnal variability. J. Geophys. Res. Oceans, 122, 77517764, https://doi.org/10.1002/2017JC012794.

    • Search Google Scholar
    • Export Citation
  • Giglio, D., S. T. Gille, B. D. Cornuelle, A. C. Subramanian, F. J. Turk, S. Hristova-Veleva, and D. Northcott, 2022: Annual modulation of diurnal winds in the tropical oceans. Remote Sens., 14, 459, https://doi.org/10.3390/rs14030459.

    • Search Google Scholar
    • Export Citation
  • Gill, A. E., 1982: Atmosphere-Ocean Dynamics. Academic Press, 662 pp.

  • Godfrey, J. S., and A. C. M. Beljaars, 1991: On the turbulent fluxes of buoyancy, heat and moisture at the air-sea interface at low wind speeds. J. Geophys. Res., 96, 22 04322 048, https://doi.org/10.1029/91JC02015.

    • Search Google Scholar
    • Export Citation
  • Hersbach, H., and Coauthors, 2020: The ERA5 global reanalysis. Quart. J. Roy. Meteor. Soc., 146, 19992049, https://doi.org/10.1002/qj.3803.

    • Search Google Scholar
    • Export Citation
  • Hourdin, F., and Coauthors, 2020: LMDZ6A: The atmospheric component of the IPSL climate model with improved and better tuned physics. J. Adv. Model. Earth Syst., 12, e2019MS001892, https://doi.org/10.1029/2019MS001892.

    • Search Google Scholar
    • Export Citation
  • Joseph, J., M. S. Girishkumar, M. J. McPhaden, and E. P. R. Rao, 2021: Diurnal variability of atmospheric cold pool events and associated air-sea interactions in the Bay of Bengal during the summer monsoon. Climate Dyn., 56, 837853, https://doi.org/10.1007/s00382-020-05506-w.

    • Search Google Scholar
    • Export Citation
  • Klaver, R., R. Haarsma, P. L. Vidale, and W. Hazeleger, 2020: Effective resolution in high resolution global atmospheric models for climate studies. Atmos. Sci. Lett., 21, e952, https://doi.org/10.1002/asl.952.

    • Search Google Scholar
    • Export Citation
  • Ledvina, D. V., G. S. Young, R. A. Miller, and C. W. Fairall, 1993: The effect of averaging on bulk estimates of heat and momentum fluxes for the tropical western Pacific Ocean. J. Geophys. Res., 98, 20 21120 217, https://doi.org/10.1029/93JC01856.

    • Search Google Scholar
    • Export Citation
  • Lee, T., and W. T. Liu, 2005: Effects of high-frequency wind sampling on simulated mixed layer depth and upper ocean temperature. J. Geophys. Res., 110, C05002, https://doi.org/10.1029/2004JC002746.

    • Search Google Scholar
    • Export Citation
  • Lin, X., X. Zhai, Z. Wang, and D. R. Munday, 2018: Mean, variability and trend of Southern Ocean wind stress: Role of wind fluctuations. J. Climate, 31, 35573573, https://doi.org/10.1175/JCLI-D-17-0481.1.

    • Search Google Scholar
    • Export Citation
  • Lin, X., X. Zhai, Z. Wang, and D. R. Munday, 2020: Southern Ocean wind stress in CMIP5 models: Role of wind fluctuations. J. Climate, 33, 12091226, https://doi.org/10.1175/JCLI-D-19-0466.1.

    • Search Google Scholar
    • Export Citation
  • Masich, J., W. S. Kessler, M. F. Cronin, and K. R. Grissom, 2021: Diurnal cycles of near-surface currents across the tropical Pacific. J. Geophys. Res. Oceans, 126, e2020JC016982, https://doi.org/10.1029/2020JC016982.

    • Search Google Scholar
    • Export Citation
  • McPhaden, M. J., and Coauthors, 2010: The global tropical moored buoy array. Proc. OceanObs’09: Sustained Ocean Observations and Information for Society, Venice, Italy, ESA Publication, WPP-306, https://doi.org/10.5270/OceanObs09.cwp.61.

  • Ouni, M., A. Mehra, and A. Harzallah, 2021: The effect of the temporal resolution of the wind forcing on a central Mediterranean Sea model. Dyn. Atmos. Oceans, 96, 101262, https://doi.org/10.1016/j.dynatmoce.2021.101262.

    • Search Google Scholar
    • Export Citation
  • Ponte, R. M., and R. D. Rosen, 2004: Nonlinear effects of variable winds on ocean stress climatologies. J. Climate, 17, 12831293, https://doi.org/10.1175/1520-0442(2004)017<1283:NEOVWO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Praveen Kumar, B., J. Vialard, M. Lengaigne, V. S. N. Murty, and M. J. McPhaden, 2012: Tropflux: Air-sea fluxes for the global tropical oceans-description and evaluation. Climate Dyn., 38, 15211543, https://doi.org/10.1007/s00382-011-1115-0.

    • Search Google Scholar
    • Export Citation
  • Praveen Kumar, B., J. Vialard, M. Lengaigne, V. Murty, M. J. McPhaden, M. F. Cronin, F. Pinsard, and K. Gopala Reddy, 2013: Tropflux wind stresses over the tropical oceans: Evaluation and comparison with other products. Climate Dyn., 40, 20492071, https://doi.org/10.1007/s00382-012-1455-4.

    • Search Google Scholar
    • Export Citation
  • Redelsperger, J.-L., F. Guichard, and S. Mondon, 2000: A parameterization of mesoscale enhancement of surface fluxes for large-scale models. J. Climate, 13, 402421, https://doi.org/10.1175/1520-0442(2000)013<0402:APOMEO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Risien, C. M., and D. B. Chelton, 2008: A global climatology of surface wind and wind stress fields from eight years of QuikSCAT scatterometer data. J. Phys. Oceanogr., 38, 23792413, https://doi.org/10.1175/2008JPO3881.1.

    • Search Google Scholar
    • Export Citation
  • Serra, Y. L., M. F. Cronin, and G. N. Kiladis, 2007: Sub-seasonal variance of surface meteorological parameters in buoy observations and reanalyses. Geophys. Res. Lett., 34, L12708, https://doi.org/10.1029/2007GL029506.

    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., 2004: Evaluating mesoscale NWP models using kinetic energy spectra. Mon. Wea. Rev., 132, 30193032, https://doi.org/10.1175/MWR2830.1.

    • Search Google Scholar
    • Export Citation
  • Takahashi, K., 2012: Thermotidal and land-heating forcing of the diurnal cycle of oceanic surface winds in the eastern tropical Pacific. Geophys. Res. Lett., 39, L04805, https://doi.org/10.1029/2011GL050692.

    • Search Google Scholar
    • Export Citation
  • Thushara, V., and P. N. Vinayachandran, 2014: Impact of diurnal forcing on intraseasonal sea surface temperature oscillations in the Bay of Bengal. J. Geophys. Res. Oceans, 119, 82218241, https://doi.org/10.1002/2013JC009746.

    • Search Google Scholar
    • Export Citation
  • Turk, F. J., S. Hristovaveleva, and D. Giglio, 2021: Examination of the daily cycle wind vector modes of variability from the constellation of microwave scatterometers and radiometers. Remote Sens., 13, 141, https://doi.org/10.3390/rs13010141.

    • Search Google Scholar
    • Export Citation
  • Ueyama, R., and C. Deser, 2008: A climatology of diurnal and semidiurnal surface wind variations over the tropical Pacific Ocean based on the tropical atmosphere ocean moored buoy array. J. Climate, 21, 593607, https://doi.org/10.1175/JCLI1666.1.

    • Search Google Scholar
    • Export Citation
  • Wu, Y., X. Zhai, and Z. Wang, 2016: Impact of synoptic atmospheric forcing on the mean ocean circulation. J. Climate, 29, 57095724, https://doi.org/10.1175/JCLI-D-15-0819.1.

    • Search Google Scholar
    • Export Citation
  • Wu, Y., Z. Wang, C. Liu, and X. Lin, 2020: Impacts of high-frequency atmospheric forcing on Southern Ocean circulation and Antarctic Sea ice. Adv. Atmos. Sci., 37, 515531, https://doi.org/10.1007/s00376-020-9203-x.

    • Search Google Scholar
    • Export Citation
  • Wunsch, C., and R. Ferrari, 2004: Vertical mixing, energy, and the general circulation of the oceans. Annu. Rev. Fluid Mech., 36, 281314, https://doi.org/10.1146/annurev.fluid.36.050802.122121.

    • Search Google Scholar
    • Export Citation
  • Yan, Y., L. Zhang, X. Song, G. Wang, and C. Chen, 2021: Diurnal variation in surface latent heat flux and the effect of diurnal variability on the climatological latent heat flux over the tropical oceans. J. Phys. Oceanogr., 51, 34013415, https://doi.org/10.1175/JPO-D-21-0128.1.

    • Search Google Scholar
    • Export Citation
  • Yan, Y., X. Song, G. Wang, and C. Chen, 2022: Importance of high-frequency (≤30-day) wind variability to the annual climatology of the surface latent heat flux inferred from the global tropical moored buoy array. J. Geophys. Res. Oceans, 127, e2021JC018094, https://doi.org/10.1029/2021JC018094.

    • Search Google Scholar
    • Export Citation
  • Yu, Y., and Coauthors, 2020: Importance of diurnal forcing on the summer salinity variability in the East China Sea. J. Phys. Oceanogr., 50, 633653, https://doi.org/10.1175/JPO-D-19-0200.1.

    • Search Google Scholar
    • Export Citation
  • Yuan, J., and R. A. Houze Jr., 2010: Global variability of mesoscale convective system anvil structure from A-Train satellite data. J. Climate, 23, 58645888, https://doi.org/10.1175/2010JCLI3671.1.

    • Search Google Scholar
    • Export Citation
  • Zeng, X., Q. Zhang, D. Johnson, and W.-K. Tao, 2002: Parameterization of wind gustiness for the computation of ocean surface fluxes at different spatial scales. Mon. Wea. Rev., 130, 21252133, https://doi.org/10.1175/1520-0493(2002)130<2125:POWGFT>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., 2013: On the wind mechanical forcing of the ocean general circulation. J. Geophys. Res. Oceans, 118, 65616577, https://doi.org/10.1002/2013JC009086.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., and C. Wunsch, 2013: On the variability of wind power input to the oceans with a focus on the subpolar North Atlantic. J. Climate, 26, 38923903, https://doi.org/10.1175/JCLI-D-12-00472.1.

    • Search Google Scholar
    • Export Citation
  • Zhai, X., H. L. Johnson, D. P. Marshall, and C. Wunsch, 2012: On the wind power input to the ocean general circulation. J. Phys. Oceanogr., 42, 13571365, https://doi.org/10.1175/JPO-D-12-09.1.

    • Search Google Scholar
    • Export Citation
  • Zhou, S., X. Zhai, and I. A. Renfrew, 2018: The impact of high-frequency weather systems on SST and surface mixed layer in the central Arabian Sea. J. Geophys. Res. Oceans, 123, 10911104, https://doi.org/10.1002/2017JC013609.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    (a) Mean surface wind stress magnitude estimated from GTMBA observations. (b),(c) Contribution from subdaily winds to the mean wind stress magnitude and the corresponding percentage. Gray arrows and contours denote the satellite-derived mean surface wind stress in (a) and precipitation rate (mm day−1) in (b) and (c), respectively. The satellite-derived mean wind stress was calculated from the Scatterometer Climatology of Ocean Winds (SCOW; Risien and Chelton 2008).

  • Fig. 2.

    All the terms in Eq. (3) estimated from the GTMBA observations: The total contribution (black curve, representing the values at all buoy sites as shown in Fig. 1b but ordered from west to the east and from south to north; the order is listed in Table S1; black dots on the x axis denote the buoys, and the corresponding labels are their longitudes), mean wind speed change term (red curve), mean drag coefficient change term (cyan curve), wind speed variance term (blue curve), and covariance term (green curve).

  • Fig. 3.

    (a) Change in the mean wind speed due to subdaily winds ( |V|¯|V¯|). (b) Variance in subdaily wind speed variability ( |V|2). (c) Kinetic energy of subdaily winds [ ρ(u2+υ2)/2]. Gray contours in (a)–(c) denote the satellite-derived mean surface precipitation rate (mm day−1).

  • Fig. 4.

    Linear relationship between the mean wind stress magnitude from subdaily winds and the kinetic energy of subdaily winds.

  • Fig. 5.

    (a) Contribution from subdaily winds to the mean wind stress magnitude at all buoy sites. (b) Change in the mean wind speed due to subdaily winds (dashed curve) and the associated contribution (solid curve). (c) Variance in subdaily wind speed variability (dashed curve) and the associated contribution (solid curve). (d) Kinetic energy of subdaily winds. Black, red, and blue curves denote values estimated from the GTMBA observations, ERA5 product, and MERRA2 product, respectively. The x axes in (a)–(d) are the same as those in Fig. 2.

  • Fig. 6.

    Kinetic energy spectrum of surface winds (variance-preserving plots, solid curve) and the ratio between the reanalysis estimates and observations (percentage, dashed curve). Black, red, and blue curves denote values estimated from the GTMBA observations, ERA5 product, and MERRA2 product, respectively.

  • Fig. 7.

    Percentage error (red curve) of the mean wind stress magnitude in the (a) ERA5 and (b) MERRA2 products based on the GTMBA observations. The black curve denotes the percentage relative error associated with inaccurate subdaily winds. The x axes are the same as those in Fig. 2.

  • Fig. 8.

    Regression model between the square of the gustiness value (G2) and precipitation rate (P) for the (left) time-independent and (right) time-dependent schemes.

  • Fig. 9.

    Effect of subdaily winds on wind stress vector (red arrows). Black arrows denote the long-term mean wind stress vector. Red arrows are amplified by 10.

  • Fig. 10.

    (a) Contribution from subdaily winds to the mean wind stress magnitude estimated from the GTMBA observation using the COARE algorithm (black curve) and from the ERA5 reanalysis using the COARE (red curve) and ECMWF (blue curve) algorithms. (b) Mean absolute percentage errors of the subdaily wind contribution (red curve), convective gustiness (black curve) and drag coefficient (blue curve) calculated using the ECMWF algorithm with respect to those calculated using the COARE algorithm.