Modeling Microseism Generation by Inhomogeneous Ocean Surface Waves in Hurricane Bonnie Using the Non-Linear Wave Equation

Abstract: It has been shown that hurricanes generate seismic noise, called microseisms, through the creation and non-linear interaction of ocean surface waves. Here we model microseisms generated by the spatially inhomogeneous waves of a hurricane using the non-linear wave equation where a second-order acoustic field is created by first-order ocean surface wave motion. We treat range-dependent waveguide environments to account for microseisms that propagate from the deep ocean to a receiver on land. We compare estimates based on the ocean surface wave field measured in hurricane Bonnie in 1998 with seismic measurements made roughly 1000 km away in Florida.

Physically the process for generating microseisms can be broken down into four steps: (1) The spatially and temporally varying hurricane winds create complicated wave patterns on the ocean surface; (2) In locations where these surface waves are in opposition, their non-linear interactions generate acoustic waves; (3) These acoustic waves propagate through the water and couple into the ocean bottom; (4) The resulting seismic waves propagate along the sea bottom where they may be detected at great distances by seismic sensors.In this paper we develop an approach for estimating the microseisms generated by a hurricane and their seismo-acoustic propagation.
Our modeling approach is based on the non-linear wave equation where a second-order acoustic field is generated by a source distribution that depends on the first-order ocean surface wave motion.The acoustic (or seismic) field at a receiver can then be expressed as the integral over the source distribution multiplied by the waveguide Green function.This approach is ideal for hurricane-generated microseisms since it can be applied to spatially inhomogeneous surface wave fields.Also, this approach may be used in range-dependent waveguide environments as is the case when a hurricane at sea generates microseisms that propagate up the continental margin to a receiver on land.
As an example, we use the ocean surface directional wave spectrum from hurricane Bonnie [12,13] to estimate the microseismic source levels.We than calculate the Green's function to describe the seismo-acoustic propagation from Bonnie's location in the North Atlantic to a seismometer in Florida.Our calculated levels are then compared against actual measurements.
Our approach builds on a long history of microseism modeling efforts.Initially, microseism generation was modeled for idealized situations [14] where the ocean surface wave field was taken to be spatially homogeneous and the ocean waveguide to be range-independent.Previous microseismic models, however, cannot be applied in typical hurricane scenarios because the surface wave fields are inhomogeneous in that the wave height spectra in different parts of the hurricane can vary both in magnitude and directionality.In some cases this is because spatial homogeneity is assumed over infinite [15][16][17][18][19][20] or very large [14,21,22] surface wave areas.One model that does account for finite microseismic source generation regions is limited by the assumption that the receiver is at the center of the source area [23,24].Also, these previous models are not applicable to the range-dependent environments typical in hurricane measurements where microseisms generated by a hurricane over the deep ocean are measured by sensor on land.Our non-linear wave equation method is shown to agree with earlier models [18,20] if we make the same simplifying assumptions that the source generation region is spatially homogeneous and that the waveguide can be modeled as an infinite half-space.
While the treatment of inhomogeneous surface wave fields and range-dependent waveguides presented here had not been previously considered, the medium non-linear wave equation used here has a parallel derivation to the wave equations or hydrodynamic equations used in earlier works.Longuet-Higgins [14], Hasselmann [17] and Brekhovskikh [15] base their derivations on perturbing Bernoulli's equation while Hughes [18] and Kibblewhite and Wu [19] begin by separately perturbing the equations of momentum, continuity and state.We show that the non-linear wave equation used here is equivalent to perturbation expansions [18] for the physical parameters typical in hurricane microseisms.Lloyd [20] and Cato [23] base their derivations on Lighthill's equation; however, Lloyd [20] shows that both perturbation and Lighthill approaches yield the same end result.We also note that the non-linear wave equation can be derived as a second-order approximation of Lighthill's equation [25].
These second-order non-linear theories should not be confused with the linear theories proposed by Banerji [26] and Bowen et al. [27].These linear theories claim that the first-order motion of a surface gravity wave creates a first order pressure fluctuation on the sea floor regardless of how deep the ocean is.This is contrary to classic surface wave theory which shows that first-order wave motion decays exponentially with depth such that, in deep water, the first-order pressure fluctuation on the sea floor goes to zero [14,[28][29][30].
Our derivation is in agreement with earlier works [14,31], which show that microseisms are generated by the non-linear interaction of ocean surface waves with roughly the same wavelength but opposing propagation directions.Measurements [13,32] and models [12] of surface directional wave height spectra in hurricane Bonnie show complex patterns with the opposing surface waves necessary to generate microseisms.Based on the wave height spectra in hurricane Bonnie, we calculate the microseismic source levels generated by the non-linear interaction of the ocean surface waves.
It has been shown both theoretically [10,14] and experimentally [4,11,[33][34][35][36][37] that microseisms propagate as Rayleigh waves along the sea floor.Here we develop a normal-mode model for the propagation of Rayleigh waves along the seafloor assuming a water depth that changes with range.The Green's function from this model is used to estimate the propagation of the microseism from deep water to a seismometer on land.
Given the wave height spectra in hurricane Bonnie, the calculation of the microseismic source level generated by the non-linear wave interaction, and the propagation of those waves along the sea floor, we estimate the noise levels at a seismometer in Florida.These estimates are then compared against measured data.
It should be noted that there are other sources of microseisms besides hurricanes.Microseisms have been associated with other non-hurricane storms, coastal regions where complex wave patterns are prevalent, and the typical ambient ocean wave spectra [11,[38][39][40][41][42][43][44].These other sources of microseisms are not considered here.

Ocean Surface Gravity Waves
Hurricanes are characterized by high winds that can vary quickly with position, both in direction and speed, as shown in Figure 1A.In addition a hurricane may move at speeds up to 15 m/s [45] so that the winds at any location can also change with time.These spatially and temporally varying winds generate complex ocean surface wave directional spectra Figure 2 with wave heights that can exceed 10 m [12,13,32].From Equation (37) we expect these locations with opposing waves to produce the greatest microseismic source levels and shown in Figure 1B.
For surface gravity waves that are homogeneous over an infinite ocean surface area we can express the complex surface wave height as the linear superposition of plane waves [46] where κ is the sea surface wavenumber and σ is the corresponding frequency where in deep water σ 2 = gκ.Throughout this paper the surface wavenumber κ is expressed either in cartesian coordinates as (κ x , κ y ) or in polar coordinates as (κ, θ).Also, a cartesian coordinate system is used for position where r = (ρ, z) = (x, y, z) where z is defined downward from the sea surface.The complex surface wave height can then be written as where A(κ) is the surface wave height spectrum.Since the complex wave height ζ(ρ, t) is produced by the contributions of many independent random physical processes, we can assume by Central Limit Theorem that the statistical distribution of ζ(ρ, t) is Gaussian [47].In addition we may define the reference position from which wave height is measured such that ζ(ρ, t) has zero mean.Given this assumption, in Appendix A.1 we show that A(κ) must also be a zero-mean complex Gaussian random variable [46,47].
Since the spectrum A(κ) is zero-mean and homogeneous over an infinite ocean surface area, the second moment of the spectrum can be written as (see Appendix A.1) [16,47,48] where <> denotes the expectation.This expression shows that, for stationary processes, the wavenumber components κ and κ decorrelate.Also, since the spectrum A(κ) is Gaussian, the fourth moment can then be written as [16] < In a hurricane the surface waves are not homogeneous over an infinite area, but rather the surface wave spectrum A(κ) changes gradually with position in the storm.To characterize the spatially varying surface wave field in the hurricane, we divide the sea surface into a grid made up of finite regions m with surface areas S m over which the wave spectrum A m (κ) can be taken to be homogeneous.The complex wave height in a particular finite region m then can be expressed as In Appendix A.2 we show that the finite size of the regions m introduces a 'windowing' effect when calculating the moments of the spectrum A m (κ).However, we also show that this windowing effect can be neglected if the dimensions L x and L y of the region are much greater than the wavelength λ g of the surface wave.For the hurricane waves of interest in this paper, the surface wavelengths λ g range from roughly 100 to 300 m requiring that the dimensions L x and L y be in the order of 1 km or greater.Provided this condition is met we can make the approximation.
where S mn (κ) is the cross-spectral density of the surface wave fields at locations m and n.
As before we can assume by Central Limit Theorem that ζm (ρ, t) is a zero-mean complex Gaussian random variable so that A m (κ) is also a zero-mean complex Gaussian random variable.Because of this we can write the fourth moment as This expression for the fourth moment of the wave height spectra will be used later in Section 2.3.
From Equation ( 4) we can write the real part of the surface wave height as We can also express the real surface wave particle velocity v as a linear superposition of plane waves such that where îx , îy and îz represent unit vectors in the x, y and z directions respectively.Experiments by Forristall et al. [49] show that the linear wave theory of Equations ( 1), ( 4), ( 7) and ( 8) is adequate to describe surface waves even in high sea states with significant wave breaking as in a hurricane.The use of linear surface wave theory can also be justified since the typical surface wave heights in a hurricane are an order of magnitude less than the wavelength [12,13,32,49].
Wright et al. [13] and Walsh et al. [32] measured the spatial variation of the directional wave spectra A m (κ) in hurricane Bonnie using an aircraft-mounted Scanning Radar Altimeter (SRA).These measurements, however, are limited to the locations and times of the aircraft flights and do not give a complete picture of the wave field in the hurricane.
To quantify the surface wave field at any location and time in hurricane Bonnie, Moon et al. [12] used the Wave Watch III (WW3) ocean surface wave model [50].These model results, shown in Figure 2, are in close agreement with measurements from [13,32] and with data from buoys and oceanographic stations.In this paper we use the surface directional wave spectra A m (κ) in hurricane Bonnie calculated using WW3 [12].

Non-Linear Wave Equation
The non-linear wave equation [51,52] describes the second-order acoustic wave field generated by first order fluid motion.For the long surface waves in a hurricane viscosity may be neglected so that the non-linear wave equation can be written as [53] ∇ where the first order velocity v 1 and pressure p 1 terms on the right-hand side generate a second-order acoustic pressure field p 2 on the left-hand side, c is the sound speed, ρ is density and Γ is the coefficient of nonlinearity [51,53].
In this work we will consider the case where the first order velocity v 1 , from Equation (8), is due to surface wave motion.The first order pressure p 1 can be found using conservation of momentum where −∇p 1 = ρ∂v 1 /∂t.Using the relationship between first order velocity and pressure we can compare the relative magnitudes of the source terms in Equation (9).Given typical parameters κ = 0.02 to 0.7 rad/m, σ = 0.4 to 0.8 rad/s, c ≈ 1500 m/s, ρ ≈ 1000 kg/m 3 , and Γ ≈ −2.6 [54] we find that the 9) exceeds the other source terms by several orders of magnitude.Dropping the lesser terms yields This equation is equivalent to the perturbation expression used in the previous microseism derivation of Hughes [18].This can be seen by taking Equation (10) of [18] and assuming irrotational flow Taking the Fourier Transform of Equation ( 10) we find the frequency domain Helmholtz equation for the second-order field We will begin by solving for the pressure field generated by a finite volume V m with a source distribution q m which depends on the fluid velocity v 1,m from Equation (8).The volume V m extends vertically from the ocean surface at z = 0 to a depth well below the surface wave region (z > λ g = 2π/κ).The volume extends over the horizontal area S m with dimensions L x and L y centered at ρ m = (x m , y m ) as defined in Section 2.1.We can solve the Helmholtz (Equation ( 11)), using the frequency domain Green function g(r r , r s , ω) such that, where r s = (x s , y s , z s ) and r r = (x r , y r , z r ) are the source and receiver locations respectively.It should be noted that, although the ocean surface is perturbed by the presence of surface waves, we choose to define the upper boundary of our volume by the quiescent surface z = 0 following the example of Lloyd [20].He argues that the Green function g "has the form g = g 0 + g 1 + g 2 + . . .where g 0 is the Green function for the quiescent ocean and where g j , j ≥ 1, is of degree j in a wave height perturbation parameter < h 2 > 1/2 .The first-order term g 1 represents Bragg reflected by the moving spatially periodic components".He shows that the terms g j for j ≥ 1 may be omitted in the excitation of p 2 .He goes on to apply the Green function for a quiescent ocean half-space and, even with this approximation, matches the result of earlier work by Hughes [18].In Appendix C we will demonstrate agreement between our end result and that of both [20] and [18].
Equation ( 12) can be rewritten as where the window function w(ρ s − ρ m ) is unity for −L x /2 < x s − x m < L x /2 and −L y /2 < y s − y m < L y /2 and zero elsewhere.The double integral over ρ s represents the field at a receiver r r due to a distribution of sources q over a finite area defined by the window function w.If the receiver r r is in the far field of the finite area y /λ a where λ a is the acoustic wavelength and where r m = (ρ m , z s ) = (x m , y m , z s )), we can make the plane wave approximation, where the acoustic wavenumber is defined in cartesian coordinates as k = (k r , k z ) = (k x , k y , k z ) where In free space this would simplify to however, we will use the more general form given in Equation (14).Inserting Equation ( 14) into Equation ( 13) yields Again the double integral over ρ s analogous to a horizontal planar array of acoustic sources with broadside directed down toward the ocean bottom.In the following analogy we can define this double integral as the array output [55] The array output in Equation ( 17) can be evaluated for any horizontal wavenumber vector k r ; however, only those wavenumber vectors where |k r | ≤ ω/c lead to a propagating acoustic field.This can be illustrated if we consider the microseismic field generated by a horizontal plane of sources near the sea surface as described in Equation (17).In previous works [14] it has been shown that microseisms generated at the ocean surface are first transmitted downward to the sea floor where they then propagate along the bottom as Rayleigh waves.This vertical propagation from the sea surface to the sea floor can be modeled to first order using the free-space Green function Assuming the acoustic pressure field given in Equation ( 16) is temporally stationary, the power spectral density may be written as [56] where <> represents the expectation and r n = (ρ n , z s ) = (x n , y n , z s ).In the next section we derive expressions for the second moment < B m (k, z s,m )B * n (k , z s,n ) > of the array output and the power spectral density S P m,n due to ocean surface waves.
The power spectral density S P m ,n in Equation (19) represents the field at a receiver from a single pair of finite volumes V m and V n .In the next section we determine the microseismic field generated by a hurricane where the surface wave field is inhomogeneous and extends over a large region hundreds of kilometers across.To do this we divide the hurricane region into finite volumes V m , as shown in Figure 3, and then sum their contributions to find the total power spectral density of the received field written as

Power Spectral Density Due to Ocean Surface Gravity Waves
Taking the ocean surface wave velocity as the first order field we can determine the source term from Equation (11) where q(r s , ω) is the Fourier transform Given Equation ( 8) for the first order velocity we find where , and κ − = (κ 1 − κ 2 ).Note that the Laplacian operator (∇ 2 ) in Equation ( 23) leads to the direction cosine functions and where ), and κ − y = (κ 1y − κ 2y ).
Substituting Equation ( 23) into the definition of q m (r s , t) from Equation ( 22) and taking the Fourier transform yields, By inserting Equation ( 26) into the definition of the array output B m (k, z s ) from Equation ( 17) we obtain Note that the wave height spectra A m and A n in Equation ( 27) can be brought out of the spatial integral since they are constant over the ocean surface area S m defined by the window function w.Integrating over ρ s then leads to where * represents the two-dimensional convolution written as, and where Note that, when integrating over ρ s , the ρ m term in the window function w(ρ s − ρ m ) of Equation ( 27) introduces a phase factor e −ik r •ρ m .This cancels the other phase term e ik r •ρ m in Equation (27).From this we see that the source array output given in Equation ( 28) has no phase dependence on location ρ m .
The delta functions in Equation (28) show that the horizontal component k r of the acoustic wavenumber vector is generated by and equal to either the sum κ + or difference κ − of the surface wavenumber.In Section 2.2 we discussed how the non-propagating field where |k r | > ω c is not significant and may be ignored.Therefore, when integrating over κ 2 we need only consider a range of integration Ω corresponding to propagating acoustic waves.For the sum terms in Equation ( 28) Ω corresponds to the integration region where |κ + | ≤ ω c and for the difference terms Ω corresponds to |κ − | ≤ ω c .Here we will define the source area dimensions L x and L y to be 'acoustically small' (L x , L y << λ a ).The acoustic wavelengths of the microseismic field in water range from roughly 6 to 11 km requiring that the dimensions L x and L y be on the order of 1 km or less.Note that this requirement on L x and L y , combined with the requirement that L x , L y >> λ g from Section 2.1, means that λ g << L x , L y << λ a .If L x , L y << λ a , the function W(k r ) ≈ L x L y in Equation ( 30) within the range of acoustically propagating wavenumbers k r ≤ ω/c.With this approximation Equation ( 28) can be written as The power spectral density of the pressure in Equation ( 19) contains the second moment or variance of the array output B m (k, z s ) which from Equation ( 31) can be expressed as, Since the wave height spectra A m and A n are uncorrelated zero-mean Gaussian random variables as discussed in Section 2.1, we can substitute Equation ( 6) into Equation ( 32) and integrate over κ 3 and Then by integrating κ 2 in Equation ( 33) over the integration region Ω, corresponding to |κ + 4S mn (κ 1 )S mn (κ 1 )δ(ω)δ(ω ) Since ω c << κ this integration leads to the approximations κ 1 ≈ −κ 2 in the sum terms and κ 1 ≈ κ 2 in the difference terms.
Substituting Equation (35) into Equation ( 19) and integrating over ω leads to the power spectral density of the pressure field where we define the 'microseismic source cross-spectral density' as Summing the contribution from all finite volumes V m gives us the total received power spectral density Equation ( 38) provides an analytic expression for microseisms generated by inhomogeneous ocean surface wave fields, as opposed to previous formulations where spatial homogeneity is assumed over infinite [15,[17][18][19][20] or very large [14,21,22] surface wave areas.Also the Green functions g(r r , r m , w) and g * (r r , r n , w) in Equation ( 38) may be calculated for any arbitrary range-dependent or range-independent ocean waveguide using standard propagation models.Later we will use a formulation for range-dependent Rayleigh wave propagation to calculate the Green functions and model microseismic propagation in a typical North Atlantic waveguide environment.
While Equation ( 38) is applicable to inhomogeneous surface wave fields and arbitrary ocean waveguides, it can also be applied to simpler homogeneous surface wave fields and Green functions.In fact, in Appendix C we find that, given an infinite half-space and assuming the surface wave spectrum is range-independent (S mn (ω/2, θ) = S(ω/2, θ)), the power spectral density of the pressure field in Equation (38) simplifies to the results derived by Hughes [18] and by Lloyd [20].

Microseismic Source Levels in Hurricane Bonnie
In this section we show how Equation (37) can be used to find the microseismic source cross-spectral density generated by the surface wave field in a realistic hurricane.In these examples we use the wave spectra from hurricane Bonnie.Through this analysis we demonstrate the relationship between the wind speeds, surface wave spectra, and microseismic source spectra in a hurricane.
Beginning on 22 August 1998, hurricane Bonnie traveled along the east coast of the Bahamas, Florida, and South Carolina [58].During this time microseisms were recorded by an Incorporated Research Institutions for Seismology (IRIS) seismometer at the Disney Wilderness Preserve in Florida (28.1 • N, 81.4 • W).Here we model the microseisms generated by hurricane Bonnie that would be measured by this seismometer.
To evaluate the microseismic source cross-spectral density S B mn of Equation ( 37) requires knowledge of the surface wave height spectrum A m .Moon et al. [12] calculate the surface wave height spectrum for the wind speeds measured in hurricane Bonnie using the Wave Watch III (WW3) wind-wave modeling program [50].They also show how their results are in close agreement with aircraft-based measurements of surface wave height spectrum by Wright et al. [13] and Walsh et al. [32].
Figure 1A shows the surface wind speed in hurricane Bonnie on 24 August 1998.Illustrated is the typical hurricane structure with the high wind speeds of the 'eye wall' surrounding the 'eye' at the hurricane's center.From these wind speeds, WW3 is used to determine the wave height spectra A m (κ) as defined in Equation ( 5) (Figure 2).In Figure 2 we see that hurricanes generate waves that propagate in many different and often opposing directions.For example in Figure 2c we see a peak in the surface wave spectrum (in red) with waves traveling North ((κ x , κ y ) ≈ (0.0, 0.04)), West ((κ x , κ y ) ≈ (−0.04, 0.0)), and South ((κ x , κ y ) ≈ (0.0, −0.04)).
From the wave height spectra calculated using WW3, the microseismic source cross-spectral density S B mn may be calculated from Equation (37).We define the 'microseismic source level' as expressed in dB re Pa 2 s/m 4 where the somewhat arbitrary length scales L x and L y have been factored out. Figure 1B shows L B mm (ω, 0) at ω = 1.26 rad/s ( f = 0.2 Hz) on 24 August 1998, the day where we found the source level to be highest.We see that the peak of L B mm at location (c) is not at the same location as the maximum wind speed (a).This is because, while location (a) has large wave heights (≈37 dB re m 2 /(rad 2 /m 2 )), the waves are all propagating primarily to the South without any opposing waves.At location (c), however, there are opposing waves (propagating both North and South) which generate microseisms as expressed in Equation (37).Also the source level is higher in the low-wind-speed eye (b) than at the maximum wind speed location (a) due to the opposition in the surface waves (Figure 2b).The arrows in Figure 1 represent the direction hurricane Bonnie was moving and we see that there is a small peak in the source level well 'behind' the hurricane at location (d) even though the wind speed there is relatively low.This is again due to the opposing waves at this location (Figure 2d).This illustrates the complex relationship between wind speed and wave spectra, where the wave spectra is a function not only of wind speed but also of hurricane geometry and translation speed [12], and the complex relationship between wave spectra and microseismic source level, where the level depends on the non-linear interaction of opposing waves.

Results
In this section we model the propagation of the microseisms generated by hurricane Bonnie through the North Atlantic waveguide.The model results are then compared with seismic data gathered in Florida.
To determine the microseisms received at the sensor in Florida the Green functions g(r r , r m , w) and g * (r r , r n , w) of Equation ( 38) are calculated using the adiabatic Rayleigh wave propagation model derived in Appendix B. The bathymetry between the hurricane, located over the roughly 5 km deep Hatteras Abyssal plain on 24 August, and the seismometer, located in Florida is characterized by a gentle upslope.We will model this up-sloping environment with the simplified geometry shown in Figure 4.The sound speeds and densities for this model waveguide are based on typical values for the deep North Atlantic Ocean [59] and bottom [60].With this environment the Green functions g(r r , r m , w) and g * (r r , r n , w) are calculated using Equation (A47).The sea-floor depth at the center of hurricane Bonnie is given in Table 1, as well as the range between the seismometer and the hurricane center.The seismometer in Florida is buried in the ground to a depth of 162 m as shown in Figure 3.
Given the source cross-spectral density from Equation (38) and Green functions from Equation (A47), the power spectral density of the field, from Equation (38), received by the seismometer in Florida is calculated.Figure 5 shows the power spectral density of the horizontal velocity (Figure 5a-d) and vertical velocity (Figure 5e-h) of the earth's crust at the Florida seismometer based on our model (the red dashed curves) and on measured data (in blue).
The dotted blue lines in Figure 5 represent portions of the measured data that are dominated by non-hurricane related ambient noise.The level of this ambient noise is determined by measuring the noise levels during the week before and the week after the hurricane passed.The un-corrupted microseism signals are taken to be those that exceed the ambient by at least 3 dB while anything below that is considered to be corrupted.

Discussion
The theoretical results and measured data show reasonable agreement with peaks in the spectrum at roughly f = 0.2 Hz and a peak level between 50 and 60 dB re (nm/s) 2 /Hz.Note that the peak of the wave height spectral in Figure 2 is at roughly κ = 0.04 rad/m which equates to a surface wave frequency σ = 0.63 rad/s.From the frequency doubling effect discussed in Section 2.3, this gives a peak in the acoustic field at ω = 1.26 rad/s or f = 0.2 Hz which we see in the received field in Figure 5.
There is a lesser peak in the model result at around 0.1 Hz that is not obvious in the measured data.This peak in the model comes from the model for the Rayleigh wave propagation and suggests perhaps a limitation of the single layer (fluid water layer over a solid half-space) approximation adopted here.
It should be noted that the horizontal-and vertical-velocity measurements are consistent with Rayleigh waves.There is no evidence that other seismic wave types play a significant role.

Conclusions
Here we present an analytic model, based on the non-linear wave equation, to describe the microseisms generated by a hurricane.This model is ideal for hurricane-generated microseisms since it can be used to calculate the acoustic field due to spatially inhomogeneous surface waves.Also, this model may be used in range-dependent waveguide environments as is the case when a hurricane at sea generates microseisms that propagate up the continental margin to a receiver on land.This modeling is useful because microseisms are a primary cause of noise in seismic measurements [7,9,10] that raise the detection threshold for monitoring earthquakes [11] and tsunamis.
Based on the ocean surface directional wave spectrum in hurricane Bonnie [12,13], we predict the microseismic source levels generated by the non-linear interaction of the ocean surface waves.
We then model propagation of the microseismic field from hurricane Bonnie in the North Atlantic to a seismometer in Florida to hindcast the measured signal.We find that these results compare reasonably well with seismic measurements.equals zero [62,64] where the dependence of k n , ξ n , η n , ζ n and H on range ρ has been suppressed to simplify notation.Equation (A38) reduces to the equation where the roots k n (ρ) correspond to the horizontal wavenumbers for each mode n.For the case of constant bathymetry H(ρ) = H, Equation (A39) is identical to those given in earlier works [10,62,64]; however, for this more general range-dependent case the roots k n (ρ) can now vary with range ρ.Unlike a 'classic' Rayleigh wave which has only one root or propagating value, 'generalized' Rayleigh waves may have multiple roots k n (ρ) depending on frequency and fluid layer thickness as shown in Equation (A39) for the range-dependent case and as shown by [10] for the range-independent case.Note that Equation (A39) has no roots for k n (ρ) > ω/v 1 .
Given the roots k n (r) from Equation (A39), we can express the potentials as

10 Figure 1 .Figure 2 .
Figure 1.(A) Wind speed in m/s and (B) microseismic source level L B mm (ω, z s = 0) at ω = 1.26 rad/s ( f = 0.2 Hz) in dB re Pa 2 s/m 4 from Equation (39) at 1200 on 24 August 1998 as a function of latitude and longitude.The arrow indicates the direction hurricane Bonnie was moving.The letters (a-d) represent features of interest; (a) indicates the location of maximum wind speed, (b) indicates the eye where wind speed is zero,and (c,d) indicate peaks in the microseismic source level.This figure shows that, while a hurricane can produce significant microseismic source levels (B), these level do not directly follow the wind speeds (A) in the hurricane.
18) where d represents the water column depth.If |k r | exceeds ω/c then k z becomes imaginary leading to an exponential decay with depth d.For example if |k r | = 2ω/c then k z = i √ 3ω/c so that, for ω = 1.3 rad/s, c = 1500 m/s and an ocean depth d = 1000 m, the field is attenuated by e −k z z r ≈ 0.2 or 13 dB.Greater depths will lead to even more attenuation.Because of this evanescent attenuation we will only consider the acoustic field where |k r | ≤ ω/c.

Figure 3 .
Figure 3. Geometry of the hurricane wave field and ocean waveguide (not to scale).In our model the waveguide environment may be range-dependent and in this paper we consider the example of upslope propagation from the deep North Atlantic to Florida.The range and ocean depth parameters R and d are given in Table1.The depth of the receiver below the ocean bottom of 162 m corresponds to the depth below the earth surface of the actual seismometer in Florida.The compression wave speeds c p are 1500 and 5200 m/s in the water and bottom respectively.The shear wave speed c s in the bottom is 3000 m/s.The densities ρ in the water and bottom are 1.0 and 2.5 g/cm respectively.

Figure 4 .
Figure 4.The ocean depth (solid line) between hurricane Bonnie and the seismometer in Florida at noon on 24 August.Also shown is the ocean depth for the idealized up-sloping environment used to calculate the Green functions in Section 3 (dashed line).The scale of the figure makes the actual slope appear to change rapidly; however, the maximum slope of the ocean floor is roughly 3 • .

Figure 5 .
Figure5.Horizontal (a-d) and vertical (e-h) velocity modeled (red dashed) and measured (blue) at the seismometer in Florida at four times on 24 August 1998.The dotted blue lines represent portions of the data that where dominated by non-hurricane related ambient noise.Note that the peak in the spectra is at roughly f = 0.2 Hz.This frequency corresponds to the peak in the wave height spectrum at κ = 0.04 rad/m seen in Figure2.

Table 1 .
Parameters for Hurricane Bonnie on 24 August 1998.