Page updated: March 20, 2021
Author: Curtis Mobley
View PDF

# Wave Variance Spectra: Examples

This page gives two examples of wave elevation variance spectra. The Pierson-Moskowitz omnidirectional spectrum describes gravity waves in a “fully developed” sea. A fully developed sea is an idealization of the statistically steady-state wave ﬁeld resulting from a steady wind blowing for an inﬁnitely long time over an inﬁnite fetch. (In practice, the duration and fetch required to achieve something close to a fully developed sea depend on the wind speed. A steady wind of $5\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$ blowing for 10 hours over a fetch of 60 km might be adequate; for hurricane winds of $35\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$, a fetch of a few thousand kilometers and a duration of several days would be required. Thus it is much easier to approach a fully developed sea at low wind speeds than at high.)

The Elfouhaily et al. directional spectrum includes both gravity and capillary waves scales. Moreover, it has a parameter that describes the wave age, so than any sea state from young (the wind has just started blowing) to fully developed can be simulated. The Elfouhaily et al. spectrum will be used to generate the two-dimensional sea surface examples on the following pages.

### The Pierson-Moskowitz Omnidirectional Gravity Wave Spectrum

The omnidirectional Pierson and Moskowitz (1964) spectrum, formulated in terms of angular spatial frequency $k$, is

 ${\mathsc{𝒮}}_{PM}\left(k\right)=\frac{\alpha }{2{k}^{3}}exp\left[-\beta {\left(\frac{g}{k}\right)}^{2}\frac{1}{{U}_{19}^{4}}\right]\phantom{\rule{2em}{0ex}}\left[{m}^{2}∕\left(rad∕m\right)\right]\phantom{\rule{0.3em}{0ex}},$ (1)

where

$\alpha =0.0081$,
$\beta =0.74$,
$g=9.82\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-2}$ is the acceleration of gravity,
${U}_{19}$ is the wind speed in $m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$ at 19.5 m above the sea surface, and
$k$ is the angular spatial frequency in $rad\phantom{\rule{2.6108pt}{0ex}}{m}^{-1}$.

The wind speed at 19.5 m can be obtained from the more commonly used wind at 10 m above the sea surface by the approximate formula

 ${U}_{19}\approx 1.026{U}_{10}\phantom{\rule{0.3em}{0ex}}.$

As has already been noted, it is often of interest to express a variance spectrum in terms of other variables, such as the wavenumber $\nu$ or the temporal angular frequency $\omega$. To change variables in a spectral density function, the key is to recall that a variance density function by deﬁnition expresses the variance per unit frequency interval. The variance contained in some interval $dk$ of the spatial angular frequency equals the variance contained in the corresponding interval $d\nu$ of the wavenumber or the interval $d\omega$ of the temporal frequency. Thus we have

 ${\mathsc{𝒮}}_{PM}\left(k\right)dk={\mathsc{𝒮}}_{PM}\left(\nu \right)d\nu ={\mathsc{𝒮}}_{PM}\left(\omega \right)d\omega \phantom{\rule{0.3em}{0ex}}.$

To change the variable from $k=2\pi \nu$ to $\nu$, the previous equation gives

 ${\mathsc{𝒮}}_{PM}\left(\nu \right)={\mathsc{𝒮}}_{PM}\left(k\right)\frac{dk}{d\nu }={\mathsc{𝒮}}_{PM}\left(k=2\pi \nu \right)2\pi \phantom{\rule{0.3em}{0ex}},$

which results in

 ${\mathsc{𝒮}}_{PM}\left(\nu \right)=\frac{\alpha }{8{\pi }^{2}{\nu }^{3}}exp\left[-\beta {\left(\frac{g}{2\pi \nu }\right)}^{2}\frac{1}{{U}_{19}^{4}}\right]\phantom{\rule{2em}{0ex}}\left[{m}^{2}∕\left(1∕m\right)\right]\phantom{\rule{0.3em}{0ex}}.$ (2)

To change variables from spatial angular frequency $k$ to temporal angular frequency $\omega$, we use the dispersion relation for gravity waves in deep water,

 ${\omega }^{2}=gk\phantom{\rule{0.3em}{0ex}},$

to evaluate

 $\frac{dk}{d\omega }=\frac{2\omega }{g}\phantom{\rule{0.3em}{0ex}},$

 ${\mathsc{𝒮}}_{PM}\left(\omega \right)=\frac{\alpha {g}^{2}}{{\omega }^{5}}exp\left[-\beta {\left(\frac{g}{\omega {U}_{19}}\right)}^{4}\right]\phantom{\rule{2em}{0ex}}\left[{m}^{2}∕\left(rad∕s\right)\right]\phantom{\rule{0.3em}{0ex}}.$ (3)

All of these formulas have units of meters squared over the appropriate frequency. (The quantities $dk∕d\nu$ and $dk∕d\omega$ seen in the conversions are the Jacobians for the one-dimensional changes of variables.) Figure 1 shows the Pierson-Moskowitz spectrum as functions of $k$ and $\omega$ for wind speeds of ${U}_{10}=5,10$, and $15\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$. Figure 1: The Pierson-Moskowitz variance spectrum as functions of $k$ and $\omega$ for wind speeds of ${U}_{10}=5,10$, and $15\phantom{\rule{1em}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$. The vertical dotted lines at $k=370\phantom{\rule{2.6108pt}{0ex}}rad∕m$ and $\omega =60.3\phantom{\rule{2.6108pt}{0ex}}rad∕s$ show the boundary between gravity and capillary waves.

This spectrum has withstood the test of time fairly well (Alves and Banner (2003)) as a description of gravity waves in fully developed seas. However, it should not be used for high-frequency (short-wavelength) gravity waves, and certainly not for capillary waves. Likewise, it does not describe young seas, which have not had the time or fetch needed to approach the state of a well developed sea.

Figure 2 shows the Pierson-Moskowitz slope spectra for three wind speeds. Note that the slope spectrum falls oﬀ much more slowly for high frequencies than does the elevation spectrum. That means that the higher frequencies contribute much more to the total slope variance than they do to the total elevation variance. Figure 2: The Pierson-Moskowitz slope spectrum for wind speeds of ${U}_{10}=5,10,15\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$. The plot uses the same ordinate scale as used for the elevation spectrum in the left panel of Fig. 1 in order to highlight the slow falloﬀ of the slope spectrum compared to the elevation spectrum. The vertical dotted line is the boundary between gravity and capillary waves.

### The Elfouhaily et al. Directional Gravity-Capillary Wave Spectrum

In order to model two-dimensional sea surfaces $z\left(x,y\right)$, we need a 2-D wave elevation variance spectrum that describes the distribution of wave variance for waves propagating in diﬀerence directions. The numerical examples of 2-D sea surfaces to be seen on the following pages use the 2-D elevation variance spectrum of Elfouhaily et al. (1997), which is described here for later reference. For brevity, this paper and their model are denoted here by “ECKV,” taken from the initials of the authors’ last names.

The ECKV spectrum has an omnidirectional variance spectrum that explicitly includes both the gravity and capillary wave scales. The boundary between gravity and capillary waves is taken to be $k=\sqrt{\rho g∕\tau }=370\phantom{\rule{2.6108pt}{0ex}}rad∕m$, the angular spatial frequency at which the restoring forces (which tend to bring a wave surface back to an unperturbed level surface) of gravity and surface tension are equal. The corresponding wavelength is $\Lambda =2\pi ∕370=0.017\phantom{\rule{2.6108pt}{0ex}}m$. ECKV then combine their omnidirectional spectrum with a spreading function to obtain their one-sided, directional variance spectrum. Using their notation, the 2-D ECKV spectrum has the form

 $\Psi \left(k,\phi \right)=\frac{1}{k}\mathsc{𝒮}\left(k\right)\Phi \left(k,\phi \right)\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}45\right]\phantom{\rule{0.3em}{0ex}}.$ (4)

Here $\mathsc{𝒮}\left(k\right)$ is the 1-D omnidirectional spectrum with units of ${m}^{2}∕\left(rad∕m\right)$, and $\Phi \left(k,\phi \right)$ is a non-dimensional spreading function. $\Psi \left(k,\phi \right)$ thus has units of ${m}^{2}∕{\left(rad∕m\right)}^{2}$. Equation labels such as [ECKV 45] give for reference the corresponding equation in the ECKV paper.

The ECKV omnidirectional spectrum is

 $\mathsc{𝒮}\left(k\right)=\frac{{B}_{l}+{B}_{h}}{{k}^{3}}\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}30\right]\phantom{\rule{0.3em}{0ex}},$ (5)

where ${B}_{l}$ is the low-frequency (long gravity wave) contribution to the variance, and ${B}_{h}$ is the high-frequency (short gravity wave to capillary wave) contribution. (The quantity ${k}^{3}\mathsc{𝒮}\left(k\right)$ is called the curvature or saturation spectrum and is of interest in physical oceanography because it is related to the rate of variance dissipation of the waves. Thus ECKV refer to ${B}_{l}$ and ${B}_{h}$ as the low and high frequency curvature spectra.) The components of the omnidirectional spectrum are given by

$\begin{array}{llll}\hfill {L}_{PM}=& exp\left[-1.25{\left({k}_{p}∕k\right)}^{2}\right]\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}2\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill \Gamma =& exp\left\{-\frac{1}{2{\sigma }^{2}}\left[{\left(\sqrt{k∕{k}_{p}}-1\right)}^{2}\right]\right\}\phantom{\rule{2em}{0ex}}\left[below\phantom{\rule{1em}{0ex}}ECKV\phantom{\rule{1em}{0ex}}3\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {J}_{p}=& {\gamma }^{\Gamma }\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}3\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {F}_{p}=& {L}_{PM}{J}_{p}exp\left\{-0.3162{\Omega }_{c}\left(\sqrt{k∕{k}_{p}}-1\right)\right\}\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}32\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {F}_{m}=& {L}_{PM}{J}_{p}exp\left[-0.25{\left(k∕{k}_{m}-1\right)}^{2}\right)\right]\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}41\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill & \left(Note:A\phantom{\rule{2.6108pt}{0ex}}typo\phantom{\rule{2.6108pt}{0ex}}in\phantom{\rule{2.6108pt}{0ex}}ECKV\phantom{\rule{2.6108pt}{0ex}}Eq.\phantom{\rule{2.6108pt}{0ex}}41\phantom{\rule{2.6108pt}{0ex}}omitted\phantom{\rule{2.6108pt}{0ex}}the\phantom{\rule{2.6108pt}{0ex}}{L}_{PM}{J}_{p}\phantom{\rule{1em}{0ex}}factor\phantom{\rule{2.6108pt}{0ex}}in\phantom{\rule{2.6108pt}{0ex}}{F}_{m}\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {B}_{l}=& 0.5{\alpha }_{p}\left({c}_{p}∕c\right){F}_{p}\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}31\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {B}_{h}=& 0.5{\alpha }_{m}\left({c}_{m}∕c\right){F}_{m}\phantom{\rule{2em}{0ex}}\left[ECKV\phantom{\rule{1em}{0ex}}40\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$

where

$\alpha =0.0081$,
$\beta =1.25$,
$g=9.82\phantom{\rule{1em}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-2}$ is the acceleration of gravity,
${U}_{10}$ is the wind speed in $m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$ at 10 m above the sea surface
$k$ is the angular spatial frequency in $rad\phantom{\rule{2.6108pt}{0ex}}{m}^{-1}$
${\Omega }_{c}$ is deﬁnes the age of the waves for the given wind speed:
= 0.84 for a fully developed sea (corresponds to Pierson-Moskowitz)
= 1 for a “mature” sea [used in ECKV Fig 8a]
= 2 to 5 for a “young” sea; the maximum allowed value is 5
$C{d}_{10N}=0.00144$ is a drag coeﬃcient [value deduced from ECKV Fig 11]
${u}^{\ast }=\sqrt{C{d}_{10N}}{U}_{10}$ is the friction velocity [using ECKV 61]
${a}_{o}=0.1733$ (ECKV 59)
${a}_{p}=4.0$
${k}_{m}=370.0\phantom{\rule{2.6108pt}{0ex}}rad∕m$
${c}_{m}=0.23\phantom{\rule{2.6108pt}{0ex}}m∕s$ is the phase speed of the wave with spatial frequency ${k}_{m}$
${a}_{m}=0.13{u}^{\ast }∕{c}_{m}$ [ECKV 59]
$\gamma =1.7\phantom{\rule{1em}{0ex}}if\phantom{\rule{1em}{0ex}}{\Omega }_{c}\le 1\phantom{\rule{1em}{0ex}}else\phantom{\rule{1em}{0ex}}\gamma =1.7+6{log}_{10}\left({\Omega }_{c}\right)$
$\sigma =0.08\left(1+4{\Omega }_{c}^{-3}\right)$
${\alpha }_{p}=0.006{\Omega }_{c}^{0.55}$ [ECKV Eq. 34]
${\alpha }_{m}=0.01\left[1+ln\left({u}^{\ast }∕{c}_{m}\right)\right]\phantom{\rule{1em}{0ex}}if\phantom{\rule{1em}{0ex}}{u}^{\ast }\le {c}_{m}\phantom{\rule{1em}{0ex}}else\phantom{\rule{1em}{0ex}}{\alpha }_{m}=0.01\left[1+3ln\left({u}^{\ast }∕{c}_{m}\right)\right]$ [ECKV 44]
${k}_{o}=g∕{U}_{10}^{2}$
${k}_{p}={k}_{o}{\Omega }_{c}^{2}$ is the spatial frequency of the maximum of the spectrum
${c}_{p}=\sqrt{g∕{k}_{p}}$ is the phase speed of the wave with spatial frequency ${k}_{p}$
$c=\sqrt{\left(g∕k\right)\left(1+{\left(k∕{k}_{m}\right)}^{2}\right)}$ is the phase speed of the wave

At the lower frequencies, the ECKV spectrum is essentially the Pierson-Moskowitz spectrum (the ${L}_{PM}$ term above) with an enhancement (the ${J}_{p}$ term) that adds more energy to the lower frequencies. The highest frequencies have a cutoﬀ due to viscous damping of the smallest capillary waves. The ECKV omnidirectional elevation and slope spectra are illustrated in Fig. 3 for the case of a fully developed sea and three wind speeds. Figure 4 shows the spectra as a function of wave age for a wind speed of ${U}_{10}=10\phantom{\rule{1em}{0ex}}m\phantom{\rule{1em}{0ex}}{s}^{-1}$. Figure 3: The omnidirectional part $\mathsc{𝒮}$ of the ECKV elevation variance spectrum (left panel) and slope spectrum ${k}^{2}\mathsc{𝒮}$ (right panel) for fully developed seas and wind speeds of ${U}_{10}=5,10,15\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$. The gray lines show the corresponding Pierson-Moskowitz spectra from Fig. 2. Figure 4: The omnidirectional part of the ECKV elevation spectrum (left panel) and slope spectrum (right panel) for a wind speed of ${U}_{10}=10\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$ and wave ages from very young (${\Omega }_{c}=5$) to mature (${\Omega }_{c}=1$) to fully developed (${\Omega }_{c}=0.84$). Compare with Fig. 3.

The ECKV spreading function is given by

$\begin{array}{llll}\hfill \Phi \left(k,\phi \right)=& \frac{1}{2\pi }\left[1+\Delta \left(k\right)cos\left(2\phi \right)\right]\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill =& \frac{1}{2\pi }\left\{1+tanh\left[{a}_{o}+{a}_{p}{\left(c∕{c}_{p}\right)}^{2.5}+{a}_{m}{\left({c}_{m}∕c\right)}^{2.5}\right]cos\left(2\phi \right)\right\}\phantom{\rule{2em}{0ex}}& \hfill \text{(6)}\end{array}$

Note that this function is symmetric about $\phi =\pi ∕2$; i.e., the function has symmetric spreading downwind and upwind. This is consistent with a symmetric variance spectrum $\Psi \left(-k\right)=\Psi \left(k\right)$ as would be obtained from the Fourier transform of a snapshot of a sea surface. This symmetry will be explained on the following pages.

A commonly used family of alternate spreading functions is given by the “cosine-2s” functions of Longuet-Higgins et al. 1963, which have the form

 $\Phi \left(k,\phi \right)={C}_{s}{cos}^{2s}\left(\phi ∕2\right)\phantom{\rule{0.3em}{0ex}},$ (7)

where the normalizing coeﬃcient is

 ${C}_{s}=\frac{1}{2\sqrt{\pi }}\frac{\Gamma \left(s+1\right)}{\Gamma \left(s+1∕2\right)}\phantom{\rule{0.3em}{0ex}},$

and $s$ is a spreading parameter that in general depends on $k$, ${U}_{10}$, and wave age. In this equation $\Gamma$ is the customary gamma function deﬁned by $\Gamma \left(p\right)\equiv {\int }_{0}^{\infty }{x}^{p-1}{e}^{-x}dx$ where $p>0$. The cosine-2s functions are asymmetric, with much stronger downwind than upwind propagation.

The ECKV and cosine-2s spreading functions are illustrated in Fig. 5. Both of these functions satisfy the normalization condition (??). Both spreading functions transition from strongly forward peaked at low spatial frequencies (long gravity waves; the red curves) to curves with signiﬁcant propagation at right angles to the wind at high frequencies (capillary waves; the blue curves). The cosine-2s curves are asymmetric in $±k$ and have at least a small amount of upwind propagation at all frequencies (except at exactly upwind, $\phi =\pi$). Not surprisingly, the real ocean is more complicated than either of these models. In particular, observations of long-wave gravity waves tend to show a bimodal spreading about the downwind direction, which transitions to a more isotropic, unimodal spreading at shorter wavelengths (Heron (2006). However, the simple models of Eqs. (6) and (7) are adequate for the present purpose of illustrating surface-generation techniques. The eﬀect of the choice of spreading function on the generated waves will be illustrated on the Spreading Function Eﬀects page. Figure 5: Example spreading functions according to the ECKV model (left) and the cosine-2s model (right) for a wind speed of $10\phantom{\rule{2.6108pt}{0ex}}m\phantom{\rule{2.6108pt}{0ex}}{s}^{-1}$. Small $s$ values correspond to large spatial frequencies $k$. Downwind is to the right, upwind is to the left.

On the following pages we will learn the important distinction between “one-sided” or “folded” spectra and the associated “two-sided” spectra. The ECKV spectrum as given above is a one-sided spectrum.