Chemistry Labs

環境・グリーン・エネルギー化学

大気化学モデル

箱モデルと三次元モデルが大気輸送・排出・化学機構・観測を組み合わせ、汚染物質や気候関連成分を説明する方法。

直観直観:移動する化学収支としての大気

モデルは分子の位置、空気による輸送・混合、化学反応による生成・除去を追跡する。単純な箱は収支台帳であり、三次元モデルは気象場に従う格子ごとにその台帳を計算する。

モードを選ぶ。箱モデルでは排出・化学・沈着・交換を追い、3D表示では地表排出と格子・鉛直層間の輸送を追跡する。

学校学校レベル:発生源・消失源・寿命

定義: 化学寿命

速度定数 kk の一次除去では、e-folding 寿命は au=1/k au=1/k である。これは特性時間であり、個々の分子の年齢とは限らない。

dCdt=P−kCCss=Pk=Pτ\frac{dC}{dt}=P-kC\qquad C_{ss}=\frac{P}{k}=P\tau

定常状態では生成と損失が釣り合う。よく混合された箱に濃度 20 ppb の成分があり、一次損失が 0.10 h⁻¹ なら寿命は 10 h、濃度を保つ生成速度は 2 ppb h⁻¹ である。

例: 2成分の箱モデル収支

A は 4 ppb h⁻¹ で排出され、0.20 h⁻¹ で除去される。A は 0.10 h⁻¹ で B に変わり、B は 0.25 h⁻¹ で除去される。定常濃度を求めよ。

解答

収支式より Ass=4/(0.20+0.10)=13.3A_{ss}=4/(0.20+0.10)=13.3 ppb。B の生成速度は 0.10A=1.330.10A=1.33 ppb h⁻¹ なので、Bss=1.33/0.25=5.33B_{ss}=1.33/0.25=5.33 ppb。変換は A の損失で B の生成であり、これを無視すると A を誤る。

大学大学:方程式と光化学

∂C∂t=−∇⋅(uC)+∇⋅(K∇C)+P−L+E−D\frac{\partial C}{\partial t}=-\nabla\cdot(\mathbf u C)+\nabla\cdot(K\nabla C)+P-L+E-D

CC は混合比または濃度、u\mathbf u は風速、KK は有効乱流拡散係数、P−LP-L は化学生成と損失の差、E−DE-D は排出と沈着の差である。単位と座標規約を統一し、境界条件で輸送問題を閉じる。

光化学箱モデルは、指定した日射・温度・圧力のもとで連成反応速度を解く。Master Chemical Mechanism(MCM)は VOC 酸化を明示的に近い形で表し、大規模な反応網の計算負荷は大きい。

例: Leighton 光定常状態

jNO2=0.008j_{NO2}=0.008 s⁻¹、kNO+O3=1.8×10−14k_{NO+O3}=1.8\times10^{-14} cm³ molecule⁻¹ s⁻¹、[NO2]=5[NO2]=5 ppb、[NO]=2[NO]=2 ppb のとき、Leighton 関係からオゾンを求めよ。

解答

光定常では j[NO2]=k[NO][O3]j[NO2]=k[NO][O3]、したがって [O3]/[NO2]=j/(k[NO])[O3]/[NO2]=j/(k[NO])。298 K、1 atm の空気は約 2.46×10192.46\times10^{19} molecules cm⁻³ で、2 ppb の NO は 4.92×10104.92\times10^{10} molecules cm⁻³。よって比は9.0、オゾンは約45 ppb。ペルオキシラジカルによる NO から NO2 への変換はこの単純関係を崩す。

上級発展:連成ソルバーと全球モデル

大気反応機構は硬い系である。ラジカル反応は非常に速い一方、貯蔵種は年単位で残ることがある。最速反応に制限される陽的刻みは非効率で、陰的 Euler は線形減衰に安定、Gear/BDF や Rosenbrock 法はヤコビアンを用いて硬い系を解く。仮定が成り立つ場合、QSSA は選択した短寿命中間体を消去する。

演算子分割では輸送・化学・放射・沈着を別々のサブステップで進める。効率的だが分割誤差を生み、演算子と時間尺度が適切なら対称 Strang 分割で形式次数を高められる。GEOS-Chem、CMAQ、WRF-Chem、CAM-chem は異なる構成と目的で気象・排出・化学・除去を結合する。

代表的な化学機構
名称用途
CB05簡略化した地域大気質化学
SAPRC光化学スモッグ機構;VOC を分類集約
MOZART全球対流圏・成層圏化学
MCM詳細な VOC 酸化;箱モデル・チャンバーで利用

排出インベントリは人為・自然起源の格子別または部門別総量を与えるが、時間プロファイル、化学種分け、煙突高度は不確かさを加える。乾性沈着は地表交換、湿性除去は雲や降水に依存する。モデル評価では単一濃度だけでなく複数の観測量を比較する。

研究研究:逆解析と評価

随伴感度とベイズ逆解析は観測から排出量や化学パラメータを推定する。OMI・TROPOMI の衛星カラムは汚染物質やメタンを制約するが、検索の平均化カーネル、雲、輸送誤差、事前仮定を考慮する必要がある。データ同化はモデル状態を更新し、逆解析は不確かな原因・パラメータを推定する。

オゾン生成効率(OPE)は、定義した気塊履歴における生成オゾン量と消費 NOx の関係を示す診断量で、普遍定数ではない。VOC 条件、希釈、沈着、集計期間が重要である。モデルと地上網・航空機・衛星の観測を対応するサンプリングと不確かさ解析で比較する。

メタン酸化は OH を消費する。メタン増加は OH を低下させ寿命を延ばす正の化学フィードバックを生みうるため、寿命式の OH 場はメタンから独立ではない。化学気候モデル相互比較(CCMI)はシナリオ下の成層圏オゾンホール回復評価にも役立つが、予測は排出と気候強制に依存する。 数値的には、kCH4+OH≈6.3×10−15k_{CH4+OH}\approx6.3\times10^{-15} cm³ molecule⁻¹ s⁻¹、質量加重平均 [OH]≈106[OH]\approx10^6 cm⁻³ では 1/(k[OH])≈51/(k[OH])\approx5 年となるが、全球収支では約9年となる。これは OH が時空間的に変動し、土壌や成層圏でも CH4 が消費されるためである。

参考文献

  • Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 3rd ed. · J. H. Seinfeld, S. N. Pandis, 2016
  • Introduction to Atmospheric Chemistry · D. J. Jacob, 1999
  • Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation · H. Bey et al., 2001
  • Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds · S. M. Saunders, M. E. Jenkin, R. G. Derwent, M. J. Pilling, 2003
  • Modeling of Atmospheric Chemistry · G. Brasseur, D. J. Jacob, 2017