本数値モデルは水平・鉛直の 2 次元モデルである. 水平方向の座標変数を 4#4, 鉛直方向の座標変数を 5#5 と表し, 時間方向の変数は 6#6 と表す.
力学的な枠組みは, 準圧縮方程式系(Klemp and Wilhelmson,1978)を用いる. この方程式系では, 予報変数を水平一様な基本場とそこからのずれに分離し, 方程式の線形化を行っている. 基本場は静水圧平衡の状態にあるものと仮定する. またガスは理想気体とみなせるものとする. 準圧縮方程式系の導出は付録 A に示す. 方程式中の変数は付録 D に示す.
Klemp and Wilhelmson(1978) では凝結性ガスや凝結物を混合比で表現している のに対し, 本モデルでは比湿で表現している. 但し本ドキュメントにおける比湿とは, 通常の気象学で用いられている比湿を拡 張し, 全密度に対する任意のガスや凝結物の密度の比を指すものとする. 主成分が凝結する惑星大気を扱う際, 非凝結ガスの密度を分母とする混合比を用 いると数値計算上の困難が生じる可能性がある. このことは様々な惑星大気を扱うことを目的とする本モデルにとって大きな問題 となりうる. そこで本モデルでは微量成分が凝結する系だけでなく, 主成分が凝結する系での 計算も行なえるよう, 凝結ガスや凝結物を比湿で表現することにする.
本モデルでは非凝結性ガス, 凝結性ガス, 雲水, 雨水(氷)の 4 つのカテゴリー を想定している.
記号 | 意味 | 内容 |
7#7 | 非凝結性ガスの比湿 | 気体の非凝結成分 |
8#8 | 凝結性ガスの比湿 | 気体の凝結成分 |
9#9 | 雲水比湿 | 落下速度がゼロである粒子で, |
大気中の雲粒に対応する. | ||
通常 100 10#10m 以下の微小な流体粒子である. | ||
11#11 | 雨水比湿 | 有意な落下速度を持つ粒子で, |
大気中の雨粒または氷粒に対応する. |
12#12, 13#13, 14#14, 15#15, 16#16 をそれぞれ非凝結成分の密
度, 凝結成分の密度, 雲水の密度, 雨水の密度, 全密度とすると, 各カテゴリー
の比湿は以下のように定義される.
17#17 | 18#18 | 19#19 | (1) |
20#20 | 18#18 | 21#21 | (2) |
22#22 | 18#18 | 23#23 | (3) |
24#24 | 18#18 | 25#25 | (4) |
各変数を基本場と擾乱場に分け, 基本場の風速 26#26, 雲水比湿と雨水比湿がゼロ
であるとみなす.
また基本場は水平一様であり, 静水圧平衡が成り立つと仮定する.
基本場の物理量に 27#27 を付し, 擾乱場の物理量に 28#28 を付すこと
にすると, 各変数は以下のように書ける.
29#29 | 18#18 | 30#30 | |
31#31 | 18#18 | 32#32 | |
33#33 | 18#18 | 34#34 | |
35#35 | 18#18 | 36#36 | |
37#37 | 18#18 | 38#38 | |
39#39 | 18#18 | 40#40 | |
41#41 | 18#18 | 42#42 | |
43#43 | 18#18 | 44#44 |
45#45 | 46#46 | 47#47 | (5) |
48#48 | 46#46 | 49#49 | (6) |
50#50 | 46#46 | 51#51 | (7) |
52#52 | 46#46 | 53#53 | (8) |
63#63 | (9) |
64#64 | 46#46 | 65#65 | (10) |
66#66 | 46#46 | 67#67 | (11) |
以下に準圧縮方程式系の時間発展方程式を一覧する. 本モデルにおける予報変数は 75#75, 76#76, 77#77, 78#78, 79#79 (80#80) である. 密度の式では乾燥成分と湿潤成分の分子量の差を考慮するが, 熱の式では考慮し ない. 凝結量は気相質量に比べて十分少ないと仮定する.
81#81 | 18#18 | 82#82 | (12) |
83#83 | 18#18 | 84#84 | |
85#85 | (13) |
86#86 | 18#18 | 87#87 | |
88#88 | |||
89#89 | |||
90#90 |
91#91 | 18#18 | 92#92 | |
93#93 | (14) |
94#94 | 18#18 | 95#95 | |
96#96 | |||
97#97 | (15) |
比湿に関して
98#98 | (16) |
106#106 | (17) |
107#107 | (18) |
本モデルでは Kessler(1969) のパラメタリゼーションと Tobie et al.(2003)
に基づくパラメタリゼーションの 2 種類が用意されている.
Kessler(1969) のパラメタリゼーションでは 4 つのカテゴリーを想定し, 微物
理素過程として以下を考慮する.
ただし, これらの量は全て正の値として定義され,
水蒸気が直接雨水に凝結する過程は無視されている.
記号 | 内容 |
108#108 | 凝結による水蒸気から雲水への変換 (condensation) |
109#109 | 蒸発による雲水から水蒸気への変換 (evaporation) |
110#110 | 蒸発による雨水から水蒸気への変換 (evaporation) |
111#111 | 併合成長による雲水から雨水への変換. 併合や水蒸気拡散により, |
雲粒子が雨粒の大きさにまで成長する (autocondensation) | |
112#112 | 衝突併合による雲水から雨水への変換. |
大水滴が小水滴を衝突併合する (collection) | |
113#113 | 雨水の重力落下に伴う雨水混合比の変化率 (precipitation) |
この微物理素過程を用いると, 生成項, 落下項, 凝結加熱項は以下のように表さ
れる.
114#114 | 18#18 | 115#115 | (19) |
116#116 | 18#18 | 117#117 | (20) |
118#118 | 18#18 | 119#119 | (21) |
120#120 | 18#18 | 121#121 | (22) |
122#122 | 18#18 | 121#121 | (23) |
123#123 | 18#18 | 124#124 | (24) |
125#125 | 18#18 | 126#126 | (25) |
def_theta, def_qa を書き直すと, 以下のようになる.
91#91 | 18#18 | 127#127 | |
128#128 | (26) | ||
129#129 | 18#18 | 130#130 | |
131#131 | (27) | ||
132#132 | 18#18 | 133#133 | |
134#134 | (28) | ||
135#135 | 18#18 | 136#136 | (29) |
141#141 | (30) |
144#144 | (31) |
145#145 | (32) |
147#147 | (33) |
149#149 | (34) |
Tobie et al.(2003) は火星大気での CO151#151 の雲物理の定式化について述べ ている. Tobie et al.(2003) では雲水を除く 3 つのカテゴリーを考える(152#152). 雲粒は拡散成長のみによって成長すると仮定し, 雲粒の併合成長は考慮しない. 微物理過程として以下を考慮する.
記号 | 内容 |
153#153 | 凝結による水蒸気から氷への変換 (condensation) |
113#113 | 氷粒の重力落下に伴う氷比湿の変化率 (precipitation) |
この微物理素過程を用いると, 生成項, 落下項, 凝結加熱項は以下のように表さ
れる.
114#114 | 18#18 | 154#154 | (35) |
118#118 | 18#18 | 155#155 | (36) |
120#120 | 18#18 | 121#121 | (37) |
123#123 | 18#18 | 124#124 | (38) |
125#125 | 18#18 | 156#156 | (39) |
def_theta, def_qa を書き直すと, 以下のようになる.
91#91 | 18#18 | 157#157 | |
128#128 | (40) | ||
129#129 | 18#18 | 158#158 | |
159#159 | (41) | ||
135#135 | 18#18 | 160#160 | (42) |
以下, 153#153 の取り扱いについて述べる.
本モデルでは単位質量の気相に含まれる凝結核の個数及び半径は空間的・時間的
に一様と仮定する.
また雲粒の半径は各格子内において空間的に一定であると仮定する.
更に雲粒は球形の凝結核を核として形成され, 雲粒自身も球形となると仮定する.
このとき
163#163 | (43) |
170#170 | (44) |
174#174 | 18#18 | 175#175 | (45) |
176#176 | 18#18 | 177#177 | (46) |
178#178 | 18#18 | 179#179 | (47) |
184#184 | (48) |
185#185 | (49) |
194#194 | (50) |
以下, 単位時間体積当たりの雲粒落下量 195#195 の取り扱いについて述べる.
Tobie et al.(2003) では雲粒落下を無視しているが, 本パラメタリゼーション
では考慮する.
195#195 は Kessler(1969) と同様に, 雲粒の終端速度 148#148 での移流
として表現する.
即ち
196#196 | (51) |
197#197 | (52) |
199#199 | (53) |
201#201 | (54) |
210#210 | (55) |
放射加熱項 100#100 は正味の上向き放射フラックス 218#218 を用いて
以下のように表される.
本モデルでは 218#218 は陽に計算せず, 100#100 は高度のみに依存する パラメタとして与える.
220#220 | 18#18 | 221#221 | |
18#18 | 222#222 | (56) |
225#225 | (57) |
227#227 | 18#18 | 228#228 | |
18#18 | 229#229 | (58) |
227#227 | 18#18 | 228#228 | |
18#18 | 230#230 | (59) |
Klemp and Wilhelmson (1978) および CReSS (坪木と榊原, 2001) と同様に,
1.5 次のクロージャーを用いて, 乱流エネルギーの時間発展方程式を表現する
.
雲が存在しない場合,
233#233 | 18#18 | 234#234 | |
235#235 | |||
236#236 | |||
237#237 | |||
238#238 | |||
239#239 | (60) |
233#233 | 18#18 | 234#234 | |
240#240 | |||
241#241 | |||
242#242 | |||
236#236 | |||
237#237 | |||
238#238 | |||
243#243 | (61) |
散逸加熱項 101#101 は, 乱流運動エネルギーの散逸項をもとに,
以下のように与える.
246#246 | (62) |