Development of Multi-model Ensemble technique and its application
Daisuke Nohara
APEC Climate Center (APCC), Busan, Korea 2007/2/21 JMA
Development of Multi-model Ensemble technique and its application - - PowerPoint PPT Presentation
Development of Multi-model Ensemble technique and its application Daisuke Nohara APEC Climate Center (APCC), Busan, Korea 2007/2/21 JMA Contents 1. Introduction of APCC 2. Seasonal forecast based on multi-model ensemble 1. Process of
APEC Climate Center (APCC), Busan, Korea 2007/2/21 JMA
To meet the Bogor Goals of free and open trade and investment in the Asia-Pacific region
1998 Oct. 1999 2003 2004 2005
Proposed APCN : 3rd APEC S&T Ministers’ Meeting
Aug.
Approved APCN : 17th APEC ISTWG Meeting
Jan.
Established at the Korea Meteorological Administration
Mar. Sep.
Supported the establishment of APCC : 27th APEC ISTWG Meeting Proposed the establishment of APCC : 4th APEC S&T Minister’s Meeting
Mar. Nov.
Endorsed the establishment of APCC : 1st APEC Senior Official’s Meeting Welcomed the establishment of APCC : 17th APEC Ministerial Meeting Opening Ceremony : 13th APEC Economic Leader’s Meeting
Meteorological Service of Canada, Environmental Canada Beijing Climate Center China Meteorological Administration International Atmospheric Physics, Chinese Academy of Sciences Central Weather Bureau Chinese Taipei International Research Institute for Climate and Society
DJF SON JJA MAM DJF 2 3 2 1 12 11 10 9 8 7 6 5 4 3 1
MAM JJA SON DJF
APCC provides 3 month forecast in each season.
Precipitation T850 Z500 U850, U200 V850, V200 SLP T2m OLR
Outlook: Interpretation and Description of Global/Regional Prediction Dissemination: - Web Information Up-load
Collection: 17 model datasets
Composite: Deterministic Forecast (4 kinds of schemes)
Probabilistic Forecast
Graphic Quality check Application:
Verification:
2 ° × 2.5 ° L34 T62L64 T42L19 T63L18 1.12 ° × 1.4 ° L28 T42L14 4 ° × 5 ° L17 T63L21 T106L21 T63L40 T42L18 4 ° × 5 ° L2 T63L16 1.875 ° × 1.875 ° L50 T47L17 Model Resolution Meteorological Service of Canada MSC Canada Bureau of Meteorology Research Centre POAMA Australia Hydrometeorological Centre of Russia HMC International Research Institute Center for Ocean-Land-Atmosphere Studies IRI COLA National Aeronautics and Space Administration NCEP Coupled Forecast System Main Geophysical Observatory Meteorological Research Institute Seoul National University Korea Meteorological Administration Japan Meteorological Agency Central Weather Bureau Institute of Atmospheric Physics National Climate Center/CMA Organization USA Russia Korea Japan Chinese Taipei China Member Economies NCEP METRI/KMA GCPS/SNU IAP NSIPP/NASA MGO GDAPS/KMA JMA CWB NCC Acronym
with equal weighting.
′ =
i i
F M P 1
′ =
i i
F M P ˆ 1
Simple composite of individual forecasts, after correction by statistical downscaling using CPPM
Optimally weighted composite of individual forecasts. The weighting coefficient is obtained by SVD based regression.
Weighted combination of statistically corrected multi model output
′ =∑
i i iF
a P
′ =
i i iF
M P ˆ 1 α
Simple composite of individual forecasts, after correction by statistical downscaling using CPPM
★
pattern which is related to the target region.
coefficient in training period.
Predictor field X(i, j, t) Predictand Y(t)
XP(t) Y(t)
Y(t) = αXP(t) + β
( model forecast) ( observation)
Actual Data Set (N) Prediction
Optimally weighted composite of individual forecasts. The weighting coefficient is obtained by SVD based regression. Regression SVD technique P O
Weighted combination of statistically corrected multi model output Actual Data Set (N) Synthetic Data Set (N) Superensemble Prediction
E(ε2) –Minimization using EOF
Combine different models
Combine : according to each model’s square root of ensemble size
Na : num. of above-normal Nn : num. of near-normal Nb : num. of below-normal
Merged 3-category distribution
(Chi-square) TEST
2
χ
O : Forecast frequencies E : Hindcast frequencies
i k i i i
E E O
2 1 2
) (
∑
=
− = χ
Under 0.05% significance level
Near-normal Below-normal Above-normal
3 3-
Categorical distribution
3
Climatological PDF Target year PDF PB: Prob. of Below-Normal PN: Prob. of Near-Normal PA: Prob. of Above-Normal
APCC Multi-Model Probabilistic Forecast Temperature at 850 hPa (2006JJA)
Taylor Diagram : 1983~2003, JJA Hindcast Taylor Diagram : 1983~2003, JJA Hindcast
(model performance)
individual model
(spatial distribution of MSSS)
(spatial distribution of MSSS) (spatial distribution of MSSS)
(regionality regionality) )
0.1 0.2 0.3 0.4 0.5 0.6 AUSTRALIA EAST ASIA NORTH AMERICA RUSSIA SOUTH AMERICA SOUTH ASIA
ACC
Anomaly Correlation Coefficient for hindcast in each region (precipitation, JJA, 1983-2003) individual model multi-model ensemble schemes
2004DJF forecast 2005DJF forecast
(ACC skill score)
individual model individual model individual model
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
CWB GCPS GDAPS_F GDAPS_O HMC IRI IRIF JMA METRI MGO NCC NCEP PMME 계열14
Observed Relative Frequency
Globe, Above Globe, Above-
Normal, Precipitation
Reliability Diagram ROC Curve
False alarm rate Forecast Probability
Forecast skill of probabilistic multi-model ensemble is better than that of individual model.
individual model
Hit rate
– Thailand, Philippine, China
– El Nino, PNA, AO, NAO, Monsoon
(figures: from NWS CPC, http://www.cpc.ncep.noaa.gov/)
Bangkok
Correlation coefficient before and after downscaling in each station
0.2 0.4 0.6 0.8 Suphan buri Lop buri Prachin buri Nakhon ratchsima Bangkok Donmuag Sattahip Chata buri Average
Before After
APCC Downscaling MME Forecast for precipitation in Thailand APCC Downscaling MME Forecast APCC Downscaling MME Forecast for precipitation in Thailand for precipitation in Thailand
Precipitation at Bangkok relates to sea level pressure over the western Pacific.
DJF 2006/07 Nino3 Forecast SON 2006 PNA Forecast
* Value : mean over the period indicated to lastest values (x) * Tendency : rate of change over the period indicated to lastest values (dx/dt)
supported by APEC countries in order to reduce economic loss by climate.
probabilistic forecast. After the verification, we choose the best deterministic forecast.
than that of individual model.
downscaling can improve forecast skill.