Estimating Air-Sea Energy and Momentum Exchanges inside Tropical Cyclones Using the Fetch- and Duration-Limited Wave Growth Properties

Transcription

1 Estimating Air-Sea Energy and Momentum Exchanges inside Tropical Cyclones Using the Fetch- and Duration-Limited Wave Growth Properties Paul A. Hwang 1 and Yalin Fan 1 Remote Sensing Division Oceanography Division U. S. Naval Research Laboratory

2 Fetch- or duration-limited wave growth Hs g 7 xg f U U U xg f U Dimensionless Fetch Graph Wave speed c/u Wave height gh/u g rms #, 4 # U 3.4 H s g 3.94 U 4 16U Fetch limited U g p H sg 8 tdg U U U td g.94 U Hwang and Wang 004 JPO Duration limited Sverdrup & Munk 1947 Fetch gx/u Wave age similarity Dimensionless Duration Graph Duration to fetch Wave speed c/u Hurricanes, buoy obs. (Young 1998, 006) Wave height gh/u Duration gt/u Hurricanes (interior obs.) Hwang and Walsh (016)

3 3 Fetch- or duration-limited growth Hs g 7 xg f U U U xg f U H s g 3.94 U 4 16U Hwang and Wang 004 JPO H sg 8 tdg U U U td g.94 U Wind-wave triplets (U, H s, T p ) H 8. U x s x T 9.8 U x p x H 1.55 U t t T 3.53 U t p s t Hwang 016 JPO U H x U 91.49T x s x p x U 39.95H t U 17.71T t s t p t

4 d de E gq d f U, H, T E Q s dt Q Q Q in nl dis t w in E s p dt t EaU ; E 0.0# # H s g Hwang and Sletten 008 JGR E 5.88 / * in # C U C U C U d a / a / a C / 1.89 # in C 1. u / g Hwang 004, 005 JO p Energy exchange U U # # E ()= M ()c(); M t =f M (U,H s,t p ) M M U ; t M a M # # # Hwang and Walsh 016 JPO (in press) U p / Momentum exchange g

5 Fetch- or duration-limited growth Hs g 7 xg f U U U xg f U H s g 3.94 U 4 16U Hwang and Wang 004 JPO H sg 8 tdg U U U td g.94 U Wind-wave triplets (U, H s, T p ) H 8. U x s U H x U 91.49T x x T 9.8 U x p x s x p x H 1.55 U t U 39.95H t U 17.71T t t T 3.53 U t p s t s t p t Hurricane fetch and duration s x 4.4 U H x p x.9 U T x Hwang 016 JPO s t 1.75 U H t p t 4.81 U T t 5

6 x r, s r I x x x x r, s r I x x x t r, s r I t t t t r, s r I t t t sand I are functions of... Hwang 016 JPO Hwang and Walsh 016 JPO (in press) 6

7 yh (km) U Hs x x r, s x r I x Tp x x r, s x r I x t t r, s t r I t t t r, s t r I t s and I are functions of and rm xh (km) Harmonics of s and I q an, q cos n bn,q sin n N n 0 q can be s x, I x, s x, I x, s t, I t, s t or I t Y p1y rm py Y represents an,q and bn,q

8 8 H s _fetch E t _fetch H s _duration E t _duration T p _fetch M t _fetch T p _duration M t _duratio n

9 Input: HRD U Output1: H sw, T pw, #, #, Output: E, E t, M, M t Output0/Input1: modeled fetch/duration 9

10 Ivan 004 Fan et al. 009 JPO Hwang and Walsh 016 JPO (in press) M C U E tuh H a 5 H ; C -0.16U 9.67U U ; tuh E a E Liu et al. 008 JPO M C U C tul L a 0.069U L E 3.5 u ; u C U tul a * * L ; Product of wind wave growth function applied to HRD wind time series

11 11 Bell et al. 01 JAS (calculated by conservation of azimuthally averaged absolute angular momentum) M U ; t M a M # # C C C a d U U a U / a /

12 Hurricane Wind Waves and Air-Sea Interaction Waves inside hurricanes are fetch- and duration-limited Full set of (U, H s, T p ) can be derived knowing only one, given fetch or duration Wind sea: energy and momentum exchange = f(u, H s, T p ), simple parameterization functions E U ; 0. ; M U ; t E a E # # t M a M # # Temporal and spatial evolution of hurricane wind wave development from time series of hurricane wind field

13 13

14 14

15 3. ESTIMATING AIR-SEA ENERGY AND MOMENTUM EXCHANGES INSIDE TROPICAL CYCLONES USING THE FETCH- AND DURATION-LIMITED WAVE GROWTH PROPERTIES Paul A. Hwang 1 and Yalin Fan 1 Remote Sensing Division, Naval Research Laboratory, Washington, DC, USA Oceanography Division, Naval Research Laboratory, Stennis Space Center, MS 1. INTRODUCTION For wind-generated waves, the wind-wave triplets (reference wind speed U, significant wave height Hs, and frequency spectral peak wave period Tp) are intimately connected through the fetch- or duration-limited wave growth functions (e.g., Hasselmann et al. 1973; Donelan et al. 1985; Kahma and Calkoen 1994; Young 1999; Hwang and Wang 004). The full set of the wind-wave triplets can be obtained knowing only one of the three variables combined with the fetch xf (duration td) information using the fetch-limited (duration-limited) wave growth functions, which constitute a pair of dimensionless equations describing the growth of wave height and wave period as a function of fetch or duration. The frequency-integrated air-sea energy and momentum exchange rates Et and Mt are functions of the wind-wave triplets and can be quantified with the wind-wave growth functions (Hwang and Sletten 008; Hwang 009; Hwang and Walsh 016). Previous studies have shown that wave development inside hurricanes follows essentially the same growth functions established for steady wind forcing conditions (Young 1988, 1998, 006; Hwang 016; Hwang and Walsh 016). Here we present the analysis of wind-wave triplets collected inside Hurricanes Bonnie 1998 and Ivan 004 at category to 4 stages using the airborne scanning radar altimeter (SRA) combined with the NOAA HRD hurricane wind velocity product (Powell et al. 1996; Wright et al. 001; Moon et al. 003; Fan et al. 009). The data yield the windwave triplets inside the hurricanes along several (between 6 and 11) transects radiating from the hurricane center. Using the wind-wave triplets from the 4 hurricane scenes (Table 1), a parameterization model is formulated for the effective fetch and duration at any location inside a hurricane. The fetch and duration model is applied to the time series of D hurricane wind fields from post analysis to investigate the detailed temporal evolution of the wave field and the associated energy and momentum exchanges over the hurricane coverage area. The spatial distributions of energy and momentum exchanges show considerable asymmetry. Referenced to the hurricane heading, the exchanges on the right half plane of the hurricane are much stronger than those on the left half plane, the right-to-left ratio of about :1 to 3:1 is common.. HARMONIC ANALYSIS OF HURRICANE FETCH AND DURATION The fetch- and duration-limited wave growth functions connect the wind-wave triplets (U, Hs, Tp). Many different formulas of growth functions have been reported in the literature, in this study we use the results obtained from least squares fitting of five field experiments conducted under quasi-steady and near-neutral stability conditions, collectively called BHDDB (Burling 1959; Hasselmann et al. 1973; Donelan et al. 1985; Dobson et al. 1989; Babanin and Soloviev 1998), as described in Hwang and Wang (004): x 7 fg x fg s U 6.19, U U U H g s 8 d d 1.7,.94 4 U U U, (1) H g t g U t g, () 16 where g is the gravitational acceleration (9.8 m/s ). 4 4 H g /16 U g / U represents the wave energy s E g w rms # rms, and U / U / cp # is dimensionless frequency or inverse wave age; w is water density. The wind and wave measurements from the hurricane hunter missions provide the necessary data to calculate the effective fetches and durations for the locations where the wind and wave data are acquired (Hwang 016; Hwang and Walsh 016). Explicitly, the fetch and duration are derived by rearranging the variables in (1) and () x 4.4 U H, x.9 U T, (3) x s x p t 1.75 U H, t 4.81 U T. (4) t s t p The analyses of Hwang (016) and Hwang and Walsh (016) show that xf and td can be represented by linear functions of the radial distance r from the hurricane center: x s r I, x s r I, (5), x x x x x x t s r I t s r I, (6) t t t t t t where is the azimuth angle referenced to the hurricane heading, positive counterclockwise (CCW). The intercepts I and slopes s of the linear functions in (5) and (6) from processing the four hurricane scenes are shown in Fig. 1; these fitting parameters can be expressed in Fourier series: N n, q cos n, q sin, (7) n0 q a n b n where q can be s x, Ix, s x, Ix, s t, It, s t or I t. The curves computed with N=1, and 3 are shown with continuous lines in Fig. 1; the data are well represented by (7) with N=3. The harmonics anq and bnq display a systematic quasilinear variation with the radius of maximal wind speed rm (Fig. ) Y p r p, (8) 1Y m Y where Y represents an,q and bn,q in (7). The fitting coefficients p1y and py are listed in Table (the headers s_ex, I_ex, s_ox, I_ox, s_et, I_et, s_ot and I_ot represent s, I, s, I, s, I, s and I x x x x t t t t. 3. DERIVATION OF WAVE PARAMETERS AND AIR-SEA EXCHANGES FROM HURRICANE WIND FIELD As discussed in Hwang (016) and Hwang and Walsh (016), for wind generated waves, (U, Hs, Tp) are

16 connected by the growth functions and the full set can be computed knowing only one of the three, coupled with the fetch or duration input. At the present, the hurricane wind field is the most available product so we illustrate the computation of wave parameters and the derived quantities such as Et and Mt from the wind input. Rewriting (1), the equations for Hs and Tp with U and fetch inputs are H 8. U x, T 9.8 U x. (9) sw x pw x Similarly, for U and duration inputs (), H 1.55 U t, T 3.53 U t. () sw t pw t The subscript w is appended to the wave variables in (9) and () to emphasize that the wind-sea portion is obtained with the fetch-limited wave growth functions. The fetch and duration in those equations are computed from the harmonic model described in section. Once (U, Hs, Tp) are available, the #, #, Mt and Et can be calculated (Hwang and Sletten 008; Hwang and Walsh 016): Et EaU ; E 0.0 # #. (11) 4.3 M U ; 0.40 t M a M # # where a is the air density. Fig. 3 shows an example (approximately the middle point of case I14, 00Z 15 Sep 004) of using the HWind U to compute the wave parameters and air-sea exchange rates. If a long time series of hurricane wind field is available, the corresponding temporal variation of the wave field and air-sea exchanges can also be computed. Fig. 4 presents an application applied to Hurricane Ivan 004 (Fan et al. 009), showing (a) the energy exchange integrated over the hurricane coverage area (here defined as the circle with a 50-km radius), (b) the relative energy exchange contributions of the left- and right half-planes, (c) the momentum exchange integrated over the hurricane coverage area, and (d) the relative momentum contributions of the left- and right half-planes. These results show the temporal evolution and the asymmetric spatial distributions of energy and momentum exchanges, the right-to-left ratio of :1 to 3:1 is common. 4. SUMMARY The wave fields inside hurricanes behave similarly to those generated by steady fetch- or duration-limited wind forcing conditions (Young 1998; 006; Hwang 016; Hwang and Walsh 016). The full set of the wind-wave triplets (U, Hs, Tp) can be calculated with the fetch- or duration-limited growth functions knowing only one of three variables and accompanied with the fetch or duration information. Using the U, Hs and Tp measured inside four different hurricane scenes, we present a fetch and duration model for any location inside a hurricane. The Et and Mt of a wind wave system are functions of the wind-wave triplets, and can be calculated as well. The procedure is very useful for investigating the spatial distribution and temporal evolution of the air-sea interaction processes inside hurricanes. References Babanin, A.V., and Y. P. Soloviev, 1998: Field investigation of transformation of the wind wave frequency spectrum with fetch and the stage of development. J. Phys. Oceanogr., 8, Black, P. G. et al., 007: Air-sea exchange in hurricanes: Synthesis of observations from the coupled boundary layer air-sea transfer experiment. Bulle. Am. Meteorol. Soc., 88, Burling, R. W., 1959: The spectrum of waves at short fetches. Dtsch. Hydrogr. Z., 1, Dobson, F., W. Perrie, and B. Toulany, 1989: On the deepwater fetch laws for wind-generated surface gravity waves. Atmos.-Ocean, 7, -36. Donelan, M. A., J. Hamilton, and W. H. Hui, 1985: Directional spectra of wind-generated waves, Phil. Trans. Roy. Soc. Lond., A315, Fan, Y., I. Ginis, T. Hara, C. W. Wright, and E. J. Walsh, 009: Numerical simulations and observations of surface wave fields under an extreme tropical cyclone. J. Phys. Oceanogr., 39, Hasselmann, K. et al., 1973: Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Deutsch. Hydrogr. Z., Suppl. A8(1), 95pp. Holthuijsen, L. H., M. D. Powell, and J. D. Pietrzak, 01: Wind and waves in extreme hurricanes. J. Geophys. Res., 117, C Hwang, P. A., 006: Duration- and fetch-limited growth functions of wind-generated waves parameterized with three different scaling wind velocities. J. Geophys. Res., 111, C0005. Hwang, P. A., 016: Fetch- and duration-limited nature of surface wave growth inside tropical cyclones: With applications to air-sea exchange and remote sensing. J. Phys. Oceanogr., 46, Hwang, P. A., and M. A. Sletten, 008: Energy dissipation of wind-generated waves and whitecap coverage. J. Geophys. Res., 113, C001. Hwang, P. A., and D. W. Wang, 004: Field measurements of duration limited growth of wind-generated ocean surface waves at young stage of development. J. Phys. Oceanogr., 34, Hwang, P. A., and E. J. Walsh, 016: Azimuthal and radial variation of wind-generated surface waves inside tropical cyclones. J. Phys. Oceanogr. (in press). Janssen, P., 004: The Interaction of Ocean Waves and Wind. Cambridge Univ. Press, Cambridge, UK, 300 pp. Kahma, K. K., and C. J. Calkoen, 1994: Growth curve observations. In Dynamics and modeling of ocean waves, (eds.) G.J. Komen, L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, and P.A.E.M. Janssen, Cambridge University Press, Komen, G. J., L. Cavaleri, M. Donelan, K. Hasselmann, S. Hasselmann, and P. A. E. M. Janssen (Eds.), 1994: Dynamics and modeling of ocean waves. Cambridge Univ. Press, Cambridge, UK, 53 pp. Moon, I.-J., I. Ginis, T. Hara, H. L. Tolman, C. W. Wright, and E J. Walsh, 003: Numerical simulation of sea surface directional wave spectra under hurricane wind forcing. J. Phys. Oceanogr., 33, Pierson, W. J., and L. Moskowitz, 1964: A proposed spectral form for full, developed wind seas based on the similarity theory of S.A. Kitaigorodskii. J. Geophys. Res., 69, Powell, M. D., S. H. Houston, and T. A. Reinhold, 1996: Hurricane Andrew s landfall in south Florida. Part I: Standardizing measurements for documentation of surface wind fields. Wea. Forecasting, 11,

17 Wright, C. W., E. J. Walsh, D. Vandemark, W. B. Krabill, A. W. Garcia, S. H. Houston, M. D. Powell, P. G. Black, and F. D. Marks, 001: Hurricane directional wave spectrum spatial variation in the open ocean. J. Phys. Oceanogr., 31, Young, I. R., 1988: Parametric hurricane wave prediction model, ASCE J. Waterway, Port, Coastal and Ocean Eng., 114, Young, I. R., 1998: Observations of the spectra of hurricane generated waves. Ocean Eng., 5, Young, I. R., 1999: Wind generated ocean waves. Elsevier, Amsterdam, the Netherlands, 88 pp. Young, I. R., 006: Directional spectra of hurricane wind waves. J. Geophys. Res., 111, C0800. Table 1. Some basic information of the four datasets, collected in hurricane Bonnie 1998 and Ivan 004, used for the analysis in this paper. File ID B4 I09 I1 I14 Start time UTC 8/4/1998 0:9 9/9/004 16:15 9/1/004 :39 9/14/004 0:09 End time UTC 8/5/1998 1:44 9/9/004 0: 9/1/004 15:41 9/14/004 :49 HRD U max (m/s) rm (km) V_a (m/s) theta_a (degn) SRA U min SRA U max SRA Hs min SRA Hs max SRA Tp min SRA Tp max # SRA spectra Table. Linear fitting coefficients of fetch and duration harmonics as a function of radius of maximum winds Y p1 Y rm p Y (8), where Y represents an,q and bn,q in (7); see text for further explanation. a0,q q s_ex I_ex s_ox I_ox s_et I_et s_ot I_ot p1-1.40e E E-0 3.E E E E E-01 p 1.57E E E E E-0 -.4E E E+00 a1,q p1-1.0e-0 1.6E E E E E E E-0 p 5.E E E-01-5.E+01.39E-0 -.3E E E+00 b1,q p1 5.1E E E E E E E E-0 p -6.87E-0 3.0E E E+0 5.5E E E E+00 a,q p1 5.55E E E-0-1.E+00 3.E E E E-0 p -4.79E-01.5E E E E-0 1.7E E-0 3.4E+00 b,q p1-3.3e e E E E E E E-0 p.11e E E E E E E E+00 a3,q p1 6.44E E E-0-9.1E E E-0 6.3E E-0 p -3.9E-01.85E E E E E E-0.18E+00 b3,q p1 1.91E E E E E E E E-0 p -.5E-01.11E E E E-0 8.7E E E+00

18 Fig. 1. Slope sq and intercept Iq of q s r I, where q can be x, x, t or t. The results obtained from datasets B4, I09, I1 and I14 (Table q q x x t t 1) are shown in the 4 quadrants, respectively upper left (UL), upper right (UR), lower left (LL) and lower right (LR). The solid, dashed and dotted curves corresponding to the first, second and third order Fourier series are also superimposed in each panel. For each dataset, s and I are plotted in (a) and (b) for xx, (c) and (d) for xx, (e) and (f) for tt, and (g) and (h) for tt. N Fig.. Harmonic representations q an, q cos nbn, q sin n of the slopes and intercepts of the linear functions n0 defining the fetch and duration, where q can be s, I, s, I, s, I, s or I x x x x t t t t. The harmonics an,q and bn,q show linear dependence on rm, n=0, 1,, 3 (b0,q=0). (UL): Slopes of fetch for wave height and wave period: (a) an,sx, (b) bn,sx, (c) an,sx, and (d) bn,sx; (UR): Intercepts of fetch for wave height and wave period: (a) an,ix, (b) bn,ix, (c) an,ix, and (d) bn,ix; (LL): Slopes of duration for wave height and wave period: (a) an,st, (b) bn,st, (c) an,st, and (d) bn,st; (LR): Intercepts of duration for wave height and wave period: (a) an,it, (b) bn,it, (c) an,it, and (d) bn,it.n.

19 Fig. 3. An example illustrating the derivation of wave and air-sea interaction parameters using the hurricane wind field as an input, shown here are (a) input U, and outputs of (b) Hsw, (c) Tpw, (d) #, (e) Et, (f) Mt; subscript w emphasizes that the wave parameters derived from the fetch- and duration-growth functions are the wind sea components.