How Uncertain Are Estimates of Marine Carbon Dioxide Removal? Insights from a Simplified Model

技術要約：海洋二酸化炭素除去（mCDR）の推定値はどの程度不確実か？簡略化されたモデルによる知見

問題提起
海洋二酸化炭素除去（mCDR）技術、特に海洋アルカリ度向上（OAE）および直接海洋除去（DOR）は、大気中のCO2を減少させるために、表層海水の二酸化炭素分圧を低下させることに依存している。理論上の最大効率（例：OAEにおける1モルのアルカリ度あたり約0.8モルのCO2）は既知であるが、炭素除去の経時的な変化は、物理的および生物地球化学的なプロセスに強く依存する。現在のmCDRの有効性の推定は、海洋モデルに大きく依存しているが、それらのモデルに内在する不確実性（特に垂直混合、海気ガス交換、および水塊の特性に関するもの）は十分に理解されていない。既存の全球モデルは、自然の変動性を過小評価する可能性のある固定されたパラメトリゼーションを採用していることが多く、その複雑さゆえに、CDR計算を駆動する特定のメカニズムを分離することが困難である。この定量化された不確実性の欠如は、MRV（モニタリング・報告・検証）の枠組み、およびmCDRプロジェクトにおける純炭素中立性のタイミングにとって、根本的な課題となっている。

手法
これらの不確実性を体系的に調査するため、本研究では一次元（1D）の非生物的な海洋モデルを用いている。この簡略化されたアプローチは、水平方向のアドベクション（移流）、サブダクション（沈み込み）、および生物学的相互作用を無視することで、垂直プロセス（乱流拡散と海気ガス交換）を分離して抽出するものである。モデルは全アルカリ度（TA）と溶存無機炭素（DIC）を時間経過とともに進化させ、海気ガス交換に基づく表面CO2フラックスを算出する。

本研究では、主に2つの分析フレームワークを利用している：

1. 感度シミュレーション： 個々のモデルパラメータ（混合層の深さ、風速、ガス移動係数、躍層の拡散率、背景となるTA/DIC、投入深度など）を、観測された自然変動と一致する値で摂動させる。これらを参照シミュレーションと比較することで、「CDR感度比（$r_{CDR}$）」を決定する。これは、感度解析におけるCDRと参照ケースにおけるCDRの比として定義される。
2. モンテカルロ・アンサンブル： 観測された海洋変動を表す統計分布からパラメータをランダムに抽出した、1,000組以上のシミュレーションペア（ベースラインと介入）を同時に実行する。このアプローチにより、物理プロセスと水塊特性の複合的な不確実性が、どのようにCDR推定値に現れるかを評価する。

シミュレーションは、不確実性の持続性に与える影響を評価するため、海洋深度200 mおよび2,000 mの設定で50年間にわたって実施される。

主な結果
本研究は、CDRの不確実性を支配する3つの明確な時間的レジームを特定している：

1. 初期の急速平衡期 ($t \ll \tau_{mix}$): 最初の数年間において、CDR推定値は、初期の海気ガス交換タイムスケール（$\tau_{air-sea}$）を制御するパラメータ（混合層の深さ、風速、ガス移動係数、および背景の炭酸塩化学系（TA, DIC, 温度））に対して非常に敏感である。この段階では、パラメータのバイアスに応じて、CDR推定値は2倍の差が生じることがある。
2. 混合層からのアルカリ度喪失 ($t \approx \tau_{mix}$): 垂直混合によってアルカリ度が表層から除去されるにつれ、不確実性は海気交換タイムスケールと混合タイムスケール（$\tau_{mix}$）の間の相互作用によって支配されるようになる。躍層の拡散率（$\kappa_{pyc}$）に関する不確実性が極めて重要となり、展開後2〜3年でCDR推定値の変動がピークに達する（最大30%）。
3. 深層混合 ($t \gg \tau_{mix}$): 長期的には、混合による深層へのタイムスケールが海気平衡のタイムスケールと同等か、あるいはそれよりも速い場合、不確率は数十年間にわたって持続する。本研究では、躍層の拡散率に関する不確実性が、長期的なCDRの不確実性を支配する主要因であることを発見した。これは、拡散率が、平衡に達したアルカリ度の割合と、未平衡の深層水が表面に戻る速度の両方を制御するためである。

定量的知見:

• 不確実性の大きさ: 最初の数年間において、CDR推定値は参照シミュレーションから最大100%（2倍の差）異なる可能性がある。20〜50年後には10〜20%に縮小するが、深海への投入については依然として重大な不確実性が残る。
• 投入深度: アルカリ度を混合層より下（たとえ20〜40 m下であっても）に投入すると、効率が著しく低下し、不確実性が高まる。一部の推定値は50年経過後も参照シミュレーションを大幅に下回ったままである。
• 純中立性のタイミング: 介入が（ライフサイクルアセスメントのコストを考慮して）純炭素マイナスになるために必要な時間は、非常に変動しやすい。参照ケースでは約1.4年とされるが、モンテカルロ・アンサンブルでは広い分布を示した。具体的には、約30%のシミュレーションは1年未満で純中立性に達するが、50%は2年以上を要し、約4%は50年以上を要する。この「ロングテール」は、主に深海への投入、または初期の平衡を妨げる急速な混合に関連している。
• 深度依存性: 2,000 mの深海では、200 mの深海と比較して不確実性が長く持続する。200 mの場合、シミュレーション期間内に水柱全体が平衡状態に達することができるためである。

意義と主張
本論文は、mCDRの推定における不確実性は単なる一時的なものではなく、内部混合のタイムスケールが海気ガス交換のタイムスケールに対して相対的に長いことにより、数十年にわたって持続し得ると主張している。著者らは、これらの結果が、支配的な物理プロセスが十分に制約されていないことが多い地域的および全球的なmCDRシミュレーションにとって、根本的な課題を特定していると述べている。

具体的には、本研究は以下の点を強調している：

• MRVへの影響: 炭素除去の経時的な推移における不確実性は、純中立性のタイミングに直接影響を与える。これは、カーボンクレジットの枠組みにおいて極めて重要である。不確実性が適切に定量化されない限り、カーボンクレジットの付与は時期尚早、あるいは不正確になる可能性がある。
• モデルの限界: 現在の海洋モデルは、乱流やガス交換のパラメトリゼーションに依存しているが、これらは大幅な空間的・時間的なバイアスを受けやすい。本研究は、mCDRにおいては、一般的な気候予測よりも、これらの特定の物理プロセス（混合とガス交換）の検証がより重要であることを示唆している。
• 展開戦略: 結果は、効率を最大化し不確実性を最小限に抑えるために、アルカリ度を厳密に表層混合層内に投入する必要性を強調している。

著者らは、自身の1Dモデルが水平方向の移流や生物学的フィードバックを除外しており、それらがさらなる不確実性を導入する可能性があることを認め、控えめな範囲設定としている。しかしながら、非生物的なシステムであっても、物理プロセスのみから生じる不確実性は十分に大きく、mCDRの介入が実際に数十年にわたって純炭素マイナスであることを保証するためには、真剣なさらなる調査が必要であると結論付けている。