２：成層圏における物質輸送

1. ２：成層圏における物質輸送 北半球夏 成層圏における物質輸送では、オゾンにからむ輸送がおおきな問題であろうから、話初めとして、オゾンの分布図から 北半球冬 =5mmの厚さ(標準状態) オゾンの鉛直分布（成層圏で大きな値）、WMO-O3 report, 2007から 全オゾンの全球分布：全オゾンでみると、中高緯度の方が多い、また冬に多いこと

2. 極での全オゾンは春が最大になる。運動が大事（主に、惑星波動によって輸送される）極での全オゾンは春が最大になる。運動が大事（主に、惑星波動によって輸送される） 東西平均した緯度高度図： オゾンホール（南半球） Ozone mixing ratioの南北−高度分布（１月、４月、７月、１０月：ここで１０月は最近の南極域オゾンホールで異なる）、ppmv ある高度でオゾンが無くなっている ＜ーその場の化学過程であろう

3. ２−１：成層圏における運動の概観 冬 夏 対流圏との運動の違いは、成層圏では水が少なく大気が基本的に安定大気であることであろう (N2= g d (lnθ )/dz が大きい ）-> 波動が卓越しやすいであろう。 西風 東風 より冷たい 西風 １月（左）と７月の東西平均した東西風、冬と夏で構造が異なる 全球平均の鉛直温度構造

4. ５０ｋｍ〜９０ｋｍは中間圏と呼ばれ、温度は高さとともに減少しているが、水はほとんどなく，温度勾配も緩やかなので，基本的には対流は起こっていない５０ｋｍ〜９０ｋｍは中間圏と呼ばれ、温度は高さとともに減少しているが、水はほとんどなく，温度勾配も緩やかなので，基本的には対流は起こっていない しかし、対流圏から重力波が伝わってきて，しかもここで振幅が大きくなり，局所的に不安定が起こっている ー＞　重力波が壊れている中間圏の重力波に伴う温度構造： また、中間圏の重力波(G)が南北循環を生み出している（Plumb, 2002, J. M.S. Japan） 外部重力波のbreakingです（北斎から） こんなimageか？ ー＞鉛直拡散係数として 局所的に乾燥断熱減率になって、乾燥対流の起きる条件をみたすことがあるー＞物質の上下拡散に重要

5. 成層圏では、惑星波動による輸送が重要である成層圏では、惑星波動による輸送が重要である 成層圏中の東西に非一様な大気波動： 冬季惑星規模の波動 東西平均温度の1979年時間変化、大きく変動している 東風 西風 図：１９７９年１月２６日、１０hPaのHeight図 (m) 。北極からみた図である。冬季であることに注意、夏は惑星波動がみえない １９７９年２月２６日の東西平均風 このような成層圏の中の大気運動に絡めて、物質輸送の問題を考える。 ー＞上の例では、東西に一様な風の上に比較的簡単な構造がのっているので、 東西平均の場に、波としての擾乱を考えて議論することが1960年ころから行われたよう

6. 惑星波動を記述する方程式： 運動方程式に現れるコリオリ項を のように βー平面近似した準地衡風方程式でのPotential Vorticity（渦位）方程式 ここでug, vgは地衡風を表し、 等温位面でのポテンシャル・渦度 qは準地衡風でのPVをあらわす。 ψは流線関数であり、 散逸、非断熱がない場合は、PVが保存される 散逸や非断熱がない時、時間的に一つの(Rossby)モードの、保存的な時間発展の式 この方程式は中・高緯度の対流圏で重要な役割をもっている傾圧不安定、いろいろな惑星波動の問題に適用される。１０００ｋｍくらいより大きいスケ−ルの運動にたいしての式であろう． PVの解析例はあとに、

7. Andrews and McIntyre, 1978, J. F. M. 南北および鉛直方向 東方向 変位 波動にともなう、流体粒子の変位を　　　　　　　とし 図を参考にして、子午面内のラグランジュ平均の流れは、 ２−２：物質輸送の１つの表し方 波が定常であれば、 １項は場所を固定したオイラー平均、２項はStokes Driftと呼ばれる。変位と速度は線形的に、 中緯度の惑星波動による輸送について、東西平均量(over bar)とそれからのずれの波動(’)に分離し、波動に伴う流体粒子の動きを考慮した、ラグランジュ的な平均を考える。 のような形になるので、 で、連続の式　　　　　　　　　　を使うと また のようになる。

8. 波が定常のときは、 東西平均した物質の変化の式は、平均子午面循環による流れと波動による輸送によるとして で、物質の変化の式は 波動による物質輸送は、保存的とすれば のようなので、対称成分も考えて、以下のような近似式が得られるであろう。 ２−２：物質輸送の１つの表し方（続） の形にかけるので、波動による輸送の項は 流体粒子の波動にともなう変位は観測が不可なので、惑星波動の熱力学の式を変形して、観測しやすいオイラー的量の表現にすると、 となる。 テンソルを対称部分と反対称部分に分けると下のようになる(Holton, 1981, JGR)。 対応した、東西平均の熱力学の式は、 前ページの、　　　　　　　　　　から は残差子午面循環と呼ばれる 運動量の式は、 はEliassen-Palm fluxと呼ばれる

9. 循環を評価した１つの例 残差子午面循環 オイラー平均の子午面循環 90N 冬半球 eq. eq. 北半球冬 成層圏下層に熱帯から上昇している循環がみえる（Brewer-Dobson循環とよばれる）、上層の方はHeatingのある夏半球が上昇流で冬半球が下降流となっている。

10. 流れ場を与えて、流体粒子をLagrange的に動かしてー＞南北、鉛直の２次元に投影してみる, Kida, 1983, JMSJ 拡散について： (A) 式で示したように、移流と拡散（波の振幅の変動で起こりうる）が起こっている。 核実験後の成層圏の物質分布が上のKidaの数値実験の結果と良く似ている (A) のXXXは初期の場所、(A)は６ヶ月後の分布図、あと１年後づつ

11. Brasseur et al. (1990)では、緯度—鉛直の方向の輸送として、下のような式が輸送の式として用いられている 重力波による鉛直拡散とRossby波による水平拡散 残差循環　　　　　　により物質が移流され、さらに南北、鉛直拡散が考慮される形 中間圏の鉛直拡散は重力波のbreakingによる運動量の拡散係数が使われている。中間圏で大きな拡散になっている 南北拡散は準地衡風近似のとき、ポテンシャル渦度の南北輸送がEP-flux 発散に等しいことで(Andrews et al., 1987)、 のように書かれ、 のような拡散の形として のように定義されている。

12. 化学反応のもとに、運動の方程式と組み合わせて、物質の輸送の問題が鉛直緯度２次元の範囲で解けることになる化学反応のもとに、運動の方程式と組み合わせて、物質の輸送の問題が鉛直緯度２次元の範囲で解けることになる 計算されたオゾン分布 ー＞　近年は、１章で説明したような、対流圏を含めた、全球３次元的運動と化学過程を直接解く３次元モデルによる研究に移行

13. ２−３：化学気候モデル モデル結果 例えば、NCARの大気大循環モデルをベースにして、成層圏オゾンを主体とした化学過程をモデルに導入して、化学気候モデルが作成される。 観測結果 Rasch et al., 1995, JGR 大雑把には再現、しかし細かいプロセスまでは表現できていないよう。 化学反応の例

14. ３DモデルによるCH4分布 Rasch et al., JGR, 1995 ３次元モデルであるNCAR GCMで再現されたCH4分布図： 観測結果 0.4 0.4 ３次元モデルで、赤道域では鉛直上昇流、中高緯度では下降流（前に述べた循環に対応）および惑星波動による水平混合など、平均的な形はおおよそ再現されている。 ３０度あたりの、低緯度と中・高緯度の境界あたりの構造などがなめらか、

15. N2OとCH4の寿命の長い物質の水平分布図：中／高緯度では成層圏の惑星波動にしたがって運動しているよう。N2OとCH4の寿命の長い物質の水平分布図：中／高緯度では成層圏の惑星波動にしたがって運動しているよう。 これまでの話しをまとめると、成層圏の惑星波動に伴う輸送は、大雑把には下図のようなimageでしょうか Roche et al., 1996, J. Geophy. Res. から、21mb(27km)の高度、CLAES ( cryogenic(低温) limb array etalon spectrometer)衛星観測 Solomon, 1999, Rev. Geophys.

16. ２−４：成層圏物質輸送のプロセスをいくつか２−４：成層圏物質輸送のプロセスをいくつか 2005年6月、100hPaにおけるCO分布、Park et al., JGR, 2009 7月、気候値の100hPaでの水蒸気分布、Gettleman et al., JGR, 2004 場所に依存していて、地域的な分布が異なる。 ー＞個々のプロセスを詳細にしらべる。

17. 1：対流圏からの流入ーインドネシア域の熱帯圏界面近傍の物質輸送ー1：対流圏からの流入ーインドネシア域の熱帯圏界面近傍の物質輸送ー Lagrange的な方法による、Hatsushika and Yamazaki, 2003, JGR 温位350K以下からの粒子から初めて、390K以上にいった粒子の軌跡。赤系統から青系統で水蒸気の量は減っていく。La Ninaの１月の条件で動かしている。全球に広がるにはtransient運動が重要とある。 対流圏から成層圏への粒子の移動の概念図 ＜ー熱帯の圏界面あたりの水平移流の重要性

18. 西風偏差で高圧偏差の構造 ２：熱帯域での成層圏から対流圏への物質移動の例 赤道ケルビン波は赤道成層圏において、いろいろの役割をしている。ケルビン波の構造のみを示しておく。この波は南北風がほとんどない波である。 東西に伝播する波として　のように仮定すると 赤道ケルビン波の水平構造 のような、南北(y=0が赤道）にガウス分布となる構造をもっている。 鉛直波長２π／m 鉛直方向には、N2＝一定大気の場合には、 下降変位 ζ＜０ のような伝播性をもち、近似的に下式の関係になる。 東西波長２π／k 東西風に相対的な振動数が小さくなると、鉛直波長は短くなる。 ケルビン波の東西鉛直構造、小さな矢は　

19. Fujiwara and Takahashi, 2001, JGR 熱帯域におけるケルビン波による輸送 34E, 8月15日、オゾンの緯度高度図、Kelvin波のガウス的南北の広がりがみえる。 Kelvin波の下方変位に伴って成層圏から対流圏の方にオゾンが輸送されているよう。 オゾンと東西風（影が東向き）の経度／高度断面図　

20. 二重の圏界面になっている Pan et al., 2009 JGRを例として ３：中緯度圏界面付近の輸送の複雑さ 200７年4月11日の色は最小のdΘ/dz、赤実線は6PVunit線、赤点線は衛星track オゾンの緯度高度図 150hPa 60N 60N 40N 30N Potential Vorticity 90W 10日のbackward trajectory, 緑60-67N, 赤50-57N,13-16km高度 1PV unit = 10-6 m2s-1K kg-1

21. ４：2002年突然昇温によるオゾン変動 10月平均オゾンホールの変化 0° 1979 1982 1985 1988 http://ozonewatch.gsfc.nasa.gov/monthly/index.html 1991 1994 1997 2000 2002 2003 2006 2009

22. 南半球（オゾンホール）の様子（1998-2003年、9月25日のみ）、全オゾン南半球（オゾンホール）の様子（1998-2003年、9月25日のみ）、全オゾン 1998 1999 2000 2002 2001 2003 基本の構造は南極で少なく、オーストラリアの南の方で多いという波数１的パターンが多い。オゾンホールの形は年によりすこしづつ異なっている。2002年はかなり形態が異なっている　ー＞　この年はmajor な突然昇温現象

23. 廣岡、森、他 (2004) から 2002年は突然昇温が起きたこと 突然昇温の現象；図は南緯６０度、５０hPaの２００２年４月から１０月までの東西に平均した温度の時間変化を示したもの。冬から春への温度変化のなかで（低温からだんだん温度が上がりつつあるとき）、たまに急に温度が上がっている。この様な突然の温度増加現象を成層圏突然昇温と呼んでいる。また極の高温は温度風の関係から東風になる可能性があるので（夏の状況）、１０ｍｂ以下で６０度から極向きに温度が増加して東風が出来るとそれをmajorwarming と呼んでいる。 惑星波動の全球東西波数s＝1の振幅 波数s＝２ MAY JUN JUL AUG SEP OCT

24. 南半球の10hPa等圧面高度図（約30kmの高度）。単位はm、等値線間隔は200mの高さの違い。南半球の10hPa等圧面高度図（約30kmの高度）。単位はm、等値線間隔は200mの高さの違い。 オゾン全量 2002年オゾンホールの急激な変動（9月19-29日） 9月19日 9月25日 波数１ー＞２ 9月23日 9月29日 オゾン偏差χ’と高度偏差の近似的な関係

25. 180° Thompson and Solomon, 2002, Science 76S, 27Wでの オゾン量 補足：オゾンホールの気候変化 0° 12-5月、地表の温度（色）と950hPaでの風のトレンド、南極で冷たく、周りの風（矢）が強化されている 11月のオゾン変化：年々で変動しつつ、減少しつつある 圧力低下 約15km 温度低下 11月 11月 温度と圧力の低下傾向をしめす

26. Austin and Wilson, 2006, JGR 温暖化（成層圏では低温化）に関わって、北極域のオゾン回復の戻りが大きいという話しも出てきている 将来どのようになっていくかの、モデルによる、オゾンホールの予測実験が行われている AMTRAC 気候センター／環境研モデル 補足：オゾンホールの予測 成層圏の低温化により、オゾン生成の効率化が上がりオゾンが増えるであろうが、南極の方ではこのように極端なことはおこっておらず、その差は北半球は循環がはやく少ないフロンが極にくるためと考えられている

27. 上部成層圏の水蒸気変化 ５：成層圏の水蒸気の気候変化Rosenrof et al., GRL, 2001 HALOEによるH2O分布例 より高い高度での変化、8mb（上）, 2mb（下）、月平均の変動 水蒸気の変化、21.5mb, Boulder, コロラド(40N, 105W)のバルーン観測　HALOEは35-45N, 95-115W　衛星観測 論文では、対流圏から成層圏への輸送の変化が一つの可能性とある

28. 気候モデルによる水蒸気変化：Austin et al., JAS, 2007 化学気候モデルを1960-2005まで走らせる。SST変化などは観測データを利用、また温暖化物質なども増加させるー＞水蒸気の気候変化をしらべる 成層圏の水蒸気が増加、要因としては、CH4の酸化と書かれてある。 のような反応（cf., Le Texier et al., 1988, QJRMS） CH4の増加は循環の強化のよう 水蒸気の時間変化、青はモデル結果、赤は衛星観測（HALOEによる）の結果、BlackはBoulderでの観測結果である。 77hPaでの上方質量輸送の変化