A Model for Air–Sea Interaction Bulk Coefficient over a Warm Mature Sea under Strong Wind

Naoto Kihara Civil Engineering Research Laboratory, Central Research Institute of Electric Power Industry, Abiko, Abiko-shi, Chiba, Japan

Search for other papers by Naoto Kihara in
Current site
Google Scholar
PubMed
Close
and
Hiromaru Hirakuchi Civil Engineering Research Laboratory, Central Research Institute of Electric Power Industry, Abiko, Abiko-shi, Chiba, Japan

Search for other papers by Hiromaru Hirakuchi in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

A boundary layer model for evaluating sensible and latent heat fluxes over a mature sea accounting for sea spray effects at wind speeds of up to 28 m s−1 is presented. Heat exchange across the ocean surface controls the development of tropical cyclones, and Emanuel’s theory suggests that the ratio of the exchange coefficient of total enthalpy to the drag coefficient should be greater than 0.75 to maintain the intensity of tropical cyclones. However, traditional bulk algorithms predict a monotonic decrease in this ratio with increasing wind speed, giving a value of less than 0.5 under tropical cyclone conditions. The present model describes a decrease in the ratio with increasing wind speed under weak to moderate winds (<20 m s−1), and a plateau at approximately 0.7 under strong winds (>20 m s−1).

Corresponding author address: Naoto Kihara, Civil Engineering Research Laboratory, Central Research Institute of Electric Power Industry, 1646 Abiko, Abiko-shi, Chiba 270-1194, Japan. kihara@criepi.denken.or.jp

Abstract

A boundary layer model for evaluating sensible and latent heat fluxes over a mature sea accounting for sea spray effects at wind speeds of up to 28 m s−1 is presented. Heat exchange across the ocean surface controls the development of tropical cyclones, and Emanuel’s theory suggests that the ratio of the exchange coefficient of total enthalpy to the drag coefficient should be greater than 0.75 to maintain the intensity of tropical cyclones. However, traditional bulk algorithms predict a monotonic decrease in this ratio with increasing wind speed, giving a value of less than 0.5 under tropical cyclone conditions. The present model describes a decrease in the ratio with increasing wind speed under weak to moderate winds (<20 m s−1), and a plateau at approximately 0.7 under strong winds (>20 m s−1).

Corresponding author address: Naoto Kihara, Civil Engineering Research Laboratory, Central Research Institute of Electric Power Industry, 1646 Abiko, Abiko-shi, Chiba 270-1194, Japan. kihara@criepi.denken.or.jp

1. Introduction

The development of tropical cyclones is governed by heat exchange between the air and the sea. Emanuel (1986) proposed a simple energy balance model for tropical cyclones in which the inner region of the tropical cyclone plays a role similar to that of a Carnot heat engine, extracting latent and sensible heat from the ocean and transferring that heat outward to overcome frictional dissipation in the boundary layer at large radii. Thus, extracting heat is seen as essential to maintain the intensity of a tropical cyclone.

In meteorological models, heat transfer from the sea surface is generally estimated using bulk algorithms, which have been investigated extensively using large datasets of observations. Fairall et al. (1996) proposed a method for flux estimation using bulk variables obtained by the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) for weak to moderate winds (0.5 ≤ u10 ≤ 10 m s−1). It was shown that the exchange coefficient for latent heat flux corrected for neutral stratification is an approximately constant value of 1.11 × 10−3, decreasing slightly with wind speed. DeCosmo et al. (1996) developed a bulk algorithm using the dataset of the Humidity Exchange over the Sea (HEXOS) program. The HEXOS data were collected near the Dutch coast at wind speeds of up to 18 m s−1 and are thus affected by sea spray due to shallow wave breaking. The exchange coefficients for sensible heat (CT) and latent heat (Cq) under neutral stratification determined in the study of DeCosmo et al. (1996) are CT = (1.12 ± 0.24) × 10−3 and Cq = (1.14 ± 0.35) × 10−3. In short, based on the large datasets collected under both of these programs (TOGA COARE and HEXOS), CT and Cq appear to take approximately constant values of 1.1 × 10−3 for the wind speed range of 0.5 ≤ u10 ≤ 18 m s−1, although the data exhibit considerable scatter. However, as there have been relatively few observations of latent heat flux over the ocean under stronger wind conditions, it remains unknown whether it is valid to extrapolate these bulk algorithms to stronger winds.

Emanuel (1995a) doubted that the bulk algorithms could be reasonably extended to tropical cyclone conditions, and he conducted a theoretical study of the effect of the exchange coefficient of total enthalpy on tropical cyclone intensity using a simple balance model (Emanuel 1995b) and a nonhydrostatic convection model (Rotunno and Emanuel 1987). It was thus shown that tropical cyclone intensity is dependent on the ratio of the exchange coefficient of total enthalpy (CG) to the drag coefficient (CD), leading to the conclusion that the ratio CG/CD in the strong wind region of intense storms must be greater than 0.75 to maintain the intensity of a tropical cyclone. Here, the total enthalpy G is defined as
i1520-0485-38-6-1313-e1
where T is the temperature, q is the specific humidity, cpa is the specific heat at a constant pressure of air, and LE is the latent heat of evaporation of water, respectively. The exchange coefficient of total enthalpy CG is defined as
i1520-0485-38-6-1313-e2
where HG is the total enthalpy flux, HtotS is the total sensible heat flux, HtotL is the total latent heat flux, ρa is the density of air, T10 is the temperature at a 10-m height, Tw is the sea surface temperature, q10 is the specific humidity at a 10-m height, and qw is the saturation-specific humidity at the sea surface temperature. In general, CD increases with wind speed (e.g., Large and Pond 1981). Extrapolating the results of the large programs (i.e., TOGA COARE and HEXOS) to the stronger winds (u10 > 18 m s−1), both the exchange coefficients of sensible heat and latent heat become approximately constant at 1.1 × 10−3 and are independent of wind speed. Thus, CG becomes approximately constant at 1.1 × 10−3 from Eq. (2). Then, the bulk algorithms predict that the ratio CG/CD will decrease with increasing wind speed and become less than 0.75 at u10 > 20 m s−1. Even using the recent findings that the drag coefficient does not increase with wind speed for strong winds (u10 > 33 m s−1) and holds a constant value of about 2.5 × 10−3 (Powell et al. 2003; Donelan et al. 2004), the ratio CG/CD becomes 0.44 at u10 > 33 m s−1, for which Emanuel’s theory predicts that the central surface pressure of a steady tropical cyclone will be greater than 980 hPa. In short, the exchange coefficient of total enthalpy must increase with wind speed under tropical cyclone conditions if Emanuel’s model is correct.

Sea spray has been considered a likely factor that may resolve this discrepancy (Andreas and Emanuel 2001). Under strong wind, sea spray is generated by wave breaking, and the droplets of sea spray drifting in the atmosphere continue to transfer sensible heat to the ambient atmosphere until the droplet temperature reaches the wet-bulb temperature, at which time the droplet evaporates and generates latent heat by absorbing sensible heat from the atmosphere. This heat exchange between sea spray droplets and the atmosphere mainly occurs within a droplet evaporation layer (DEL).

The possible effect of sea spray on tropical cyclone intensity has been investigated by a number of groups. For example, using a coupled ocean–atmosphere model of a tropical cyclone, Kepert et al. (1999) concluded that although sea spray had little effect on the total enthalpy flux, it could increase storm intensity by cooling and moistening the boundary layer. Wang et al. (2001) came to a similar conclusion using the parameterization of Fairall et al. (1994), but also showed that sea spray contributes directly to the total enthalpy flux if the parameterization of Andreas and DeCosmo (1999) is assumed. These early studies thus suggest that sea spray may increase the total enthalpy flux from the ocean and also cool and moisten the boundary layer.

Sensible and latent heat are transferred by turbulence in the airflow. As the airflow immediately above the ocean is disturbed by ocean waves, forming a wave boundary layer (WBL), the turbulent transfer in the WBL must be affected by ocean waves (Makin and Mastenbroek 1996; Makin 1999). A number of WBL models have been proposed to predict the drag coefficient (e.g., Janssen 1989; Makin and Kudryavtsev 1999), and most consider momentum conservation in the WBL. These studies indicate that turbulent momentum flux in the WBL is reduced by wave-induced momentum flux related to periodic airflow perturbations generated by ocean waves. Recently, Hara and Belcher (2004) proposed a WBL model based on the conservation of both momentum and kinetic energy over ocean waves, and Moon et al. (2004a) combined the model of Hara and Belcher with a wavenumber spectrum calculated by the WAVEWATCH III wave model (Tolman 2002) to investigate the momentum exchange process for young and mature seas. Application of that model under conditions of hurricane wind forcing (Moon et al. 2004b, c) revealed the same plateauing of the drag coefficient at u10 ≥ 35 m s−1 as indicated by the observations of Powell et al. (2003). This result suggests that it may soon be possible to solve many of the unanswered questions related to transfer mechanisms in the boundary layer over the ocean.

Moon et al. (2007) proposed a new roughness length of z0 for momentum, based on a sequence of their results using the WBL model (Moon et al. 2004b, c) under conditions of hurricane wind forcing, and they investigated the effects of the use of the new roughness length for momentum on the intensity of tropical cyclones. They showed that the ratio CG/CD lies in the range of 0.7–1.4 and increases as wind speed increases, by using the relation of
i1520-0485-38-6-1313-e3
i1520-0485-38-6-1313-e4
where κ is the von Kármán constant. However, because they focused on the impact of the new roughness length z0 for momentum, they simply used a Charnock-type roughness length (Charnock 1955) for temperature and humidity:
i1520-0485-38-6-1313-e5
In the case that the bulk algorithm of Zeng et al. (1998), which is based on the results of the abovementioned large programs, is adopted for the estimation of zT [shown in Eq. (49)], the ratio CG/CD using z0 of Moon et al. (2007) monotonically decreases as wind speed increases, and the ratio becomes less than 0.5 at u10 > 35 m s−1, as also shown by the traditional bulk algorithms. In their conclusion, they pointed out the need for further studies of heat transfer in strong wind conditions.

Makin and Mastenbroek (1996) and Makin (1999) evaluated heat transfer in the WBL using their WBL model.

In the present study, a boundary layer model for evaluating sensible and latent heat fluxes including sea spray–mediated heat fluxes over a mature sea is developed. The WBL model of Hara and Belcher (2004) is adopted as a basis, the wavenumber spectrum for mature sea presented by Elfouhaily et al. (1997) is used for the low-wavenumber regime, and that proposed by Hara and Belcher (2002) is used for the high-wavenumber regime. Using the information on the airflow structure obtained from the WBL model, turbulent heat transfer in the WBL is calculated using the approach of Makin and Mastenbroek (1996). The sea spray–mediated heat flux in the DEL is incorporated by introducing the parameterization of Fairall et al. (1994). The results of the developed model are used to evaluate the increase in total enthalpy under strong wind conditions.

Reduction of momentum flux in the airflow due to sea spray is not considered. The momentum of airflow is extracted by the spray droplets drifting in the atmosphere due to the acceleration of the droplets in the direction of the wind. Wu (1973) found that such momentum exchange is negligible under moderate winds, and Andreas and Emanuel (2001) showed that the contribution of this momentum exchange is three orders of magnitude less than the interfacial stress for moderate winds, although the interfacial and spray stresses become comparable for strong winds (u* ≥ 3 m s−1). The friction velocity (u* = 3 m s−1) corresponds to u10 = 60 m s−1 when CD = 0.0025. In the present study, heat transfer is investigated at wind speeds of up to u10 = 28 m s−1, at which the top of DEL reaches z = 10 m (the reference height of bulk variables). Thus, the sea spray–mediated momentum flux is considered to be negligible in the present study.

This paper is organized as follows: in section 2, the model and solving procedures are introduced. In section 3, the sensible and latent heat transfer over a mature sea without sea spray effects is evaluated, and the influence of sea spray is then discussed. The growth of the exchange coefficient of total enthalpy under tropical storm conditions is then evaluated. The main results of the present study and future work are finally presented in section 4.

2. Model description

a. Wave boundary layer model

The momentum flux over the ocean waves is estimated by applying the WBL model proposed by Hara and Belcher (2004), which is reviewed briefly here. The WBL is within a constant flux layer, and thus the total momentum flux τtot is constant with height. The downward total momentum flux τtot is composed of turbulent stress τt(z) and wave-induced stress τw(z), and thus the budget of the momentum flux is expressed as
i1520-0485-38-6-1313-e6
where u* is the friction velocity of airflow. The origin of the z axis is the sea surface. The wave-induced stress is a momentum flux due to periodic motion in the airflow caused by wind waves. The wave-induced stress at the wave surface τw(0) governs the growth and decay of wind waves and is expressed as (Makin and Mastenbroek 1996)
i1520-0485-38-6-1313-e7
where ψ(k, θ) is the wavenumber spectrum; and k, σ, θ, βg(k, θ), and ρw are the wavenumber, wave angular frequency, wave propagation direction relative to the mean wind direction, wave growth rate, and density of water, respectively.
The wave-induced stress tends to decay with height. Hara and Belcher (2004) employed a step function to represent this decay behavior:
i1520-0485-38-6-1313-e8
Here, zi(k) is the depth of the inner region in the airflow over waves of wavenumber k. Hara and Belcher (2004) set kzi(k) = δ, where δ is a constant of O(0.05) used as a fitting parameter. Using this decay function, the wave-induced stress at a height z is given by
i1520-0485-38-6-1313-e9
The local turbulent stress τt(z = δ/k) at the top of the inner region for waves of wavenumber k is given by
i1520-0485-38-6-1313-e10
where ul*[z = zi(k)] is the local friction velocity. The local turbulent stress τt(z = δ/k) forces waves of wavenumber k (Belcher and Hunt 1993); thus, the growth rate βg(k, θ) for waves of wavenumber k is estimated using the local turbulent stress. Solving the momentum flux budget above affords the vertical distribution of the local turbulent stress.
To calculate the vertical distribution of the mean wind velocity, Hara and Belcher (2004) considered the total kinetic energy conservation in the WBL, that is, the sum of the kinetic energy of the mean motions, wave-induced motion, and turbulence, as expressed by
i1520-0485-38-6-1313-e11
where ε is the viscous dissipation of turbulent kinetic energy, and Π and Π′ are the vertical transport of kinetic energy of wave-induced motion and the vertical transport of turbulent kinetic energy, respectively. It was assumed that the largest contribution to Π is from pressure transport, while the decay of Π is governed by the same decay function as the wave-induced stress:
i1520-0485-38-6-1313-e12
where g is the acceleration of gravity. It was further assumed that the divergence of turbulent transport within the WBL is smaller than the other terms in Eq. (11), because the WBL is a homogeneous rough-wall boundary layer, and can thus be neglected. The viscous dissipation ε(z) is related to the local friction velocity ul*(z) at each height and is approximated by ε = ul*(z)3/κz. Solving the kinetic energy budget allows the vertical distribution of the mean wind velocity to be estimated.

For fast-moving waves [k > km = 0.072g/(ul*)2], the growth rate βg(k, θ) is negligible. Thus, as Eq. (9) shows, the effect of wind waves does not extend above zm = δ/km, where the height zm is defined as the top of the WBL, above which τt is constant. For detailed procedures for solving the WBL model, refer to Hara and Belcher (2004) or Moon et al. (2004a).

b. Scalar transfer in WBL

Scalar transfer in the WBL is evaluated following the approach of Makin and Mastenbroek (1996). In the WBL, the total scalar flux is constant with height, as expressed by
i1520-0485-38-6-1313-e13
where HtΘ and HwΘ denote the turbulent scalar flux and the wave-induced scalar flux of a scalar Θ. It is assumed that the wave-induced scalar flux HwΘ is negligible, because it is much smaller than the turbulent scalar flux. This assumption is considered acceptable because the wave-induced scalar fluxes at the sea surface and the top of the WBL are zero and the depth of the WBL estimated by solving the WBL model of Hara and Belcher (2004) is thin, as shown by Moon et al. (2004a), who show that the depth is 1.4 m for the mature sea state under u10 = 45 m s−1 (cf. their Fig. 9).
Using the eddy viscosity assumption, the relationship between scalar flux and the vertical distribution of the concentration of Θ is determined by relating the local turbulent stress to the vertical gradient of the mean velocity:
i1520-0485-38-6-1313-e14
where Km is the eddy viscosity. Similarly, the turbulent scalar flux is expressed as
i1520-0485-38-6-1313-e15
Here, ScΘt = Km/KΘ is the turbulent Schmidt number of Θ, where KΘ is the eddy diffusivity of Θ. If Θ denotes the temperature T, then the turbulent Schmidt number denotes the turbulent Prandtl number Prt = Km/Ks, where Ks is the eddy diffusivity of sensible heat, and the upward turbulent sensible heat flux Hst is expressed as
i1520-0485-38-6-1313-e16
In the case of turbulent latent heat flux HtL,
i1520-0485-38-6-1313-e17
where Sct is the Schmidt number of water vapor. The boundary condition for q at the sea surface is qw. As the exchange bulk coefficients for sensible heat flux and latent heat flux are almost identical in many sets of observational results, Sct for water vapor can be assumed to be approximately equal to Prt. The turbulent sensible and latent heat fluxes [ HtS(0), HtL(0)] represent interfacial heat transfer.

It should be noted that although the wave-induced scalar flux HwΘ in Eq. (13) is ignored, the effects of wind waves are indirectly reflected in the scalar transfer through the use of the local turbulent stress τt in Eq. (14), which is derived from the WBL model.

c. Sea spray–mediated heat transfer

The effect of sea spray on heat transfer is evaluated by applying the sea spray–mediated heat flux parameterization proposed by Fairall et al. (1994). This parameterization allows the effects of the evaporation of oceanic spray droplets on the sensible and latent heat fluxes to be determined using bulk meteorological quantities such as wind velocity, temperature, and relative humidity (RH) at a height of 10 m and at sea surface temperature.

Fairall et al.’s (1994) parameterization is based on the droplet microphysics model of Andreas (1990, 1992). In Andreas (1990), three time scales were introduced for droplets in relation to sea spray–mediated heat transfer; the time scales of temperature evolution (τT), evaporation (τr), and atmospheric resistance for a droplet (τf). The time scale of temperature evolution determines the time necessary for a droplet to reach its equilibrium temperature, while τr represents the time necessary for a droplet to reach its equilibrium radius and is two to three orders of magnitude longer than τT. Droplets thus exchange sensible heat much more rapidly than latent heat (Andreas 1992). The time scale of atmospheric resistance for a droplet determines the duration the droplet exchanges sensible and latent heat with the ambient atmosphere. These three time scales are dependent on the initial radius of droplets. Fairall et al. (1994) integrated the heat exchanges among initial radii r0, satisfying the relation
i1520-0485-38-6-1313-e18
The sensible and latent heat exchange for droplets occurs within a DEL, the height of which can be approximated by the significant wave amplitude as (i.e., half the significant wave height hs). According to Eq. (18), the droplets transfer all excess sensible heat to the ambient atmosphere in the DEL. The sea spray–mediated sensible heat flux in the DEL is parameterized as
i1520-0485-38-6-1313-e19
where ρw and cpw denote the density and specific heat of water. Here, W(u10) is the whitecap area fraction, which is related to the sea spray–generation function dF/dr0 as
i1520-0485-38-6-1313-e20
Here, f (r0) is the source spectrum per unit area of whitecap. As the sea spray–generation function proposed by Andreas (1992) applies to an upper limit of wind velocity of 20 m s−1, application of the parameterization of Fairall et al. (1994) to strong wind conditions such as those for tropical cyclones requires a new generation function corresponding to stronger winds. Here, the generation function proposed by Andreas (1998) is employed, which can be applied to wind conditions of up to 32 m s−1. Corresponding to the new generation function, new values of W(u10) and Sυ in Eq. (19) are used because Sυ = ∫0 (4π/3)r3f (r) dr. After some arithmetic, we obtain
i1520-0485-38-6-1313-e21
Here,
i1520-0485-38-6-1313-e22
and
i1520-0485-38-6-1313-e23
where u14 is the wind velocity at a height of 14 m, which can be estimated by solving the present model.
The evaporation of droplets is assumed in Fairall et al. (1994) to occur constantly in the DEL. The sea spray–mediated latent heat flux is thus parameterized as
i1520-0485-38-6-1313-e24
where qs(Tz=as) is the saturation-specific humidity for the temperature at the height as, and β is given by
i1520-0485-38-6-1313-e25
where Rd is the gas constant for dry air. As the new generation function of Andreas (1998) is used here, the value of W(u10)Sa must be estimated. Using the value estimated based on the definition of Fairall et al. (1994) to evaluate the latent heat flux mediated by sea spray under the same meteorological conditions as considered in Andreas (1992) over mature seas, the latent heat fluxes given by the present model are one order of magnitude smaller than those reported by Andreas (1992). The value of W(u10)Sa is thus modified by multiplying by a constant, as given by
i1520-0485-38-6-1313-e26
where γ1 = 1.10 × 10−9 and γ2 = 2.05 × 10−6.
As the sensible and latent heat fluxes are constant in the constant flux layer, the total sensible and latent heat fluxes are expressed as
i1520-0485-38-6-1313-e27
i1520-0485-38-6-1313-e28
Here, the sea spray–mediated heat fluxes are assumed to be constant values estimated from Eqs. (19) and (24) in the DEL and to vanish above the DEL. This assumption means that sea spray droplets are generated at the top of the DEL and the heat exchange between the sea spray droplets and the ambient atmosphere occurs only at the top of the DEL. Although heat exchange will actually occur everywhere in the DEL, information on the heat exchange in the DEL is not sufficient to be able to determine the vertical profiles of QS(z) and QL(z). Thus, the validity of this assumption remains to be verified in future work.
The upward turbulent sensible and latent heat fluxes at a reference height higher than the DEL is expressed as
i1520-0485-38-6-1313-e29
i1520-0485-38-6-1313-e30
Because the WBL is much thinner than the DEL, the turbulent sensible and latent heat fluxes are constant above the DEL. Furthermore, because the wave-induced scalar flux is assumed to be negligible, the turbulent heat fluxes are also constant below the DEL. Thus, Eqs. (29) and (30) can be rewritten as
i1520-0485-38-6-1313-e31
i1520-0485-38-6-1313-e32
Here, HtSzas and HtLzas denote the interfacial sensible and latent heat fluxes, and HtSz>as and HtLz>as denote the total sensible and latent heat fluxes.

d. Effects of stable/unstable stratification

Using the bulk meteorological quantities at a height of 10 m and at sea surface (i.e., u10, T10, q10, Tw, qw), the fluxes are estimated by solving the WBL model and sea spray model. The profiles of the vertical gradients of wind velocity, temperature, and specific humidity [i.e., du/dz = fu(z), dT/dz = fT(z), and dq/dz = fq(z)] are also obtained. Note that the obtained values of fu(z), fT(z), and fq(z) are under a neutral condition. As bulk meteorological quantities contain stratification effects, the buoyancy effects on turbulent transfer should be taken into account. Following the Monin–Obukhov similarity, the vertical gradients are corrected in the cases of stably/unstably stratified airflow as
i1520-0485-38-6-1313-e33
i1520-0485-38-6-1313-e34
i1520-0485-38-6-1313-e35
where L is the Monin–Obukhov length defined by
i1520-0485-38-6-1313-e36
Note that the turbulent sensible and latent heat fluxes (i.e., HtS and HtL) are distributed like a step function: they are constant below and above z = as, but vary rapidly at z = as. Thus, the value of L below z = as is different from that above z = as.
The stratification corrections for the unstable condition (i.e., L < 0) are given by (Kondo 1994):
i1520-0485-38-6-1313-e37
and
i1520-0485-38-6-1313-e38
The stratification corrections for the stable condition (i.e., L > 0) are given by
i1520-0485-38-6-1313-e39

In the present study, these stratification effects are accounted for by performing iterative calculations. First, the WBL model and sea spray model are solved under a neutral condition (i.e., ϕu = ϕT = ϕq = 1), and we obtain the turbulent momentum, sensible and latent heat fluxes, and the value of fu(z), fT(z), and fq(z). Then, L and the stratification corrections (i.e., ϕu, ϕT, and ϕq) are calculated by using the turbulent fluxes. We substitute the values of L and the stratification corrections into Eqs. (33)(35), and we solve the WBL model and sea spray model again. The iteration continues until the fluxes reach appropriate convergence limits.

e. Solving procedure

Using bulk meteorological quantities and the wavenumber spectrum ψ(k, θ) for a mature wave state, the sensible and latent heat transfer under the effects of surface waves and sea spray is calculated. The bulk meteorological quantities include the wind velocity, temperature, and relative humidity at a height of 10 m and at the temperature and atmospheric pressure at the sea surface.

A wavenumber spectrum near the peak frequency proposed by Elfouhaily et al. (1997) is adopted here. This wavenumber spectrum is based on the wavenumber spectra proposed by Pierson and Moskowitz (1964) and the Joint North Sea Wave Project (JONSWAP) (Hasselmann et al. 1973) for a low-wavenumber regime, and the wavenumber spectra of Phillips (1985) and Kitaigorodskii (1973) for a high-wavenumber regime. The wavenumber spectrum ψ(k, θ) of Elfouhaily et al. have a relatively simple closed form and are consistent with past observational results. Further, ψ(k, θ) can be considered to represent a saturation value based on the wave age (cP/u10) or fetch, where cP is the wave velocity of the peak frequency. In the present study, the mature sea state is considered, for which cP/u10 is set at 1.19 (Elfouhaily et al. 1997).

In the high-wavenumber regime, which is considered to be an equilibrium range (Phillips 1977), the wavenumber spectrum of Hara and Belcher (2002) is employed. An equilibrium wave state is maintained by the balance of energy among the three sources of wave action: input from wind, dissipation due to wave breaking, and divergence of wave action flux by nonlinear wave interactions (Phillips 1985). The analytical model for ψ(k, θ) presented by Hara and Belcher (2002) is developed by introducing energy input sources obtained from the WBL model of Hara and Belcher (2004). Thus, ψ(k, θ) in the equilibrium range of Hara and Belcher (2002) is regarded as appropriate for solving the WBL model of Hara and Belcher (2004).

Hara and Belcher (2002) identify a sheltering wavenumber (ks) at which the local turbulent stress begins to be affected by sheltering by longer waves. The wavenumber spectrum in the equilibrium range [ψeq(k, θ)] is thus expressed as
i1520-0485-38-6-1313-e40
Here, h(θ) = cos2(θ) for −π/2 < θ < π/2, h(θ) = 0 for |θ| ≥ π/2, and cθ = ∫ππ [h(θ)]3/2 cosθdθ = 3π/8. The sheltering wavenumber is determined by continuously joining the omnidirectional spectrum S(k) = ∫ππ ψ(k, θ) kdθ of the equilibrium range with that of the lower wavenumber at the cutoff wavenumber (kc). Following Moon et al. (2004a), kc is set at 9.0 × kpeak, where kpeak is the wavenumber of the peak frequency. The sheltering wavenumber is then given by
i1520-0485-38-6-1313-e41

Figure 1 shows the distribution of the omnidirectional spectrum S(k) for a mature sea at several wind velocities. This omnidirectional spectrum includes the peak frequency regime of Elfouhaily et al. (1997) and the equilibrium regime of Hara and Belcher (2002). The figure shows that these two regimes are smoothly connected.

The WBL model of Hara and Belcher (2004) is then solved following Moon et al. (2004a). The upper bound of wavenumber (k1) is set at 400 rad m−1, above which Hara and Belcher (2004) showed that the calculated roughness length of momentum does not depend on the value of k1. The boundary condition at the sea surface is defined by u(z = δ/k1) = 0 at z = 0.1υa/ul*(k1). The value of δ is set at 0.008, such that the drag coefficient for the mature sea state is consistent with the parameterization of Zeng et al. (1998) for weak to moderate winds (4 ≤ u10 ≤ 18 ms−1), as shown in Fig. 2. Solving the WBL model, we can obtain the total momentum flux τtot, or friction velocity u*, and the vertical distributions of local turbulent stress τt(z) and wind velocity u(z).

Next, using the vertical distributions of τt(z) and u(z), the relationship between the turbulent heat fluxes (constant in the WBL) and the vertical distributions of temperature or specific humidity is obtained from Eqs. (14), (16), and (17). Thus, the relationship between the turbulent heat fluxes and the temperature/specific humidity at the top of the WBL (z = zm) can be expressed as
i1520-0485-38-6-1313-e42
and
i1520-0485-38-6-1313-e43
where σz=zm = ∫zmz0T Km(z)−1 dz is a value obtained by solving the WBL model. For the boundary condition at the sea surface, T(z0T) = Tw is set at z0T = 0.18υa/ul*(k1) (Liu et al. 1979), and q(z0q) = qw is set at z0q = z0T. The values of Prt and Sct are set at 1.1, such that the exchange coefficients of sensible and latent heat fluxes without sea spray effects agree with the parameterization of Zeng et al. (1998; see Fig. 2).
The sea spray–mediated heat fluxes are then determined by solving Eqs. (19), (24),
i1520-0485-38-6-1313-e44
i1520-0485-38-6-1313-e45
i1520-0485-38-6-1313-e46
and
i1520-0485-38-6-1313-e47
Equations (44)(47) are derived from Eqs. (16), (17), (31), and (32). As these six equations contain six unknown quantities (i.e., Tz=as, qz=as, HtSzas, HtLzas, QS, QL), the equations can be solved. The effects of stratification are accounted for by solving these equations iteratively.

3. Results and discussion

a. Heat transfer over a mature sea without sea spray

The sensible and latent heat fluxes affected by fully developed ocean waves without sea spray effects (i.e., for the case that heat is transferred only from the ocean interface) are compared here with the results of the bulk algorithm of Zeng et al. (1998). Their algorithm is based on data obtained by TOGA COARE under weak to moderate wind conditions (<12 m s−1) (Fairall et al. 1996), and in the HEXOS experiment under stronger wind conditions (<18 m s−1) (DeCosmo et al. 1996). The roughness length of momentum (z0), temperature (zT), and humidity (zq) are given by
i1520-0485-38-6-1313-e48
and
i1520-0485-38-6-1313-e49
where Re* = u*z0/υa is the roughness Reynolds number. Figure 2 shows the exchange coefficients (CD, CT, Cq) for mature sea states obtained in the present study along with those reported by Zeng et al. (1998). Here, the coefficients are defined as
i1520-0485-38-6-1313-e50
i1520-0485-38-6-1313-e51
and
i1520-0485-38-6-1313-e52
For this calculation, Tw = 300 K, T10 = 298 K, RH is set at 80%, and the atmospheric pressure (P0) is 1000 hPa. The values of δ and Prt are set such that the coefficients match those of Zeng et al. (1998), as described in the previous section. The exchange coefficient of latent heat flux coincides with that of sensible heat flux if sea spray effects are not incorporated, because Prt is assumed to agree with Sct in Eqs. (16) and (17). Figure 2 shows that although the dependence of CD on wind velocity is consistent among the results shown, slight differences can be observed among the results for CT. Such differences can be attributed to the wide scattering of heat flux data in TOGA COARE and HEXOS observations (see Fairall et al. 1996; DeCosmo et al. 1996). The CT and Cq parameterized from the HEXOS data of DeCosmo et al. (1996) exhibit ranges of approximately 60%, and the present results fall well within this scattering range of data. Thus, CT and CD obtained in the present study are considered acceptable for weak to moderate winds (4 ≤ u10 ≤ 18 m s−1).

b. Heat transfer over a mature sea with sea spray

To evaluate the sensible and latent heat fluxes over mature sea states considering the effects of sea spray, the depth of the DEL (as) is required. Here, as is estimated using ψ(k, θ) by (Komen et al. 1994)
i1520-0485-38-6-1313-e53

In Table 1, the present results are compared with those determined using the microphysics of Andreas (1992). The table lists the sea spray–mediated sensible and latent heat fluxes under the six conditions tabulated in Table 3 of Andreas (1992). The wind speed is fixed at u10 = 20 m s−1, while the temperatures (T10, Tw) and RH vary among cases. As can be seen in the table, the present results are generally in agreement with those of Andreas (1992), although the sensible heat flux at the lowest low in Table 1 is underestimated in the present model. This may indicate that the value of W(u10)Sυ in Eq. (19) need to be better regulated to agree with the results of Andreas (1992).

At the condition of the lowest low in Table 1, the interfacial turbulent sensible heat flux was 63.3 W m−2, which is much smaller than QL (=183.2 W m−2), and thus the total sensible heat flux is downward, even though T10 is lower than Tw. Similarly, in Andreas (1992), the corresponding interfacial turbulent sensible heat flux is 46 W m−2 using a bulk algorithm (Large and Pond 1981; Andreas and Murphy 1986), and it is also much smaller than QL (=170 W m−2). In this case, the sea spray from the warmer ocean makes the cool atmosphere cooler, and this cooling and moistening effect would influence the tropical cyclone intensity (e.g., Kepert et al. 1999; Lighthill et al. 1994).

Figure 3 shows the effects of sea spray and surface waves on the vertical profiles of temperature and specific humidity under moderate (u10 = 7 m s−1) and strong (u10 = 25 m s−1) winds for a stably stratified condition with Tw = 300 K, T10 = 302 K, P0 = 1000 hPa, and RH set at 80%. The calculated heat fluxes for the moderate wind case are QS = −0.08 W m−2, QL = 0.35 W m−2, HtotS = −12.28 W m−2, and HtotL = 27.55 W m−2, while those for the strong wind case are QS = −14.4 W m−2, QL = 855.4 W m−2, HtotS = −921.2 W m−2, and HtotL = 1006.8 W m−2. For the moderate wind conditions, the magnitudes of sea spray–mediated heat fluxes are much smaller than those of the total heat fluxes, and the sea spray would have little effect on the heat transfer mechanism. Under strong wind conditions, on the other hand, the sea spray–mediated latent heat flux is comparable in magnitude to the total heat fluxes; thus, the sea spray effect would be strong.

Figure 3 shows that the WBL is much thinner than the DEL, but becomes thicker with increasing wind speed. Although the temperature and specific humidity are approximately distributed logarithmically above the WBL, the distributions in the WBL are not logarithmic. The vertical gradient of temperature just below the top of the WBL (at approximately z = 0.1 m at 25 m s−1) is larger than that extrapolated from the logarithmic layer (at approximately z = 1 m at 25 m s−1), attributable to the wave-induced kinetic energy in Eq. (11), which decreases with height [i.e., dΠ/dz in Eq. (11) is negative]. Thus, the vertical gradient of velocity and eddy viscosity in the WBL become larger than those above the WBL, resulting in a corresponding increase in the modulus of the vertical gradient of temperature and humidity. If the vertical transport of kinetic energy of the wave-induced motion is not considered in Eq. (11), the vertical gradient of temperature just below the WBL becomes smaller than that in the logarithmic layer, because the modulus of the vertical gradient of temperature in this case is expressed as
i1520-0485-38-6-1313-e54

For moderate winds, the contributions of sea spray–mediated fluxes are much smaller than those of turbulent fluxes at the interface, and thus almost no changes are observed in the vertical gradients of T and q at the top of the DEL (i.e., at z = as). However, because of the effect of stable stratification, the vertical gradients are misaligned with a logarithmic low near z ≅ 10 m. For strong winds, on the other hand, the vertical gradient of T and q changes quickly at z = as. Strong winds generate substantial volumes of sea spray, and the evaporation of droplets cool the atmosphere in the DEL. The temperature below the DEL thus becomes lower than that without sea spray effects.

Figure 4 shows the vertical profiles of temperature and specific humidity for an unstably stratified condition with Tw = 300 K, T10 = 298 K, P0 = 1000 hPa, and RH set at 80%. The values obtained for moderate wind are QS = 0.09 W m−2, QL = 0.35 W m−2, HtotS = 14.83 W m−2, and HtotL = 113.81 W m−2, and those for strong wind are QS = 27.5 W m−2, QL = 821.7 W m−2, HtotS = −693.8 W m−2, and HtotL = 1398.6 W m−2. As is the case with the stably stratified condition, the magnitudes of sea spray–mediated heat fluxes are much smaller than those of the total heat fluxes under moderate wind conditions, whereas the magnitude of sea spray–mediated latent heat flux is comparable to that of the total heat fluxes under strong wind conditions. Thus, under moderate wind conditions, the contributions of sea spray–mediated fluxes are much smaller than those of turbulent fluxes at the interface, and almost no changes are observed in the vertical gradients of T and q at the top of the DEL. For strong wind, the sign of the vertical gradient of T changes at z = as, because the turbulent sensible heat flux above the DEL is negative while that below the DEL is positive. In this case, the temperature at z = as is lower than T10, resulting in larger sea spray–mediated sensible and latent heat fluxes than that estimated by Fairall’s original parameterization, in which the temperature and specific humidity at z = as in Eqs. (19) and (24) are approximated from T10 by
i1520-0485-38-6-1313-e55
and
i1520-0485-38-6-1313-e56
The value of (TwTz=as) in Eq. (19) is 1.96 K in Fairall’s original parameterization for u10 = 25 m s−1, while the value given by the present model is 2.46 K, resulting in larger QS and QL. Andreas and DeCosmo (2002) incorporated this cooling effect by reducing the feedback of QL to Hstot, effectively adjusting the heat fluxes by tuning the respective contributions.

The cooling effect also increases the interfacial sensible heat fluxes (i.e., HtSz<as), because the turbulent sensible heat flux depends on the vertical gradient of temperature. This is an indirect contribution of sea spray evaporation to the total heat fluxes and would facilitate the development of tropical cyclones.

It should be noted that under such unstably stratified and strong wind conditions, the total enthalpy flux is upward but the total sensible heat flux is downward. The total enthalpy flux HG is given by
i1520-0485-38-6-1313-e57
Emanuel (1995a) argued that the development of tropical cyclones is dependent on the total enthalpy flux. The exchange coefficients of total enthalpy (CG) under stably and unstably stratified conditions are plotted in Fig. 5 for the wind speed range of 4 ≤ u10 ≤ 28 m s−1. Although the spray-generation function adopted here is applicable at wind speeds of up to u10 = 32 m s−1, the top of the DEL under stronger winds (u10 > 28 m s−1) exceeds 10 m (not shown). Under the stably stratified condition, Tw = 300 K, T10 = 302 K, P0 = 1000 hPa, and RH is set at 80%. Under the unstably stratified condition, Tw = 300 K, T10 = 298 K, P0 = 1000 hPa, and RH is set at 80%.

For weak to moderate winds (4 ≤ u10 ≤ 12 m s−1), the buoyancy effect on turbulence causes a difference in CG among the two cases. In the stably stratified airflow, the turbulence becomes weak, causing CG to become smaller than that in the neutral case, while under unstably stratified conditions, turbulence increases, resulting in higher CG. In the case of strong winds (u10 > 20 m s−1), the sea spray–mediated heat fluxes are responsible for the difference between the two cases. Table 2 lists the interfacial and sea spray–mediated heat fluxes for both cases under the strongest wind condition (u10 = 28 m s−1), as well as the interfacial heat fluxes without sea spray effects. Under the stably stratified condition, the total enthalpy flux becomes positive, even though the interfacial and sea spray–mediated sensible heat fluxes are negative. Although the sea spray indirectly increases the interfacial heat fluxes by 7 W m−2 because of cooling in the DEL, the sea spray–mediated heat flux directly reduces the total enthalpy flux by 36 W m−2, resulting in an overall decrease in total enthalpy flux. Under the unstably stratified condition, on the other hand, the large sea spray–mediated latent heat flux cools the DEL, and the interfacial heat fluxes are indirectly increased by 44 W m−2. In conjunction with the direct increase in heat flux of 49 W m−2, the combined effects increase the total enthalpy flux by approximately 10%. The numerical results of Wang et al. (2001) indicate that the parameterization of Fairall et al. (1994) results in a smaller direct contribution of sea spray to the total enthalpy flux than its indirect contribution in the tropical cyclone model. This is considered to be caused by the omission of the cooling effect below the reference height of bulk variables (i.e., z = 10 m), resulting in an underestimation of sea spray–mediated heat flux. Thus, the present model would give a larger total contribution by sea spray evaporation in the tropical cyclone model.

In the unstably stratified case, the increasing rate of CG with wind speed under strong wind (u10 ≥ 20 m s−1) is consistent with that for CD. In other words, the dependence of CG on wind speed is comparable to that for CD under strong wind conditions. Figure 6 shows the ratio of CG to CD under the unstably stratified condition. The parameterization of Zeng et al. (1998) for a weak to moderate wind regime (4 ≤ u10 ≤ 18 m s−1) and the extrapolation to a strong wind regime (18 ≤ u10 ≤ 28 m s−1) are also plotted. While CG/CD given by the parameterization of Zeng et al. decreases with increasing wind speed, CG/CD in the present study reaches a plateau at 0.7 for strong winds (u10 ≥ 20 m s−1). This result supports the theory of Emanuel (1995a); that is, the ratio CG/CD should be greater than 0.75 to maintain the intensity of a tropical cyclone.

The slightly lower value CG/CD = 0.7 (<0.75) in the present study may be due to an underestimation of the sea spray–mediated sensible heat flux QS. The lowest low of Table 1, which is a similar condition to Fig. 6, shows that the sea spray–mediated sensible heat flux QS in the present study is approximately half that reported by Andreas (1992). If the value of W(u10)Sυ in Eq. (19) is multiplied by a constant value of 2.2 to agree with our QS of the lowest low of Table 1 with that of Andreas (1992), the ratio increases with wind speed, reaching 0.75 at u10 = 27 m s−1 (see Fig. 6). Note that although this regulation makes the prediction of QS under the condition of Fig. 6 better, it is not appropriate for all the conditions. Thus, we have not used the regulated QS in the previous discussions.

Figure 7 shows the variation in CG/CD under various conditions. These are calculated using the value of QS without the abovementioned regulation. The temperature difference TwT10 in these examples varies from 0.5 to 3.0 K, and RH varies from 70% to 90%. Under the strongest wind conditions (u10 = 28 m s−1), CG/CD increases with increasing temperature difference or decreasing relative humidity. It is interesting to note that in every case, the ratio levels off at approximately 0.7 with winds exceeding 20 m s−1.

4. Conclusions

A boundary layer model over a mature sea was developed for the evaluation of sensible and latent heat fluxes at wind speeds of up to 28 m s−1. To estimate the turbulent stress, we solved the wave boundary layer model of Hara and Belcher (2004). Sea spray effects were incorporated using the parameterization of Fairall et al. (1994). It was confirmed that for weak to moderate winds (4 ≤ u10 ≤ 18 m s−1), the drag coefficient and the heat exchange coefficient without sea spray effects are consistent with a bulk algorithm developed using TOGA COARE and HEXOS observations.

The main result of the present study is that while the ratio of the exchange coefficient of total enthalpy to the drag coefficient (CG/CD) decreases with increasing wind speed below 20 m s−1, the ratio plateaus at approximately 0.7 for stronger winds (20 ≤ u10 ≤ 28 m s−1) under tropical cyclone conditions (wet air blowing over a warm mature sea). It has not been possible to obtain this result using traditional bulk algorithms, which afford a monotonic decrease in CG/CD with increasing wind speed. Our result supports the argument of Emanuel (1995a), who showed that the intensity of a tropical cyclone can only be maintained if this ratio remains greater than 0.75 in the strong wind region of intense storms.

The wave boundary layer model adopted here assumes that the perturbations due to interfacial waves reach the inner region depth, which is expressed simply as zi = δ/k. This expression is applicable for low or high wave ages, but it is not appropriate for intermediate wave ages, at which the inner region depth is thicker and dependent on wave age. The characteristics of flow perturbations differ among these cases, as shown in the previous study of Kihara et al. (2007). Although the wave-age range of intermediate wave age is not large and its effects on the momentum and heat transfer mechanism would be weak, we should confirm the effects of intermediate wave age.

Under tropical cyclones, the local wave state is determined not only by the local wind speed but also by the past wave state; that is, it is not a mature sea state, although the proposed boundary layer model in the present study is developed using wave spectrum under mature sea states. In applying our proposed model to the tropical cyclone conditions, the local wave spectrum should be used instead of the wave spectrum under a mature sea. To incorporate the local wave spectrum, it is necessary to couple our model with a wave model, such as WAVEWATCH III or WAM cycle 4, as carried out in Moon et al. (2004a, b,c).

Acknowledgments

I thank two anonymous referees for their useful comments and suggestions.

REFERENCES

  • Andreas, E. L. 1990. Time constants for the evolution of sea spray droplets. Tellus 42B:481497.

  • Andreas, E. L. 1992. Sea spray and the turbulent air-sea heat fluxes. J. Geophys. Res. 97:1142911441.

  • Andreas, E. L. 1998. A new sea spray generation function for wind speeds up to 32 m s−1. J. Phys. Oceanogr. 28:21752184.

  • Andreas, E. L. and B. Murphy. 1986. Bulk transfer coefficients for heat and momentum over leads and polynyas. J. Phys. Oceanogr. 16:18751883.

    • Search Google Scholar
    • Export Citation
  • Andreas, E. L. and J. K. DeCosmo. 1999. Sea spray production and influence on air-sea heat and moisture fluxes over the open ocean. Air-Sea Exchange: Physics, Chemistry and Dynamics, G. L. Geernaert, Ed., Kluwer Academic, 363–409.

    • Search Google Scholar
    • Export Citation
  • Andreas, E. L. and K. A. Emanuel. 2001. Effects of sea spray on tropical cyclone intensity. J. Atmos. Sci. 58:37413751.

  • Andreas, E. L. and J. K. DeCosmo. 2002. The signature of sea spray in the HEXOS turbulent heat flux data. Bound.-Layer Meteor. 103:303333.

    • Search Google Scholar
    • Export Citation
  • Belcher, S. E. and J. C. R. Hunt. 1993. Turbulent shear flow over slowly moving waves. J. Fluid Mech. 251:109148.

  • Charnock, H. 1955. Wind stress on a water surface. Quart. J. Roy. Meteor. Soc. 81:639640.

  • 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:1200112016.

    • Search Google Scholar
    • Export Citation
  • 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.

    • Search Google Scholar
    • Export Citation
  • Elfouhaily, T., B. Chapron, K. Katsaros, and D. Vandemark. 1997. A unified directional spectrum for long and short wind-driven waves. J. Geophys. Res. 102:1578115796.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1986. An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. J. Atmos. Sci. 43:585604.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1995a. Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics. J. Atmos. Sci. 52:39693976.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1995b. The behavior of a simple hurricane model using a convective scheme based on subcloud-layer entropy equilibrium. J. Atmos. Sci. 52:39603968.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., J. D. Kepert, and G. J. Holland. 1994. The effect of sea spray on surface energy transports over the ocean. Global Atmos. Ocean Syst. 2:121142.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., E. F. Bradley, D. P. Rogers, J. B. Edson, and G. S. Young. 1996. Bulk parameterization of air-sea fluxes for Tropical Ocean-Global Atmosphere Coupled-Ocean Atmosphere Response Experiment. J. Geophys. Res. 101:37473764.

    • Search Google Scholar
    • Export Citation
  • Hara, T. and S. E. Belcher. 2002. Wind forcing in the equilibrium range of wind-wave spectra. J. Fluid Mech. 470:223245.

  • Hara, T. and S. E. Belcher. 2004. Wind profile and drag coefficient over mature ocean surface wave spectra. J. Phys. Oceanogr. 34:23452358.

    • Search Google Scholar
    • Export Citation
  • Coauthors 1973. Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Dtsch. Hydrogh. Z. Suppl. A 8:195.

    • Search Google Scholar
    • Export Citation
  • Janssen, P. S. E. M. 1989. Wave-induced stress and the drag of air flow over sea waves. J. Phys. Oceanogr. 19:745754.

  • Kepert, J. D., C. W. Fairall, and J-W. Bao. 1999. Modelling the interaction between the atmospheric boundary layer and evaporating sea spray droplets. Air-Sea Exchange: Physics, Chemistry and Dynamics, G. L. Geernaert, Ed., Kluwer Academic, 363–409.

    • Search Google Scholar
    • Export Citation
  • Kihara, N., H. Hanazaki, T. Mizuya, and H. Ueda. 2007. Relationship between the airflow at the critical height and momentum transfer to the traveling waves. Phys. Fluids 19.015 102, doi:10.1063/1.2409736.

    • Search Google Scholar
    • Export Citation
  • Kitaigorodskii, S. A. 1973. The Physics of Air-Sea Interactions. Keter Press, 273 pp.

  • Komen, G. J., L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, and P. A. E. M. Janssen. 1994. Dynamics and Modeling of Ocean Waves. Cambridge University Press, 520 pp.

    • Search Google Scholar
    • Export Citation
  • Kondo, J. 1994. Meteorology of the Water Environment: Water and Heat Balance of the Earth’s Surface. (in Japanese). Asakura Shoten Press, 348 pp.

    • Search Google Scholar
    • Export Citation
  • Large, W. G. and S. Pond. 1981. Open ocean momentum flux measurements in moderate to strong winds. J. Phys. Oceanogr. 11:324336.

  • Lighthill, J., G. Holland, W. Gray, C. Landsea, G. Craig, J. Evans, Y. Kurihara, and C. Guard. 1994. Global climate change and tropical cyclones. Bull. Amer. Meteor. Soc. 75:21472157.

    • Search Google Scholar
    • Export Citation
  • 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 constrains at the interface. J. Atmos. Sci. 36:17221735.

    • Search Google Scholar
    • Export Citation
  • Makin, V. K. 1999. A note on wind speed and sea state dependence of the heat exchange coefficient. Bound.-Layer Meteor. 91:127134.

  • Makin, V. K. and C. Mastenbroek. 1996. Impact of waves on air-sea exchange of sensible heat and momentum. Bound.-Layer Meteor. 79:301306.

    • Search Google Scholar
    • Export Citation
  • Makin, V. K. and V. N. Kudryavtsev. 1999. Coupled sea surface-atmosphere model 1. Wind over waves coupling. J. Geophys. Res. 104:76137623.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., T. Hara, I. Ginis, S. E. Belcher, and H. Tolman. 2004a. Effect of surface waves on air–sea momentum exchange. Part I: Effect of mature and growing seas. J. Atmos. Sci. 61:23212333.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, and T. Hara. 2004b. Effect of surface waves on air–sea momentum exchange. Part II: Behavior of drag coefficient under tropical cyclones. J. Atmos. Sci. 61:23342348.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, and T. Hara. 2004c. Effect of surface waves on Charnock coefficient under tropical cyclones. Geophys. Res. Lett. 31.L20302, doi:10.1029/2004GL020988.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, T. Hara, and B. Thomas. 2007. Physics-based parameterization of air–sea momentum flux at high wind speeds and its impact on hurricane intensity predictions. Mon. Wea. Rev. 135:28692878.

    • Search Google Scholar
    • Export Citation
  • Phillips, O. M. 1977. The Dynamics of the Upper Ocean. 2nd ed. Cambridge University Press, 336 pp.

  • Phillips, O. M. 1985. Spectral and statistical properties of the equilibrium range in the wind-generated gravity waves. J. Fluid Mech. 156:505531.

    • Search Google Scholar
    • Export Citation
  • Pierson, W. J. and L. Moskowitz. 1964. A proposed spectral form for fully developed wind sea based on the similarity theory of S. A. Kitaigorodskii. J. Geophys. Res. 69:51815190.

    • Search Google Scholar
    • Export Citation
  • Powell, M. D., P. J. Vickery, and T. A. Reinhold. 2003. Reduced drag coefficient for high wind speeds in tropical cyclones. Nature 422:279283.

    • Search Google Scholar
    • Export Citation
  • Rotunno, R. and K. A. Emanuel. 1987. An air–sea interaction theory for tropical cyclones. Part II. J. Atmos. Sci. 44:542561.

  • Tolman, H. L. 2002. User manual and system documentation of WAVEWATCH-III, version 2.22. NOAA/NWS/NCEP/OMB Tech. Note 222, 133 pp.

  • Wang, Y., J. D. Kepert, and J. Holland. 2001. The effect of sea spray evaporation on tropical cyclone boundary layer structure and intensity. Mon. Wea. Rev. 129:24812500.

    • Search Google Scholar
    • Export Citation
  • Wu, J. 1973. Spray in the atmospheric surface layer: Laboratory study. J. Geophys. Res. 78:511519.

  • Zeng, X., M. Zhao, and R. E. Dickinson. 1998. Intercomparizon of bulk aerodynamic algorithms for the computation of sea surface fluxes using TOGA COARE and TAO data. J. Climate 11:26282644.

    • Search Google Scholar
    • Export Citation

Fig. 1.
Fig. 1.

Distributions of the omnidirectional spectrum S(k) for a mature sea for u10 = 5, 10, 15, 20, and 25 m s−1.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 2.
Fig. 2.

Wind speed dependence of exchange coefficients CD, CT, and Cq for mature sea states obtained in the present study. Those reported by Zeng et al. (1998) are also shown. The solid line denotes CT and Cq in the present study; the dashed line is CD in the present study; the chain line is CT and Cq of Zeng et al.; and the chain double-dashed line is CD of Zeng et al.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 3.
Fig. 3.

(a) Vertical profiles of temperature and (b) specific humidity for u10 = 7 and 25 m s−1 under a stably stratified condition, with Tw = 300 K, t10 = 30 2K, RH = 80%, and P0 = 1000 hPa. A black circle denotes the height of the top of the wave boundary layer; an open circle denotes the height of the top of the drop evaporation layer.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 4.
Fig. 4.

(a) Vertical profiles of temperature and (b) specific humidity for u10 = 7 and 25 m s−1 under an unstably stratified condition, with Tw = 300 K, T10 = 298 K, RH = 80%, and P0 = 1000 hPa. A black dot denotes the height of the top of the wave boundary layer; an open circle denotes the height of the top of the drop evaporation layer.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 5.
Fig. 5.

Wind speed dependence of the exchange coefficient of total enthalpy CG and the drag coefficient CD under stably and unstably stratified conditions. The solid line denotes CG under the unstably stratified condition; the dashed line is CG under the stably stratified condition; the chain line is CD under the unstably stratified condition; and the chain double-dashed line is CD under the stably stratified condition.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 6.
Fig. 6.

Wind velocity dependence of CG/CD under an unstably stratified condition. The solid line denotes our result; the dashed line is the parameterization of Zeng et al. (1998); the dotted line is the extrapolation of the parameterization of Zeng et al.; and the chain line is our regulated result.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Fig. 7.
Fig. 7.

Wind velocity dependence of CG/CD under unstably stratified conditions. The temperature difference TwT10 varies from 0.5 to 3.0 K, and the relative humidity varies from 0.7 to 0.9.

Citation: Journal of Physical Oceanography 38, 6; 10.1175/2007JPO3828.1

Table 1.

Comparison of sea spray–mediated sensible and latent heat fluxes obtained in the present study (fourth and fifth columns) with those of Andreas (1992) (sixth and seventh columns). For every case, the wind velocity is 20 m s−1.

Table 1.
Table 2.

The interfacial and sea spray–mediated heat fluxes under both stably and unstably stratified cases for the strongest wind condition (u10 = 28 m s−1). The interfacial heat fluxes without sea spray effects are also shown. Under the stably stratified condition, we set Tw = 300 K, T10 = 302 K, RH = 80%, and P0 = 1000 hPa. On the other hand, under the unstably stratified condition, we set Tw = 300 K, T10 = 298 K, RH = 80%, and P0 = 1000 hPa.

Table 2.
Save
  • Andreas, E. L. 1990. Time constants for the evolution of sea spray droplets. Tellus 42B:481497.

  • Andreas, E. L. 1992. Sea spray and the turbulent air-sea heat fluxes. J. Geophys. Res. 97:1142911441.

  • Andreas, E. L. 1998. A new sea spray generation function for wind speeds up to 32 m s−1. J. Phys. Oceanogr. 28:21752184.

  • Andreas, E. L. and B. Murphy. 1986. Bulk transfer coefficients for heat and momentum over leads and polynyas. J. Phys. Oceanogr. 16:18751883.

    • Search Google Scholar
    • Export Citation
  • Andreas, E. L. and J. K. DeCosmo. 1999. Sea spray production and influence on air-sea heat and moisture fluxes over the open ocean. Air-Sea Exchange: Physics, Chemistry and Dynamics, G. L. Geernaert, Ed., Kluwer Academic, 363–409.

    • Search Google Scholar
    • Export Citation
  • Andreas, E. L. and K. A. Emanuel. 2001. Effects of sea spray on tropical cyclone intensity. J. Atmos. Sci. 58:37413751.

  • Andreas, E. L. and J. K. DeCosmo. 2002. The signature of sea spray in the HEXOS turbulent heat flux data. Bound.-Layer Meteor. 103:303333.

    • Search Google Scholar
    • Export Citation
  • Belcher, S. E. and J. C. R. Hunt. 1993. Turbulent shear flow over slowly moving waves. J. Fluid Mech. 251:109148.

  • Charnock, H. 1955. Wind stress on a water surface. Quart. J. Roy. Meteor. Soc. 81:639640.

  • 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:1200112016.

    • Search Google Scholar
    • Export Citation
  • 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.

    • Search Google Scholar
    • Export Citation
  • Elfouhaily, T., B. Chapron, K. Katsaros, and D. Vandemark. 1997. A unified directional spectrum for long and short wind-driven waves. J. Geophys. Res. 102:1578115796.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1986. An air-sea interaction theory for tropical cyclones. Part I: Steady-state maintenance. J. Atmos. Sci. 43:585604.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1995a. Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics. J. Atmos. Sci. 52:39693976.

    • Search Google Scholar
    • Export Citation
  • Emanuel, K. A. 1995b. The behavior of a simple hurricane model using a convective scheme based on subcloud-layer entropy equilibrium. J. Atmos. Sci. 52:39603968.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., J. D. Kepert, and G. J. Holland. 1994. The effect of sea spray on surface energy transports over the ocean. Global Atmos. Ocean Syst. 2:121142.

    • Search Google Scholar
    • Export Citation
  • Fairall, C. W., E. F. Bradley, D. P. Rogers, J. B. Edson, and G. S. Young. 1996. Bulk parameterization of air-sea fluxes for Tropical Ocean-Global Atmosphere Coupled-Ocean Atmosphere Response Experiment. J. Geophys. Res. 101:37473764.

    • Search Google Scholar
    • Export Citation
  • Hara, T. and S. E. Belcher. 2002. Wind forcing in the equilibrium range of wind-wave spectra. J. Fluid Mech. 470:223245.

  • Hara, T. and S. E. Belcher. 2004. Wind profile and drag coefficient over mature ocean surface wave spectra. J. Phys. Oceanogr. 34:23452358.

    • Search Google Scholar
    • Export Citation
  • Coauthors 1973. Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Dtsch. Hydrogh. Z. Suppl. A 8:195.

    • Search Google Scholar
    • Export Citation
  • Janssen, P. S. E. M. 1989. Wave-induced stress and the drag of air flow over sea waves. J. Phys. Oceanogr. 19:745754.

  • Kepert, J. D., C. W. Fairall, and J-W. Bao. 1999. Modelling the interaction between the atmospheric boundary layer and evaporating sea spray droplets. Air-Sea Exchange: Physics, Chemistry and Dynamics, G. L. Geernaert, Ed., Kluwer Academic, 363–409.

    • Search Google Scholar
    • Export Citation
  • Kihara, N., H. Hanazaki, T. Mizuya, and H. Ueda. 2007. Relationship between the airflow at the critical height and momentum transfer to the traveling waves. Phys. Fluids 19.015 102, doi:10.1063/1.2409736.

    • Search Google Scholar
    • Export Citation
  • Kitaigorodskii, S. A. 1973. The Physics of Air-Sea Interactions. Keter Press, 273 pp.

  • Komen, G. J., L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, and P. A. E. M. Janssen. 1994. Dynamics and Modeling of Ocean Waves. Cambridge University Press, 520 pp.

    • Search Google Scholar
    • Export Citation
  • Kondo, J. 1994. Meteorology of the Water Environment: Water and Heat Balance of the Earth’s Surface. (in Japanese). Asakura Shoten Press, 348 pp.

    • Search Google Scholar
    • Export Citation
  • Large, W. G. and S. Pond. 1981. Open ocean momentum flux measurements in moderate to strong winds. J. Phys. Oceanogr. 11:324336.

  • Lighthill, J., G. Holland, W. Gray, C. Landsea, G. Craig, J. Evans, Y. Kurihara, and C. Guard. 1994. Global climate change and tropical cyclones. Bull. Amer. Meteor. Soc. 75:21472157.

    • Search Google Scholar
    • Export Citation
  • 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 constrains at the interface. J. Atmos. Sci. 36:17221735.

    • Search Google Scholar
    • Export Citation
  • Makin, V. K. 1999. A note on wind speed and sea state dependence of the heat exchange coefficient. Bound.-Layer Meteor. 91:127134.

  • Makin, V. K. and C. Mastenbroek. 1996. Impact of waves on air-sea exchange of sensible heat and momentum. Bound.-Layer Meteor. 79:301306.

    • Search Google Scholar
    • Export Citation
  • Makin, V. K. and V. N. Kudryavtsev. 1999. Coupled sea surface-atmosphere model 1. Wind over waves coupling. J. Geophys. Res. 104:76137623.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., T. Hara, I. Ginis, S. E. Belcher, and H. Tolman. 2004a. Effect of surface waves on air–sea momentum exchange. Part I: Effect of mature and growing seas. J. Atmos. Sci. 61:23212333.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, and T. Hara. 2004b. Effect of surface waves on air–sea momentum exchange. Part II: Behavior of drag coefficient under tropical cyclones. J. Atmos. Sci. 61:23342348.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, and T. Hara. 2004c. Effect of surface waves on Charnock coefficient under tropical cyclones. Geophys. Res. Lett. 31.L20302, doi:10.1029/2004GL020988.

    • Search Google Scholar
    • Export Citation
  • Moon, I-J., I. Ginis, T. Hara, and B. Thomas. 2007. Physics-based parameterization of air–sea momentum flux at high wind speeds and its impact on hurricane intensity predictions. Mon. Wea. Rev. 135:28692878.

    • Search Google Scholar
    • Export Citation
  • Phillips, O. M. 1977. The Dynamics of the Upper Ocean. 2nd ed. Cambridge University Press, 336 pp.

  • Phillips, O. M. 1985. Spectral and statistical properties of the equilibrium range in the wind-generated gravity waves. J. Fluid Mech. 156:505531.

    • Search Google Scholar
    • Export Citation
  • Pierson, W. J. and L. Moskowitz. 1964. A proposed spectral form for fully developed wind sea based on the similarity theory of S. A. Kitaigorodskii. J. Geophys. Res. 69:51815190.

    • Search Google Scholar
    • Export Citation
  • Powell, M. D., P. J. Vickery, and T. A. Reinhold. 2003. Reduced drag coefficient for high wind speeds in tropical cyclones. Nature 422:279283.

    • Search Google Scholar
    • Export Citation
  • Rotunno, R. and K. A. Emanuel. 1987. An air–sea interaction theory for tropical cyclones. Part II. J. Atmos. Sci. 44:542561.

  • Tolman, H. L. 2002. User manual and system documentation of WAVEWATCH-III, version 2.22. NOAA/NWS/NCEP/OMB Tech. Note 222, 133 pp.

  • Wang, Y., J. D. Kepert, and J. Holland. 2001. The effect of sea spray evaporation on tropical cyclone boundary layer structure and intensity. Mon. Wea. Rev. 129:24812500.

    • Search Google Scholar
    • Export Citation
  • Wu, J. 1973. Spray in the atmospheric surface layer: Laboratory study. J. Geophys. Res. 78:511519.

  • Zeng, X., M. Zhao, and R. E. Dickinson. 1998. Intercomparizon of bulk aerodynamic algorithms for the computation of sea surface fluxes using TOGA COARE and TAO data. J. Climate 11:26282644.

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

    Distributions of the omnidirectional spectrum S(k) for a mature sea for u10 = 5, 10, 15, 20, and 25 m s−1.

  • Fig. 2.

    Wind speed dependence of exchange coefficients CD, CT, and Cq for mature sea states obtained in the present study. Those reported by Zeng et al. (1998) are also shown. The solid line denotes CT and Cq in the present study; the dashed line is CD in the present study; the chain line is CT and Cq of Zeng et al.; and the chain double-dashed line is CD of Zeng et al.

  • Fig. 3.

    (a) Vertical profiles of temperature and (b) specific humidity for u10 = 7 and 25 m s−1 under a stably stratified condition, with Tw = 300 K, t10 = 30 2K, RH = 80%, and P0 = 1000 hPa. A black circle denotes the height of the top of the wave boundary layer; an open circle denotes the height of the top of the drop evaporation layer.

  • Fig. 4.

    (a) Vertical profiles of temperature and (b) specific humidity for u10 = 7 and 25 m s−1 under an unstably stratified condition, with Tw = 300 K, T10 = 298 K, RH = 80%, and P0 = 1000 hPa. A black dot denotes the height of the top of the wave boundary layer; an open circle denotes the height of the top of the drop evaporation layer.

  • Fig. 5.

    Wind speed dependence of the exchange coefficient of total enthalpy CG and the drag coefficient CD under stably and unstably stratified conditions. The solid line denotes CG under the unstably stratified condition; the dashed line is CG under the stably stratified condition; the chain line is CD under the unstably stratified condition; and the chain double-dashed line is CD under the stably stratified condition.

  • Fig. 6.

    Wind velocity dependence of CG/CD under an unstably stratified condition. The solid line denotes our result; the dashed line is the parameterization of Zeng et al. (1998); the dotted line is the extrapolation of the parameterization of Zeng et al.; and the chain line is our regulated result.

  • Fig. 7.

    Wind velocity dependence of CG/CD under unstably stratified conditions. The temperature difference TwT10 varies from 0.5 to 3.0 K, and the relative humidity varies from 0.7 to 0.9.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 116 38 16
PDF Downloads 66 16 5