590 likes | 943 Views
Ensemble Forecasting. Yuejian Zhu Environmental Modeling Center December 6 th 2011 Acknowledgment for: Members of Ensemble & Probabilistic Guidance Team. Outlines. Responsibility of ensemble team Include all ensemble systems Current available data and products
E N D
Ensemble Forecasting Yuejian Zhu Environmental Modeling Center December 6th 2011 Acknowledgment for: Members of Ensemble & Probabilistic Guidance Team
Outlines • Responsibility of ensemble team • Include all ensemble systems • Current available data and products • Data access through all available resource • Digital probabilistic products • Web-based products • Implementations for past year • Include all major and minor • Update information from other centers • On going implementations • Include all major and minor • Reforecasting for future ensemble post process • Future plans • Development of major system • Post process and calibration • How to satisfy user requests?
Responsibilities of Ensemble Team - Assess, model, communicate uncertainty in numerical forecasts • Present uncertainty in numerical forecasting • Tasks • Design, implement, maintain, and continuously improve ensemble systems • Topics • Initial value related uncertainty • Model related forecast uncertainty • Ensemble systems • Global – GEFS / NAEFS / NUOPC • Regional – SREF / HREF / NARRE-TL / HWAF ensemble • Climate – Contributions to future coupling CFS configuration • NAEFS/GEFS downscaled • Ocean wave ensemble (MMA/EMC) • Statistical correction of ensemble forecasts • Tasks • Correct for systematic errors on model grid • Downscale information to fine resolution grid (NDFD) • Combine all forecast info into single ensemble/probabilistic guidance • Probabilistic product generation / user applications • Contribute to design of probabilistic products • Support use of ensembles by • Internal users (NCEP Service Center, WFOs, OHD/RFC forecasters and et al.) • External users (research, development, and applications)
Web-based Probabilistic Products • NCEP web-site • NCO supported all ensemble probabilistic products (main page) • EMC (experimental) developed probabilistic products • PQPF for various thresholds (include precipitation types) • RMOP for 500hPa height • Tropical storm track forecast • Extra-tropical storm track forecast (global tracking) • NAEFS web-site • CPC developed extended probabilistic forecast • Temperature, 500hPa height and precipitation • NCO supported NAEFS products • Spaghetti, • EMC (experimental) NAEFS products • PQPF for various thresholds • Anomaly forecast • CMC supported NAEFS products • Metagrams for NA major cities (example) • NUOPC web-site • EMC (experimental) NUOPC products • PQPF for various thresholds (include precipitation types) • TS track forecast for multi-model ensembles • FNMOC’s NUE probabilistic products (example)
Implementations for past year • Alaska downscaling - Dec. 7th 2010 • 6km probabilistic guidance for Alaska region (improvement for Max/Min temperature, wind speed/direction) • NUOPC – IOC Jan. 18 2011 • Getting FNMOC ensemble data (both raw and bias corrected) in NCEP and NOMADS for public • Ceremony of NUOPC IOC • NAEFS products upgrade – March 1st 2011 • Adding more variables for exchange, more variables for bias correction • 5th Ensemble User Workshop – May 9-11 2011 • In Laurel Maryland (highlight) • NAEFS data upgrade – May 24th 2011 • Receive CMC’s grib2 bias corrected forecast directly; • CCPA upgrade – July 26th 2011 • Adding 3-hrly analysisfor regional application • CMC’s GEFS upgrade – August 17th 2011 (highlight) • FNMOC’s GEFS upgrade – September 14th 2011 (highlight,stats) • NMME (August) and IMME (December) implementation in CPC • CMC’s GEFS data on NOMADS – November 2011
On going implementations • Major GEFS upgrade – Jan. 2012 • Highlight the changes • What do we expect from this change? • Overall skills, TC tracks, bias correction • More stats: http://www.emc.ncep.noaa.gov/gmb/yzhu/html/imp/201109_imp.html • NAEFS product upgrade – April. 2012 • Extend variables for CONUS downscaling • Temperature, winds, humidity and dew point • Precipitation calibrations (CONUS, RFC) • Anomaly forecast from CFSRR climatology • 6th NAEFS workshop – May 1-3 2012 • At Monterey CA – FNMOC will be a localhost • GEFS reforecast – August 2010 – • Leading by Tom Hamill (ESRL) • Bias correction for TS forecast • CCPA update • Use more historical data for calibration • Additional supports • CSTAR (the Collaborative Science, Technology, and Applied Research) program - ensemble sensitivity analysis
Future plans • GEFS • Plan for FY2013-2015 (EMC-ESRL collaboration) • NAEFS • Increasing products resolution and output frequencies • Improving calibration method • Apply GEFS reforecast for calibration • EMC-MDL collaboration for high moment adjustment • EMC-OHD collaboration for improving temperature and precipitation uncertainty forecast • NUOPC – adding FNMOC ensemble to NAEFS • Project information and highlights • Evaluation metrics • Tentative Schedule • Extended range forecast • Coupling ensemble system • Extended to 45 days to cover week-3 and week-4 • Seamless forecast system • Exchange extended forecast data with CMC • Support IMME and NMME to improve MJO prediction • User requests • Mean sea level pressure – more useful for users • Relative humidity or surface dew-point – more useful for users • Anomaly forecast – extreme weather index • Clustering (questionable for single model ensemble?)
http://mag.ncep.noaa.gov/NCOMAGWEB/appcontroller Global ensemble NAEFS Regional ensemble Return
http://www.cpc.ncep.noaa.gov/products/predictions/short_range/NAEFS/Outlook_D264.00.phphttp://www.cpc.ncep.noaa.gov/products/predictions/short_range/NAEFS/Outlook_D264.00.php Example of temperature forecast Upper tercile Return Lower tercile
Map of PQPF and Precipitation Types: every 6 hours, 4 different thresholds Return
Example: The tracker puts out 4 time a day for all cyclones (Northern Hemisphere) Return
Application for Alaska region and HPC Alaska desk 10-m U Max Temp Solid – RMS error Dash - spread Bias (absolute value) 10-m wind speed Max Temp Bias (absolute value) CRPS – small is better Return
5th NCEP Ensemble User Workshop • Logistics • Workshop organized by EMC/NCEP and DTC/NCAR (co-organizer) • May 10-12 2011, Laurel, MD, 90+ participants • NWS Regions (6), Headquarters (17), NCEP (44) • OAR (5), other government agencies (4), private (2), academic (5) & international (11) • For further info, see: http://www.dtcenter.org/events/workshops11/det_11/ • Main Theme • How to support NWS in its transition from single value to probabilistic forecasting • Goal is to convey forecast uncertainty in user relevant form • 46 presentations • Covering all ensemble forecast systems • SREF, GEFS/NAEFS, Wave ensemble, CFS and NMME • Reports from NCEP Service Centres and Regions (WFOs) • E.g., first numerical ensemble-based 2-day tornado, week 3-4, monthly MJO outlook • Working groups • Ensemble configurations - Ensemble forecasting • Statistic post processing - Reforecast/hindcast generation • Probabilistic product generation - Forecaster’s role and training • Ensemble data depository / access - Database interrogation / forecaster tools • Outcome / Recommendations • Prepared report for NWS roadmap reference • Plan for immediate steps (interim solution to be implemented in 2-3 years) • Outline for long term solution and resource requirements (5-10 years) • All activities to be coordinated under NWS Forecast Uncertainty Program (NFUSE) Return
CMC’s GEFS Implementation • Modification to EnKFanalysis configuration • Use 192 ensemble members instead of 96 • Has more satellite data • Upgrade GEM version • Use 4.2 version (vertical staggering) instead of version 3.0 • Increase model top to 2 hPa from 10 hPa • Resolutions • Horizontal: 600x300 (66km) from 400x200 (100km) • Vertical: 40 levels from 28 levels Return
FNOMC’s GEFS Implementation and Plan • 9 latitude band Ensemble Transform initialization instead of 5-banded • T159L42 instead of T119L30 in horizontal and vertical resolutions • Plan for T239L42 in operations for June 2012 • Implement the bias correction for selected variables (NAEFS algorithm) • Implement forecast vs observation verification system • NUOPC (FNMOC+NCEP+CMC) products Return
Tropical Cyclone Track Error T119, T159, T239 Homogenous NHTC track forecast error (km), for G119, G159, and G239 ensemble mean tracks as denoted in key. Also shown is the average forecast error of the T239L30 NOGAPS operational deterministic forecast (CTL). The numbers of verifying forecasts are shown below the x axis. The differences between G119 and G159 are statistically significant at the 95% level out to 96 h. The differences between G159 and G239 are statistically significant at 48 and 72 h. All used global ET. (Fig 3 from Impact of Resolution and Design on the U.S. Navy Global Ensemble Performance in the Tropics, Reynolds, et al., MWR, July 2011, p 2145-2155.) Return
Proposal Changes Model and initialization Using GFS V9.01 (current operational GFS) instead of GFS V8.00 Improved Ensemble Transform with Rescaling (ETR) initialization Improved Stochastic Total Tendency Perturbation (STTP) Configurations T254 (55km) horizontal resolution for 0-192 hours (from T190 – 70km) T190 (70km horizontal resolution for 192-384 hours (same as current opr) L42 vertical levels for 0-384 hours (from L28) Add Sunshine duration for TIGGE data exchange Part of products will be delayed by approximately 20 minutes Due to limit CCS resources 40-42 nodes for 70 minutes (start +4:35 end: +5:45) Unchanged: 20+1 members per cycle, 4 cycles per day pgrb file output at 1*1 degree every 6 hours GEFS and NAEFS post process output data format Why do we make this configurations? Considering the limited resources and resolution makes difference What do we expect from this implementation? Improve general probabilistic forecast skill overall Significant improvement of tropical storm tracks (especially for Atlantic basin) Return
Anomaly Correlation Winter 2 months 11.00d Skillful line 10.25d SH 500hPa height NH 500hPa height GFS V8.0 .vs V9.0 NH 850hPa temperature SH 850hPa temperature
Atlantic, AL01~19 (06/01~11/30/2011) GEFSx runs once per day before Oct. GEFSo---GEFS T190 (operational run) GEFSx---GEFS T254 (parallel run) GFS ------GFS T574 (operational run) 11% Improvement 20% Track error(NM) 22% 12% #CASES 309 279 251 227 202 162 125 88 Return Forecast hours
CRPS for NH 850hPa Temperature CRPS for NH 10m U-wind CRPS for NH 2-meter Temperature CRPS for NH 10m V-wind
Tmax Tmin CRPS CRPS • Latest evaluation for CONUS temperature forecast by apply : • Bias correction at 1*1 degree for NCEP GFS/GEFS, CMC/GEFS • Bybrid bias corrected NCEP GFS and GEFS • Apply statistical downscaling for all bias corrected forecast • Combined all forecasts at 5*5 km (NDGD) grid with adjustment - NAEFS Temperature CRPS Return
V 10m U 10m CRPS CRPS Wind speed Wind direction CRPS CRPS Return
Dew point T RH RMS & Spread RMS & Spread Dew Point T RH CRPS CRPS Return
Precipitation calibration for 2009-2010 winter season (CONUS only) Comparison for GFS and ensemble control (raw and bias corrected) ETS for all lead-time ETS for 0-6hr fcst BIAS for all lead-time BIAS for 0-6hr fcst Perfect bias = 1.0 Return Courtesy of Yan Luo The probabilistic scores (CRPS -not show here) is much improved as well. We are still working on the different weights, different RFC regions, downscaled to 5km as well. More results will come in soon. Plan for implementation: Q4FY11
Significantly reduced bias for CONUS and each RFC CONUS NWRFC NCRFC MBRFC NERFC • Consistently effective along with leading time MARFC CNRFC OHRFC Bias=1.0 CBRFC ABRFC LMRFC SERFC WGRFC • More effective on lower amount precip 1 * 1 deg CBRFC OHRFC before bias correction before bias correction after bias correction 1 * 1 deg after bias correction 1 * 1 deg
Experimental maps to support CSTAR program for winter season
GEFS Implementation Plan for FY13-15 Hybrid data assimilation based GEFS initializations Using 6-hr EnKF forecast to combine improved ETR (without cycling) (Schematic diagram) Improving ETR (still in discussion and investigation) Adaptive modification of initial and stochastic model perturbation variances Based on recursive average monitoring of forecast errors and ensemble spread Avoid having to tune perturbation size after each analysis/model/ensemble changes Improving performance and easy maintenance Real time generation of hind-casts (pending on resource) Make control forecast once every ~5th day (6 runs for each cycle) T254L42 (0-192) and T190L42 (192-384), and use new reanalysis (~30y) Increasing sample of analysis – forecast pairs for statistic corrections Improving bias correction beyond 5-d Potential for regime/situation dependent bias correction Coupled ocean-land-atmosphere ensemble Couple MOM4/HYCOM with land-atmosphere component using ESMF Depending on skill, extend integration to 35 days Merge forecasts with CFS ensemble for seamless weather climate interface Land perturbation and surface perturbation (later) Explore predictability in intra-seasonal time scale Potential skill beyond 15 days Hydro-meteorological (river flow) ensemble forecasting Pending on operational LDAS/GLDAS, and RFC application Return
Development of Statistical Post-Processing for NAEFS • Opportunities for improving the post-processor • Utilization of additional input information • More ensemble, high resolution control forecasts (hybrid?) • Using reforecast information to improve week-2 forecast and precipitation • Analysis field (such as RTMA and etc..) • Improving calibration technique • Calibration of higher moments (especially spread) • Use of objective weighting in input fields combination • Processing of additional variables with non-Gaussian distribution • Improve downscaling methods Future Configuration of EMC Ensemble Post-Processor Return
Project Information and Highlights • Evaluate the value added for current NAEFS inclusion of FNMOC ensembles • Current NAEFS products attached • Period: December 1st 2011 – May 31st 2012 • Cover winter and spring seasons • Available data for each participating centers • NCEP • Raw and bias corrected NCEP, CMC and FNMOC ensembles • CMC • Raw and bias corrected NCEP, CMC and FNMOC ensembles • FNMOC • Raw NCEP, CMC and FNMOC ensemble data only • AFWA • Raw NCEP, CMC and FNMOC ensemble data only Return
Evaluation metrics • NUOPC evaluation metrics • RMSE (MAE) and spread of ensemble mean, CRPS, Brier score for selected thresholds • Targeting evaluation parameters: • 2m T, 10m winds, 250hPa winds, 700hPa RH, 500hPa Z (targeting for next NCEP implementation) • TS tracks, precipitation (potential for future consideration) • Significant wave high, total cloud cover • NCEP is expected to have most evaluations against own analysis (GSI and RTMA) and observations (to connect with NCEP users) • NAEFS evaluation metrics • RMSE and spread of ensemble mean, CRPS (resolution and reliability) and etc… • Targeting evaluation parameters: • 2m T, 10m winds, 250 and 850 winds, 500hPa, 1000hPa Z, 850 T. (Targeting for next NCEP implementation) • Total precipitation (raw) Return
Scheduling • Current – November 30: preparation • December 1st 2011 – May 31st 2012: Data collection and evaluation for all participate centers • Mid-term performance review through NAEFS and UEO meetings • Another performance review by NAEFS workshop (May 1-3, Monterey CA) • Late of June • Fully evaluations from all participated centers • Possibility to have one day meeting at Silver Springs/Camp Springs or NUOPC workshop in June (???) • Mid of July • Decision will be made to recommend for NCEP implementation (or not) • July 1st – deadline for EMC RFCs for implementation • August – September 2012: NCO test and real time parallel for NCEP users evaluations (see additional slide for details) • September 25 2012: targeting for implementation Return
Future seamless forecast system NCEP/GEFS will plan for T254L42 (2011 GFS version) resolution with tuned ETR initial perturbations and adjusted STTP scheme for 21 ensemble members, forecast out to 16 days and 4 cycles per day. Extended to 45 days at T126L28/42 resolution, 00UTC only (coupling is still a issue?) NAEFS will include FNMOC ensemble in 2011, with improving post process which include bias correction, dual resolution and down scaling Main products: ENSO predictions??? Seasonal forecast??? Main event MJO GEFS/NAEFS service CFS service week-1 week-2 one month Weather/Climate linkage SEAMLESS • Main products: • Probabilistic forecasts for every 6-hr out to 16 days, 4 times per day: 10%, 50%, 90%, ensemble mean, mode and spread. • D6-10, week-2 temperature and precipitation probabilistic mean forecasts for above, below normal and normal forecast • MJO forecast (week 3 & 4 … ) Operational CFS has been implemented in Q2FY2011 with T126L64 atmospheric model resolution (CFSv2, 2010version) which is fully coupled with land, ocean and atmosphere (GFS+MOM4+NOAH), 4 members per day (using CFS reanalysis as initial conditions, one day older?), integrate out to 9 months. Future: initial perturbed CFS
Flow Chart for Hybrid Variation and Ensemble Data Assimilation System (HVEDAS) - concept Lower resolution EnKF assimilation t=j EnKF assimilation t=j+1 Ensemble fcst (1) t=j-1 j Ensemble fcst t=j, j+1 Ensemble initialization Ensemble fcst (2) t=j 16 days Replace Ensemble Mean Estimated Background Error Covariance from Ensemble Forecast (6 hours) Estimated Background Error Covariance from Ensemble Forecast (6 hours) Two-way hybrid GSI/3DVAR t=j GSI/3DVAR t=j+1 Hybrid Analysis? Higher resolution
Atlantic, AL01~17 (06/01~09/30/2011) GEFSo --- GEFS T190 (operational run) GEFSx --- GEFS T254 (parallel run) GFS ------ GFS T574 (operational run) 17% Improvement 24% Track error(NM) 26% 13% #CASES 235 213 194 178 159 133 103 75 Forecast hours
NAEFS downscaling parameters and products Last update: May 1st 2010 (NDGD resolutions) All products at 1*1 (lat/lon) degree globally Ensemble mean, spread, 10%, 50%, 90% and mode back
NEXT NAEFS pgrba_bc files (bias correction) back
2-m temp 10/90 probability forecast verificationNorthern Hem, period of Dec. 2007 – Feb. 2008 90% 3-month verifications 10% Top: 2-m temperature probabilistic forecast (10% and 90%) verification red: perfect, blue: raw, green: NAEFS Left: example of probabilistic forecasts (meteogram) for Washington DC, every 6-hr out to 16 days from 2008042300 back
Optimal Threshold = 15% Decision Theory Example Forecast? YES NO Critical Event: sfc winds > 50kt Cost (of protecting): $150K Loss (if damage ): $1M Hit False Alarm Miss Correct Rejection YES NO $150K $1000K Observed? $150K $0K back
Overall temperature forecasts: Average over past 30 days: (20080929-20081028) MAE Bias >10 err <3 err off. rank Best G. 2nd G. Worst G. 1 12-hr 2.44 0.7 0.1% 67.3% 1 out of 7 NAM40 65.4% NAM12 60.1% NGM80 44.4% 2 24-hr 2.84 1.0 0.3% 59.1% 2 out of 7 NAM40 60.3% NAM12 56.9% SREF 47.0% 3 36-hr 2.94 0.8 0.3% 57.8% 1 out of 7 NAM40 55.9% NAM12 52.6% NGM80 44.0% 4 48-hr 3.36 1.6 2.1% 52.8% 1 out of 7 MOSGd 48.9% NAM40 48.3% NGM80 12.9% 5 60-hr 3.26 1.0 1.7% 54.8% 1 out of 6 MOSGd 50.1% NAM12 48.8% NAM40 6.2% 6 72-hr 3.35 1.3 2.1% 53.1% 1 out of 5 MOSGd 49.9% NAM12 49.5% SREF 44.0% 7 84-hr 3.80 0.6 4.7% 49.0% 1 out of 5 NAEFS 48.6% SREF 44.5% NAM12 2.6% 8 96-hr 3.96 0.7 4.0% 44.4% 2 out of 4 NAEFS 46.2% HPCGd 42.6% MOSGd 40.6% 9 108-hr 4.43 0.9 5.5% 38.5% 2 out of 3 NAEFS 41.7% MOSGd 37.7% MOSGd 37.7% 10 120-hr 4.57 1.0 5.9% 36.6% 2 out of 4 NAEFS 40.9% HPCGd 36.5% MOSGd 36.3% 11 132-hr 4.83 0.7 7.8% 34.7% 1 out of 3 NAEFS 34.5% MOSGd 34.4% MOSGd 34.4% 12 144-hr 4.83 0.5 7.4% 34.7% 3 out of 4 HPCGd 36.4% NAEFS 35.5% MOSGd 33.3% 13 156-hr 5.43 0.1 11.9% 30.3% 3 out of 3 NAEFS 32.1% MOSGd 30.8% MOSGd 30.8% 14 168-hr 5.74 0.3 14.4% 27.7% 2 out of 4 HPCGd 27.7% MOSGd 26.9% NAEFS 26.1 Minimum temperature forecast: Average over past 30 days: (20080929-20081028) 1 12-hr 3.17 -1.2 1.0% 53.4% 3 out of 7 NAEFS 59.7% SREF 57.1% NGM80 21.8% 2 24-hr 3.03 -0.9 0.6% 55.5% 2 out of 7 SREF 57.2% NAEFS 54.2% NGM80 24.9% 3 36-hr 3.25 -0.8 0.9% 51.6% 3 out of 7 NAEFS 54.2% SREF 53.9% NGM80 23.2% 4 48-hr 3.94 -1.1 2.9% 43.2% 3 out of 7 NAEFS 51.9% SREF 45.8% NGM80 6.2% 5 60-hr 4.30 -0.4 4.4% 39.1% 4 out of 6 NAEFS 49.2% SREF 43.0% NAM40 8.9% 6 72-hr 4.76 0.1 6.4% 33.7% 5 out of 5 NAEFS 42.9% SREF 40.1% NAM12 35.2% 7 84-hr 4.85 0.3 7.5% 34.7% 2 out of 6 NAEFS 40.0% MOSGd 33.4% NAM12 8.9% 8 96-hr 5.24 0.4 13.0% 33.1% 1 out of 3 NAEFS 32.7% MOSGd 29.9% MOSGd 29.9% 9 108-hr 5.11 0.8 12.8% 35.4% 1 out of 4 HPCGd 34.5% NAEFS 32.1% MOSGd 30.5% 10 120-hr 5.31 0.7 12.0% 31.9% 1 out of 3 MOSGd 31.6% NAEFS 24.8% NAEFS 24.8% 11 132-hr 4.97 0.7 9.9% 35.1% 2 out of 4 HPCGd 38.0% MOSGd 30.9% NAEFS 27.2% 12 144-hr 5.42 0.6 15.0% 35.0% 1 out of 3 MOSGd 31.3% NAEFS 29.0% NAEFS 29.0% 13 156-hr 5.40 0.5 14.9% 35.7% 1 out of 4 HPCGd 32.9% MOSGd 32.7% NAEFS 23.4% 14 168-hr 5.46 1.1 17.7% 38.1% 1 out of 3 MOSGd 35.6% NAEFS 28.4% NAEFS 28.4% Official Guidance: NGM80, NAM40, SREF, NAM12, MOSGd, HPCGd, NAEFS Contributed by Richard Grumm (WFO)
Ocean Wave Ensemble System- Hendrik Tolman • Configuration of ocean wave ensemble system Wave ensemble has been running since 2008 • Running on 1°×1° wave model grid as the control. • 20 wave members generated through GEFS using ETR method • Cycling initial conditions for individual members to introduce uncertainty in swell results. • 10 day forecast using the GEFS bias corrected 10m wind (future operation) • Improving forecast uncertainties through • Introducing ensemble initial perturbations from previous model cycle • Introducing bias corrected ensembles as external forcing. • Example of comparison (wave heights) • Plans • Work towards a combined NCEP-FNMOC ensemble • Analyze the role of swell played in the wave ensemble
Comparison of the ensemble systems old cycle Old ensemble setup, ensemble with cycling of initial conditions and wind bias correction (BC). Mean wave height (contours) and spread (shading) 2008/03/28 t06z 120h forecast cycle, BC back
50th (median) and mean are best back Courtesy of Dave Novak
EMC-MDL COLLABORATION • Compare quality of current operational / experimental products • Gridded MOS vs. Downscaled NAEFS • Ongoing • Kathy Gilbert, Val Dragostano – Zoltan Toth, Bo Cui, Yuejian Zhu • Proxy for truth issue unresolved • Need observations independent of MOS • MDL experimental ensemble guidance vs. Downscaled NAEFS • 10/50/90 percentiles to be evaluated • Matt Peroutka & Zoltan Toth • Proxy for truth issue • Proxy for truth? • Agree on best proxy for truth • Collaborate on • Improving RTMA, including bias correction for FG • Creating best CONUS precipitation analysis & archive • Joint research into best downscaling methods? • Climate, regime, case dependent methods • Addition of fine temporal/spatial variability into ensemble
SREF Probability of STP Ingredients: Time Trends 24 hr SREF Forecast Valid 21 UTC 7 April 2006 Prob (MLCAPE > 1000 Jkg-1) X Prob (6 km Shear > 40 kt) X Prob (0-1 km SRH > 100 m2s-2) X Prob (MLLCL < 1000 m) X Prob (3h conv. Pcpn > 0.01 in) Shaded Area Prob > 5% Max 50%