「津波の心配はありません」のあと
2022年1月15日19時1分、気象庁は「津波の心配はありません」と伝えた。若干の海面変動はありうるという見立てである。13時ごろに南太平洋のフンガ・トンガ-フンガ・ハアパイ火山が大規模に噴火してから6時間がたっていたが、日本へ向かう途中の潮位観測点に大きな変化は出ていなかった。

異変は20時台に始まる。
父島をはじめ各地で、地震の津波なら届くはずの時刻より2時間以上も早く潮位が動き出した。同じころ、国内の気圧計は2 hPaほどの変化を記録している。潮位の変化はその後も大きくなり、翌16日0時15分、気象庁は奄美群島・トカラ列島に津波警報を、太平洋沿岸などに津波注意報を出した。津波の高さの測り方で読むと、奄美市小湊で134 cm、岩手県の久慈港で107 cmである。
まだ来ないはずの波が、なぜ来たのか。
答えは海の上の空気にあった。噴火で生まれた「ラム波」という気圧の波が秒速約300 mで太平洋を渡り、海面を持ち上げながら進んだ。ふつうの津波は、太平洋の水深では秒速200〜220 mでしか進めない。だから、気圧の波に連れられた水位の変化が先に届く。

持ち上げる向きも、直観に反する。高潮で習うとおり、気圧が1 hPa下がれば海面は約1 cm上がり、気圧が上がれば海面は下がるはずだ。ところが久保田達矢らの計算では、この夜の気圧の山の下で、沖の海面は上がる。
クラカタウの気圧計から、26日後の原稿まで
遠い噴火が気圧計に届いたのは、これが初めてではない。1883年8月のクラカタウ噴火のあと、英国のR. H. スコットは、8月27日から31日にかけてヨーロッパを通った一連の気圧の乱れを王立協会に報告した。気象庁の検討会も、気圧波に伴う津波の例にクラカタウを挙げ、こうした潮位変化の頻度を「約150年で2~3回程度」と書く。

ラム波という名は、英国の応用数学者ホレス・ラムにちなむとみられる。ラムは1911年に論文「大気の振動について」を書き、教科書『流体力学』の版を重ねた。久保田らの原稿も、ラム波を説明する箇所でこの教科書を引く。

海の側の式は1929年に出た。リバプール潮汐研究所を率いたジョゼフ・プラウドマンの論文「大気圧の変化が海に及ぼす影響」である。動く気圧が海の長波に近い速さで進むと、気圧の傾きが同じ波に仕事をし続け、波が育つ。今はプラウドマン共鳴と呼ばれ、長崎の「あびき」のような気象津波を説明するときに使う。

1967年には地球物理学者のD. ハークライダーとF. プレスが、クラカタウの波を大気と海が結合した系のパルスとして論じた。ただ、クラカタウは現代の観測網ができるはるか前の出来事である。
2022年の噴火は違った。世界の海底に水圧計が並んでいた。防災科学技術研究所の久保田らは2021年にも、東北沖の海底水圧計で気象津波を調べる論文を出している。噴火の26日後の2月10日には、同じ研究所の齊藤竜彦、東京大学地震研究所の西田究と書いた原稿がEarthArXivに載った。査読を経た論文は5月12日(米国東部時間)に『サイエンス』のサイトに載り、7月の号に収められた。
久保田らは秒速300 mの気圧波で海の波を計算し、世界の海底水圧計に残った第一波を再現した。防災科研の発表は、この波について「伝播と共に増大していないことから典型的な共振現象とはいえず、むしろラム波による海面の強制振動によって引き起こされたと考えられる」と書く。一方、気象庁の津波予測技術に関する勉強会は同じ年の4月、共鳴には「10,000m近い水深が必要」としながら、それより浅くても「大きく増幅されうる」と書き、効果の大きさは今後の調査に委ねた。
押したのか、共鳴したのか。式を立てれば、この2つを同じ1行で読める。
気圧の山が海面を持ち上げる理由
津波の波長は水深よりずっと長く、沖で水位が変わる幅は数cmで、水深の数kmより桁違いに小さい。そこでナビエ・ストークス方程式から粘性と非線形の項を落とし、水深方向に平均した線形の浅水方程式を使う。海面の高さを、水平の流速を、一定の水深を、海水の密度を、気圧の変化をとする。
-
水の体積の保存を書く。
-
水平の運動方程式を書く。力は海面の傾きと気圧の傾きの2つである。
-
1をで、2をで微分し、を消す。左辺は速さの波動方程式で、右辺が気圧の駆動である。
-
気圧の山が形を変えず速さで進むとし、、を代入する。
-
について解く。はフルード数である。
同じ気圧の山でも、ゆっくり動くと山の下で海面が下がり、津波より速く動くと山の下で海面が上がる。 -
気圧がゆっくり動く極限をとる。高潮の「1 hPaで約1 cm」がここで出る。
-
気圧が津波より速いときは符号が変わる。
-
噴火の時刻に海が静止していたとする。解は強制波と、左右へ進む2つの自由波の和で、と書く。
-
とをで課す。
-
とを解く。ならで、強制波を沈降の波が追う。
-
の極限をとる。発散は消え、距離に比例して育つ項が残る。
なぜ符号が変わるのか。海の1点に立って考えると見える。山の前面が近づくと、気圧の傾きが水を進行方向へ押し、流速が増えていく。頂が真上に来たとき流速は最大になり、後ろ側で減って止まる。つまり山の前面の下では、後ろの速い水が前の遅い水に追いつき、水が詰まって海面が上がる。頂が過ぎると水は離れ、海面は元に戻る。海面の山が気圧の山の真下に来るのは、このためだ。
山がゆっくりなら話は逆になる。水には山の下から逃げる時間があり、海面の傾きが気圧の傾きとつり合うまで海面が下がる。速い山は、水が逃げる前に通り過ぎる。
手順9のは水の体積の話である。持ち上げた分の水は、どこかから来なければならない。その穴が発表にある「海面沈降の波」で、重力で崩れて秒速200〜220 mで後を追う。
手順11のは気圧の傾きそのものだ。防災科研の発表が用語の説明で書く「気圧の水平勾配による仕事」が、距離に比例する項として見える。
を0から動かすと、倍率はから下がり続け、を境に正へ飛び移る。、なら1.20倍。にするとの分母はほぼ0になる。
もっとも、ここまでは平らな海底の1次元の話である。実際の波は火山から円く広がり、海底の起伏で散乱もする。
2 hPaから海底の4 hPaへ
沖の海面は、どれだけ上がったのか。久保田らの査読前の原稿によると、世界の気圧計で見た気圧の山は約2 hPaである。秒速300 mで、重力加速度を9.8 m/s²、海水の密度を1025 kg/m³とする。
-
2 hPa = 200 Paを水の高さに直す
-
水深5,000 mの津波の速さを出す
-
フルード数の2乗を出す
-
倍率を出す
-
海面の盛り上がりを出す
-
海底の水圧は気圧と水の重さの和である
-
水深5,000 mの値を入れる
-
水深4,000 mの値を入れる
上空の2 hPaが、海底では3.5〜4.4 hPaになる。原稿によると、海底水圧計の観測は約3〜4 hPa、地形を入れた計算は約4 hPaだった。原稿の補足資料も同じ関係から、水深約4,000 mで1.8倍と見積もっている。海底地形を入れない1行の式でここまで近づくのだから、少し驚く。
では、共鳴はどの深さで起きるのか。
-
をについて解く
気圧の速さごとに、津波の速さと等しくなる水深を棒で示す。ラム波の速さでは9 km以上、大気重力波の速さではトンガと日本の間の典型的な水深に近い。
ラム波の速さは、測り方で幅がある。勉強会の報告書は到達時刻に合う値として312 m/s、音速として340 m/sを挙げ、Wrightらの解析は地表近くで318.2±6 m/sとした。同じ式で312 m/sなら9.9 km、340 m/sなら11.8 kmで、後者は報告書の11,796 mと一致する。勉強会の「10,000m近い」も、この範囲に入る。しかし勉強会が典型的な値に置いた水深5,000 mでは、倍率は1.2倍にとどまる。発散からは遠く、山と自由波も早く離れる(例題5)ので、波は距離とともに育たない。久保田らが共鳴ではないとした理由と、この数字は矛盾しない。
ラム波の速さが音速に近い理由も、短い計算で出る。鉛直の動きが無い音波を考える。
-
水平の運動と連続の式を書く
-
断熱の関係を入れる
-
とを消す
-
理想気体の音速を出す。、とする。
-
鉛直は静水圧のつり合いで決まる(は一定温度の大気のスケールハイト)
気圧の変化は上へ行くほど弱まり、エネルギーは地表の近くに集まる。だから波は上へ逃げず、水平に広がって遠くまで届く。勉強会の報告書にある「大気と地面や海面との境界に捕捉されて」という性質である。実際の大気は上空ほど冷たく、風の影響も受けるので、地表の音速340 m/sより遅くなる。
秒速300 mで到達を予想する
時刻の差を手で出す。 火山から東京までは、大円で約7,840 kmある。噴火を勉強会の報告書の計算と同じ13時10分とすれば、300 m/sの波は20時26分に届く。気象衛星ひまわりの画像から、気圧波は20時40分ごろに本州付近へ達していた。200〜220 m/sの津波なら23時4分から0時3分で、差は2.6〜3.6時間になる。
資料の数字は少しずつ違う。見ている観測点と比べ方が違うからだ。『サイエンス』の論文の要旨は2時間超、防災科研の発表は世界の水圧計について「2~3時間以上」、気象庁の検討会は日本の潮位観測点について「2時間以上」と書く。水深を入れた伝播計算と潮位を地域ごとに比べた勉強会では、東北から東海で3〜3.5時間、四国から沖縄で3.5〜4時間だった。
地球を一回りした時間から速さを出す。 海底の水圧計S-netは、1月17日9時ごろ、気圧波が地球を一回りしてきたのに対応する変化もとらえた。20時40分からの約36時間20分で4万km進んだなら、およそ306 m/sになる。
到達の情報に使う。 気象庁は2022年4月7日から、火山噴火のときはラム波の速さで求めた到達予想時刻を「遠地地震に関する情報」に載せることにした。同じ年の検討会は、どの気圧波が大きな潮位変化を起こすかは現段階で不明だとして、最も速い波の時刻を知らせるのが適当とした。情報の中では「津波」と呼び、津波警報・注意報の仕組みを使うのも適切としている。
水圧計の読みを直す。 上の計算の手順7にある4.4 hPaをそのまま水の高さと読むと約4.4 cmで、海面の2.4 cmの2倍近くになる。地震の津波では気圧の分を無視できたが、今回は引かないと波高を見誤る。
沖の計算は2.4 cmだが、岸の観測は134 cmに達する。この開きは、1行の式では埋まらない。大陸棚で浅くなると波は高くなり、港にはそれぞれ揺れやすい周期がある。勉強会が全国の観測点で比べると、最大の高さのときの周期は港の固有周期とおおむね一致した。それでも、最大の高さを生んだ主な要因とその割合は、勉強会の時点で見積もれていない。
この潮位変化で、日本では避難中の転倒による軽傷者が2人出て、高知、徳島、三重、宮城の4県で転覆や沈没などの船が30隻にのぼった。国外では、気象庁の月報によるとトンガで3人、ペルーで2人が亡くなっている。津波警報が出たときの行動は気象庁の案内にある。
水深と速さから手で出す
例題 1基本速さがそろう水深
大気重力波の速さを215 m/sとすると、津波と速さがそろう水深はいくらか。300 m/sのラム波と比べよ。
解説を見る
に入れる。
215 m/sは太平洋の典型的な水深に近く、300 m/sはその約2倍の深さを要する。
例題 2入試水深6,000 mの強制波
水深6,000 mの海の上を、高さ1 hPaの気圧の山が300 m/sで進む。海面の盛り上がりと海底の水圧の変化を求めよ。 kg/m³とする。
解説を見る
水の高さに直す。
海面は cm上がる。海底の水圧は倍で、約2.9 hPaである。
例題 3大学沈降の波の大きさ
水深5,000 m、2 hPa、300 m/sのとき、証明の手順10の, , をcmで求め、和が0になることを確かめよ。
解説を見る
m/s、 cmなので cmである。
となる。沈降の波は隆起より少し大きく、後ろへ進む小さな隆起と合わせて体積が保たれる。
例題 4大学岸へ近づくと何倍か
エネルギーの流れが保たれるとして、振幅と水深の関係を導け。水深5,000 mで2.4 cmの波が水深100 mへ進むと何cmになるか。
解説を見る
から、次の関係(グリーンの法則)が出る。
なので、約6.4 cmである。Villalongaらの論文も、岸へ近づく増幅を2〜4倍としている。134 cmには、これだけでは届かない。
例題 5大学計算が壊れない刻みと、山と波が離れる距離
(1)格子の幅 km、水深12 kmで線形の浅水方程式を解く。CFL条件を満たすの上限を求めよ。
(2)幅 kmの山が300 m/sで進むとき、山が自由波からだけ離れるまでに進む距離を、水深5 kmと8 kmで求めよ。
解説を見る
(1) m/sなので、次のとおりである。
この記事のシミュレーションは10 sを使い、比は0.34になる。
(2)離れるまでの時間はである。山が進む距離はこれにを掛けてになる。
浅いと共鳴は早く切れ、深いと長く続く。勉強会の「長く維持されれば」は、この距離の話である。
気圧の競走と2段の波
水深を9 kmへ近づける
線形の浅水方程式を、気圧の山を動かしながら実際に解いている。水深5 kmでは、盛り上がりは山に付いて走るが高さは一定で、その後ろに沈降の波が置いていかれる。浅くするほど盛り上がりは小さくなる。水深を9.2 kmへ近づけると、山の前が上がり後ろが下がる形で、距離とともに波が育つ。9.2 kmより深くすると自由波が山を追い越し、山の下の海面は下がる側へ回る。読み出しの「定常解」は山と自由波が離れたあとに落ち着く値で、9 km付近ではこの距離のうちに届かない。
「式の解を重ねる」を入れると、証明の手順10と11の解が、破線で少し上にずれて重なる。太い線の数値解と形が合うかを見てほしい。
同じ分母が現れる場所
- 地球科学 潮位の変動が高潮か気象津波かを見分けるとき、気圧の速さと津波の速さの比を見る。なら静水圧のつり合いで1 hPaあたり約1 cm、が1に近いと共鳴で気圧の変化以上に潮位が動く。
- 物理学 動く外力で波を駆動する問題を解くときに、共振の分母が同じ形で現れる。波数の成分で、とすれば、は強制振動のと同じ形である。
- 工学 港や湾の長周期の揺れを調べるとき、沖から入る水位変動の周期と港の固有周期の近さを比べる。沖の入力はの側で決まり、港の中の増幅は固有周期との比で決まるので、2つを分けて考えられる。
押したのか、育てたのか
- 2022年Omiraら(ポルトガル海洋大気研究所、リスボン大学など): 観測と数値モデルから、音響重力波という動く波源が海にエネルギーを渡すと考え、早い到達と長い継続を説明した。トンガ・ケルマデック海溝の深い海では、共鳴で波が増幅するという。が1に近づく場所を重く見る立場である。(出典: Global Tonga tsunami explained by a fast-moving atmospheric source, Nature)
- 2022年SekizawaとKohyama(東京大学、お茶の水女子大学): 1次元の浅水モデルで水深の分布を変えて実験し、日本で観測された時刻と振幅に合う結果を得た。深い海盆での共鳴に近い増幅と大陸斜面での浅水変形が効き、海溝のように深く狭い地形は二次的だという。同じでも、海溝の役割は論文によって見方が分かれる。(出典: SOLA掲載論文の原稿)
- 2022年Nishikawaら(高知工科大学、九州大学など): 全国に置いた微気圧計と海底水圧計の記録から、ラム波の310 m/sのあとに215 m/sの大気重力波を見つけ、この波が共鳴で津波を励起したとみる。太平洋の典型的な水深ではになる速さである。(出典: Scientific Reports)
- 2023年Villalongaら(スペイン海洋研究所バレアレス海洋センター、バレアレス諸島大学など): 世界の検潮所の記録を比べ、大陸棚と港の固有モードが、沖で生まれた波の特定の周期だけを増幅すると論じた。式のは沖の入力で、岸の高さはその先で決まる。(出典: Scientific Reports)
参考文献
- T. Kubota, T. Saito, K. Nishida, "Global fast-traveling tsunamis driven by atmospheric Lamb waves on the 2022 Tonga eruption", Science 377 (2022), 91-94. doi:10.1126/science.abo4364(要旨はEurope PMC、数値の一部は査読前の原稿EarthArXivで確認)
- 防災科学技術研究所・東京大学地震研究所「2022年1月トンガ噴火に伴う地球規模の津波発生と伝播メカニズムを解明」(2022年5月13日)
- 津波予測技術に関する勉強会「フンガ・トンガ-フンガ・ハアパイ火山の噴火により発生した潮位変化に関する報告書」(2022年4月)、火山噴火等による潮位変化に関する情報のあり方検討会「火山噴火等による潮位変化に関する情報のあり方(報告書)」(2022年7月)
- 気象庁「地震・火山月報(防災編)令和4年1月 特集」(2022年)
- J. Proudman, "The effects on the sea of changes in atmospheric pressure", Geophys. Suppl. Mon. Not. R. Astron. Soc. 2 (1929), 197-209.
- H. Lamb, "On atmospheric oscillations", Proc. R. Soc. Lond. A 84 (1911), 551-572. 同Hydrodynamics, 6th ed. (1932).
- R. H. Scott, "Note on a series of barometrical disturbances which passed over Europe between the 27th and the 31st of August, 1883", Proc. R. Soc. Lond. 36 (1883), 139-143.
- D. Harkrider, F. Press, "The Krakatoa air-sea waves: an example of pulse propagation in coupled systems", Geophys. J. Int. 13 (1967), 149-159.