Earth:Ekman layer

From HandWiki
Short description: Force equilibrium layer in a liquid
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. Ekman layers occur both in the atmosphere and in the ocean.

There are two types of Ekman layers. The first type occurs at the surface of the ocean and is forced by surface winds, which act as a drag on the surface of the ocean. The second type occurs at the bottom of the atmosphere and ocean, where frictional forces are associated with flow over rough surfaces.

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 begins by assuming a neutrally stratified fluid, a balance between the forces of pressure gradient, Coriolis and turbulent drag.

[math]\displaystyle{ \begin{align} -fv &= -\frac{1}{\rho_o} \frac{\partial p}{\partial x}+K_m \frac{\partial^2 u}{\partial z^2}, \\[5pt] fu &= -\frac{1}{\rho_o} \frac{\partial p}{\partial y}+K_m \frac{\partial^2 v}{\partial z^2}, \\[5pt] 0 &= -\frac{1}{\rho_o} \frac{\partial p}{\partial z}, \end{align} }[/math]

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

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. Each of these situations can be accounted for through the boundary conditions applied to the resulting system of ordinary differential equations. The separate cases of top and bottom boundary layers are shown below.

Ekman layer at the ocean (or free) surface

We will consider boundary conditions of the Ekman layer in the upper ocean:[2]

[math]\displaystyle{ \text{at } z = 0 : \quad A \frac{\partial u}{\partial z} = \tau^x \quad \text{and} \quad A \frac{\partial v}{\partial z} = \tau^y, }[/math]

where [math]\displaystyle{ \ \tau^x }[/math] and [math]\displaystyle{ \ \tau^y }[/math] are the components of the surface stress, [math]\displaystyle{ \ \tau }[/math], of the wind field or ice layer at the top of the ocean, and [math]\displaystyle{ \ A \equiv \rho K_m }[/math] is the dynamic viscosity.

For the boundary condition on the other side, as [math]\displaystyle{ \ z \to -\infty : u \to u_g, v \to v_g }[/math], where [math]\displaystyle{ \ u_g }[/math] and [math]\displaystyle{ \ v_g }[/math] are the geostrophic flows in the [math]\displaystyle{ \ x }[/math] and [math]\displaystyle{ \ y }[/math] directions.

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:

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

The value [math]\displaystyle{ d }[/math] is called the Ekman layer depth, and gives an indication of the penetration depth of wind-induced turbulent mixing in the ocean. Note that it varies on two parameters: the turbulent diffusivity [math]\displaystyle{ K_m }[/math], and the latitude, as encapsulated by [math]\displaystyle{ f }[/math]. For a typical [math]\displaystyle{ K_m=0.1 }[/math] m[math]\displaystyle{ ^2 }[/math]/s, and at 45° latitude ([math]\displaystyle{ f=10^{-4} }[/math] s[math]\displaystyle{ ^{-1} }[/math]), then [math]\displaystyle{ d }[/math] is approximately 45 meters. This Ekman depth prediction does not always agree precisely with observations.

This variation of horizontal velocity with depth ([math]\displaystyle{ -z }[/math]) 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

[math]\displaystyle{ 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 ]. }[/math]

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.

Ekman layer at the bottom of the ocean and atmosphere

The traditional development of Ekman layers bounded below by a surface utilizes two boundary conditions:

  • A no-slip condition at the surface;
  • The Ekman velocities approaching the geostrophic velocities as [math]\displaystyle{ z }[/math] goes to infinity.

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 planetary 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 in the ocean were made in the Arctic Ocean from a drifting ice floe in 1958.[8] More recent observations include (not an exhaustive list):

  • 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]
  • In the eastern tropical Pacific, at 2°N, 140°W, using 5 current meters between 5 and 25 meters depth.[13] This study noted that the geostrophic shear associated with tropical stability waves modified the Ekman spiral relative to what is expected with horizontally uniform density.
  • North of the Kerguelen Plateau during the 2008 SOFINE experiment [14]

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][14]

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 978-0-13-353301-9. 
  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 978-0-521-84969-2. 
  3. ↑ Ekman, V.W. (1905). "On the influence of the earth's rotation on ocean currents". Ark. Mat. Astron. Fys. 2 (11): 1–52. 
  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 978-0-12-354015-7. 
  5. ↑ Santala, M. J.; Terray, E. A. (1992). "A technique for making unbiased estimates of current shear from a wave-follower". Deep-Sea Research 39 (3–4): 607–622. doi:10.1016/0198-0149(92)90091-7. Bibcode: 1992DSRA...39..607S. 
  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). 
  7. ↑ Weller, R.A.; Davis, R.E. (1980). "A vector-measuring current meter". Deep-Sea Research 27 (7): 565–582. doi:10.1016/0198-0149(80)90041-2. Bibcode: 1980DSRA...27..565W. 
  8. ↑ Hunkins, K. (1966). "Ekman drift currents in the Arctic Ocean". Deep-Sea Research 13 (4): 607–620. doi:10.1016/0011-7471(66)90592-4. Bibcode: 1966DSRA...13..607H. 
  9. ↑ Davis, R.E.; de Szoeke, R.; Niiler., P. (1981). "Part II: Modelling the mixed layer response". Deep-Sea Research 28 (12): 1453–1475. doi:10.1016/0198-0149(81)90092-3. Bibcode: 1981DSRA...28.1453D. 
  10. ↑ 10.0 10.1 Price, J.F.; Weller, R.A.; Schudlich, R.R. (1987). "Wind-Driven Ocean Currents and Ekman Transport". Science 238 (4833): 1534–1538. doi:10.1126/science.238.4833.1534. PMID 17784291. Bibcode: 1987Sci...238.1534P. 
  11. ↑ 11.0 11.1 Chereskin, T.K. (1995). "Direct evidence for an Ekman balance in the California Current". Journal of Geophysical Research 100 (C9): 18261–18269. doi:10.1029/95JC02182. Bibcode: 1995JGR...10018261C. 
  12. ↑ 12.0 12.1 Lenn, Y; Chereskin, T.K. (2009). "Observation of Ekman Currents in the Southern Ocean". Journal of Physical Oceanography 39 (3): 768–779. doi:10.1175/2008jpo3943.1. Bibcode: 2009JPO....39..768L. 
  13. ↑ Cronin, M.F.; Kessler, W.S. (2009). "Near-Surface Shear Flow in the Tropical Pacific Cold Tongue Front". Journal of Physical Oceanography 39 (5): 1200–1215. doi:10.1175/2008JPO4064.1. Bibcode: 2009JPO....39.1200C. 
  14. ↑ 14.0 14.1 Roach, C.J.; Phillips, H.E.; Bindoff, N.L.; Rintoul, S.R. (2015). "Detecting and Characterizing Ekman Currents in the Southern Ocean". Journal of Physical Oceanography 45 (5): 1205–1223. doi:10.1175/JPO-D-14-0115.1. Bibcode: 2015JPO....45.1205R. http://ecite.utas.edu.au/99058. 

External links

fr:Spirale d'Ekman#Couche d'Ekman



Retrieved from "https://handwiki.org/wiki/index.php?title=Earth:Ekman_layer&oldid=2945699"