<p>It is helpful to understand the meaning of Fourier coefficients in the directional wave spectrum analysis.</p>

<p>A pitch-and-roll buoy measures time series of water elevation $\eta$, and slopes in both Cartesian directions, $\eta_x$ and $\eta_y$. The Fourier components $a_0$, $a_1$, $b_1$, $a_2$, $b_2$, are related to cross-spectra of elevation and slope time series (introduced by Longuet-Higgins 1963):</p>

<p>$$
a_0 = \dfrac{C_{11}}{\pi}
$$</p>

<p>$$
a_1 = \dfrac{Q_{12}}{\pi}
$$</p>

<p>$$
b_1 = \dfrac{Q_{13}}{\pi}
$$</p>

<p>$$
a_2 = \dfrac{C_{22}-C_{33}}{k^2\pi}
$$</p>

<p>$$
b_2 = \dfrac{2C_{23}}{k^2\pi}
$$</p>

<p>where the cross-spectra are:</p>

<p>$$
C_{11}(f) = C[\eta(f)\eta(f)] = \int_0^{2\pi} E(f,\theta)d\theta
$$</p>

<p>$$
Q_{12}(f) = -Q[\eta(f)\eta_x(f)]k^{-1} = -\int_0^{2\pi} E(f,\theta)\cos{\theta}d\theta
$$</p>

<p>$$
Q_{13}(f) = -Q[\eta(f)\eta_y(f)]k^{-1} = -\int_0^{2\pi} E(f,\theta)\sin{\theta}d\theta
$$</p>

<p>$$
C_{22}(f) = C[\eta_x(f)\eta_x(f)]k^{-2} = \int_0^{2\pi} E(f,\theta)\cos^2{\theta} d\theta
$$</p>

<p>$$
C_{23}(f) = C[\eta_x(f)\eta_y(f)]k^{-2} = \int_0^{2\pi} E(f,\theta)\sin{\theta}\cos{\theta} d\theta
$$</p>

<p>$$
C_{33}(f) = C[\eta_y(f)\eta_y(f)]k^{-2} = \int_0^{2\pi} E(f,\theta)\sin^2{\theta} d\theta
$$</p>

<p>From here you can see that the mean direction, as defined in the paper that you cite, is:</p>

<p>$$
\theta_0 = \arctan\left({\dfrac{b_1}{a_1}}\right) = \arctan\left({\dfrac{Q_{13}}{Q_{12}}}\right) = \arctan\left({\dfrac{\int_0^{2\pi}E(f,\theta)\sin{\theta}d\theta}{\int_0^{2\pi}E(f,\theta)\cos{\theta}d\theta}}\right)
$$</p>

<p>Because $Q_{12}$ and $Q_{13}$ are proportional to the integrated energy in the zonal and meridional directions, respectively, the arctangent of their ratio yields mean wave direction.</p>

<p>The authors could use higher moments to calculate direction, but this would have different physical meaning. In this specific case, using $\theta = \arctan\left({b_2/a_2}\right)$ would yield peak (dominant) wave direction in the context of pitch-and-roll buoys. </p>

<p>You see that the above analysis is inherently limited by the quantities that are measured, specifically $\eta$, $\eta_x$, and $\eta_y$. If higher order quantities, say $\eta_{xx}$, $\eta_{xy}$, $\eta_{yy}$ were available, the directional spectrum could be described at a higher accuracy even using just Fourier decomposition.</p>

<p>Further reading:</p>

<ul>
<li><p><a href="http://journals.ametsoc.org/doi/abs/10.1175/1520-0485%281984%29014%3C1800%3AADAOWD%3E2.0.CO%3B2" rel="nofollow noreferrer">Oltman-Shay and Guza (1984)</a></p></li>
<li><p><a href="http://www.ndbc.noaa.gov/wavemeas.pdf" rel="nofollow noreferrer">Nondirectional and directional wave data analysis procedures at NDBC</a></p></li>
</ul>

CC BY-SA 3.0
2017-01-30T21:43:24.290Z
Thank you again - you answered another question of mine a few days ago, you're very helpful!
RH_data_maths
CC BY-SA 3.0
2017-01-30T20:25:11.230Z
2017-01-30T21:06:04.393Z
2017-01-30T21:06:04.393Z
milancurcic
<p>When modelling ocean waves, a directional distribution $D_f(\theta)$ is used together with a frequency spectrum $S(f)$ to describe the energy of waves at a particular frequency $f$ and angle $\theta$. </p>

<p>I understand that the directional distribution can be written as a Fourier series i.e.
$$D_f(\theta) = \frac{1}{2\pi}\left[ 1 + 2\sum_{n=1}^{\infty}\{a_n\cos(n\theta) + b_n\sin(n\theta) \} \right] $$
where $a_n = \int_0^{2\pi} D_f(\theta)\cos(\theta)\,d\theta$ and $b_n = \int_0^{2\pi} D_f(\theta)\sin(\theta)\,d\theta$.</p>

<p>In <a href="http://journals.ametsoc.org/doi/pdf/10.1175/1520-0485(1988)018%3C1020%3AAMFTRA%3E2.0.CO%3B2" rel="nofollow noreferrer">Kuik (1988)</a>, the mean wave direction, $\theta_0$, is found by calculating 
$$\theta_0 = \arctan\left(\frac{b_1}{a_1}\right)$$
where $b_1$ and $a_1$ are the first order Fourier coefficients. </p>

<p>Alongside this definition, the author refers the reader to Borgman (1969) but I can't find this paper on the web. </p>

<p>My question is why is it only the first order Fourier coefficients $a_1$ and $b_1$ used in this calculation?</p>

<p>EDIT: 
After giving this more thought, I think that the fact they are Fourier coefficients is somewhat of a coincidence. </p>

<p>If the directional distribution is seen as the PDF (as the integral of it is equal to 1) then $a_1$ and $b_1$ are more like the expected values that the cosine and sine of the angle $\theta$ take. The average values can then be used in the $atan2$ function to determine the mean angle.</p>

CC BY-SA 3.0
2017-01-30T18:24:43.250Z
I imagine it is a typo and they are referring to Borgman (1979) "Directional Wave Spectra from Wave Sensors" http://link.springer.com/chapter/10.1007/978-1-4684-3399-9_12 I think that article dealt with the formulation.
arkaia
<ocean><waves>
The mean direction of waves in a directional distribution
RH_data_maths
# The mean direction of waves in a directional distribution

(Asked by: username_0 on 2017-01-30)

When modelling ocean waves, a directional distribution $D\_f(\theta)$ is used together with a frequency spectrum $S(f)$ to describe the energy of waves at a particular frequency $f$ and angle $\theta$.

I understand that the directional distribution can be written as a Fourier series i.e.
$$D\_f(\theta) = \frac{1}{2\pi}\left[ 1 + 2\sum\_{n=1}^{\infty}\{a\_n\cos(n\theta) + b\_n\sin(n\theta) \} \right] $$
where $a\_n = \int\_0^{2\pi} D\_f(\theta)\cos(\theta)\,d\theta$ and $b\_n = \int\_0^{2\pi} D\_f(\theta)\sin(\theta)\,d\theta$.

In [Kuik (1988)](http://journals.ametsoc.org/doi/pdf/10.1175/1520-0485(1988)018%3C1020%3AAMFTRA%3E2.0.CO%3B2), the mean wave direction, $\theta\_0$, is found by calculating
$$\theta\_0 = \arctan\left(\frac{b\_1}{a\_1}\right)$$
where $b\_1$ and $a\_1$ are the first order Fourier coefficients.

Alongside this definition, the author refers the reader to Borgman (1969) but I can't find this paper on the web.

My question is why is it only the first order Fourier coefficients $a\_1$ and $b\_1$ used in this calculation?

EDIT:
After giving this more thought, I think that the fact they are Fourier coefficients is somewhat of a coincidence.

If the directional distribution is seen as the PDF (as the integral of it is equal to 1) then $a\_1$ and $b\_1$ are more like the expected values that the cosine and sine of the angle $\theta$ take. The average values can then be used in the $atan2$ function to determine the mean angle.

**Tags:** ocean, waves

> **Comments:**
> **username_1**: I imagine it is a typo and they are referring to Borgman (1979) "Directional Wave Spectra from Wave Sensors" http://link.springer.com/chapter/10.1007/978-1-4684-3399-9_12 I think that article dealt with the formulation.  

---

## 1 answers

### Answer by username_2 on 2017-01-30. Score: 4 | Accepted answer

It is helpful to understand the meaning of Fourier coefficients in the directional wave spectrum analysis.

A pitch-and-roll buoy measures time series of water elevation $\eta$, and slopes in both Cartesian directions, $\eta\_x$ and $\eta\_y$. The Fourier components $a\_0$, $a\_1$, $b\_1$, $a\_2$, $b\_2$, are related to cross-spectra of elevation and slope time series (introduced by Longuet-Higgins 1963):

$$
a\_0 = \dfrac{C\_{11}}{\pi}
$$

$$
a\_1 = \dfrac{Q\_{12}}{\pi}
$$

$$
b\_1 = \dfrac{Q\_{13}}{\pi}
$$

$$
a\_2 = \dfrac{C\_{22}-C\_{33}}{k^2\pi}
$$

$$
b\_2 = \dfrac{2C\_{23}}{k^2\pi}
$$

where the cross-spectra are:

$$
C\_{11}(f) = C[\eta(f)\eta(f)] = \int\_0^{2\pi} E(f,\theta)d\theta
$$

$$
Q\_{12}(f) = -Q[\eta(f)\eta\_x(f)]k^{-1} = -\int\_0^{2\pi} E(f,\theta)\cos{\theta}d\theta
$$

$$
Q\_{13}(f) = -Q[\eta(f)\eta\_y(f)]k^{-1} = -\int\_0^{2\pi} E(f,\theta)\sin{\theta}d\theta
$$

$$
C\_{22}(f) = C[\eta\_x(f)\eta\_x(f)]k^{-2} = \int\_0^{2\pi} E(f,\theta)\cos^2{\theta} d\theta
$$

$$
C\_{23}(f) = C[\eta\_x(f)\eta\_y(f)]k^{-2} = \int\_0^{2\pi} E(f,\theta)\sin{\theta}\cos{\theta} d\theta
$$

$$
C\_{33}(f) = C[\eta\_y(f)\eta\_y(f)]k^{-2} = \int\_0^{2\pi} E(f,\theta)\sin^2{\theta} d\theta
$$

From here you can see that the mean direction, as defined in the paper that you cite, is:

$$
\theta\_0 = \arctan\left({\dfrac{b\_1}{a\_1}}\right) = \arctan\left({\dfrac{Q\_{13}}{Q\_{12}}}\right) = \arctan\left({\dfrac{\int\_0^{2\pi}E(f,\theta)\sin{\theta}d\theta}{\int\_0^{2\pi}E(f,\theta)\cos{\theta}d\theta}}\right)
$$

Because $Q\_{12}$ and $Q\_{13}$ are proportional to the integrated energy in the zonal and meridional directions, respectively, the arctangent of their ratio yields mean wave direction.

The authors could use higher moments to calculate direction, but this would have different physical meaning. In this specific case, using $\theta = \arctan\left({b\_2/a\_2}\right)$ would yield peak (dominant) wave direction in the context of pitch-and-roll buoys.

You see that the above analysis is inherently limited by the quantities that are measured, specifically $\eta$, $\eta\_x$, and $\eta\_y$. If higher order quantities, say $\eta\_{xx}$, $\eta\_{xy}$, $\eta\_{yy}$ were available, the directional spectrum could be described at a higher accuracy even using just Fourier decomposition.

Further reading:

* [Oltman-Shay and Guza (1984)](http://journals.ametsoc.org/doi/abs/10.1175/1520-0485%281984%29014%3C1800%3AADAOWD%3E2.0.CO%3B2)
* [Nondirectional and directional wave data analysis procedures at NDBC](http://www.ndbc.noaa.gov/wavemeas.pdf)

> **Comments:**
> **username_0**: Thank you again - you answered another question of mine a few days ago, you're very helpful!  

---

