雲・エアロゾル・放射フィードバック精密評価 (kuba_02.12.18

共生課題２
２．温暖化・大気組成変化相互作用モデル開発
温暖化 - 雲・エアロゾル・放射フィードバック精密評価
久芳奈遠美
KUBA Naomi
地球フロンティア研究システム
水循環予測研究領域・雲降水グループ
IPCC 2001
対流圏エアロゾルの間接放射強制力の見積もりは
いまだに不確定性の大きい問題であり、その主な
要因はエアロゾルと雲の関係が不確定であること
による。
雲凝結核（ＣＣＮ）が雲の光学特性および降水形成効率に及ぼす影響
数値実験
微物理学的雲モデル
活性化凝結 / 蒸発併合重力落下
CCN (エアロゾル)
粒径分布
化学組成
上昇流速度
雲粒
気候変動
粒径分布
放射収支
光学特性
水循環
降水形成効率
観測
観測
ＣＣＮカウンター
２DC FSSP 雲レーダ
ドップラーライダ
衛星データレーダ
長期・全球観測
雲凝結核（ＣＣＮ）が雲の光学特性および降水形成効率に及ぼす影響
数値実験
微物理学的雲モデル
活性化凝結 / 蒸発併合重力落下
CCN (エアロゾル)
粒径分布
化学組成
上昇流速度
雲粒
気候変動
粒径分布
放射収支
光学特性
水循環
降水形成効率
観測
観測
ＣＣＮカウンター
２DC FSSP 雲レーダ
ドップラーライダ
衛星データレーダ
長期・全球観測
GCM におけるエアロゾルの間接放射強制力の評価
1. CCSR/NIES の場合
放射計算に用いる有効半径 re を雲粒数密度 Nd から計算する。
Nd = e Na Nm / (e Na + Nm )
Na ⇒ 0 : Nd ⇒ e Na
Na ⇒ ∞ : Nd ⇒ Nm (constant)
Nd：雲粒数密度
Na：総エアロゾル数密度(sulfate, carbon, sea salt)
2.
ECHAM GCMの場合
Nd = w Na / ( w + a Na )
Ghan et al. (1993)
w = wm + c TKE ½
Lohmann et al. (1999)
w : 雲内上昇流速度
wm : grid 内平均上昇流速度
TKE ：turbulent kinetic energy
FRSGC での研究成果
詳細雲モデルを開発し、数値実験によって雲の光学的性質
や雲粒数密度を予測するパラメタリゼーションを開発。
Kuba et al., 2003 :
Parameterization of the Effect of Cloud Condensation
Nuclei on Optical Properties of a Non-precipitating
Water Layer Cloud
J. Meteorol. Soc. Japan
Smax の予測式：Nc(X%)の関数による近似
Smax (%)
1
0.1
S max = H N c (0.075%) I
Vbase
m s-1-0.230
S max = 0.710
N =0.24
c (0.075%)
H = 1.16V base 0.344
-1-0.262
S max = 0.559
N=0.12
(0.075%)
Vbase
m
s
c
I = - 0.176V base -0.187
-0.298
S max = 0.441
N c=0.06
(0.075%)
Vbase
m s-1
0.01
1
10
100
-3
N c (0.075%) (cm )
1000
1000
N d = min[a{N c ( S=S max ) - b )( N c ( S = S max ) - c } , N c ( S = S max )]
Nd をＣＣＮスペクトルから予測する近似式
a=-0.0231V
-0.0108V
-0.00180
2
base
base
b=70.0V base -12.9
c=8420V
c maxbase +278
Nd = min [ E(N (S
) - F)(Nc(Smax) - G) , Nc(Smax) ]
for N c (S=S max ) < 0.5(b+c)
Vbase= 0.24 m s-1
Nd
Vbase= 0.12 m s-1
500
Vbase= 0.06 m s-1
Nc(Smax) が大きい時
Vbase が小さい時
Nd = Nc(Smax)
にならない
0
0
500
N c (S=S max )
Nc(Smax) (cm-3)
1000
雲の光学的厚さの近似式
100
t =tA=aN
Ndd Bb
a=0.121LWP
0.702
Optical Thickness t
0.702
A = 0.121 LWP
0.0538
b=0.274LWP
B = 0.274 LWP 0. 0538
10
Nd と LWP の関数として表す
LWP = 100 g m-2
(300 m)
LWP = LWP=100g
70 g m-2 m-2
(250 m)(300 m)
-2
LWP = LWP=
44 g m70g
m-2
(200 m)(250 m)
LWP= 44g m-2
(200 m)
Nd および LWP がわかれば、
τを予測できる。
τ および LWP がわかれば、
Nd を推測できる。
1
10
100
1000
Number m
concentration
of cloud droplets
雲底から100
での雲粒数密度
Nd (cm-3)
at 100 m above cloud base N d (cm-3)
雲粒の有効半径の近似式
Nd の方が Re より指標として有効
100
100
-D
Re = C N
R ed= aN d -b
Effective Radius Re (mm)
a=6.41Z 0.380
C = 6.41b=0.288Z
Z 0.380 0.0254
D = 0.288 Z 0. 0254
雲底からの高度
10
10
Z=275 m
Z=225 m
Z=175 m
Z=125 m
Z= 75 m
Z= 25 m
11
10
10
100100
Number Concentration of Cloud Droplets
-3
雲底から100
m での雲粒数密度
at 100 m above
cloud base N d (cm ) Nd
1000 1000
(cm-3)
CCN数密度の近似式
Nd と Vbase の関数として表す
1000
b
K
NcN(0.2%)
= J dN
c (S=0.2%)=aN
d
a=2.19V base
1.64
Nc(S=0.2%)
1.64
-0.359
J =b=0.714V
2.19 Vbase
base
K = 0.714 Vbase-0.359
Vbase=0.24 m s-1
Vbase=0.12 m s-1
Vbase=0.06 m s-1
100
Nc(X%) と Vbase がわかれば、
Nd を予測できる。
Nd と Vbase がわかれば、
Nc(X%) を推定できる。
10
1000
1010
100100
1000
-3
Number Concentration of Cloud Droplets at 100 m above Cloud Base N -3
d (cm )
雲底から100 m での雲粒数密度 Nd (cm )
開発した手法 (パラメタリゼーション)
Revised
CCN 過飽和度スペクトル
雲底付近の
上昇流速度
Vbase
Nc(0.075%)
New
Nc(0.2%)
Nc(Smax)
Smax
鉛直積算雲水量
LWP
雲粒数密度
Nd
雲の光学的厚さ・雲粒の有効半径
t, Re
研究方針
ＧＣＭに組み込むための、エアロゾル（ＣＣＮ）の気候への
影響を評価できるをパラメタリゼーションを開発する。
既存のGCMにおけるパラメタリゼーションの検討
手持ちの成果
を発展させて反映させる
非降水性の層雲の光学的厚さ t を
雲粒数密度 Nd と鉛直積算雲水量 LWP から導く式
Nd をＣＣＮの情報と上昇流速度から予測する式
問題点：雲内の上昇流、雲の LWP を格子点値からどう
やって算出するか？

Download Report