Skip to main content

Acoustic-gravity waves in the atmosphere generated by infragravity waves in the ocean

Abstract

Infragravity waves are surface gravity waves in the ocean with periods longer than approximately 30 s. Infragravity waves propagate transoceanic distances and, because of their long wavelengths, provide a mechanism for coupling wave processes in the ocean, atmosphere, and the solid Earth. Here, we present a strict physical justification for the hypothesis that background ocean waves may generate waves in the upper atmosphere. We show that, at frequencies below a certain transition frequency of about 3 mHz, infragravity waves continuously radiate their energy into the upper atmosphere in the form of acoustic-gravity waves. Based on ionospheric observations and estimates of the fluxes of the mechanical energy and momentum from the deep ocean, we conclude that acoustic-gravity waves of oceanic origin may have an observable impact on the upper atmosphere. We anticipate our work to be a starting point for a detailed analysis of global manifestations of the ocean-generated background acoustic-gravity waves.

Findings

Introduction

Ionospheric manifestations of tsunamis (Artru et al. 2005; Makela et al. 2011; Galvan et al. 2012; Komjathy et al. 2012; Occhipinti et al. 2013) provided dramatic evidence of coupling between wave processes in the ocean and in the atmosphere. But tsunamis are only a transient feature of a more general phenomenon: infragravity waves (IGWs) (Webb et al. 1991; Herbers et al. 1995; Aucan and Ardhuin 2013; Godin et al. 2014). IGWs are permanently present surface gravity waves in the ocean with periods longer than the longest periods (approximately 30 s) of wind-generated waves. IGWs propagate transoceanic distances and, because of their long wavelengths (from approximately 1 km to hundreds of km), offer a mechanism for coupling wave processes in the ocean, atmosphere, and the solid Earth (Hines 1972; Crawford et al. 1998; Rhie and Romanowicz 2004; Bromirski and Stephen 2012). The notion that background ocean waves may generate waves in the atmosphere has existed for a long time (Hines 1972; Peltier and Hines 1976; Livneh et al. 2007), but no quantitative theory of this effect has been suggested. Here, we provide a strict physical justification for this hypothesis. Using a coupled atmosphere-ocean model taking into account both compressibility of the air and the gravity, we show that there exist two distinct regimes of IGW penetration into the atmosphere. At frequencies above a transition frequency, there are evanescent waves in the atmosphere propagating horizontally along the ocean surface. At frequencies below the transition frequency, IGWs continuously radiate their energy into the upper atmosphere in the form of acoustic-gravity waves (AGWs). The transition frequency is a function of ocean depth; it varies slowly near 3 mHz for typical depth values and drops sharply to zero only for extreme depth values. Using a semi-empirical model of the IGW power spectrum (Godin et al. 2013), we derive an estimate of the flux of the mechanical energy and mechanical momentum from the deep ocean into the atmosphere due to IGWs and predict specific forcing on the atmosphere in coastal regions. We conclude that AGWs of oceanic origin may have an observable impact on the upper atmosphere.

Theory

Infragravity waves and atmospheric waves

IGWs are usually modeled as surface gravity waves in an incompressible ocean with a free surface and a rigid seafloor. In this approximation, IGWs have the dispersion equation

$$ {\omega}^2=gk \tan hkH, $$
(1)

where ω and k are the wave frequency and wavenumber, g ≈ 9.8 m/s2 is the acceleration due to gravity, and H is the ocean depth (Gill 1982). IGW phase speed c = ω/k does not exceed (gH)1/2 and is much smaller than sound speed in water and compressional and shear wave speeds in the ocean bottom. Therefore, corrections to the dispersion equation due to water compressibility and seafloor compliance are rather small (Tsai et al. 2013; Watada et al. 2014) and do not change the nature of IGWs as horizontally propagating surface waves, where both k and ω are real-valued (Brekhovskikh and Godin 1998).

The free-surface boundary condition implies that the wave-induced motion of the air has no effect on oceanic waves. For a simple description of coupling between waves in the ocean and atmosphere, we model the latter as an isothermal gas half-space. Continuous waves with horizontal wave vector k on the sea surface generate plane waves in the half-space. Horizontal and vertical components of the wave vector of the atmospheric waves are k and m, where

$$ {m}^2=\frac{\omega^2}{c_s^2}-{k}^2-\frac{\gamma^2{g}^2}{4{c}_s^4}+\frac{\left(\gamma -1\right){g}^2{k}^2}{\omega^2{c}_s^2}, $$
(2)

k 2 = k∙k, c s is the sound speed in air, and γ ≈ 1.4 is the ratio of specific heats at constant pressure and constant volume (Gill 1982). The atmospheric waves have a finite impedance, i.e., the ratio of pressure perturbation in a fluid parcel to its vertical velocity. All these quantities must be continuous on the air-sea interface (Gill 1982). For oceanic waves, the impedance is small as long as the ratio R ~ 1.3 × 10−3 of the mass densities of the air and water at the water-air interface is small. We use a perturbation theory to quantify the coupling between oceanic and atmospheric waves.

Perturbation theory

Our theoretical analysis of the effects of the atmospheric load on infragravity waves is based on modeling the ocean surface as a boundary with finite impedance and application of the recently developed perturbation theory (Godin 2012) for normal modes of acoustic-gravity waves. The perturbation theory predicts the change in the dispersion relation of normal modes of acoustic-gravity waves in generic waveguides in horizontally stratified fluids, which results from small perturbations of the waveguide’s parameters.

Introduce a Cartesian coordinate system with horizontal coordinates x and y and a vertical coordinate z increasing upwards. Consider a fluid layer with a rigid boundary at z = −H and a compliant boundary that coincides with the horizontal plane z = 0 in the absence of waves. The fluid is stationary in the absence of waves. Let W 0(z, ω)exp(ik(ω)x − iωt) and P 0(z, ω)exp(ik(ω)x − iωt) describe the wave-induced vertical displacement of fluid parcels and pressure variation, which satisfy linearized (with respect to the wave amplitude) governing equations and boundary conditions W(z, ω) = 0, z = − H on the rigid boundary and P(z, ω) = ρgW(z, ω), z = 0 on the free surface (Gill 1982). Here, t is time, ω is wave frequency, k is wavenumber, and ρ is fluid density. The functions W 0(z, ω) and P 0(z, ω) describe a normal mode with the dispersion relation k(ω) in an unperturbed waveguide and are assumed to be known.

Now, consider a waveguide that differs from an unperturbed waveguide by having a different upper boundary. For waves with harmonic dependence exp(ikx − iωt) on horizontal coordinates and time, the boundary condition can be written as

$$ P\left(z,\omega \right)-\rho gW\left(z,\omega \right)=-i\omega Z\left(k,\omega \right)W\left(z,\omega \right),\kern0.48em z=0. $$
(3)

Here, P − ρgW and Z have the physical meaning, respectively, of the wave-induced pressure variation in a moving fluid parcel (i.e., Lagrangian pressure perturbation) and the boundary impedance (Godin 2012). Impedance of the free boundary equals zero. When Z is small, there exists a normal mode with the dispersion relation k 1(ω) = k(ω) + δk(ω) (Godin 2012), where

$$ \delta k=\frac{i{\omega}^3}{2k}{W}_0^2\left(0,\omega \right){\left[{\displaystyle \underset{-H}{\overset{0}{\int }}}\frac{dz}{\rho }{P}_0^2\left(z,\omega \right)\right]}^{-1}Z\left(k\left(\omega \right),\omega \right)+O\left({Z}^2\right)\;. $$
(4)

Note that the impedance in the right-hand side of Equation (4) is calculated for the unperturbed value of the normal mode wavenumber.

IGWs are usually modeled as surface gravity waves in an incompressible ocean with a constant water density, a free surface, and a rigid seafloor. In this approximation, one has (Gill 1982)

$$ P\left(z,\omega \right)=P\left(0,\omega \right)\frac{ \cosh {k}_0\left(z+H\right)}{ \cosh {k}_0H},\kern0.48em W\left(z,\omega \right)=\frac{P\left(0,\omega \right)}{\rho g}\frac{ \sinh {k}_0\left(z+H\right)}{ \sinh {k}_0H}, $$
(5)

and the unperturbed dispersion relation (1).

Lagrangian pressure perturbations and normal displacement are continuous at fluid-fluid interfaces (Gill 1982). Therefore, the impedance Z(k, ω) is also continuous and can be calculated from the ratio of the Lagrangian pressure perturbation and the vertical displacement in an atmospheric wave with the frequency ω and horizontal wavenumber k. The ocean surface impedance depends on atmospheric conditions. The simplest atmospheric model is that of a stationary, isothermal, perfect gas of constant composition occupying the half-space z > 0. In the isothermal atmosphere, acoustic-gravity waves (AGWs) have the dispersion relation (2). The lower boundary of the isothermal gas half-space has the impedance (Godin 2012; Godin and Fuks 2012)

$$ Z\left(k,\omega \right)=\frac{R\rho (0)}{\omega \left({\omega}^2{c}_s^{-2}-{k}^2\right)}\left[{\omega}^2m\left(k,\omega \right)+ig\left({k}^2-0.5{\omega}^2\gamma {c}_s^{-2}\right)\right], $$
(6)

where ρ(0) and R ~ 1.3 × 10−3 are the water density and the ratio of the air and water densities at the ocean surface. In Equation (6), the sign of the vertical component m(k, ω) of the wave vector should be chosen as follows: m = i|m|, when m 2 < 0; \( m=\left|m\right|\mathrm{s}\mathrm{g}\mathrm{n}\left({\omega}^2-\left(\gamma -1\right){g}^2{c}_s^{-2}\right), \) when m 2 ≥ 0 (Godin and Fuks 2012).

From Equations (1) and (4) to (6), we find the perturbation to the IGW dispersion relation (1) caused by the atmospheric loading:

$$ \delta k=\frac{Rk\left(2{k}^2{c}_s^2-\gamma {\omega}^2-2i{\omega}^2{c}_s^2m/g\right) \sinh 2kH}{2\left({k}^2{c}_s^2-{\omega}^2\right)\left( \sinh 2kH+2kH\right)}. $$
(7)

Here, only terms of the first order in R are retained, and ω and k satisfy the unperturbed dispersion relation (1). Equation (7) reduces to known results in the previously considered case of an infinitely deep ocean (Godin and Fuks 2012).

An equivalent of Equation (7) can be obtained, although with much more effort, without using the perturbation theory and by deriving instead an exact dispersion equation for acoustic-gravity waves in a coupled system consisting of a liquid layer and an isothermal gas half-space.

Equation (7) indicates that there are two distinct wave-propagation regimes, which are separated by a transition frequency such that m = 0. When m 2 < 0 in Equation (2), waves in the air are evanescent, their energy density decreasing with height z above the ocean surface as exp(−2z|m|). Then, correction δk (7) is real, and the IGW remains a true surface wave. When m 2 > 0, waves in air propagate away from the ocean surface and carry energy and momentum into the upper atmosphere. For propagating waves, the amplitude of the velocity oscillations increases with height as [ρ a (0)/ρ a (z)]1/2, where ρ a (z) is the air density and ρ a (0) refers to the mean sea level. Since the IGW phase speed is always smaller than sound speed in air, m 2 (2) can be positive only at frequencies ω smaller than the buoyancy frequency N = (γ − 1)1/2 g/c s . At ω < N, it follows from the radiation condition at z → +∞ that m < 0 when m 2 > 0 (Godin and Fuks 2012). Then, δk (7) is complex, and the imaginary part Im δk > 0. This means that, due to the atmospheric loading, the IGW slowly attenuates in the direction of its horizontal propagation and thus becomes a leaky wave (Brekhovskikh and Godin 1998). This is a natural corollary of the energy conservation and radiation of IGW energy into the atmosphere. Radiation into the atmosphere is only one of the mechanisms of IGW decay. As discussed in the section ‘IGW attenuation’ below, comparison with observed tsunami decay rates (Rabinovich et al. 2013) shows that the decay rate predicted by Equation (7) is smaller than measured decay rates.

Radiation of atmospheric waves

An analysis of Equations (1) and (2) shows that, within our simple environmental model, radiation of IGW energy into the atmosphere occurs at frequencies 0 < ω < ω b (H), where the transition frequency ω b steadily decreases with depth from ω b (0) = N in shallow water (about 3.0 mHz for c s  = 330 m/s) to zero at \( H={H}_b\equiv 4\left(\gamma -1\right){\gamma}^{-2}{c}_s^2/g\approx 0.82{c}_s^2/g \) (about 9,100 m for c s  = 330 m/s; see Figure 1a). When a more realistic environmental model is used, ω b (H) and H b also depend on the temperature and wind-velocity profiles in the atmosphere. For most of the ocean, the transition frequency is between 2.5 and 3 mHz (Figure 1a). Only small areas of the ocean with the greatest depths, such as the Mariana, Kermadec, and Puerto Rico trenches, have transition frequencies below 2.25 mHz.

Figure 1
figure 1

Geography of atmospheric wave radiation by infragravity waves in the ocean. (a) The transition frequency ω b from surface to leaky infragravity waves (in mHz) is shown by color for various geographical locations. At frequencies below the transition frequency, infragravity waves continuously radiate acoustic-gravity waves into the upper atmosphere. (b) The surface density I of the vertical flux of mechanical power radiated into the atmosphere (in W/m2) is shown by color on the logarithmic scale. The transition frequency and the power flux density are calculated for an isothermal atmosphere with sound speed c s  = 330 m/s.

Prediction of the upper cutoff frequency ω b is consistent with an earlier analysis (Peltier and Hines 1976) of the limiting case of non-dispersive tsunami waves, where kH 1. Indirect experimental support for our predictions is provided by the observations at the Infrasound Laboratory on the Island of Hawaii that the noise level, i.e., the spectral level of ambient pressure fluctuations in air, increases sharply at frequencies below approximately 4 mHz (Garces et al. 2011).

When investigating the conditions for IGW radiation into the atmosphere, it is imperative to consider the atmospheric waves as acoustic-gravity waves, i.e., to take into account effects due to buoyancy and compressibility simultaneously. When surface waves are considered as a potential motion of an incompressible fluid, which is an excellent approximation for wind seas (Gill 1982), orbital velocities in the wave exponentially decrease with distance from the sea surface at the same rate in the water and the air, and no radiation of wave energy into the atmosphere is predicted. In the opposite case, when buoyancy is neglected, we have acoustic waves in air, and radiation into the atmosphere occurs when the phase speed c of surface waves exceeds the sound speed (Brekhovskikh and Godin 1998). That would require the ocean depth \( H>{c}_s^2/g \) (Eckart 1953). Perhaps surprisingly, an account of the buoyancy effects results in converse dependence on the ocean depth. As discussed above, there is no radiation into the atmosphere from a horizontally homogeneous ocean with depth \( H>{H}_b\approx 0.82{c}_s^2/g. \) The reason for the qualitatively different behavior of acoustic and acoustic-gravity waves lies in the fact that m 2 decreases with increasing k 2, when g = 0, but increases when g ≠ 0 and ω < N, see Equation (2).

Observations of ionospheric manifestations of recent catastrophic tsunamis (Artru et al. 2005; Makela et al. 2011; Galvan et al. 2012) serve as a spectacular manifestation of radiation of acoustic-gravity atmospheric waves by IGWs in the ocean. Background IGWs, which are believed to be generated by the nonlinear interaction of sea swell with bathymetry close to shore, have much smaller amplitude than destructive tsunamis but are continuously present throughout the ocean (Webb et al. 1991; Herbers et al. 1995; Aucan and Ardhuin 2013; Godin et al. 2014). Over extended time periods, the contribution of the background IGWs into the flux of mechanical energy from the ocean into the upper atmosphere is expected to dominate over the contributions of tsunamis. For a rough estimate, we assume that powerful tsunamis with ocean surface displacement of 0.2 m account for half of the tsunami energy, occur on average every 5 years, last 2 h at open ocean locations, and affect 10% of the ocean surface. For RMS surface elevations in background IGWs, we conservatively take 3 mm (Webb et al. 1991; Aucan and Ardhuin 2013; Godin et al. 2013). Then, the time-averaged tsunami energy is about 2% of the background IGW energy.

The mechanical power radiated into the atmosphere can be calculated as an integral over the vertical component of the power flux density on the ocean surface (Godin and Fuks 2012). Random IGWs with power spectrum S of the ocean surface displacement radiate into the upper atmosphere the power

$$ I={\rho}_a(0){\displaystyle \underset{\omega_{\min }}{\overset{\omega_b(H)}{\int }}}\frac{\omega^3S\left(\omega, H\right)\left|m\left(\omega \right)\right|d\omega }{m^2\left(\omega \right)-{\omega}^2{c}_s^{-2}} $$
(8)

per unit area of the ocean surface. Here ω min is the lowest frequency of IGWs present and ρ a (0) is the air density at sea level; I = 0, when ω min ≥ ω b (H). When AGWs dissipate in the upper atmosphere, they exert a force on the mean flow (wind), which can be calculated as the vertical flux of the AGW pseudomomentum (Fritts and Alexander 2003). AGW forcing plays an important role in large-scale circulation of the middle and upper atmosphere (Fritts and Alexander 2003; Geller et al. 2013). The horizontal force per unit area of the ocean surface is given by the same integral as in Equation (8) but with ω 2k(ω)〉 substituted for ω 3 in the numerator in the integrand. Here, 〈k(ω)〉 is the horizontal wave vector of IGWs averaged over all propagation directions. The horizontal force is non-zero only for anisotropic IGW fields. The magnitude of the horizontal force reaches its upper bound, when all IGWs propagate in the same direction; then, |〈k(ω)〉| equals the wavenumber k(ω, H) of IGWs from the dispersion equation (1). This limiting case can be thought of as corresponding to IGWs generated near and propagating away from a straight shoreline. A force of the same order of magnitude in the direction perpendicular to the shoreline is expected to be applied to the upper atmosphere above the shore due to AGW radiation by IGWs with generic directivity. In the latter case, the requisite anisotropy arises from the absence of IGWs on the shore side of the coastline.

We apply the recently developed semi-empirical model S(ω, H) = B tanh−3/2 kH, where B = 2.1 × 10−5 m2/Hz, of the IGW power spectrum (Godin et al. 2013) to estimate the mechanical energy and momentum radiated by the ocean into the upper atmosphere. The spectrum was obtained for IGWs with frequencies 0.3 mHz and higher. The efficiency of AGW radiation decreases rapidly at frequencies below approximately 0.5 mHz (Figure 2a,b), which makes the surface densities of radiated power [Equation (8)] and pseudomomentum insensitive to unknown IGW spectral levels below 0.3 mHz. The IGW spectrum and the IGW and AGW dispersion relations (1) to (7) combine to make AGW radiation more efficient at shallower ocean depths and at frequencies close to but smaller than the transition frequency ω b (Figure 2a,b). The effect of the ocean depth is particularly pronounced for the momentum flux. Atmospheric waves have the unusual property that the horizontal and vertical components of their group velocity depend on frequency in rather different ways. The horizontal component is insensitive to frequency. In contrast, the vertical component (Figure 2c) and the direction of wave propagation (Figure 2d) vary strongly with frequency. The maximum vertical velocity exceeds 80 m/s and is attained at frequencies close to 0.8 ω b (Figure 2c). Geographical distribution of the mechanical energy radiation is illustrated by Figure 1b. The total power radiated into the middle and upper atmosphere by IGWs globally is estimated to be approximately 0.3 MW. The sound speed value c s  = 330 m/s, which is used in our model calculations, corresponds to air temperature of about 271 K. The modeling results presented in Figures 1, 2, 3, and 4 are insensitive to the particular choice of c s within the range of possible sound speed values above the ocean surface.

Figure 2
figure 2

Spectra of atmospheric waves radiated by infragravity waves in the ocean. (a) Power spectrum of the surface density of the mechanical energy flux is shown on the logarithmic scale relative 1 W/m2 Hz. (b) Power spectrum of the vertical flux of horizontal pseudomomentum is shown for the straight coast model on the logarithmic scale relative 1 Pa/Hz. AGW power spectra are calculated for various ocean depths using a semi-empirical model (Godin et al. 2013) of the IGW power spectrum. (c) Vertical component of the group velocity of AGWs in m/s. (d) Angle (in degrees) between the ocean surface and the direction of AGW propagation.

Figure 3
figure 3

Attenuation of infragravity waves in the ocean due to radiation of atmospheric waves. Decay time T of infragravity waves is shown as a function of wave frequency for ocean depths H = 2,000, 4,000, and 6,000 m. IGW attenuation is calculated assuming an isothermal atmosphere with sound speed c s  = 330 m/s.

Figure 4
figure 4

Measured power spectra. Measured power spectra of velocities’ fluctuations along the magnetic field line at an altitude of 150 km are compared with predicted contributions due to atmospheric waves radiated by infragravity waves in the ocean. The spectra are measured by radio sounding of the ionosphere using an onshore Dynasonde (Wright and Pitteway 1999; Zabotin et al. 2006) HF radar system located at NASA’s Wallops Island Flight Facility, Virginia, USA. The five panels refer to five distinct seasons in 2013 to 2014, and each panel contains four experimental spectral estimates obtained from the data series, specific dates and times of which are listed in the text. Blue dots show the results of spectral measurements. Red curves represent the results of smoothing the experimental spectra using a 0.5-mHz sliding window. Black curves show contributions of the atmospheric waves of oceanic origin, which are estimated for various ocean depths (200, 1,000, and 5,000 m, top to bottom) using a semi-empirical model (Godin et al. 2013) of the IGW spectra.

IGW attenuation

Atmospheric loading results in the perturbation δk, which is given by Equation (7), to the dispersion relation (1) of IGWs in an ocean with a free surface. From Equation (7) it follows that, at frequencies 0 < ω < ω b (H), IGW amplitude attenuates with propagation distance L as exp[−LIm(δk)] due to continuous radiation of mechanical energy into the atmosphere. Since R 1, an IGW needs to propagate a long distance of the order of R −1 of its wavelengths to radiate half of its energy into the atmosphere. As a function of time t, the energy E of free IGWs decreases exponentially: E(t) = E(0)exp(−t/T), where the e-folding time T = [2c g Im(δk)]−1 is usually referred to as the decay time. From Equations (1) and (7), for the IGW group speed c g and decay time, we find

$$ {c}_g={\left(\frac{\partial \omega }{\partial k}\right)}_H=\frac{\omega }{k}\cdot \frac{2kH+ \sinh 2kH}{2 \sinh 2kH} $$
(9)

and

$$ T=\frac{g\left({k}^2-{\omega}^2{c}_s^{-2}\right)}{R{\omega}^3\left|\mathrm{R}\mathrm{e}(m)\right|}. $$
(10)

Here, k is the real-valued IGW wavenumber; small relative corrections of the order of R are neglected in Equations (7) and (10). IGW attenuation due to AGW radiation is small over the wave period, and the decay time is usually large since the air-to-water density ratio R 1 (Figure 3). As expected, the decay time tends to infinity at the transition frequency ω b , where m = 0.

In addition to radiation of atmospheric waves, which takes place only at frequencies ω < ω b (H), many other processes (viscous attenuation, bottom friction and percolation, scattering at small-scale bathymetric features, nonlinear interaction with other ocean wave types) can potentially contribute to IGW attenuation. Energy loss at reflection from coastlines is usually considered to be the dominant decay mechanism, at least for long-period IGWs such as tsunamis (Munk et al. 1962; Munk 1963; Van Dorn 1984, 1987). Recent measurements by Rabinovich et al. (2013) give decay rates of 18.5 to 20.0 h and 20.5 to 21.5 h for the 2010 Chile tsunami at wave periods 6 to 20 min (wave frequency 0.8 to 2.8 mHz) and 20 to 60 min (wave frequency 0.3 to 0.8 mHz), respectively. For the 2011 Tohoku tsunami, the measured decay rates (Rabinovich et al. 2013) are 17.0 to 19.0 h and 21.3 to 25.0 h at wave periods 6 to 20 min and 20 to 60 min, respectively. Similar decay rates were derived previously by Munk (1963) and Van Dorn (1984) from much scarcer data sets than those employed in the recent study by Rabinovich et al. (2013).

Comparison of the predicted contribution of the atmospheric wave radiation (Figure 3) and measured (Munk 1963; Van Dorn 1984; Rabinovich et al. 2013) IGW decay rates shows that (i) IGW energy loss due to AGW radiation is consistent with available experimental data; (ii) radiation of atmospheric waves is a non-negligible mechanism of attenuation of IGWs away from shore; and (iii) depending on the average ocean depth, AGW radiation can account for up to one third of the observed decay of IGWs with periods from 6 to 20 min.

Experiment

Methods

We use Dynasonde techniques to measure spectra of vertical velocity disturbances in the upper atmosphere. Dynasonde is a method of precision ionospheric radio sounding, which is based on detection of radio echoes, i.e., radio wave reflections from the inhomogeneous ionospheric plasma, in a broad frequency range of 1 to 20 MHz. Various analysis routines specific for Dynasonde applications perform inversions of measured parameters (range, angles of arrival, polarization, line-of-sight Doppler speed, amplitude) of the echoes to retrieve a number of ionospheric parameters (vertical profile of the plasma density, vector velocities of the plasma structures, spectra of small-scale irregularities, etc.) attributed to real altitudes (Wright and Pitteway 1999; Zabotin et al. 2006). The transition from the echo ranges measured in the so-called virtual space (where the difference between the group speed of the high-frequency (HF) electromagnetic waves in plasma and the vacuum speed of light is neglected) to actual spatial coordinates is essential for accuracy of our analysis. A component of Dynasonde data analysis software, the 3-D plasma density inversion procedure NeXtYZ (Zabotin et al. 2006), finds 3-D spatial location of the reflection point of each echo. The inversion procedure possesses the accuracy and sensitivity, which are necessary to visualize and measure wave activity in the upper atmosphere. Averaging of the vertical component of the Doppler velocities for all echoes reflected within a narrow altitude range during a sounding session produces a single data point for our analysis.

We have used data obtained with the Wallops Island VIPIR/Dynasonde HF radar system in 2013 to 2014 in the regime where the sounding sessions were repeated every 2 min. The data screening has been performed using web-based visualization tools both provided by the standard Dynasonde software suite (Zabotin et al. 2005) and developed specifically for this project. Figure 5 presents an example of the standard output of the visualization software: a day-long scan of the vertical Doppler speed values shown by color as a function of the time of the day and of the real altitude. One can see a typical transition from night to day conditions near 13 UT and a transition from day to night conditions near 21 UT. The slightly inclined color strips are indications of the traveling ionospheric disturbances with their phase fronts propagating downward and therefore most likely caused by the acoustic-gravity waves propagating upward. Several criteria were used for the data selection. First of all, errors in ionogram trace identification and screening made by autonomous analysis must be infrequent in the data sets to not influence the statistics. Second, data series used for the spectral calculations must be long enough (at least several hours) because the frequencies of the wave disturbances of interest are as low as 0.1 mHz. We have adopted the coupling model between ionized and neutral components of the atmosphere according to which the velocity of the ion motion is estimated as a projection of the AGW-induced velocity perturbations onto the magnetic field line direction (Nicolls et al. 2014). This means that our analysis requires data from the upper ionosphere (above the E region). On the other hand, we cannot go too high because we would like to minimize the uncertainties and complications associated with the AGW absorption due to air viscosity and thermal conductivity. Only daytime ionospheric sounding data can provide information about the desired 130- to 180-km altitude range. Ionospheric reflections from the altitudes below 150 km are often blanketed by the underlying regular and sporadic E layers. The gray lacuna in Figure 5 between approximately 14 and 20 UT and between roughly 110 and 140 km is an example of manifestation of such blanketing. One can tolerate infrequently missing data points due to, for example, occasional high radio interference at some narrow frequency bands or due to short-lasting sporadic layers. But long interruptions in data series should be avoided. That is why summer seasons, when distortions to the data caused by sporadic E layers are common, are not suitable for this analysis.

Figure 5
figure 5

An example of the standard output of the Dynasonde analysis visualization software. A day-long scan of the vertical Doppler speed values is shown by color as a function of the time of the day (in UT hours) and of the real altitude (in km) for Wallops Island, Virginia, USA on 15 December 2013.

Combination of all the factors stated above makes 150 km the optimal altitude for data collection. Specifically, the following 20 time intervals, each one between 6 and 11 h long, were selected: 12 to 23 UT on 8 May 2013, 11 to 21 UT on 9 May 2013, 12 to 22 UT on 14 May 2013, 11 to 23 UT on 18 May 2013, 12 to 22 UT on 29 September 2013, 12 to 21 UT on 24 October 2013, 14 to 20 UT on 17 November 2013, 13 to 22 UT on 30 November 2013, 13 to 21 UT on 15 December 2013, 13 to 21 UT on 22 December 2013, 13 to 21 UT on 27 December 2013, 13 to 21 UT on 29 December 2013, 12 to 22 UT on 10 January 2014, 12 to 22 UT on 15 March 2014, 12 to 22 UT on 18 March 2014, 12 to 22 UT on 22 April 2014, 13 to 20 UT on 26 April 2014,12 to 21 UT on 26 September 2014, 13 to 22 UT on 29 September 2014, and 12 to 22 UT on 30 September 2014.

Every data point was obtained by averaging the vertical projections of the line-of-sight Doppler speed data for the 150- to 151-km altitude range. All the time intervals are affected by the diurnal variations of the Doppler speed; there is a background motion of the ionization layers downward after sunrise and upward before sunset. To suppress the spectral harmonics related to the diurnal variations, the Hann window was used for every data segment, with subsequent application of restoration coefficient 2 to the resulting power spectra. Some of the data points are missing, because of the reasons stated above, so the data series are inherently non-uniform. To account for this, the spectrum evaluation was performed with Hocke’s (1998) modification of the Lomb-Scargle algorithm.

The experimental Doppler speed spectra were converted into the spectra of the ion movements by applying the factor sin−2 φ, where φ is the magnetic dip angle. These are compared to the power spectra of the IGW-induced velocities parallel to the magnetic field line at an altitude of 150 km, which are calculated using the semi-empirical model spectra S(ω, H) (Godin et al. 2013) of the ocean surface displacements caused by the infragravity waves. The altitude chosen for the comparison is low enough (Hines 1960) to neglect effects of absorption of IGW-radiated acoustic-gravity waves, which have frequencies below approximately 3 mHz and horizontal spatial scales of tens and hundreds of kilometers. From AGW polarization relations (e.g., Godin 2012), we find the relation

$$ \left\langle {v}_m^2\right\rangle =\left\langle {v}_3^2\right\rangle \left[{ \sin}^2\varphi +\frac{\left(\gamma -1\right){g}^2{\omega}^{-2}-{c}_s^2}{2\left({c}_s^2-{\omega}^2{k}^{-2}\right)}{ \cos}^2\varphi \right] $$
(11)

between the wave-induced fluctuations v m of the velocity along the magnetic line and the vertical velocity v 3 in a stationary atmosphere. Here, angular brackets denote statistical averaging. The first and second terms in the brackets in Equation (11) come from the wave-induced vertical and horizontal motions, respectively. We assume here that the IGW field is isotropic. The power spectrum S v (ω, z) of the vertical velocity fluctuations v 3 at height z in the atmosphere is related to the spectrum of the ocean surface elevations as follows: S v (ω, z) = ω 2[ρ a (0)/ρ a (z)]S(ω, H). The power spectra of the IGW-induced fluctuations v m of the velocity along the magnetic line, which are calculated from S v (ω, z) using Equation (11), are utilized below for comparisons with experimental v m power spectra that are derived from Doppler speed measurements (Figure 4). Note that in a vertically stratified atmosphere with the horizontal winds, the total vertical power flux carried by an acoustic-gravity wave remains unchanged until the wave reaches a turning point, and existence and altitude of the turning points depend on the direction of wave propagation in the horizontal plane (Fritts and Alexander 2003). Therefore, the details of the vertical stratification are much less significant for the integral quantities considered in this paper than, for example, for calculating the spatial distribution of the wave field.

Observations

Atmospheric waves are known to be generated by a number of physical processes, associated most notably with convection, wind shear, and air-flow interaction with orography (Fritts and Alexander 2003; Geller et al. 2013). To assess the significance of the IGW contribution to the atmospheric wave field in the middle and upper atmosphere above oceans, we compare statistics of measured fluctuations of velocity with the velocity fluctuations due to random AGWs radiated by background IGWs in the ocean (Figure 4). To our knowledge, detailed measurements of velocity fluctuations in the middle and/or upper atmosphere above the ocean are not available at this time. As a proxy for the atmosphere above the ocean, we use experimental results obtained with HF radar systems at a coastal location in the North Atlantic.

Specifically, we use the data obtained with the VIPIR/Dynasonde HF radar system located at NASA’s Wallops Flight Facility in Virginia, USA. Projections of the velocity fluctuations onto the magnetic field line direction at ionospheric altitudes and their power spectra, which are shown in Figure 4, are retrieved from Doppler shifts of radio echoes reflected by the ionosphere in the real altitude range of 150 to 151 km (Wright and Pitteway 1999). The choice of this altitude range is determined by specifics of the ionospheric radio sounding as explained in the previous section. Figure 2d shows that AGWs of oceanic origin propagate obliquely away from the ocean surface. AGWs observed by the radar would be radiated by free IGWs at different distances from shore depending on the wave frequency and bathymetry. For comparison with the measured spectra, in Figure 4, we present predictions for oceans with three different depths H, which encompass the range of the relevant H values in the vicinity of Wallops. The IGW contribution to velocity fluctuations in the atmosphere is expected to be a weighted average of spectra corresponding to different H.

Although measured velocity fluctuations include contributions of both atmospheric waves of various origins and turbulence, Figure 4 shows that the expected velocity fluctuations due to the waves of oceanic origin are within the range covered by the natural variability of the measured spectra. Comparison of the smoothed measured spectra (red lines in Figure 4) and theoretical predictions (black lines in Figure 4) indicates that there are extended periods of time when the bulk of the velocity variance at the thermospheric altitudes observed at Wallops can be attributed to AGWs of oceanic origin. These observations suggest the feasibility of direct measurements that would confirm correlation (coherence) between the wave phenomena in both media. A project of this kind is underway with a support from the Basic Research Challenge program of the US Office of Naval Research (Zabotin et al. 2013).

From tropospheric to mesospheric heights, wave processes are believed to be dominated by AGWs that are generated by moist convection, shear flow instabilities, jets, atmospheric fronts, and interaction of winds with the earth’s topography (Gossard and Hooke 1975; Fritts and Alexander 2003; Geller et al. 2013). The power flux in these AGWs rapidly decreases with height due to viscous dissipation, wave breaking, and blocking at the critical levels (Gossard and Hooke 1975; Fritts and Alexander 2003; Geller et al. 2013). Only a very small (and unknown) fraction of the energy of these AGWs reaches ionospheric heights. AGWs of oceanic origin, and especially the AGWs radiated by IGWs in deep water, are characterized by considerably higher horizontal phase speeds, which can exceed the maximum wind speeds. This leads to larger horizontal spatial scales and, correspondingly, to weaker viscous absorption (Godin 2014). Fast AGWs do not have critical levels and escape the critical level filtering (Hines 1960; Gossard and Hooke 1975). Therefore, it should not be surprising that the AGWs of oceanic origin, which have small amplitudes and are hardly observable in the troposphere and stratosphere, play a much larger role at ionospheric heights.

Conclusion

Atmospheric loading leads to a qualitative change in properties of infragravity waves in the ocean. These waves cease to be surface ones and become leaky waves at frequencies below a certain depth-dependent transition frequency of the order of 3 mHz. The leaky waves continuously radiate a small fraction of their energy into the atmosphere as acoustic-gravity waves. The acoustic-gravity waves of oceanic origin propagate into the middle and upper atmosphere, where they appear to account for a significant part of the observed velocity fluctuations. As the atmospheric waves dissipate at high altitudes, they impact the atmospheric dynamics through deposition of mechanical energy and momentum transported from the ocean. Further research is needed to refine the simple model presented in this paper by account of geographic and seasonal variability of IGWs in the ocean and three-dimensional inhomogeneity of the atmosphere, in order to quantify the energy and momentum transport due to radiation of long atmospheric waves by background infragravity waves.

Abbreviations

AGW:

acoustic-gravity wave

IGW:

infragravity wave

UT:

universal time

References

  • Artru J, Ducic V, Kanamori H, Lognonné P, Murakami M (2005) Ionospheric detection of gravity waves induced by tsunamis. Geophys J Int 160:840–848

    Article  Google Scholar 

  • Aucan J, Ardhuin F (2013) Infragravity waves in the deep ocean: an upward revision. Geophys Res Lett 40:3435–3439

    Article  Google Scholar 

  • Brekhovskikh LM, Godin OA (1998) Acoustics of layered media. 1: plane and quasi-plane waves. Springer, Heidelberg

    Google Scholar 

  • Bromirski PD, Stephen RA (2012) Response of the Ross Ice Shelf, Antarctica, to ocean gravity-wave forcing. Ann Glaciol 53:163–172

    Article  Google Scholar 

  • Crawford WC, Webb SC, Hildebrand JA (1998) Estimating shear velocities in the oceanic crust from compliance measurements by two-dimensional finite difference modeling. J Geophys Res 103:9895–9916

    Article  Google Scholar 

  • Eckart C (1953) The theory of noise in continuous media. J Acoust Soc Am 25:195–199

    Article  Google Scholar 

  • Fritts DC, Alexander MJ (2003) Gravity wave dynamics and effects in the middle atmosphere. Rev Geophys 41:1003

    Article  Google Scholar 

  • Galvan DA, Komjathy A, Hickey MP, Stephens P, Snively J, Song YT, Butala MD, Mannucci AJ (2012) Ionospheric signatures of Tohoku-Oki tsunami of March 11, 2011: model comparisons near the epicenter. Radio Sci 47:RS4003

    Article  Google Scholar 

  • Garces MA, Perttu A, Williams B (2011) Infrasonic troposphere-ionosphere coupling in Hawaii, Paper SA23A-1897 presented at the 2011 AGU Fall Meeting, San Francisco, CA, USA, 5–9 December 2011

    Google Scholar 

  • Geller MA, Alexander MJ, Love PT, Bacmeister J, Ern M, Hertzog A, Manzini E, Preusse P, Sato K, Scaife AA, Zhou T (2013) A comparison between gravity wave momentum fluxes in observations and climate models. J Clim 26:6383–6405

    Article  Google Scholar 

  • Gill AE (1982) Atmosphere–ocean dynamics. Academic, New York

    Google Scholar 

  • Godin OA (2012) Acoustic-gravity waves in atmospheric and oceanic waveguides. J Acoust Soc Am 132:657–669

    Article  Google Scholar 

  • Godin OA (2014) Dissipation of acoustic-gravity waves: an asymptotic approach. J Acoust Soc Am 136:EL411–EL417

    Article  Google Scholar 

  • Godin OA, Fuks IM (2012) Transmission of acoustic-gravity waves through gas–liquid interfaces. J Fluid Mech 709:313–340

    Article  Google Scholar 

  • Godin OA, Zabotin NA, Sheehan AF, Yang Z, Collins JA (2013) Power spectra of infragravity waves in a deep ocean. Geophys Res Lett 40:2159–2165

    Article  Google Scholar 

  • Godin OA, Zabotin NA, Sheehan AF, Collins JA (2014) Interferometry of infragravity waves off New Zealand. J Geophys Res Oceans 118:1103–1122

    Article  Google Scholar 

  • Gossard EE, Hooke WH (1975) Waves in the atmosphere. Elsevier, Amsterdam

    Google Scholar 

  • Herbers THC, Elgar S, Guza RT (1995) Generation and propagation of infragravity waves. J Geophys Res 100:24863–24872

    Article  Google Scholar 

  • Hines CO (1960) Internal atmospheric gravity waves at ionospheric heights. Can J Phys 38:1441–1481

    Article  Google Scholar 

  • Hines CO (1972) Gravity waves in the atmosphere. Nature 239:73–78

    Article  Google Scholar 

  • Hocke K (1998) Phase estimation with the Lomb-Scargle periodogram method. Ann Geophys 16:356–358

    Google Scholar 

  • Komjathy A, Galvan DA, Stephens P, Butala MD, Akopian V, Wilson VO, Mannucci AJ, Hickey M (2012) Detecting ionospheric TEC perturbations caused by natural hazards using a global network of GPS receivers: the Tohoku case study. Earth Planets Space 64:1287–1294

    Article  Google Scholar 

  • Livneh DJ, Seker I, Djuth FT, Mathews JD (2007) Continuous quasiperiodic thermospheric waves over Arecibo. J Geophys Res 112, A07313

    Article  Google Scholar 

  • Makela JJ, Lognonné P, Hébert H, Gehrels T, Rolland L, Allgeyer S, Kherani A, Occhipinti G, Astafyeva E, Coïsson P, Loevenbruck A, Clévédé E, Kelley MC, Lamouroux J (2011) Imaging and modeling the ionospheric airglow response over Hawaii to the tsunami generated by the Tohoku earthquake of 11 March 2011. Geophys Res Lett 38:L13305

    Article  Google Scholar 

  • Munk WH (1963) Some comments regarding diffusion and absorption of tsunamis. In: Proceedings of the Tsunami Meeting, X Pacific Science Congress, IUGG Monographs, 24th edn. IUGG, Paris, pp 53–72

    Google Scholar 

  • Munk WH, Miller GR, Snodgrass FE (1962) Long-period waves over California’s continental borderland, III. The decay of tsunamis and the dissipation of tides. J Mar Res 20:119–120

    Google Scholar 

  • Nicolls MJ, Vadas SL, Aponte N, Sulzer MP (2014) Horizontal parameters of daytime thermospheric gravity waves and E region neutral winds over Puerto Rico. J Geophys Res 119:575–600

    Article  Google Scholar 

  • Occhipinti G, Rolland L, Lognonné P, Watada S (2013) From Sumatra 2004 to Tohoku‐Oki 2011: the systematic GPS detection of the ionospheric signature induced by tsunamigenic earthquakes. J Geophys Res 119:3626–3636

    Article  Google Scholar 

  • Peltier WR, Hines CO (1976) On the possible detection of tsunamis by monitoring of the ionosphere. J Geophys Res 81:1995–2000

    Article  Google Scholar 

  • Rabinovich AB, Candella RN, Thomson RE (2013) The open ocean energy decay of three recent trans-Pacific tsunamis. Geophys Res Lett 40:3157–3162

    Article  Google Scholar 

  • Rhie J, Romanowicz B (2004) Excitation of Earth’s continuous free oscillations by atmosphere–ocean-seafloor coupling. Nature 431:552–556

    Article  Google Scholar 

  • Tsai VC, Ampuero J-P, Kanamori H, Stevenson DJ (2013) Estimating the effect of Earth elasticity and variable water density on tsunami speeds. Geophys Res Lett 40:492–496

    Article  Google Scholar 

  • Van Dorn WG (1984) Some tsunami characteristics deducible from tide records. J Phys Oceanogr 14:353–363

    Article  Google Scholar 

  • Van Dorn WG (1987) Tide gage response to tsunamis. Part II: other oceans and smaller seas. J Phys Oceanogr 17:1507–1516

    Article  Google Scholar 

  • Watada S, Kusumoto S, Satake K (2014) Traveltime delay and initial phase reversal of distant tsunamis coupled with the self-gravitating elastic Earth. J Geophys Res 119:4287–4310

    Article  Google Scholar 

  • Webb SC, Zhang X, Crawford W (1991) Infragravity waves in the deep ocean. J Geophys Res 96:2723–2736

    Article  Google Scholar 

  • Wright JW, Pitteway MLV (1999) A new data acquisition concept for digital ionosondes: phase-based echo recognition and real-time parameter estimation. Radio Sci 34:871–882

    Article  Google Scholar 

  • Zabotin NA, Wright JW, Bullett TW, Ye Zabotina L (2005) Dynasonde 21 principles of data processing, transmission, storage and web service. In: Proceedings of the ionospheric effects symposium 2005. Alexandria, VA, USA., pp 7B3-1–7B3-3

    Google Scholar 

  • Zabotin NA, Wright JW, Zhbankov GA (2006) NeXtYZ: three-dimensional electron density inversion for dynasonde ionograms. Radio Sci 41:RS6S32

    Article  Google Scholar 

  • Zabotin NA, Godin OA, Bullet TW (2013) Tracing acoustic-gravity waves from the ocean into the ionosphere: ONR BRC Project., University of Colorado Boulder. http://surf.colorado.edu. Accessed 29 Dec 2014

    Google Scholar 

Download references

Acknowledgements

The authors are grateful to L. Zabotina for development of web-based visualization software for results of the Dynasonde data analysis. Insightful comments by Nozomu Nishitani and two anonymous referees helped to improve presentation. The work reported in this paper was supported by the Office of Naval Research through award N00014-13-1-0348. This project utilized the Janus supercomputer, which is supported by the National Science Foundation (award number CNS-0821794) and operated by the University of Colorado Boulder. The VIPIR/Dynasonde system is hosted by NASA’s Wallops Flight Facility.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Oleg A Godin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

OG and NZ conceived this project and wrote the paper. OG developed the analytical theory. NZ performed numerical calculations and Dynasonde data analysis. TB performed radar experiments and collected data. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Godin, O.A., Zabotin, N.A. & Bullett, T.W. Acoustic-gravity waves in the atmosphere generated by infragravity waves in the ocean. Earth Planets Space 67, 47 (2015). https://doi.org/10.1186/s40623-015-0212-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40623-015-0212-4

Keywords