エクマン境界層

エクマン境界層 は、コリオリ力と渦粘性がつりあう流体に見られる層である。ヴァン・ヴァルフリート・エクマン によって初めて記述された。

エクマン境界層の一般解

二次元回転平面流体において、圧力傾度力、コリオリ力、渦粘性が釣り合っている状況を考える。すなわち。平面流れ(u,v)が次の微分方程式系 に従うと考える。

${\begin{aligned}-fv&=-{\frac {1}{\rho _{o}}}{\frac {\partial p}{\partial x}}+K_{z}{\frac {\partial ^{2}u}{\partial z^{2}}},\\fu&=-{\frac {1}{\rho _{o}}}{\frac {\partial p}{\partial y}}+K_{z}{\frac {\partial ^{2}v}{\partial z^{2}}},\\\end{aligned}}$

ここで地衡流$(u_{g},v_{g})$を考えると、

${\begin{aligned}-fv_{g}&=-{\frac {1}{\rho _{o}}}{\frac {\partial p}{\partial x}}\\fu_{g}&=-{\frac {1}{\rho _{o}}}{\frac {\partial p}{\partial y}}\\\end{aligned}}$

非地衡流成分$u^{\prime }=u-u_{g},v^{\prime }=v-v_{g}$ が従う方程式は

${\begin{aligned}-fv^{\prime }&=K_{z}{\frac {\partial ^{2}u^{\prime }}{\partial z^{2}}},\\fu^{\prime }&=K_{z}{\frac {\partial ^{2}v^{\prime }}{\partial z^{2}}},\\\end{aligned}}$

以下では,非地衡流成分$u^{\prime },v^{\prime }$を省略して$u,v$と記す。

ここで複素速度$V=u+iv,(i={\sqrt {-1}})$ を用いると、次の二階連立偏微分方程式にまとめることができる。

$K_{z}{\frac {\partial ^{2}V}{\partial z^{2}}}=ifV$

この微分方程式の一般解は

$V=A\exp \left[\left({\frac {if}{K_{z}}}\right)^{\frac {1}{2}}z\right]+B\exp \left[-\left({\frac {if}{K_{z}}}\right)^{\frac {1}{2}}z\right]$

ここで$h_{E}\equiv {\sqrt {2K_{z}/f}}$ と置くと,${\sqrt {i}}=(1+i)/{\sqrt {2}}$ であるから

$V=A\exp \left[(1+i){\frac {z}{h_{E}}}\right]+B\exp \left[-(1+i){\frac {z}{h_{E}}}\right]$

この$h_{E}$は、エクマン境界層の厚さの目安のひとつになっている。

海上を吹く風によって形成されるエクマン境界層

海上を吹く風が海面に及ぼす風応力によって海面に運動がおこる。この運動が渦粘性によって次第に下層まで伝わり、エクマン境界層が生じる。 風が長期間一定であると仮定したときのエクマン境界層について記述する。 これは上記の一般解に2つの境界条件を課すと解くことができる。

境界条件のひとつは。海洋表層における力が風応力による摩擦によって生じるという仮定から導かれる。

$\rho K_{z}{\frac {\partial V}{\partial z}}{\Big |}_{z=0}=\tau ,\tau =\tau _{x}+i\tau _{y}$

もうひとつは、いま考えている層が海面に沿った境界層であることから導かれる。海はエクマン層に比べて 無限に深いとし、海底ではエクマン層に由来する流れが消えていると仮定する。

$V^{\prime }(-\infty )=0$

2つ目の境界条件から$B=0$ がわかる。次に1つ目の境界条件から、

${\begin{aligned}\rho K_{z}{\frac {\partial V}{\partial z}}{\Big |}_{z=0}=A\rho K_{z}{\frac {1+i}{h_{E}}}=\tau ,\\A={\frac {\tau h_{E}}{\rho K_{z}(1+i)}}={\frac {\tau h_{E}(1-i)}{2\rho K_{z}}},\\\end{aligned}}$

これを一般解に代入して整理し、実部と虚部に分けると、

${\begin{aligned}u=u_{g}+{\frac {h_{E}}{2\rho K_{z}}}\left[(\tau _{x}+\tau _{y})\cos {\frac {z}{h_{E}}}+(\tau _{x}-\tau _{y})\sin {\frac {z}{h_{E}}}\right],\\v=v_{g}+{\frac {h_{E}}{2\rho K_{z}}}\left[(-\tau _{x}+\tau _{y})\cos {\frac {z}{h_{E}}}+(\tau _{x}+\tau _{y})\sin {\frac {z}{h_{E}}}\right],\\\end{aligned}}$

特に$z=0$において,地衡流を除いた風応力のみによる吹送流をみると,

$u_{0}={\frac {h_{E}}{2\rho K_{z}}}(\tau _{x}+\tau _{y})\qquad ,v_{0}={\frac {h_{E}}{2\rho K_{z}}}(-\tau _{x}+\tau _{y})$

すなわち海表面における流れの速度ベクトルは,風応力ベクトルの右45度の角をなす.(北半球)

関連項目

• エクマン螺旋
• エクマン輸送

参考文献

気象力学通論　小倉義光