A Training Course on Quantitative Precipitation Estimation/Forecasting (QPE/QPF) Crowne Plaza Manila Galleria, Quezon City, Philippines 27-30 March 2012
Frequently used QPF techniques and developments Extrapolation, Climatology, Similarity
Kuo-Chen Lu Central Weather Bureau, Chinese Taipei
Training Course on QPE/QPF Manila 2012 1 contents
Introduction Philosophy Synoptic forcing Frontal system Development of forecast guide Typhoon Climatology approach
Rain gauge and Grid base Dynamic model approach
Development of forecast guide
Training Course on QPE/QPF Manila 2012 2 philosophy
QPF : is the expected amount of melted precipitation accumulated over a specified time period over a specified area. Amount: Maximum value, or the mean value ? Time period : Hourly, Daily or, A Storm lifespan. Specified Area : A rain gauge, a river basin, A township, A high risk area. How they evaluate QPF ? What’s the reporter said on the newspaper ? They evaluated the QPF is based on the rain gauges. So the forecaster should forecast the value to meet the reporter need, that is based on the rain gauges. Besides, they show the maximum value first. What’s the type of QPF that the user need ? The type of the value should be in the range or in probability. Daily or hourly, the spatial resolution and the duration of request.
Training Course on QPE/QPF Manila 2012 3 Forecasting QPF
Must determine Where When How Much rainfall will occur area. Must understand the processes that determine the size, scale and intensity of an area of precipitation (synoptic, mesoscale, and even microscale meteorology) Must ..Possess Good Pattern Recognition Skills and understand what gives the pattern the potential to produce significant rainfall Must Possess a working Knowledge of Local Climatology Understand numerical models especially model biases and why they occur
Training Course on QPE/QPF Manila 2012 4 Subjectively QPF
Analyze situation looking at the synoptic and mesoscale environment. (use the current data and model output to assess situation Does the environment favor high rainfall rates. Use model guidance as a first guess but understand model limitations and biases. Based on model output, radar and satellite imagery and conventional upper air data, try to figure out where rainfall will be most intense for the longest period. This is where the heaviest rainfall will occur. Modify model guidance based on your understanding of the physics that determines how much rainfall will fall. Calibrate forecasts through verification. Verification is very important for manual and computer generated forecast. But how do you verify a forecast.
Training Course on QPE/QPF Manila 2012 5 Philosophy Synoptic Checking Large Scale checking on a specific area, with regard to the itemized synoptic forcing. Dynamic guidance Numerical Model may perform well on synoptic forecast, but not on the convective scale, especially on the heavy rain fall. It is worth to evaluate the synoptic forcing by NWP. Decision Making A check list of synoptic forcing by the NWP might help us to expect when and where the heavy rain will occur. Training Course on QPE/QPF Manila 2012 6 Synoptic Checking of Frontal System
1. Surface front 2. Low level horizontal wind shear 3. Sub synoptic system 4. Moist Field 5. Low level Jet 6. Cold tongue 7. Pressure depression 8. Upper level diference 9. Instibility Composited chart 0000 UTC June 1, 2000
Training Course on QPE/QPF Manila 2012 7 Terrential Rainfall Check list Table for 12 48 hrs during the Mei-Yu Season
Initial time Year______Mon.______Day______Hr.______Z NWP: WRF JMA
Items of check list OBJ 12H r 24Hr 36Hr 48Hr remark
1. Surface front location within 20 N 28 N 118 E 124 E ( ) ( ) ( ) ( ) ( )
Taipei near 100Km or passed 200Km ( ) ( ) ( ) ( ) ( ) KuoHsung near 200Km ( ) ( ) ( ) ( ) ( )
2. Horizontal wind shear 22 N 28 N 114 E 127 E
850/700hPa ( ) ( ) ( ) ( ) ( )
3. Sub synoptic South of Taiwan or east of 114 E surface/850hPa Meso-low ( ) ( ) ( ) ( ) ( ) 700/500hPa Short wave trough ( ) ( ) ( ) ( ) ( )
4. Moist fields 850hPa Td 15 ( ) ( ) ( ) ( ) ( )
850hPa e axi cross Taiwan ( ) ( ) ( ) ( ) ( ) 700hPa T-Td 3 ( ) ( ) ( ) ( ) ( )
5. Low level Jet 18 N 26 N 115 E 125 E surface 10 20kts southwestly wind 850hPa 25kts southwestly wind ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 700hPa 30kts southwestly wind ( ) ( ) ( ) ( ) ( ) 850hPa 10kts south southwestly within north of ( ) ( ) ( ) ( ) ( ) South China Sea 15 N
6. Thermal Fields 850hPa cold tongue/warm tongue arount the wind shear ( ) ( ) ( ) ( ) ( ) 700hPa cold tongue/warm tongue arount the wind shear ( ) ( ) ( ) ( ) ( )
7. Pressure depression (around Taiwan) Surface pressure P 24hr -3hPa 850hPa H 24hr -15gpm ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Sea level pressure 1005hPa ( ) ( ) ( ) ( ) ( ) Low pressure around Hai-Nan Low Pressure ( ) ( ) ( ) ( ) ( ) 1004hPa
8. Upper level wind difluence
300/200hPa angle of difluence 45 ( ) ( ) ( ) ( ) ( )
9. Stability Total Index 40 ( ) ( ) ( ) ( ) ( ) K-Index 35 ( ) ( ) ( ) ( ) ( )
Summary Yes ( ) 22
Training Course on QPE/QPF Manila 2012 8 Terrential Rainfall Check list Table for 12 48 hrs during the Mei-Yu Season
Initial time Year______Mon.______Day______Hr.______Z NWP: WRF JMA
Items of check list OBJ 12H r 24Hr 36Hr 48Hr remark
1. Surface front location within 20 N 28 N 118 E 124 E ( ) ( ) ( ) ( ) ( )
Taipei near 100Km or passed 200Km ( ) ( ) ( ) ( ) ( ) KuoHsung near 200Km ( ) ( ) ( ) ( ) ( )
2. Horizontal wind shear 22 N 28 N 114 E 127 E
850/700hPa ( ) ( ) ( ) ( ) ( )
3. Sub synoptic South of Taiwan or east of 114 E surface/850hPa Meso-low ( ) ( ) ( ) ( ) ( ) 700/500hPa Short wave trough ( ) ( ) ( ) ( ) ( )
4. Moist fields 850hPa Td 15 ( ) ( ) ( ) ( ) ( )
850hPa e axi cross Taiwan ( ) ( ) ( ) ( ) ( ) 700hPa T-Td 3 ( ) ( ) ( ) ( ) ( )
5. Low level Jet 18 N 26 N 115 E 125 E surface 10 20kts southwestly wind 850hPa 25kts southwestly wind ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 700hPa 30kts southwestly wind ( ) ( ) ( ) ( ) ( ) 850hPa 10kts south southwestly within north of ( ) ( ) ( ) ( ) ( ) South China Sea 15 N Training Course on QPE/QPF Manila 2012 9
Items of check list OBJ 12H r 24Hr 36Hr 48Hr remark
6. Thermal Fields 850hPa cold tongue/warm tongue arount the wind shear ( ) ( ) ( ) ( ) ( ) 700hPa cold tongue/warm tongue arount the wind shear ( ) ( ) ( ) ( ) ( )
7. Pressure depression (around Taiwan) Surface pressure P 24 hr -3hPa 850hPa H 24 hr -15gpm ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Sea level pressure 1005hPa ( ) ( ) ( ) ( ) ( ) Low pressure around Hai-Nan Low Pressure ( ) ( ) ( ) ( ) ( ) 1004hPa
8. Upper level wind difluence
300/200hPa angle of difluence 45 ( ) ( ) ( ) ( ) ( )
9. Stability Total Index 40 ( ) ( ) ( ) ( ) ( ) K-Index 35 ( ) ( ) ( ) ( ) ( )
Summary Yes ( ) 22
Heavy raainfall
(50mm < PQF <130mm a day) resuilts Terrential Rainfall
Training(QPF Course >= on 130mm QPE/QPF a day)Manila 2012 10
EX: Surface Fronts 20°N~28°N,118°E~124°E ( x ) Taipei near 100Km or passed 200Km ( x ) KuoHsung near 200Km ( - )
Training Course on QPE/QPF Manila 2012 11 EX: Low level Jet(18°N~26°N,115°E~125°E)
surface 10~20kts southwestly wind 850hPa >25kts southwestly wind 700hPa >30kts southwestly wind
Training Course on QPE/QPF Manila 2012 12 Composited chart 0000 UTC June 14, 1997
QPE 15 items selected
Training Course on QPE/QPF Manila 2012 13 Composited chart 0000 UTC May 18, 1998
QPE 9 items selected
Training Course on QPE/QPF Manila 2012 14 Composited chart 0000 UTC May 19, 1999
QPE 11 items selected
Training Course on QPE/QPF Manila 2012 15 Composited chart 0000 UTC June 12, 2000
QPE 15 items selected
Training Course on QPE/QPF Manila 2012 16 Distribution of the number of items for the synoptic forecing
15.0% 14.3% 13.9%
12.0% percentage
9.0% 8.3% 7.8% 7.8% 7.4% 6.5% 百分比 6.0% 5.7% 5.7% 5.7% Rare events 4.3%
3.9%
3.0% 2.6% 2.6% 1.7% 1.3% 0.4% 0.0% 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Number項目數 of items
Training Course on QPE/QPF Manila 2012 17 Distribution of the number of items for the frequency of heavy rainfall
100% 2500
90% Blue 大雨發生頻率bar : freq. of
80% heavy發生總雨量 rain 2000 black line the is total rain fall for stations every
70%
frequency heavyrain 60% 1500
50%
1000
40% 總雨量值 發生大雨頻率
30%
20% 500
10%
0% 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 Number勾選因子數目 of items
Training Course on QPE/QPF Manila 2012 18 the correlation among the synoptic items
Pressure Upper Synopic Level Sub Moist Low Thermo Average Front Depressio level Instability items Shear System Fields Level Jet Fields of items n difluence Average 17% 8% 6% 23% 14% 7% 4% 4% 17% 5.7 Rain 16% 8% 7% 22% 15% 6% 3% 4% 19% 4.1 Heavy 19% 7% 6% 24% 12% 8% 5% 4% 15% 8.6 rain Terrenti 20% 6% 6% 23% 13% 6% 6% 5% 15% 10.6 al Rain Correlat ion 0.55 0.35 0.26 0.60 0.31 0.29 0.25 0.35 0.31 0.65 Coef.
Training Course on QPE/QPF Manila 2012 19 Verificaion of the number of synoptic items (forecast v.s. Analyis )
21 21
18 12 hrs 18 24 hrs 十 二 二 15 十 15 小 四 時 12 小 12 預 時 測 9 預 9 項 測 目 6 項 6 數 目 3 數 3
0 0 3 6 9 12 15 18 21 0 0 3 6 9 12 15 18 21 客觀分析項目數 客觀分析項目數 21 36 hrs 21 t s a c e r o F 18 18 48hrs 三 四 十 十 15 15 六 八 小 小 12 12 時 時 預 預 9 9 測 測 項 項 6 6 目 目
數 數 3 3
0 0 0 3 6 9 12 15 18 21 0 3 6 9 12 15 18 21 客觀分析項目數 客觀分析項目數 a n a lTraining ys i sCourse on QPE/QPF Manila 2012 20 Verificaion for differenct valied time
Training Course on QPE/QPF Manila 2012 21 Summary of Frontal System
There are good relation between the rain event and Synoptic forcing, if the forcing can be quantitatively define. Find an appropriated way to use NWP can help the forecast of heavy rain event subjectively.
Training Course on QPE/QPF Manila 2012 22 Development of Computer Interface
Training Course on QPE/QPF Manila 2012 23 Historical data base of QPE
Training Course on QPE/QPF Manila 2012 24 Samples of difference combinaiton of synoptic items.
Training Course on QPE/QPF Manila 2012 25 THE FLOW CHART OF SUBJECTIVE QPF
Interface for
Environment Environment Historical Synoptic Synoptic Database similarity Checking Analysis + NWP
Fine tune
Rainfall Decision Database Making
Training Course on QPE/QPF Manila 2012 26 Philosophy Climatology guidance --- for day 3, day 4 forecast Large Uncertainty for the track and the structure of Typhoon Dynamic guidance --- for day2 and day 0 forecast High resolution model (WRF Ensemble) is available Similarity or analogue --- for the day 1 and day 0 Radar echo and QPE is available for calibration
Training Course on QPE/QPF Manila 2012 27 The mechanism of typhoon precipitation
1. convergence of typhoon circulation (radius Spiral band dependent),
2. Terrain slope lifting of typhoon circulation (location Inner rain ban dependent), Topographic lifting 3. convections within eye wall Eye Wall and rain bands (transient)
4. interactions with environmental flows (e.g. NE or SW monsoon).
Training Course on QPE/QPF Manila 2012 28 The terrain slope lifting (topographically locked) – key of typhoon rainfall climatology model
Training Course on QPE/QPF Manila 2012 29 Training Course on QPE/QPF Manila 2012 30 Training Course on QPE/QPF Manila 2012 31 Collection of Historical case for climatological QPF model. Locations of 371 rain gauges
CWB: 340 ; WRB: 8 ; Reservoir: 23
To study the rainfall characteristics in Taiwan during typhoon period, we analyzed all hourly rainfall data taken at 371 rain gauges for 91 typhoons in 1989~2007. Training Course on QPE/QPF Manila 2012 32 115 120 125 130 135 140 35 35
30 30 Typhoon Nielson (1985) 08/23
25 08/24 08/21 25 08/19 08/22 08/20
20 08/18 20 33 Training Course on QPE/QPF Manila 2012
115 120 125 130 135 140 120 125 130 135 140
30 30
10/03 Typhoon Longwang (2005) 10/02 25 25 09/29 09/30 09/27
10/01 09/28 20 20 34 Training Course on QPE/QPF Manila 2012 09/26
15 15
120 125 130 135 140 110 115 120 125 130 135
30 30
Typhoon Dujan (2003) 25 25 09/03
09/01 09/02 08/31 20 20 35 Training Course on QPE/QPF Manila 2012
08/30 15 15 110 115 120 125 130 135 Climatology by Track Type
The track pattern is the key point for the topographic block of rainfall
Training Course on QPE/QPF Manila 2012 36 Training Course on QPE/QPF Manila 2012 37 Track Type 1, Moving Westward, and passed the sea of Northern Taiwan
Average rainfall on the 25 Synoptic weather sites around Taiwan Training Course on QPE/QPF Manila 2012 38 Track Type 2, Moving Westward and Landing on Taiwan
Average rainfall on the 25 Synoptic weather sites around
Training Course on QPE/QPF ManilaTaiwan 2012 39 Track Type 6,7 Moving Northward and passed Eastern/Western sea of Taiwan
Training Course on QPE/QPF Manila 2012 40 Collection of Historical case for climatological QPF model. Locations of 371 rain gauges
CWB: 340 ; WRB: 8 ; Reservoir: 23
To study the rainfall characteristics in Taiwan during typhoon period, we analyzed all hourly rainfall data taken at 371 rain gauges for 91 typhoons in 1989~2007. Training Course on QPE/QPF Manila 2012 41 The climatology model contains a set of 371 typhoon rainfall climatology maps, one for each rain gauge station
Maximum Data counts
Taipei (4669 Mean Standard 2) deviation The rainfalls at a station 0.5o x 0.5o grid box (e.g. 46692, Taipei) were averaged15 mm/h when the typhoon 50 mm/h centers were located within each 0.5o x 0.5o grid box (256 grid boxes) to represent Use Barnes (1973) scheme the climatology at that box. to interpolate to 0.1o x 0.1o → Typhoon rainfall climatology map
Training Course on QPE/QPF Manila 2012 42 Climatology Approach
66 1o x 1o Mean rainfall rate
Training Course on QPE/QPF Manila 2012 43 Climatology
29 31 Approach
26 47 55 57
36 68 66 66
1o x 1o Mean rainfall rate Training Course on QPE/QPF Manila 2012 44 Climatology
29 31 37 29 Approach
26 47 55 57
36 68 66 66
54 79 103 59
36 29 55 80 1o x 1o Mean
9 23 65 59 rainfall rate Training Course on QPE/QPF Manila 2012 45 Rainfall estimated based on climatology rainfall forecast and CWB 120-h track forecast 2009080603UTC
2009080612UTC
Training Course on QPE/QPF Manila 2012 From Lee et al 46 Rainfall estimated based on climatology rainfall forecast and CWB 120-h track forecast 2009080700UTC
2009080712UTC
Training Course on QPE/QPF Manila 2012 47 Model can provide hourly and cumulative rainfalls when integrating along the typhoon track.
Observed 26.9mm481.8 25 443.5450 21.3mm 400 20 350 300 15 Seth (1994) 250 10 200 150 5 100 1.009 50
0 0
1008 08 1008 10 1008 12 1008 14 1008 16 1008 18 1008 20 1008 22 1008 24 1008 02 1009 04 1009 06 1009 08 1009 10 1009 12 1009 14 1009 16 1009 18 1009 20 1009 22 1009 24 1009 02 1010 04 1010 06 1010 08 1010 10 1010 12 1010 14 1010 16 1010 18 1010 20 1010
Training Course on QPE/QPF Manila 2012 48 46693 Lynn (1987) 200 175 3 hr rainfall
150 (mm)
125
100
75
50
25
0 2300 2312 2400 2412 2500 2512 2600
Training Course on QPE/QPF Manila 2012 49 Climatology model underestimated seriously where the interactions taking place ( Lynn, 1987).
Observations Rainfall (blue) vs. model results (red)
Taitong station Chu-tze-hu station
Training Course on QPE/QPF Manila 2012 50 Summary of Climatology Model
Motivation For a given typhoon, the distribution and amount of rainfall over Taiwan is highly correlated to the relation between typhoon position and Taiwan island, it is so called topographic locking. It is workable to provide a statistical based rainfall estimate according to a given typhoon track forecast
Weakness The limited historical typhoon cases Difficult to resolve the case dependent variability of environmental situation and mesoscale precipitation structure.
Lee et al. 2006 51 Training Course on QPE/QPF Manila 2012 52 Numerical model forecast Strength If track forecast is OK, the model QPF is worth to be a guidance Is able to represent the interaction between the environment and typhoon circulation Depend on case, has the potential to capture the mesoscale precipitation process Weakness Uncertainty to the track forecast Official track forecast still better than the model forecast Uncertainty from the initial condition and the model physical process Limitation to the model resolution Hard to configure a model that performs THE BEST all the time
How to maximally take advantage of the strength and well handle the uncertainty of the weakness from the model QPF?
Training Course on QPE/QPF Manila 2012 53 What can we do to extract the useful model TY QPF information from ensemble forecasts? Training Course on QPE/QPF Manila 2012 54 • Select the model QPF cases from ensemble members according to the prior estimate of the typhoon position, then produce the composite rain map and probability products based on the selected samples. • Maximally use the ensemble QPF based on the optimal track forecast • More members are required to provide the enough spread of the forecast typhoons. • “Realtime” is not a critical issue, so that can join all the model resources from Taiwan community, e.g. the ensemble typhoon forecast experiment conducted by TTFRI in 2010. • Climatology Model ensemble typhoon QPF (ETQPF) • Using the near-realtime ensemble typhoons instead of the real typhoons. Training Course on QPE/QPF Manila 2012 55 OBS
Examples of the member QPFs with similar typhoon location Training Course on QPE/QPF Manila 2012 56 Observation
Ty Megi (2010)
Strong interaction between the typhoon circulation and north-east wind Training Course on QPE/QPF Manila Ensemble2012 based Typhoon QPF 57 Ty Fanapi (2010)
Climatology Observations ETQPF 12-hr accumulated rainfall official track : 2010/09/18 00Z valid time from 09/18 12Z Training Course on QPE/QPF Manila 2012 to 09/1958 00Z Ty Fanapi (2010)
Climatology Station Obs ETQPFs 12-hr accumulated rainfall official track : 2010/09/18 06Z valid time from 09/18 18Z Training Course on QPE/QPF Manila 2012 to 09/1959 06Z Ty Fanapi (2010)
Climatology Station Obs ETQPFs 12-hr accumulated rainfall official track : 2010/09/18 06Z valid time from 09/19 06Z Training Course on QPE/QPF Manila 2012 to 09/1960 18Z Ty Fanapi (2010)
Climatology Station Obs ETQPFs 12-hr accumulated rainfall official track : 2010/09/18 12Z valid time from 09/18 12Z Training Course on QPE/QPF Manila 2012 to 09/1961 00Z Ty Fanapi (2010)
Climatology Station Obs ETQPFs 12-hr accumulated rainfall official track : 2010/09/18 12Z valid time from 09/19 00Z Training Course on QPE/QPF Manila 2012 to 09/1962 12Z Observation
Better estimate of the typhoon track Better QPF
Training Course on QPE/QPF Manila 2012 63 Systematic errors of ETQPF based on the best track
Megi (O-F) Fanapi (O-F)
Over prediction over mountain range
Under prediction: Smooth out due to the composite process Training Course on QPE/QPF Manila 2012 64 How to promote the ETQPF performance in advance? How to include the uncertainty of the track forecast?
• Improve the ensemble prediction system. • Develop a Human–Machine Interaction information system to maximally handle the forecast uncertainty and perform the risk assessment for decision making.
Training Course on QPE/QPF Manila 2012 65 • Apply certain criterion to filter out the undesired model forecast based on • The selection radius, model ID, model resolution, initial time, forecast time, rainfall intensity, typhoon strength, moving speed, vertical shear, …
To optimally estimate the typhoon QPF
Training Course on QPE/QPF Manila 2012 66 ETQPFs
A user friendly GUI interface to manually modify the forecast typhoon track
To include the uncertainty of the track forecast
Training Course on QPE/QPF Manila 2012 67 Optimal ETQPF guide
Track improve
1. Best track 2. CWB official track 3. WRF ensemble track 4. Man Machine Mix
Training Course on QPE/QPF Manila 2012 68 ETQPF based on the track uncertanity
Training Course on QPE/QPF Manila 2012 69 ETQPF Summary An Ensemble based Typhoon QPF (ETQPF) technique is able to provide trustable model QPF products over Taiwan island. Ensemble mean, standard deviation, probabilistic QPF products based on the prior estimate of the typhoon track forecast from official or ensemble forecast. A Human–Machine Interaction information system (ETQPFs) has been implemented to maximally handle the forecast uncertainty and perform the risk assessment for decision making. Taking the advantage of information technology is essential and has great potential to well handle the forecast guidance with large uncertainty.
Training Course on QPE/QPF Manila 2012 70 The track of Typhoon ( location and speed) The Structure ( Intensity and Size) The Mesoscale convection ( Instability )
Training Course on QPE/QPF Manila 2012 71 There’s always something you don’t know
Blue line is Forecast Track.
Red line is traced by radar. It took longer and closer than Synoptic view.’
Nowcasting is crucial.
Training Course on QPE/QPF Manila 2012 72 QPF is highly depend on QPE for two purposes, one is for climatological, one is for the calibration.
Both climatology and dynamic model are necessary for the approach of QPF.
QPF is strongly depend on the process of Man Machine Mix for two reasons, one is that guidance is too complicated, the other one is that the it need perform in precisely and efficiently.
Training Course on QPE/QPF Manila 2012 73 Training Course on QPE/QPF Manila 2012 74