1 対流のパラメータ化の位置づけ
1.1 数値モデルにおける解像度とサブグリッド過程
大気や海洋の数値モデルでは格子間隔が有限であり、対流が生む細かな雲の発達・崩壊、局所的な上昇下降、降水粒子の生成分裂といった現象を、そのまま格子で逐一計算することは難しい。その結果、格子平均の場を得るために、格子スケールより小さい運動の寄与を“サブグリッド過程”として別途表現する必要が生じる。対流のパラメータ化は、この要請に対応するため、未解像の上昇下降運動やそれに伴う水分変換、乱流混合、降水生成を、格子スケールの予報変数へ平均的に反映させる手続きである。
1.2 対流が予報・予測へ与える影響
対流は熱と水分の鉛直輸送を通じて、温度構造と湿潤分布を急速に変える。さらに雲による潜熱放出や蒸発による冷却は、安定度や風の鉛直分布に影響し、結果として降水のタイミングや分布の見積り精度にも直結する。加えて、対流は境界層における混合や放射への間接的な作用をもたらすため、日々の天気予報だけでなく、季節〜気候の長期積分でも誤差の伝播経路の中心に位置づけられることが多い。したがって、対流の扱いは“計算上の便利さ”ではなく、予測の物理妥当性を左右する要素となる。
1.3 パラメータ化の基本目標
基本目標は、格子平均量に対して、対流が与える熱・水分・運動量の平均フラックスを、適切な大きさと鉛直位置で供給することである。加えて、対流が発生しやすい条件と消えやすい条件を表現し、発生の頻度や持続の見積りを過度に偏らせないよう設計する。最終的には、モデルがもつ水平解像度・鉛直格子・力学的時間ステップの違いに依存して振る舞いが大きく変わらないよう、整合的な“閉じ込み”を達成することが求められる。
2 対流パラメータ化の理論的枠組み
2.1 格子平均量の扱い
2.1.1 熱・水分・運動量の保存と閉じ込み
対流パラメータ化は、格子セル内の平均場に対して保存則の形を保つことが重要である。具体的には、熱に対応する温位やエネルギー、総水量(あるいは水蒸気・雲水などの成分)、運動量に関わる運動の鉛直拡散・輸送を、パラメータ化成分としてモデルの更新式へ組み込む。閉じ込みの観点では、未解像の相関(例:上昇運動と水分の同時変動)が格子平均に現れるが、それらを直接計算できないため、平均フラックスとして置き換える。これにより、格子内で起きる対流の効果を“平均的に整合した形”で反映しようとする。
2.2 1次元/多次元の簡略化
2.2.1 仮想的な上昇流・下降流の表現
対流を直接追跡する代わりに、代表的な気塊運動を仮想的な上昇流(アップドラフト)と下降流(ダウンドラフト)として表す方法が広く用いられる。上昇流は凝結や凍結などの相変化を伴い、潜熱の放出を通じて周囲との熱差を変える。下降流は蒸発冷却や沈降により、周囲の安定度にフィードバックしうる。1次元の枠組みでは、これらの流れを鉛直方向の輸送として表現し、気塊の速度や水分量、エントレインメント(周囲の混合)といった量をパラメータ化して閉じ込める。多次元に近い表現を採る場合でも、基本思想としては代表気塊を通じた平均輸送を再現しようとする。
2.3 時間発展の記述方法
2.3.1 発生から消滅までの定式化
対流は発生し、成長し、やがて弱まって消滅するという時間発展を持つ。パラメータ化では、この過程を連続的な時間積分に組み込むため、発生条件を判定し、その後の輸送や水変換を時間発展方程式として更新する設計が採られることが多い。発生から消滅までの定式化には、対流の“寿命”をどのように見積もるか、上昇・下降双方の強度をどのように増減させるか、供給される不安定エネルギーや水分が尽きたときにどのように減衰させるかが含まれる。結果は、モデルが持つ時間解像度と整合するよう、更新頻度や安定性の観点から工夫される。
2.4 放射・降水・境界層との結合
対流は雲を形成し、放射過程(長波・短波)の吸収・散乱を変えるため、放射と相互作用する。さらに降水が落下すれば、地表付近の蒸発や冷却により境界層の乱流混合や安定度が変わる。パラメータ化は、雲水・氷相・降水粒子の分布に関わる変数を通じて、これらの物理過程に影響を及ぼす必要がある。結合の設計では、対流が更新する時間スケールと放射計算や境界層過程の時間更新との整合が問われるため、サブステップ化や反復手法などの実装上の工夫がしばしば用いられる。
3 主要パラメータとその設定
3.1 発生条件(トリガー)のパラメータ
3.1.1 不安定度・湿潤性・収束の指標
対流は任意に起きるのではなく、周囲に対して上昇が有利になる条件が必要である。ここで用いられる代表的な指標には、不安定度(鉛直方向の安定性が破れつつある度合い)、湿潤性(凝結開始や潜熱放出が起こりやすい水分状態)、収束(運動の集中を通じて上昇を誘導しやすい力学状況)がある。これらはしばしば組み合わせて、ある閾値を超えた場合に発生が“許可”されるよう設計される。指標の選択と閾値の置き方は、対流の初期頻度や地域差に影響し、結果の現実性を左右するため慎重な調整対象となる。
3.2 鉛直輸送を決めるパラメータ
3.2.1 エントレインメント(混合)係数
エントレインメントは、上昇する気塊が周囲の空気を取り込む割合や混合の強さを表す概念である。この混合が強いほど、気塊の熱・水分特性は周囲に近づきやすくなり、上昇が持続しにくくなる場合がある。逆に混合が弱いと気塊はより“自律的”になり、強い上昇や大きな水輸送を生みやすい。係数は鉛直高度依存や条件依存として定義されることが多く、気塊の安定度や湿度構造、あるいは乱流の見積りと結び付けられる。
3.2.2 輸送率・フラックスの計算方法
輸送率は、上昇流と下降流それぞれの速度や気塊の断面規模に対応する量から決まる。フラックスは、輸送率と気塊が運ぶ諸量(熱量や水分混合比、運動量成分など)との積として表され、格子平均量の更新に用いられる。計算では、気塊の速度をどのように与えるか、気塊内の水分や温度が周囲とどう相対的に変化するか、混合による効果をどこに含めるかがポイントになる。また、フラックスの符号や大きさが不自然に振れないよう、物理制約(負の水量を避ける等)を反映する設計も重要である。
3.3 濡れたプロセスのパラメータ
3.3.1 蒸発・凝結・凍結の扱い
対流内では水の相変化が進行し、熱収支と水収支を同時に決める。凝結は上昇で冷却された空気が水蒸気を雲水へ変える過程として表れ、凍結は温度低下によって氷相が現れる場面で関与する。下降流では蒸発が起こりやすく、これが周囲を冷やして安定度を変える。パラメータ化ではこれらのプロセス速度を、雲微物理の詳細を省略した形で速度係数、しきい値、効率のパラメータとして表現することが多い。相変化の扱いは降水生成とも結びつき、結果として雲粒子の量や寿命推定を通じて降水量へ影響する。
3.4 降水の生成と落下のパラメータ
3.4.1 雨滴・雪・霰に関する簡略化
降水形成は、雲中での衝突・合体、凍結過程、成長に関する微視的機構を含むが、パラメータ化では粒径分布を直接追わないことが多い。その代わり、雨滴・雪・霰といったカテゴリに分類し、それぞれの生成・移行・降下速度を簡略化した形で扱う。落下速度は重力沈降により大きくなるため、軌道が短時間で高度を移ると地表の雨に対する寄与が変わる。分類の仕方や変換効率の置き方は、同じ総水量でも“どの高度から降るか”の違いを通じて降水の鉛直構造に反響する。
3.5 停止条件(マダニング)のパラメータ
3.5.1 限界不安定度・餌切れの評価
対流が持続するには、上昇のためのエネルギーと水分の供給が必要である。停止条件(マダニング)は、これらの供給が不足したときに対流を弱めるための規則を指す。具体的には、ある程度以上に安定化した状態が形成されると上昇が維持できないという“限界不安定度”の考え方や、気塊が運べる水分の余裕が失われ“餌切れ”になるという見積もりが用いられることがある。停止のタイミングは、対流の平均寿命や再発生までの間隔を決め、降水の時間変動の振る舞いに波及する。
4 濃淡度とモデル挙動の制御
4.1 風の鉛直シアと組織化
4.1.1 複数セル/セル群の表現
対流は単発の気塊としてだけでなく、組織化した構造として現れる場合がある。風の鉛直シアは上昇流や下降流の相対配置を変え、セルが連鎖するような形や、複数の領域が時間的・空間的にまとまる形に影響する。パラメータ化では、解像できない組織化を平均量として表現するため、複数の代表気塊や“セル群”の概念を導入し、セル間の相互影響を簡略化して反映する場合がある。これにより、水平スケールと鉛直輸送の結びつきが自然な範囲に収まるよう工夫される。
4.2 対流抑制と抑制解除
4.2.1 乾燥域・安定成層の影響
下降流による蒸発冷却は、周辺に乾燥層や安定層を形成し、次の対流発生を抑えることがある。さらに大規模な安定成層が存在すると、上昇が潜熱放出で補える前に持ち上がりにくくなる。抑制解除は、乱流混合や収束、熱力学的な再不安定化によって可能になる。パラメータ化では、抑制を単なる一括停止ではなく、場の状態に応じた発生許可の変化として表す工夫が重要である。これにより、短時間の抑制とその後の再発生が、平均的に妥当な頻度で起こる。
4.3 パラメータ感度と調整(キャリブレーション)
4.3.1 設定値がもたらす系統的偏り
対流パラメータは、発生閾値、混合強度、相変化効率、停止規則など多数の要素からなる。その設定を変えると、降水の総量、降水の開始時刻、雲量、鉛直温度・湿度プロファイルに系統的偏りが生じうる。たとえば混合を強めると対流が弱まりやすくなり、降水の頻度が減る一方で大雨イベントが相対的に強調されるなどの形が現れることがある。したがって、キャリブレーションでは単一の評価指標に頼らず、複数の観測制約やモデル診断を同時に考慮する必要がある。調整は地域ごとの気候特性と結びつくため、過度な一般化は避ける設計が望ましい。
4.4 計算安定性とタイムステップ依存性
対流パラメータ化は、時間積分の中で強い非線形性を持つため、タイムステップの大きさにより挙動が変わる場合がある。発生判定やフラックスの更新がタイムステップに対して遅れると、対流の開始が遅れたり、逆に短時間で過大な輸送が発生して不安定化したりすることがある。安定性確保のため、対流成分に対してサブタイムステップや制限(上限・下限)を設ける設計、あるいは反復的な結合手順を導入する場合がある。結果として、物理妥当性と計算効率の両立が課題となる。
5 検証・評価と改善の道筋
5.1 観測データによる評価
5.1 雲量・降水・鉛直プロファイル
対流パラメータ化の妥当性は、雲量や降水量の統計だけでなく、鉛直温度・湿度分布、雲の高度構造、降水の生成高度などの“中間的診断”によって評価される。観測には地上雨量、衛星による雲・降水の推定、レーダーによる鉛直断面、ゾンデによる鉛直プロファイルなどがある。特定の指標が一致しても他の制約と矛盾する場合があるため、多面的な検証が必要となる。さらに季節や気候帯ごとに評価を分けることで、調整が偏った領域に依存していないかを確認しやすくなる。
5.2 再解析データ・モデル比較
5.2.1 同化と対流表現の関係
再解析データは観測同化によって生成されるため、対流に関する誤差が同化の仕方を通じてモデル構造に影響する場合がある。そのため、再解析を“真値”として扱うと、評価が歪むことがある。たとえば同化の観測種が雲や降水に敏感であれば、対流が強い場面では同化が特定の物理表現へ引っ張られる可能性がある。評価の際は、同化システムの特徴と、比較対象モデルの対流評価軸を慎重に整合させる必要がある。
5.3 大規模LESやクラウド解像モデルとの比較
5.3.1 サブグリッド統計の導出
大規模LESやクラウド解像モデルは、対流の多くをより直接的に表現できるため、パラメータ化の検証に有用である。比較では、上昇下降の頻度、鉛直輸送の統計、混合の強さ、降水形成の鉛直分布など、パラメータ化が暗黙に仮定する量に対応する統計量を取り出す。こうしたサブグリッド統計は、係数の設定根拠として使われたり、パラメータ化の仮定が現実の統計と整合しているかを点検するために利用される。目的は“置き換え”ではなく、パラメータ化の設計思想を現象に照らして修正することである。
5.4 改良の代表的アプローチ
5.4.1 物理の拡張、パラメータ化の階層化、学習の利用(一般論)
改良は複数の方向で進む。第一に、雲微物理や混合、相変化の簡略化が不十分な場合には物理過程の表現を拡張する。第二に、単一の係数で全てを説明するのではなく、状況に応じて異なるスケールや状態の寄与を階層的に結合する設計が検討される。第三に、統計的手法や機械学習を一般論として利用し、観測や高解像度モデルから得られる関係を補助的に導入することもある。ただし、学習に依存するほど過去のデータ範囲外での外挿が難しくなるため、物理制約を併用する設計が重要になる。
6 記述体系と代表的なモデルの分類
6.1 直接フラックス型と式ベース型の違い
分類の一つは、フラックスを直接与える方式か、気塊の運動や相変化を含む式に基づいて結果としてフラックスを得る方式か、という観点である。直接フラックス型では、格子平均へ追加される熱・水分のフラックスを統計的に見積もり、保存則を満たす形で適用する。式ベース型では、上昇流・下降流の速度や混合、相変化の進行を計算し、その積として輸送が決まる。どちらが優れているというより、目的と計算制約に応じた設計選好として理解される。
6.2 単一カラム型と複数カラム型の違い
単一カラム型では、格子セル内の平均的な代表状態を用い、対流の鉛直構造をそのカラムに対して計算する。複数カラム型では、同一セル内に複数の代表的気塊や状態を並列に扱い、確率的な多様性をより強く反映させることがある。これにより、同じ平均条件でも実際には異なる対流モードが共存しうるという現象を、より表現しやすくなる。ただし計算コストが増える傾向があるため、実装上の折り合いが課題となる。
6.3 「準定常」仮定と「時間依存」仮定
準定常仮定は、対流の変化がある時間内では局所的な平衡に近いという考え方に基づく。これにより、対流の内部過程を簡略な関係式として更新できる場合がある。一方、時間依存仮定では、対流の成長・減衰を時間発展として追跡し、非平衡な段階の影響を反映する。実際の現象は両者の中間的性格を持つため、モデル設計では時間スケールの整理と、計算上の手間とのバランスが重要になる。
6.4 自己調整型・統計調整型の考え方
自己調整型では、対流パラメータ化が与えるフラックスが環境を変え、それが次の対流発生や停止に直接影響することを通じて、長期的にある程度の統計特性が安定化することを狙う。統計調整型では、観測や高解像度計算に合うように、特定の統計量が目標範囲に入るように調整を行う考え方が中心になる。いずれも“現象の一部を平均として閉じ込める”点は共通であり、調整の主体がどこにあるかが違いとして現れる。
7 実装上の注意点
7.1 入力変数の品質管理
対流パラメータ化の発生判定や輸送計算は、格子平均の温度・湿度・風・収束など入力に強く依存する。したがって、入力変数が数値的にノイズを含んでいる、あるいは鉛直配置の不整合があると、発生の閾値判定が頻繁に揺れたり、フラックスが不連続になることがある。実装では、鉛直補間、物理量の範囲制限、診断に用いる平均化の手順など、入力側の品質確保が結果の再現性に直結する。
7.2 パラメータの単位・スケーリング
パラメータの定義はモデル座標系や物理単位系と密接である。係数が次元を持つ場合、スケーリングの置き方(高度依存、密度依存、温度依存など)を誤ると輸送量や相変化が桁違いに変化する恐れがある。さらに、水平解像度や鉛直格子間隔が変わったとき、パラメータがどのように再解釈されるべきかも重要である。単位チェックと整合性テストは、運用段階での“静かなバグ”を防ぐために欠かせない。
7.3 他物理過程との整合
7.3.1 境界層乱流・放射との干渉
対流は境界層乱流と結びつき、雲による放射加熱・冷却は温度構造を変え、乱流の生成にも影響する。したがって、対流更新と境界層・放射更新の順序、時間同期の取り方が結果へ影響する。たとえば放射が先に温度を変えるのか、対流が先に雲水量を更新してから放射計算を呼ぶのかで、微妙な差が蓄積される可能性がある。整合性は計算結果の安定性と物理一貫性の両面に関わるため、結合スキームの設計とテストが不可欠である。