ここで、ラグランジュ微分を表すDと平均水深を表すDを混同しないでください。次節の渦位保存則の証明には、こちらのより厳密な方程式を用います。簡略化された(1’)式の両辺をyで偏微分、(2’)式の両辺をxで偏微分するととなる。(2”)式から(1”)式を辺々引き算するとが得られる。おけるケルビン波」1.(2)の(3)式と比較してみられたし。そこの[補足説明2]で注意したように(∂u/∂x+∂v/∂y)はベクトルの“発散”を表しています。そこで説明したように今考えている場所の柱状領域に周囲から流体が流入すれば渦度は増えます。その柱状領域を横切るy方向の流れvが存在してもその場所の渦度は変化します。それは座標系の回転角速度が緯度依存性を持っているので流体が高緯度へ移動すると、渦度がより小さい低緯度の流体要素でもってその場所の柱状領域が占めなおされるからです。新たに付け加わったβvの項がそのことを表しています。この場合も、オイラー形式の微分方程式ですから、座標空間に固定された微少柱状領域での話です。(5)、(5’)式は、運動方程式(1’)、(2’)と連続の式(3)式以上の内容を含んでいるわけではないのですが、この様な表現にすると渦度が変化するメカニズムがよくわかります。座標空間に固定された座標値を変数とするオイラー形式に乗っ取った厳密な形で書いてみる。そのために、左辺の時間微分∂v/∂tを厳密なラグランジュ微分[Dv/Dt=∂v/∂t+(v・grad)v]にした運動方程式を用いる。厳密な(A)式の両辺をyで偏微分、(B)式の両辺をxで偏微分するとが得られる。(B’)式から(A’)式を辺々引き算するとこれら三つの方程式は、その形は微妙に異なりますが、いずれも(x,y)座標が固定された場所の微少柱状領域に存在する流体の渦度の変化を表す微分方程式です。厳密な“渦度方程式”です。3.渦位保存則ここで、厳密な“連続の式”がであることに注意して、上式を書き直してみる。これを“渦位保存則”と言う。vorticity)”と“渦度(vorticity)”ζ=(∂v/∂x-∂u/∂y)を混同しないでください。両者は全く別の量です。力学」の3.(3)2.ラグランジュの渦定理や5.(4)2.ヘルムホルツの渦定理の証明と比較してみられたし。そこでは[渦度ζ]÷[密度ρ]や[渦度ω]×[渦糸の断面積dσ]が保存されていたが、ここでは([渦度ζ]+[コリオリ力項f])÷[その点での流体の厚さD+h]が保存される。このとき、(∂v/∂x+∂u/∂y)=0の非発散性(つまり波高が変化しない)流れの場合、上記の渦位保存則はに簡単化される。vorticity)”は、もともと有限の長さの渦糸に対して、渦糸の渦度を長さで割ったものです。そのとき渦糸の長さが変化するとき渦度は変化するのですが、渦度を渦糸の長さで割ったものは変化しなかった。それは角運動量保存則が働くからです。アイススケートのスピンで手を広げたり狭めたりしたとき回転の角速度は変わるが角運動量は保存されたのと同じです。渦糸の長さが長くなれば、太さは細くなります。このあたりを説明する定理が流体力学におけるラグランジュの渦定理・ヘルムホルツの渦定理・ケルビンの循環定理です。そのとき用いた運動方程式には回転を引き起こすような項は存在しなかったことを思い出してください。ここでは回転座標系上での話なのでf0+βyの項も加わっています。そのため、1.(2)1.で説明したように、渦糸の長さが変化(太さが変化)する場合、f0+βyの項が効いてきて渦度の大きさが変化すると同時に普通の意味のポテンシャル渦度も変化します。単純な角運動量保存則は成り立たないのです。そのため、ポテンシャル渦度をf0+βyの項も考慮した上記の形にしないと渦位保存則は成り立たない。大規模な大気や海洋の流れは、ほぼ水平流なので、鉛直軸をもった渦糸の集合とみなせる。このため、その時間変化は上記のポテンシャル渦度“渦位”の保存則によってよく記述できる。大気や海洋の厚さはせいぜい数km〜数十kmであるのに対して、今から考察する波動の波長や広がりは数百km〜数千kmであることは常に意識しておかねばなりません。[補足説明2]本によっては、(f0+βy)のことを[宇宙空間からみた地球の回転の事だから]惑星渦度、(∂v/∂x-∂u/∂y)のことを[地面に相対的な渦度だから]相対渦度と言い、(∂v/∂x-∂u/∂y)+(f0+βy)のことを絶対渦度と言っているものがあります。絶対静止の宇宙空間からみて惑星渦度と相対渦度の和が保存されることを絶対渦度保存則と呼んでいる。確かに絶対静止の空間から見ればその通りなのだから誤りではないのですが、ここでの議論の様に地球表面上に設置したデカルト平面座標(β平面)上で論じている場合には、そのような言い方は誤解を招きます。デカルト平面座標の上で論じているのは、別稿でも説明したように、仮想的な力である慣性力(コリオリ力)を導入すれば、議論がとても簡単になるからです。して初めて(∂v/∂x+∂u/∂y)の効果(波高の変化による渦度の変化)なども旨く取り込むことができます。[補足説明1]で説明した通りです。偏西風の蛇行は地球規模の現象だから、絶対静止の宇宙空間から見る立場でその挙動を説明したくなるのは山々ですが、今は地球表面上に設定したデカルト平面座標(β平面)上で論じているのですから、その誘惑は押さえておかねばなりません。渦位保存則は、流体の実質部分に付属する“渦位”が時間的に変化しない事を示している。流体の運動が運動方程式と連続の式を満足している限りに於いて、それが波の伝播のようなものであろうと、定常的な流れであろうと、任意の運動状態にある流体の実質部分[具体的には場所(x,y)の微少柱状領域を占める流体]が持つ渦位は、場所的(二次元的に考えている)にも時間的にも保存されて一定です。ある地点(x,y)に存在する微少柱状領域の波高が、波の伝播に伴って高くなったとする。そのとき渦位が変化しないと言うことは、その微少柱状領域の流体の深さD+hが増大するのだから渦度ζは増大しなければならない事を言っている。これはコリオリ力のためにそうなる。別稿の[補足説明2]で説明したようにケルビン波において、そのことがまさに実現されている。ある地点(x,y)の微少柱状領域を占めていた流体がΔyだけyの正方向に移動したとする。その場合、(f0+βy)が増大するのだから、渦位が保存されるためには渦度ζは減少しなければならない。これはコリオリ力の緯度依存性(β効果)のためにそうなるのです。このことが、後ほど説明するロスビー波の蛇行が保持され、波が西進することの理由です。これらは1.(2)1.ですでに説明した事柄です。それを保存則の形で述べているにすぎない。8.文献(3)まとめ議論の出発点となる連立方程式をまとめておきます。厳密な式(A)(B)(C)は非線形方程式だから解析的に解けないので、今後は最初に述べた微少変位の場合に成り立つ簡略化された方程式を用いる。議論て必要になる概念を、いくつかの特別な場合で説明します。(1)時間的変化が無視できる流れ(地衡流)最初に議論する特別な例は運動が時間的に変化しない定常的な場合です。方程式の時間の偏微分項が無視できる場合です。時間的な変化が無視できる場合は、前述の方程式の∂/∂tの項がすべてゼロとおけるのでとなる。ここで利用した圧力傾度と水面傾度の関係式は別稿で説明した。(6)、(7)式は、“定常的”流れでは[水面傾度(圧力傾度)]と[コリオリ力]が常に釣り合ている事を示している。この場合に実現されるのは、波高が時間的に変化しない“非発散性”の流れであることを示している。波高hは時間的には変化しないが、場所的には変化している。currents)と言う。これは“地球回転と平衡した流れ”という意味で名付けられたようです。2.水面傾度(圧力傾度)と流速ベクトル“地衡流”では[水面傾度(圧力傾度)を表すベクトル]と[流速ベクトル)]は必ず直交します。実際、(6)式にuを乗じたものと、(7)式にvを乗じたものはとなる。これらの式を辺々足し算するとが得られる。これらはベクトルの内積の定義から明らかなようにが直交することを示している。そのとき、(6)、(7)式の符号は北半球の場合です。uとvのどちらか一方のみがゼロで無い特別な場合を考えれば解るように、北半球では水面の高い所(圧力の高い所)を右に見ながら流れる。定常的な流れに対しては、“流線関数(streamfunction)Ψ”を導入することができて、流れの様子を視覚的・感覚的に解りやすく表現できます。力学」2.(2)で説明しましたが、重要な概念なのでもう一度説明します。ある瞬間に流体の微小部分がもつ速度の方向を連ねた曲線を流線と言う。各瞬間に於ける流体微少部分が持つ速度ベクトルはその点に於ける流線の接線方向です。流線と速度ベクトルは平行で同じ方向を向いている。ここで導入される流線こそ流体力学を特徴づける最重要な概念で、それは連続の方程式と密接に関係している。速度ベクトルは流線に沿った線分の方向を向く。二次元流における流線の線分ベクトルを(dx、dy)とすると流線の定義よりが成り立つ。ところで定常流(水面高の時間的変化が無い)の場合には連続の式(3)が(8)式の様に簡単化される。(8’)式は、流線の定義式における(−v)dx+udyなる微分形式(9)式がある一つの関数(−Ψ)の完全微分になる必要十分条件であることを表している。[この言葉の意味は別稿「絶対温度とは何か(積分因子とは何か)」5.(2)を参照されたし]そのため、この−Ψ(x,y)はxy平面上の一つの連続な関数曲面z=−Ψ(x,y)を形づくりと表現できる。力学」2.(2)の場合と逆にしているのは、流線関数を(北半球での)水面の波高と対応させるためです。それ以上の意味は有りません。そのとき、速度成分u、-vはある関数−Ψ(x,y)の勾配で表すことができる。実際(10)式が成り立てばとなって(8’)式が常に満たされていることが解る。−Ψ(x,y)を用いると流線が満たすべき(9)式はとなる。ところで、−dΨ=0を満たす線分ベクトル(dx,dy)の方向は関数曲面上のΨ=一定の曲線をxy平面に投影した等高線の方向です。そのため、(11)式は[等高線の線分ベクトル(dx,dy)]と[関数曲面Ψ(x,y)の最大勾配の方向ベクトル(∂Ψ/∂x,∂Ψ/∂y)]の内積がゼロで互いに直交していることを表している。(10)式はとなるので[速度ベクトル(u,v)]と[勾配ベクトル(∂Ψ/∂x,∂Ψ/∂y)]の内積がゼロで直交している事を意味している。そのため曲面z=Ψ(x,y)の等高線と速度ベクトル(u,v)は平行になり、曲面z=Ψ(x,y)の等高線が流線になる。そのときΨがyの正方向に高くなっていればx方向の速度成分uが負、Ψがxの正方向に高くなっていればy方向の速度成分vが正となるので、流体は常にΨの値の高い所を右に見ながら、等高線に沿って流れる。等高線の密集している所[つまり曲面Ψ(x,y)の勾配が大きいところ]ほど高速で流れることを示している。この関数Ψ(x,y,t)を“流れ関数(streamfunction)”あるいは“流線関数”という。ここで、さらに以下のことに注目してほしい。コリオリ力の場所的な変化を表すβ項が小さくてf=f0+βy≒f0と近似できるときには、運動方程式(6)、(7)式と(10)式を比較してみれば明らかなようにとおける。波高や圧力が、ほぼそのまま流線関数に比例する量だと見なせる。そのとき流体は、波高(圧力)の高い所を(北半球では)右側に見ながら、h(x,y)の等高線あるいはp(x,y)の等圧線にほぼ沿って流れる。4.地衡流(地衡風)流体の運動を規制している最大の制約が“連続の式”で、これこそが流体力学を一つの学問領域として成り立たせている中心原理です。連続の式から流線は、xy平面上のz軸方向に分布する連続でなめらかな関数曲面流線関数z=Ψ(x,y)の等高線を(x,y)平面に射影した曲線になる。そのとき、さらに、“運動方程式”が[コリオリ力]と[圧力傾度]が釣り合うような形をしていて、時間的な変化が無視できる定常的な流れでβ項が微少なら、流線関数Ψを波高分布関数z=h(x,y)や圧力分布関数z=p(x,y)で近似的に置き換えることができる。流体はh、pの等高線に沿って流れ、h、pの勾配が急な所ほど速く流れる。h、pの勾配が急なところは等高線[流線]が密集している所ですから、流れの幅が狭まり流速は増大する。流れが速いほど右向き(北半球)のコリオリ力は強まる。北半球で流体が動くと、コリオリ力のために常に右向きの力を受けるが、流れが定常的になるためには、その力を打ち消す反対向きの力が常に働いていなければならない。それが進行方向の右に高く積もった水面によって生じる圧力傾度力です。圧力傾度力は水面の勾配が大きいほど強くなる。最初に掲げた“定常流”で成り立つ(6)、(7)式では[水面傾度(圧力傾度)]と[コリオリ力]が両辺で互いに釣り合ていた。それがまさにこのことを示している。currents)と言う。地衡流は別稿「台風のメカニズム」2.などで説明したが、大規模な大気や海洋の定常的な流動運動では良い近似解となる。β項が小さくてf=f0+βy≒f0と近似できるときには、(10)式から明らかなように等圧線に沿って吹く風の速度Vは、nを等圧線に対して直角方向の距離とすると次のように表される。海洋では「吹送流の西岸強化」で説明した大洋を循環する大規模海流がまさに地衡流の例です。そこでは表面を吹く風の駆動力と海底の摩擦力が釣り合って定常流が実現されていたが、本来そのような駆動力が無くても、平べったい容器に最初からそういった流れが存在して、内部摩擦が無い理想流体であればその流れは存在し続けることができます。この場合にも南北方向の圧力傾度を生じるメカニズムや偏西風を加速するメカニズムが別に有ります。大気が内部摩擦の無い理想流体と見なせて最初からそういった流れが存在するのならば、本来そのようなメカニズムが無くても地衡流は存在し解となります。後で説明するロスビー波はまさにこの地衡風(偏西風)に生じるものです。温帯低気圧や温帯高気圧の等圧線に沿って吹く風も基本的には地衡風です。あるとき、A点での風速|V|を求めてみる。おけるケルビン波」3.(3)で述べた解は、本来x方向に関して非定常な波動ですが、そこの[補足説明1]で述べたように各瞬間には地衡流と見なせます。8.文献(2)水面傾度(圧力傾度)が無視できる流れ(慣性振動)二番目の簡単化の例として運動方程式(6)、(7)式の右辺の水面高(圧力)の空間的な変化が無視できる場合を考察する。流体の運動を規定する力がコリオリ力のみの場合です。圧力項(水面高の場所的な変化)が無視できる場合には、運動方程式は次のように簡単化される。このとき空間微分項が無視できるなら、もとの浅水方程式系で連続の式の波高の時間的変化も無視できる。そのため、連続の式を考慮する必要はなくなる。近似できる場合にはとなる。(12)式の両辺をtで偏微分したものに(13)式を代入、あるいは(13)式の両辺ををtで偏微分したものに(12)式を代入するととなる。角振動数がf0の振動方程式が得られる。このとき、波高hの変化は無視しているので、“連続の式”や“渦位保存則”は意味を持たない。前項の振動方程式は、いわゆる単振動の式だから簡単に解ける。最初の式からが得られるが、これを(12)式に代入して積分するとが得られる。uとvの位相が90°異なっていることから解るように、これは流体が回転運動することを表している。実際、tを媒介変数としたuv平面上のグラフを描いてみると、速度ベクトル(u,v)の先端は時計回り(北半球)の円を描く。速度解を時間tで積分して距離に直す。積分定数を適当に調整するとが得られる。これらの式は、tを媒介変数としたxy平面上のグラフを描いてみれば明らかなように、時計回り(北半球)の円運動を表している。代入してみるととなるので、円運動の半径=(V0/f0)であることが解る。回転と言っても空間的に分布する流体のことですから、流体全体が同時に平行移動するように回転する運動です。コリオリ力が働く回転座標系の上で、物体がある初速度を持って動き始めると回転運動をして元の位置に戻ることを別稿「コリオリ力」3.で説明した。それと同じ事が流体全体に対して起こるのです。これらの運動については以下の事が言える。初速度V0が大きければ大きな半径の、小さければ小さな半径の回転運動を行う。上で求めた円の半径(V0/f0)がそのことを端的に表している。一回転に要する時間Tはすべて同じT=2π/f0となる。円の半径が初速度V0に比例するのだからそれは当然です。極では半日程度であり、それ以外の緯度ではそれをsinθ0で割った値となる。そのため赤道に近づくにつれて振動周期は長くなる。oscillation)と呼ばれる。説明した“極偏東風”や“貿易風”の様に、南方へ速度V0で全体が動き出した流体が、コリオリ力のために全体的にそろって風向きを変える様な現象です。最初にそのような速度を得るメカニズムは、太陽輻射による南北の不均等加熱なのですが、そういった初速が生じる原因についてはここでは立ち入っていません。伴う回転運動とは、そのメカニズムが異なります。低気圧や高気圧の風は前節で説明した、[圧力傾度]と[コリオリ力]が釣り合って定常流が実現されている地衡風ですが、ここではコリオリ力と初速度でのみ生じるものです。慣性振動の初期の観測例としてスエーデンのグスタフソン(T,Gustafson)とクレンベリイ(B.Kullinberg)が1933年8月17〜24日にバルト海で行ったものがある。強風による一定方向への海水移動に加わって周期約14hの円運動が見られた。この周期は観測点の緯度から計算れさる周期14h8mに良く一致している。中心ベクトル図と進行ベクトル図を下に示す。下図は緯度N13°付近の深さ30mの海水の運動のスペクトル分析です。緯度N13°での慣性振動の周期は53時間ですが、半日周期の潮汐とともに、慣性振動周期の運動に多くのエネルギーが分配されているのが読み取れる。このように慣性振動は様々な動きのノイズのなかに埋もれて存在する場合が多い。慣性振動は海面に強い風が吹いた後の海中でよく観測される。下図は深さ1000mで浮力が釣り合うようにした中立浮きの動きを追跡したものです。図から中立浮きが円を描きつつ、その場の一様流に流されている様子が読み取れる。そのとき、緯度N31°での慣性振動周期は23.3時間ですが、回転周期からまさしく慣性振動による運動であることが解る。8.文献(3)コリオリ力が無視できる流れ(重力波)三番目の簡単化の例は運動方程式(6)、(7)式の左辺のコリオリ力項が無視できる場合です。流体に働く力が水面傾度(圧力傾度)のみによって規定される場合です。1.運動方程式・連続の式・渦位保存則コリオリ力項が無視できる場合には、運動方程式は次のように簡単化される。このとき連続の式は前節の様に省略することはできません。波高の時間偏微分項を省略することはできません。ここで、速度成分(u,v)の波動方程式を導くだけならば、簡単化した渦度方程式(5)に於いてf=f0+βy=0と於いた式を時間積分して得られるを用いることができる。渦位保存則としては、h→0の場合にのみ成り立つ近似的なもので正確ではありません。正しい渦位保存則を得るには1.(2)1.で説明した渦度方程式(5)式でコリオリ力=0としたのでは簡単化しすぎています。1.(2)2.で説明した厳密な渦度方程式に於いてコリオリ力=0としなければなりません。その場合にはとなり、コリオリ力が存在しないときの正確な“渦位保存則”が得られます。“連続の式”(3)の両辺をtで偏微分した式に、“運動方程式”(14)の両辺ををxで偏微分したものと、(15)の両辺をyで偏微分したものを代入するとが得られる。さにら、“連続の式”(3)の両辺をxで偏微分した式に、“運動方程式”(14)の両辺をtで偏微分したものを代入し、さらに“渦度方程式”(16)をyで偏微分したものを代入するとが得られる。同様に、“連続の式”(3)の両辺をyで偏微分した式に、“運動方程式”(15)の両辺をtで偏微分したものを代入し、さらに“渦度方程式”(16)をxで偏微分したものを代入するするとが得られる。これらは単純な長波の理論でお馴染みのものです。コリオリ力項を含んだもっと一般的な場合は3.(1)で説明します。これらの別々の波動方程式から導かれる波高h解と流速(u,v)解は、運動方程式(14)と(15)式によって定まる関係も同時に満足しなければならない。3.波動方程式の解これらの波動方程式はお馴染みのもので、最も簡単な解を例示するとの様になる。速度成分(u,v)は(20)式を(14)、(15)式に代入して時間tで積分すれば求まる。これらは速度成分に関する波動方程式(18)、(19)式を満足しています。(20)式を(17)式に代入してみれば明らかなようにが満足されれねばならない。wave)と呼ばれています。front)”と言う。frequency)”あるいは“角周波数”と言う。普通の本では平均水深をH、波高をη、角振動数をωと表しているが、この稿では平均水深をDとしており、ωは地球の自転角速度、ω’はコリオリ力項の表現に用いている。そのため、この稿では波高をh、波高振幅をH、波数を(k,l)、角振動数をσで表している。k、l、σでは、解りにくいので別稿「波動方程式と一般解」3.(3)で説明した波長λ、周期Tの形に書き換える。このときそれぞれの量は下記の意味を持つ。これらの関係式は後に繰り返し利用します。ここで、H=1、λx=6、λy=8、T=12とした波高解の様子をアニメーション表示する。アニメーションから明らかなようにベクトル(k,l)は“波面”に垂直なベクトルで波が伝播していく方向を示す。このベクトルを“波数ベクトル”と言う。これらはすでに別稿「偏光とは何か」2.などで説明した。波数ベクトル(k,l)と波長λの関係は下図の様になる。図中にできる直角三角形はすべて互いに相似になる。それらの三角形について辺の比例関係を用いるととなる。“重力波”の角振動数σは波数ベクトルの大きさに比例する。そのため、伝播速度(位相速度)は、水深と重力加速度のみに依存し波長に依存しない。このように波長に依存しないことを分散性が無いと言う。“重力波”では流速ベクトル(u,v)の方向は波数ベクトル(k,l)の方向と一致する。流体粒子は波の伝播する方向と平行に運動する。後に説明する“慣性重力波”や“ロスビー波”ではそうなりません。流速は波高に比例し、流速ベクトルの方向は、h>0の部分で波の伝播方向を向き、h<0の部分で反対方向を向く。流速ベクトルの大きさはとなる。波高との関係は一次元の長波理論でお馴染みのものです。重力波の代表例がコリオリ力が無視できる場合の“津波”や“潮汐波”です。これらは地球規模の現象なので、実際はコリオリ力項を無視できません。ここでは津波や潮汐波が生じる原因については全く考慮していません。そのような波が存在したならば、ここで説明した方程式の解であると言うだけです。対流圏大気の成層中には周期10分程度で上下に振動する重力波が存在する。フィンランドのヴァイサラは20世紀の初め頃に気温観測などの目的で気球観測をするとき、気球の上昇する速さがいつも10分程度の周期で増減する事から大気中の重力波を発見した。重力波”、「天気」1973年11月号http://www.metsoc.jp/tenki/pdf/1973/1973_11_0571.pdf風が山を越えるときに、山岳の風下側にできる波動も重力波です。山岳によって強制的に持ち上げられた重い大気は周囲の軽い空気の中に存在する事になるので、、その密度差により引き続く状況で下方へ沈み込みます。この浮き沈みが波動となって下流に伝播していく。一様流中の内部波なので方程式系はここで説明したものより少し込み入ってきますが重力波であることには間違い有りません。風下波)と呼ばれ、グライダーの飛行や山岳地帯を飛ぶ飛行機に重要な影響をおよぼす。山岳波の位相は山に固定されています。密度の内部境界面に生じる波なのでエネルギーは上空に伝播していきます。1960年代におけるロケット観測(バラシュートでラジオゾンデを落として気温や風の高度分布を測定)から下図の様な情報が得られた。これらのデータを波動力学的に解析することから成層圏・中間圏は水平波長数十〜数百kmの重力波に満たされていることが解ってきた。その後のさらなる観測から、それらは山岳波、前線や台風、熱帯の積雲対流などの影響で対流圏中に生じた重力波が成層圏・中間圏に伝播したもので、中層大気の大きな変動(プラネタリー波)に本質的な影響を与えていることが解ってきた。ブラネタリー波とは地球規模で生じる慣性重力波やロスビー波のことです。波動の規模が大きくなってくると、必然的にコリオリ力の影響が出てきます。ここではまずコリオリ力は無視できないが、緯度変化に伴う影響(β効果)は無視できる限られた緯度・領域で生じる波動を考察します。ここで、すべての項を考慮した浅水方程式系[1.(3)参照]を考察する。β項を導入すると話は難しくなるので、ここではまずβ項が無視できてf=f0+βy≒f0=一定と置ける場合を議論します。コリオリ力と圧力傾度力の両方を考慮するが、コリオリ力の緯度依存性は無視できる場合です。f0+βy≒f0=一定と置くととなる。2.渦位保存則このとき(23)、(24)式から得られる渦度方程式は(5’)式に於いて、f=f0+βy≒f0=一定と置いたものになり、両辺を時間積分できます。なり、これは波動方程式を導くときに使えます。これは渦位保存則としては近似的なものです。正確な渦位保存則を導くには1.(2)2.で求めた正確な渦度方程式を時間積分する必要があります。それは1.(2)3.の保存則でβ=0と置いたもので、となります。これが正確な“渦位保存則”です。このときhが小さいときには(5”)になることは簡単な計算で確かめられます。運動方程式(23)と(24)の両辺をtで偏微分するととなる。ここで、波高hと流速(u,v)の間に成り立つ関係式を求めておきます。(23’)式に“運動方程式”(24)を代入するととなる。(24’)式に“運動方程式”(23)を代入するととなる。これらは波高解hから流速解(u,v)を導くときに用います。8.文献(2)波動方程式(23)(24)(3)は、3つの未知関数h、u、vが互いに入り交じった連立方程式なので、一つの未知関数のみを含む単一の式に変換する。おけるケルビン波」2.で求めており、そこの2ω’を現在の記号f0に置き換えればよい。繰り返しになりますが大事な所なのでもう一度説明します。運動方程式(23)の両辺をxで偏微分するととなる。運動方程式(24)の両辺をyで偏微分するととなる。ここで、“連続の式”(3)の両辺をtで偏微分したものに、“運動方程式”の両辺をxやyで偏微分した(23”)と(24”)を代入し、さらに“渦運動を表す式”(5”)の積分定数をゼロと置いたものを用いるとが得られる。“渦運動を表す式”(5”)の両辺をyで偏微分するととなる。(5”)の両辺をxで偏微分するととなる。ここで、“連続の式”(3)の両辺をxで偏微分したものに、“運動方程式”(23)の両辺をtで偏微分した(23’)と“運動方程式”(24)を代入し、さらに“渦運動を表す式”(5”)をyで偏微分して得られた(5*)用いるとが得られる。同様にして、“連続の式”(3)の両辺をyで偏微分したものに、f=f0=一定と置いた“運動方程式”(2’)の両辺をtで偏微分した(2”’)とf=f0=一定と置いた“運動方程式”(1’)を代入し、さらに“渦運動を表す式”(5”)をxで偏微分して得られた(5†)用いるとが得られる。結局、三つの量h、u、vは全て同じ形の“2階の同次波動方程式”(28)、(29)、(30)を満足する。これらは2.(3)で説明した“重力波の波動方程式”に、コリオリ力の影響を表すf0項が付け加わった形をしている。これが“慣性重力波の波動方程式”です。[補足説明]本来ならば、この式中に流体の運動を引き起こす力の項が導入されていてしかるべきです。いまは、それらの効果はすべて考えられていません。”を表す波動方程式の議論です。おけるケルビン波」3.(1)で説明したようにと表せます。ここでは見通しを良くするために、最も簡単な実数解を取り上げる。この波の様子は3.(3)で取り上げた重力波の解(20)式と同じですからそこのアニメーションを参照されてください。そのためk、l、σと波長λ、振動周期T、伝播速度μの関係も全く同じです。この波について以下の事が言える。(31)式を波動方程式(28)に代入するとが得られる。これは、(31)式が波動方程式(28)の解であるとき、σが満足しなければならない関係を表している。(32)式においてコリオリ力が無視できる[f0→0]場合には、重力波の関係式に一致する。重力波では波長が長くなる[(k2+l2)0.5→0]につれて、σ→0となり振動周期T=2π/σはいくらでも長くなれた。慣性重力波のσにはそれ以上小さくなれない最小値f0が存在する。振動周期Tは慣性振動周期2π/f0よりも長くなれない。(31)式を波長λと周期Tで表すととなる。(32)式が成り立つとき、波の伝播速度(位相速度)はとなる。これは、伝播速度(位相速度)が波長λ=2π/(k2+l2)0.5に依存して変化する事を示している。このように波の位相速度が波長(波数)に依存するとき、波は分散性を持つ(dispersive)と言う。deformation)”と言いますが、λ>Rdになると、位相速度の値が重力波のそれと顕著に異なってくる。波長が同じ場合、慣性重力波の位相速度は重力波の位相速度よりも常に速くなる。伝播速度が重力波の伝播速度(gD)0.5より速くなるのはコリオリ力が働くためです。コリオリ力の影響が無視できる(f0=0)場合には普通の重力波の伝播速度(gD)0.5になる。このとき、位相速度は波数ベクトルの絶対値の大きさ(波長)には依存するが、波数ベクトルの方向には依存しない。後で説明する“ロスビー波”は方向にも依存する。wave)と言います。wave)と呼ぶ場合があります。おけるケルビン波」の1.(2)[補足説明2]や1.(4)で説明したように(∂u/∂x+∂v/∂y)はベクトルの“発散”を表しており、考えている場所の柱状領域に周囲から流体が流入すれば波高の増大が生じるとともに渦度も変化します。それが波の復元力に慣性力(コリオリ力)が関わっている理由です。そのことを通じて位相速度に関係してくる。波高変化に伴って渦度が変化することは、次節で表示する流速ベクトルのアニメーションを見れば解ります。波高の変化に伴って隣り合う流体間の速度成分がずれあうように変化するさまが読み取れます。2.流速解h(x,y,t)に対応するu(x,y,t)やv(x,y,t)の解をを求めよう。そのためには、(31)式を(26)、(27)式に適用すればよい。するとやとなるのでが得られる。満足する(x,y)座標上に存在する流体粒子は、すべて同じ位相でそろって同じ方向、同じ速度で動く。波数ベクトル(k,l)と位置ベクトル(x,y)の内積がゼロだから互いに直交していることを意味する。一つの波面(同一波高の点を結んだライン)上の点がすべてラインダンスをするようにそろって同じ動きをする。入り交じった複雑な形をしているが、慣性重力波の位相速度は波数ベクトル(k,l)の絶対値だけに依存して、その方向にはよらない。波動方程式には元々その方向による依存性は無い。そのため波数ベクトル(k,l)=(k,0)として、波の伝播方向をx軸と一致させても解の一般性を失わない。x軸方向に伝播する波は、(35)、(36)式でl=0と置いた式になるが、この様に簡単化しても一般性を失わない。このとき当然α=β=0となるが、さらに(1/k)(H/D)は、波の波長λ(=2π/k)、振幅H、水深Dが定まると、一定値と見なすことができるのでC(=定数)と置ける。理解しやすくするために(x,y)=(0,0)での(u,v)の時間的変化を考察する。その場合、上式はさらに簡単化できてとなる。f0Cの楕円を描くことを示している。この運動を、波数ベクトルが(k,l)となる座標系(x’,y’)から見たときの流速ベクトル(u’,v’)に変換する。(u’,v’)は、元の座標系での流速ベクトル(u,v)と以下の関係で結び付けられる。これらの関係式を用いるととなり、(35)、(36)式と一致することが証明できる。2.流速解の意味(35’)、(36’)式の意味するところをまとめると以下のようになる。流速は波高の振幅に比例し波高が大きな所ほど速い。波高が負になれば流速は反対向きに速くなる。コリオリ力が無視できる(f0→0)場合には、普通の重力波の流速表現に一致する。流速の変化幅は、波長λ(=2π/k)が長くなると大きくなる。波長が長いと、同じ周期で同一の波高変化幅を実現するためには、流体はより急速に移動しなければならないからです。波長λ→∞の場合、伝播速度の式のルート記号の中の1は省略できて伝播速度μ〜f0λ/2π→∞となる。これでは意味をなさないので、元のσの式に帰るととなる。λ→∞(振幅の変化がない)場合、流体の運動は単なる慣性振動になることを示している。波長が無限大でも速度成分を持つと言うことは、元々の流体が全体的に平行移動していたと言うことです。波の伝播方向に平行な速度成分uがσに比例することは、σが周期T(=2π/σ)に関係することから了解できます。周期が短い(σが大きい)ほどuは大きくならないと、同じ波高振幅Hの変化を実現できないからです。波の伝播方向に垂直な速度成分vがコリオリ力効果f0に比例して生じることが、慣性重力波で新たに出てきた事柄です。コリオリ力が存在するために流体は波の伝播に伴って楕円運動することになる。場合をアニメーション表示する。波高は黄色(高)青色(低)で表現し、流速ベクトルは赤矢印で示す。このアニメーションから以下のことが読み取れる。x方向の速度成分uの加速度は、水面勾配(圧力勾配)に関係する。y方向の速度成分vの加速・減速には表面勾配(圧力勾配)関係していません。y方向の速度成分vの変化ははコリオリ力によって生じます。特定の地点の速度ベクトルの時間的変化に着目すれば、柱状流体要素は波の伝播に伴って右回り(北半球)の楕円運動をすることが解る。波の伝播速度に比べて流体速度は遙かに小さな値です。柱状流体要素はほぼ同じ位置で小さな楕円運動をする。一回転に要する時間は波の周期に同期している。慣性重力波の周期はですから、常に“慣性振動”の周期(2π/f0)より短くなります。ちなみに、慣性振動の周期はT’=2π/σ=2π/f0=12h/sinθ0で、緯度によって定まっていました。“慣性重力波”の回転運動は波高の変化に伴って生じるものですが、元々の流体が全体的に平行移動する速度を持っていたら、“慣性振動”の回転運動が重なって生じます。速度ベクトルの変化を観察すれば解るように隣り合う領域の速度ベクトルとの関係は、単なる平行流の関係ではありません。隣り合う部分との速度差が増減して互いにずれあう様な関係が生じる。これによって波高の変化と共に柱状領域の流体要素が持つ“渦度”が変化する。“渦位保存則”(25)が成り立つためには、“渦位”のコリオリ力項f0が一定でも分母のD+hが変化するので、“渦度”(∂v/∂x-∂u/∂y)が変化しなければならない。着目する柱状領域の“渦度”は波高の最高点で最大、最低点で最小です。渦度の変化率は表面勾配が最も大きな点で最大です。f0を単位として表した角振動数を縦軸に、f0/(gD)0.5を単位として表した波数を横軸に取ったグラフで、“慣性振動”、“重力波”、“慣性重力波”の関係を示す。f0/(gD)0.5は次節で説明する“ロスビーの変形半径”の逆数です。以下は、この図についての補足説明です。“慣性重力波”の角振動数は、波の伝播に伴って柱状流体要素が水平楕円運動をする角振動数の意味ですが、それは水面が上下動する角振動数でもあります。“慣性振動”が慣性重力波の曲線の中に存在することの意味はなかなか解りにくいところですが、波数→0(波長→∞)で水面の上下動は無いが柱状流体要素は水平方向の回転運動していると言うことです。元々の流体が全体として平行移動する速度成分を持っていればそれは可能です。“重力波”の直線は、波長と周期からきまる波の伝播速度μ=λ/Tが常に一定であることを示している。波長が変化しも伝播速度は常に一定で分散性が無いことを示している。“慣性重力波”の曲線はコリオリ力項f0→0とともに“重力波”の直線に近づきます。そのとき慣性振動の振動周期もT→∞となりますので当然σ=2π/T→0となり、“慣性振動”を表す赤丸もσ/f0=1を保ったまま原点に近づきます。グラフから明らかなように、慣性重力波の振動周期T=2π/σは同じ波長の重力波の振動周期よりも必ず短くなる。これはコリオリ力が復元力の一部となって働くためです。振動周期に最大値Tmax=2π/f0があることも慣性重力波の特徴です。負の波数というのは解りにくい所ですが、負の方向に伝播する波の波数であると言うだけです。[補足説明]ここでは、波長や周期を定めるメカニズムについては何も議論していません。潮汐波ならば月や太陽の起潮力の周期が波の伝播速度(gD)0.5と一緒になって波長を決めました。波動方程式の中にはそういった波長や周期を決める原因になるような項はどこにもありませんから、その点に付いては議論のしようが無いからです。仮に何らかの原因で海面や圏界面に特定の波長・周期の擾乱・攪乱が生じて波が発生したら、その波動は(30)式で与えられる位相速度を持って伝播すると言うだけです。性質から入って行かざるをえないからなのでしょう。おけるケルビン波」で用いたものと全く同じです。ケルビン波も“慣性重力波”の一種です。慣性重力波の波動方程式は、本来その解の中に慣性振動を含んでいますから、流体粒子の運動方向は曲げられます。ところが、側壁境界面がある海洋の場合、境界が運動をブロックするために慣性振動(流体の回転運動)ができなくなる。x軸に沿った側壁が有る特殊な境界条件のときには、f=f0+βy≒f0=一定と置いた方程式系(23)、(24)、(3)で、最初からv=0と仮定した方程式系に簡略化できる。それが別項「ケルビン波」3.(3)3.[補足説明1]に示した(A)、(B)、(C)式です。このとき(B)式は側壁に沿った方向の流れuに関して地衡流が成立していることを示している。それは特殊な境界条件(v=0)のために出てきたものです。そのため、側壁方向に関しては分散性の無い“重力波”が解となる。慣性重力波は、本来その解の中に重力波も含んでいるのですから、この様な重力波と見なせる解が生じるのも頷けます。そのときコリオリ力は慣性振動的な回転運動を生じるのではなくて、別稿「ケルビン波」3.(3)3.[補足説明2]で述べたようにy方向の波高変化を生じさせ、それとともに渦度を変化させるメカニズムとして効いている。“慣性重力波”が湾や海峡に進入していくとき、それが側壁境界の影響を受けてどのようにして“ケルビン波”に変化していくかについて、スヴェズドラップが文献1.のp61〜62で説明していますのでご覧ください。2.ケルビン波の特徴ケルビン波は別稿で詳しく論じましたが、ここのテーマに沿ってその特徴をもう一度まとめます。側壁(x軸に平行)の場所にとどまらずすべての領域でv=0となり、壁に平行な速度成分uしか存在しない。波高振幅は壁の位置で最大で、壁から離れるにつれて指数関数的に減少する。流速成分uも壁から離れるにつれて指数関数的に減少する。波高hや、流速uの大きさは、壁から距離y=(gD)0.5/f0だけ離れると、壁での値の1/e倍に減少する。deformation)”という。この中にコリオリ力項f0が含まれていることから明らかなように、y方向の振幅や流速の変化はコリオリ力が決めている。コリオリ力は、流体粒子に慣性振動的な回転運動を生じさせるのではなくて、波高や流速をy方向で変化させる所に働いている。8.文献4.赤道ケルビン波ここで丁度良い機会なので、近年注目されている“赤道ケルビン波”について説明します。それは、f=f0+βyに於いてf0=0と置けるがβ項が無視できない場合です。コリオリ力の緯度依存性と圧力傾度力の両方を考慮するのだが、コリオリ力f0は無視できる場合です。赤道にケルビン波が存在することは、ケルビン波のメカニズムを考察すれば推測できる。別稿で説明したケルビン波はx軸(東西)方向に平行な壁がある場合の解でした。この場合にはv=0という解が得られた。地球の赤道上では、赤道を挟んでコリオリ力が働く方向が逆転するために、壁が無くても、赤道をはさんで南北両半球で対称(すなわちv=0)のケルビン波が存在できそうです。事実そのような解が存在することは次のようにして示すことができる。低緯度に行くほど小さくなり赤道上では0になってしまうことです。それで赤道付近の運動を論じるには、1.(1)2.で説明したコリオリ力の緯度依存性を含んだβ平面近似を用いる必要がある。この様な近似を赤道β平面近似という。f0とβの緯度による変化の様子は1.(1)2.を復習してください。ここでは、赤道上にx軸を取る。この近似を使い、かつ初めからv=0とおくと、赤道付近のケルビン波を表す方程式としてが得られる。おけるケルビン波」の3.(3)3.[補足説明1]、あるいは前章3.(5)1.、で示した(A)、(B)、(C)式に似ていますが、(38)式の左辺の係数がyの関数になっている所が違ってる。そのため解の様子が普通のケルビン波と異なります。(37)式の両辺をtで偏微分した式に、(39)式の両辺をxで偏微分した式を適用するとが得られる。これはx方向の一次元波動方程式で、x方向(赤道方向)に関しては位相速度がμ=±(gD)0.5で伝播する“重力波”となることを示している。これは普通のケルビン波と同じです。2.南北方向の変化(40)式を考慮すると、振幅がy方向でのみ変化するの形の解が考えられる。これらを(37)式と(38)式に代入すれば、H(y)とU(y)の具体的な形を決めることができる。ここでH0はy=0(すなわち赤道)におけるh(x,y,t)の振幅H(0)です。y→∞でh→0となるためには、位相速度μ=σ/kは正でなければならない。赤道ケルビン波は常に東に向かって進む。これはコリオリ力と水面傾度(圧力傾度)の釣り合いを考えると当然の結果です。この結論を用いると求める解は以下のようになる。波は、赤道に沿って常に東に向かって伝わり、位相速度は重力波と同じ(gD)0.5で分散性は無い。赤道から高緯度へ離れるにつれて、波高振幅や流速振幅は減少する。赤道での値の1/e倍になる特徴的な距離はで与えられる。赤道付近の大洋深度をD=4000mとして、上記の特徴的スケールを計算すると4180≒約4000km程度になる。大気の場合には赤道圏界面の高度D=17000mとすると、特徴的スケールは約6000km程度になる。おけるケルビン波」5.(2)4.[補足説明2]で、潮汐現象に伴う赤道ケルビン波の存在の可能性を述べましたが、この場合には潮汐周期(12.42h)が水深Dと相まって波長9000kmを決めます。赤道の成層圏大気に存在すると考えられている赤道ケルビン波(位相速度約25m/s、周期10〜20日、波長約30000〜40000km)については、また別のメカニズムが必要です。Sci.,35,p714〜722,1978年、ネットからダウンロード可]。[補足説明2]今日エルニーニョ現象に絡んで存在の可能性が指摘されている赤道ケルビン波は、海洋表層水の温度分布変化に伴う内部密度境界層に生じるものです。そのため、gやDは内部波の理論に基づいたものに置き換える必要があります。位相速度や特徴的スケールもそれに応じた値になります。そのため伝播速度は3m/s程度となり、約3ヶ月かけて太平洋を西から東へ横断すると考えられている。この赤道ケルビン波の伝播に引き続いて発生する赤道ロスビー波[6.(3)2.参照]の太平洋横断(東から西へ)期間は約9ヶ月と考えられています。これと前記の3ヶ月を足し合わせたものを二倍して得られるエルニーニョ周期(遅延振動子理論)は2年程度となります。ところが、この周期は実際のエルニーニョの2〜6年周期と合わない。いう新しい理論がより広い支持を受けているそうです。[参考文献11.第10章参照]2.赤道ケルビン波のアニメーション話を解りやすくするためにk、σではなくて、波長λ、周期Tに書き直した式にする。したをアニメーション表示する。表示の中心線が赤道で、波は西から東へ向かって伝播する。上図に速度ベクトル(赤矢印)を重ね合わせると下図になる。伝播速度に比べて流体速度は遙かに小さな値です。柱状流体要素はほぼ同じ位置で赤道方向に沿った小さな往復運動をします。このとき、波高の変化に伴って隣り合う速度成分との間で互いの大きさが変動する。これが波高の変化に伴って渦度が変化することを示しています。赤道に沿う方向の流れでは地衡流が実現されています。これらは普通のケルビン波と同じです。波動の規模が地球規模になってくると、地球が球形をしていることによるコリオリ力の緯度変化(β効果)が無視できなくなります。ここで行うのはまさに地球規模の波動(プラネタリー波)の考察です。3.章の波動はf=f0=一定の場合、4章の波動はf=βyと表される場合でした。ここでいよいよ、中緯度地帯で最も重要な“ロスビー波”について説明します。それは、f=f0+βyに於いてf0もβ項も無視できない場合です。コリオリ力と圧力傾度力の両方を考慮し、コリオリ力の緯度依存性も考慮する場合です。(1)方程式1.方程式系ここで取り扱う方程式系をあらためて書いておく。前記の方程式系では、係数が変数yの関数となる項が(1’)(2’)(5’)のあちらこちらに含まれている。一般に線形偏微分方程式は、その係数がすべて定数であると解が求めやすいが、そうでない場合には容易に解けない。そのためβ効果に関する効果として(5’)式の真ん中の項βvだけ残して、それ以外のf=f0+βyはすべてf=f0=一定で置き換える。このとき、赤道付近ではf0とβyが同じオーダーとなるので、この近似は使えず、新たに考え直さなければならない。高緯度ではβ→0となるので、やはりこの近似は成り立たない。そのためここでの議論はあくまで中緯度に限定される。赤道地帯に生じる波については6.赤道波で説明します。このように近似すると方程式の係数はすべて定数となり、解はたいへん求めやすくなる。結局のところが解くべき連立方程式です。一つの未知関数のみを含む、単一の式に変換する。“連続の式”(46)の両辺をtで偏微分したものに、“運動方程式”(44)の両辺をxで、(45)の両辺をyで偏微分したものを代入する。この両辺をtで偏微分して、“渦度方程式”(48)を代入する。“運動方程式”(44)式の左辺の時間微分項が無視できる地衡流近似を用いるとが得られる。この中の左辺第3項のf02の項は自由表面の渦管の伸縮の効果を現している。第4項がβ効果に依存する項です。“連続の式”(46)の両辺をxで偏微分したものに、“運動方程式”(44)の両辺をtで偏微分したものと、“運動方程式”(45)を代入する。この両辺をtで偏微分したものに、“渦度方程式”(48)の両辺をyで偏微分したものを代入する。“連続の式”(46)式の右辺の時間微分項を無視した式を用いる。これは波高の時間的変化が無視できるとする地衡流近似です。同様にして、“連続の式”(46)の両辺をyで偏微分したものに、“運動方程式”(45)の両辺をtで偏微分したものと、“運動方程式”(44)を代入する。この両辺をtで偏微分したものに、“渦度方程式”(48)の両辺をxで偏微分したものを代入するとが得られる。このときコリオリ力項をf=f0に近似した式ではなくて、f=f0+βyをそのまま用いても上記の式変形は可能です。運動方程式のtによる偏微分や渦度方程式のxによる偏微分にはその違いは影響しません。最後の変形に於いて地衡流近似を用いていません。vに関する方程式(52)は、[5.(1)2.で述べた近似]や[地衡流近似]を用いなくても導けます。(52)のf2の項が定数係数で無くなるので、実際に解くときには結局5.(1)2.の近似を用いなければなりません。結局、三つの量h、u、vは全て同じ形の“定数係数の3階線形同次方程式”(50)、(51)、(52)を満足する。これらは3.(2)で説明した“慣性重力波の波動方程式”に、コリオリ力の緯度依存性を表すβ項が付け加わった形をしている。equation)と呼ばれている。3.波高hと流速(u,v)の関係ここで、波高hと流速(u,v)の間に成り立つ関係式を求めておきます。運動方程式(44)と(45)の両辺をtで偏微分するととなる。(44’)式に“運動方程式”(45)を代入するととなる。(45’)式に“運動方程式”(44)を代入するととなる。今用いている近似の範囲では慣性重力波の場合と同じ関係式が使える。これらは波高解hから流速解(u,v)を導くときに使いますが、“ロスビー波”も、波の伝播に伴った楕円運動をすることが解る。5.(5)3.で説明するように“非発散性ロスピー波”の平面波解の場合には波面に平行な往復運動になる。最初に最も簡単な形の波高解を考察する。これが“準地衡流方程式”(50)を満たしていることは代入してみれば直ちに確認できる。そのときσはを満足する場合に限定される。これはσに関する三次の代数方程式だから可能な解は3つある。三つの解の様子は、例えばY(X)=−X3+20X+5のグラフの様子を考えれば解る。このグラフはY=−X3のグラフにY=20Xの一次関数のグラフを足して、それをさらにYの正方向へ5だけずらしたグラフとなる。−X3+20X+5=0の解は、X=0付近に一つ[20X+5=0からほぼ求まる]、X軸の正負にほぼ同程度離れた二カ所に二つ[−X2+20=0からほぼ求まる]の合計三つの解を持つ。今の場合に当てはめると、σ≪k(gD)0.5の場合と、σ≫β/kの場合を仮定して代数方程式を解けばよい。2.σ≫β/kの場合これはβ項の効果が無視できる場合です。このとき(56)式のβgDkの項が省略できてとなる。これは3.(3)1.で求めた慣性重力波の角振動数の表現式(32)式と同じです。これはコリオリ力の緯度依存性が省略できる(β→0)場合には、解は“慣性重力波”になることを表している。β→0の場合には、準地衡流方程式は3.(2)で説明した慣性重力波の波動方程式となるわけですから、この様な解が含まれているのは当然です。今は中緯度に於ける波動を考えているので必ずしもβ→0では有りません。βが小さくなくてもkが有る程度大きければσ≫β/kを満足するσ領域がありますから、その様な状況の慣性重力波が存在するということです。そのとき、同じkに対してσ≪k(gD)0.5を満足するσ領域には次項で説明するロスビー波の解も同時に存在できます。これら条件不等式の意味は解りにくいところですから、5.(4)5.[補足説明2]のグラフをご覧になりながら、その意味を検討してみてください。k→0となるとσは重力波から離れて原点から有る程度離れたσminに収束しますから、σ≫β/kを満足できるのです。このとき、波面(位相が一定の場所)が進む速度(位相速度)は3.(3)1.の(34)式で求めたようにとなり、重力波の位相速度(gD)0.5より常に大きくなりますが、その大きさ自体は同程度です。緯度45°における深さ4000mの海洋の場合には、f0=10-4rad/s、程度ですからとなります。波長λが例えば潮汐波の10000kmでも、この項は0.6程度ですから、μは重力波の位相速度(gD)0.5〜200m/sとそんなに違いません。“慣性重力波”の伝播の様子や流速の様子は、3章を復習されてください。3.σ≪k(gD)0.5の場合これは位相速度が重力波のそれに比較して小さい場合です。このとき、(56)式の−σ3の項が省略できてとなる。このとき必ずσ<0となるので、ロスビー波の位相は必ず西進します。波のエネルギーが常に西進するとは限りません。ここはロスビー波の議論で最も解りにくいところですから5.(4)節で詳しく説明します。ロスビー波の位相速度はとなる。ロスビー波は位相速度が波長に依存する分散性の波です。ここで注意すべきは慣性重力波と違ってロスビー波の位相速度を決めているのはコリオリ力の緯度変化に関係するβ項です。慣性重力波の場合は重力に関係する(gD)0.5が大半を決めていた。コリオリ力項f0と波長λはその偏差に関係しただけです。ロスビー波では、位相速度にf0や波長λは関係してはいますが、その大半を決めているのはβ項です。そのためロスピー波と慣性重力波の位相速度は全く異なります。実際、緯度45°における深さ4000mの海洋の場合、f0=10-4rad/s、β=1.6×10-11rad/m・s、Rd〜2000km程度ですからとなる。慣性重力波の200m/sに対して25m/sですから、ロスビー波の方が一桁小さい。このことから最初に用いた近似が良い精度で成り立っていることが確認できる。一つ注意すべきは、ロスビー波の位相速度μは波長λが同じでも、波面の伝播する方向によって変化することです。位相速度μの表現式に掛かっているは波面の伝播方向によって位相速度が変化することを意味する。このようになったのはσの表現式にkが入っているためですが、それは元々コリオリ力が緯度と共に変化することに起因する。コリオリ力の緯度依存性β項のために位相速度が伝播方向によって変化する。これはロスビー波の大きな特徴です。4.中緯度に存在する波ここでも波を生じるメカニズムについては何も論じていませんが、中緯度ではf0とβ項は常に存在しますので、方程式の線形性から明らかなように慣性重力波とロスビー波は両方同時に重なって存在できます。そのときσの違いから慣性重力波とロスビー波は明瞭に区別できる。l=0の場合には[東進する慣性重力波]、[西進する慣性重力波]、[西進するロスビー波]の三つの波が同時に存在します。l=0の場合について上記の関係をグラフ表示してみる。ロスビー波の場合、波数ベクトルの大きさは同じでも方向が違えばσの値が変化しますので、3.(4)3.で示した慣性重力波の様に波数ベクトルの絶対値(k2+l2)0.5を横軸に取ることはできません。(k,l)=(k,0)の特別な場合のグラフをkを横軸に取って表示しています。このグラフのσ<0は位相速度が西向きになることをあらわしているに過ぎません。σを常に正として、kの正負で東進・西進の方向を表すこともできます。すると、上記のグラフはの様になる。これが良く見かけるグラフです。これらのグラフから明らかなように、中緯度に於いては位相速度が(gD)0.5の重力波を境にして、速い位相速度を持つ慣性重力波と、遅い位相速度を持つロスビー波がはっきり区別できる。このとき、流速解についても慣性重力波とロスビー波に対応するものが存在する。両者が重なって存在すると考えればよい。“慣性重力波”に伴う流速解は3章で詳しく調べました。“ロスピー波”に伴う流速解は5.(2)3.で説明したように基本的に慣性重力波と同じですが、σが異なることから流速変動の振幅と周期が異なっています。“非発散性ロスピー波”の場合は5.(5)で検討します。8.文献(4)ロスビー波の位相速度と群速度1.位相速度と波数ベクトルの関係ロスビー波の[位相速度が波数ベクトルの方向に依存して変化する]ことや、[位相速度の西向きの成分は波数ベクトルの絶対値が同じなら波数ベクトルの方向によらず一定になる]ことは解りにくいところなので、アニメーションで説明します。(55)式に(57)式を代入するととなる。ここで位相速度μや周期Tが負になっているのは、波が波数ベクトルの方向と逆向きに伝播するからで、それ以上の意味は有りません。(C)(k,l)=(1/3,80.5/3)このとき波数ベクトルと波面との関係は下図の様になっている。これらのアニメーションから、[波数ベクトルの方向が変わると位相速度がkに比例して変化する]ことと[ロスビー波は波数ベクトルの向きとは逆に進む]ことを読み取ってください。このときさらに、波面がx軸上を西に進む速度はどの場合にも同じになります。[ロスビー波の波面がx軸に対してどのように傾いていても波数ベクトルの絶対値が同じで有る限り西向きの速度は同じ]になります。なることは、上図の赤半円に内接する直角三角形のkが位相速度に比例する量であることに注意すれば明らかです。この様な特別の性質を持つのは、ロスビー波がβ効果に伴う波だからです。2.角振動数σと波数kの関係ロスビー波の理論で最も解りにくい所は、位相が常に西に向かって進むと言うことです。このことを正しく理解するには、波のエネルギーの伝播とは何かを理解する必要があります。それは波の位相速度と群速度の違いを理解することですが、それを説明するために5.(3)4.で説明したl=0の場合についてのσとkのグラフをもう一度取り上げます。そこのグラフのロスビー波の部分だけを取り出して、σを拡大して描くとのようになる。このとき以下の事柄が言える。ロスビー波の角振動数には上限σmax(振動周期に下限Tmin=2π/σmax)があり、振動数がそれ以上高い波は存在できない。ちなみに、慣性重力波については振動数に下限が有る。角振動数が極値を取るのは、波数kがロスビーの変形半径Rdの逆数に一致するときです。波長λがロスビーの変形半径の2π倍のときです。角振動数の最大値σmax=βRd/2は、Rd=(gD)0.5/f0ですから、β値に比例しf0値に反比例します。波長がロスビーの変形半径の2π倍を境にして、それよりも波長が長い場合には位相速度ベクトルと群速度ベクトルは同一の方向を向くが、それよりも短い場合には、位相速度ベクトルと群速度ベクトルは互いに反対方向を向く。この様なことが生じるのはロスビー波だけてあって、重力波や慣性重力波では生じない。今までの議論では波を発生させるメカニズムについては何も話していませんでしたが、ここで以下のような状況でロスビー波を発生させる。β平面のx=0の位置にy軸に平行な垂直板を置いて、その板をx軸方向に角振動数σ0でゆっくりと振動させる。角振動数が定まった波をx=0の位置で発生させる。そのとき、振動板で作られた波は当然x軸の正方向と負方向にエネルギーを伝えていくはずです。ところが、今まで述べてきたように波の位相は必ず西向きにしか進みません。それではエネルギーが東と西の両方に伝わるとはどういうことなのでしょうか。その鍵は上記のグラフの中にあります。上記のグラフでσが定まった波を発生させた場合、発生する波の波数には二つの解が存在します。図中のk1とk2です。このとき、いずれのkの値のときもσ<0ですから、今まで説明したようにこの二つの波の位相速度は負となり必ず西に進みます。このとき、波の波数が定まっていると言っても実際に生じる波は、例えばk1の波ならばεを微少な値としてk1付近の(k1−ε)〜(k1+ε)の波数の波が重ね合わされたものと考えるべきです。k2の波についても同様です。重ね合わせる波の振幅Hは当然kの関数ですからH(k)と書かれています。実際に発生する波は非常に小さな幅2εを有し、平均の波数がk1ork2の多数の波が重ね合わされたものであると考えるべきです。そのときと置くと、前記の式はとになります。]括弧内の被積分項目は、(k−ki)、(σ−σ0)が小さな値なので、正弦・余弦関数の振幅が空間的にも時間的にも非常にゆっくり変化します。これはちょうどアナログのラジオ放送電波の“搬送波”とそれに重ね合わされた“変調波”の関係に似ています。電波の場合は波長が変わっても伝播速度(光速)は変わらない非分散性の波なので、搬送波と変調波は同じ速度(光速)で伝わりますが、今は、波長により伝播速度が異なる分散性の波なので、“変調波”の伝播速度は“搬送波”の伝播速度とは異なってきます。そのとき、波のエネルギー輸送を担っているのは変調波の伝播速度です。その様に、搬送波の振幅の変化(変調波)が移動する速度を“群速度”といいで表される。搬送波そのものが伝わる速度は“位相速度”といわれで表される。[補足説明]位相速度と群速度の概念は解りにくい所です。参考文献13の§26に解りやすい説明がありますが、より直観的な理解には別稿「ファインマン物理第�U巻の§6-3の説明が有益です。[光の波]と[流体を伝わる重力波]の違いはありますが本質は同じです。流体を構成する部分要素が(光波の場合と同様に)入射してきた一次波によって二次波を生じ、それが重ね合わされたものが進行波となる。そのとき二次波の位相の遅れに関係するのがコリオリ力の効果です。その効果の関わり方の違いが分散性を生じ、位相速度と群速度の違いを生み出す。波のエネルギーが位相速度ではなくて群速度で伝わるというのはなかなか解りにくいところなので、前記のグラフのk1とk2の近傍の二カ所の波を重ね合わせた合成波を用いて説明します。k1の近傍の二カ所(k1’,σ1)と(k1,σ0)の波を重ね合わせる。k2の近傍の二カ所(k2’,σ2)と(k2,σ0)の波を重ね合わせたものも全く同様にしてとなる。ここで具体的な数値を下記の様に取る。これら数値を用いると、重ね合わせた波高解は以下のようになる。k1の波はk2の波はのようになる。これらのグラフの横軸はそれぞれ原点からの距離(正が東側、負が西側)を表しますが、アニメーションが見やすいようにk1波とk2波では縮尺スケールが変えてあります。これらの例から以下のことを読み取れる。搬送波(波長小)と変調波(波長大)の伝播する速度が違うが、波のエネルギーの伝わる速度は変調波の伝わる速度(群速度)であることを波の伝わる様子から読み取れる。変調波の節の移動速度を超えて波のエネルギーは移動できないのですから。k1の波は原点の造波装置から発生して西に向かうロスビー波を表している。位相速度も群速度もいずれも西に進んでいる。k1の波については西進する位相速度の方が、西進する群速度よりも少し速い。それらの値はすでに示したグラフ中の傾きに対応する。そのとき、搬送波の波長が長くなる(k→0)につれて位相速度は大きくなるのですが、慣性重力波の場合と違って有限の値にとどまる。k→0とともに群速度と位相速度の差はだんだんなくなる。これらはσ−kグラフの傾きの変化から読み取れます。k2の波は原点の造波装置から発生して東に向かうロスビー波を表している。そのとき位相速度は西進するが、エネルギーを伝える群速度は東進している。位相速度と群速度が逆方向であることはとても解りにくいところですが、アニメーションから了解できる。そのとき、搬送波の波長が短くなる(k→∞)につりて位相速度は小さくなるのですが、そのとき群速度はいったん増大した後にだんだん遅くなることがσ−kグラフの傾きの変化から読み取れます。k1とk2のいずれの場合も、搬送波の位相速度は常に西向き(ロスビー波の特徴)だが、エネルギーを伝える変調波はそれぞれ造波装置から外向きに向かっている。そのとき、k1とk2がσが極値(σ=σmax)をとるときの値1/Rdに近づくにつれて、群速度はゼロになる。同じ振動数の造波装置から作られた波であるにもかかわらず、西と東では波長の異なった波が生じる。普通は搬送波の位相速度と変調波の群速度は、例え両者の速度が異なっていても、同じ方向を向いています。重力波や慣性重力波がそうでした。ところがロスビー波では位相速度の方向と群速度の方向が反対方向を向く場合があります。それはσとkのグラフで位相速度を表す直線と群速度を表す直線の傾きが正負に分かれる場合があるからです。k>1/Rdの領域がそうです。慣性重力波では、位相速度を表す直線と群速度を表す直線の傾きが正負に分かれることが生じないことを5.(3)4.のグラフから読み取ってください。Phillipsは、この様な波が実際に生じることを巧みな実験で示した。次項でその内容を説明します。それぞれ異なっており、それは理論の予想する通りであることを、次のような模型実験で示した[参考文献3.と4.]。彼らは図に示す様な、同心円筒の間に流体をいれ、共通軸の回りに回転させた。このとき、容器の底は平坦ではなくて、半径が増大すると半径方向に水深が深くなるようにしてある。このようにすると、1.(2)3.で述べた渦位保存則に於いてコリオリ力項f=f0=一定であるが、水深項DがD−β’yと変化しているために、β効果が擬似的に実現できる[詳細は5.(6)参照]。この円形水路は直径2m程度、水路幅0.3mのかなり巨大なものです。水深は内側円筒壁部で6cm、外側円筒壁部で12cmです。この水路は4.2秒で1回転する。造波装置の板の振幅は5cm程度で、振動周期は10s〜30s程度の範囲で変化させた。造波装置によって発生する波の波長や位相速度を測定するために、表面張力を押さえる薬剤を表面に滴下したのちにアルミ粉末を浮かべてある。波は水路と同期して回転するように水路上方に設置された、カメラと動画撮影用の映写機により記録された。1965年にPhillipsが得た結論を文献3.から引用するとの様になった。グラフ中のsが波数で、波数s=40がおよそ波長λ=11cmに対応する。1967年の論文(文献4.)から引用するととなった。得られた結果は、振動版の東(右)側にできる波も、西(左)側にできる波もいずれも西進(左進)した。造波板の東西で波長の違う波が観測された。side)が長く、それらは理論的に予測されるものと一致した。5.二次元波数空間に於ける位相速度と群速度5.(4)2.や3.はl=0の一次元の波数空間kに関する議論でしたが、二次元波数空間(k,l)に拡張した場合を説明します。(k,l)座標上にσ曲面を描くと下図の様になる。この図に置いて、三角形OCDに於ける線分ODの傾きが波速ベクトル(k,l)方向の位相速度を表す。同様に三角形BCDに於ける線分BDの傾きがx方向の位相速度、三角形ACDに於ける線分ADの傾きがy方向の位相速度を表す。そのとき、位相を表す波面は線分ベクトルOCに垂直で、波面の進行方向は線分ベクトルOCの逆向き(線分ODの傾きは負だから)となる。その位相波の群速度はどうなるのだろうか。前々項の議論から明らかなように、D点に於けるσ曲面の最大勾配の大きさが群速度を表す。群速度ベクトルの方向はσ曲面の最大勾配の方向です。群速度ベクトルはσ曲面のD点に於けるσの等高線(図の赤半円)に垂直です。そのとき、エネルギーの伝播に関係する変調波の波面は群速度ベクトルの方向に垂直で、変調波の波面は群速度ベクトルの方向(D点の勾配ベクトルは正方向に向く)に進んでいく。位相速度ベクトルと群速度ベクトルの方向が異なることは理解しにくい所なので、具体的な波の例について時間的変化を追跡してみる。(k,l,σ)での話だから非常に解りにくい。上図のD点の状況をxy座標・波高hの実空間(x,y,h)の値で解釈し直す。搬送波の波長λp、振動周期Tp、位相速度μpはとなる。実空間でのエネルギー輸送速度を表す変調波の伝播速度(群速度)のx成分とy成分は、で与えられた。これは変調波の波面の実空間での伝播速度成分を表すので、実際の群速度の大きさはで与えられる。実空間において、変調波の波面は搬送波の約0.065/0.058=1.1倍の速度で移動する。そのあたりはグラフ曲面の傾きからも読み取れます。(k,l,σ)に於ける変調波の波数ベクトル(Δk,Δl)で表されるが、その波数と変調波の角振動数Δσは、搬送波と同様な関係式を満たさなければならない。このときΔσと波数ベクトルの大きさ(Δk2+Δl2)0.5は、両者の比率が上記の関係を満たせば任意に取れる。アニメーション表示したときに解りやすいのは、変調波の波長λgが搬送波の波長λpの10倍程度の場合です。その様にするにはΔσをいくらにしなければならないか見積もってみる。変調波の波長を搬送波の波長の10倍とすると、変調波の角振動数Δσは程度にすればよいことが解る。このときΔσを0.0234ではなくて0.021とすると、丁度この10倍が搬送波の角振動数0.21に一致する。このような整数倍にする必要性など全く無いのだが、アニメーションファイルを構成するときに整数比にしておいた方が都合が良いのでΔσを0.021として、波長の方を調整することにする。Δσ=0.021とすると波長はとすればよいことが解る。ここで変調波の波面が群速度ベクトルに直交するように0.3231をΔkとΔlに分割するにはとすればよいことが解る。搬送波の位相速度、変調波の群速度が上記の値となる波は以下の式によって表される。この波の様子を以下でアニメーション表示する。変調波の波長λgは19.45となりますが、アニメーション中の変調波の山と山の間隔はその半分の9.725であることに注意してください。λgx、λgyについても同様です。なることは5.(4)3.のアニメーションを復習すると解ります。そこの灰色カーブが変調波の波形を表しています。アニメーションから、この例の位相速度ベクトルは波数ベクトルとは逆向きであること、群速度ベクトルが波数ベクトルの方向と異なっていることなどがわかる。群速度が位相速度の約1.1倍になっていることも読み取れます。[補足説明1]5.(4)1.で、[ロスビー波の波面がx軸に対してどのように傾いていても波数ベクトルの絶対値が同じである限り西向きの速度は同じ]になることを説明しましたが、σ曲面の上で考えるとこのことが直ちに理解できます。下図中の赤ラインの傾きがすべて同じであることが西向きの位相速度が同じになることを示しています。[補足説明2]5.(3)4.で説明した角振動数σと波数(k,l)の関係を三次元グラフで示す。簡単化のためにf0=gD=Rd=1、β=0.5と置いている。それぞれの波のσ曲面は下図の様になる。これらの曲面について、ロスビー波のσ曲面と同様な考察をすると、以下のことが直ちに了解できる。“重力波”は、波数ベクトルの大きさ・方向が変化しても位相速度ベクトルと群速度ベクトルの大きさ・方向は変化せずに常に一定値となる。位相速度ベクトルと群速度ベクトルは常に同じ大きさであり両者の方向も常に一致する。“慣性重力波”は、波数ベクトルがゼロのときには位相速度は無限大ですが、波数ベクトルの大きさが増大するにつれて位相速度は無限大からだんだん減少してきて重力波の位相速度(=群速度)に近づいていく。群速度の大きさは波数ベクトルがゼロのときにはゼロだが、波数ベクトルの大きさが増大するにつれてやはり重力波の群速度(=位相速度)に近づいていく。そのとき、σ曲面はσ軸を中心軸として回転対称だから、これらの傾向は波数ベクトルの方向に依らない。σ曲面を比較してみると、“ロスビー波”の特異な性質はすべてのそのσ曲面の特殊な形状に由来していることが解る。8.文献(5)非発散性ロスビー波ロスビー波が重力波や慣性重力波と根本的に違う所は、ロスビー波は完全に二次元の流体でも存在できることです。海面や圏界面に蓋をしてしまって境界面の上下動を無くしてもロスビー波は存在できる。この場合には、H=0だから方程式(50)やその解(55)は意味をなさない。流速場に関する(51)、(52)式や、(50)式のhを圧力場と見なした式はそのまま成り立つ。流速場や圧力場がロスピー波となって伝わると考えればよいのです。wave)と呼びます。このとき慣性重力波が重なって存在する事はなくてロスピー波のみが存在する。その為、ロスビー波の性質を調べるには特に適した状況です。1.非発散性ロスビー波の方程式系この場合の方程式系は5.(1)1.の(1’)、(2’)、(3)式でhが絡む項をすべて省略したものとなる。解くべき連立方程式です。ここでは(44’)と(45’)の右辺に圧力項を残していますが、今はこの圧力項が波高に比例するわけではありません。(u,v)の水平方向での変動を実現するために残してあるだけです。(44’)と(45’)の右辺の圧力項をゼロで近似してしまえば、これは、2.(2)1.で説明した慣性振動を表す運動方程式にβ項が付け加わった形ですが、全体的な平行流があるとき慣性振動が重要になるのであって、慣性振動が解として注目さているわけではありません。このとき(44’)と(45’)の両辺をyとxでクロス微分したものを片々引き算して右辺の圧力項を消去した式を得る。それに連続の式(46’)を適用するとが得られる。このとき得られる渦度方程式(48’)では、非発散性のためにf0に関係する項が無くなっている。慣性重力波の場合コリオリ力f0は発散性を通して運動に効いていたのですが、非発散性の場合、f0は慣性振動には効いてきますが、波動運動(渦度方程式)には見かけ上効いてこない。コリオリ力の緯度依存性(β項)のみが効いてきます。連続の式(46’)が成り立つ場合、2.(1)3.で説明した“流線関数”Ψを導入する事ができます。2.(1)3.の議論は定常流についてのもので、波高は時間的に変化しないが場所的には変化しているものでした。ここでは波高は時間的・場所的に変化しないが、流速分布は時間的に変化する非定常な場合を考察します。非定常な場合でも、非発散性(波高が変化しない)が成り立てば流線関数は定義できます。ここでは流線関数と波高の関連性は有りません。いずれにしても、流線関数Ψ(x,y)を用いると流速成分(u,v)はと表される。これらを(48’)式に代入するとが得られる。これは5.(2)2.で求めた流速成分(u,v)が満たす準地衡流方程式の三階時間偏微分項とf0項を省略した形と同じです。今用いている近似の場合、元の流速成分(u,v)もこの形の方程式を満たしていると考えられる。話を簡単にするために、平面波解を考察する。この形の解が方程式(59)を満たすには、角振動数σはでなければならない。これは5.(3)3.の(57)式の分母の(1/Rd2)の項が無い場合に相当します。(61)式と(10)、(60)式を用いると流線関数Ψと流速解(u,v)はとなる。これらの波の位相速度は波長とともに変化するので分散性の波です。これらの波は波数ベクトルとは逆方向に伝播します。位相速度は波数ベクトルの方向に依存して変化することや、位相速度の西向きの成分は波数ベクトルの方向によらず一定になること、また群速度ベクトルの方向は位相速度ベクトルの方向と異なることなどは5.(4)で述べた発散性ロスピー波の場合と同様です。[補足説明1]5.(3)3.の(57)式に於いてならば(61)式となりますが、“非発散性ロスピー波”の波長λが“ロスビーの変形半径Rd”に比較して特に短いわけではありません。そのあたりは5.(5)4.で示す三次元表示したσ曲面でl≠0の領域を見れば解るように波長の長いものもあります。(1/Rd2)の項は∂h/∂tに伴う水平発散の効果(f0が関係)を現していたのであって、ここでは発散が無いのでこの項が無くなったのです。この項が無くなるとσ曲面がどのように変化するかは5.(4)5.のσ曲面と5.(5)4.のσ曲面を比較すれば解ります。[補足説明2](59)式はσに関して一つの解しか持たない。これは言うまでも、なく海面に蓋をしてしまったために慣性重波波が発生できなくなったからです。そのとき、角振動数σと波数ベクトル(k,l)の関係がどのようになるのかは、5.(5)4.で示す三次元表示したσ曲面上で検討してみてください。前項の平面波解をアニメーション表示してみる。ここでも以下の具体的な例を示す。流線関数を黄色(高)青色(低)で表現し、流速ベクトルは赤矢印で示す。柱状流体要素はほぼ同じ位置で往復運動をします。(A)(k,l)=(1,0)(B)(k,l)=(2/3,50.5/3)(C)(k,l)=(1/3,80.5/3)アニメーションの色の濃淡は流線関数の高低を表している。速度ベクトルは流線関数の等高線に沿っており、柱状流体要素は流線関数の高い部分を右に見る方向に動く。流線関数の等高線間隔が密なところほど速度は大きくなる。このときk→0となると平面波の波面はx軸に平行になりますが、平面波の位相速度もゼロになりますから、それは非定常な波動ではなくなります。流速場の変動が波となって伝わるときの復元力はβ効果です。β効果が復元力として働くメカニズムはなかなか理解しにくいところですが、(48’)式や(59)式がそのメカニズムを表しています。波動解(60)はそれらの微分方程式を満足しますから、実際にこのようなメカニズムの波が存在するはずです。その様な波の存在を指摘したのがロスビーです。“非発散性ロスビー波”の場合、角振動数σを表す(61)式を変形するととなりますが、これは次のような解釈ができます。前項で示した三つのアニメーションの様に、波数ベクトルの長さ(k2+l2)0.5をすべて共通にして、角振動数σの値をσA:σB:σC=3:2:1(周期をTA:TB:TC=1/3:1/2:1/1)で変えた場合。(k2+l2)0.5を固定すれば、上記の式は、波数ベクトルの先端が、中心位置(β/2σ+1/2Rd2,0)と半径(β/2σ+1/2Rd2)を三通りに変えた、三つの円の上にあることを示している。変えてアニメーション表示する場合。σを固定すれば、上記の式は、波数ベクトルの先端が、波数空間で(β/2σ,0)を中心とした半径β/2σの、一つの円の上に有ることを示している。この二つのやり方の違いを図示するとの様な関係になります。多くの参考書では二番目のやり方で説明していますので、この稿の図と混同しないでください。ロスピー波が発生するとき、波長を決めるメカニズムが有るのか、周期を決めるメカニズムが有るのかによって両者を使い分ければ良い。[ただし、“発散性ロスビー波”の場合は、1/Rd2の項が有るために、上記の様な簡単な関係にはなりません。]二番目の方法で角振動数を同じになるようにして波数ベクトルの方向を変えると位相速度がどのように変化するか見てみよう。上図赤矢印の様に波数ベクトルを変化させてみるわけです。このあたりの関係は5.(4)5.と同様に、非発散性ロスビー波のσ曲面を書いてみれば直ちに了解できます。図中の青色線分OA、OB、OCの傾きが位相速度を表していますが、角振動数σが同じ場合には波数ベクトル(k,l)の大きさ(k2+l2)0.5に反比例することが見取れます。5.非発散性ロスビー波の位相速度と群速度ここで、角振動数σを一定にして波数ベクトルの大きさと方向が変われば、群速度ベクトルの大きさと方向がどのように変わるかを示す。具体的な数値例として前項で示したβ=1、σ=−0.333・・の場合を取り上げて、波数ベクトルkA、kB、kCで表される三つの波を考える。位相速度がすべて負であることは、波数ベクトルとは逆向きに搬送波の波面が進むことを示している。東西方向に関しては、位相速度は必ず西に向かう。群速度に関しては、波数ベクトルにより、様々な方向を向き、その大きさも変化する。東西方向に関しては西に進む場合も有れば東に進む場合もある。このあたりは発散性ロスビー波の場合と同じです。偏西風中に生じるロスビー波を考えるときには、一般流の効果が無視できない。ここでは非発散性で一般流中に生じるロスビー波を考える。一般流がx軸の正方向を向いており速度Uを持つとする。この場合の運動方程式は、ラグランジュ微分D/Dtの中のu∂/∂xの項を残した形の近似式を用いる。(63)式の両辺をxで偏微分したものから、(62)式の両辺をyで偏微分したものを引算して、さらに連続の式(46’)を適用するとなる渦度方程式が得られる。非発散性の二次元流では、流線関数が導入できる。これを用いれば渦度方程式(64)は単一の関数の偏微分方程式に帰着できる。これは(59)式にUが関係する項が付け加わった形をしている。ここで、見通しを良くするために5.(5)2.と同じ平面波解(60)式を仮定する。これを(65)式に代入するとが得られる。(64)式が、ロスビーの論文[参考文献2.]中の(11)式に相当します。ロスビーはこの式を渦位保存則から導いているのですが、渦位保存則はもともと渦度方程式から得られるものですから、ここで導いたものと同じです。ロスビーはこの式に、幾つかの仮定を適用して論文中(17)式を導いていますが、これは上記(66)式の両辺をkで割って位相速度度にした式に相当します。ロスビーはこれらの式を元にして、帯状流(偏西風)の蛇行との関係を論じています。そのとき用いられる主な考え方は、偏西風の蛇行に伴って帯状流中の流体部分が南北に移動すればβ効果により帯状流中での渦度が変化し、帯状流の蛇行が保持され、克つその蛇行が西進するであろうと言うものです[例えは文献9.p12〜13、文献11.p34等を参照]。ロスビーは、さらに今日ロスビー波の特徴として知られている性質を論じています。確かに、そのようなメカニズムが働いて偏西風の蛇行や西進が説明できそうです。そこで行われている仮定や実際の偏西風の蛇行との関係の説明は、かなり曖昧なところがあり、私はロスビーの論文を読んでもよく理解できません。偏西風にロスビー波は何らかの形で関係しているのは確かですが、こんにち偏西風の形成やその蛇行そのものはもっと違う原因に依ると考えられています。[7.(2)参照]そのため、ここではそのあたりの議論はやりません。ロスビーの考察の詳細については文献2.をお読みください。8.文献(6)地形性ロスビー波1.方程式系ここでは緯度変化に伴うコリオリ力変化(β効果)が無視できるf0平面近似で考える。ここで、Dは水深を表しておりβ’は斜面の勾配ではなくて浅くなる割合を表しているのだから、yの正方向に向かって水深が浅くなるためには、β’y項の前の符号は負でなければならないことに注意。この場合、出発点となる運動方程式は5.(1)1.の(1’)、(2’)式においてf=f0(すなわちβ=0)とおいた(44)、(45)式です。ここでは連続の式(46)は使えない。水深Dが場所的に変化する場合にも正しいを用いなければならない。(45)式の両辺をxで偏微分したものから、(44)式の両辺をyで偏微分したものを引算して、さらに連続の式(67)を適用するとなり、5.(1)1.の渦度方程式(48)と全く同じ形の渦度方程式が得られる。地面(海底)の傾斜に依って擬似的にβ効果を生み出すことができる。(44)、(45)、(48)式は5.(1)2.で説明した方程式系と、全く同じですから、そこでの議論がそのまま成り立ちます。wave)といいます。β効果に関係する南北(北に向かってy軸の正方向を取る)の関係を、海岸線に垂直な方向(深海から海岸・浅海に向かってy軸の正方向を取る)の関係に置き換えないといけない。このとき、ロスビー波の位相速度ベクトルは北半球でも南半球でも西向きでした。南半球でのy軸は南極の方向を向いており、f=f0+βy<0ですから、ロスビー波はy軸の正方向に対して右向きに進むからです。これを地形性ロスビー波での疑似β効果に焼き直すと、南半球では深海から海岸(浅海)に向かうy軸の正方向に対して右向きに位相は進むことになります。北半球では地形性ロスビー波はβ効果が増大する方向(海岸側)を右に見ながら進む。南半球ではβ効果が増大する方向(海岸側)を左に見ながら進む。[補足説明]ここで100km沖合に出ると水深が100m深くなる大陸棚について、β値を見積もってみると程度になる。緯度45度の深さ100m程度の大陸棚を伝播するロスビー波の位相速度はとなり、それなりの値となる。2.地形性ロスビー波[模型実験]水深変化は模型実験でβ効果を出すために良く使われる手法です。実験は、まさしくこの地形による擬似的なβ効果を利用したものです。1962年)は、水位が周期4日と9日で規則的に上下動する事を発見した。その後の詳しい観測から波長3500km程度、位相速度4m/s程度でオーストラリア東海岸を大陸中央上空から見て反時計回りに(つまり海岸を左に見て)伝播していることが解った。wave)と名付けた。この水位の変動は、低気圧の移動に伴って吹く強風によって励起されるもので、波長は数千kmになる。周期も4〜5日程度で慣性振動周期(中緯度で1日程度)よりも長い。陸棚波の特徴をまとめると下記のようになる。波高は10〜20cm、波長は数千km、振動周期は4〜5日程度で、伝播速度は数m/s程度です。位相速度がこの程度になるのは、コリオリ力のf0成分に関係するケルビン波と違って地形のβ効果によるロスビー波だからです。波は傾斜した斜面状の海底(大陸棚)を持つ海に生じ、波は等深線に沿って伝播する。伝播の方向は浅くなる海岸線を北半球では右手に、南半球では左手に見ながら進む。伝播の方向はケルビン波と同じですが、位相速度はケルビン波に比べてかなり遅い。ちなみに、水深100m程度の大陸棚を伝播するケルビン波の位相速度は30m/s程度です。沿っても、同じような水位の変化が報告されている。北半球なので北米大陸中央上空からみて時計回りに(つまり海岸を右に見て)伝播する。気象庁(http://www.data.kishou.go.jp/kaiyou/db/tide/knowledge/tide/ijochoi.html)によると、本州の太平洋沿岸や日本海の沿岸では、台風や低気圧の風によって海洋中に励起された陸棚波が非常にゆっくりとした速度で伝搬し、数日経過してから潮位が上昇して浸水被害等が発生することがあるそうです。太平洋沿岸では東から西に、日本海沿岸では西から東に陸を右に見る方向に伝播し順次、潮位偏差の大きな地域が現れる。1971年秋に本州南岸で大規模な異常潮位が発生したが、この要因の一つが台風によって励起された陸棚波(地形性ロスビー波)だと考えられているそうです。次図は福島県の沖合5km(水深20m)に入れた流速計の記録から観測された陸棚波です。wave)の全体像は、下記の方針に従った考察から松野太郎氏によって明らかにされた。ここでは、その内容をできるだけ解りやすく説明しますが、詳細は1966年に松野太郎氏によって報告された参考文献5.をご覧ください。前章5.(1)2.で注意したように、赤道付近では5章のロスビー波の議論で用いた近似は使えない。4章で説明した赤道β平面近似を使わなければならない。ここではx座標を赤道に沿って、y座標を北極方向を正として経線(南北線)に沿って取るとする。f=f0+βyに於いてf0=0と置けるがβ項が無視できない場合です。4章の“赤道ケルビン波”との違いはv≠0の場合を含むもっと一般的な解を求めることです。コリオリ力の緯度依存性と圧力傾度力の両方を考慮するのだが、f0項は無視できる場合の一般的議論です。方程式系は4.(1)1.で説明したものと同じです。赤道をx軸を取るととなります。ここでは赤道ケルビン波の場合と違ってv≠0を含めたもっと一般的な解を考察しますので、vが関係する項も残されています。ここで、(69)式の両辺をxで偏微分したものから、(68)式の両辺をyで偏微分したものを片々引き算してが得られる。連立方程式を、一つの未知関数のみを含む単一の式に変換する。このときyを含む項があるために、hやuに関する式に変形するのは難しいが、vに関する方程式には帰着できる。このあたりは5.(2)2.で注意したのと同じ事情があるからです。“運動方程式”(1’)、(2’)の両辺をtで偏微分し、“連続の式”(70)の両辺をxおよびyで偏微分するとが得られる。ここで、(68’)と(70’)からhの項を消去して、さらにβyを乗じるととなる。同様に(69’)と(70”)からhの項を消去して、更に時間tで偏微分するととなる。“渦度方程式”(71)の両辺をxで偏微分して、更にgDを乗じるとを得る。(72)(73)(74)式を加えあわせるとv(c,y,t)に関する偏微分方程式が得られる。これは5.(2)2.で導いた(52)式のf02がβ2y2に置き換わった形になっている。赤道地帯に於いては中緯度におけるf0の効果をβyが担っているのですが、この係数が変数yの関数であることが重要です。赤道地帯ではこの項のためにy方向に特異な変化が生じる。3.流速成分vと波高h、流速成分uとの関係ここで、流速成分vと波高h、流速成分uとの間に成り立つ関係式を求めておきます。“運動方程式”(68)の両辺をtで偏微分したものと、“連続の式”(70)の両辺をxで偏微分したものからが得られる。同様に、“運動方程式”(68)の両辺をxで偏微分したものと、“連続の式”(70)の両辺をtで偏微分したものからが得られる。解v(x,y,t)が求まると、これらの式に代入する事で、波高解h(x,y,t)と流速成分解u(x,y,t)を決定する事ができます。前節の(75)式は、第4項にyについて変化する係数が含まれているので、y方向に関する境界条件を満足する解は簡単な三角関数型とはならない。そのため、y方向に関しては未知関数Φ(y)を仮定して、次のような試行解を考察する。これを(75)式に代入するとが得られる。これがΦ(y)を決定するための、変数yに関する2階常微分方程式です。ここでと置くことにすると(79)式はとなる。ここでαは長さの二乗の逆数の次元を持つ。実際そうであることはβ、g、Dの次元を考察すればすぐに確認できます。そこで変数yの変わりに次元の無いなる変数ξに置き換えると、上式はとなる。(80)式は、量子力学の調和振動子の固有値問題(一般にはスツルム・リュービル型の固有値問題と言われる)でお馴染みのものです。ここのβ2y2は、ちょうど調和振動子のポテンシャル曲線の働きをしている。function)と呼ばれる特殊関数で記述できる。境界条件から固有値が決定される所の意味は解りにくいところですが、朝永先生が参考文献14.の§41(�B)および付録�\で説明されています。そのまま引用しておきますので、その解説をお読みください。朝永先生の説明から解るように、方程式(80)の解Φ(ξ)がξ→±∞でΦ→0となる境界条件を満たすたすためには、λ/αは勝手な値であることは許されず、必ず奇数の自然数でなければならない。条件がみたされる場合のみ解の存在が許される。そのとき、方程式(80)の解はnの値に応じてと表される。これはn番目の固有振動に於けるξ方向(y方向)の波動の形を表している。ここのAnは任意の定数です。Hn(ξ)は“Hermiteの多項式”と言われるものでn=0,1,2,3,4・・・・・に対してのようなn次の多項式です。この多項式は一般にで表すことができる。この関係式の証明は前記朝永引用文献の(�\・22)式以降を参照されたし。Φn(ξ)の具体的な形をグラフ表示すると下記のようになる。Anに関してはとしている。この様に置くと、(81’)の関数系Φn’は(−∞,+∞)において2乗可積分な関数空間で完全正規直交系をつくる。完全正規直交系の意味は朝永引用文献の(41・69)、(41・70)式を参照されたし。ξをyに戻した(79)式を(76)式に代入すると最終的な解としてが得られる。ここのkとσは、固有値の関係式を元のk、σで表したときの関係式を満足しなければならない。固有値nが0,1,2,・・・の場合の解の具体的な形を順番に見てゆこう。ここで解の性質を調べる見通しを良くするために角振動数σとx方向の波数kに関してなる置き換えをする。要するに、σ*、K*なる無次元量で置き換える。以下で出て来る(σ*,k*)のグラフは、角振動数σを、を単位として、波数kを、を単位として、測定したものと見なすわけです。このような置き換えをすると(83)式はとなる。1.n=0の場合(83’)、(83”)にn=0を代入した式のグラフは下図の様になる。σに関する三次の代数方程式は以下のように因数分解できる。それは上図右側のグラフから明らかです。それぞれの解を表す式でσig-*=σR*と置けば求まります。この点を境にしてロスビー波と慣性重力波の解釈を下記のように繋ぎ替えるべきです。そのあたりの事情は上右グラフ中に記している通りです。σ*とk*の関係式をσとkの関係式に戻して(82)式に代入すれば、n=0に対するv(x,y,t)の解が求まる。H0(ξ)=1であることを考慮するとが得られる。vの添え字0はn=0に対する解であることを意味する。ここで、以後の議論を簡単にするために、6.(3)でyに対して行った無次元化を、x、t、h、u、vに対しても行う。t*、σ*、k*で表わすととなる。νの添え字0はn=0に対する解であることを表している。これらの解ν0(χ,ξ,t*)を6.(1)3.の(76)、(77)式に代入して、偏微分方程式を解けば波高解η(χ,ξ,t*)とx方向の流速成分解μ(χ,ξ,t*)が求まる。“慣性重力波”と“ロスビー波”の違いは、σ*の所だけで、ν0の関数形そのものは両者で同じです。ηとμを求める手順はすべて同様に取り扱える。ここで、(76)、(77)式を無次元量を用いた偏微分方程式に変換しておく。無次元量に変換するとき、2階偏微分項の上側微分は1乗、下側微分は2乗で変換式を適用する。上記の変換式を用いて整理し、両辺の共通項を割り算消去するとおよびとなる。解ν0(χ,ξ,t*)を(76’)に代入すると同じく、解ν0を(77’)に代入してが得られる。以上をまとめるととなる。境にて性質が大きく変化するが、解の形から以下のことに注意しなけれはならない。解の形から明らかなようにk*=|σ*|となる領域ではη0とμ0の解が発散してしまう。領域の“西向き慣性重力波”は解として成り立たない。その領域のそれらの解は存在しないと考えるべきです。領域の“西向き慣性重力波”(群速度は東向き)に繋ぎ変わると考えなければならない。wave)と呼ばれる。式を見ただけでは解りにくいので、グラフ表示してみる。以下のグラフでξ=0の線が赤道です。χの正方向が東で、負方向が西です。上記の事情により、“東向き慣性重力波”と“ロスビー波”についてはk*=1.0の場合を、“西向き慣性重力波”についてはk*=0.5の場合を表示することにする。この“ロスビー波”と“西向き慣性重力波”をあわせたものが“混合ロスビー-重力波”です。このとき、さらにH=1、A0=0.7511・・・とする。J.Meteor.Soc.Japan,45,p391〜408,1967年、これらはJ.Metior.soc.Japanのサイトからダウンロード可]。2年振動に於けると赤道波と慣性内部重力波の役割”、「天気」2012年09月号http://www.metsoc.jp/tenki/pdf/2012/2012_09_0807.pdf2.n=1の場合(83’)、(83”)にn=1を代入した式のグラフは下図の様になる。この場合も、正確なσを得るには三次方程式を解く必要がある。ここでは、5.(3)1.で採用した方法で求まる近似的な値を用いることにする。σ≫β/kの場合の“慣性重力波”に対する解としてが、またσ≪k(gD)0.5の場合の“ロスビー波”に対する解としてが得られる。これらをまとめるととなる。H1(ξ)=2ξで表した解ν1(χ,ξ,t*)を(76’)に代入すると同じく、解ν1を(77’)に代入してが得られる。以上をまとめるとが得られる。k*=1の場合をグラフ表示してみる。このとき、さらにH=1、A1=0.5311・・・とする。[東へ伝播する慣性重力波(群速度も東向き)](n=1,k*=1.0,σ*=2.115)[西へ伝播するロスビー波(k*=1なので群速度は西向き)](n=1,k*=1.0,σ*=-0.254)これがエルニーニョ現象を説明する遅延振動子理論で存在の可能性が指摘されている赤道ロスビー波です。それは温度が変化する海洋の内部密度境界層に生じるものです。gやDは内部波の理論に基づいたものに置き換える必要があります。そのため伝播速度は1m/s以下となり、約9ヶ月かけて太平洋を東から西へ伝播すると考えられています。そのとき、4.(1)で説明した赤道ケルビン波と相まって、この遅延振動子理論が予想するエルニーニョ周期はだいたい2年程度になり、実際の2〜6年周期を旨く説明できないようです。いう新しい理論がより広い支持を受けているそうです。[参考文献11.第10章参照][西へ伝播する慣性重力波(群速度も西向き)](n=1,k*=1.0,σ*=-1.861)3.n=2の場合(83’)、(83”)にn=2を代入した式のグラフは下図の様になる。この場合も、正確なσを得るには三次方程式を解く必要がある。ここでも、5.(3)1.で採用した方法で求まる近似的な値を用いることにする。σ≫β/kの場合の“慣性重力波”に対する解としてが、またσ≪k(gD)0.5の場合の“ロスビー波”に対する解としてが得られる。これらをまとめるととなる。表した解ν2(χ,ξ,t*)を(76’)に代入すると同じく、解ν2を(77’)に代入してが得られる。以上をまとめるとが得られる。k*=1の場合をグラフ表示してみる。このとき、さらにH=1、A2=0.2655・・・とする。赤道波の全体は下図の様にまとめることができる。図には4章で考察した“赤道ケルビン波“(空色破線)も附記してある。置いた方程式(83’)を解いた解に相当します。実際の現象では、どういった原因で波動が生じるのかというメカニズムの問題やx方向の境界条件に依存するが、これらの波が重ね合わされた形で存在すると考えて良いであろう。8.文献7.編集後記(1)お断りこの稿は、ただ一層からなる非圧縮性流体の浅い層に生じる波の議論です。そのとき、コリオリ力(f0効果)とその緯度依存性(β効果)が波動の性質にどのように関係するかを説明するものです。現実の大気や大洋は、温度(塩分濃度)成層に伴う密度成層をしておりその内部境界面を伝わる波が重要になります。そのとき波の伝播速度は高度(水深)によって変化し、波動は水平方向のみならず鉛直方向にも伝播します。ここでの議論は、そういった事柄をすべて省略して単純化したものです。コリオリ力効果、β効果に起因する複雑な性質を理解するには、上記の簡単化はやむを得ないところかもしれません。一つお断りしなければならないのは、ここでは、波が発生するメカニズムについて何も説明していないことです。この稿の内容は、何らかの原因で適当な波長・振動数の波が存在したら、その波はどのような性質を持っており、どのように振る舞うのかを説明するものです。波の波長や振動数を決めるメカニズムは有るのか、波動を励起する物理的過程は何なのか、等々・・・について何も論じていないのは申し訳ないところです。ロスビー波は大気や大洋の中に確かに存在します。比較して位相速度はかなり遅く周期も長いために、様々なノイズの中に埋もれていて、時間的・場所的に平均化しないとその存在を認知するのは難しいようです。持って西進するロスビー波は確かに存在して、大気や大洋の振る舞いに大きな影響を与えているようです。(1999年刊)第7章“発見と立証(2つの研究例)”この本は気象解析の思想書のような本でとても面白い。(2)偏西風とロスビー波ロスビーの論文(文献2.)は「大気の帯状循環流の強さの変動と作用の準恒久的中心の変位との関係」というとても長い題が付いています。これはもともと移動性高・低気圧を運ぶ偏西風帯の蛇行とその変化の特徴を明らかにするためになされたもので、“ロスビー波”の論文として知られています。私自身この論文を注意深く読んでも偏西風との関係が今ひとつ良く理解できません。ロスビー波と偏西風との関係はロスビーが当初予想したような簡単なものではないようです。1949年))によると考えられています。中緯度に発生する寒冷な高気圧、温暖な低気圧が偏西風を作り出し・強め・加速し・保持し、その蛇行を生み出すと考えられています。5.中緯度における循環を参照]そのとき、偏西風の存在と蛇行は高気圧、低気圧の発達・減衰にフィードバック的に影響を与えます。偏西風が北上すればβ効果により時計回り(北半球で)の回転が強まりますから高気圧性の渦を加速して強めます。南下すれば同じくβ効果により反時計回り(北半球)の回転が強まりますから低気圧性の渦を加速して強めます。偏西風によって西から東へ一日あたり約1000km流されて行きます。このように、偏西風波動と移動性高・低気圧は互いに相互作用しており、単純に切り離して論じる事ができません。そのため偏西風の蛇行をロスビー波で単純に説明することは難しいようです。確かに、偏西風の蛇行は複雑ですが、部分的にロスビー波で旨く説明できる現象もあるようです。山岳や海陸の温度差などによって励起された内部ロスビー波は鉛直上方に伝播して、成層圏・中間圏の大気循環に顕著な効果を及ぼすと考えられているようです。8.文献8.参考文献ケルビン波のページを作るために浅水方程式を復習したので、ロスビー波を説明するページも作ろうと思い立ちました。私にとって気象学の論文は難しくて、ロスビー波も何のことか良く理解できていなかったのですが、今度こそ理解したくて下記の文献などを参考にして試行錯誤しながら作りました。ここでは、ロスビー波の発生原因について何も考察していませんが、私自身この稿を作ることで浅水方程式にコリオリ力が関わる所の全体像がやっと理解できました。全体的に眺めることで、ロスビー波の性質を決めるf0効果とβ効果の働きどころも解ったような気がします。p1〜75,1926年(Oslo)慣性重力波による潮汐現象を報告した歴史的論文です。下記URLよりダウンロード可。1939年これは“ロスビー波の論文”と言われているものです。下記URLよりダウンロード可。http://www.aos.princeton.edu/WWWPUBLIC/gkv/history/Rossby_collab39.pdfまた、下記論文も参照されたし。下記URLよりダウンロード可。赤道波の全体像を論じた有名な論文です。下記URLよりダウンロード可。(1978年刊)数式の取り扱いが正確で、その説明も解りやすい。Oceanography”下記のURLよりダウンロード可。2013年のテキストです。(1983年刊)面白い本ですが、少し難しいです。力学」(2011年著)下記URLよりダウンロード可。2000年(http://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1152-14.pdf)大気波動力学の全体像が良く解りました。感謝!付記URLよりダウンロード可。function)の解りやすい丁寧な説明があります。この本は若い頃かなり気合いをいれて勉強したのですが、今はその内容をすっかり忘れてしまいました。(1999年刊)などを倉庫から引っ張り出して読み直してみました。廣田先生の熱い語りに導かれて、大気中の風・波動が引き起こす大気現象を研究することの面白さにあらためて引き込まれてしまいました。この中で紹介されている論文をネットで探してみると、今日その大半をpdfファイルとして無料で入手できますね。この稿で立ち入れなかったメカニズムを説明する論文も入手可能です。我々のようなアマチュアに取って面白い時代に成りました。