Ekman layer

From Infogalactic: the planetary knowledge core
Jump to: navigation, search
The Ekman layer is the layer in a fluid where the flow is the result of a balance between pressure gradient, Coriolis and turbulent drag forces. In the picture above, the wind blowing North creates a surface stress and a resulting Ekman spiral is found below it in the column of water.

The Ekman layer is the layer in a fluid where there is a force balance between pressure gradient force, Coriolis force and turbulent drag. It was first described by Vagn Walfrid Ekman.

History

Ekman developed the theory of the Ekman layer after Fridtjof Nansen observed that ice drifts at an angle of 20°-40° to the right of the prevailing wind direction while on an Arctic expedition aboard the Fram. Nansen asked his colleague, Vilhelm Bjerknes to set one of his students upon study of the problem. Bjerknes tapped Ekman, who presented his results in 1902 as his doctoral thesis.[1]

Mathematical formulation

The mathematical formulation of the Ekman layer can be found by assuming a neutrally stratified fluid, with horizontal momentum in balance between the forces of pressure gradient, Coriolis and turbulent drag.


\begin{align}
  -fv &= -\frac{1}{\rho_o} \frac{\part p}{\part x}+K_m \frac{\part^2 u}{\part z^2}, \\
   fu &= -\frac{1}{\rho_o} \frac{\part p}{\part y}+K_m \frac{\part^2 v}{\part z^2}, \\
    0 &= -\frac{1}{\rho_o} \frac{\part p}{\part z},
\end{align}

where \ u and \ v are the velocities in the \ x and \ y directions, respectively, \ f is the local Coriolis parameter, and \ K_m is the diffusive eddy viscosity, which can be derived using mixing length theory. Note that p is a modified pressure: we have incorporated the hydrostatic of the pressure, to take account of gravity.

Boundary conditions

There are many regions where an Ekman layer is theoretically plausible; they include the bottom of the atmosphere, near the surface of the earth and ocean, the bottom of the ocean, near the sea floor and at the top of the ocean, near the air-water interface. Different boundary conditions are appropriate for each of these different situations. We will consider boundary conditions of the Ekman layer in the upper ocean:[2]


  \text{at } z = 0 : 
  \quad A \frac{\part u}{\part z} = \tau^x \quad \text{and} 
  \quad A \frac{\part v}{\part z} = \tau^y,

where \ \tau^x and \ \tau^y are the components of the surface stress, \ \tau , of the wind field or ice layer at the top of the ocean and \ u_g and \ v_g are the geostrophic flows in the \ x and \ y directions – as \ z \to \infty : u \to u_g, v \to v_g.

In the other situations, other boundary conditions, such as the no-slip condition, may be appropriate instead.

Solution

Three views of the wind-driven Ekman layer at the surface of the ocean in the Northern Hemisphere. The geostrophic velocity is zero in this example.

These differential equations can be solved to find:


\begin{align} 
  u &= u_g + \frac{\sqrt{2}}{fd}e^{z/d}\left [\tau^x \cos(z/d - \pi/4) - \tau^y \sin(z/d - \pi/4)\right ],
  \\
  v &= v_g + \frac{\sqrt{2}}{fd}e^{z/d}\left [\tau^x \sin(z/d - \pi/4) + \tau^y \cos(z/d - \pi/4)\right ],
  \\
d &= \sqrt{2 K_m/f}.

\end{align}

This variation of horizontal velocity with depth (-z) is referred to as the Ekman spiral, diagrammed above and at right.

By applying the continuity equation we can have the vertical velocity as following

w = \frac{1}{f\rho_o}\left [-\left (\frac{\partial \tau^x}{\partial x} + \frac{\partial \tau^y}{\partial y} \right )e^{z/d}\sin(z/d) + \left (\frac{\partial \tau^y}{\partial x} - \frac{\partial \tau^x}{\partial y} \right )(1-e^{z/d}\cos(z/d))\right ].

Note that when vertically integrated the volume transport associated with the Ekman spiral is to the right of the wind direction in the Northern Hemisphere.

Experimental observations of the Ekman layer

There is much difficulty associated with observing the Ekman layer for two main reasons: the theory is too simplistic as it assumes a constant eddy viscosity, which Ekman himself anticipated,[3] saying

and because it is difficult to design instruments with great enough sensitivity to observe the velocity profile in the ocean.

Laboratory demonstrations

The bottom Ekman layer can readily be observed in a rotating cylindrical tank of water by dropping in dye and changing the rotation rate slightly.[1] Surface Ekman layers can also be observed in rotating tanks.[2]

In the atmosphere

In the atmosphere, the Ekman solution generally overstates the magnitude of the horizontal wind field because it does not account for the velocity shear in the surface layer. Splitting the boundary layer into the surface layer and the Ekman layer generally yields more accurate results.[4]

In the ocean

The Ekman layer, with its distinguishing feature the Ekman spiral, is rarely observed in the ocean. The Ekman layer near the surface of the ocean extends only about 10 – 20 meters deep,[4] and instrumentation sensitive enough to observe a velocity profile in such a shallow depth has only been available since around 1980.[2] Also, wind waves modify the flow near the surface, and make observations close to the surface rather difficult.[5]

Instrumentation

Observations of the Ekman layer have only been possible since the development of robust surface moorings and sensitive current meters. Ekman himself developed a current meter to observe the spiral that bears his name, but was not successful.[6] The Vector Measuring Current Meter [7] and the Acoustic Doppler Current Profiler are both used to measure current.

Observations

The first documented observations of an Ekman-like spiral were made in the Arctic Ocean from a drifting ice flow in 1958.[8] More recent observations include:

  • The 1980 Mixed Layer Experiment[9]
  • Within the Sargasso Sea during the 1982 Long Term Upper Ocean Study [10]
  • Within the California Current during the 1993 Eastern Boundary Current experiment [11]
  • Within the Drake Passage region of the Southern Ocean [12]
  • North of the Kerguelan Plateau during the 2008 SOFINE experiment [13]

Common to several of these observations spirals were found to be 'compressed', displaying larger estimates of eddy viscosity when considering the rate of rotation with depth than the eddy viscosity derived from considering the rate of decay of speed.[10][11][12][13]

See also

References

  1. Cushman-Roisin, Benoit (1994). "Chapter 5 - The Ekman Layer". Introduction to Geophysical Fluid Dynamics (1st ed.). Prentice Hall. pp. 76–77. ISBN 0-13-353301-8.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  2. 2.0 2.1 Vallis, Geoffrey K. (2006). "Chapter 2 – Effects of Rotation and Stratification". Atmospheric and Oceanic Fluid Dynamics (1st ed.). Cambridge, UK: Cambridge University Press. pp. 112–113. ISBN 0-521-84969-1.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  3. Ekman, V.W. (1905). "On the influence of the earth's rotation on ocean currents". Ark. Mat. Astron. Fys. 2 (11): 1–52.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  4. 4.0 4.1 Holton, James R. (2004). "Chapter 5 – The Planetary Boundary Layer". Dynamic Meteorology. International Geophysics Series. 88 (4th ed.). Burlington, MA: Elsevier Academic Press. pp. 129–130. ISBN 0-12-354015-1.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  5. Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  6. Rudnick, Daniel (2003). "Observations of Momentum Transfer in the Upper Ocean: Did Ekman Get It Right?". Near-Boundary Processes and their Parameterization. Manoa, Hawaii: School of Ocean and Earth Science and Technology.<templatestyles src="Module:Citation/CS1/styles.css"></templatestyles>
  7. Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  8. Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  9. Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  10. 10.0 10.1 Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  11. 11.0 11.1 Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  12. 12.0 12.1 Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).
  13. 13.0 13.1 Lua error in Module:Citation/CS1/Identifiers at line 47: attempt to index field 'wikibase' (a nil value).

External links