Solution of the Advection-Diffusion Equation in Two Dimensions using the Deposition Velocity

技術要約：沈着速度を用いた二次元移流拡散方程式の解法

問題提起
本研究は、不安定な大気条件下における二次元定常状態移流拡散方程式（ADE）の解析解を扱うものである。主な目的は、汚染物質の拡散、特に中心線における風向横方向積分正規濃度に焦点を当てることである。本モデルは、地上レベルにおける沈着速度境界条件を組み込み、垂直方向の渦拡散係数および風速の両方に、垂直高度の関数としてのべき乗則関数を用いている。本研究は、この解析的手法を、デンマークのコペンハーゲン・フィールド実験から得られた硫黄ヘキサフルオライド（SF6）の観測濃度データに対して検証することを目的としている。

手法
著者らは、支配的なADEを解くために変数分離法を採用している。数学的枠組みは以下のように定義される：

• 支配方程式： 二次元ADEは、風速 $u(z)$ および垂直渦拡散係数 $k_z(z)$ が垂直高度のべき乗（それぞれ $z^p$ および $z^q$）として表現される形式で定式化されている。
• 境界条件：
  1. 上端境界（混合層高度 $h$）： ゼロフラックス条件（$k_z \frac{\partial C}{\partial z} = 0$）。
  2. 地表面境界 ($z=0$)： 沈着フラックス条件。ここで、フラックスは沈着速度（$V_d$）と濃度の積に等しい（$k_z \frac{\partial C}{\partial z} = -V_d C$）。
  3. ソース項： スタック高さ $h_s$ に位置する連続点源であり、ディラックのデルタ関数によって表される。
• 解析的導出：
  • 方程式は空間成分 $X(x)$ と $Z(z)$ に分離される。
  • 垂直成分は、可変係数を持つシュトゥルム＝リウヴィル型の微分方程式を導く。
  • 指数の特定値（例：$q=0.85$）を仮定し、変数変換を用いることで、方程式は一般化ベッセル方程式へと簡約される。
  • 垂直方向の濃度プロファイルに対する解は、第一種および第二種ベッセル関数を用いて導出され、これらは半次ベッセル関数として三角関数（正弦および余弦）へと簡略化される。
  • 固有値および固有関数は、境界条件を適用することによって決定され、濃度 $C_y(x, z)$ の級数解へとつながる。
  • 係数は、風速プロファイルから導かれる直交条件を用いて算出される。

主な貢献

• 解析的枠組み： 本論文は、不安定条件下における2D ADEに対する、沈着速度境界条件（これは拡散モデリングにおいて複雑な要因となることが多い）を明示的に組み込んだ、閉じた形式の解析解を提示している。
• べき乗則パラメータ化： モデルは風速と渦拡散係数のべき乗則表現を利用しており、これにより異なる安定度クラス（指数 $p$ および $q$ によって定義される）にわたる柔軟性を実現している。
• 検証： 解析モデルは、様々な安定度クラス（A, B, C, D）および下流距離をカバーする、よく知られたコペンハーゲンSF6データセットに対して厳密にテストされている。

結果
本研究では、予測された風向横方向積分正規濃度を、コペンハーゲン実験における観測されたSF6濃度と比較している。

• データ比較： 表2および図1は、様々な下流距離における複数のラン（Run 1からRun 9）について、観測値と予測値の直接的な比較を示している。
• 統計的性能： モデルの性能は、標準的な統計指標（Hanna, 1989）を用いて評価されている：
  • 正規平均二乗誤差 (NMSE)： 0.07（低い誤差を示す）。
  • 分数バイアス (FB)： 0.003（無視できるバイアスを示す）。
  • 相関係数 (COR)： 0.82。
  • 2倍以内の範囲に入る割合 (FAC2)： 0.997。
• 視覚的一致： 図1は、観測された濃度プロファイルと予測値との間の高い一致度を示している。図2は、予測値が観測データの2倍以内の範囲に収まっていることを裏付けている。

意義および主張
著者らは、提案された解析モデルが、観測されたフィールドデータに対して「高度な閉じた一致（high degree of closed agreement）」を達成していると主張している。具体的には、FAC2の値が0.997であり、バイアスがほぼゼロであることから、予測された濃度は中心線における観測濃度と**100%**一致したと述べている。本研究は、変数分離法と風および拡散のべき乗則パラメータ化を組み合わせた手法が、複雑な数値シミュレーションを必要とせずに、不安定な大気条件下での汚染物質拡散を予測するための堅牢な解析ツールを提供することを結論付けている。本研究は、不安定境界層における汚染物質輸送の本質的な物理現象を捉える上で、変換ベースおよび分離ベースの解析解の有用性を強化するものである。