Improvements of the Eastward Propagation of the MJO in MIROC6

Presentation on theme: "Improvements of the Eastward Propagation of the MJO in MIROC6"— Presentation transcript:

1 Improvements of the Eastward Propagation of the MJO in MIROC6
Nagio Hirota1, T. Ogura1, H. Tatebe2, H. Shiogama1, M. Watanabe3, M. Kimoto3 NIES JAMSTEC AORI, UTokyo 14 Dec. 2016 Hirota, N., T. Ogura, H. Tatebe, H. Shiogama, M. Watanabe, M. Kimoto, Improvements of the Eastward Propagation of the MJO in MIROC6, in preparation. 20min presentation / 10min discussion T. Kataoka？

2 (Madden & Julian, 1972) Introduction Madden Julian Oscillation: The dominant intraseasonal variations over the tropics with period of 20-90days (Madden & Julian 1971) Large impacts on ENSO, monsoons, tropical cyclones, and … (e.g. Lau & Waliser 2011) ➡ Important for Seasonal & future prediction Africa Indonesia America

3 No CMIP3 models simulate realistic MJO (Lin et al. 2006)
MJO in climate models No CMIP3 models simulate realistic MJO (Lin et al. 2006) Some CMIP5 models, including MIROC5, reproduce the MJO spectral peak qualitatively, but its amplitude is underestimated & the eastward propagation over the Pacific is poorly simulated (e. g. Kim et al. 2009) Lag regression wrt Prcp at (0°,85°E) (Hung et al. 2013) Observation　　　 　MRI3　　　 　 MIROC5 Hungの図は30-70日フィルターがかかっている。 MRIもMIROC5もMJO振幅が観測の半分以上はあるいいモデルとして説明されている MRIもMIROC5もインド洋では東進・西進の振幅比2程度で比較的良い。太平洋で東進・西進の振幅比が比較的よいのは、CNRMとFGOALのみ。 MRIはequivalent depthが浅すぎる; MIROC5は比較的よい。などなど

4 Important factors for the MJO representation
Moisture mode: coupling of convection and free trop. moisture (e.g. Sobel & Maloney 2012; 2013) Cold pool dynamics (e.g. Sugiyama 2009) Conv momentum transport (e.g. Miyakawa & Takayabu 2012) Cloud radiative feedback (e.g. Kim et al. 2015) Coupling between convective activities & free trop. moisture (e.g. Sobel & Maloney 2012; 2013) RH composite [%] decay mature dvlp supress Moistening from lower to upper troposphere (Kikuchi & Takayabu, 2004) Kim 2011 Fig. 3 MJO東進メトリック：観測、tokiokaナシ、tokiakaアリはそれぞれ3.5, 0.9, 3.2ぐらい？ day prcpの30S-30N平均：2.3, 2.1,2.6ぐらい？ 観測(MJOmet=3.5)　　 ConvTriggerナシ(0.9)　 ConvTriggerアリ(3.2) May-Oct mean Prcp [mm/day]

5 【Objective】 Improve representation of the MJO amplitude & eastward propagation in MIROC6. 【Data】 NOAA-OLR( ) ERAI( ) 【Model】 MIROC5 (for CMIP5) MIROC6 (for CMIP6) A shallow convective scheme is newly implemented

6 Wavenumber-frequency spectra (Wheeler & Kiladis, 1999)
OLR(15S-15N) AVHRR MIROC5 MIROC6 9.4 (1<ss(pn).and.ss(pn)<3).and. & ! MJO signal Kim et al. (2014) 30<ff(pt).and.ff(pt)<60) 14.7 ← ← ↑ 11.3 East power (Kim et al. 2014)

7 MJO composites OLR(color) & Velocity potential(cont)
(Wheeler & Hendon, 2004) Events = amplitude (√PC12+PC22) > & maximum in a 15-day-window AVHRR/ERAI (43.9%) u200 OLR u850 OLR u850 u200 (20-100day,15S-15N) [W/m2]

8 OLR(color) & Velocity potential (cont; 106m2/s)
MJO propagation OLR(color) & Velocity potential (cont; 106m2/s) AVHRR/ERAI (43.9%) MIROC5 (26%) MIROC6 (32.2%) ← phase ← E 120E W 60W 0 E 120E W 60W 0 E 120E W 60W 0 20oS-20oN [W/m2] Improvements in MIROC6 MJO continues to propagate over the Pacific The explained variance is enhanced

9 Time evolutions of the MJO amplitude (phase4)
amplitude (√PC12+PC22) (life time = days amp > 1) Normalized PC2 Well simulated in MIROC6 Normalized PC1

10 Improved over the Indian & Pacific Ocn. (phase2-4)
AVHRR/ERAI MIROC6 MIROC5 For each phase phase phase phase 3 phase4 phase phase phase 7 phase8 Improved over the Indian & Pacific Ocn. (phase2-4)

11 Specific humidity(color) & ω(contour;Pa/hr)
in the MJO life cycles over the W. Pac. ( E,10S-10N) AVHRR/ERAI AVHRR/JRA55 MIROC5 MIROC6 kfp03 kfp07 → phase → → phase → → phase → [g/kg] WH04b-v15p0-r120 Moistening from lower troposphere is more realistic in MIROC6 Still too much moistening around 600hPa

12 Dry bias in the mean field is mitigated by shallow convection
Specific humidity (contour) & its bias (color) MIROC5 MIROC6 800hPa [g/kg] 800hPa Moistening [g/kg day] by shallow convection 20S-20N [g/kg]

13 Sensitivity experiments
kfp: a factor of TKE used for shallow convection (MIROC6 default value = 0.5) OLR(color) & Velocity potential(cont) MJO amplitude (phase 4) kfp=0.3 (28.8%) less active shallow kfp=0.7 (32.8%) more active shallow [W/m2] Importance of shallow convection is verified

14 Discussion: Coupling of convection & trop. humidity
Prcp. percentile (abscissa) vs RH profile (color) (e.g. Thayer-Calder & Randall 2009; Del Genio et al. 2012; Kim et al. 2014) ERAI MIROC5 MIROC6 pressure(hPa) kfp03 kfp07 (0-360E,20S-20N,Ocn) For the intense prcp (percentile=80~100) MIROC6 shows more realistic RH profile with more moisture at 800hPa compared to MIROC5

15 Improvements of MJO in MIROC6: The amplitude is enhanced
Summary Improvements of MJO in MIROC6: The amplitude is enhanced Eastward propagation & its lifetime (especially over the Pacific) Moistening from lower trop. is more realistic The shallow convective scheme, transporting the boundary layer moisture to the free troposphere more effectively, is responsible for the improvements

16 Structure of MJO Cloud image Convection (shade)
& circulation (vector, cont) Gill (1980)

17 Lag regression wrt Prcp at (0°,85°E) (Hung et al. 2013)
Obs MIROC5 MRI3

18 6/23 SuPAM SPCZの表現は?オーストラリアモンスーンとか？(木本) BSISOは？(建部) フィリピン、MJOの改善と浅い積雲の役割(木本・シンドウ) 3ステージの表現の図とか(木本)？ 東太平洋のMJOも良い？(塩竈、小倉) MJOの寿命は(木本)？ 海が細かいMIROC6も見る?(建部) 冬のモンスーンはWPパターン? ストームトラックは? (森) 7/12 気候感度 MRIも対流圏を湿らせてMJOをだしている(川合) MJOを良くするために雲表現は必要という認識？(佐藤) MJOが必ず良くなったという流れじゃない方がよい？(木本) 7/27 pr/prwの図は話の流れがわるい(木本) 境界層が湿っていれば、カラムでなくてもということ？種別降水は？(塩竈) 段階的な雲の発達を示したい(木本) 言葉遣いが変、ASはエントレインメントない？(木本) 加熱プロファイル見たら？(建部、木本) の加湿が重要といわれている(しんどう) 積雲のパラメータは？(しんどう) 11/28 下層が湿るのは浅い対流？フィードバック？(三井) 浅い対流スキームはどんなもの？(McGregor) 位相速度は?(中野さん) 12/6 位相速度は?(前田さん) 降水持続性の解析は?(木本先生) 聴衆が聞きたいことを意識した方がよい(渡部さん)

19

20

21 MIROC5 (26.0%) MIROC6 (32.2%) OLR[W/m2] MIROC5は14+11=25%
Kim et al. 2009によれば SPCAM(AGCM): =34.6% ECHAM4/OPYC(flux adjustment): =35.5% OLR[W/m2]

22 Wavenumber-frequency spectra (Wheeler & Kiladis, 1999)
OLR(15S-15N) kfp=0.3 MIROC6(kfp=0.5) kfp=0.7 (1<ss(pn).and.ss(pn)<3).and. & ! MJO signal Kim et al. (2014) 30<ff(pt).and.ff(pt)<60)

23 MJO metric

24 平均場の加湿 赤：下層雲 青：中層雲 オレンジ：上層雲 浅い積雲(色)と乱流スキーム(線) 800hPa 20S-20N MIROC5
[g/(kg day)] ISCCP 下層雲量 バイアス [%]

25 平均場のドライバイアス(0-360E) MIROC5 MIROC6 平均場比湿(線)と そのバイアス(色) 浅い対流(色)&
(refはERAI) [g/kg] 浅い対流(色)& 乱流スキーム(線) による加湿 [g/(kg day)]

26 平均場蒸発量 雲の高さの比較とか

27 浅い積雲で増えた水蒸気が貿易風で西に流れるから？
議論1：なぜ西太平洋で降る様になった？ 浅い積雲で増えた水蒸気が貿易風で西に流れるから？ 降水量 乱流と浅い積雲スキーム による加湿率(色) 水蒸気フラックス(矢) GPCP MIROC5 MIROC5 MIROC6 MIROC6 800hPa

28 90-100percentileと、0-20percentileのRh700-850の差
Kim et al. 2014 MIROC6 MJO東進メトリック MIROC5 Pivs_amp_yr: hPa, X=10, Y=20, 0-360E,30S-30N,ocn ERAI-t85.txt: metric: n= GPCP1DD-t85.txt: metric: n= JRA55-t85.txt: metric: n= RAD6s.txt: metric: n= kfp03.txt: metric: n= kfp07.txt: metric: n= pictl.txt: metric: n= wk99_amp: ERAI-t85_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= JRA55-t85_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= NOAA-OLR-t85_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= RAD6s_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= kfp03_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= kfp07_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= pictl_rlnt_yr_b090_ns014_trg096_olp060.txt: amp0= ERAI-t85_pr_yr_b090_ns014_trg096_olp060.txt: amp0= GPCP1DD-t85_pr_yr_b090_ns014_trg096_olp060.txt: amp0= JRA55-t85_pr_yr_b090_ns014_trg096_olp060.txt: amp0= RAD6s_pr_yr_b090_ns014_trg096_olp060.txt: amp0= kfp03_pr_yr_b090_ns014_trg096_olp060.txt: amp0= kfp07_pr_yr_b090_ns014_trg096_olp060.txt: amp0= pictl_pr_yr_b090_ns014_trg096_olp060.txt: amp0= 対流と大気湿度の結合メトリック 90-100percentileと、0-20percentileのRh の差

29 古い！議論2：対流と自由大気湿度の関係は弱化 Prwd_land 30S-30N海上
TRMMPRv7 GPCP1DD ERAI MIROC6 MIROC5 rhCol_ocn 降水量[mm/day] rhCol_ocn jc, j0, ja, je 可降水量[mm] rhCol_land

30 ISCCP cld MJO偏差 観測 MIROC5 MIROC6 西太平洋 東インド洋

31 上層、中層（、下層）の雲量バイアスは全球平均でいずれも低減。
ISCCP雲量バイアスの特徴 (tau>0.3, 全球・年平均) MIROC5 MIROC6 /oraid01/ogura/moon/work/miroc6.0/CTL/shell/plot{35,36}.sh [%] “Too thick” の傾向は悪化。 上層、中層（、下層）の雲量バイアスは全球平均でいずれも低減。

32 ISCCP cld: red(shallow); blue(middle); orange(high)
( E, 20S-20N) 観測 MIROC5 MIROC6 (偏差)

33 ISCCP下層雲量（tau>0.3, 年平均）
地理分布の改善点 ISCCP下層雲量（tau>0.3, 年平均） CFMIP-OBS [%] MIROC5 MIROC6 /oraid01/ogura/moon/work/miroc6.0/CTL/shell/plot{17,18,19}.sh 亜熱帯海上の東西コントラストが改善。

34 MJO: 数1000kmの対流活発域が赤道を20-90日周期で東進
Wheeler & Hendon (2004)に従ってCEOF解析で抽出(OLR, u200, u850) 観測 (45%) MIROC5 (25%) MIROC6 (33%) phase2 phase2 phase2 OLR [w/m2] 20oS-20oN →day→ [説明]Madden Julian Ocsillation と呼ばれる現象の表現が改善した。MJOは数1000kmスケールの対流活発域が、赤道を20-90日の周期で東進する現象で、熱帯気候のみならず、日本気候にも大きく影響することが知られている。 Wheeler & Hendon (2004)に従って、EOF解析でMJOモードを抽出した。図は対流活動度の指標であるOLRを示す。左から観測、MIROC5、MIROC6。いずれにおいても数1000kmの対流活発域がインド洋にあり、それが東進するMJOの様子が確認できる。 観測データにおいてはこのMJOモードが熱帯対流変動(20-90日成分)の45%程度を説明する。一方、MIROC5では25%程度でMIROC5はMJO振幅を過小評価している。それがMIROC6では33%程度に改善した。また、MIROC5は対流活発域の東進が少し遅いが、その伝播特性もMIROC6では改善した。MJOシグナルの東進が太平洋でも確認できることになったことも重要な改善点。 一般に、MJOの気候モデルによる再現は難しいことが知られており、CMIP5の世界中の気候モデルのMJO比較において、MIROC5は他のモデルに比べてMJO表現が比較的良いことが知られていた(Kim et al. 2009; Hung et al. 2013)。それがMIROC6ではさらに良くなった。改善理由は調査中だが、MIROC6で新たに導入した浅い積雲スキームによって、モデルの雲表現が改善したことが関係していると考えられる。 [メモ] 観測はNOAA-OLR(AVHRR)とERAI。WheelerHendon(2004)に従って、20-100日フィルターを掛けたdailyのolr, u200, u850の15S-15N平均(季節は全て)を結合EOFして、第1と第2を合わせて伝播性のMJOを抽出した。第1と第2の大小で8つのフェーズを定義した。図示しているのはxy(水平面)がphase2でxtがphase3; MIROC6の印象で選んだ; 揃えた方がよければ、スライド下のphase3のxyの図を使ってください。 phase3 phase3 phase3 OLR [w/m2] 東進特性と振幅が改善 phase3 phase3 phase3