Ensemble Forecasting of Typhoon Rainfall and Floods over a Mountainous Watershed in Taiwan Hsiao, L.-F., M.-J. Yang, et all, 2013: Ensemble forecasting.

Slides:



Advertisements
Similar presentations
應用聯合分析法探討婦女對子宮頸癌篩檢方案之偏 好 子宮頸癌是國人常見之婦女癌症,民國 95 年的發生率為台灣婦女癌症之第二 位,死亡率為第六位。許多研究證實子宮頸癌篩檢為預防子宮頸癌最有效之 方法,我國婦女的子宮頸癌篩檢率低是導致子宮頸癌發生率及死亡率偏高的 重要原因,顯示如何提升篩檢率為一重要議題。因此本研究擬瞭解婦女受檢.
Advertisements

2011年度汇报 科技部973项目 《日地空间天气预报的物理基础与模式研究》 第六课题组:空间天气预报方法和技术的应用与集成研究
Measures of location and dispersion
牙齒共振頻率之臨床探討 論 文 摘 要 論文名稱:牙齒共振頻率之臨床探討 私立台北醫學院口腔復健醫學研究所 研究生姓名:王茂生 畢業時間:八十八學年度第二學期 指導教授:李勝揚 博士 林哲堂 博士 在口腔醫學的臨床診斷上,到目前為止仍缺乏有效的設備或方法可以評估或檢測牙周之邊界狀態。臨床上有關牙周病的檢查及其病變之診斷工具,
B型肝炎帶原之肝細胞癌患者接受肝動脈栓塞治療後血液中DNA之定量分析
Chapter 8 Liner Regression and Correlation 第八章 直线回归和相关
Chaoping Li, Zhejiang University
綠色創意伙伴Green Creative Partner
版權所有 翻印必究 指導教授:林克默 博士 報告學生:許博淳 報告日期: 2011/10/24. 版權所有 翻印必究 Results and discussion The crystalline peak at 33° corresponds to the diffraction of the (200)
《大洋环流和海气相互作用的数值模拟》 第九讲 中高纬度海气相互作用 周天军
-Artificial Neural Network- Adaline & Madaline
题目 第一作者1,2,第二作者1,3, 及第三作者等 1,4* 1,大学,部门,城市
袁 星 谢正辉,梁妙玲 中国科学院大气物理研究所
Thinking of Instrumentation Survivability Under Severe Accident
Population proportion and sample proportion
地球科學領域開發高中教材及高中生研習活動
模式识别 Pattern Recognition
INTRODUCTION TO TRENCHLESS TECHNOLOGY 管道铺设非开挖技术的介绍
華爾街的物理學家 混沌碰上華爾街.
《大洋环流和海气相互作用的数值模拟》 第九讲 中高纬度海气相互作用 周天军
Creating Animated Apps (I) 靜宜大學資管系 楊子青
Sampling Theory and Some Important Sampling Distributions
Watershed Management--10
Digital Terrain Modeling
HLA - Time Management 陳昱豪.
第一章 敘述統計學.
Coupling TRIGRS and TOPMODEL in shallow landslide Prediction
報告者:紀瑋欣(Wei-Hsin Chi) 指導老師:楊明仁老師(Prof. Ming-Jen Yang)
The role of leverage in cross-border mergers and acquisitions
Reporter: Prudence Chien
Interval Estimation區間估計
參加2006 SAE年會-與會心得報告 臺灣大學機械工程系所 黃元茂教授
Demand Forecasting in a Supply Chain
A Revised Approach to Ice Microphysical Process for the Bulk Parameterization of Cloud and Precipitation SONG-YOU HONG, JIMY DUDHIA, SHU-HUA CHEN January2004,
塑膠材料的種類 塑膠在模具內的流動模式 流動性質的影響 溫度性質的影響
辐射带 1958年:探险者一号、探险者三号和苏联的卫星三号等科学卫星被发射后科学家出乎意料地发现了地球周围强烈的、被地磁场束缚的范艾伦辐射带(内辐射带)。 这个辐射带由能量在10至100MeV的质子组成,这些质子是由于宇宙线与地球大气上层撞击导致的中子衰变产生的,其中心在赤道离地球中心约1.5地球半径。
The Nature and Scope of Econometrics
Word Revision. Word Revision 1.灾难 2.洪水 3.经历 4.引起,导致 5.猛烈的,激烈的 6.埋葬 7.发生 8.袭击 旋转的,循环的 飓风 雷暴 龙卷风 柱状物 热带的 墓地 棺材.
校園地震預警系統的建置與應用 林沛暘.
村镇污水处理的政策导向 住房和城乡建设部村镇建设司 副司长 赵晖.
準確性(Accuracy) 誤差種類 儀器準確度 時間因素 儀器參數.
Chapter 04 流程能力與績效分析.
中国科学技术大学计算机系 陈香兰 2013Fall 第七讲 存储器管理 中国科学技术大学计算机系 陈香兰 2013Fall.
Simple Regression (簡單迴歸分析)
Dual-Aircraft Investigation of the Inner Core of Hurricane Nobert
L. Ruby Leung and Yun Qian 2003 J.Hydrometeor.,4, 報告:陳心穎
Inter-band calibration for atmosphere
An Evaluation of a High-resolution Hydrometeorological Modeling System for Prediction of a Cool-season Flood Event in a Coastal Mountainous Watershed Kenneth.
第九章 摘要.
Journal of Applied Meteorology, 39,
島嶼地形對颱風路徑偏折之影響 Huang, Y.-H., C.-C. Wu and Y. Wang 2011: The
Yang, Y., and Y.-L. Chen, 2008: Mon Wea. Rev., 136,
An Efficient MSB Prediction-based Method for High-capacity Reversible Data Hiding in Encrypted Images 基于有效MSB预测的加密图像大容量可逆数据隐藏方法。 本文目的: 做到既有较高的藏量(1bpp),
Intercomparison of bulk microphysics schemes in model simulations of polar lows Prudence Chien Wu, L., and G. W. Petty, 2010: Intercomparison.
納莉(2001)颱風之數值模擬研究 A Modeling Study of Typhoon Nari (2001)
An organizational learning approach to information systems development
Review of Statistics.
Sensitivity of Orographic Precipitation to Changing Ambient Conditions and Terrain Geometries: An Idealized Modeling Perspective BRIAN A. COLLE 2004 ,JOURNAL.
Dual-Doppler radar analysis of a near-shore line-shaped convective system on 27 July 2011, Korea: a case study J-T. Lee et al. (2014) Tellus Paper Review.
Yang, M.-J., D.-L. Zhang, X.-D. Tang, and Y. Zhang, J. Geophys. Res.
1-2-2 影響 (a) precipitation There is discussion about during LaNina
Class imbalance in Classification
Monsoonal influence on Typhoon Morakot (2009). Part II: Numerical study. Liang, J., L. Wu*, X. Ge, and C.-C. Wu, 2011: Monsoonal influence on Typhoon   Morakot.
利用循环Hybrid ETKF-3DVAR 改进黄海海雾数值模拟的初始场
Reference: Bryan, G. H., and R. Rotunno, 2009, The maximum intensity
Fei Chen and Jimy Dudhia April 2001 (Monthly Weather Review) 報告:陳心穎
簡單迴歸分析與相關分析 莊文忠 副教授 世新大學行政管理學系 計量分析一(莊文忠副教授) 2019/8/3.
Rainfall Simulation Associated with Typhoon Herb( 1996 )near Taiwan
Principle and application of optical information technology
Gaussian Process Ruohua Shi Meeting
Presentation transcript:

Ensemble Forecasting of Typhoon Rainfall and Floods over a Mountainous Watershed in Taiwan Hsiao, L.-F., M.-J. Yang, et all, 2013: Ensemble forecasting of typhoon rainfall and floods over a mountainous watershed in Taiwan. J. Hydrology, doi: http://10.1016/j.jhydrol.2013.08.046, in press. Keywords: Ensemble forecast ; runoff prediction

Ensemble forecast

By Bob Gell ,NOAA

Runoff prediction WASH123D (Yeh et al.) 1-D Stream-River Network 2-D Overland Regime 3-D Subsurface

Outline Introduction Data and methods Meteorological verification Hydrological verification Conclusions

Introduction Regional scale ensemble prediction systems have been developed to address the need for detailed and high-impact weather forecasting with higher spatial resolution (Du et al., 2009, Yamaguchi et al., 2009 and Clark et al., 2010). 2010 results find that cumulus scheme can effectively provide physics perturbations (25 % track error difference in 36-km WRF) The one-way coupled hydrometeorological approach with rainfall forcing from an ensemble mesoscale modeling system was used in this study to predict rainfall and flooding during the landfall of Typhoon Nanmadol (2011). An ensemble forecast that explicitly represents these uncertainties would provide useful quantitative information regarding the probability of the weather systems (Murphy, 1990). Kain et al., 2008 and Weisman et al., 2008, and Clark et al. (2010) indicated that the convection-allowing NWP models with fine horizontal grid spacing provide value-added predictions for severe convective storms and their associated heavy rainfall. Regional scale ensemble prediction systems have been developed in research and operational modes to address the need for detailed and high-impact weather forecasting with higher spatial resolution (Du et al., 2009, Yamaguchi et al., 2009 and Clark et al., 2010). Yamaguchi et al. (2009) showed that the ensemble mean track forecasts for typhoons in the western North Pacific in 2007 had a 40-km error reduction in the 5-day forecasts compared to the deterministic model forecast.

The hydrological responses of most watersheds in Taiwan are fast and complicated due to the steep slopes of the Central Mountain Range (CMR). In this study, the Lanyang creek basin was selected as the target area for watershed modeling . 圖片來源 : Environmental Protection Administration Executive Yuan, R.O.C

Nanmadol became a tropical storm at 1200 UTC 23 August 2011 , and then moved north–northwestward before making landfall in southeastern Taiwan on 28 August. After Nanmadol passed over Taiwan, it rapidly weakened before dissipating over the Taiwan Strait. Typhoon Nanmadol produced heavy rainfall that resulted in agricultural and industry damage and the loss of many lives. Nanmadol became a tropical storm at 1200 UTC 23 August 2011 as it moved westward to northwestward along the southern edge of the subtropical high. Following landfall in the northeastern Philippines at 0000 UTC 27 August, its intensity was reduced from category 3 to category 2 [based on the Saffir–Simpson hurricane scale (Simpson, 1974)]. Nanmadol then moved north–northwestward due to the westward extension of the subtropical high before making landfall in southeastern Taiwan on 28 August (Fig. 1). After Nanmadol passed over Taiwan, it rapidly weakened before dissipating over the Taiwan Strait. From 1200 UTC 27 August to 0000 UTC 30 August, the Central Weather Bureau (CWB) of Taiwan issued typhoon warnings for heavy rainfall and strong winds. In Taiwan, a total property loss of 100 million Taiwan dollars (3.3 million US dollars) resulted from Typhoon Nanmadol. 圖片來源 : 中央氣象局CWB

Data and methods-Observations Lanyang stream watershed 512 automatic rain-gauge stations Rainfall forecast interpolated to each stations using the Kringing technique (Bras and Rodriguez-Iturbe, 1985)

Data and methods-Model setups Three nested domains with 51 vertical levels 18 ensemble members in WRF and MM5 221*127 150*180 183*195 三層巢狀網格 垂直51層以高解析度解析邊界層

解釋 Cold start :以GFS的大尺度預報場做初始化 ,做在解析時間前做12小時cold stat以及六小時資料同化循環 CV3 CV5 : 使用3DVAR同化時,使用的error covariance matrices OL : 使用3DVAR時的outer loop,藉由將算好的Xa帶回去計算以穩定非線性的H Bogus vortices : 給予假渦旋做初始場 MM5的模擬使用4DVAR和noda的模擬 LBCs : provided every 6 hours from NCEP GFS and CWB gfs 積雲參數化主要使用前兩個domain Grell-Devenyi 對於深對流和環境的閉合,為對每個grib都用不同的閉合 Grell 3D ensemble 同Grell 只是對3維每個grib都給予不同的閉合 Betts-Miller-Janjic 藉由柱狀的溼調整來調整成完全混和的profile Kain-Fritsch 對於次網格的深淺對流系統 使用質量通量的向下傳遞和CAPE 來移除時間尺度的影響 Grell scheme 為MM5中較舊版本的Grell scheme 主要是要看基雲參數話對颱風路徑的影響 微物理參數化為了清楚解系對流,在5km也使用 使用Goddard 和 WSM5 Goddard 主要是對冰 雪 grapel有高解析度模擬 WESM5則是對混相位過程和過冷水有較多著墨 PBL使用YSU和MRF 系集結果每六小時輸出一次

Data and methods- Skill score Observed Yes no Forecast Hits False alarms Misses Correct negatives Threat score (TS) Equitable threat score (ETS) H : Hits F : Forecast yes O : Observed yes

Data and methods- Skill score Observed Yes no Forecast Hits False alarms Misses Correct negatives Bias score (BS) False alarm rate (FAR) Standard deviations (SD) H : Hits F : Forecast yes O : Observed yes

Meteorological verification To establish the veracity of the track forecast ensemble system,219 forecasts from 21 typhoons in 2011 were verified relative to the observed (CWB best-track analysis) TC positions. The ensemble track forecasts of Nanmadol were better than the average for the 21 typhoons in 2011. F4 The ensemble mean track errors for the 21 typhoons were 93, 180, 295 km at 24, 48, and 72 h forecasts, respectively. These track error values were superior to the Navy Operational Global Atmospheric Prediction System (NOGAPS) values of 140, 216, and316 km. In the next section, these ensemble rainfall forecasts for Typhoon Nanmadol and the flood simulations from the WASH123D hydrological model are evaluated.

Rainfall forecast skill parameters Fig. 5. Box-and-whisker plot of (a) threat score (TS), (b) bias score (BS), (c) equitable threat score (ETS), and (d) false alarm rate (FAR) for 24-h accumulated rainfall forecast of the 18 individual members at the 130-mm threshold for Typhoon Nanmadol from 1200 UTC 27 to 0000 UTC 30 August during which the CWB issued typhoon warnings. Values of the ensemble mean of rainfall forecast are indicated with dots. The box-and-whisker plot is interpreted as follows: the middle line shows the median value; the top and bottom of the box show the upper and lower quartiles (i.e., 75th and 25th percentile values); and the whiskers show the minimum and maximum values. 可以看到系集平均比單一member的預報分數來的高 BS分數隨著颱風登陸台灣(1800UTC AUG28) 而上升並在離開後下降(0000UTC AUG29) 系及預報的降水大部分是高估的,主要是因為颱風接近時迎風面地形和颱風環流的交互作用的不確定性 颱風的離開帶動了FAR的下降 以及ETS的上升但隨後FAR就更多 因為台灣離開陸地後開始減弱且結構被破壞,使得模擬的不確定性增加,以及海上觀測少 Extremely heavy rain (豪雨) : 24-hour accumulated rainfall exceeds 130 millimeters

0-24 h fcst rainfall by each member Obs F8a The 18 individual ensemble members for the 0–24 h accumulated rainfall forecasts are provided in Fig. 8a. In these ensemble forecasts, the maximum 24-h rainfall occurs over eastern Taiwan, except for the M02 and M05 members. except for the M02 and M05 members. The forecast tracks for these two members were two outliers 2. The two MM5 ensemble members (M17 and M18) predicted that more rainfall would occur in eastern Taiwan than the two members (M14 and M15) that were based on the WRF model with similar 0–24 h forecast tracks. =>This early rainfall forecast from the two MM5 members resulted from a more rapid translation speed. Thus, the interaction between the typhoon circulation and the Central Mountain Range occurred earlier Ensemble mean probability distribution

24-48 h fcst rainfall by each member Obs F8b 因為系集TC的移動緩慢,造成中央山脈東側累積降雨高估 南部的降雨低估 且主要降雨區及中在東部 而非屏東的488mm Ensemble mean probability distribution

Forecasted tracks by ensemble members and the observed track Compared two cases : Initiated at 1200 UTC 27 August Initiated at 1200 UTC 28 August ∵rainfall variability are sensitive to typhoon track F6 To examine the improvements in the ensemble rainfall forecasts and runoff simulations, two cases for the forecasts that were initiated at 1200 UTC 27 August and at 1200 UTC 28 August (Fig. 6) were selected. These forecasts contained large and small rainfall variabilities that likely result from large and small typhoon track forecasts variabilities, respectively.

Initiated at 1200 UTC 27 August (a)Obs. 24-h rainfall (b)Fcst. 24-h rainfall by ensemble mean (c)Probability of 24-h rainfall >130mm (a,b,c) 0-24h (d,e,f) 24-48h F7 Fig. 7. The 0–24-h accumulated rainfall from the forecast initiated at 1200 UTC 27 August: (a) observed rainfall, (b) ensemble mean from 18 members, and (c) the rainfall probability distribution (%) exceeding the threshold of 130 mm for 18 ensemble members. The observed rainfall at the 130-mm threshold is shown by the black solid lines. (d,e, and f) as in panels (a, b, and c), except for the 24–48-h accumulated rainfall. 0-24 h 觀測到的豪大雨在東台灣地區 而24-48h 則在東及南台灣 從左圖c中,看到可能性50%以上的0-24h豪大雨降水都及中在東部地區 而對於24-48h預測 則70%以上豪大雨降水都在東部,只有30%在南部 這表示台灣的降雨預測受地形效應影響很深 而比較游28日開始模擬的系集結果,因為路徑預報較為準確的關係,對於各時段的降雨pattern上掌握的準確, 看到圖c的豪大雨區基本上都跟觀測相符 但是雨量則高估許多

Time series of 3-h rainfall for (a) three basins over southern Taiwan and (b) Lanyang basin Initiated at 1200 UTC 27 August Initiated at 1200 UTC 28 August F9 aug27 Time series of areal-average 3-h rainfall (in units of mm) for (a) three basins over southern Taiwan and (b) Lanyang basin from the ensemble members with minimum, lower quartile, median, upper quartile, and maximum depicted by boxand- whiskers plot from the forecast initiated at 1200 UTC 27 August, and the ensemble mean (MEAN; gray solid line), the rainfall observations (OBS; black solid line), and standard deviation (SD; black dash line). 看到27日開始模擬的降雨 在南部地區三個集水區的雨量上,28 1800和29 0300過小外 其他的趨勢都有抓到, 而對蘭陽溪集水區的降雨則略有所高估 而28日開始模擬的降雨,在28日2100和30日0000的高值沒有抓得很精準 雖然前一個峰直有抓到,但第二個卻完全沒有 以系集平均來說再第二個peak也過於低估 而在蘭陽溪流域的模擬,在第一個峰值的部分延遲了六了小時,但若降雨的峰值還是有抓到

Horizontal distribution of radar reflectivity at 00 UTC 29 August Initiated at 1200 UTC 28 August (12h ) F11 Aug28 Fig. 11. Spatial distribution of radar reflectivity (dBZ) from observations (OBS) at 0000 UTC 29 August (left panel) for the 18 ensemble members (right panels) at 12 h in the forecast initiated at 1200 UTC 28 August. 從雷達迴波圖可以發現,大部分member的回波形狀都相似,但相較於觀測回波則略偏東北 但因為中央山脈地形的關係而在贏風處產生了系統性的高估 MM5的迴波更是明顯

Horizontal distribution of radar reflectivity at 00 UTC 30 August Initiated at 1200 UTC 28 August (36h ) F12 aug28 Fig. 12. As in Fig. 11, except for the observations (OBS) at 0000 UTC 30 August (left panel) and for the 18 ensemble members at 36 h in the forecast initiated at 1200 UTC 28 August. 從30號的迴波可以清楚地看到,在中央山脈南側地區有一條狹長的雨帶, 但是在模擬中都沒辦法抓到這一條回波, 因為grid size最小網格為5公里無法解析的關係

Data and methods-Hydrological model WASH123D (Yeh et al.1998) Finite-element approach Terrain spatial resolution 400m*400m Finer grids : 40m*40m Interpolated 5-km rainfall from atmospheric model using nearest neighbor interpolation. River and overland: Diffusive wave equations Infiltration : Green–Ampt model Coastal inundations : semi-Lagrangian and Galerkin finite-element methods An integrated watershed simulation that includes groundwater calculations for flood forecasting has not been considered as a practical alternative in Taiwan due to its steep terrain and the resulting short hydraulic response times. Therefore, groundwater routing was ignored in this WASH123D hydrological model.

Hydrological verification Initiated at 1200 UTC 27 August During Typhoon Nanmadol, most of the rainwater was presumed to have infiltrated into the groundwater. Thus, excess overland flow was assumed to move slowly due to the dry soil conditions. Rain gauge observations Initiated at 1200 UTC 28 August F14 Comparison of water stages (m) between the measurements and the WASH123D simulation driven by rain gauge observations starting from 1200 UTC 27 August. 看到觀測水位實際上是緩慢上升,而預報水位則是快速上升後,就變得低估 F16 Hourly time series of the areal-averaged water stage (in units of m) for the Lanyang basin estimated from the ensemble members with minimum, lower quartile, median, upper quartile, and maximum depicted by box-and-whiskers plots for the ensemble forecasts initiated at (a) 1200 UTC 27 August and (b) 1200 UTC 28 August, and hourly water stage from ensemble mean (MEAN; gray solid line), the observation (OBS; black closed circles), and standard deviation (SD; black dash line). Just as the WASH123D model over-predicted the water stage when driven by the rain gauge data, river runoff was over-forecasted when driven by the ensemble rainfall forecast. When the river water stage was high, more precipitation directly affects rainfall because the soil is moist. In general, this integrated hydrometeorology modeling system is useful for predicting (albeit a likely over-forecast) the occurrence of extreme floods during typhoon events in the mountainous watersheds on the windward side of Taiwan. This result can be used in other mountainous watersheds by using hydrological models that are familiar based on local soil conditions.

The 48-h simulations by 18 ensemble members Initiated at 1200 UTC 28 August Initiated at 1200 UTC 27 August F15a water stage was over-predicted by the ensemble members. F15b water stages were lower and were closer to the observed water stages.

Conclusions In 2011, the ensemble provide a better track prediction than those of operational centers. 90% probability that accumulated rainfall exceeded 130mm for 0-24 h forecast at 1200 UTC 28 August is in good agreement with the distribution of observed 130-mm rainfall. Ensemble forecasting system adequately estimated the topographic locations where rainfall may occur.

In this case, the river runoff patterns were reasonably predicted despite the mismatch between the runoff maximum and the actual time and quantity of flooding. The omission of a ground water routing component in the watershed model contributed to the over-prediction of river runoff. Despite the systematic over-prediction of rainfall and water stage in the watershed on the windward side of Taiwan, the coupled hydrometeorological modeling system can potentially improve the accuracy and timing of flood predictions.

References Hsiao, L.-F., M.-J. Yang, et all, 2013: Ensemble forecasting of typhoon rainfall and floods over a mountainous watershed in Taiwan. J. Hydrology, doi: http://10.1016/j.jhydrol.2013.08.046, in press. Environmental Protection Administration Executive Yuan, R.O.C Nasrollahi, Nasrin, Amir AghaKouchak, Jialun Li, Xiaogang Gao, Kuolin Hsu, Soroosh Sorooshian, 2012: Assessing the Impacts of Different WRF Precipitation Physics in Hurricane Simulations. Wea. Forecasting, 27, 1003–1016. Wandishin, Matthew S., Steven L. Mullen, David J. Stensrud, Harold E. Brooks, 2001: Evaluation of a Short-Range Multimodel Ensemble System.Mon. Wea. Rev., 129, 729–747. Clark, Adam J., and Coauthors, 2011: Probabilistic Precipitation Forecast Skill as a Function of Ensemble Size and Spatial Scale in a Convection-Allowing Ensemble. Mon. Wea. Rev., 139, 1410–1418. Hamill, Thomas M., 1999: Hypothesis Tests for Evaluating Numerical Precipitation Forecasts. Wea. Forecasting, 14, 155–167. 《李天浩,2009:應用克利金法建立高解析度網格》

Thanks for your attention

Outer loop