## 1. Introduction

Turbulence research in the atmosphere has built on the groundwork laid by fluid mechanics research in the laboratory. But, atmospheric research may also have suffered by assuming too much similarity with laboratory fluid mechanics. As an example, we consider here drag parameterizations at the air–sea interface.

*P*is the fluid pressure,

*ρ*is the fluid density,

*U*is the flow speed,

*g*is the acceleration of gravity, and

*z*is the height above some arbitrary reference level. This equation essentially states that the energy per unit volume is constant along a streamline in a fluid flow.

*A*is the frontal area of the obstacle and

*C*is its drag coefficient. In the context of (1.1),

_{D}*D*/

*A*can be thought of as the change in the constant along a streamline that results from frictional losses to the obstacle; the ½ in (1.2) emphasizes the concept that the drag is a change in kinetic energy per unit volume in the fluid. Because energy can be changed only by adding or subtracting momentum,

*D*/

*A*can also be thought of as a momentum flux.

*τ*is the drag per unit area of sea surface (also called the surface stress or the momentum flux),

*ρ*is the air density,

_{a}*U*is the wind speed at some reference height

_{r}*r*above the sea, and

*C*is the dimensionless drag coefficient appropriate for

_{Dr}*r*. Although the ½ appeared in some early atmospheric versions of (1.3) (e.g., Sutton 1953, p. 232; von Arx 1967, p. 113) to emphasize its derivation from (1.2), modern version are like (1.3), with the ½ absorbed into

*C*.

_{Dr}*C*became a theoretical—not just empirical—coefficient (e.g., Garratt 1992, 52–55):

_{Dr}*k*(=0.40) is the von Kármán constant;

*z*

_{0}, the roughness length; and

*ψ*, an empirical function of the stratification parameter

_{m}*L*, the Obukhov length.

*z*is the height above the surface and

*u*is the friction velocity such that

_{*}*z*

_{0}is the artificial height at which the wind speed is zero and is presumably a fundamental aerodynamic property of the surface (e.g., Wieringa 1993). Hence,

*z*

_{0}and

*C*imply the same information.

_{Dr}*C*more meaningful, we usually eliminate the stability dependence in (1.4) and choose 10 m as the standard reference height. Then (1.4) becomes an expression for the neutral-stability, 10-m drag coefficient:

_{Dr}*z*

_{0}is expressed in meters. Likewise, we often plot

*C*

_{DN}_{10}versus the neutral-stability wind speed at 10 m, which derives from (1.5) with the stability term ignored:

Over 50 years of research to develop a unified parameterization for *C _{DN}*

_{10}has, however, not narrowed the range of reported

*C*

_{DN}_{10}values or satisfactorily explained that range. Reviews repeatedly show plots of widely spread

*C*

_{DN}_{10}values at any given wind speed (e.g., Kraus 1968; Garratt 1977; Blanc 1985; Geernaert 1990; Banner et al. 1999; Toba et al. 2001; Drennan et al. 2005). It has long been suspected that

*C*

_{DN}_{10}responds to other variables than just wind speed, including to sea state; see Jones and Toba (2001) for a review. Ignoring these dependencies was presumed to explain the scatter.

We have some other ideas, however, on why *C _{DN}*

_{10}has been so hard to pin down. First, consider the fundamental uncertainty in

*C*

_{DN}_{10}as computed from (1.9). Measurements of

*u*over the sea have a minimum uncertainty of about ±10%. [See, for instance, Fairall et al. (1996, Table 1) for typical uncertainties.] When

_{*}*u*is small, this uncertainty can approach ±100%. Although

_{*}*U*may be measured at sea on a ship, from a tower or buoy, or from an aircraft with an uncertainty of, say, ±5%, obtaining

_{r}*U*

_{N}_{10}necessitates a stability correction involving

*L*. Because

*L*requires

*C*

_{DN}_{10}(i.e., if we assume that

*U*

_{N}_{10}has the same uncertainty as

*U*) is ±30%. Other phenomena that are presumed to affect

_{r}*C*

_{DN}_{10}—such as sea state, swell, and the relative wind direction between wind and waves—likely make contributions smaller than this in moderate and high winds (cf. Janssen 1997); these effects are thus hidden in the uncertainty of the basic measurements.

Our original datasets come from tabulations in the cited references or were provided by the cited authors. The “Number of runs” gives the number of *u _{*}*–

*U*

_{N}_{10}pairs in the dataset. The cited wind speed range is for

*U*

_{N}_{10}.

In addition to these fundamental considerations of uncertainty, in stable stratification, the magnitude of *ψ _{m}* in (1.4) and (1.5) can become comparable to ln(

*z*/

*z*

_{0}). Hence, in light winds,

*C*can become unrealistically small and is dominated by the uncertainty in

_{Dr}*L*. At the same time,

*U*

_{N}_{10}can also become very small—and, at times, negative [see (2.1) below]. Clearly

*C*

_{DN}_{10}is problematic in such conditions.

*C*

_{DN}_{10}by suggesting that the definition of the drag coefficient [which has its roots in (1.2)] is fundamentally flawed. Namely, we could reasonably infer from (1.9) that

*u*is proportional to

_{*}*U*

_{N}_{10}with a proportionality constant of

*u*versus

_{*}*U*

_{N}_{10}, however, they obtained

*u*and

_{*}*U*

_{N}_{10}are in meters per second,

*a*= 0.051, and

*b*= −0.14 m s

^{−1}. Equation (1.10) shows that

*u*is linearly related to

_{*}*U*

_{N}_{10}but is not proportional to it.

*U*

_{N}_{10}, as in most formulations of

*C*

_{DN}_{10}(e.g., Garratt 1977; Smith 1980; Geernaert 1990; Smith et al. 1992), here

*C*

_{DN}_{10}is more complex. Moreover, because

*b*is negative,

*C*

_{DN}_{10}rises, rolls off, and asymptotes to

*a*

^{2}at high wind speed.

The hurricane community has been searching for behavior like this in *C _{DN}*

_{10}since Emanuel (1995) reported that hurricane models could not produce storms with enough intensity if their drag parameterization was simply an extrapolation of results from moderate wind speeds, which had

*C*

_{DN}_{10}increasing linearly with

*U*

_{N}_{10}, without bounds. Modern hurricane and ocean mixed-layer models, on the other hand, have had some success in predicting storm intensity and ocean response by limiting the value of

*C*

_{DN}_{10}in high winds (Jarosz et al. 2007; Moon et al. 2007; Sanford et al. 2007; Chiang et al. 2011).

Equation (1.10) has features aligned with our own philosophy of air–sea interaction: 1) The experimental coefficient *a* has only half the experimental uncertainty of *C _{DN}*

_{10}and is, thus, more reliably measured; 2) a plot of

*u*versus

_{*}*U*

_{N}_{10}does not have pathological behavior when

*U*

_{N}_{10}is near zero, as do plots of

*C*

_{DN}_{10}; 3) (1.10) minimizes reliance on Monin–Obukhov similarity theory and thereby suffers little from the fictitious correlation typical of these types of analyses (e.g., Mahrt et al. 2003; Klipp and Mahrt 2004; Grachev et al. 2007a,b; Andreas 2011b); and 4) (1.11) produces a natural limit to

*C*

_{DN}_{10}.

*u*versus

_{*}*U*

_{N}_{10}for data that we had on hand. Figure 1 shows that our results corroborate those of Foreman and Emeis. We find

As such, we add over 6000 more values measured by low-flying aircraft in winds up to 27 m s^{−1} to the 778 points shown in Fig. 1. This aircraft set also shows a straight-line relation between *u _{*}* and

*U*

_{N}_{10}in the aerodynamically rough regime, and the fitting coefficients are not statistically different from those in (1.12).

Both datasets also suggest that *u _{*}* follows the prediction for aerodynamically smooth flow for low

*U*

_{N}_{10}. Consequently, we devise a continuous drag relation for all

*U*

_{N}_{10}by smoothly combining this aerodynamically smooth regime with (1.10) for the aerodynamically rough regime.

On extrapolating this relation to hurricane-strength winds, we find that it predicts the rolloff in *C _{DN}*

_{10}that hurricane models seem to require. Moreover, the straight-line behavior of

*u*with

_{*}*U*

_{N}_{10}—even in high winds—and the rolloff in

*C*

_{DN}_{10}are compatible with theoretical models by Moon et al. (2007) and Mueller and Veron (2009) that compute the air–sea drag as resulting from just skin friction and the form drag from flow separation over the waves. In other words, our analysis suggests that there is no need to invoke exotic processes, such as sea spray loading or the disintegration of the air–sea interface, to explain the roll off in

*C*

_{DN}_{10}with increasing wind speed. Wind–wave coupling suffices.

## 2. Datasets

Table 1 summarizes the data with which we made our first test of the Foreman and Emeis (2010) approach (i.e., our Fig. 1). We will refer to this as the “original” dataset. Most of these sets are available as tabulations in the cited references. We obtained the Fronts and Atlantic Storm Track Experiment (FASTEX) and Greenland Flow Distortion Experiment (GFDex) data, however, as electronic files from the scientists referenced for these sets.

We did not include the Coupled Boundary Layers Air–Sea Transfer (CBLAST) hurricane dataset mentioned in Table 1 in this analysis because these data are not consistent with our other data: The *u _{*}* values tend to be low, as we will show later. We suspect that this bias resulted because these aircraft data were obtained at flight levels that were never below 70 m, were as high as 383 m, and had a median level of 193 m, while the depth of the boundary layer for these flights during Hurricanes Fabian and Isabel was 350–550 m (Zhang et al. 2009). That is, because the stress is known to decrease with height through the boundary layer (e.g., Caughey et al. 1979; Nicholls and Readings 1979; Zhang et al. 2009; Wyngaard 2010, 244–247, 286–287), the measured flight-level stress was less than the surface stress. Although French et al. (2007) tried to correct for this flux divergence, their reported values of

*u*remain low.

_{*}Table 2 summarizes a second set of data that we use in this study. Because all of these data come from low-flying aircraft, we will refer to this as our “aircraft” dataset.

Our “aircraft” dataset consists of 4-km flight segments that we processed ourselves; see Mahrt et al. (2012) for additional details. “Number of runs” here is the number of such 4-km legs. The “Altitude range” gives the aircraft flight level; the wind speed noted is the range of measured wind speeds at those levels. The cited references give more details on the measurements.

Four different aircraft collected these data: the National Oceanic and Atmospheric Administration’s Long-EZ, the C-130 and Electra from the National Center for Atmospheric Research (NCAR), and the Twin Otter from the Naval Postgraduate School’s Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) (Khelif et al. 2005).

To measure the turbulent wind vector required for computing the momentum flux, the Twin Otter, C-130, and Electra used five-port radomes on the nose of the aircraft. Each port had a pressure sensor that sampled at 20–25 Hz; Lenschow (1986) describes the principles of obtaining the wind vector from aircraft pressure measurements. Each of these three aircraft used the Global Positioning System (GPS) to correct the aircraft’s inertial navigation system to find true ground speed (Khelif et al. 1999).

The Long-EZ used the Best Atmospheric Turbulence Probe (BAT) for measuring the wind vector (Crawford and Dobosy 1992; Garman et al. 2006). This is a baseball-bat-shaped device with nine pressure ports in its thicker end; it protruded from the Long-EZ into the undisturbed free stream ahead of the aircraft. A pressure sensor in each port sampled at 50 Hz; again, the aircraft’s inertial navigation system was corrected with GPS positioning to obtain the true wind vector with respect to the ground.

Regardless of aircraft, each flux value in the aircraft set is the average over a 4-km flight segment that we processed ourselves from the raw data. We use the flight-level momentum flux as

We did only minimal initial screening of the aircraft data for quality control (cf. Mahrt et al. 2012). In low winds, wave effects or uncertainty in the aircraft turbulence measurements can produce a stress that appears upward, contrary to boundary layer theory. We screened for such spurious stress measurements and eliminated 662 cases from the initial 6080 flight legs summarized in Table 2. Over 95% of these questionable measurements occurred with flight-level winds of less than 8 m s^{−1}.

On the other hand, the original authors of the Janssen, Risø Air–Sea Experiment (RASEX), Southern Ocean Waves Experiment (SOWEX), and Surface Wave Dynamics Experiment (SWADE) data in Table 1 probably screened these datasets more strictly before publishing them and reported only “high quality” fluxes. Although the FASTEX, GFDex, and Humidity Exchange over the Sea (HEXOS) sets had been processed when we received them, we suspect that these sets had been screened only for instrument malfunctions—not for stationarity, homogeneity, wave characteristics, or other relevant quality metrics.

We further screened the SWADE set ourselves. The original set has 126 records; and Donelan et al. (1997) had identified whether each record represented conditions of wind sea or whether there was swell with wind in the same direction, in the opposing direction, or at right angles. In Fig. 1, we include only the 27 SWADE cases with wind sea; later we will discuss the remaining SWADE data that had swell present.

A distinct feature of all the data represented in Tables 1 and 2 is that they come from eddy-covariance measurements of the momentum and heat fluxes. We eschewed available datasets that were based on inertial-dissipation estimates of the fluxes because such fluxes rely heavily on Monin–Obukhov similarity theory. We want to minimize our reliance on similarity theory.

Initially in our analysis, we will treat the original and aircraft datasets separately. In effect, we are using the aircraft dataset to validate our analysis of the original dataset, or vice versa.

The method used for estimating *U _{N}*

_{10}for this work is crucial. One approach, obviously, is to use (1.8) and the measured

*u*. The required

_{*}*z*

_{0}could come from the corresponding measured value, from a lookup table or a comparable standard (e.g., Panofsky and Dutton 1984, 121–123; Stull 1988, p. 380; Wieringa 1993), or from a parameterization such as the Charnock relation. With that approach, however, the dependent variable in the analysis, the measured

*u*, will be very well (and artificially) correlated with the independent variable

_{*}*U*

_{N}_{10}.

*U*

_{N}_{10},

*U*

_{N}_{10}is most sensitive to the actual wind measurement,

*U*(

*z*), and has generally only a modest built-in dependence on the measured

*u*. Furthermore, if the measurement height is close to 10 m [when ln(

_{*}*z*/10) ~ 0] and if the stratification is near neutral [when

*ψ*(

_{m}*z*/

*L*) ~ 0], the

*U*

_{N}_{10}obtained from (2.1) has very weak built-in correlation with

*u*, and

_{*}*U*

_{N}_{10}has nearly the same uncertainty as

*U*(

*z*). A

*U*

_{N}_{10}estimated from (1.8), in contrast, always has an uncertainty no smaller than the uncertainty in

*u*.

_{*}When we could, we estimated *U _{N}*

_{10}from (2.1). For all the aircraft data in Table 2, this was the case. For

*ψ*we used the function from Paulson (1970) in unstable stratification and the function from Grachev et al. (2007a) in stable stratification. For the FASTEX, GFDex, and RASEX data in Table 1, we had enough information to also calculate

_{m}*U*

_{N}_{10}according to (2.1). In the CBLAST-hurricane set,

*U*

_{N}_{10}was estimated with a stepped-frequency microwave radiometer (Drennan et al. 2007) and thus has no built-in dependence on

*u*. In the HEXOS set, DeCosmo (1991) reported only

_{*}*U*

_{N}_{10}and did not explain how she obtained this value. In the SOWEX set, Banner et al. (1999) obtained

*U*

_{N}_{10}from (1.8). In the SWADE set, Donelan et al. (1997) reported

*U*

_{N}_{10}and explained that they obtained it from (2.1).

Finally, Janssen (1997) reported only a variable denoted *U*_{10} but did not explain how this was obtained or whether it is the neutral-stability value. He did, however, report two simultaneous, independent sets of measurements: *u _{*}* and

*U*

_{10}were measured with both a pressure anemometer and a sonic anemometer. Under the assumption that

*U*

_{10}is

*U*

_{N}_{10}but to avoid the built-in correlation in case Janssen estimated

*U*

_{N}_{10}from (1.8), we switched the pressure anemometer and sonic measurements of

*u*in our analysis. In other words, in Fig. 1 and subsequent figures, we plot

_{*}*u*from the sonic against the corresponding

_{*}*U*

_{N}_{10}value from the pressure anemometer and

*u*from the pressure anemometer against the sonic

_{*}*U*

_{N}_{10}. Now, in the Janssen set, the

*u*and

_{*}*U*

_{N}_{10}pairs have no direct built-in correlation from a shared

*u*.

_{*}## 3. Results

Figure 2 is a plot like Fig. 1 but for the aircraft data summarized in Table 2. Unlike the original set, a high percentage of the aircraft data were collected in stable stratification—1123 of the 5418 data records left after our screening the stress. While we are not concerned about the flux divergence in unstable stratification for fluxes measured at the aircraft altitudes in this set (up to 49 m), we worry about possible flux divergence in stable stratification because of the generally shallower boundary layer. (Remember, we obtain the *u _{*}* surface value from the uncorrected momentum flux measured at flight level.)

To avoid biasing our analysis with *u _{*}* values biased low because of vertical flux divergence, we made several plots and analyses as in Fig. 2. Table 3 summarizes the calculations. First, we considered the aircraft data collected in unstable stratification and, separately, the data collected in stable stratification. Admittedly, the stable cases constituted only 26% as much data as in the unstable cases; still, as expected because of the flux divergence, the

*a*value for the data collected in stable stratification is significantly less than the

*a*value for the unstable cases.

Fits to the model *r*, the parameter *a*, the 95% confidence interval on *a* [from Bendat and Piersol (1971), p. 131], parameter *b*, and the 95% confidence interval on *b*.

Moreover, when we further segregated the stable data into cases with *z* is the aircraft altitude and *L* is the measured Obukhov length, the set with the more stable conditions had a smaller *a* value than for the weakly stable set (Table 3). Finally, when we included just this weakly stable set (i.e., *a* and *b* values were indistinguishable from the *a* and *b* values for just the unstable data. Hence, as our analyzed “aircraft” dataset here, we use all the aircraft data collected in unstable stratification and the data from weakly stable stratification when

In both Figs. 1 and 2, the data clouds change character in the *U _{N}*

_{10}range 8–10 m s

^{−1}. Below this range, the points have a shallower slope than above it. This tendency is compatible with aerodynamically smooth flow at low wind speeds and the transition to aerodynamically rough flow as the wind speed increases.

*ν*is the kinematic viscosity of air. For

*α*is the Charnock parameter,

Figures 3 and 4 show our evaluation of *R _{*}* for the two datasets. Each figure includes the individual values, bin averages in 1 m s

^{−1}bins in

*U*

_{N}_{10}, and bin medians. Because

*R*values are approximately lognormally distributed, the proper bin average is computed as the geometric mean—as opposed to the arithmetic mean. If the

_{*}*R*values were perfectly lognormal within a bin, the median would be the same as the geometric mean. Notice that in both figures the medians and geometric means are close.

_{*}*U*

_{N}_{10}is greater than 8 m s

^{−1}. As a conservative estimate, we therefore used only the data for which

*u*and

_{*}*U*

_{N}_{10}are in meters per second. The correlation coefficient is 0.835.

Bendat and Piersol (1971, p. 131) explain that both the slope and intercept in (1.12) and (3.3) follow Student *t* distributions. We therefore calculated 95% confidence intervals on the slopes and intercepts fitted to the original and aircraft data in Figs. 1 and 2 (Table 3). For the slope in (1.12), the 95% confidence interval is [0.0564, 0.0599]; for the intercept, [−0.239, −0.189]. For the slope in (3.3), the 95% confidence interval is [0.0566, 0.0600]; for the intercept, [−0.266, −0.221].

Because both slope and intercept intervals coincide well, Figs. 1 and 2 give us essentially the same result. In effect, we validate the fitting line in Fig. 1 with Fig. 2, and vice versa. Henceforth, we will use the coefficients in (3.3) as our main result because they come from the larger dataset.

*R*values at small

_{*}*U*

_{N}_{10}in Figs. 3 and 4 are related to the “pathological behavior” in

*C*

_{DN}_{10}for small

*U*

_{N}_{10}that we mentioned earlier. The

*z*

_{0}values used to create these figures came from the left-hand side of (2.1) [or, alternatively, from (1.5)]:

*U*

_{N}_{10}and

*u*are small, their uncertainties often cause

_{*}*z*

_{0}is unrealistically large, and so is

*R*.

_{*}The slope in (1.10) that Foreman and Emeis (2010) reported (0.051) is smaller than our values, and their intercept (−0.14 m s^{−1}) is larger. We suspect that, because they used *U _{N}*

_{10}= 8 m s

^{−1}as the lower limit for aerodynamically rough flow in their analysis, they may have retained some data reflecting aerodynamic transition. Notice in Figs. 1 and 2 how the

*u*values at lower winds turn up and away from (1.12) and (3.3). By including such data in their calculations, Foreman and Emeis would have inadvertently decreased

_{*}*a*and increased

*b*from what the data in truly rough flow suggest.

At *U _{N}*

_{10}= 9 m s

^{−1}, (3.3) gives

*u*= 0.28 m s

_{*}^{−1}. In our analysis, this is the friction velocity at the transition to aerodynamically rough flow. For comparison, Wu (1969) concluded that this transition is at

*U*

_{N}_{10}= 7 m s

^{−1}, whereas Wu (1980) obtained

*u*= 0.263 m s

_{*}^{−1}, although he assumed

*R*= 2.33 at the transition to aerodynamically rough flow. From his data analysis, Melville (1977) concluded that

_{*}*u*was in the range 0.15–0.30 m s

_{*}^{−1}at the onset of aerodynamically rough flow, although he also used for the transition an

*R*limit (=2) lower than ours. On invoking the Charnock relation, Kraus and Businger (1994, p. 145) and Foreman and Emeis (2010) estimated

_{*}*u*was 0.29 and 0.28 m s

_{*}^{−1}, respectively, at the transition to aerodynamically rough flow.

Hence, our estimate that *u _{*}* = 0.28 m s

^{−1}at the transition to aerodynamically rough flow agrees with most previous estimates; but our result that

*U*

_{N}_{10}= 9 m s

^{−1}at this transition is a bit higher than previous estimates.

## 4. Discussion

### a. Consistency of the results

Thoughtful readers might suspect that the data clouds in Figs. 1 and 2 obscure differences in behavior among the different datasets that are typical in plots of *C _{DN}*

_{10}versus

*U*

_{N}_{10}. Then our fitting lines in Figs. 1 and 2 would just be average results that ignore true differences in drag relations among the sets. To allay these worries, we created Figs. 5 and 6.

These show the individual datasets that went into Figs. 1 and 2. Figure 5 shows our original data; Fig. 6, the aircraft data. Reassuringly, 17 of the 18 datasets individually either lie along our aircraft fit, (3.3), suggest aerodynamically smooth scaling at low wind speed, or do both. That is, the individual datasets are not biased high or low such that, when we fitted (1.10) to the two consolidated datasets, the fitting line simply split the difference between systematically high and systematically low values.

The one exception to this behavior is the CBLAST hurricane dataset (Fig. 5). These *u _{*}* values seem to be too low—probably for the reasons we discussed earlier. We have therefore not included them in our least squares fittings.

### b. A unified roughness parameterization

*U*

_{N}_{10}= 0 to where they intersect our aerodynamically rough results, (1.12) and (3.3), we fitted them with

*U*

_{N}_{10}in [0, 8.76 m s

^{−1}].

*u*with a smooth, differentiable function of

_{*}*U*

_{N}_{10}by combining (3.3) and (4.2) in a hyperbola. The result that best fits our data for all

*U*

_{N}_{10}is

*u*and

_{*}*U*

_{N}_{10}are both in meters per second.

Figure 7 shows how well this expression fits the bin-averaged *u _{*}* values from the combined original and aircraft datasets. Only for

*U*

_{N}_{10}< 3 m s

^{−1}do the data in Fig. 7 deviate significantly from (4.3). Instead of highlighting missing physics, these three large

*u*values reveal how difficult measuring

_{*}*u*is in low winds.

_{*}Equation (4.3) now becomes a core component of a new bulk flux algorithm that we are developing. The appendix sketches how we use (4.3) in this algorithm.

For readers used to looking at flux algorithms in terms of *C _{DN}*

_{10}, we can insert (4.3) into (1.9) to obtain an expression for

*C*

_{DN}_{10}for all wind speeds. Figure 8 shows this result and how it fits the bin-averaged

*C*

_{DN}_{10}values in our combined original and aircraft datasets.

Figures 7 and 8 also reiterate some of the advantages of a drag relation based on *u _{*}* over one based on

*C*

_{DN}_{10}that we described in the introduction. Although the averaged

*u*values in the three lowest

_{*}*U*

_{N}_{10}bins in Fig. 7 do not follow aerodynamically smooth scaling, at least they are well behaved and have some of the smallest error bars in the plot. The

*C*

_{DN}_{10}values in the two lowest

*U*

_{N}_{10}bins in Fig. 8, in contrast, are above the upper limit of the plot and thus do not show up at all. Moreover, the errors bars on the

*C*

_{DN}_{10}values for small

*U*

_{N}_{10}are generally the largest on the plot and even encompass negative

*C*

_{DN}_{10}values.

Figure 8 also suggests that the distribution of individual *C _{DN}*

_{10}values within bins is skewed toward larger values: For

*U*

_{N}_{10}< 9 m s

^{−1}, the averages are larger than the medians. In Fig. 7, (4.3) fits the bin-averaged

*u*values very well. In contrast, the bin-averaged

_{*}*C*

_{DN}_{10}values for

*U*

_{N}_{10}< 9 m s

^{−1}in Fig. 8 are above the

*C*

_{DN}_{10}curve derived from (4.3) although the same data as in Fig. 7 went into this plot. All of these features are evidence of what we termed pathological behavior in

*C*

_{DN}_{10}.

### c. Drag relations in high winds

Figure 9 shows (4.3) extrapolated to hurricane-strength winds. The figure also shows both the original and aircraft datasets to emphasize how consistent they are. Furthermore, Fig. 9 includes the CBLAST hurricane data to demonstrate that they are generally below reliable data measured at similar wind speeds.

The main features of Fig. 9, however, are the curves attributed to Moon et al. (2007) and Mueller and Veron (2009). These are theoretical results in which both sets of authors modeled the total wind stress on the sea as a combination of the viscous stress (or skin friction or tangential stress), a wave-induced stress from form drag, and the reduction of the viscous and wave-induced stresses by sheltering (or flow separation).

Moon et al. (2007) obtained their results [summarized in their (4) and (5)] by simulating the surface stress in 10 Atlantic hurricanes with wind speeds up to 70 m s^{−1} using the hurricane model of Moon et al. (2004). They then inferred *z*_{0} from this modeled stress through similarity theory. Their *z*_{0} values are thus statistical averages in wind speed bins. The Mueller and Veron (2009) model, on the other hand, simply provides a deterministic prediction of the surface stress for the given input conditions. [F. Veron (2011, personal communication) reran the Mueller and Veron model for our specified conditions for wind speeds up to 80 m s^{−1}.] Neither group evidently realized that its model yielded a nearly straight-line relation between *u _{*}* and

*U*

_{N}_{10}for

*U*

_{N}_{10}above 20–25 m s

^{−1}(Fig. 9). We use this theoretical behavior in these two models to justify extrapolating our own result, (4.3), to hurricane-strength winds.

We cast our results in the familiar form of a drag coefficient in Fig. 10. Remember, because of (1.11), our drag coefficient rolls off and approaches an asymptotic limit of *a*^{2} in high winds: *α* = 0.0185, Andreas et al. (2008)], to get a unified expression for *z*_{0} (e.g., Zilitinkevich 1969; Smith 1988; Fairall et al. 1996).

In models of the oceanic mixed layer under hurricanes, both Sanford et al. (2007) and Chiang et al. (2011) based parameterizations for *C _{DN}*

_{10}on the observations reported by Powell et al. (2003). While we do not endorse the Powell et al. results for several reasons, we show in Fig. 10 the Sanford et al. and Chiang et al.

*C*

_{DN}_{10}parameterizations because they are continuous functions, like the other curves in Fig. 10, and because, we feel, they represent the lowest reasonable drag coefficient possible in high winds.

Finally, Fig. 10 also includes the drag coefficients that Bell (2010) obtained by using dropsondes launched in Hurricanes Fabian and Isabel in 2003 to estimate the angular momentum budget under the assumption that the storms were axisymmetric. These are the most recent determinations of air–sea drag in high winds that are available. Bell computed drag coefficients for six separate aircraft missions and used 72 control volumes per mission for flux calculations. The error bars on his points in Fig. 10 are thus (probably) standard deviations around the means of the 72 samples per mission. That is, they are not uncertainty estimates but rather indicators of the scatter in the individual values.

Bell’s (2010) results, unfortunately, do not help us decide which of the candidate drag parameterizations in Fig. 10 is the most realistic. His data range from below the lowest realistic parameterization to above our parameterization, which we suggest gives the greatest upper bound.

Ingel (2011) refers to the rolloff in *C _{DN}*

_{10}depicted in Fig. 10 as a “drag crisis,” invoking the terminology of classical fluid mechanics when laminar flow transitions to turbulence and the drag coefficients of cylinders and spheres drop suddenly by a factor of 5 with increasing Reynolds number (Monin and Yaglom 1971, 82–83; Faber 1995, 266–267). That is, according to Ingel, something fundamental about the air–sea coupling has changed. Kudryavtsev (2006) believes that this rolloff signals a saturation of the surface stress: the stress no longer increases with increasing wind speed. Neither of these inferences is true.

If *u _{*}* is a

*linear*function of

*U*

_{N}_{10}for moderate and high wind speeds—as our data and the theories of Moon et al. (2007) and Mueller and Veron (2009) suggest—nothing fundamental changes about the way air and sea couple as the wind speed increases. There is certainly no drag crisis in the classical sense. Nor does the surface stress saturate: (3.3) confirms that

*u*and thus the surface stress increase with wind speed for all wind speeds.

_{*}The roughness length does saturate, however, as Donelan et al. (2004) suggest. From (3.3) and (3.4), we see that *z*_{0} approaches a limiting value of *k*/*a*)] at large *U _{N}*

_{10}.

Because Moon et al. (2007) and Mueller and Veron (2009) account for the behavior that we see in the data by modeling just wind–wave coupling, invoking more exotic processes to explain the rolloff in *C _{DN}*

_{10}seems unnecessary. Many of these attempts to explain the rolloff in

*C*

_{DN}_{10}involve injecting sea spray into the near-surface air in copious amounts (e.g., Makin 2005; Kudryavtsev 2006; Soloviev and Lukas 2010; Ingel 2011; Bianco et al. 2011). Such spray loading may stabilize the near-surface air and thus reduce the momentum transfer somewhat in very high winds.

Shpund et al. (2011) recently suggested, however, that such spray loading may not be as important as has been estimated from one-dimensional, eddy-diffusivity models (e.g., Lighthill 1999; Makin 2005; Kudryavtsev 2006; Ingel 2011; Bianco et al. 2011). When Shpund et al. introduced large eddies into their two-dimensional Lagrangian model, these eddies carried spray that was generated at the sea surface to higher levels in the marine boundary layer, thereby reducing the spray loading near the surface.

Still, we cannot rule out the possibility that, with increasing wind speed, spray loading may cause the actual drag relation to fall slightly below our result in Fig. 10. Nevertheless, wind–wave coupling appears to be the dominant mechanism causing the drag coefficient to roll off.

### d. Effects of sea state

The idea that the roughness of the sea surface *z*_{0} may depend on wave age *C _{p}*/

*u*, where

_{*}*C*is the phase speed of the dominant waves, goes back at least to Kitaigorodskii and Volkov (1965). Many authors have investigated whether this and other sea-state parameters influence air–sea drag (e.g., Smith et al. 1992; Janssen 1997; Bourassa et al. 2001; Oost et al. 2002; Drennan et al. 2005). Jones and Toba (2001) dedicated an entire book to this subject of how to parameterize the drag; in that book, Toba et al. (2001) highlighted the fact that “there is still lively debate” on this issue of how sea state influences the air–sea drag.

_{p}Briefly, the question is whether an air–sea drag parameterization can be successful as a function of just wind speed or whether it must also include sea-state variables such as wave age, wave height, wave period, or swell height and direction. Although our drag parameterization for aerodynamically rough flow requires only *U _{N}*

_{10}, we cannot conclude that sea state is unimportant. After all, our extrapolated drag relation for aerodynamically rough flow, (4.3), corroborates the model predictions of Moon et al. (2007) and Mueller and Veron (2009) (Figs. 9 and 10), which involved integrating the ocean wave spectrum. Because these models depend implicitly on wave age, wave height, and wave steepness through the wave spectrum, our results, in turn, imply that these parameters are important for air–sea coupling.

The SWADE dataset (Table 1), which includes observations of swell, also has implications for our parameterization. Remember that, in the SWADE set, Donelan et al. (1997) noted whether their measurements of surface stress were made over a wind sea or in the presence of swell that was across wind, counter to the wind, or following the wind. Figure 11 shows the full SWADE dataset, plotted as *u _{*}* versus

*U*

_{N}_{10}, with these four regimes identified.

Although we add only 99 points to our analysis, Fig. 11 suggests how robust our *u _{*}*-versus-

*U*

_{N}_{10}analysis is and how weakly swell affects the drag. For most of the

*U*

_{N}_{10}range, the swell cases on average have

*u*values only 10%–15% higher than the wind sea cases. At the upper end of the

_{*}*U*

_{N}_{10}range, however, the swell and wind sea cases are essentially coincident (cf. Drennan et al. 2003). Only for the lowest

*U*

_{N}_{10}values—well below the aerodynamically rough regime (

*U*

_{N}_{10}< 9 m s

^{−1})—do swell cases deviate significantly from the wind sea cases and from the central tendency of our analysis.

We hypothesize that, in moderate and high winds, the sea surface is so strongly forced by the wind that swell has a low-order effect on the surface stress. In light winds, on the other hand, the relationship between wind stress and swell becomes significant.

Still, the observed erratic behavior of the drag in light winds—expressed as the drag coefficient *C _{DN}*

_{10}, the roughness length

*z*

_{0}, or

*u*(Figs. 7, 8, and 11; e.g., Rutgersson et al. 2001; Drennan et al. 2003)—may simply reflect how difficult defining a meaningful reference frame and then measuring stress is in light winds with swell present. In fact, with swell opposing a light wind, the ocean may actually drag the atmosphere (e.g., Smedman et al. 1994; Grachev and Fairall 2001; Grachev et al. 2003).

_{*}## 5. Conclusions

Despite many measurements, the drag coefficient formulated as *C _{DN}*

_{10}is naturally prone to such variability. Here, we therefore evaluated an alternative formulation of the air–sea drag relation, following the suggestion by Foreman and Emeis (2010).

*u*, measured over the ocean in aerodynamically rough flow increases linearly with

_{*}*U*

_{N}_{10}, the 10-m, neutral-stability wind speed. We find

*C*

_{DN}_{10}.

We suggest that wind–wave coupling explains (5.1). Theoretical models for the surface stress by Moon et al. (2007) and Mueller and Veron (2009) include terms for only the skin friction, form drag on the waves, and flow sheltering. Yet, these produce nearly straight-line relations between *u _{*}* and

*U*

_{N}_{10}up to winds of major hurricane strength. (Neither group evidently realized this behavior.) Furthermore, both model predictions are very close to our (5.1). As a result, we conclude that known processes involving wind–wave coupling may be enough to explain the behavior of the air–sea drag for all wind speeds. These theoretical results also motivate our extrapolating (5.1) to hurricane-strength winds.

The literature contains data-based estimates that suggest *C _{DN}*

_{10}can be as low as

^{−1}winds (e.g., Powell et al. 2003). We believe that this estimate is the smallest lower bound on the drag coefficient in hurricane-strength winds. On the other hand, one way to view our result (5.2) is as the greatest upper bound on the drag coefficient.

Processes that the models of Moon et al. (2007) and Mueller and Veron (2009) did not include—such as spray loading—may reduce the drag coefficient from what (5.2) predicts. We hypothesize, however, that any such effects will be second order, reducing *C _{DN}*

_{10}from the level in (5.2) by, perhaps, 10%. Because the sea surface is so strongly forced in high winds, we also hypothesize that swell will negligibly affect air–sea drag for wind speeds above ~15 m s

^{−1}.

From the behavior of the roughness Reynolds number in our two datasets, we also estimated the wind speed and the friction velocity at which the sea surface becomes aerodynamically rough. Although the roughness Reynolds numbers are scattered, we have enough data to reliably determine mean behavior. We conclude that the sea surface becomes aerodynamically rough for *U _{N}*

_{10}between 8 and 10 m s

^{−1}; as an operational estimate, we use

*U*

_{N}_{10}= 9 m s

^{−1}as the wind speed at transition. From (5.1), this wind speed gives

*u*= 0.28 m s

_{*}^{−1}as the friction velocity when the sea surface becomes aerodynamically rough.

Although flux measurements at sea in light winds have larger uncertainty, our data suggest that *u _{*}* follows aerodynamically smooth scaling at low

*U*

_{N}_{10}, where the roughness length is presumed to obey

*U*

_{N}_{10}< 9 m s

^{−1}, this expression produces nearly straight-line behavior in plots of

*u*versus

_{*}*U*

_{N}_{10}. We thus fitted this smooth region with a straight line, (4.2), to complement our straight-line result in rough flow, (5.1).

*u*and

_{*}*U*

_{N10}are in meters per second.

## Acknowledgments

We thank the following colleagues for providing datasets: Will Drennan and Jun Zhang for the CBLAST-hurricane set; Ola Persson, Jeff Hare, Chris Fairall, and Bill Otto for the FASTEX set; Nína Petersen and Ian Renfrew for the GFDex set; and Janice DeCosmo for help with the HEXOS set. We thank Fabrice Veron for rerunning the Mueller and Veron (2009) model to produce the curves we show in Figs. 9 and 10 and Emily Moynihan for preparing Figs. 5 and 6. Two anonymous reviewers made insightful comments that helped us focus the manuscript better. The U.S. Office of Naval Research supported this work with Award N00014-11-1-0073 to NorthWest Research Associates. ONR also supported ELA through the National Ocean Partnership Program with Award N00014-10-1-0154 to the University of Rhode Island, for whom he is a subcontractor.

## APPENDIX

### Bulk Flux Algorithm

*u*parameterization (4.3) was to use it in a bulk flux algorithm for high-wind, spray conditions. Sea spray enhances the air–sea exchange of sensible and latent heat; consequently, in such an algorithm, the total sensible (

_{*}*H*

_{s,T}) and latent (

*H*

_{L,T}) heat fluxes each include two terms, an

*interfacial*part (

*H*and

_{s}*H*) and a

_{L}*spray-mediated*part (

*Q*

_{s,}

_{sp}and

*Q*

_{L,sp}) (Andreas et al. 2008; Andreas 2011a). With the surface stress or momentum flux,

*τ*, the bulk flux algorithm thus comprises three main equations:

*τ*once

*U*

_{N}_{10}is calculated from (2.1). Andreas et al. (2008) give the equations for computing

*Q*

_{s,sp}and

*Q*

_{L,sp}, which also depend on

*u*. Our specifying

_{*}*u*directly from

_{*}*U*

_{N}_{10}, however, allows us to deviate from standard methods for calculating the interfacial fluxes

*H*and

_{s}*H*.

_{L}*c*is the specific heat of air at constant pressure,

_{p}*L*is the latent heat of vaporization, and

_{υ}*θ*and

_{*}*q*are temperature and specific humidity flux scales analogous to

_{*}*u*. These

_{*}*θ*and

_{*}*q*, in turn, appear in the Monin–Obukhov similarity expressions for the profiles of average potential temperature Θ and specific humidity

_{*}*Q*in the atmospheric surface layer (Garratt 1992, 52–54; Andreas et al. 2008):

*r*is a reference height; Θ

_{r}and

*Q*are thus the average potential temperature and specific humidity at height

_{r}*r*; Θ

_{s}and

*Q*are the corresponding values at the sea surface. Also,

_{s}*L*is the Obukhov length—a function of

*u*,

_{*}*θ*, and

_{*}*q*. For the stability correction

_{*}*ψ*, we use the same sources that we did for the

_{h}*ψ*function in (2.1): Paulson (1970) for unstable stratification, and Grachev et al. (2007a) for stable stratification. Andreas et al. (2008) explain how we estimate the roughness lengths

_{m}*z*and

_{T}*z*.

_{Q}*U*

_{N}_{10}, (2.1), are coupled through the Obukhov length.

Notice that this new algorithm requires no parameterization for *z*_{0}, which is one of the most variable and uncertain quantities in surface layer similarity theory.

## REFERENCES

Anderson, K., and Coauthors, 2004: The RED Experiment: An assessment of boundary layer effects in a trade winds regime on microwave and infrared propagation over the sea.

,*Bull. Amer. Meteor. Soc.***85**, 1355–1365.Andreas, E. L, 2011a: Fallacies of the enthalpy transfer coefficient over the ocean in high winds.

,*J. Atmos. Sci.***68**, 1435–1445.Andreas, E. L, 2011b: The fallacy of drifting snow.

,*Bound.-Layer Meteor.***141**, 333–347.Andreas, E. L, and G. Treviño, 2000: Comments on “A physical interpretation of von Kármán’s constant based on asymptotic considerations—A new value.”

,*J. Atmos. Sci.***57**, 1189–1192.Andreas, E. L, P. O. G. Persson, and J. E. Hare, 2008: A bulk turbulent air–sea flux algorithm for high-wind, spray conditions.

,*J. Phys. Oceanogr.***38**, 1581–1596.Andreas, E. L, P. O. G. Persson, R. E. Jordan, T. W. Horst, P. S. Guest, A. A. Grachev, and C. W. Fairall, 2010: Parameterizing turbulent exchange over sea ice in winter.

,*J. Hydrometeor.***11**, 87–104.Banner, M. L., W. Chen, E. J. Walsh, J. B. Jensen, S. Lee, and C. Fandry, 1999: The Southern Ocean Waves Experiment. Part I: Overview and mean results.

,*J. Phys. Oceanogr.***29**, 2130–2145.Batchelor, G. K., 1970:

*An Introduction to Fluid Dynamics.*Cambridge University Press, 615 pp.Bell, M. M., 2010: Air–sea enthalpy and momentum exchange at major hurricane wind speeds. Ph.D. dissertation, Naval Postgraduate School, 133 pp.

Bendat, J. S., and A. G. Piersol, 1971:

*Random Data: Analysis and Measurement Procedures.*Wiley-Interscience, 407 pp.Bianco, L., J.-W. Bao, C. W. Fairall, and S. A. Michelson, 2011: Impact of sea-spray on the atmospheric surface layer.

,*Bound.-Layer Meteor.***140**, 361–381.Blanc, T. V., 1985: Variation of bulk-derived surface flux, stability, and roughness results due to the use of different transfer coefficient schemes.

,*J. Phys. Oceanogr.***15**, 650–669.Bourassa, M. A., D. G. Vincent, and W. L. Wood, 2001: A sea state parameterization with nonarbitrary wave age applicable to low and moderate wind speeds.

,*J. Phys. Oceanogr.***31**, 2840–2851.Caughey, S. J., J. C. Wyngaard, and J. C. Kaimal, 1979: Turbulence in the evolving stable boundary layer.

,*J. Atmos. Sci.***36**, 1041–1052.Chiang, T.-L., C.-R. Wu, and L.-Y. Oey, 2011: Typhoon Kai-Tak: An ocean’s perfect storm.

,*J. Phys. Oceanogr.***41**, 221–233.Crawford, T. L., and R. J. Dobosy, 1992: A sensitive fast-response probe to measure turbulence and heat flux from any airplane.

,*Bound.-Layer Meteor.***59**, 257–278.DeCosmo, J., 1991: Air–sea exchange of momentum, heat and water vapor over whitecap sea states. Ph.D. dissertation, University of Washington, 212 pp.

Donelan, M. A., W. M. Drennan, and K. B. Katsaros, 1997: The air–sea momentum flux in conditions of wind sea and swell.

,*J. Phys. Oceanogr.***27**, 2087–2099.Donelan, M. A., B. K. Haus, N. Reul, W. J. Plant, M. Stiassnie, H. C. Graber, O. B. Brown, and E. S. Saltzman, 2004: On the limiting aerodynamic roughness of the ocean in very strong winds.

,*Geophys. Res. Lett.***31**, L18306, doi:10.1029/2004GL019460.Drennan, W. M., H. C. Graber, D. Hauser, and C. Quentin, 2003: On the wave age dependence of wind stress over pure wind seas.

,*J. Geophys. Res.***108**, 8062, doi:10.1029/2000JC000715.Drennan, W. M., P. K. Taylor, and M. J. Yelland, 2005: Parameterizing the sea surface roughness.

,*J. Phys. Oceanogr.***35**, 835–848.Drennan, W. M., J. A. Zhang, J. R. French, C. McCormick, and P. G. Black, 2007: Turbulent fluxes in the hurricane boundary layer. Part II: Latent heat flux.

,*J. Atmos. Sci.***64**, 1103–1115.Edson, J., and Coauthors, 2007: The Coupled Boundary Layers and Air–Sea Transfer Experiment in low winds.

,*Bull. Amer. Meteor. Soc.***88**, 341–356.Emanuel, K. A., 1995: Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics.

,*J. Atmos. Sci.***52**, 3969–3976.Faber, T. E., 1995:

*Fluid Dynamics for Physicists.*Cambridge University Press, 440 pp.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**, 3747–3764.Foreman, R. J., and S. Emeis, 2010: Revisiting the definition of the drag coefficient in the marine atmospheric boundary layer.

,*J. Phys. Oceanogr.***40**, 2325–2332.Francis, J. R. D., 1954: Wind stress on a water surface.

,*Quart. J. Roy. Meteor. Soc.***80**, 438–443.French, J. R., W. M. Drennan, J. A. Zhang, and P. G. Black, 2007: Turbulent fluxes in the hurricane boundary layer. Part I: Momentum flux.

,*J. Atmos. Sci.***64**, 1089–1102.Garman, K. E., and Coauthors, 2006: An airborne and wind tunnel evaluation of a wind turbulence measurement system for aircraft-based flux measurements.

,*J. Atmos. Oceanic Technol.***23**, 1696–1708.Garratt, J. R., 1977: Review of drag coefficients over oceans and continents.

,*Mon. Wea. Rev.***105**, 915–929.Garratt, J. R., 1992:

*The Atmospheric Boundary Layer.*Cambridge University Press, 316 pp.Geernaert, G. L., 1990: Bulk parameterizations for the wind stress and heat flux.

*Surface Waves and Fluxes,*Vol. 1, G. L. Geernaert and W. L. Plant, Eds., Kluwer, 91–172.Grachev, A. A., and C. W. Fairall, 2001: Upward momentum transfer in the marine boundary layer.

,*J. Phys. Oceanogr.***31**, 1698–1711.Grachev, A. A., C. W. Fairall, J. E. Hare, J. B. Edson, and S. D. Miller, 2003: Wind stress vector over ocean waves.

,*J. Phys. Oceanogr.***33**, 2408–2429.Grachev, A. A., E. L Andreas, C. W. Fairall, P. S. Guest, and P. O. G. Persson, 2007a: SHEBA flux-profile relationships in the stable atmospheric boundary layer.

,*Bound.-Layer Meteor.***124**, 315–333.Grachev, A. A., E. L Andreas, C. W. Fairall, P. S. Guest, and P. O. G. Persson, 2007b: On the turbulent Prandtl number in the stable atmospheric boundary layer.

,*Bound.-Layer Meteor.***125**, 329–341.Ingel, L. Kh., 2011: On the effect of spray on the dynamics of the marine atmospheric surface layer in strong winds.

,*Izv. Atmos. Ocean. Phys.***47**, 119–127.Janssen, J. A. M., 1997: Does wind stress depend on sea-state or not?—A statistical error analysis of HEXMAX data.

,*Bound.-Layer Meteor.***83**, 479–503.Jarosz, E., D. A. Mitchell, D. W. Wang, and W. J. Teague, 2007: Bottom-up determination of air–sea momentum exchange under a major tropical cyclone.

,*Science***315**, 1707–1709.Johnson, H. K., J. Højstrup, H. J. Vested, and S. E. Larsen, 1998: On the dependence of sea surface roughness on wind waves.

,*J. Phys. Oceanogr.***28**, 1702–1716.Jones, I. S. F., and Y. Toba, Eds., 2001:

*Wind Stress over the Ocean.*Cambridge University Press, 307 pp.Khelif, D., S. P. Burns, and C. A. Friehe, 1999: Improved wind measurements on research aircraft.

,*J. Atmos. Oceanic Technol.***16**, 860–875.Khelif, D., C. A. Friehe, H. Jonsson, Q. Wang, and K. Rados, 2005: Wintertime boundary-layer structure and air–sea interaction over the Japan/East Sea.

,*Deep-Sea Res. II***52**, 1525–1546.Kitaigorodskii, S. A., and Yu. A. Volkov, 1965: On the roughness parameter of the sea surface and the calculation of momentum flux in the near-water layer of the atmosphere.

,*Izv. Atmos. Ocean. Phys.***1**, 566–574.Klipp, C. L., and L. Mahrt, 2004: Flux-gradient relationship, self-correlation and intermittency in the stable boundary layer.

,*Quart. J. Roy. Meteor. Soc.***130**, 2087–2103.Kraus, E. B., 1968: What we do not know about the sea-surface wind stress.

,*Bull. Amer. Meteor. Soc.***49**, 247–253.Kraus, E. B., and J. A. Businger, 1994:

*Atmosphere–Ocean Interaction.*2nd ed. Oxford University Press, 362 pp.Kudryavtsev, V. N., 2006: On the effect of sea drops on the atmospheric boundary layer.

,*J. Geophys. Res.***111**, C07020, doi:10.1029/2005JC002970.Lenschow, D. H., 1986: Aircraft measurements in the boundary layer.

*Probing the Atmospheric Boundary Layer,*D. H. Lenschow, Ed., Amer. Meteor. Soc., 39–55.Lighthill, J., 1999: Ocean spray and the thermodynamics of tropical cyclones.

,*J. Eng. Math.***35**, 11–42.Mahrt, L., and D. Khelif, 2010: Heat fluxes over weak SST heterogeneity.

,*J. Geophys. Res.***115**, D11103, doi:10.1029/2009JD013161.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.Mahrt, L., D. Vickers, E. L Andreas, and D. Khelif, 2012: Sensible heat flux in near-neutral conditions over the sea.

,*J. Phys. Oceanogr.***42**, 1134–1142.Makin, V. K., 2005: A note on the drag of the sea surface at hurricane winds.

,*Bound.-Layer Meteor.***115**, 169–176.Melville, W. K., 1977: Wind stress and roughness length over breaking waves.

,*J. Phys. Oceanogr.***7**, 702–710.Monin, A. S., and A. M. Yaglom, 1971:

*Statistical Fluid Mechanics: Mechanics of Turbulence.*Vol. 1. MIT Press, 769 pp.Moon, I.-J., T. Hara, I. Ginis, S. E. Belcher, and H. L. Tolman, 2004: Effect of surface waves on air–sea momentum exchange. Part I: Effect of mature and growing seas.

,*J. Atmos. Sci.***61**, 2321–2333.Moon, I.-J., I. Ginis, T. Hara, and B. Thomas, 2007: A physics-based parameterization of air–sea momentum flux at high wind speeds and its impact on hurricane intensity predictions.

,*Mon. Wea. Rev.***135**, 2869–2878.Mueller, J. A., and F. Veron, 2009: Nonlinear formulation of the bulk surface stress over breaking waves: Feedback mechanisms from air-flow separation.

,*Bound.-Layer Meteor.***130**, 117–134.Neumann, G., 1956: Wind stress on water surfaces.

,*Bull. Amer. Meteor. Soc.***37**, 211–217.Neumann, G., and W. J. Pierson Jr., 1966:

*Principles of Physical Oceanography.*Prentice-Hall, 545 pp.Nicholls, S., and C. J. Readings, 1979: Aircraft observations of the structure of the lower boundary layer over the sea.

,*Quart. J. Roy. Meteor. Soc.***105**, 785–802.Oost, W. A., G. J. Komen, C. M. J. Jacobs, and C. Van Oort, 2002: New evidence for a relation between wind stress and wave age from measurements during ASGAMAGE.

,*Bound.-Layer Meteor.***103**, 409–438.Panofsky, H. A., and J. A. Dutton, 1984:

*Atmospheric Turbulence: Models and Methods for Engineering Applications.*John Wiley and Sons, 397 pp.Paulson, C. A., 1970: The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer.

,*J. Appl. Meteor.***9**, 857–861.Persson, P. O. G., J. E. Hare, C. W. Fairall, and W. D. Otto, 2005: Air-sea interaction processes in warm and cold sectors of extratropical cyclonic storms observed during FASTEX.

,*Quart. J. Roy. Meteor. Soc.***131**, 877–912.Petersen, G. N., and I. A. Renfrew, 2009: Aircraft-based observations of air–sea fluxes over Denmark Strait and the Irminger Sea during high wind speed conditions.

,*Quart. J. Roy. Meteor. Soc.***135**, 2030–2045.Powell, M. D., P. J. Vickery, and T. A. Reinhold, 2003: Reduced drag coefficient for high wind speeds in tropical cyclones.

,*Nature***422**, 279–283.Roll, H. U., 1965:

*Physics of the Marine Atmosphere.*Academic Press, 426 pp.Romero, L., and W. K. Melville, 2010: Airborne observations of fetch-limited waves in the Gulf of Tehuantepec.

,*J. Phys. Oceanogr.***40**, 441–465.Rutgersson, A., A.-S. Smedman, and U. Högström, 2001: Use of conventional stability parameters during swell.

,*J. Geophys. Res.***106**, 27 117–27 134.Sanford, T. B., J. F. Price, J. B. Girton, and D. C. Webb, 2007: Highly resolved observations and simulations of the ocean response to a hurricane.

,*Geophys. Res. Lett.***34**, L13604, doi:10.1029/2007GL029679.Shpund, J., M. Pinsky, and A. Khain, 2011: Microphysical structure of the marine boundary layer under strong wind and spray formation as seen from simulations using a 2D explicit microphysical model. Part I: The impact of large eddies.

,*J. Atmos. Sci.***68**, 2366–2384.Smedman, A.-S., M. Tjernström, and U. Högström, 1994: The near-neutral marine atmospheric boundary layer with no surface shearing stress: A case study.

,*J. Atmos. Sci.***51**, 3399–3411.Smith, S. D., 1980: Wind stress and heat flux over the ocean in gale force winds.

,*J. Phys. Oceanogr.***10**, 709–726.Smith, S. D., 1988: Coefficients for sea surface wind stress, heat flux, and wind profiles as a function of wind speed and temperature.

,*J. Geophys. Res.***93**, 15 467–15 472.Smith, S. D., and Coauthors, 1992: Sea surface wind stress and drag coefficients: The HEXOS results.

,*Bound.-Layer Meteor.***60**, 109–142.Soloviev, A., and R. Lukas, 2010: Effects of bubbles and sea spray on air-sea exchange in hurricane conditions.

,*Bound.-Layer Meteor.***136**, 365–376.Stull, R. B., 1988:

*An Introduction to Boundary Layer Meteorology.*Kluwer, 666 pp.Sun, J., J. F. Howell, S. K. Esbensen, L. Mahrt, C. M. Greb, R. Grossman, and M. A. LeMone, 1996: Scale dependence of air–sea fluxes over the western equatorial Pacific.

,*J. Atmos. Sci.***53**, 2997–3012.Sun, J., D. Vandemark, L. Mahrt, D. Vickers, T. Crawford, and C. Vogel, 2001: Momentum transfer over the coastal zone.

,*J. Geophys. Res.***106**, 12 437–12 448.Sutton, O. G., 1953:

*Micrometeorology.*McGraw-Hill, 333 pp.Sverdrup, H. U., M. W. Johnson, and R. H. Fleming, 1942:

*The Oceans.*Prentice-Hall, 1087 pp.Toba, Y., S. D. Smith, and N. Ebuchi, 2001: Historical drag expressions.

*Wind Stress over the Ocean,*I. S. F. Jones and Y. Toba, Eds., Cambridge University Press, 35–53.von Arx, W. S., 1967:

*An Introduction to Physical Oceanography.*Addison-Wesley, 422 pp.Wieringa, J., 1993: Representative roughness parameters for homogeneous terrain.

,*Bound.-Layer Meteor.***63**, 323–363.Wilson, B. W., 1960: Note on surface wind stress over water at low and high wind speeds.

,*J. Geophys. Res.***65**, 3377–3382.Wu, J., 1969: Wind stress and surface roughness at air-sea interface.

,*J. Geophys. Res.***74**, 444–455.Wu, J., 1980: Wind-stress coefficients over sea surface near neutral conditions—a revisit.

,*J. Phys. Oceanogr.***10**, 727–740.Wyngaard, J. C., 2010:

*Turbulence in the Atmosphere.*Cambridge University Press, 393 pp.Zhang, J. A., W. M. Drennan, P. G. Black, and J. R. French, 2009: Turbulence structure of the hurricane boundary layer between the outer rainbands.

,*J. Atmos. Sci.***66**, 2455–2467.Zilitinkevich, S. S., 1969: On the computation of the basic parameters of the interaction between the atmosphere and the ocean.

,*Tellus***21**, 17–24.