10章

第10章:成層圏突然昇温など
--惑星波動による平均東西風変化の例-10−1: EP-フラックスによる解析
Andrews McIntyre(1976)によって導入された、変換された(transformed)オイラ—平均で惑星波動の振る舞いをみることが
よくなされる。
準地衡風系では(cf. Edmon et al., 1980, JAS)
r
u
1
 f 0v * 
F  X
t

*
v
1 

w *  0
y
 z
T
H

N 2 w*  J / c p
t
R
熱力学の式で、擾乱の効果が見えないこと(非断熱が鉛直循環を直接駆動する形)、運動の方程式において東西平均風

の加速の項が、波の2次量であるEP
flux (Eliassen-Palm flux)の発散によって表現される(惑星波動ではu’w’とか dU/dz
等の項は落ちる)。
Fy    u' v'
Fz   f0
R
v' T' / N 2
H
この量は惑星波に伴う東風運動量を南北、鉛直に運ぶ指標で、psudo-運動量フラックスとも呼ばれる。
運動量フラックスの発散によって、風(運動量)が直接変化することを示していて、物理的に理解しやすい表現になってい
る。また残差子午面内循環は、近似的に重心の平均的な南北、鉛直の運動状態を記述していると考えられる。
v*  v 
w
*
1 R 
(  v' T' / N 2 )
 H z
 w
R 
2
(v' T' / N )
H y
惑星波の
Eliassen-Palm
フラックスを図
にのせておこう。
有効な手法で、
解析でよく使わ
れている.
ー>波の振幅
と位相の表現
から、波の運動
量の流れという
考え方か
図:定常
planetary wave
のEliassenPalm flux。196
3年から1969
年までの1月で
波数1である。
Sato(1980, J. M.
S. J. )より、年
により非常に異
なる
10−2:成層圏突然昇温について
西風
東風
冬の成層圏の基本的な平均東西風は西風である
ー>東風に変わるときがある
西風中 U(y,z) の定常惑星波動の伝播の様子、モ
デル中での波のエネルギの流れ
(  'v'  U u'v' ,  ' w'  U
f2
v'  z ')
N2
突然昇温の現象について例示しよう。図は北緯80度、10hPaの1978年10月から1979年5月までの東西
に平均した温度の時間変化を示したものである。冬から春への温度変化のなかで(低温からだんだん温度が
上がりつつあるとき)、時々急に温度が上がっている。この年は3度起こっている。この様な突然の温度増加
現象を成層圏突然昇温と呼んでいる。英語ではstratospheric sudden warming 。また極の高温は温度風の関
係から東風になる可能性があるので(夏の状況)、10mb以下で60度から極向きに温度が増加して東風が
出来るとそれを major stratospheric warming 、と呼んでいる。かなり不規則で(北半球では2年に一回程度)、
どの年にmajor warmingが起こるかわかっていない。対流圏の年々の状況にもよるであろうし(惑星波の生成
問題?)、また赤道下部成層圏の準2年振動と関係があるともいわれてはいるが(これは波の伝播問題と関
わるであろう?)
1978
1979
表:majorな突然昇温の起こった年、Andrews et al.(1987)
より
1979年突然昇温の平均東西風時間変化、Andrews et al.(1987)の
教科書
12月8日/78年
1月25日
2月26日
2月6日
前ページ図に対応したときの平均東西
風変化の様子、それぞれ12月8日、1
月25日、2月6日、2月26日、3月3
日である。12月8日は冬のはじめで
西風が強い。1月25日および2月6日
は温度が上がっており、それにともな
い極域に東風が吹いているが10mb
では東風になっていないので minorと
している。また2月26日には10mbで
東風になっているのでこれはmajor
warmingとなっている。
3月3日
s=1の振幅
s=2
時間ー緯度
断面図:
2/6
2/26
1/25
1/25
2月
2/26
major
warming
20N
warm
cold
40N
180
major warmingに対応した、10mbでの Planetary
wave の振舞いを図に示す。日にちはそれぞれ2
月17日、2月19日、2月21日、2月26日、3月1
日、3月5日の温度(5度おき、dashed curve)とハ
イト(0.2kmおき)を示す。はじめ気圧場の水平構
造は極渦が引き延ばされて、楕円のような構造に
なっていて、渦の中心が少しpoleから離れている。
それにアリューシャン高気圧が付随している。(b)
では(a)のような定常・惑星波の構造が少し変形
しつつある。ー>次へ
図:1979年の突然昇温のときの10mbの温度
と高度の分布。Andrews et al.(1987)の教科書よ
り。
0
a : 2月17日
b : 2月19日
c : 2月21日/
1979年
warm
(c)で大きな変
化が起こって
いる。低気圧
の渦が2つに
分離されたよ
うな形になり、
極が高温にな
りつつある。
(d)では極が
高温になり、ま
た極が高気圧
になっている。
極の高気圧に
ともない東風
が吹く。そして
しばらく時間
(数日、放射の
緩和時間)が
たった後また
冬の状態(完
全ではない
が)に戻る。
180
0
240
cold
e : 3月1日
warm
high
f : 3月5日
d : 2月26日
このような現象を、Matsuno (1971)によるPlanetary wave の鉛直伝播と、その波と平均東西流(および平均温度場)との
相互作用の観点から見てみよう。概略を述べると以下のようになるであろうか。あるとき対流圏においてPlanetary wave
が増幅される。この増幅の機構は対流圏の現象と関係があるらしいがまだ明確になっていない。ともかく惑星波が強まっ
てその波が鉛直へ伝播していく。上方に伝播し波の振幅は密度factorによりさらに強められる。そのときtransientな波の
非線形効果により平均流を変化させる。波が定常であればEliassen-Palmの定理により何の変化ももたらさない。しかしい
まは波が急に増幅したので、上の定理は破綻して、東西平均流は変化していく。
数値実験による説明:波動に関する時間発展の式(形は線形)

  1  cos '
1
2 '
p ' 
2 2 
(   )
(
)
 4 a
(
)
t
 cos   sin 2  
cos2  sin 2  2
pz N 2 z 


q 1 '
0
 cos  
 
u
a cos 
Zonal mean equation:(東西平均場が時間的に変化する式)

u
1

 2sin  v  
(u' v' cos 2  )
2
t
a cos  
 
1
 '
2
(
)N w  
(
v' cos  )
t z
a cos   z


( pv cos  ) 
( pw cos  )  0
a
z
のような式をcoupleして解いてある。擾乱が東西平均場を変え(下の方の式)、変わった平均場を擾乱が感じて(上の方の
式)… のように発展していく。
数値実験:
モデル振幅
高
度
点線は観測
西風
下部境界での惑星波動の振幅変動、t=0から波を
強制する
波数1の振る舞い
ゼロwind-lineが
下方に伝播
東風
緯度
初期 t=0(初期条件)における平均東西風
計算された波動振幅の時間変化、β平面モデルで初
期は一定の風(33m/s)の場合、shadeの部分が東風、
30度と90度に壁で、60度での様子
時間変動の様子:
高
度
時間
60Nの平均東西風の時間変化、西風であったところ
から東風が生成されている。波数1の強制
波数2の場合の波の振幅の時間変化 -> 下図
と対応
温度下降
高
度
初期条件からの極の温度の時間変化、成
層圏は温度が上昇、中間圏は温度下降して
いる、波数1の場合。
平均東西風の時間変化、波数2の場合は低
い高度で東風に
水平の構造:
t=0で波を
forcing
10日後
緯
度
平均東西風の時間変化、波数2の場合
緯
度
西風
東風
水平パターンの時間的変化の様子(30km)、波数2の場
合、極の低気圧が高気圧に変わっている。
波の振幅の時間変化、線形の為か
critical level が動き過ぎのよう
上方伝播に限ったオイラー平均的説明:惑星波に伴って、熱輸送がある。北側で上昇流が作られ、
Eliassen-Palmの定理の破れのために連続の式から南北風は北風、それにコリオリが働いて、東
風をつくる。
u
 fv  0
t

 (  )  N2 w    v   
t z
y z
波
の
鉛
直
伝
播
熱輸送の効果の方が
勝って温度上昇
北側の温度変化
対応した東西風変化
成層圏warmingのLagrange平均的な見方について:
成層圏warmingを惑星波動が鉛直に伝播して平均東西風と相互作用をすることで説明した。それを波にともなう流体粒子
の変位に伴った平均をするLagrangian 平均の立場で説明する。
Andrews and McIntyre(1978)から引用する:
r r

( x,t)
とすると、Lagrange
的平均は
上図のような波の変位を考慮した平均をする。この粒子変位を
r
r
r r
 ( x, t)   ( x  ( x, t), t)
L

のように定義される。
定常の惑星波が鉛直に伝わり、平均を東西方向にえらべば、Lagrange平均は

 (y,z)    (x  , y  ,z   )dx /  dx
L
のようになり(Matsuno and Nakamura, 1979, J. Atmos. Sci.)、平均の南北風はTaylor展開することで
v
L
w
L
 w
r
  v
r
   w
 v 
のようになる。右辺1項はこれまでのEuler的平均を意味する。右辺2項をStokes Driftとよぶ。

いまの場合、変位は
u
v


 ez / H
(e z / H w)  0
x
y
z
を時間積分すれば
  


 ez / H
(e z / H )  0
x
y
z
<いいのであろう

なる連続の式をみたす。Lagrange的子午面循環は実質的な流体の重心の平均運動と考えられていて物質粒子の運動に
適用される。これは近似的に、* のついたresidual circulationである(AM76参照)
惑星波が臨界層(今の場合はU=0に対応)に伝播しつつあるとき(線形の定常波のとき波は吸収される)のLagrange的な
子午面循環は以下の図のようになる。Critical Levelでは北向きとなっている。
東
風
定常惑星波動
にとっての臨
界層(U=0)
U=0
西
風
下
降
流
南
北
対応した平均東西風と温度の時間的変化
オイラー平均(場所に固定した東西平均)における子午
面循環(流線関数と鉛直流)
極側
今の場合、
T
H

N 2w L  0
t
R

昇温
EP フラックスのLagrange的解釈:鉛直伝播する定常な惑星波動では(Matsuno and Nakamura, 1979)
u
1 
f R
 fv * 
( 2
v'T')
t
 z
N H
1 
1
1 
1


(
p' w') 
(
p' u0 (z)
)
 z u0 (z)
 z u0 (z)
x
1 


p'
 z
x
西向き
波はここま
できていな
い  0
東向き+
のような形になる。

最後の項は、惑星波の鉛直変位
しめす。
にともなう応力を
西風
平均東西風として、西風が吹いている。図のA点で、矢

羽はEの方(西風)をむいている。下のほうから定常惑
星波が伝播している。波にともなって流体粒子面は凸
凹している(図のB点に対応しており、定常惑星波動に
ともなう鉛直変位のx微分が+のとき圧力偏差は+に
なっているので、その積は+となる)。その鉛直微分は
Aでは波がなく、Bで波が伝播しているとすればz-微
分は - (負)となり、力として -加速(東風加速、西風を
減速)のようになっている。
南風
p'

0
x


ただ、変位は普通は、解析できないので、EP-flux
Divergence、残差循環で評価される。

0
x


0
x
東風加速
波数1が主
EP-フラックスによる解析:図は前
述の1979年major warmingのと
きのもので矢印は前に述べた
Eliassen-Palmのフラックス。
加速は
 u  fv*  X  (  a cos  )  1  F  D
0
F
t
簡単には右辺のEliassen-Palm
flux の発散が平均東西風を変
化させると思う。図には収束に
よる加速 ものっている。時間的
に非常に複雑な変化を示してい
る。21日あたりは波が収束的
で東風加速になっている。一方、
28日ではEP-fluxは発散になっ
ており、西風加速になっている。
太い矢羽根は
波2の寄与
2002年の南半球は特別の年:南半球で観測史上
初めてのMajor Warmingがおこった
01年(は平年的)と比較する
2002年
1/0ct15/Oct
太い実線は1979-2002平均の東西風の季節変化(影は標準
偏差)、Hio and Yoden, 2005, JASから
東西風
2002年
16/Aug-30/Sep
鉛直EP-fluxと平均東西風の関係の各年
のscatter plot、2002年は鉛直EP-fluxが
強く、平均東西風が弱い
細線が2002年、鉛直EP-fluxが10月に大きな値をもっ
ている
近年の南半球(オゾンホール)の様子(1998-2003年、 9月25日のみ)、全オゾン
1998
2001
1999
2002
基本の構造は南極で少なく、オーストラリアの南の方で多いという波数1的パターンが多い。
オゾンホールの形は年によりすこしづつ異なっている。
2002年はかなり形態が異なっている ー> この年に major warming
2000
2003
2002年オゾンホールの急激な変動(9月19-29日)
オゾン全量
9月19日
9月23日
9月25日
左図に対応した、南半球の10hPa等圧面高度図(約
30kmの高度) 。単位はm、等値線間隔は200mの高さ
の違い。
図は廣岡、森、他 (2004) から
9月29日
波数1から、2が卓越している
補足:大気変動の仕方は年により違う
2002年と2001年の南緯60度における東西に平均した東西風の時間と高度(縦軸)の図、2002年
で9月の終わり頃(波の形態がものすごく変形した時)に東風が下方に伝わっている様子をしめす。
それに比べ、2001年は比較的ゆっくりした季節変動をしている
高
度
月
2002年(左図)9月の終わり頃(波の形態がものすごく変形した時)、極でオゾンが増大している時期ー>極の方が温度が
高温になり、西風が東風に変わっている(Majorの突然昇温になっている)様子で、10hPa(約30kmの高度)における東西
に平均した温度、東西風の時間変化(5月 -10月)と、惑星波動の振幅変動を示す。
一方、右図は、2001年での同じ図であるが、比較的ゆっくりした季節変動
緯
度
東西風
が東風
にならず
に、西風
が長い
期間吹
いている。
波動の強
さは2002
年ほどに
は強くな
い
波数k=1振幅
波数2の
振幅は小
さい
波数k=2
MAY
JUN
2002
JUL
AUG
SEP
OCT
2001
2002年と2001年のEP flux(矢羽根の長さ)が異なるー>大きな変動をおこす
2002年
2001年
波数k=1
波数k=2
波数k=3
2002年は波の活動が強く、成層圏の中にまで侵入している様子がみえる、上から東西波数が小さいもの
(k=1, 2, 3)から並べている(南緯50−70度平均)、色は波のEP-flux 発散の大きさをしめす、矢の右向き
は極向きを意味する。
対流圏の様子:予測実験との違いから
実況
9月19日−21日の対
流圏の様子
k=2のEP-flux
傾圧波動( k=4-6 )によるEP-flux、
大きな活動度および減速がある
予測を外している例
東西風が異なっている
9月13日を初期値にした
ときのモデル結果の対流
圏のパターン
9月13日を初期値にしたときの
モデル結果のEP-flux
10−3:Downward controlについて
Haynes et al., JAS, 1991
成層圏の中である外力が働いた時の大気の応答についての議論で、時間がt=∞たった時を考える
1

1 
(v * cos  ) 
( w * )  0
a cos  
 z
のような残差循環についての連続の式が成り立つから、流線関数を導入する。

1

1

v*  
w* 
 cos  z
a cos  
一方、定常状態での運動方程式は
 

1

u
v * 
(u cos  )  2sin  w *
 F
z
a cos  

角運動量mを導入する
m  acos  (u  acos  )

上式の運動方程式は

 (, m )
 m
 m


 a 2 F cos 2 
 (, z)
 z
z 
m
 (m , z)
divided
by


 (, z)
 (, m )
 (, z)  (  )  a 2 F cos 2 
 (m , z)
z m
m /
 (, z)

Local なForcingを与えた数値実験の例
境界条件は上端の無限で、
u
流れ関数
w *   0 z  
   0 z  
a 2 F cos2 
 (, z)   z {
}dz
m /
2
2
 a F cos 
1

*
w 
 { m / }dz
 cos   z

u
流れ関数
ある高度の鉛直流はそれより上に
あるFの分布だけできまる。外力
は下方のみに影響を及ぼす。
極限の式では
中緯度forcingにより、赤道では上昇流
(a), (b)放射減衰がないときのuと流れ関数の応答、(c),(d)は4日放
射減衰で10日後、(e),(f)は20日、30日後の流れ関数
熱帯圏界面の温度変化について:
赤道上の温度の1年振動:北半球の冬に低温化、夏に高温化、
ECMWF ERA-15 data、1979-2001の平均
Kerr-Munslow and Norton, JAS, 2006, pp1410-1419
90hPa, 10N-10Sでの温度変化、実線が温位変化、
dashが鉛直移流の効果、dotは水平移流、dot-dash
は非断熱、triple dot-dashがその他
Q


t 
wave heating


残差鉛直移流によって温度が変動、
w * N2
w * mm/s
V*の季節変化について:

net  EPD
v'T'  term

夏はw*が弱くなる、南半球冬がdownward
controlが弱いと思っていいか?


(u')v'
y

((u')w')
z

EPDがV*の変化に対して主、その年変動に対して、3点
dot-dashが対応:u’w’項である。

90hPa, 10S-10N平均
赤道Rossby波による
ー>downward controlではないと言っている
10−4:成層圏のArctic
Oscillationについて
時間的変動、下方伝播のように見える、赤がweak,
warm vortexである、赤□はmajorまたはearly final昇温、
Cはカナダwarming、90day low-pass filterがかかってい
る
AOの高度別パターン
東西風の南北dipole構造
下方伝播の様子
極で低圧偏差のとき、中緯
度では高圧パターン
赤のweak vortexは突然昇温の弱い風に対応している
と思っていいであろう、
Baldwin and Dunkerton, 1999か
ら
10-dayのlow-pass filterの場合
Itoh and Harada, 2008, JCでは、第1 mode は
同時的
一方、第2モード(成層圏ではs=1の惑星波
動)は対流圏から成層圏への上方伝播的
に見える
対流圏での先行が目立つように感じる
補足:AOへの水平運動量輸送の役割
Limpasuvan and Hartmann, 2000, JC
対流圏西風が極にshift 波が赤道に屈折、極にeddy momentum fluxを輸送し、高緯度に西風をさらに強化、
北半球Annular Modeは惑星波動が卓越、一方南半球Annular Modeについてはtransient波が卓越と言ってい
る。
high
順圧 forcing difference
low
high low
全forcing difference、矢はEP-flux
difference
jetの軸は55N
北
北半球の k=1屈折率差, 差
として負(外部波的)が6070度にある
別GCM中のAOとEP-flux ( Kornich et al., GRL, 2003 )
GCM中のAOの高度別パターン
300hPa(上)と10hPa(下)におけるAO
のパターン
東西平均した東西風anomalyとEP-flux anomaly
簡単モデルにおけるAO (Eichelberger and Holton, 2002,
JGR
EP-flux とその発散、t =590 〜 650までの間の時間変動、
t=605で大きな減速、t = 635で変形してる:sin2lyを含める
ことでdiple 構造を表現
東西風および惑星波動を、南北には2つのモードを取
り出して議論するモデル
15m/s間隔
類似のパ
ターンが得
られている。
5m/s間隔
東西風の sinly 第1モードanomalyとsin2ly 第2モード(南北
反対)の時間変動、影は西風anomalyに対応
ー>AOと似たような振動が再現されている
対応した東西風の変動、10m/s間隔