A Simple Model for the Skewness of Global Sea Surface Winds

Adam Hugh Monahan School of Earth and Ocean Sciences, University of Victoria, Victoria, British Columbia, and Earth Systems Evolution Program, Canadian Institute for Advanced Research, Toronto, Ontario, Canada

Search for other papers by Adam Hugh Monahan in
Current site
Google Scholar
PubMed
Close
Full access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

A strong linear relationship between the mean and skewness of global sea surface winds (both zonal and meridional) is shown to exist, such that where the wind component is on average positive, it is negatively skewed (and vice versa). This relationship is observed in reanalysis, satellite, and buoy data. This relationship between the mean and skewness fields of sea surface winds follows from the nonlinear surface drag predicted for a turbulent boundary layer by Monin–Obukhov similarity theory since forcing perturbations speeding the wind up are subject to a stronger drag force than perturbations slowing it down. Furthermore, it is demonstrated that the results of an empirical fit of observed surface winds to a stochastic differential equation presented in a recent study by Sura are consistent with the white-noise limit of the momentum equations for a turbulent boundary layer subject to fluctuating forcing, albeit with a somewhat different physical interpretation.

Corresponding author address: Adam Hugh Monahan, School of Earth and Ocean Sciences, University of Victoria, P.O. Box 3055, STN CSC, Victoria, BC V8P 5C2, Canada. Email: monahana@uvic.ca

Abstract

A strong linear relationship between the mean and skewness of global sea surface winds (both zonal and meridional) is shown to exist, such that where the wind component is on average positive, it is negatively skewed (and vice versa). This relationship is observed in reanalysis, satellite, and buoy data. This relationship between the mean and skewness fields of sea surface winds follows from the nonlinear surface drag predicted for a turbulent boundary layer by Monin–Obukhov similarity theory since forcing perturbations speeding the wind up are subject to a stronger drag force than perturbations slowing it down. Furthermore, it is demonstrated that the results of an empirical fit of observed surface winds to a stochastic differential equation presented in a recent study by Sura are consistent with the white-noise limit of the momentum equations for a turbulent boundary layer subject to fluctuating forcing, albeit with a somewhat different physical interpretation.

Corresponding author address: Adam Hugh Monahan, School of Earth and Ocean Sciences, University of Victoria, P.O. Box 3055, STN CSC, Victoria, BC V8P 5C2, Canada. Email: monahana@uvic.ca

1. Introduction

Interactions between the ocean and the atmosphere are mediated through their respective boundary layers, such that the character of the turbulence in these boundary layers determines the rate of exchange of momentum, of energy, and of material substances such as water vapor, aerosols, and gases. Classical similarity theory predicts that the character of the turbulence in the atmospheric boundary layer can be simply related to the density stratification and the vertical shear of the wind, each averaged over a typical eddy time scale of 1 h or so (Stull 1997). This relationship can then be used to develop bulk parameterizations of turbulent air–sea fluxes as functions of the eddy-averaged winds at a given altitude (typically 10 m). Denoting the eddy-averaging of the variable x by an overbar with an e, xe, standard bulk parameterizations for the turbulent flux of a quantity Λ across the ocean–atmosphere interface take the form
i1520-0469-61-16-2037-e1
where |U| is the wind speed and F is a (typically nonlinear) function which will generally depend on the stability of the boundary layer and on the sea state (Taylor 2000); this is discussed in more detail in section 4.
It is quite often the case, however, that only observations of boundary layer variables averaged over time scales longer than the eddy time scale are available for estimating air–sea fluxes, in which case a time-averaged version of Eq. (1) must be used. Denoting averaging of the variable x over the time scale T (assumed to be longer than the eddy-averaging time scale) by 〈xT, the T-averaged flux will not equal the flux corresponding to the T-averaged variables if the function F is nonlinear (which it typically is): in general,
i1520-0469-61-16-2037-e2
The importance of time-averaging on estimates of air– sea momentum fluxes has been demonstrated by (for example) Wright and Thompson (1983), Mahrt et al. (1996), and Ponte and Rosen (2004), and for air–sea CO2 fluxes by Bates and Merlivat (2001).
The T-averaged eddy fluxes can, however, be represented in terms of the T-averaged boundary layer quantities plus correction terms involving higher-order moments of fluctuations on time scales shorter than T. Using a Taylor series expansion, we obtain
i1520-0469-61-16-2037-e3
where the remaining terms represent higher-order moments of |U|e, moments of Λe, and cross moments of |U|e and Λe. If F is a polynomial in |U|e and Λe, then an expansion of (3) will be exact with a finite number of terms; otherwise, it will be approximate. If the higher-order moments in the expansion (3) can be related to the T-averaged quantities, then workable closures for the eddy fluxes in terms of the T-averaged variables can be derived. This paper demonstrates the existence of a strong linear relationship between the mean and the skewness (normalized third-order moment) of wind components (both zonal and meridional) in the marine boundary layer, such that the more positive the mean wind component is at a given location, the more negative is the skewness (and vice versa). Furthermore, it will be shown that this relationship follows naturally from the nonlinear surface drag predicted by Monin–Obukhov similarity theory for a turbulent boundary layer.

Estimation of eddy fluxes is not the only motivation for studying the statistics of sea surface winds. Recent years have seen considerable interest in the effect on oceanic variability of fluctuations in atmospheric forcing, for example, in the context of the dynamics of El Niño–Southern Oscillation (e.g., Penland 1996; Kleeman and Moore 1997; Thompson and Battisti 2000), of midlatitude gyre variability (e.g., Moore 1999; Sura et al. 2001), of the stability of the meridional overturning circulation (e.g., Aeberhardt et al. 2000; Monahan 2002a,b), and of variability in open-ocean deep convection (e.g., Kuhlbrodt et al. 2001; Kuhlbrodt and Monahan 2003). In a recent study, Sura (2003) fit observed surface winds from the Southern and Pacific Oceans to a scalar stochastic differential equation, and found evidence that the dynamics is characterized by multiplicative noise [noise with intensity dependent on the state variable; see Penland (2003a,b) for a discussion]. This result was interpreted as reflecting increasing gustiness at stronger wind speeds. We will demonstrate that the results of Sura (2003) can largely be accounted for by boundary layer momentum dynamics subject to fluctuating forcing, in an appropriate limit, and that the multiplicative noise can be understood as arising from turbulent fluctuations in the surface drag coefficient.

A description of the data used in this study is presented in section 2. The relationship between the mean and the skewness fields of both the zonal and meridional winds is characterized in section 3. Section 4 describes the boundary layer momentum equations subject to fluctuating forcing, and provides a minimal theory for the observed relationship between the mean and skewness of wind components. In section 5, the white-noise limit of the model presented in section 4 is compared to the results of Sura (2003). Finally, a discussion and conclusions are presented in section 6.

2. Data

Five datasets were analyzed in this study:

  1. Six-hourly National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis sea level pressure, surface air temperature, and 995-hPa (the lowest sigma level of the reanalysis model) zonal and meridional winds, available on a 2.5° × 2.5° grid from 1 January 1948 to 31 December 2002 (Kalnay et al. 1996). This is the primary dataset considered in this study. [Data available online at the National Oceanic and Atmospheric Administration–Cooperative Institute for Research in Environmental Sciences (NOAA–CIRES) Climate Diagnostics Center, http://www.cdc.noaa.gov/.]

  2. Six-hourly European Centre for Medium-Range Weather Forecasts (ECMWF) 40-yr reanalysis period (ERA-40) 10-m zonal and meridional winds, available on a 2.5° × 2.5° grid from 1 September 1957 to 31 August 2002 (Simmons and Gibson 2000). (Data available online at http://data.ecmwf.int/data/d/era40/.)

  3. Six-hourly 10-m zonal and meridional wind data from a blend of National Aeronautics and Space Administration's (NASA's Quick Scatterometer (QuikSCAT) scatterometer observations with NCEP analyses, available on a 0.5° × 0.5° grid from 19 July 1999 to 31 March 2003 (Chin et al. 1998; Milliff et al. 1999). [Data available online at the NCAR Data Support Section (DSS), http://dss.ucar.edu/datasets/ ds744.4/.]

  4. Variationally gridded six-hourly Special Sensor Microwave Imager (SSM/I) 10-m zonal and meridional wind observations, available on a 1° × 1° grid from 1 July 1987 to 31 December 2001 (Atlas et al. 1996). [Data available online at the NASA Jet Propulsion Laboratory (JPL) Distributed Active Archive Center, http://podaac.jpl.nasa.gov.]

  5. Hourly 4-m zonal and meridional winds from 65 buoys in the Tropical Atmosphere–Ocean (TAO) array (McPhaden et al. 1998). (Data available online at the TAO Project Office at http://www.pmel.noaa.gov/tao/data_deliv/.) Wind data were reported from these buoys intermittently between 1990 and 2001; the number of hourly measurements at individual buoys varied from a minimum of 4000 to more than 80 000.

No preprocessing, such as filtering or removing the annual cycle, was carried out on any of these datasets.

3. Moments of surface ocean winds

The mean, standard deviation, and skewness fields of the sea surface zonal and meridional wind components from the NCEP–NCAR reanalysis are presented in Figs. 1 and 2, respectively. The skewness of a random variable x is the normalized third-order moment,
i1520-0469-61-16-2037-e4
where angle brackets 〈 · 〉 denote long-term time averaging. The mean sea surface zonal wind field displays the familiar tropical easterlies and midlatitude westerlies, while the meridional wind field shows strong flow only on the eastern flanks of the subtropical highs. The standard deviation fields of both zonal and meridional winds display variability minima in the Tropics and subtropics, and variability maxima in the storm tracks. The skewness field of the zonal wind is generally positive in the Tropics and negative in the midlatitudes, while that of the meridional wind is generally small except along the eastern flanks of the subtropical highs, where skewness is positive in the Northern Hemisphere and negative in the Southern Hemisphere.

A linear relationship between the mean and skewness fields of both the zonal and meridional sea surface wind components is evident from Figs. 1 and 2: positive (negative) mean wind components are associated with negative (positive) skewness. This spatial anticorrelation is further illustrated in Fig. 3, which presents scatterplots over all ocean grid points of the mean versus the skewness of the zonal and meridional winds for the NCEP– NCAR reanalysis, the blended QuikSCAT–NCEP analysis, and the TAO buoy datasets. Similar plots for the ERA-40 and SSM/I surface winds (not shown) display the same relationship between the mean and skewness fields. The spatial correlation coefficients between the mean and skewness fields of the surface winds for the different datasets are displayed in Table 1; these correlation coefficients display essentially no seasonal variability.

Zonal wind skewness is seen to range between approximately −1 and 1.5; skewness toward the upper and lower end of this range is strong enough to be manifest in the probability density functions. Kernel density estimates (Silverman 1986) of the probability density functions of the zonal wind from the NCEP– NCAR reanalysis at three representative grid points [(50°S, 110°W), (20°N, 160°W), and (50°N, 60°W)] are presented in Fig. 4. The negative and positive skewness characteristic of the zonal winds in the Southern Ocean and the tropical Pacific respectively are evident in this figure, as is the relatively weak skewness in the central North Atlantic. It is interesting to note the relatively weak skewness of the wind components at Ocean Station Papa (50°N, 145°E) evident in Figs. 1 and 2, which is consistent with the observation in Monahan and Denman (2004) that the zonal and meridional winds at this location are nearly Gaussian.

The negative skewness of zonal winds in the Southern Ocean was noted by Sura (2003), who attributed the skewness to the presence of multiplicative noise (that is, noise with state-dependent intensity) in an empirical stochastic equation for the zonal wind. In the next two sections, we will demonstrate that while the dynamics of the winds may involve multiplicative noise (in an appropriate limit, to be defined in section 5), the skewness of the wind components and the relationship between the skewness and the mean do not require multiplicative noise, but can be understood as a consequence of the nonlinear drag law characteristic of turbulent boundary layers.

4. A minimal theory for the skewness of sea surface wind components

Neglecting advection terms, the horizontal momentum equations averaged over a layer from the surface to the altitude h read
i1520-0469-61-16-2037-e5
where u and υ are respectively the depth-averaged zonal and meridional components of the wind, ϕ and λ are respectively the latitude and the longitude, a is the radius of the earth, ρ is the air density, p is the pressure, f is the Coriolis parameter, and τx and τy are respectively the zonal and meridional components of the stress vector. Monin–Obukhov similarity theory for the boundary layer (Haltiner and Williams 1980) implies that these equations may be expressed as
i1520-0469-61-16-2037-e7
where
i1520-0469-61-16-2037-e9
is the drag coefficient. The drag coefficient is a function of the reference altitude za (typically 10 m), of the momentum roughness length z0, and of the stability of the boundary layer through the term
i1520-0469-61-16-2037-e10
where L is the Obukhov length and Φm is a dimensionless function characterizing buoyant effects on boundary layer turbulence. The Obukhov length L is a function of the local wind speed |u| = (u2 + υ;th2)1/2 and the surface buoyancy flux. The von Karman constant k0 is a fixed number. Note that it has been assumed that the wind at the reference level za is approximately equal to the wind averaged from z = 0 to z = h.

The sea surface momentum roughness length z0 is determined by the sea state, which is affected both by remotely generated swell and by surface waves generated by local winds. Consequently, z0 is a function of the local wind speed |u|, and the drag coefficient will depend on |u| even in the neutral (L → ∞) limit in which buoyant effects on turbulence are negligible. The precise form of the dependence of the neutral drag coefficient on swell and local winds is a field of active research, and a number of different formulations have been proposed (e.g., Mahrt et al. 1996; Rieder 1997; Taylor 2000; Grachev and Fairall 2001; Taylor and Yelland 2001; Fairall et al. 2003; Mahrt et al. 2003). These formulations differ in detail, but there is a general consensus that the drag coefficient increases with |u|, except perhaps for very light winds.

The simplest parameterization for (τx, τy) at z = h assumes an eddy viscosity κ so that
i1520-0469-61-16-2037-e11
where (u*, υ*) is the (nondepth-averaged) wind vector. Equation (11) can be approximated using finite differences as
i1520-0469-61-16-2037-e12
where U and V are, respectively, the zonal and meridional winds in the atmospheric layer above z = h.

Evidently, fluctuations in u and υ can result from fluctuations in

  • the pressure gradient forces,

  • the drag coefficient cD,

  • the wind components U and V above the surface layer, and

  • the size of the viscosity κ.

Variability in U, V, and the pressure gradients will reflect variations in weather above the surface layer. The drag coefficient will fluctuate with changes in stability of the boundary layer (as measured by the Obukhov length L) and the surface wave field (as measured by the roughness length z0). Similarly, κ will vary depending on the character and intensity of the turbulence in the boundary layer.

For convenience, we define
i1520-0469-61-16-2037-e13
so that
i1520-0469-61-16-2037-e16
where
i1520-0469-61-16-2037-e18
As already noted, the statistics of u and υ will be determined by the fluctuations in Πu, Πυ, Γ, and K; we will model these fluctuations as stochastic processes. The simplest stochastic representation of (16), (17) is obtained by representing Πu, Πυ, Γ, and K as independent Ornstein–Uhlenbeck (red noise) processes:
i1520-0469-61-16-2037-e19
where
i1520-0469-61-16-2037-e23
is the neutral (L → ∞) drag coefficient (which may depend on |u| through z0) divided by h. The fluctuations are assumed to be governed by the stochastic differential equations (SDEs)
i1520-0469-61-16-2037-e24
where (τu, τυ, τγ, τμ) are relaxation time scales, (σu, συ, σγ, σμ) tune the strength of the fluctuations, and (1, 2, 3, 4) are independent white-noise processes:
it1jt2δijδt1t2
where δij is the Kronecker delta and δ(t1t2) is the Dirac delta function. A basic introduction to stochastic differential equations is given in Penland (2003a); more detailed expositions appear in Penland (2003b) and Gardiner (1997). Note that the Ornstein–Uhlenbeck process z described by the stochastic differential equation
i1520-0469-61-16-2037-e29
will (after initial transients have died out) have a mean of 0, the stationary autocovariance function
i1520-0469-61-16-2037-e30
and will converge to the white-noise process σẆ in the limit that τ → 0.
Together, Eqs. (16)–(27) define a stochastic differential equation for the six-dimensional Markov process X = (u, υ, Πu, Πυ, γ, μ):
AX
where
i1520-0469-61-16-2037-e32
In the terminology of SDEs, the “deterministic” vector A(X) is known as the drift and the matrix 𝗕 multiplying the white-noise vector is known as the diffusion. As the matrix 𝗕 is independent of the state variable X, the noise is said to be additive.
To demonstrate that the relationship between the mean and skewness of sea surface wind components described in the previous section follows from Eq. (31), the statistics of Πu and Πυ were estimated from the NCEP–NCAR reanalysis data as follows. First, the pressure gradient fields were estimated using centered finite-difference approximations to the derivatives, and surface air density ρ was calculated from sea level pressure and surface air temperature. Next, the Πu and Πυ fields were calculated using Eqs. (14) and (15) and the NCEP– NCAR 995-hPa u and υ fields. At each grid point; the sample mean and lag autocovariance functions of Πu and Πυ were estimated, yielding Πu, τu, Πυ, and τυ. The noise strength parameters σu and συ were then estimated from the standard deviation fields (std) of Πu and Πυ using Eq. (30) with t1 = t2:
i1520-0469-61-16-2037-e35
Note that it follows from Eqs. (35) and (36) that the units of σu and συ are m s−3/2. As it was found that the autocorrelation times τu and τυ were generally on the order of 1 day, a constant value of
τuτυ
was used in these and all subsequent calculations.
Because finite differences of the NCEP–NCAR reanalysis surface pressure field are coarse approximations of the derivatives, estimates of the mean forcing fields Πu and Πυ were noisy and uncertain. Estimates of σu and συ were more robust. Upper and lower limits of the estimates of these parameters are listed in Table 2. For the sake of simplicity, the effect of downward mixing of momentum was ignored by setting
K
and a constant drag coefficient was assumed:
i1520-0469-61-16-2037-e39
where the layer depth h = 80 m corresponds to the average thickness of the lowest σ level of the NCEP– NCAR reanalysis model. The numerical integration of Eq. (31) was carried out using a forward-Euler scheme, with parameters Πu, Πυ, σu, and συ sampled evenly through the ranges indicated in Table 2. The numerical integration of stochastic differential equations is described in Kloeden and Platen (1992); in brief, for a timestep δt, the forward-Euler numerical recursion is
i1520-0469-61-16-2037-e41
where ξk is a vector of mutually independent zero-mean, unit-variance Gaussian random variables generated independently at each time step:
ξiξTjδij
where 𝗜 is the 4 × 4 identity matrix.

Figure 5 displays scatterplots of the mean versus the skewness of the simulated wind components for different parameter values. Clearly, the mean and skewness of the simulated wind components are anticorrelated, with a linear dependence that matches favorably that illustrated in Fig. 3. Note that the skewness is a result of nonlinearities in the drift term, and not of multiplicative noise as suggested by Sura (2003); by assumption, the intensities of fluctuations in Πu and Πυ in this model are independent of the state of the boundary layer.

In this analysis, we have neglected the dependence of the neutral drag coefficient on the wind speed. Several numerical integrations of Eq. (31) were carried out in which Γ was taken to depend on |u| using various formulations described in Taylor (2000) and Fairall et al. (2003). In none of these integrations was the relationship between the mean and skewness fields of the surface wind components found to be significantly different from that displayed in Fig. 5 for a drag coefficient independent of wind speed. The surface wind speed dependence of the neutral drag coefficient does not appear to play a leading-order role in the relationship between the mean and skewness fields of u and υ.

An intuitive understanding of the relationship between the mean and skewness of the surface wind components is straightforward. Consider the atmospheric surface layer subject to zonal forcing Πu with zero mean and with fluctuations that are equally as likely to be positive as negative. By symmetry, the mean and skewness of the resulting zonal wind will both be zero. Now suppose that the zonal forcing Πu has a non-zero mean, which we will take (without loss of generality) to be positive, so that the mean u will be positive. Because of the nonlinear surface drag, positive anomalies in u will be subject to stronger friction than negative anomalies. Consequently, a positive perturbation in the forcing will produce a weaker response than a negative perturbation of the same magnitude, and the symmetric fluctuations will produce an asymmetric response. In particular, a tail toward slower zonal winds will develop in the distribution, so the zonal wind will be negatively skewed. As the mean of Πu increases, the asymmetry in drag between positive and negative u anomalies will also increase, and so the skewness will become more strongly negative. By the same token, negative mean Πu will yield negative mean zonal winds that are positively skewed, with the skewness becoming more positive as the mean of u becomes more negative. In this manner, the nonlinear surface drag results in the anticorrelation of the mean and skewness fields of zonal surface winds. A corollary of this argument is that for a given mean zonal forcing, the skewness of u should decrease as the variance of u decreases, because the difference in drag between positive and negative fluctuations of typical (e.g., unit standard deviation) magnitude will decease. In fact, this behavior is observed in Fig. 5. The scatter of points with relatively low skewness for relatively high wind speed correspond to relatively low standard deviations. This effect explains the greater scatter in skewness for easterly zonal winds than for westerly zonal winds evident in Fig. 3. In some parts of the Tropics (e.g., the eastern tropical Pacific), surface zonal winds are subject to mean forcings of about the same magnitude as surface winds in the midlatitudes, but of considerably weaker variability. Consequently, these are regions of relatively weak skewness. The preceding arguments hold of course for meridional winds as well.

It is admittedly naive to assume that Πu, Πυ, Γ, and K should be independent Ornstein–Uhlenbeck processes. First, it cannot be expected that Πu, Πυ, Γ, and K should have Gaussian probability distributions and exponential temporal autocorrelation functions. In particular, Πu and Πυ estimated from the NCEP–NCAR reanalyses are themselves non-Gaussian. However, the spatial distribution of the skewness of u bears no obvious relation to that of Πu (and similarly for υ and Πυ), so the skewness of the wind components is not simply a result of skewed forcing fluctuations. Furthermore, the Obukhov length L and the eddy viscosity K are to some extent determined by the boundary layer wind speed, which itself is affected by pressure gradient forces. Consequently, fluctuations in Γ and K should not be independent of those in Πu and Πυ. The choice to represent these fluctuations by independent Ornstein– Uhlenbeck processes was not based on any physical motivation, but on the desire to have this representation be as simple as possible. The simplicity of the system (31) and its success in characterizing the relationship between the mean and skewness of surface wind components indicates that the approximations lead to a model that is qualitatively useful.

The strong linear relationship between the mean and skewness of the surface wind components has been shown to be a consequence of the nonlinearity of the surface drag, which follows from well-established similarity theories for turbulent boundary layers. A number of previous studies of the momentum balance of monthly averaged surface winds (e.g., Deser 1993; Ward and Hoskins 1996; Chiang and Zebiak 2000) have employed a surface drag law in which the surface stress varies linearly with wind speed to represent the combined effects of surface friction and downward mixing of momentum into the surface layer. While such a linear drag law may be appropriate for the equilibrium balance of spatially and temporally averaged winds, it is clearly inadequate for understanding the character of nonequilibrium fluctuations in the winds on the daily and shorter time scales that play important roles in turbulent transfers between the ocean and the atmosphere.

5. Stochastic dynamics of sea surface winds

A recent effort to characterize the variability of sea surface winds is the study of Sura (2003), in which the drift and diffusion functions of SDEs for the zonal and meridional winds were estimated empirically from 6-hourly blended QuikSCAT–NCEP analysis data. Sura assumed that the zonal and meridional components of the wind can be described individually by one-dimensional stochastic differential equations, that is, that the 6-hourly zonal and meridional winds are individually Markov processes. However, the discussion in the previous section suggests that this assumption is inconsistent with dynamical balances in the surface layer when the forcings Πu and Πυ decorrelate on time scales of the same order as the decorrelation time scale of the winds. Given that marine boundary layers evolve slowly relative to terrestrial boundary layers, it is also to be expected that the decorrelation time scales of Γ and K will not be much less than 6 h. In fact, Rieder (1997) found that changes in the drag coefficient follow changes in surface wind speed with a lag of about 4 h, suggesting that the turbulent marine boundary layer has an adjustment scale on the order of several hours. Furthermore, u and υ are coupled frictionally, and even in the limit of white-noise fluctuations in Πu, Πυ, Γ, and K, neither u nor υ is individually Markov. Nevertheless, it is instructive to consider the limit of Eq. (31) in which these fluctuations become white relative to those of the winds, and to compare the resulting system to the results of the empirical analysis of Sura (2003). For the sake of simplicity, we will neglect the dependence of the drag coefficient on the wind speed.

In the limit that τu, τυ, τγ, τμ → 0, Πu, Πυ, γ, and μ become white-noise processes, as described in section 4. By the Wong–Zakai theorem (Penland 2003b), the six-dimensional stochastic differential equation (31) converges to a stochastic differential equation for the two-dimensional process u = (u, υ):
CUU
where
i1520-0469-61-16-2037-e44
Note that in this limit the SDE involves products of state variables and white-noise processes; the corresponding noise terms are said to be multiplicative (Gardiner 1997). The open circle in Eq. (43) indicates that the multiplicative noise terms are to be interpreted in the Stratonovich sense, as is appropriate for the white-noise limit of a process with a non-zero autocorrelation time (Gardiner 1997; Penland 2003b). The convergence of Eq. (31) to (43) is a special case of the general convergence problem considered in Majda et al. (2001, 2003) for the simple situation in which the dynamics of the “fast” variables (Πu, Πυ, γ, μ) do not depend on the “slow” variables u and υ (although of course the dynamics of the slow variables depend on the fast variables). In the present simple case, the theorem of Wong and Zakai is sufficient to characterize the white-noise limit of the SDE (31).
The equivalent Itô drift
i1520-0469-61-16-2037-e46
is obtained through the standard transformation
i1520-0469-61-16-2037-e47
where ∂k for k = 1 and 2 denotes partial differentiation with respect to u and υ, respectively. The difference between Itô and Stratonovich representations of SDEs arises in the presence of multiplicative noise, for which the state dependence of fluctuations can change the average dynamics of the system. Penland (2003a) presents an intuitive description of the origins of this “noise-induced drift;” technical details are given in Penland (2003b), Gardiner (1997), and Kloeden and Platen (1992). We need to consider the Itô drift CI(U) because, as we will see, the drift and diffusion functions calculated by the algorithm used in Sura (2003) are those of an Itô SDE.
In the limit that fluctuations in υ are significantly smaller than those in u, to a good approximation u itself can be described by a Markov process described by the Itô SDE
i1520-0469-61-16-2037-e48
with the associated Fokker–Planck equation for the probability density function p(u):
i1520-0469-61-16-2037-e49
where
i1520-0469-61-16-2037-e50
Sura (2003) assumes that the zonal wind u is a Markov process governed by an Itô SDE of the form
αuβuẆ,
where it has been assumed that only a single white-noise process enters the SDE, so the diffusion β is a scalar function of u. This equation is similar in form to the scalar equation (48) for the evolution of u, which we have shown follows from the momentum equations in a turbulent boundary layer in the limit of white-noise forcing and vanishing variability in υ. Equations (48) and (52) differ, however, in that three independent white-noise processes are present in (48) so the diffusion is a 1 × 3 matrix.
The algorithm used in Sura (2003) to estimate the parameters of the SDE (52) for the zonal wind takes advantage of the fact that the Itô SDE
Eyy
(where y is an n-dimensional vector, is an m-dimensional white-noise vector, and 𝗙(y) is an n × m matrix function of y) is associated with the Fokker–Planck equation for the joint probability density function p(y):
i1520-0469-61-16-2037-e54
and thus that the product 𝗙𝗙T can be estimated from the limit
i1520-0469-61-16-2037-e55
Thus, in Sura's methodology, the diffusion matrix 𝗙 is not estimated directly, but only through the product 𝗙𝗙T. Sura notes the nonuniqueness of 𝗙 satisfying Eq. (55) for n > 1, and asserts that this nonuniqueness vanishes (up to a sign) for the scalar case, n = 1. This is true if m = 1, that is, if only a single white-noise process enters the SDE for the scalar process, but is certainly not true if m > 1, which the preceding analysis suggests is the case relevant to the stochastic dynamics of sea surface winds. In fact, the estimation algorithm presented in Sura (2003) does not provide a methodology to estimate the number m of white-noise processes entering the SDE (53), and thus the matrix 𝗙 cannot be uniquely determined.
Note that an SDE with m = 1 and diffusion
i1520-0469-61-16-2037-e56
which is approximately quadratic in u, would produce the same Fokker–Planck equation as the SDE (48). Consequently, Δ(u) would be the diffusion estimated by the method of Sura (2003) from a dataset produced by the SDE (48). Figure 6 shows plots of the drift CI1 and effective diffusion Δ as functions of u, calculated using characteristic values Πu = 0.3 × 10−4 m s−2, σu = 0.15 m s−3/2, K = 0, υ = 3 m s−1, σγ = 30 s1/2 (corresponding to a 20% standard deviation for fluctuations of Γ with a decorrelation time scale of 3 h). The drift and effective diffusion functions displayed in Fig. 6 bear a striking resemblance to those estimated in Sura (2003). In particular, the effective diffusion varies approximately quadratically with u and the drift function decreases monotonically with u and displays inflection points on either side of the zonal wind value ueq which balances the pressure gradient, Coriolis, and frictional forces. The physical interpretation of the drift and diffusion coefficients shown in Fig. 6 is that the “deterministic” component of the dynamics acts to drive the zonal wind to the balanced value ueq while the “stochastic” component drives the zonal wind away from this equilibrium value, such that typical fluctuations become larger with increasing |u|. Interestingly, the diffusion functions estimated in Sura (2003) have minima at nonzero values of u; in contrast, the minimum of the effective diffusion term Δ(u) is necessarily at u = 0. This discrepancy may be a consequence of the simplifying assumptions we have used to obtain the SDE (48), or it may be an artifact of Sura's estimation algorithm. In particular, the removal of the seasonal cycle from the winds in Sura (2003) before estimating the drift and diffusion coefficients could lead to biases in the estimated drift and diffusion functions, as the drag is a function of the total wind speed and not its standardized anomaly. Nevertheless, while the agreement between the drift and diffusion functions estimated by Sura (2003) and those arising from boundary layer dynamics is not perfect, it is remarkably good.

Sura (2003) attributes the quadratic dependence of Δ(u) on u to burstiness in the turbulent boundary layer, suggesting that as midlatitude wind speeds increase, so too does the variability. In fact, there is no evidence of this behavior in the 6-hourly winds. Figure 7 displays a scatterplot of the absolute value of the 6-hourly zonal wind increment, |u(t + 6 h) − u(t)|, as a function of [u(t + 6 h) + u(t)]/2 for winds at the representative point (57°S, 0°). It is evident that instead of the 6-h wind increments being larger at larger wind values, they are generally smaller. Inspection of similar plots at different grid points indicates that this relationship is independent of latitude and longitude. The increase in wind variability at higher wind speeds may be a feature of the winds at shorter time scales, but there is no evidence of this behavior at the 6-hourly time scales used in Sura (2003) to estimate the SDE parameters. The previous analysis suggests that multiplicative noise in the τγ → 0 limit arises not from burstiness in the winds, but from turbulent fluctuations in the magnitude of the drag coefficient.

6. Conclusions

This study has demonstrated the existence of a strong spatial anticorrelation between the mean and skewness fields of surface wind components over the ocean, and has shown that this relationship can be understood as a consequence of the nonlinear surface drag predicted by Monin–Obukhov similarity theory for a turbulent boundary layer. As this relationship is evident in the NCEP–NCAR and ERA-40 reanalysis products, in blended QuikSCAT scatterometer–NCEP analysis data, in SSM/I satellite data, and in TAO buoy data, it is not simply an artifact of the observing system or the reanalysis algorithm. Furthermore, it has been demonstrated that the surface layer momentum dynamics subject to fluctuation forcing are consistent with the results of the empirical study by Sura (2003), although the interpretation of the results is different. In particular, it is suggested that the presence of multiplicative noise does not arise from increasing gustiness with higher wind speeds, as suggested by Sura, but from fluctuations in the turbulent drag.

A number of approximations were used in this study that cannot be expected to hold in the real world. First, the downward mixing of momentum into the surface layer, which has been ignored in this study, has been demonstrated to be a potentially significant component of the surface layer momentum budget in the Tropics (e.g., Deser 1993; Chiang and Zebiak 2000). Second, advection terms in the momentum budget have been ignored; while this may be justifiable for the equilibrium momentum budget of monthly averaged winds (Deser 1993), it may be important during periods of strong winds on shorter time scales. Third, Monin–Obukhov similarity theory predicts that buoyancy effects in a nonneutral boundary layer introduce Richardson number– dependent correction terms to the drag coefficient. Fourth, it has been assumed that all fluctuating quantities in the momentum equations (Πu, Πυ, Γ, and K) can be represented as independent Gaussian random variables; as is described in section 4, there are compelling reasons why this should not be true. More complex models of boundary layer turbulence could be used to develop a better understanding of the fluctuations in the drag coefficient and the eddy viscosity, and the relationship of these distributions to the surface winds. Such an extension of the present study would be necessary for the detailed quantitative analysis of the relationships between the moments of sea surface wind components. Finally, the statistics of Πu, Πυ, γ, and μ have been assumed to be stationary, when in fact they will display both diurnal and annual variability (Deser 1994). Nonetheless, the model presented in this study is able to reproduce the essential features of the relationship between the mean and skewness of surface winds, and is broadly consistent with the empirical results of Sura (2003). An investigation of the reasons for those differences that do exist between the results of the present study and those of Sura (2003) are an interesting direction of future study. Thus, while the simplifying approximations may limit the model's quantitative accuracy, they do not seem to have significantly degraded its qualitative utility.

Previous studies of the oceanic response to fluctuating atmospheric forcing (e.g., Moore 1999; Sura et al. 2001; Monahan 2002a; Kuhlbrodt and Monahan 2003) have assumed Gaussian fluctuations, while the results presented in the present study demonstrate that in the Tropics and midlatitudes the surface wind components are generally skewed. In particular, a leading theory of El Niño–Southern Oscillation variability posits that the tropical Pacific coupled atmosphere–ocean system acts as a damped linear oscillator subject to Gaussian perturbations (e.g., Penland 1996; Kleeman and Moore 1997; Thompson and Battisti 2000). Such linear models cannot account for the marked asymmetry of the magnitude and structure of sea surface temperature anomalies between El Niño and La Niña events that has been noted in observations (e.g., Monahan 2001; Hannachi et al. 2003; An and Jin 2004; Monahan and Dai 2004). This asymmetry is generally ascribed to nonlinearities in ocean dynamics, but may also reflect the strong positive skewness of the zonal wind in the western equatorial Pacific. An interesting question is the extent to which the strength of El Niño events relative to La Niña events simply reflects the fact that westerly wind bursts in the western equatorial Pacific are not matched by easterly wind bursts.

A natural extension of the present study would be an investigation of the dependence on time scale of the relationship between the mean and skewness of surface wind components (e.g., the relationship between monthly averaged winds and the skewness of submonthly fluctuations), and the incorporation of this information into eddy flux closure schemes. In particular, such a study could provide important constraints on the development of stochastic parameterization schemes, which have been proposed as improvements to the representation of subgrid-scale phenomena in models of the climate system (Palmer 2001; Imkeller and Monahan 2002). Errors in the representation of turbulent fluxes between the ocean and the atmosphere are a potentially significant source of systematic errors in global climate models, and the development of more accurate parameterization schemes is an important aspect of improving our understanding of this complex coupled system.

Acknowledgments

The author acknowledges support from the Natural Sciences and Engineering Research Council of Canada, by the Canadian Foundation for Climate and Atmospheric Sciences, and by the Canadian Institute for Advanced Research Earth System Evolution Program. The author is grateful to Norm McFarlane, Ken Denman, and Alexandra Guerrero-Martinez for helpful discussions. The author would like to thank the three anonymous referees whose thoughtful comments considerably improved this manuscript.

REFERENCES

  • Aeberhardt, M., M. Blatter, and T. F. Stocker, 2000: Variability on the century time scale and regime changes in a stochastically forced zonally averaged ocean–atmosphere model. Geophys. Res. Lett, 27 , 13031306.

    • Search Google Scholar
    • Export Citation
  • An, S-I., and F-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17 , 23992412.

  • Atlas, R., R. Hoffman, S. Bloom, J. Jusem, and J. Ardizzone, 1996: A multiyear global surface wind velocity dataset using SSM/I wind observations. Bull. Amer. Meteor. Soc, 77 , 869882.

    • Search Google Scholar
    • Export Citation
  • Bates, N. R., and L. Merlivat, 2001: The influence of short-term wind variability on air–sea CO2 exchange. Geophys. Res. Lett, 28 , 32813284.

    • Search Google Scholar
    • Export Citation
  • Chiang, J. C., and S. E. Zebiak, 2000: Surface wind over tropical oceans: Diagnosis of the momentum balance, and modeling the linear friction coefficient. J. Climate, 13 , 17331747.

    • Search Google Scholar
    • Export Citation
  • Chin, T., R. Milliff, and W. Large, 1998: Basin-scale, high-wavenumber sea surface wind fields from a multiresolution analysis of scatterometer data. J. Atmos. Oceanic Technol, 15 , 741763.

    • Search Google Scholar
    • Export Citation
  • Deser, C., 1993: Diagnosis of the surface momentum balance over the tropical Pacific Ocean. J. Climate, 6 , 6474.

  • Deser, C., 1994: Daily surface wind variations over the equatorial Pacific Ocean. J. Geophys. Res, 99 (D11) 2307123078.

  • Fairall, C. W., E. F. Bradley, J. E. Hare, A. A. Grachev, and J. B. Edson, 2003: Bulk parameterization of air–sea fluxes: Updates and verification for the COARE algorithm. J. Climate, 16 , 571591.

    • Search Google Scholar
    • Export Citation
  • Gardiner, C. W., 1997: Handbook of Stochastic Methods for Physics, Chemistry, and the Natural Sciences. Springer, 442 pp.

  • Grachev, A. A., and C. W. Fairall, 2001: Upward momentum transfer in the marine boundary layer. J. Phys. Oceanogr, 31 , 16981711.

  • Haltiner, G. J., and R. T. Williams, 1980: Numerical Prediction and Dynamic Meteorology. Wiley, 477 pp.

  • Hannachi, A., D. Stephenson, and K. Sperber, 2003: Probability-based methods for quantifying nonlinearity in the ENSO. Climate Dyn, 20 , 241256.

    • Search Google Scholar
    • Export Citation
  • Imkeller, P., and A. H. Monahan, 2002: Conceptual stochastic climate models. Stochast. Dyn, 2 , 311326.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc, 77 , 437471.

  • Kleeman, R., and A. Moore, 1997: A theory for the limitation of ENSO predictability due to stochastic atmospheric transients. J. Atmos. Sci, 54 , 753767.

    • Search Google Scholar
    • Export Citation
  • Kloeden, P. E., and E. Platen, 1992: Numerical Solution of Stochastic Differential Equations. Springer-Verlag, 632 pp.

  • Kuhlbrodt, T., and A. H. Monahan, 2003: Stochastic stability of open-ocean deep convection. J. Phys. Oceanogr, 33 , 27642780.

  • Kuhlbrodt, T., S. Titz, U. Feudel, and S. Rahmstorf, 2001: A simple model of seasonal open ocean convection. Part II: Labrador Sea stability and stochastic forcing. Ocean Dyn, 52 , 3649.

    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, J. Howell, J. Højstrup, J. M. Wilczak, J. Edson, and J. Hare, 1996: Sea surface drag coefficients in the Risø Air Sea Experiment. J. Geophys. Res, 101 , 1432714335.

    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, P. Frederickson, K. Davidson, and A-S. Smedman, 2003: Sea-surface aerodynamic roughness. J. Geophys. Res.,108, 3171, doi:10.1029/2002JC001383.

    • Search Google Scholar
    • Export Citation
  • Majda, A. J., I. Timofeyev, and E. Vanden-Eijnden, 2001: A mathematical framework for stochastic climate models. Commun. Pure Appl. Math, 54 , 891974.

    • Search Google Scholar
    • Export Citation
  • Majda, A. J., I. Timofeyev, and E. Vanden-Eijnden, 2003: Systematic strategies for stochastic mode reduction in climate. J. Atmos. Sci, 60 , 17051722.

    • Search Google Scholar
    • Export Citation
  • McPhaden, M., and Coauthors, 1998: The Tropical Ocean–Global Atmosphere (TOGA) observing system: A decade of progress. J. Geophys. Res, 103 , 1416914240.

    • Search Google Scholar
    • Export Citation
  • Milliff, R., W. Large, J. Morzel, and G. Danabasoglu, 1999: Ocean general circulation model sensitivity to forcing from scatterometer winds. J. Geophys. Res, 104C , 1133711358.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2001: Nonlinear principal component analysis: Tropical Indo–Pacific sea surface temperature and sea level pressure. J. Climate, 14 , 219233.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2002a: Lyapunov exponents of a simple stochastic model of the thermally and wind-driven ocean circulation. Dyn. Atmos. Oceans, 35 , 363388.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2002b: Stabilization of climate regimes by noise in a simple model of the thermohaline circulation. J. Phys. Oceanogr, 32 , 20722085.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., and A. Dai, 2004: The spatial and temporal structure of ENSO nonlinearity. J. Climate, 17 , 30263036.

  • Monahan, A. H., and K. L. Denman, 2004: Impacts of atmospheric variability on a coupled upper-ocean/ecosystem model of the subarctic Northeast Pacific. Global Biogeochem. Cycles, in press.

    • Search Google Scholar
    • Export Citation
  • Moore, A. M., 1999: Wind-induced variability of ocean gyres. Dyn. Atmos. Ocean, 29 , 335364.

  • Palmer, T. N., 2001: A nonlinear dynamical perspective on model error: A proposal for nonlocal stochastic–dynamic parameterisation in weather and climate prediction models. Quart. J. Roy. Meteor. Soc, 127 , 279304.

    • Search Google Scholar
    • Export Citation
  • Penland, C., 1996: A stochastic model of IndoPacific sea surface temperature anomalies. Physica D, 98 , 534558.

  • Penland, C., 2003a: Noise out of chaos and why it won't go away. Bull. Amer. Meteor. Soc, 84 , 921925.

  • Penland, C., 2003b: A stochastic approach to nonlinear dynamics: A review. Bull. Amer. Meteor. Soc, 84 , ES43ES52.

  • Ponte, R. M., and R. D. Rosen, 2004: Nonlinear effects of variable winds on ocean stress climatologies. J. Climate, 17 , 12831293.

  • Rieder, K. F., 1997: Analysis of sea-surface drag parameterizations in open ocean conditions. Bound.-Layer Meteor, 82 , 355377.

  • Silverman, B. W., 1986: Density Estimation for Statistics and Data Analysis. Chapman and Hall, 175 pp.

  • Simmons, A., and J. Gibson, 2000: The ERA-40 Project Plan. ERA-40 Project Report Series No. 1, ECMWF, Reading, United Kingdom, 63 pp.

  • Stull, R. B., 1997: An Introduction to Boundary Layer Meteorology. Kluwer, 670 pp.

  • Sura, P., 2003: Stochastic analysis of Southern and Pacific Ocean sea surface winds. J. Atmos. Sci, 60 , 654666.

  • Sura, P., K. Fraedrich, and F. Lunkheit, 2001: Regime transitions in a stochastically forced double-gyre model. J. Phys. Oceanogr, 31 , 411426.

    • Search Google Scholar
    • Export Citation
  • Taylor, P. K., Ed.,. 2000: Intercomparison and validation of ocean– atmosphere energy flux fields. Joint WCRP/SCOR Working Group on Air–Sea Fluxes Final Report, WMO Tech. Doc. 1036, 306 pp.

    • Search Google Scholar
    • Export Citation
  • Taylor, P. K., and M. J. Yelland, 2001: The dependence of sea surface roughness on the height and steepness of the waves. J. Phys. Oceanogr, 31 , 572590.

    • Search Google Scholar
    • Export Citation
  • Thompson, C., and D. Battisti, 2000: A linear stochastic dynamical model of ENSO. Part I: Development. J. Climate, 13 , 28182832.

  • Ward, M. N., and B. J. Hoskins, 1996: Near-surface wind over the global ocean 1949–1988. J. Climate, 9 , 18771895.

  • Wright, D. G., and K. R. Thompson, 1983: Time-averaged forms of the nonlinear stress law. J. Phys. Oceanogr, 13 , 341345.

Fig. 1.
Fig. 1.

Mean, std dev, and skewness fields of NCEP–NCAR reanalysis zonal winds

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 2.
Fig. 2.

Mean, std dev, and skewness fields of NCEP–NCAR reanalysis meridional winds

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 3.
Fig. 3.

Grid point by grid point scatterplots of the mean vs the skewness of the zonal and meridional sea surface wind components for the NCEP–NCAR reanalysis data (red dots), the blended scatterometer–NCEP analysis data (black dots), and the TAO buoy data (green circles)

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 4.
Fig. 4.

Kernel density estimates of the probability density functions of the zonal wind, u, at 50°S, 110°W; 20°N, 160°W; and 50°N, 60°W

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 5.
Fig. 5.

Scatterplots of the mean vs the skewness of the zonal and meridional components of the simulated wind components.

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 6.
Fig. 6.

Drift and diffusion functions corresponding to the system Eq. (48) that would be estimated using the methodology of Sura (2003)

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Fig. 7.
Fig. 7.

Scatterplot of the absolute 6-h increment of the zonal wind, |u(t + 6 h) − u(t)|, as a function of the zonal wind strength [u(t + 6 h) + u(t)]/2 at 57°S, 0°

Citation: Journal of the Atmospheric Sciences 61, 16; 10.1175/1520-0469(2004)061<2037:ASMFTS>2.0.CO;2

Table 1.

Spatial correlation coefficients between the mean and skewness fields of the zonal and meridional winds from the NCEP– NCAR reanalysis data, the ERA-40 reanalysis, the blended Quik SCAT–NCEP analysis data, the SSM/I passive microwave data, and the TAO buoy data

Table 1.
Table 2.

Stochastic boundary layer model parameter ranges estimated from NCEP–NCAR reanalysis data

Table 2.
Save
  • Aeberhardt, M., M. Blatter, and T. F. Stocker, 2000: Variability on the century time scale and regime changes in a stochastically forced zonally averaged ocean–atmosphere model. Geophys. Res. Lett, 27 , 13031306.

    • Search Google Scholar
    • Export Citation
  • An, S-I., and F-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17 , 23992412.

  • Atlas, R., R. Hoffman, S. Bloom, J. Jusem, and J. Ardizzone, 1996: A multiyear global surface wind velocity dataset using SSM/I wind observations. Bull. Amer. Meteor. Soc, 77 , 869882.

    • Search Google Scholar
    • Export Citation
  • Bates, N. R., and L. Merlivat, 2001: The influence of short-term wind variability on air–sea CO2 exchange. Geophys. Res. Lett, 28 , 32813284.

    • Search Google Scholar
    • Export Citation
  • Chiang, J. C., and S. E. Zebiak, 2000: Surface wind over tropical oceans: Diagnosis of the momentum balance, and modeling the linear friction coefficient. J. Climate, 13 , 17331747.

    • Search Google Scholar
    • Export Citation
  • Chin, T., R. Milliff, and W. Large, 1998: Basin-scale, high-wavenumber sea surface wind fields from a multiresolution analysis of scatterometer data. J. Atmos. Oceanic Technol, 15 , 741763.

    • Search Google Scholar
    • Export Citation
  • Deser, C., 1993: Diagnosis of the surface momentum balance over the tropical Pacific Ocean. J. Climate, 6 , 6474.

  • Deser, C., 1994: Daily surface wind variations over the equatorial Pacific Ocean. J. Geophys. Res, 99 (D11) 2307123078.

  • Fairall, C. W., E. F. Bradley, J. E. Hare, A. A. Grachev, and J. B. Edson, 2003: Bulk parameterization of air–sea fluxes: Updates and verification for the COARE algorithm. J. Climate, 16 , 571591.

    • Search Google Scholar
    • Export Citation
  • Gardiner, C. W., 1997: Handbook of Stochastic Methods for Physics, Chemistry, and the Natural Sciences. Springer, 442 pp.

  • Grachev, A. A., and C. W. Fairall, 2001: Upward momentum transfer in the marine boundary layer. J. Phys. Oceanogr, 31 , 16981711.

  • Haltiner, G. J., and R. T. Williams, 1980: Numerical Prediction and Dynamic Meteorology. Wiley, 477 pp.

  • Hannachi, A., D. Stephenson, and K. Sperber, 2003: Probability-based methods for quantifying nonlinearity in the ENSO. Climate Dyn, 20 , 241256.

    • Search Google Scholar
    • Export Citation
  • Imkeller, P., and A. H. Monahan, 2002: Conceptual stochastic climate models. Stochast. Dyn, 2 , 311326.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc, 77 , 437471.

  • Kleeman, R., and A. Moore, 1997: A theory for the limitation of ENSO predictability due to stochastic atmospheric transients. J. Atmos. Sci, 54 , 753767.

    • Search Google Scholar
    • Export Citation
  • Kloeden, P. E., and E. Platen, 1992: Numerical Solution of Stochastic Differential Equations. Springer-Verlag, 632 pp.

  • Kuhlbrodt, T., and A. H. Monahan, 2003: Stochastic stability of open-ocean deep convection. J. Phys. Oceanogr, 33 , 27642780.

  • Kuhlbrodt, T., S. Titz, U. Feudel, and S. Rahmstorf, 2001: A simple model of seasonal open ocean convection. Part II: Labrador Sea stability and stochastic forcing. Ocean Dyn, 52 , 3649.

    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, J. Howell, J. Højstrup, J. M. Wilczak, J. Edson, and J. Hare, 1996: Sea surface drag coefficients in the Risø Air Sea Experiment. J. Geophys. Res, 101 , 1432714335.

    • Search Google Scholar
    • Export Citation
  • Mahrt, L., D. Vickers, P. Frederickson, K. Davidson, and A-S. Smedman, 2003: Sea-surface aerodynamic roughness. J. Geophys. Res.,108, 3171, doi:10.1029/2002JC001383.

    • Search Google Scholar
    • Export Citation
  • Majda, A. J., I. Timofeyev, and E. Vanden-Eijnden, 2001: A mathematical framework for stochastic climate models. Commun. Pure Appl. Math, 54 , 891974.

    • Search Google Scholar
    • Export Citation
  • Majda, A. J., I. Timofeyev, and E. Vanden-Eijnden, 2003: Systematic strategies for stochastic mode reduction in climate. J. Atmos. Sci, 60 , 17051722.

    • Search Google Scholar
    • Export Citation
  • McPhaden, M., and Coauthors, 1998: The Tropical Ocean–Global Atmosphere (TOGA) observing system: A decade of progress. J. Geophys. Res, 103 , 1416914240.

    • Search Google Scholar
    • Export Citation
  • Milliff, R., W. Large, J. Morzel, and G. Danabasoglu, 1999: Ocean general circulation model sensitivity to forcing from scatterometer winds. J. Geophys. Res, 104C , 1133711358.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2001: Nonlinear principal component analysis: Tropical Indo–Pacific sea surface temperature and sea level pressure. J. Climate, 14 , 219233.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2002a: Lyapunov exponents of a simple stochastic model of the thermally and wind-driven ocean circulation. Dyn. Atmos. Oceans, 35 , 363388.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., 2002b: Stabilization of climate regimes by noise in a simple model of the thermohaline circulation. J. Phys. Oceanogr, 32 , 20722085.

    • Search Google Scholar
    • Export Citation
  • Monahan, A. H., and A. Dai, 2004: The spatial and temporal structure of ENSO nonlinearity. J. Climate, 17 , 30263036.

  • Monahan, A. H., and K. L. Denman, 2004: Impacts of atmospheric variability on a coupled upper-ocean/ecosystem model of the subarctic Northeast Pacific. Global Biogeochem. Cycles, in press.

    • Search Google Scholar
    • Export Citation
  • Moore, A. M., 1999: Wind-induced variability of ocean gyres. Dyn. Atmos. Ocean, 29 , 335364.

  • Palmer, T. N., 2001: A nonlinear dynamical perspective on model error: A proposal for nonlocal stochastic–dynamic parameterisation in weather and climate prediction models. Quart. J. Roy. Meteor. Soc, 127 , 279304.

    • Search Google Scholar
    • Export Citation
  • Penland, C., 1996: A stochastic model of IndoPacific sea surface temperature anomalies. Physica D, 98 , 534558.

  • Penland, C., 2003a: Noise out of chaos and why it won't go away. Bull. Amer. Meteor. Soc, 84 , 921925.

  • Penland, C., 2003b: A stochastic approach to nonlinear dynamics: A review. Bull. Amer. Meteor. Soc, 84 , ES43ES52.

  • Ponte, R. M., and R. D. Rosen, 2004: Nonlinear effects of variable winds on ocean stress climatologies. J. Climate, 17 , 12831293.

  • Rieder, K. F., 1997: Analysis of sea-surface drag parameterizations in open ocean conditions. Bound.-Layer Meteor, 82 , 355377.

  • Silverman, B. W., 1986: Density Estimation for Statistics and Data Analysis. Chapman and Hall, 175 pp.

  • Simmons, A., and J. Gibson, 2000: The ERA-40 Project Plan. ERA-40 Project Report Series No. 1, ECMWF, Reading, United Kingdom, 63 pp.

  • Stull, R. B., 1997: An Introduction to Boundary Layer Meteorology. Kluwer, 670 pp.

  • Sura, P., 2003: Stochastic analysis of Southern and Pacific Ocean sea surface winds. J. Atmos. Sci, 60 , 654666.

  • Sura, P., K. Fraedrich, and F. Lunkheit, 2001: Regime transitions in a stochastically forced double-gyre model. J. Phys. Oceanogr, 31 , 411426.

    • Search Google Scholar
    • Export Citation
  • Taylor, P. K., Ed.,. 2000: Intercomparison and validation of ocean– atmosphere energy flux fields. Joint WCRP/SCOR Working Group on Air–Sea Fluxes Final Report, WMO Tech. Doc. 1036, 306 pp.

    • Search Google Scholar
    • Export Citation
  • Taylor, P. K., and M. J. Yelland, 2001: The dependence of sea surface roughness on the height and steepness of the waves. J. Phys. Oceanogr, 31 , 572590.

    • Search Google Scholar
    • Export Citation
  • Thompson, C., and D. Battisti, 2000: A linear stochastic dynamical model of ENSO. Part I: Development. J. Climate, 13 , 28182832.

  • Ward, M. N., and B. J. Hoskins, 1996: Near-surface wind over the global ocean 1949–1988. J. Climate, 9 , 18771895.

  • Wright, D. G., and K. R. Thompson, 1983: Time-averaged forms of the nonlinear stress law. J. Phys. Oceanogr, 13 , 341345.

  • Fig. 1.

    Mean, std dev, and skewness fields of NCEP–NCAR reanalysis zonal winds

  • Fig. 2.

    Mean, std dev, and skewness fields of NCEP–NCAR reanalysis meridional winds

  • Fig. 3.

    Grid point by grid point scatterplots of the mean vs the skewness of the zonal and meridional sea surface wind components for the NCEP–NCAR reanalysis data (red dots), the blended scatterometer–NCEP analysis data (black dots), and the TAO buoy data (green circles)

  • Fig. 4.

    Kernel density estimates of the probability density functions of the zonal wind, u, at 50°S, 110°W; 20°N, 160°W; and 50°N, 60°W

  • Fig. 5.

    Scatterplots of the mean vs the skewness of the zonal and meridional components of the simulated wind components.

  • Fig. 6.

    Drift and diffusion functions corresponding to the system Eq. (48) that would be estimated using the methodology of Sura (2003)

  • Fig. 7.

    Scatterplot of the absolute 6-h increment of the zonal wind, |u(t + 6 h) − u(t)|, as a function of the zonal wind strength [u(t + 6 h) + u(t)]/2 at 57°S, 0°

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 499 230 57
PDF Downloads 183 50 8