２０２３年12月4日
第１４回歴遊会かねさは

１８５６年安政江戸台風（安政3年の大風災）による高潮、高波の試計算

福谷 陽
関東学院大学理工学部 准教授

安政江戸台風の概要■台風通過日1856年9月23日から9月24日（安政3年8月25日から26日）に関東周辺を通過し、江戸湾で大規模な高潮被害が発生■経路（平野(2017)による推測）
静岡県伊豆半島付近に上陸し、江戸のすぐ西を通過。
その後、関東地方北部、東北地方へ進んだと推定（ただし、坂崎ら(2015)は伊豆半島東側を北上、と推測）【出典】平野順平(2017)歴史イベントと気候との関わりをどう教えるか─歴史気候学からの視点，歴史と地理708：
地理の研究197号

安政江戸台風の概要◼ 日記に書かれた風向変化や天気（晴や大雨など）の記述が、台風の経路を推定するのに役立つ

参考文献【出典】平野(2017)

高潮被害の概要◼ 一説には、死者約10万人と推定◼ 伊豆半島西側、三浦半島東側～東京湾の西側（芝・品川・本所・深川など）にかけての広い範囲で床上三尺（約1m）の水深が推定（坂崎ら(2015)の記述）◼ 東京湾では最大水位2.5m～3.2mに達する高潮が発生し、水路や川が氾濫（鎌ヶ谷市HPの記述）◼ 高潮被害により、江戸や浦安、市川、船橋で浸水が発生し海のようだった◼ 木造建築が多く、風に弱いため倒壊した建物から火災も発生『安政風聞集』の被害記録マッピングの図。
赤点が高潮の記述があった地点。
（坂崎ら, 2015）

高潮被害の概要

現在の江東区大島にある本所羅漢寺付近の高潮の様子（安政風聞集より、作者不明）※この付近には荒川があり、高潮に加えて、洪水氾濫により家が流され、町全体が海のようになったのではないかと推測される。
（平山先生より情報提供）

江戸の様子（『安政風聞集』より）
https://edo-g.com/blog/2019/10/typhoon.html

高潮被害の概要

江戸の様子（『安政風聞集』より）
https://edo-g.com/blog/2019/10/typhoon.html

安政風聞集（あんせいふうぶんしゅう）◼ 金屯道人 （きんとんどうじん ）こと仮名垣魯文による風水害ルポ 。
安政4年（1857）刊◼ 倒壊した築地西本願寺の本堂や半壊した永代橋など 、すさまじい風
雨の爪あとを多色刷りの図と文章で紹介すると共に 、洪水のなか幼児を背負い大木にしがみついて夜を明かし九死に一生を得た男の話
ほか災害が生んだ奇談悲話も載せています。
◼ 全3冊

天候記録の記述から推定した台風の通過時刻◼ 台風が関東地方を通過した時間を1856年9月23日19時～24日4時と推定◼ 横浜や立川では、22時頃に次第に風が強くなっていることが分かり、東風や南風が吹いていると推定される場所にプロット◼ 玉村では、4時頃に北風が吹いていることが分かっているため、玉村の東側にプロット

4時
3時

2時
1時
0時
23時

当時の日記に記録された天候（平野ら, 2017)）
22時

地名

9月23日

9月24日
21時

静岡県下田市
神奈川県横浜市

22時頃（亥の刻）前に激しい暴風雨

晴れ

東京都立川市22時頃（亥の刻）に大雨が降り出す

昨夜に東風が吹き出し、徐々に和らぐ群馬県玉村町

昨夜に北東風が吹き出し、4時頃（七つ）北風に変わる。

20時
19時

設定した台風情報◼ 台風の中心気圧はデータが得られないため、まずは、２０１９年台風１５号（上陸時955Pa程度）の気圧を当てはめた◼ 以下の台風情報を設定して高潮計算に使用
日付

時間

緯度

経度

気圧の低下量（ｈPa）

最大風速半径（km）

移動速度（km/h）

1856年
9月23日

19時

34.525138.583

９６０

21.30

30.4

20時

34.796

138.630

９６０

21.30

30.4

21時

35.117138.793

９６０

21.30

38.6

22時

35.377138.916

９６０

21.30

30.9

23時

35.594139.084

９６０

21.30

28.5

0時

35.715139.135

９６０

21.30

14.2

1時

35.872

139.203

９６０

21.30

18.5

2時

36.025139.319

９６０

21.30

19.9

3時

36.174139.487

９7０

19.17

22.4

4時

36.288139.593

９7０

19.17

15.8

1856年
9月24日

高潮推算の基礎方程式(SuWAT)
連続の式𝜕𝜂 𝜕𝑀 𝜕𝑁++=0
𝜕𝑡 𝜕𝑥 𝜕𝑦

X軸方向の運動方程式

𝜕𝑀 𝜕 𝑀2
𝜕 𝑀𝑁
𝜕𝜂 𝐷 𝜕𝑝𝑜 𝜏𝑠𝑥 − 𝜏𝑏𝑥𝜕2𝑀 𝜕2𝑁
1 𝜕𝑆𝑥𝑥 𝜕𝑆𝑦𝑦++= 𝑓𝑁 − 𝑔𝐷−++ Aℎ+−+𝜕𝑡 𝜕𝑥 𝐷𝜕𝑦 𝐷𝜕𝑥 𝜌𝑤 𝜕𝑥𝜌𝑤
𝜕𝑥 2 𝜕𝑦 2
𝜌𝑤 𝜕𝑥
𝜕𝑦

𝜕 𝑁2
𝜕𝜂 𝐷 𝜕𝑝𝑜 𝜏𝑠𝑦 − 𝜏𝑏𝑦𝜕2𝑀 𝜕2𝑁
1 𝜕𝑆𝑥𝑥 𝜕𝑆𝑦𝑦
Y軸方向の 𝜕𝑀 𝜕 𝑀𝑁++= 𝑓𝑀 − 𝑔𝐷−++ Aℎ+−+
運動方程式 𝜕𝑡 𝜕𝑥 𝐷𝜕𝑦 𝐷
𝜕𝑦 𝜌𝑤 𝜕𝑦𝜌𝑤
𝜕𝑥 2 𝜕𝑦 2𝜌𝑤 𝜕𝑥
𝜕𝑦
速度変化項（非定常項）

移流項

コリオリ項 重力項 気圧傾度項

せん断応力項

拡散項

ラディエーション応力項（Radiation Stress）※波の存在により流体内に生じる応力

𝑝𝑜 ：海面での大気圧
𝜂：静水面からの海面上昇量
𝑀： x軸方向の流量フラックス（海面から海底までの流速の積分量） 𝜌𝑤 ：海水密度
𝑁： y軸方向の流量フラックス（海面から海底までの流速の積分量） 𝜏𝑠𝑥 , 𝜏𝑠𝑦 ：海面せん断応力（風が海面に及ぼす接線応力）
𝜏𝑏𝑥, 𝜏𝑏𝑦 ：海底面せん断応力（海底からの摩擦力）
𝐷：海面から海底までの水深
Aℎ ：水平渦動粘性係数
𝑓：コリオリ係数
𝑆𝑥𝑥 , 𝑆𝑦𝑦 ：ラディエーション応力テンソル項
10

高波推算の基礎方程式(SWAN)ത有義波高𝐻1/3 ・卓越周期𝐸の算出

波作用量平衡方程式

𝜕𝐸 𝜕𝐶𝑥 𝐸 𝜕𝐶𝑦 𝐸 𝜕𝐶𝜃 𝐸+++= S𝑛𝑒𝑡
𝜕𝑡
𝜕𝑥
𝜕𝑦
𝜕𝜃
移流項

屈折項

エネルギー
外力関数

𝐻1/3 = 4.0 ඵ 𝐸 𝑓, 𝜃 𝑑𝑓𝑑𝜃
全方位、全周波数の波のエネルギーを積分し、その平方根に4.0を乗じる

S𝑛𝑒𝑡 = S𝑖𝑛 + S𝑛𝑙 + S𝑑𝑠 + S𝑏𝑡𝑚
𝐸(𝑓, 𝜃)： 波浪2次元スペクトル𝐶𝑥 𝑓 , 𝐶𝑦 𝑓 : 各空間での群速度（波エネルギーの進む速さ）
2𝜋
S𝑛𝑒𝑡 ： 波浪スペクトルに対するエネルギーの入出力を表現する 𝐸ത 𝑓 = න 𝐸 𝑓, 𝜃 𝑑𝜃0外力関数
周波数𝑓ごとに、全方位の波のエネルギーを
S𝑖𝑛 ：エネルギー入力（風による風浪の発生・発達）
積分
S𝑛𝑙 ：非線形エネルギー輸送（波と波の相互作用、3波・4波共鳴）
S𝑑𝑠 ：エネルギー散逸（逆風や砕波による波浪の減衰）
S𝑏𝑡𝑚 ：海底摩擦11高潮・高波の計算結果
台風中心◼ 東京湾奥で最大約1m超、金沢沖では最大約0.5mの潮位偏差◼ 金沢沖の潮位偏差は、２０１９年台風１５号の潮位偏差よりは低くなった

1856年9月23日22時

1856年9月23日23時

1856年9月24日0時

1856年9月24日1時
潮位偏差(m)３０ｍメッシュ標高

水深（ｍ）

１０ｍメッシュ標高

水深（ｍ）

久保田ら(2022)
オランダ軍艦の航海日誌に基づいて推定した安政江戸台風の大きさ◼ 安政江戸台風の1000hPaの半径は359kmと推定
14

参考文献
1. 久保田ら(2022)オランダ軍艦の航海日誌に基づいて推定した安政江戸台風の大きさ、2022年度日本地理学会春季学術大会、セッションID: 516
2. 平野淳平（２０１７）歴史イベントと気候との関わりをどう教えるか─歴史気候学からの視点，歴史と地理708：地理の研究197号
3. 坂崎貴俊ら（２０１５）安政 江戸台風 (1856)の被害と当時の気象場推定，生存圏研究，
11:64-70.

15