永久凍土融解・活動層深さシミュレーター

永久凍土（パーマフロスト）の夏季融解で生じる活動層の深さを、ステファン方程式で予測するツールです。年平均気温・融解度日・土壌種・氷含有量・植生・積雪・温暖化シナリオを変えると、現在と将来の活動層深さ、地温、永久凍土から放出される炭素量（CO₂eq）までリアルタイムで見えます。

パラメータ設定

年平均気温 MAAT

°C

永久凍土帯の指標。-2 °C 以下で連続永久凍土

融解度日 TDD

°C·day

夏季の正の度日積算。亜寒帯で 500〜2000

土壌タイプ

熱伝導率 k と比熱・密度を自動設定

氷含有量 θ_ice

単位体積あたりの氷の体積分率

地表植生被覆

苔・地衣類・低木の被覆率（夏季の断熱）

冬期積雪

冬の地温を上げる断熱層として作用

ΔT 気温上昇

°C

IPCC AR6: SSP2-4.5 で +2 °C 程度

予測年数

year

何年後の活動層を予測するか

計算結果

—

現状活動層深 (cm)

—

平均地温 (°C)

—

50年後活動層深 (cm)

—

活動層深変化 (cm)

—

炭素排出ポテンシャル (kg C/m²)

—

CO₂eq 排出 (kg/m²)

—

永久凍土断面 — 季節変化アニメーション

夏は地表から活動層が融解し、冬は再凍結します。緑線は植生、白い帯は積雪、青→濃青のグラデーションは地中温度。CH₄ の気泡が融解部から立ち上ります。

活動層深さ vs 年平均気温 MAAT

土壌種別の熱物性比較

理論・主要公式

$$h_{active} = \sqrt{\frac{2\,k_{soil}\cdot \text{TDD}}{L_v}},\quad L_v = \theta_{ice}\,\rho_{ice}\,L_f$$

k_soil：土壌の熱伝導率 [W/m/K]、TDD：融解度日（K·s）、L_v：単位体積融解潜熱 [J/m³]（氷含有量に依存）。Stefan の凍結融解問題の解析解で、活動層深さは融解熱と熱伝導のバランスで決まる。

$$h' = h_{active}\cdot(1-0.3\,V_{cov}),\quad T_g = T_{air}+5\cdot 0.4\cdot \frac{d_{snow}}{300}$$

V_cov：植生被覆率（夏季最大 30% の遮熱）、d_snow：積雪深 [cm]（冬季最大 40% の断熱で地温を引き上げる）。表層プロセスを補正項として組み込む簡易モデル。

$$\Delta C = C_{perm}\cdot \frac{\Delta h}{1\,\text{m}}\cdot \eta_{decomp},\quad \text{CO}_2\text{eq} = \Delta C\cdot \text{GWP}_{\text{CH}_4}$$

C_perm ≈ 50 kg C/m²（永久凍土帯の典型有機炭素）、η_decomp ≈ 10%、GWP_CH₄ ≈ 28（100年）。永久凍土炭素フィードバックの一次近似。

永久凍土の融解と活動層深さ — 気候変動応答

🙋

「永久凍土」って、北極のずっと凍ったままの土ですよね？でも「活動層」って言葉、初めて聞きました。何が「活動」してるんですか？

🎓

いい質問だね。永久凍土（パーマフロスト）の定義は「2年以上連続して 0 °C 以下にとどまる土壌」のこと。シベリア、アラスカ、カナダ北部、チベット高原で北半球の約 24% を占めるんだ。でも、その表面が一年中ガチガチに凍ってるわけじゃない。夏になると地表から数十センチ〜数メートルだけ融け、冬にまた再凍結する。この「夏に融け、冬に凍る」表層を「活動層」と呼ぶんだ。植物が生育するのも、微生物が有機物を分解するのも、全部この活動層でだけ起きるんだよ。

🙋

なるほど！じゃあ、活動層の深さってどうやって決まるんですか？ツールを見ると、年平均気温じゃなくて「融解度日」っていうので計算してますけど…。

🎓

そう、ここが永久凍土工学の核心だ。鍵を握るのが「ステファン方程式」、h = sqrt(2·k·TDD / L_v) という驚くほどシンプルな式なんだ。TDD（Thawing Degree Days、融解度日）は「夏に気温が 0 °C を超えた日数 × その超過温度」を積算した量で、夏の累積暑さを表す。k は土壌の熱伝導率、L_v は氷を融かすのに必要な単位体積あたりの潜熱。要するに「与えられた熱量で、どれだけ深くまで融けるか」を熱伝導と潜熱のバランスで解く問題だ。だから年平均気温そのものより、夏の暖かさの積算と、土の中の氷の量が効くんだよ。

🙋

「土壌タイプ」を泥炭にすると活動層が浅くなって、岩盤だと深くなりますね。これは熱伝導率の差ですか？

🎓

その通り。泥炭は有機質で空隙が多く、k がたった 0.5 W/m/K しかない。一方、緻密な岩盤は k=3.0 と 6 倍も大きい。同じ夏の暑さを与えても、岩盤の方が下まで熱が伝わって深くまで融ける。だからシベリアの泥炭湿地（バイダラフ）は活動層 30〜50 cm しかないのに、アラスカの礫質斜面では 1.5〜2 m に達することもある。さらに植生被覆と積雪も大きく効く。スライダーで植生 100% にすると最大 30% 浅くなり、積雪 300 cm にすると冬の冷気が遮断されて年平均地温が上がる。地表の苔と冬の雪布団がパーマフロストの守護神なんだ。

🙋

右下の「CO₂eq 排出」って何ですか？活動層が深くなると CO₂ が出るんですか？

🎓

これが「永久凍土炭素フィードバック」と呼ばれる、気候科学で最も心配されている tipping point の一つだ。永久凍土には全世界で約 1,600 Gt、土壌炭素全体の半分にあたる有機物が氷漬けのまま眠っている。マンモスの遺骸も含めてね。温暖化で活動層が深くなると、それまで凍ってた有機物が初めて微生物に分解され、CO₂ と CH₄ になって出てくる。CH₄ は CO₂ の 28 倍の温室効果を持つから、放出が増えるとさらに温暖化が進み、活動層がもっと深くなる…という正のフィードバックループだ。ツールでは ΔT を 4 °C に上げてみて。50 年後の活動層が一気に深くなり、CO₂eq 排出も急増するのが見えるはずだ。

🙋

気候だけじゃなく、北極圏の道路や建物にも影響しそうですね。実際どんな被害があるんですか？

🎓

深刻だよ。活動層が深くなる、あるいは永久凍土そのものが融けると、土壌中の氷が水になって体積が減り、地面が陥没する。これを「サーモカルスト（thermokarst）」と呼ぶ。アラスカ縦断パイプラインは熱が伝わらないよう数千本のヒートパイプ付きパイルで持ち上げてるけど、それでもサポート脚の沈下が頻発してる。ロシアでは 2020 年にノリリスクで燃料タンクが永久凍土沈下で破断し、ディーゼル 21,000 トンが河川に流出する大惨事になった。シベリアのヤクーツク、アラスカのフェアバンクス、ヌナブト準州の集落では、家が傾き、道路がうねり、空港滑走路が割れている。IPCC AR6 では永久凍土融解は気候系の「臨界要素（tipping element）」の一つとして明示的にリストされているんだ。

よくある質問

本ツールは Stefan 方程式 h = sqrt(2·k·TDD / L_v) を使用します。h は活動層深さ [m]、k は土壌の熱伝導率 [W/m/K]、TDD は融解度日（Thawing Degree Days）を秒単位に直した値 [K·s]、L_v は単位体積あたりの融解潜熱 [J/m³] です。L_v は氷含有量 θ_ice と氷の密度・融解潜熱から L_v = θ_ice·ρ_ice·L_f で求めます。さらに植生被覆と積雪による断熱補正を掛け、温暖化シナリオでは TDD と平均気温を引き上げて将来値を予測します。

永久凍土には全世界で約 1,600 Gt（土壌炭素全体の約半分）の有機炭素が長期間凍結されたまま蓄積されています。温暖化で活動層が深くなると、それまで凍結していた有機物が微生物分解を受け、CO₂ と CH₄ として大気中に放出されます。CH₄ は CO₂ の約 28 倍（100 年 GWP）の温室効果を持つため、放出量が増えるとさらに温暖化が進み、活動層がさらに深くなる正のフィードバックを形成します。本ツールではこのループの一次近似として、活動層拡大分から排出炭素量と CO₂eq を試算します。

夏季の活動層深さは、地表に届く熱が土壌へどれだけ伝わるかで決まります。苔・地衣類・低木などの植生被覆は日射を遮り、地表面温度を下げて土壌への熱流入を減らすため、活動層を浅く保ちます。本ツールでは植生被覆 100% で最大 30% の融解抑制を仮定します。一方、冬期の積雪は逆方向に働きます。雪は熱伝導率が低いため毛布のように冷気を遮断し、冬の凍結を弱めて年平均地温を上昇させます。本ツールでは積雪 300 cm で最大 40% の冬期断熱効果を与え、平均地温を引き上げる形でモデル化しています。

ステファン方程式は一次元・準定常・潜熱支配を仮定した古典的な解析解で、シベリアやアラスカの観測値（活動層 30〜200 cm）と桁オーダーで一致することが知られています。一方で、地形による日射差、土壌の含水率の鉛直分布、年々の気象変動、地下水流動、地盤の不均一性などはモデルに含まれません。インフラ設計では本ツールのような概算で初期検討を行い、現地モニタリング（thermistor、サーモカルスト測量）と CryoGrid・GIPL2 などの非定常 1〜3次元モデルで詳細評価するのが標準的な流れです。

実世界での応用

北極圏インフラ設計：アラスカ縦断パイプライン（Trans-Alaska Pipeline、800 マイル）は、活動層と永久凍土の境界が温度変動で動かないよう、数千本のヒートパイプ付き垂直支持パイルで持ち上げて建設されました。本ツールのような活動層予測は、パイルの埋設深さ・基礎の浮き上がり対策・道路の断熱層厚さの初期検討に使われます。鉄道（バイカル・アムール鉄道）、空港滑走路、ガス田開発（ヤマル半島）でも同様の検討が不可欠です。

気候モデルと炭素循環研究：CMIP6 / IPCC AR6 の地球システムモデル（ESM）では、永久凍土炭素フィードバックがまだ十分に組み込まれていません。CryoGrid、GIPL2、CLM-Permafrost などの 1〜3 次元モデルで活動層の動態と炭素分解を計算し、衛星観測（GRACE による地下水量変化、Sentinel-1 SAR による地表変位）と組み合わせて、将来の CO₂・CH₄ 放出量を推定します。本ツールはその基礎となるステファン方程式の感度を直感的に理解するための入口です。

サーモカルスト・地盤沈下リスク評価：シベリア・ヤマル半島では 2014 年以降、メタン爆発と思われる巨大クレーター（直径 30 m 級）が次々発見されています。2020 年のノリリスク燃料タンク事故（21,000 トン軽油流出）も永久凍土沈下が原因でした。集落・パイプライン・原油タンクの周囲では、本ツールで将来の活動層深さを推定し、サーモカルストや地盤の不均一沈下を予測した上で、ヒートパイプ・断熱層・基礎強化などの対策を選定します。

先住民コミュニティと食料安全保障：カナダ・ヌナブト準州やシベリアの先住民集落では、永久凍土に掘った地下貯蔵庫（ice cellar）で長年トナカイ肉や魚を保存してきました。活動層の拡大で貯蔵庫内温度が上昇し、伝統食の保存が困難になっています。また、トナカイの越冬地形変化、湖沼の急速排水、海岸侵食（北極沿岸で年 5〜20 m）など、生活基盤そのものを脅かす変化が起きています。本ツールのような概算でも、地域別の長期見通しに貢献できます。

よくある誤解と注意点

まず最も多い誤解が、「永久凍土＝表面まで一年中凍結」というイメージです。実際は地表から数十センチ〜数メートルは毎年融け（活動層）、その下に永久凍土層が広がります。連続永久凍土帯（北極海沿岸）では永久凍土層の厚さが数百メートル〜1500 m に達することもありますが、不連続永久凍土帯（南端）ではパッチ状で 10〜50 m と薄く、温暖化で先に消失します。本ツールが計算するのはあくまで活動層であり、永久凍土層全体の厚さや消失時期を予測するには、より長期間（数百〜数千年）の非定常熱伝導解析が必要です。

次に、「Stefan 方程式は精密な予測ツール」という誤解。本式は一次元・準定常・潜熱支配の理想化であり、現実の永久凍土には地下水流動、不均一な含氷量、季節の気温変動、地形による日射差、植生フィードバック（湿地化）など多くの過程が絡みます。シベリアの観測サイトでは、観測値とステファン解の差が ±30% に達することも珍しくありません。本ツールの値は「概算」と「感度の理解」のためであり、実プロジェクトでは現地ボーリング・サーミスタ列・衛星 InSAR モニタリングと CryoGrid 級モデルを併用してください。

最後に、「永久凍土が融けても問題は地元だけ」という誤解。実際には永久凍土炭素フィードバックは全球規模の問題です。Schuur ら（Nature 2015）の推計では、4 °C 温暖化シナリオで 2100 年までに 130〜160 Gt CO₂eq が放出される可能性があり、これは現在の年間化石燃料排出（約 40 Gt）の数年分に相当します。CH₄ は短寿命だが GWP が高く、突発的なクレーター爆発や thermokarst lake からの放出は局所的にも検出可能なほど大量です。さらに永久凍土には炭素だけでなく、古いウイルス・細菌（2014 年に Pithovirus が再活性化）、水銀（北極水銀ストックの約半分）も封じ込められており、これらの放出リスクも近年急速に注目されています。

永久凍土融解・活動層深さシミュレーター

永久凍土の融解と活動層深さ — 気候変動応答

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点

永久凍土 融解・活動層深さシミュレーター

永久凍土の融解と活動層深さ — 気候変動応答

よくある質問

実世界での応用

よくある誤解と注意点

使い方ガイド

具体的な計算例

実務での注意点

永久凍土融解・活動層深さシミュレーター