スカイラジオメーターによる東京及び小笠原に おける対流圏エアロゾルの光学的特性 に関する研究 東京商船大学 商船学研究科 交通電子機械工学専攻 齋藤 泰治

1. スカイラジオメーターによる東京及び小笠原に　　　おける対流圏エアロゾルの光学的特性　　に関する研究スカイラジオメーターによる東京及び小笠原に　　　おける対流圏エアロゾルの光学的特性　　に関する研究 東京商船大学　　商船学研究科 交通電子機械工学専攻 齋藤　泰治

2. 1.　研究の背景及び目的 2.　観測装置及び観測原理 　　・skyradiometerの観測動作 　　・観測場所 　　・観測から得られる 　　　各種パラメーター 　　・観測原理 3.　解析結果及び考察 ・光学的厚さの波長依存性についての考察 　　・本学上空における観測 　　・父島上空における観測 　　・2003年3,4月の黄砂と5月のシベリア森林火災による影響 4.　まとめ 発表内容

3. 研究の背景及び目的 ・地球環境、大気環境への関心が高まり、人間活動による排出物や、気候変動による砂塵嵐、森林火災の多発、拡大化によって、大気中を浮遊する微小粒子（エアロゾル）による影響が懸念されている。 ・エアロゾルの光学的な特性や粒径分布などの特徴、季節変化、経年変化を知ることで、その気象、気候に与える影響を知る基礎となる。 ・本研究では、東京海洋大学、越中島キャンパス（以下、本学）及び、小笠原諸島、父島での大気エアロゾルの観測から、その光学的特性、粒径分布などの解析を行った

4. ・本学は、東京都心に位置し、ここでのエアロゾルの光学的特　性やその挙動を知ることは、大都市域においてエアロゾルが　気象、気候に与える影響を知る基礎となる。・本学は、東京都心に位置し、ここでのエアロゾルの光学的特　性やその挙動を知ることは、大都市域においてエアロゾルが　気象、気候に与える影響を知る基礎となる。 ・父島における観測は、国立極地研究所によって行われ、　　　　平成16年度、共同一般研究 「スカイラジオメーターによる極域及び中低緯度域における　　エアロゾルの光学特性に関する研究」 の一環として、観測データの提供をいただき、解析を行った。

5. ・父島は、東京から南方約1000kmに浮かぶ亜熱帯性気候の　島で、中緯度偏西風帯と貿易風帯の境に位置するため、アジ　ア起源の気団と太平洋中心部の清浄な気団に、季節によって交代で支配される。・父島は、東京から南方約1000kmに浮かぶ亜熱帯性気候の　島で、中緯度偏西風帯と貿易風帯の境に位置するため、アジ　ア起源の気団と太平洋中心部の清浄な気団に、季節によって交代で支配される。 ・そのため、大陸性（黄砂や人為汚染）のエアロゾルや、海洋　性バッククグラウンドの特徴を捉えることが期待される。 ・2ヶ所で、同時期に観測、解析を行うことで、観測点間のエア　ロゾルの光学的特性や、挙動の差異を比較することができる。 ・また、エアロゾルの広域における影響や輸送過程に伴うエアロゾルの量的、光学的特性の変化なども捉えると期待される。

6. 2.観測装置及び観測原理 ・スカイラジオメーターは、太陽直達光と、角度別散乱光の放射輝度を測定する、分光放射計である。 ・設定エアマス（0.25）毎もしくは設定時間（10分）毎に自動測定を行う。 ・散乱角度毎の全天放射輝度分布の測定から、エアロゾルの散乱の寄与、粒径分布など多くの情報を得ることができる。 本学におけるSkyradiometer観測風景

7. 観測装置仕様 （株）プリード社製　 　　　Skyradiometer POM-01 視野半角　　　　　　　　　 0.5° 測定散乱角度　　　　　　 3,4,5,6,7,10,15,20,25,30,40,50,60,70,80,90,100, 110,120,130,140,150,160° 観測波長　　　　　　　　 315,400,500,670,870,940,1040 [nm]　（本学） 315,400,500,675,870,940,1020 [nm]　（父島） 波長選択　　　　　　　　 フィルターホイール方式 光検出器　　　　　　　　 シリコンフォトダイオード 駆動方式　　　　　　　　 パルスモーター駆動方式 動作角度　　　　　　　　　 350°（方位角） 150°（天頂角） 角度分解能　　　　　　　 0.0036° サンセンサー　　　　　　　4 素子シリコンセンサー

8. 光学系（分光干渉フィルター、検出器） 鏡筒 駆動系（太陽追尾装置） 降雨センサー サンセンサー PC（制御及びデータ蓄積用） Skyradiometer本体 制御ユニット Skyradiometer 構成概念図

9. Almucantar plane （水平スキャン） f q0 Q f =180 f =0 天頂角q > 15° エアマスm < 3.0 散乱角Q 180°まで エアマスm ≧ 3.0散乱角Q 30°まで Skyradiometer の観測動作

10. 天頂角q < 15° 散乱角Q 150°まで Q q Principal plane （垂直スキャン） Skyradiometer の観測動作

11. 観測場所 ・東京海洋大学、越中島キャンパスでの観測 観測場所　：　東京海洋大学　越中島キャンパス内　 　　　　　　　　　　二号館屋上（海抜約35m） 　　　　　　　　　　（ 北緯35.664°東経139.796°） 　　　　　　　　　　気圧計測 　　　　　　　　　　同内　第三実験棟屋上 　　　　　　　　　　（緯度､経度､高度情報は、輝度観測点と同　　　　　　　　　　　　一とする） 観測期間(解析対象期間）　：　2002年 8月～2004年12月

12. ・小笠原諸島、父島での観測 観測場所　：　東京都小笠原村父島字桑ノ木山　 JAXA　小笠原追跡所　管制棟 （海抜約220m） 　　　　　　　　　　（ 北緯27.067°東経142.217°） 　　　　　　　　　　気圧計測 　　　　　　　　　　東京都小笠原村父島西町 　　　　　　　　　　気象庁　父島気象観測所 （海抜3m） 　　　　　　　　　　（ 北緯27.001°東経142.183°） 　　　　　　　　　　*（財）気象業務支援センターにより一般に提供 観測期間(解析対象期間）　：　2003年 2月～2004年11月

13. 気象庁　父島気象観測所 JAXA　小笠原追跡所 観測地点の位置関係

14. ・複素屈折率 　　物質の屈折率で、nreが散乱成分、nimは吸収成分を表す。 • ・光学的厚さ　t • 　光が密度rの物質中を距離dz通過する際の消散の割合。 • 　地上0mから高度z までの大気の光学的厚さは次式で定義される。 （2.1） 観測から得られる各種パラメーター

15. ・体積粒径分布v(r) 　　単位面積における大気気柱中のエアロゾルの内、ある粒径範囲dlnrの範疇にあるエアロゾルの体積量 ・単一散乱アルベド　w (l) 　　エアロゾルによる、太陽光の消散の内、散乱による割合。 　　消散は波長と物質により決まる定数 → 消散係数 kext 　　散乱部分 → 散乱係数 kscat、 吸収部分 → 吸収係数 kabs （2.2）

16. ・Ångstrom指数 エアロゾルの光学的厚さは、波長に依存し、Ångstrom(1961,1964)の経験式 で示される。ここで、aを波長指数（ Ångstrom指数）、bを混濁係数という。 Ångstrom指数は、エアロゾルの粒径分布と関係付けられ、 一般に微小粒子が多いと大きく、粗大粒子が多いと小さくなる。 （2.3） （2.4）

17. 観測原理 各散乱角Qにおける放射強度 E(Q) と直達太陽放射強度 F から DW;　装置の視野立体角 m;　大気路程（エアマス）、1/cosq0で近似できる。 と規格化された相対的強度が定義できる(Nakajima et al, 1994)。 (2.5)

18. 一次散乱寄与であるb は また、光学的厚さは粒径分布 n(r) を持つエアロゾル粒子の消散の効果の和として、 （K は散乱光の振幅、位相等を決定する関数、Qextはその　粒子による散乱の効率を示す散乱効率因子である。厳密　な計算はMie散乱理論に基づき求められる。） (2.6) (2.7)

19. Input data DW; solid view angle A; ground albedo ; complex refractive function rmin,rMAX ; Radius interval b (n)⇒v(n)(r) v(n)(r) ⇒t(n)(l) , R(n)(l,Q) Up date b e < 0.1% or n > 20 no yes Output v(r) ; volume size distribution t(l) ; Optical thickness P(l,Q) ; single scattering phase function w (l) ; single scattering albedo Skyrad.pack コード(Nakajima et al,1994) における解析フローチャート

20. 3.結果及び考察光学的厚さの波長依存性に関する考察3.結果及び考察光学的厚さの波長依存性に関する考察 ・観測結果の議論の前に、 Ångstrom指数に関して考察した。 微小粒子や粗大粒子が顕著に卓越する場合、Ångstromの経験式（1次式）からの逸脱が大きくなり、十分に近似が成されない(Eck at el, 1999)。そこで、 と表される2次曲線による近似を考えることとする。 （3.1）

21. 微小粒子の卓越時(森林火災;’03.05.21) 粗大粒子の卓越時(黄砂 ; ‘04.04.17) の体積粒径分布　　　　　(日平均)

22. 微小粒子および、 粗大粒子の卓越時 での光学的厚さの波長依存性 ln(t) vs ln(l)上で、 微小粒子が多いと近似曲線は、凸 粗大粒子が多いと近似曲線は、凹 となる傾向を示す。

23. ここで、式（3.1）において a* = -a1,g = -a2 を定義する（Schmid at el,2003）。g は一次曲線からの逸脱を示し、g = 0　ならば、式（2.4）＝式（3.1）（1次式）である。 ln(t) vs ln(l) における波長依存性の振る舞いにおいて、 　　微小粒子が多いと近似曲線は、凸 　　粗大粒子が多いと近似曲線は、凹 であったから、gは 　　微小粒子が多いと大きい 　　粗大粒子が多いと小さい

24. 2003年に観測された微小粒子の卓越時と粗大粒子の卓越時における、a* とg の相関 微小粒子が多い 　　　　　→　 a* 、g 大 粗大粒子が多い 　　　　　→　 a* 、g 小

25. a*、gの導入により、 ・光学的厚さ（絶対量）に左右されず、粒径分布のみを議論す　ることができる。 ・ Ångstromの経験式からの逸脱が大きい場合においても、　　十分に波長依存性を議論できる。 ・これまでと同様の簡便さで、 Ångstrom指数aより多くの情　　報を得ることができると考えられる。 従って、本研究において、これまでのÅngstrom指数aに加え、 a*、gも用いて議論した。 a*、g により推定される結果は、粒径分布やÅngstrom指数aにおける解析結果とも矛盾が無かった。

26. 本学上空における観測 2002年8月から2004年12月におけるt,a,wの時系列変化

27. 2003年及び2004年におけるaとt の相関の季節変化

28. 2003年及び2004年におけるa*とgの相関の季節変化2003年及び2004年におけるa*とgの相関の季節変化

29. 父島における観測 2003年3月から2004年11月におけるt,a,w の時系列変化

30. 2003年及び2004年におけるaとt の相関の季節変化

31. 2003年及び2004年におけるa*とgの相関の季節変化2003年及び2004年におけるa*とgの相関の季節変化

32. 2003年3,4月の黄砂と5月のシベリア森林火災による影響2003年3,4月の黄砂と5月のシベリア森林火災による影響 本学及び父島における2003年春季のtの時系列変化

33. 黄砂飛来時（2003年3月26日及び4月20日）における黄砂飛来時（2003年3月26日及び4月20日）における NOAA HYSPLIT MODELによる後方粒跡線解析

34. 日本に飛来する黄砂 2003年4月16日 （JAXA；EORC)

35. 本学及び父島における2003年3月の黄砂飛来時の体積粒径分布本学及び父島における2003年3月の黄砂飛来時の体積粒径分布

36. シベリア森林火災による、煙粒子飛来時（2003年5月22日及26日）におけるNOAA HYSPLIT MODELによる後方粒跡線解析

37. 2003/05/19 at 03 :45 UTC,Fires and smoke in Russia and China Courtesy of MODIS Rapid Response Project at NASA/GSFC

38. 日本に飛来する煙粒子　2003年10月17日　（JAXA：EORC)日本に飛来する煙粒子　2003年10月17日　（JAXA：EORC)

39. 本学及び父島における2003年5月（本学は及び6月）のシベリア森林火災による煙粒子飛来時の体積粒径分布本学及び父島における2003年5月（本学は及び6月）のシベリア森林火災による煙粒子飛来時の体積粒径分布

40. 本学及び父島における、黄砂時と森林火災時のa*とgの相関本学及び父島における、黄砂時と森林火災時のa*とgの相関

41. 本学及び父島における、黄砂時と森林火災時の体積粒径分布の比較本学及び父島における、黄砂時と森林火災時の体積粒径分布の比較

42. まとめ 今回の解析で、以下のことが認められた。 ・本学（東京）では、春から夏にかけて光学的厚さが高くなり、 　　Ångstrom指数は比較的低く、逆相関を示す。一方、秋から冬　　にかけては、光学的厚さは小さくなり、 Ångstrom指数は高い。　特に冬において、光学的厚さとÅngstrom指数は正相関を示す。 ・単一散乱アルベドは、夏に大きく、冬に小さい。 ・父島においては、春と秋に光学的厚さが高くなり、 夏は小さい。 ・これらの変化は、東アジア・西部太平洋域の気候的季節変化に対応しており、　特に、春と秋には黄砂などの自然起源エアロゾルに加え、中国　大陸からの汚染大気の流入も認められた。

43. ・2003年春に発生した、黄砂及び森林火災に関して、本学及び　　父島の両観測点で捉えられ、その広範囲にわたる影響と輸送　が認められた。・2003年春に発生した、黄砂及び森林火災に関して、本学及び　　父島の両観測点で捉えられ、その広範囲にわたる影響と輸送　が認められた。 ・スカイラジオメーターの観測から、日本本州上空を通過した後、　2~3日後に父島に達したことが窺え、後方粒跡線解析の結果　　からも、これが認められた。 ・シベリア森林火災に関して、父島では本学ほど微小粒子の卓　　越が認められず、輸送過程におけるエアロゾルの拡散、沈着　　や海塩粒子との凝集による粒子の成長も予想される。 ・ Ångstromの経験式の2次式拡張に関して、a*、gはa 同様にエ　アロゾルの粒径分布に関して十分に議論でき、より多くの情報　を与えうることが解析により示唆できた。

44. 今後の課題 ・本研究によって、グローバルな季節的気圧配置、気流の変化に　対応した、エアロゾルの光学的特性や体積粒径分布の変化は　捕らえられたものの、その中において、観測点近傍からの局所　的な原因によるものとの判断は困難であった。地表から大気上　端までの気柱積算量だけではなく、水平、鉛直方向に対する情　報など、より厳密な観測、解析が必要である。 ・ Ångstromの経験式の2次式拡張に関して、 a*、g の粒径分布　における意味は、その定義と解析によって示唆されたが、数値　的、物理的な裏づけと、厳密な議論に関して、Mie散乱理論に　　基づく数値シミュレーションによる検証を行う。