TMCnet News
Fallacies of the Enthalpy Transfer Coefficient over the Ocean in High Winds [Journal of the Atmospheric Sciences](Journal of the Atmospheric Sciences Via Acquire Media NewsEdge) ABSTRACT Mesoscale and large-scale atmospheric models use a bulk surface flux algorithm to compute the turbulent flux boundary conditions at the bottom of the atmosphere from modeled mean meteorological quantities such as wind speed, temperature, and humidity. This study, on the other hand, uses a state-of-the-art bulk air-sea flux algorithm in stand-alone mode to compute the surface fluxes of momentum, sensible and latent heat, and enthalpy for a wide range of typical (though randomly generated) meteorological conditions over the open ocean. The flux algorithm treats both interfacial transfer (controlled by molecular processes right at the air-sea interface) and transfer mediated by sea spray. Because these two transfer routes obey different scaling laws, neutral-stability, 10-m transfer coefficients for enthalpy CKN10, latent heat CEN10, and sensible heat CHN10 are quite varied when calculated from the artificial flux data under the assumption of only interfacial transfer-the assumption in almost all analyses of measured air-sea fluxes. That variability increases with wind speed because of increasing spray-mediated transfer and also depends on surface temperature and atmospheric stratification. The analysis thereby reveals as fallacious several assumptions that are common in air-sea interaction research-especially in high winds. For instance, CKN10, CEN10, and CHN10 are not constants; they are not even single-valued functions of wind speed, nor must they increase monotonically with wind speed if spray-mediated transfer is important. Moreover, the ratio CKN10/CDN10, where CDN10 is the neutral-stability, 10-m drag coefficient, does not need to be greater than 0.75 at all wind speeds, as many have inferred from Emanuel's seminal paper in this journal. Data from the literature and from the Coupled Boundary Layers and Air-Sea Transfer (CBLAST) hurricane experiment tend to corroborate these results. (ProQuest: ... denotes formulae omitted.) 1. Introduction Sixteen years ago on these pages, Emanuel (1995) reported calculations of hurricane intensity based on his simple, balanced axisymmetric model of tropical cyclones. His specific objective was to build on the work by Ooyama (1969) and Rosenthal (1971) and investigate how sensitive storm intensity was to the surface exchange coefficients for enthalpy CK and momentum CD. He summarized his findings in the following key statements: The direct implication of this result ... is that the ratio CK/CD is most likely to lie in the range 1.2-1.5 in the high wind region of these intense storms. In no event are the results from either model consistent with values of CK/CD less than about three-fourths. This paper was pivotal for air-sea interaction research because it focused attention on the need to better understand the physics of air-sea exchange in high winds: The data and air-sea flux parameterizations available in 1995, which treated winds up to only about 20 m s-1, were not compatible with Emanuel's (1995) constraint on CK/CD. Theory and observations in moderate winds put CK/CD at values much less than 0.75 and thus spurred researchers to understand the discrepancy. Emanuel's (1995) paper, however, also had some unintended adverse effects. Although models had employed surface transfer coefficients for many years (e.g., Ooyama 1969; Benoit et al. 1989; Fantini and Buzzi 1993), some hurricane modelers took Emanuel's paper as further evidence that using an enthalpy transfer coefficient is the best way to model air-sea exchange even in high winds and, furthermore, presumed that the enthalpy transfer coefficient is constant - or, perhaps, a single-valued function of wind speed. These inferences have impeded progress on storm modeling because all are fallacious. Here is a list of misconceptions that have been inferred from Emanuel's (1995) paper and that otherwise exist in the air-sea interaction community when the discussion involves high winds (e.g., Braun and Tao 2000; Bao et al. 2002; Smith 2003; Zedier et al. 2009). Notice that some of these assumptions are mutually exclusive, an observation that emphasizes the persistent uncertainty in aspects of modeling air-sea fluxes. I will address all of these misconceptions. * CK (and CH, CE) is constant (CH and CE are comparable transfer coefficients for sensible and latent heat). * CK (and CH, CE) is a single-valued function of wind speed. * CK (and CH, CE) should increase monotonically with wind speed when sea spray becomes a significant agent for air-sea exchange. * CK/CD is constant. * CK/CD is a single-valued function of wind speed. * Observations suggest that Emanuel's (1995) guideline, CK/CD = 0.75, is invalid. In high winds, when copious amounts of sea spray fill the near-surface air, two routes exist by which enthalpy crosses the air-sea interface. One is through exchange controlled by molecular processes right at the air-sea interface (the interfacial route); a second is through spray-mediated processes (the spray route). Because the two transfer processes scale differently with wind speed and other mean meteorological variables, the enthalpy transfer coefficient is not single-valued and thus cannot be used alone to predict the air-sea enthalpy flux in high winds. Here, I demonstrate these ideas by using a state-ofthe-art bulk flux algorithm to compute total air-sea enthalpy fluxes for a wide range of typical environmental conditions that I generate from random numbers. That algorithm (Andreas et al. 2008; Andreas 2010) treats both the interfacial and spray routes by which enthalpy crosses the air-sea interface; hence, the total fluxes referred to above are the sums of both contributions. The CK values inferred from the computed fluxes are not single-valued functions of wind speed but rather range over a factor of 2 at the lower limits of hurricane-strength winds because of a dependence on surface temperature. The ratio CK/CD likewise has a spread in values for any given wind speed, and that spread increases with increasing wind speed. This CK/CD ratio is below Emanuel's (1995) limit of 0.75 for wind speeds between 12 and 33 m s-1 but rises above 0.75 and heads toward Emanuel's preferred range of 1.2-1.5 for the high winds in intense storms, as he suggests. Finally, these calculations are compatible with aircraft fluxes measured during the Coupled Boundary Layers and Air-Sea Transfer experiment (CBLAST; e.g., Black et al. 2007) and thus put in context measurements that Drennan et al. (2007) and Zhang et al. (2008) reported. 2. Approach a. Flux algorithm For wind speeds above about 5 m s-1, spray droplets are always present over the sea; their numbers increase roughly as the third power of the wind speed (e.g., Wu 1979; Monahan and Muircheartaigh 1980; Monahan et al. 1986). These spray droplets provide another route by which sensible heat, water vapor, and enthalpy cross the air-sea interface. Figure 1 shows my conceptual picture of spray-mediated transfer. Sensible and latent heat are always crossing the sea surface right at the interface. I call these the interfacial fluxes and denote them as Hs and HL, respectively. In my convention, positive fluxes are upward-from sea to air. Spray droplets provide another transfer route. The droplets originate with the same temperature as the sea surface but cool rapidly to an equilibrium temperature Teq below the air temperature. In so doing, they give up sensible heat. The droplets also give up water vapor as they evaporate; but because this is a slower process, the droplets are at temperature Teq for most of their evaporation. As a result, they convert sensible heat from the air to the latent heat of evaporation. These spray processes generally take place near the sea surface in the droplet evaporation layer (Smith 1989; Andreas et al. 1995). A rough estimate is that this layer extends about one significant wave height above mean sea level. The total sensible and latent heat fluxes at the top of this layer - Hs,T and HL,T, respectively - include the contributions from both interfacial and spray transfer. The sum of these fluxes is the total enthalpy flux Qen,T (e.g., Businger 1982), while the sum of Hs and HL is the interfacial enthalpy flux Qen,int - The quantities Hs,T, HL,T, and Qen,T can serve as the lower flux boundary conditions for atmospheric models. Andreas et al. (2008) and Andreas (2010) developed an algorithm to estimate all of these fluxes and tuned algorithm with data from the experiment to study Humidity Exchange over the Sea (HEXOS) and from the Fronts and Atlantic Storm Track Experiment (FASTEX) In simple form, the algorithm is ... (2-1a) ... (2.1b) ... (2.1c) Here, QS,sp, QL,sp, and Qen,sp are the spray-mediated contributions to the sensible heat, latent heat, and enthalpy fluxes. Also, ... (2.2) Because Qen,sp is tuned separately with the HEXOS and FASTEX data, it is close to but not always equal to QS,sp + QL,sp. The flux algorithm computes the interfacial fluxes in (2.1a) and (2.1b) as (Andreas et al. 2008) ... (2-3a) ... (2-3b) Here, ?a is the air density; cp, the specific heat of air at constant pressure; and Lv, the latent heat of vaporization. Also in (2.3), Sr, Tr, and Qr are the effective wind speed, potential temperature, and specific humidity, respectively, at reference height r (assumed to be above the droplet evaporation layer); and Ts and Qs are the sea surface values of temperature and specific humidity. The transfer coefficients for sensible and latent heat in (2.3), respectively CHr and CEr, come essentially from the Coupled Ocean-Atmosphere Response Experiment (COARE) version 2.6 algorithm (Fairall et al. 1996), with a few small modifications as discussed in Andreas and Emanuel (2001), Andreas and DeCosmo (2002), Perrie et al. (2005), and Andreas et al. (2008). Andreas et al. (2008) deduced simple parameterizations for QS,sp and QL,sp in (2.1a) and (2.1b): ... (2.4a) ... (2.4b) In these equations, ?s is the density of seawater, cps is the specific heat of seawater at constant pressure, Teq,100 is the equilibrium temperature of spray droplets that have a radius of 100 µm when they form, and rd (tf50) is the radius when they fall back into the sea of droplets that form with 50-µm radius. This radius comes from ... (2.5) where req50 is the equilibrium radius of droplets that form with radius 50 µm, tr,50 is the e-folding time to reach that equilibrium radius through evaporation, and tf50 is the residence time in the droplet evaporation layer of droplets formed with 50-µm radius. See Andreas (2005a) for equations to calculate the quantities Teq,100 req,50, and tr,50. Lastly, in (2.4), u* is the friction velocity, which the algorithm computes as ... (2.6) Here, Sr and CDr, the drag coefficient for height r, come from the COARE version 2.6 algorithm, with the small modifications mentioned above. In (2.3) and (2.4), the fluxes are in watts per square meter when the other quantities have MKS units; in (2.4b) and (2.5), however, the radii are in micrometers. In (2.6), u* and Sr are in meters per second. From the HEXOS and FASTEX data, Andreas (2010) derived an equation similar in form to (2.4a) for the spray enthalpy flux in (2.1c): ... (2.7) This result is consistent with Andreas and Emanuel's (2001) conclusion that, if sea spray affects the net air-sea enthalpy flux, the mechanism must be through an exchange of sensible heat. b. Implied transfer coefficients If the total enthalpy flux, (2.1c), were parameterized in terms of CK (i.e., under the assumption of only interfacial exchange), the expression would be (cf. Emanuel 1995) ... (2.8) Here, cpd is the specific heat of dry air at constant pressure. Because the Andreas (2010) algorithm computes Qen,T for specified mean meteorological conditions such as Ts, Tr, Qs, Qr, and wind speed Ur, I can invert (2.8) to compute the implied enthalpy transfer coefficient: ... (2.9) Ultimately, however, comparing the so-called neutralstability transfer coefficients is more meaningful. For a standard reference height of 10 m, these values come from (e.g., Andreas and Murphy 1986) ... (2.10) Here, k (=0.40) is the von Kármán constant; CDr is the same value used in (2.6); CKr comes from (2.9); ?h (r/L) is a stratification correction that appears in the near-surface temperature and humidity profiles, where L is the Obukhov length; and ... (2.11) Here, Z0 is the roughness length for wind speed from the COARE version 2.6 algorithm: ... (2.12) where a (=0.0185) is the Charnock constant, v is the kinematic viscosity of air (m2 s-1), and g is the acceleration of gravity (m s-2). In (2.10)-(2.12), r and z0 must be in meters. The same set of mean meteorological variables yields Hs,T and HL,T from (2.1a) and (2.1b). Hence, I can likewise obtain implied values of the neutral-stability, 10-m transfer coefficients for sensible and latent heat, CHN10 and CEN10, under the assumption of only interfacial transfer by inverting (2.3) after first replacing Hs and HL with Hs,T and HL,T, respectively. The conversions from CHr to CHN10 and from CEr to CEN10 also follow (2.10). c. Parameter space The final step in my approach is to generate sets of mean meteorological variables to use as input for the flux algorithm. I select these meteorological variables by using a random number generator to create uniform distributions of five key variables within a broad range of marine conditions. Without losing generality, and to eliminate the need to convert from height r to 10 m in (2.10), I first assume that wind speed U10, air temperature T10, and relative humidity (RH) are all for a reference height of 10 m. The Andreas et al. (2008) and Andreas (2010) algorithm is theoretically based and tuned with data for wind speeds up to 20 m s-1; therefore, I presume that it is reasonably accurate for U10 values between 5 and 40 m s-1. Nevertheless, I assign U10 by generating a random number between 0 and 1 and linearly mapping this to a wind speed between 5 and 50 m s-1; but I subsequently ignore the cases with wind speeds from 40 to 50 m s-1. Some of the "data" selection functions described below led to spurious fluxes at the upper end of the wind speed range. Letting the allowable wind speed be above the 40 m s-1 limit in my algorithm but ignoring cases with winds above that limit eliminated these end effects. Likewise, I generate another random number between 0 and 1 and linearly map it to a sea surface temperature Ts between 0° and 30°C. I next tie the 10-m air temperature to this surface temperature by assuming that Ts - T10 will usually be in the interval ?Tmin = -1°C to dTmax = 3°C (e.g., Jackson and Wick 2010). Certainly, ?T can be outside this range-under cold-air advection off continents or off regions covered with sea ice, for instance. But this limited range is sufficient to demonstrate my points. A wider range would have produced even more variability in the figures I present soon. I also assume that the Ts - T10 interval narrows with increasing wind speed and calculate the actual sea-air temperature difference as ... (2.13) where ? is another random number between 0 and 1. Notice that when U10 = 50 m s-1, the sea-air temperature difference is 0°C. In essence, (2.13) is a way to parameterize the feedback between the surface fluxes and the air temperature. I also adjust the minimum sea-air temperature difference for wind speed according to ... (2.14) This step is necessary to make both the minimum and maximum Ts - T10 limits approach 0°C as U10 approaches 50 m s-1 because now ... (2.15) Last, I convert T10 to the potential temperature at 10 m T10. The surface salinity of the open ocean typically is 32-37 psu. Again, I generate another random number between 0 and 1 and map this linearly to a salinity S from 32 to 37 psu. With Ts and S now specified, I compute Qs as the saturation specific humidity for seawater with temperature Ts and salinity S. The Andreas et al. (2008) and Andreas (2010) algorithm is formulated with relative humidity, rather than specific humidity, as the input humidity variable. Once RH is specified, however, the specific humidity at 10 m Q10 is easy to obtain. I assume that relative humidity increases with increasing wind speed. As with temperature, this assumption is a simple way to treat the feedback between the latent heat flux and the near-surface humidity. The lowest relative humidity that I consider (RHmin) is 75%; the highest is constrained by saturation with seawater of salinity 5 (e.g., Roll 1965, p. 262): ... (2.16) The mean relative humidity for U10 is thus ... (2.17) I further assume the maximum variability about this mean humidity is ±3%. Hence, I choose the actual maximum deviation about RHmean by finding another random number ? between 0 and 1 and calculating ... (2.18) Equation (2.18) centers the deviation in relative humidity on 0% . Next, I assume that this deviation will also decrease with increasing wind speed and calculate the relative humidity at 10 m as ... (2.19) Last, I check that RH is not less than RHmin and not greater than RHmax. If it is, I set it to the respective limit. The Andreas et al. (2008) and Andreas (2010) algorithm also requires the air pressure P and the water depth D. The depth occurs only in the module that computes the significant wave height, which the algorithm uses in turn to calculate the residence time tf50 in (2.4b) and (2.5). Because the significant wave height is not very sensitive to water depth for typical open ocean depths (Andreas and Wang 2007), I set D = 3000 m for all calculations. I likewise fix P at 1000 mb because barometric pressure only weakly influences specific humidity and air density in these calculations. Table 1 summarizes the nominal range of all of these input variables. Andreas (2005b) gives the equations I use to evaluate the dynamic and thermodynamic quantities such as ?a, cp, Lv, v, Qs, and Q10. The steps that I have just described produce one set of input variables for forcing the Andreas et al. (2008) and Andreas (2010) algorithm and thus yield one set of flux estimates. I repeated this sequence 3000 times to produce a reliable and representative set of fluxes and the resulting transfer coefficients. After I ignored cases with UN10 above 40 m s-1, 2327 flux sets remained. 3. Results Figure 2 shows the neutral-stability enthalpy transfer coefficient for a reference height of 10 m CKN10 that the 3000 calculations described in the last section yielded. For comparison, the figure also shows the curve that represents interfacial transfer in the Andreas et al. (2008) and Andreas (2010) algorithm. The interfacial transfer coefficient for enthalpy in neutral stability is a singlevalued function of wind speed, but the implied enthalpy transfer coefficient computed from the total flux is not: CKN10 is a function of surface temperature and wind speed (and, less obviously, relative humidity and surface salinity). The CKN10 values in Fig. 2 are approximately constant for wind speeds up to about 20 m s-1. Such behavior in actual measurements of CKN10 (or CHN10 and CEN10, the comparable neutral-stability transfer coefficients for sensible and latent heat) is usually interpreted to mean that sea spray has no effect on air-sea heat transfer (e.g., Smith 1989; DeCosmo et al. 1996; Smith et al. 1996). This misconception seems to go back at least to Bortkovskii (1987), who suggested in his Fig. 3.10 that CHN10 and CEN10 increase dramatically for winds above 15 m s-1 because of enhanced transfer mediated by sea spray. Figure 2, however, suggests a different scenario. The COARE version 2.6 algorithm, which uses the surface renewal theory of Liu et al. (1979) to estimate the interfacial transfer coefficients, predicts that CKN10, CHN10, and CEN10 all decrease with increasing wind speed beyond their maxima at about 6 m s-1. Thus, the transfer coefficients need not increase to reflect spray effects; the nearly constant CKN10 values for 5 = UN10 = 20 m s-1 in Fig. 2 are above the interfacial transfer limit because of spray-mediated exchange (cf. Andreas and DeCosmo 2002; Andreas et al. 2008; Andreas 2010). Figure 3 shows what the Andreas algorithm says about the ratio CKN10/CDN10. Although CDN10 is a singlevalued function of wind speed in this algorithm, CKN10/ CDN10 is not because of how spray-mediated transfer affects CKN10. The figure also shows Emanuel's (1995) limit - that CKN10/CDN10 should be larger than 0.75 - and what strict interfacial transfer predicts for the value of that ratio. This interfacial curve explains why Emanuel's limit caused such concern: The CKN10/CDN10 ratio for interfacial transfer is below 0.75 for all wind speeds above about 10 m s-1 and is roughly half this limit for hurricanestrength winds. After Emanuel's (1995) paper was published, Bister and Emanuel (1998) identified dissipative heating in the atmospheric boundary layer in the high-wind core of tropical cyclones as an additional source of heat that increases the intensity of these storms (cf. Businger and Businger 2001; Smith 2003). Emanuel (2010, personal communication) thus speculates that his limiting CKN10 CDN10 ratio could actually be significantly lower than 0.75. Bryan and Rotunno (2009) recently upgraded the Rotunno and Emanuel (1987) nonhydrostatic, axisymmetric model, including adding dissipative heating. This Rotunno-Emanuel model is one of the hurricane models that Emanuel (1995) used in his sensitivity study. Bryan and Rotunno found that, in their simulations, storm intensity still depended on the CKN10/CDN10 ratio but was also very sensitive to the parameterization for horizontal turbulent diffusion. With their preferred horizontal turbulent mixing length of 1500 m, however, CKN10/CDN10 still needed to be above 0.75 to produce a category-5 hurricane. Hence, I will continue to take 0.75 as "Emanuel's limit." Even with spray-mediated transfer, the CKN10/CDN10 ratio in Fig. 3 is less than 0.75 for wind speeds from about 12 up to 36 m s-1 for some surface temperatures. Beyond this upper limit, though, CKN10/CDN10 is greater than 0.75 - a result that is entirely compatible with Emanuel's (1995) recommendation: CKN10/CDN10 should be above 0.75 in "the high wind region" of "intense storms." For any wind speed in Fig. 3, the CKN10.CDN10 ratio generally decreases with increasing surface temperature. Using similarity arguments and the assumption that spray processes dominate both air-sea enthalpy and momentum exchange in very high winds, Emanuel (2003) made a similar prediction. In his theory, the temperature variable is air temperature and neither Ck nor CD depends on wind speed in his self-similar regime; hence, the predictions in his Fig. 3 and the results in my Fig. 3 are not exactly comparable. Nevertheless, to demonstrate that Emanuel's similarity arguments are compatible with my theory of spray-mediated transfer, I read from Emanuel's Fig. 3 values for Ck/Cd of 1.27 and 0.85 for temperatures of 20° and 30°C, respectively, to compare with CKN10/CDN10 values of 1.02 and 0.80 at 40 m s-1 for the same temperatures in my Fig. 3. Figure 4 focuses more closely on tropical cyclones. This figure repeats Fig. 3 but shows only the surface temperature range typical of tropical storms. The plot also includes aircraft measurements of CKN10/CDN10 during the CBLAST hurricane experiment in 2003 (Zhang et al. 2008; see also French et al. 2007; Drennan et al. 2007). Three features stand out in Fig. 4. First, most of the aircraft data are above the curve for strict interfacial transfer; thus, the CBLAST data probably document significant spray-mediated transfer. Second, predictions from Andreas's (2010) algorithm neatly cut through the CBLAST data and thus seem to capture the magnitude and wind speed dependencies of CKN10 and CDN10 Third, most of the CBLAST data lie below Emanuel's (1995) limiting value of 0.75. Zhang et al. (2008) therefore concluded that Emanuel's model needs to be "revisited." But the Andreas (2010) algorithm demonstrates that CKN10/CDN10 values less than 0.75 are compatible with theory in the wind speed range of the CBLAST data. Figure 5 shows algorithm predictions of the implied neutral-stability latent heat transfer coefficient for a reference height of 10 m CEN10. Again, CEN10 is basically constant for wind speeds up to 20 m s-1, as in the HEXOS data (DeCosmo et al. 1996). Even in this range, though, CEN10 is not single-valued: Cenio spreads out in response to the surface temperature (among other variables). Furthermore, that spread increases with increasing wind speed. Figure 6 reproduces from Fig. 5 just the cases for which the surface temperature is typical of tropical cyclones because the figure also shows the CBLAST CEN10 values measured by aircraft in Hurricanes Fabian and Isabel in 2003 (French et al. 2007; Drennan et al. 2007). These CBLAST values tend to cluster around the line representing strict interfacial transfer - a result that is a bit curious in light of Fig. 4. The average of the CBLAST points in Fig. 6 is CEN10 = 1-18 × 10-3, and Drennan et al. (2007) ascribe an uncertainty to this average of ±25%. As a result, the data are not precise enough to rule out spray contributions. In the wind speed range for most of the CBLAST data, 18 =UN10 <25 m s-1, the CEN10 values in Fig. 6 from the flux algorithm average 1.49 × 10-3. This average is not statistically different for the average of the CBLAST data in this range, 1.19 × 10-3 ± 25%. Another reason for the different implications of Figs. 4 and 6 is that the CBLAST values of CDN10 used in calculating the CBLAST CKN10/CDN10 values in Fig. 4 are generally smaller than CDN10 in the Andreas et al. (2008) and Andreas (2010) algorithm. This algorithm calculates Conio in (2.11) from the Charnock-like relation in the COARE version 2.6 algorithm, (2.12). Figure 7 is the final plot in this series. It shows values of the neutral-stability, 10-m coefficient for sensible heat transfer CHN10 implied by the dataset generated in section 2. The data segregate according to the sea-air potential temperature difference ... and their spread increases with increasing wind speed. Because of the way I selected the air temperature in section 2 - where ?T| got smaller with increasing wind speed - the midrange orange symbols ... in Fig. 7 extend only to about 34 m s-1, and the high-range red symbols ... extend only up to about 20 m s-1. Because spray droplets always cool to a temperature less than the air temperature, spray processes almost always transfer sensible heat from sea to air. The blue symbols in Fig. 7 are, thus, quite scattered. These symbols represent stable stratification ... which implies a negative interfacial sensible heat flux. The spray-mediated transfer can still be positive, however. Consequently, with spraymediated transfer, the total sensible heat flux can actually be countergrathent - a net positive heat flux despite a negative sea-air temperature difference. The variability in CHN10 in Fig. 7 caused by spraymediated exchange may explain why measurements of CHN10 over the ocean are typically more scattered than measurements of CEN10, especially in winds above 10 m s-1 (e.g., Large and Pond 1982; DeCosmo et al. 1996; Dupuis et al. 2003; Persson et al. 2005). Moreover, a main conclusion in Large and Pond (1982) is that CHN10 is significantly larger in unstable stratification than in stable stratification (cf. Smith 1980). Heretofore, this result has been unexplained because interfacial transfer theory predicts only minor differences between CHN10 and CEN10 and no dependence of either on atmospheric stratification. In Fig. 7, however, the average value of CHN10 is obviously smaller in stable stratification than it is in unstable stratification because of the spraymediated transfer. 4. Discussion and conclusions Because the nature of spray-mediated air-sea transfer has only recently emerged, many earlier concepts about parameterizing air-sea heat exchange now appear fallacious. Here I have revealed several such fallacies by generating an artificial dataset of mean meteorological quantities and then using a bulk flux algorithm (Andreas et al. 2008; Andreas 2010) to estimate the air-sea fluxes of momentum, sensible heat, latent heat, and enthalpy from these data. From these input variables and the resulting fluxes, I could calculate neutral-stability, 10-m values of the drag coefficient and the transfer coefficients for sensible heat, latent heat, and enthalpy - CDN10, CHN10, CCEN10, and CN10, respectively - using interfacial scaling. This is a process that most experimentalists follow once they have measured the fluxes. The bulk flux algorithm used to generate the fluxes parameterizes both interfacial transfer and spray-mediated transfer, however. Consequently, the resulting CEK10, CEN10, and CEN10 values do not exhibit strict interfacial scaling but rather are multivalued because spray-mediated heat exchange does not scale the same as interfacial exchange. To see this, compare (2.2) and (2.3) with (2.4) and (2.7). These analyses led to my conclusions that several common assumptions regarding CKN10, CEN10, and CHN10 are fallacious. Although CKN10 and CEN10 may be measured as nearly constant for 10-m, neutral-stability winds UN10 less than 20 m s-1, this behavior cannot be extrapolated to higher winds nor is it evidence that sea spray is playing no role in the air-sea exchanges of enthalpy and latent heat. Theory predicts that both CKN10 and CEN10 should decrease monotonically as wind speed increases from 6 to 20 m s-1 if interfacial transfer is the only exchange mechanism. The fact that observations of CKN10 and CEN10 do not decrease in this wind speed range is thus evidence of spray-mediated transfer that complements the interfacial transfer. Moreover, above UN10 = 20 m s-1, I predict that both CKN10 and CEN10 increase dramatically with wind speed because of the ever-increasing importance of spray-mediated transfer. The behavior of the inferred CHN10 is more complex because spray-mediated exchange almost always leads to a sea-to-air flux of sensible heat. Interfacial sensible heat transfer, on the other hand, is always down the potential temperature grathent. Hence, sometimes the interfacial and spray fluxes are complementary, but sometimes they are in opposite directions. The two transfer routes thus lead to CHN10 values that are much more variable than CKN10 and CEN10. In fact, the opposing fluxes can lead to negative values of CHN10 - that is, to countergrathent fluxes - which is a clear violation of strict interfacial transfer theory. As a result, I suspect that measurements of negative CHN10 (which are true) are often not reported because they fail to pass quality controls that (mistakenly) require sensible heat fluxes to be down the potential temperature grathent. Many reported measurements of CEN10 and CHN10 over the ocean have yielded plots with some of the general features of Figs. 5 and 7. Even in modest winds (where spray is still, nevertheless, making a contribution), measured CHN10 values are often more scattered than measured CEN10 values. The explanation until now has usually been that the sensible heat flux is a weaker signal over the ocean than is the latent heat flux; the transfer coefficient for sensible heat thus would have a larger uncertainty and more scatter. An equally plausible explanation now is that spray-mediated transfer confounds analyses that presume only interfacial transfer (i.e., of CHN10). When Large and Pond (1982) segregated their CHN10 data into measurements in stable and unstable stratification, they found CHN10 to be significantly larger in unstable stratification - a result that, for them, was "neither predicted nor expected." Figure 7, however, now predicts just this result in the wind speed range up to 20 m s-1, where most of Large and Pond's data fell. Meanwhile, Large and Pond and Anderson and Smith (1981) did not find a comparable dependence of CEN10 on stratification (cf. Fig. 5). Many have inferred from Emanuel's (1995) modeling that CKN10/CDN10 can be taken as a constant in hurricane models or, alternatively, that this ratio is a single-valued function of wind speed. Furthermore, many often assume that CKN10/CDN10 must be greater than 0.75 for all wind speeds though Emanuel specifically gave this limit as a constraint only "in the high wind region" of tropical cyclones. My analyses have likewise showed these three inferences all to be fallacious. Figure 3 shows that the theory and data tuning on which my flux algorithm (Andreas et al. 2008; Andreas 2010) is based predict CKN10/CDN10 values that are neither constant nor single-valued functions of wind speed. This ratio spreads out with increasing wind speed because interfacial and spray-mediated enthalpy transfers scale differently with wind speed, air temperature, surface temperature, relative humidity, and surface salinity. In Fig. 3, the ratios group and spread according to surface temperature, but these other environmental conditions also contribute to the spread. Figure 4 focuses on Emanuel's (1995) constraint of 0.75. My theory predicts that CKN10/Conio is les CDN10 s than 0.75 in gale-force winds (UN10 between 12 and 30 m s-1) for surface temperatures that favor tropical cyclones. But this ratio rises above 0.75 at the lower limit of hurricanestrength winds (CN10 of 30-36 m s-1); my theory is thus compatible with Emanuel's guidance. Aircraft flux measurements from CBLAST in Hurricanes Fabian and Isabel and resulting evaluations of CDN10, CEN10, and CKN10 (French et al. 2007; Drennan et al. 2007; Zhang et al. 2008) seemed incompatible with spray-mediated transfer and with Emanuel's (1995) prediction for CKN10/CDN10 when viewed alone. But when I presented the CBLAST data in the context of my predictions for CKN10/CFN10 and CEN10 in Figs. 4 and 6, those data refute neither my theory of spray-mediated transfer nor Emanuel's guidance that CKN10/CDN10 must be larger than 0.75 in hurricane winds. In brief, the CBLAST data do not include measurements in hurricane-strength winds but fall in the wind speed range 16-30 m s-1, where my theory suggests that CKN10/CDN10 should, indeed, be below 0.75, as the CBLAST data generally are. We now seem to know enough about air-sea exchange to move beyond the practice of parameterizing heat fluxes in terms of transfer coefficients that obey interfacial scaling - the gravest fallacy of all. Spray-mediated transfer is significant even in winds of only 12 m s-1 and does not obey interfacial scaling. Thus, no heat fluxes parameterized with CH, CE, and CK can be generally accurate, especially as the wind speed increases. The interfacial and spray transfer routes must be treated separately. Free FORTRAN code, based on Andreas et al. (2008) and Andreas (2010), is available online (at www.nwra.com/resumes/andreas/software.php) for making these calculations. Acknowledgments. I thank Jun Zhang and Will Drennan for generously providing their CBLAST data, Kerry Emanuel for valuable comments on an early version of the manuscript, two anonymous reviewers for their editorial guidance, and Emily Moynihan for creating Fig. 1. The U.S. Office of Naval Research supported this work with awards N00014-08-1-0411 and N00014-11-1-0073 to NorthWest Research Associates and through the National Ocean Partnership Program with award N00014-10-1-0154 to the University of Rhode Island, for whom I am a subcontractor. REFERENCES Anderson, R. J., and S. D. Smith, 1981: Evaporation coefficient for the sea surface from eddy flux measurements. J. Geophys. Res., 86, 449-456. Andreas, E. L, 2005a: Approximation formulas for the microphysical properties of saline droplets. Atmos. Res., 75, 323-345. _____ , 2005b: Handbook of Physical Constants and Functions for Use in Atmospheric Boundary Layer Studies. ERDC/CRREL Monograph M-05-1, U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, NH, 42 pp. _____ , 2010: Spray-mediated enthalpy flux to the atmosphere and salt flux to the ocean in high winds. J. Phys. Oceanogr, 40, 608-619. _____ , and B. Murphy, 1986: Bulk transfer coefficients for heat and momentum over leads and polynyas. J. Phys. Oceanogr., 16, 1875-1883. _____ , and K. A. Emanuel, 2001: Effects of sea spray on tropical cyclone intensity. J. Atmos. Sci., 58, 3741-3751. _____ , and J. DeCosmo, 2002: The signature of sea spray in the HEXOS turbulent heat flux data. Bound.-Layer Meteor., 103, 303-333. _____ , and S. Wang, 2007: Predicting significant wave height off the northeast coast of the United States. Ocean Eng., 34, 1328-1335. _____ , J. B. Edson, E. C. Monahan, M. P. Rouault, and S. D. Smith, 1995: The spray contribution to net evaporation from the sea: A review of recent progress. Bound. -Layer Meteor., 72, 3-52. _____ , P. O. G. Persson, and J. E. Hare, 2008: A bulk turbulent airsea flux algorithm for high-wind, spray conditions. J. Phys. Oceanogr., 38, 1581-1596. Bao, J.-W., S. A. Michelson, and J. M. Wilczak, 2002: Sensitivity of numerical simulations to parameterizations of roughness for surface heat fluxes at high winds over the sea. Mon. Wea. Rev., 130, 1926-1932. Benoit, R., J. Côté, and J. Mailhot, 1989: Inclusion of a TKE boundary layer parameterization in the Canadian regional finite-element model. Mon. Wea. Rev., 117, 1726-1750. Bister, M., and K. A. Emanuel, 1998: Dissipative heating and hurricane intensity. Meteor. Atmos. Phys., 65, 233-240. Black, P. G., and Coauthors, 2007: Air-sea exchange in hurricanes: Synthesis of observations from the Coupled Boundary Layer Air-Sea Transfer Experiment. Bull. Amer. Meteor. Soc., 88, 357-374. Bortkovskii, R. S., 1987: Air-Sea Exchange of Heat and Moisture during Storms. D. Reidel, 194 pp. Braun, S. A., and W.-K. Tao, 2000: Sensitivity of high-resolution simulations of Hurricane Bob (1991) to planetary boundary layer parameterizations. Mon. Wea. Rev., 128, 3941-3961. Bryan, G. H., and R. Rotunno, 2009: The maximum intensity of tropical cyclones in axisymmetric numerical model simulations. Mon. Wea. Rev., 137, 1770-1789. Businger, J. A., 1982: The fluxes of specific enthalpy, sensible heat and latent heat near the earth's surface. J. Atmos. Sci., 39, 1889-1892. Businger, S., and J. A. Businger, 2001: Viscous dissipation of turbulence kinetic energy in storms. J. Atmos. Sci., 58, 3793-3796. DeCosmo, J., K. B. Katsaros, S. D. Smith, R. J. Anderson, W. A. Oost, K. Bumke, and H. Chadwick, 1996: Air-sea exchange of water vapor and sensible heat: The Humidity Exchange over the Sea (HEXOS) results. J. Geophys. Res., 101, 12 001-12 016. 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. Dupuis, H., C. Guerin, D. Hauser, A. Weill, P. Nacass, W. M. Drennan, S. Cloché, and H. C. Graber, 2003: Impact of flow distortion corrections on turbulent fluxes estimated by the inertial dissipation method during the FETCH experiment on R/V L'Atalante. J. Geophys. Res., 108, 8064, doi:10.1029/ 2001JC001075. 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. _____ , 2003: A similarity hypothesis for air-sea exchange at extreme wind speeds. J. Atmos. Sci., 60, 1420-1428. 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. Fantini, M., and A. Buzzi, 1993: Numerical experiments on a possible mechanism of cyclogenesis in the Antarctic region. Tellus, 45A, 99-113. 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. Jackson, D. L., and G. A. Wick, 2010: Near-surface air temperature retrieval derived from AMSU-A and sea surface temperature observations. J. Atmos. Oceanic Technol., 27, 1769-1776. Large, W. G., and S. Pond, 1982: Sensible and latent heat flux measurements over the ocean. J. Phys. Oceanogr., 12, 464-482. Liu, W. T., K. B. Katsaros, and J. A. Businger, 1979: Bulk parameterization of air-sea exchanges of heat and water vapor including the molecular constraints at the interface. J. Atmos. Sci., 36, 1722-1735. Monahan, E. C., and I. Ó. Muircheartaigh, 1980: Optimal powerlaw description of oceanic whitecap coverage dependence on wind speed. J. Phys. Oceanogr., 10, 2094-2099. _____ , D. E. Spiel, and K. L. Davidson, 1986: A model of marine aerosol generation via whitecaps and wave disruption. Oceanic Whitecaps and Their Role in Air-Sea Exchange Processes, E. C. Monahan and G. Mac Niocaill, Eds., D. Reidel, 167-174. Ooyama, K., 1969: Numerical simulation of the life cycle of tropical cyclones. J. Atmos. Sci., 26, 3-40. Perrie, W., E. L Andreas, W. Zhang, W. Li, J. Gyakum, and R. McTaggart-Cowan, 2005: Sea spray impacts on intensifying midlatitude cyclones. J. Atmos. Sci., 62, 1867-1883. 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. Roll, H. U., 1965: Physics of the Marine Atmosphere. Academic Press, 426 pp. Rosenthal, S. L., 1971: The response of a tropical cyclone model to variations in boundary layer parameters, initial conditions, lateral boundary conditions, and domain size. Mon. Wea. Rev., 99, 767-777. Rotunno, R., and K. A. Emanuel, 1987: An air-sea interaction theory for tropical cyclones. Part II: Evolutionary study using a nonhydrostatic axisymmetric numerical model. J. Atmos. Sci., 44, 542-561. Smith, R. K., 2003: A simple model of the hurricane boundary layer. Quart. J Roy. Meteor. Soc., 129, 1007-1027. Smith, S. D., 1980: Wind stress and heat flux over the ocean in gale force winds. J. Phys. Oceanogr., 10, 709-726. _____ , 1989: Water vapor flux at the sea surface. Bound.-Layer Meteor., 47, 277-293. _____ , K. B. Katsaros, W. A. Oost, and P. G Mestayer, 1996: The impact of the HEXOS programme. Bound.-Layer Meteor., 78, 121-141. Wu, J., 1979: Oceanic whitecaps and sea state. J. Phys. Oceanogr., 9, 1064-1068. Zedier, S. E., P. P. Niiler, D. Stammer, E. Terrill, and J. Morzel, 2009: Ocean's response to Hurricane Frances and its implications for drag coefficient parameterization at high wind speeds. J. Geophys. Res., 114, C04016, doi:10.1029/2008JC005205. Zhang, J. A., P. G. Black, J. R. French, and W. M. Drennan, 2008: First direct measurements of enthalpy flux in the hurricane boundary layer: The CBLAST results. Geophys. Res. Lett, 35, L14813, doi:10.1029/2008GL034374. EDGAR L ANDREAS NorthWest Research Associates, Inc., Lebanon, New Hampshire (Manuscript received 25 October 2010, in final form 28 February 2011) Corresponding author address: Dr. Edgar L Andreas, NorthWest Research Associates, Inc. (Seattle Division), 25 Eagle Ridge, Lebanon, NH 03766-1900. E-mail: [email protected] (c) 2011 American Meteorological Society |