1. Introduction


Note that the wind power input to the large-scale geostrophic circulation, which is the focus of the present study, is only part of the total wind power input to the ocean. However, the wind power input to surface waves and surface ageostrophic currents, although of much larger magnitude, tends to be dissipated within the surface layer and is therefore not available to the deep ocean (e.g., von Storch et al. 2007; Zhai et al. 2009).
Over most of the ocean the speed of ocean surface currents is at least one order of magnitude smaller than that of the 10-m wind, and the wind stress is thus often computed using the 10-m wind alone, neglecting the contribution from the surface ocean currents. Recently, a few studies (e.g., Duhaut and Straub 2006; Zhai and Greatbatch 2007; Hughes and Wilson 2008; Scott and Xu 2009) have found a positive bias in calculations of wind power input (about 20%–30%) if the relative air–sea velocities are not accounted for in the stress calculation. Not accounting for ocean surface velocity in the stress calculation is hereafter referred to as the “resting ocean approximation.”
Attempts to interpret physically the cause of the positive bias have, in the past, emphasized the smaller spatial scales of ocean currents and the vortex structure of ocean eddies (Zhai and Greatbatch 2007; Hughes and Wilson 2008), schematically illustrated in Fig. 1a. When the wind blows over the eddy shown, the stress is smaller on the northern side because wind and current are aligned and larger on the southern side where they oppose each other. The net effect, when integrated spatially, takes energy out of the eddy: that is, the wind mechanically damps the eddy.

Schematic illustrating the mechanical damping effect by the wind. The dashed line denotes the wind; the solid line denotes the surface ocean current; and the open arrow denotes the surface stress, which depends on the relative motion between the air and the surface ocean. (a) The wind blows over an eddy, where the stress is smaller on the northern side because wind and current are aligned and larger on the southern side because they oppose each other. The net effect, when integrated spatially, takes energy out of the eddy: that is, the wind mechanically damps the eddy. (b) An oscillating wind blows over a steady current. When the wind is aligned with the current, the stress is smaller than in the motionless ocean case and hence the wind does less positive work. When the wind opposes the current, the stress is larger and hence does more negative work. Integrated over time, the energy is removed from the current, damping the flow. Both damping effects illustrated in the schematic are termed the wind mechanical damping effect.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

Schematic illustrating the mechanical damping effect by the wind. The dashed line denotes the wind; the solid line denotes the surface ocean current; and the open arrow denotes the surface stress, which depends on the relative motion between the air and the surface ocean. (a) The wind blows over an eddy, where the stress is smaller on the northern side because wind and current are aligned and larger on the southern side because they oppose each other. The net effect, when integrated spatially, takes energy out of the eddy: that is, the wind mechanically damps the eddy. (b) An oscillating wind blows over a steady current. When the wind is aligned with the current, the stress is smaller than in the motionless ocean case and hence the wind does less positive work. When the wind opposes the current, the stress is larger and hence does more negative work. Integrated over time, the energy is removed from the current, damping the flow. Both damping effects illustrated in the schematic are termed the wind mechanical damping effect.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
Schematic illustrating the mechanical damping effect by the wind. The dashed line denotes the wind; the solid line denotes the surface ocean current; and the open arrow denotes the surface stress, which depends on the relative motion between the air and the surface ocean. (a) The wind blows over an eddy, where the stress is smaller on the northern side because wind and current are aligned and larger on the southern side because they oppose each other. The net effect, when integrated spatially, takes energy out of the eddy: that is, the wind mechanically damps the eddy. (b) An oscillating wind blows over a steady current. When the wind is aligned with the current, the stress is smaller than in the motionless ocean case and hence the wind does less positive work. When the wind opposes the current, the stress is larger and hence does more negative work. Integrated over time, the energy is removed from the current, damping the flow. Both damping effects illustrated in the schematic are termed the wind mechanical damping effect.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
An overlooked aspect, however, is that the atmospheric wind tends to vary on faster time scales than do surface geostrophic currents. At a fixed location in space, the synoptically varying wind tends to damp the underlying ocean currents, regardless of their spatial structure. This situation is illustrated schematically in Fig. 1b for the case of an oscillating wind blowing over a steady current. When the wind is aligned with the current, the stress is smaller than in the motionless ocean case and hence the wind does less positive work. When the wind opposes the current, the stress is larger and hence does more negative work. Integrated over time, the energy is removed from the current, damping the flow. Both damping effects illustrated in Fig. 1 are termed here the “wind mechanical damping effect.”
In this paper, the wind power input to the ocean general circulation is reexamined using available observations. Of particular interest is the role of the synoptic wind in supplying energy to the ocean general circulation and in taking energy out of the ocean when the ocean surface velocity is taken into account in the wind stress calculation.
2. Theory
a. With the resting ocean approximation












The time-averaged wind stress
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The time-averaged wind stress
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
The time-averaged wind stress
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
b. Without the resting ocean approximation










In the next section, the above ideas are tested using available observations.
3. Data
Following Hughes and Wilson (2008), the absolute sea surface height (SSH) for the period from January 1995 to December 2008 is obtained by combining the ocean mean dynamic height from the Maximenko and Niiler (2005) product and the SSH anomaly product compiled by the Collecte Localisation Satellites (CLS) Space Oceanographic Division of Toulouse, France. The SSH anomaly values result from merging the Ocean Topography Experiment (TOPEX)/Poseidon and European Remote Sensing Satellite-1 (ERS-1)/ERS-2 along-track SSH measurements for a temporal gridding of 7 days on a ⅓° Mercator grid (Le Traon et al. 1998). The Maximenko and Niiler (2005) product, which integrates information from surface drifters, satellite altimetry, surface winds, and the Gravity Recovery and Climate Experiment (GRACE) gravity mission for the period from 1992 to 2002, is interpolated from a 1/2° latitude–longitude grid to the same grid as the SSH anomalies. Surface currents ug are then computed through geostrophy from the absolute SSH with temporal resolution of 7 days. Readers are referred to Hughes and Wilson (2008) and Scott and Xu (2009) for a detailed discussion of errors associated with each product.
The 6-hourly, daily, and monthly 10-m wind fields are taken from the National Centers for Environmental Prediction (NCEP) reanalysis product (Kalnay et al. 1996) and interpolated to the same grid as ug. Surface wind stress is then computed from 10- m wind using the Large et al. (1994) formula for the drag coefficient and ρa = 1.223 kg m−3. There are subtle issues associated with the drag coefficient (e.g., whether the same drag coefficient should be used for both the resting and nonresting ocean cases), but these issues are not considered in the present study. The wind power input to the ocean general circulation is finally computed as
4. Results
a. With the resting ocean approximation
Figures 3a,b show the average rate of power input to the surface geostrophic currents by the NCEP monthly and 6-hourly winds, respectively. The spatial pattern in both cases is very similar to that found in previous studies (Wunsch 1998; Hughes and Wilson 2008; Scott and Xu 2009; Roquet et al. 2011), with the majority of the wind power input entering in the Southern Ocean. However, regions of both positive and negative power input become more pronounced when the 6-hourly wind is used. This change in magnitude can be clearly seen in Fig. 3c, which shows the power input by the 6-hourly wind minus that by the monthly wind. Positive power input is strongly enhanced in the Southern Ocean, Gulf Stream, and Kuroshio Extension when the synoptic values are used. In Figs. 3a,c, most of the wind power input north of 30°N in the North Atlantic and North Pacific is seen owing to wind periods between 6 h and a month, instead of the monthly- or climatological-mean values. The significant increase of wind power input to the surface geostrophic currents in the storm-track regions is consistent with the argument presented in section 2a because in these regions the synoptic wind variability is particularly strong. Figures 4a,b show the atmospheric kinetic energy

Power input (W m−2) by (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind with the resting ocean approximation. (c) Shown is (b) minus (a). The color bar is saturated. The maximum value in the Southern Ocean in (b) is about 0.04 W m−2.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

Power input (W m−2) by (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind with the resting ocean approximation. (c) Shown is (b) minus (a). The color bar is saturated. The maximum value in the Southern Ocean in (b) is about 0.04 W m−2.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
Power input (W m−2) by (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind with the resting ocean approximation. (c) Shown is (b) minus (a). The color bar is saturated. The maximum value in the Southern Ocean in (b) is about 0.04 W m−2.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The mean atmospheric kinetic energy (
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The mean atmospheric kinetic energy (
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
The mean atmospheric kinetic energy (
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
Integrated globally, the power input to the surface geostrophic currents by the monthly-mean wind is about 0.42 TW (1 TW = 1012 W), whereas that by the 6-hourly wind is about 0.72 TW (see Table 1), an increase of over 70%. The explanation for this significant increase in wind power input, as outlined in section 2, lies in the quadratic dependence of wind stress on wind itself, such that the high-frequency wind contributes to the time-averaged wind stress. Wind power input to the ocean general circulation depends on the wind energy integrated over the whole spectrum. Power input by the daily wind integrates to about 0.65 TW, slightly less than that by the 6-hourly wind, and the difference is again concentrated in the storm-track regions (not shown). Readers are referred to Scott and Xu (2009) for an in-depth discussion of error estimates, where the authors found the uncertainty is about 10% of the mean using a range of wind and ocean current products.
The wind power input to the ocean general circulation by the monthly-mean NCEP wind stress taken directly from the reanalysis product and stresses computed from the 6-hourly, daily-mean, and monthly-mean NCEP 10-m winds. All numbers are globally integrated values in TW. In theory, estimates without the resting ocean approximation are more accurate. The most trustable estimates are highlighted in bold.


It is possible that the time-dependent wind stress associated with the 6-hourly wind projects more effectively onto the time-dependent ocean surface velocities than does the stress associated with the monthly wind, and it could lead to an increase in wind power input. This possibility was tested by calculating the projection of time-dependent winds onto time-dependent currents, but no significant effect was found (see appendix B for a detailed discussion).
Figure 5 shows the time-mean wind stress by averaging the monthly NCEP wind stress taken directly from the reanalysis product, the time-mean wind stress computed from the 6-hourly wind using (1), and the difference between them. Note again that the monthly NCEP wind stress is a monthly average of the instantaneous surface wind stress at every 20-min NCEP model time step, and it therefore includes contributions from wind variability at all periods longer than 40 min. As a result, the wind stress computed using the 6-hourly wind still underestimates the NCEP modeled wind stress (Fig. 5c). This difference explains why the present estimate of wind power input to the ocean general circulation using the 6-hourly wind is somewhat less than previous estimates of ~0.9 TW when the resting ocean approximation is used (Wunsch 1998; Hughes and Wilson 2008; Scott and Xu 2009). Indeed, when the wind power input is computed directly using the monthly-mean NCEP wind stress, a value of 0.87 TW is obtained, closer to previous estimates. We will discuss this issue further in section 5.

The time-mean wind stress (N m−2) by averaging (a) the monthly-mean NCEP wind stress taken directly from the reanalysis product and (b) the stress computed from the 6-hourly NCEP wind using (1) with uo = 0. (c) Shown is (a) minus (b).
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The time-mean wind stress (N m−2) by averaging (a) the monthly-mean NCEP wind stress taken directly from the reanalysis product and (b) the stress computed from the 6-hourly NCEP wind using (1) with uo = 0. (c) Shown is (a) minus (b).
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
The time-mean wind stress (N m−2) by averaging (a) the monthly-mean NCEP wind stress taken directly from the reanalysis product and (b) the stress computed from the 6-hourly NCEP wind using (1) with uo = 0. (c) Shown is (a) minus (b).
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
b. Without the resting ocean approximation
In section 2b, it is argued that the synoptically varying wind can significantly reduce the energy of the surface geostrophic currents through the direct wind damping effect when the resting ocean approximation is not used. Figure 6 shows the reduction in power input by the monthly and 6-hourly NCEP winds because of the mechanical damping effect illustrated in Fig. 1. The spatial pattern is, again, very similar to that found in previous studies (Hughes and Wilson 2008; Scott and Xu 2009). When integrated globally, the reduction in power input by the 6-hourly wind is about 0.25 TW, almost double the reduction by the monthly wind (0.14 TW; see Table 1). Hughes and Wilson (2008) estimated the reduction in power input to be ~0.19 TW using a weekly-mean wind, which fits well with the values in Table 1. One may anticipate that the mechanical damping effect by the wind will be even greater if the wind output from every NCEP time step is used.

The reduction in power input (W m−2) associated with (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind, after removing the resting ocean approximation. The color bar is as in Fig. 3 for easy comparison.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The reduction in power input (W m−2) associated with (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind, after removing the resting ocean approximation. The color bar is as in Fig. 3 for easy comparison.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
The reduction in power input (W m−2) associated with (a) the monthly NCEP wind and (b) the 6-hourly NCEP wind, after removing the resting ocean approximation. The color bar is as in Fig. 3 for easy comparison.
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
5. Discussion
The quadratic dependence of the stress law on the wind speed produces a qualitative change in the calculation of power input to the ocean when synoptic weather systems are present. This result also amplifies the process by which the wind field damps eddy motions when the ocean surface velocity is accounted for. By comparing calculations from the NCEP 6-hourly, daily-mean, and monthly-mean estimated winds, the following is found
Power input to the ocean general circulation is increased by roughly 70% when 6-hourly winds are used instead of monthly winds.
With the resting ocean approximation, the power increase is from 0.42 to 0.72 TW (an increase of 71%). Much of the increase occurs in the storm-track regions of the Southern Ocean, Gulf Stream, and Kuroshio Extension.
Without the resting ocean approximation, the power increase is from 0.28 to 0.47 TW (an increase of 68%).
With the resting ocean approximation, the power input by the 6-hourly wind is found to be about 0.2 TW less than that by the monthly-mean NCEP wind stress, which can be explained by the contribution of wind with periods less than 6 h to the monthly NCEP wind stress. However, this 0.2-TW difference raises a serious question: is the wind variability at such high frequencies (less than a few hours) really important for inputting energy to the geostrophic ocean circulation, which eventually finds its way to feed the deep-ocean mixing? Or does it simply generate shear, mixing, and waves (e.g., near-inertial waves) in the surface layer, with the majority of its energy input being dissipated there?1 If the latter is true, the wind power input to the ocean general circulation that is relevant to deep-ocean mixing may be less than previously thought (~0.5 TW instead of ~1 TW).
Pathways by which the energy input by wind working on the sea surface enters the interior ocean circulation remain the subject of considerable uncertainty, partly because of the complex turbulent structure of the near-surface boundary layer (e.g., von Storch et al. 2007). As Roquet et al. (2011) show, the assumption of Ekman layer physics implies that the energy is pumped into the interior, sometimes far from the region of surface working. The accuracy of the wind products remains obscure, and previous estimates of possible factor of two errors in the total power input as calculated by these methods are probably still appropriate.
Finally, forcing ocean-only models is problematic. If ocean models are forced directly by the NCEP wind stress, they will be forced too strongly since the ocean surface velocity is not accounted for in the wind stress calculation. On the other hand, if ocean models are forced by the 6-hourly wind using (1), the energy input by the higher-frequency wind will be missed. Neither of these two caveats are trivial. Our study points toward the importance of coupled atmosphere–ocean models.
Acknowledgments
XZ thanks Benjamin Grandey and Chris Wilson for helpful discussions. Financial support was provided by the U.K. Natural Environment Research Council. HLJ is funded by a Royal Society University Research Fellowship, for which she is grateful. CW was supported by the George Eastman Visiting Professorship of Balliol College, Oxford. DPM is partially supported by the Oxford Martin School. We thank two anonymous reviewers for their constructive and insightful comments that led to significant improvements.
APPENDIX A
Mean Wind Stress Formula with Stochastic Variability
a. With the resting ocean approximation















b. Without the resting ocean approximation












APPENDIX B
The Time-Dependent Wind Power Input
It is possible that the time-dependent wind stress associated with the 6-hourly wind projects more effectively onto the time-dependent ocean surface velocities than does the stress associated with the monthly wind, and it could lead to an increase in wind power input. This possibility is tested here by considering the time-dependent component of the power input: that is,

The time-dependent component of the power input
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1

The time-dependent component of the power input
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
The time-dependent component of the power input
Citation: Journal of Physical Oceanography 42, 8; 10.1175/JPO-D-12-09.1
REFERENCES
Dewar, W. K., and P. D. Killworth, 1995: Do fast gravity waves interact with geostrophic motions? Deep-Sea Res., 42, 1063–1081.
Duhaut, T. H., and D. N. Straub, 2006: Wind stress dependence on ocean surface velocity: Implications for mechanical energy input to ocean circulation. J. Phys. Oceanogr., 36, 202–211.
Ferrari, R., and C. Wunsch, 2009: Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annu. Rev. Fluid Mech., 41, 253–282.
Hughes, C. W., and C. Wilson, 2008: Wind work on the geostrophic ocean circulation: An observational study on the effect of small scales in the wind stress. J. Geophys. Res., 113, doi:10.1029/2007JC004371.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471.
Large, W. G., J. C. McWilliams, and S. C. Doney, 1994: Oceanic vertical mixing: A review and a model with a nonlocal boundary layer parameterization. Rev. Geophys., 32, 363–403.
Le Traon, P. Y., F. Nadal, and N. Ducet, 1998: An improved mapping method of multisatellite altimeter record. J. Atmos. Oceanic Technol., 15, 522–534.
Maximenko, N. A., and P. P. Niiler, 2005: Hybrid decade-mean global sea level with mesoscale resolution. Recent Advances in Marine Science and Technology, N. Saxena, Ed., PACON International, 55–59.
Roquet, F., C. Wunsch, and G. Madec, 2011: On the patterns of wind-power input to the ocean circulation. J. Phys. Oceanogr., 41, 2328–2342.
Scott, R. B., and Y. Xu, 2009: An update on the wind power input to the surface geostrophic flow of the World Ocean. Deep-Sea Res., 56, 295–304.
Thompson, K. R., R. F. Marsden, and D. G. Wright, 1983: Estimation of low-frequency wind stress fluctuations over the open ocean. J. Phys. Oceanogr., 13, 1003–1011.
von Storch, J. S., H. Sasaki, and J. Marotzke, 2007: Wind-generated power input to the deep ocean: An estimate using a 1/10 general circulation model. J. Phys. Oceanogr., 37, 657–672.
Wunsch, C., 1998: The work done by the wind on the oceanic general circulation. J. Phys. Oceanogr., 28, 2332–2340.
Zhai, X., and R. J. Greatbatch, 2007: Wind work in a model of the northwest Atlantic Ocean. Geophys. Res. Lett., 34, L04606, doi:10.1029/2006GL028907.
Zhai, X., R. J. Greatbatch, C. Eden, and T. Hibiya, 2009: On the loss of wind-induced near-inertial energy to turbulent mixing in the upper ocean. J. Phys. Oceanogr., 39, 3040–3045.
To leading order, the geostrophic flow does not exchange energy with the internal gravity wave field (Dewar and Killworth 1995).