Model Post Processing

1. Model Post Processing

2. Model Output Can Usually Be Improved with Post Processing • Can remove systematic bias • Can produce probabilistic information from deterministic information • Can provide forecasts for parameters that the model incapable of modeling successfully due to resolution or physics issues (e.g., shallow fog)

3. Post Processing • Model Output Statistics was the first post-processing method used by the NWS (1969) • Based on multiple linear regression. • Essentially unchanged in 40 years. • Does not consider non-linear relationships between predictors and predictands. • Does take out much of systematic bias.

4. Day 2 (30-h) GFS MOS Max Temp Equation for KSLC (Cool Season – 0000 UTC cycle)

5. Day 2 (42-h) GFS MOS Max Temp Equation for KUNV (Warm Season – 1200 UTC cycle)

6. Day 2 (30-h) GFS MOS Min Temp Equation for KDCA (Cool Season - 1200 UTC cycle) Develop / Evaluate

8. MOS Developed by and Run at the NWS Meteorological Development Lab (MDL) • Full range of products available at: http://www.nws.noaa.gov/mdl/synop/index.php

10. Global Ensemble MOS • Ensemble MOS forecasts are based on the 0000 UTC run of the GFS Global model ensemble system. These runs include the operational GFS, a control version of the GFS (run at lower resolution), and 20 additional runs. • Older operational GFS MOS prediction equations are applied to the output from each of the ensemble runs to produce 21 separate sets of alphanumeric bulletins in the same format as the operational MEX message.

11. Gridded MOS • The NWS needs MOS on a grid for many reasons, including for use in their IFPS analysis/forecasting system. • The problem is that MOS is only available at station locations. • To deal with this, NWS created Gridded MOS. • Takes MOS at individual stations and spreads it out based on proximity and height differences. Also does a topogaphic correction dependent on reasonable lapse rate.

14. Current “Operational” Gridded MOS

15. Localized Aviation MOS Program(LAMP) • Hourly updated statistical product • Like MOS but combines: • MOS guidance • the most recent surface observations • simple local models run hourly • GFS output

16. Practical Example of Solving a LAMP Temperature Equation Y = b + a1x1 + a2x2 + a3x3 + a4x4 Y = LAMP temperature forecast Equation Constant b = -6.99456 Predictor x1 = observed temperature at cycle issuance time (value 66.0) Predictor x2 = observed dewpoint at cycle issuance time (value 58.0) Predictor x3 = GFS MOS temperature (value 64.4) Predictor x4 = GFS MOS dewpoint (value 53.0)

17. Theoretical Model Forecast Performance of LAMP, MOS, and Persistence

19. Verification of LAMP 2-m Temperature Forecasts

21. MOS Performance • MOS significantly improves on the skill of model output. • National Weather Service verification statistics have shown a narrowing gap between human and MOS forecasts.

22. Cool Season Mi. Temp – 12 UTC Cycle Average Over 80 US stations

23. Prob. Of Precip.– Cool Season(0000/1200 UTC Cycles Combined)

24. MOS Won the Department Forecast Contest in 2003 For the First Time!

28. Average or Composite MOS • There has been some evidence that an average or consensus MOS is even more skillful than individual MOS output. • Vislocky and Fritsch (1997), using 1990-1992 data, found that an average of two or more MOS’s (CMOS) outperformed individual MOS’s and many human forecasters in a forecasting competition.

29. Some Questions • How does the current MOS performance…driven by far superior models… compare with NWS forecasters around the country. • How skillful is a composite MOS, particularly if one weights the members by past performance? • How does relative human/MOS performance vary by forecast projection, region, large one-day variation, or when conditions vary greatly from climatology? • Considering the results, what should be the role of human forecasters?

30. This Study • August 1 2003 – August 1 2004 (12 months). • 29 stations, all at major NWS Weather Forecast Office (WFO) sites. • Evaluated MOS predictions of maximum and minimum temperature, and probability of precipitation (POP).

31. National Weather Service locations used in the study.

32. Forecasts Evaluated • NWS Forecast by real, live humans • EMOS: Eta MOS • NMOS: NGM MOS • GMOS: GFS MOS • CMOS: Average of the above three MOSs • WMOS: Weighted MOS, each member is weighted by its performance during a previous training period (ranging from 10-30 days, depending on each station). • CMOS-GE: A simple average of the two best MOS forecasts: GMOS and EMOS

33. The Approach: Give the NWS the Advantage! • 08-10Z-issued forecast from NWS matched against previous 00Z forecast from models/MOS. • NWS has 00Z model data available, and has added advantage of watching conditions develop since 00Z. • Models of course can’t look at NWS, but NWS looks at models. • NWS Forecasts going out 48 (model out 60) hours, so in the analysis there are: • Two maximum temperatures (MAX-T), • Two minimum temperatures (MIN-T), and • Four 12-hr POP forecasts.

34. Temperature Comparisons

35. Temperature MAE (F) for the seven forecast types for all stations, all time periods, 1 August 2003 – 1 August 2004.

36. Precipitation Comparisons

37. Brier Scores for Precipitation for all stations for the entire study period.

38. Brier Score for all stations, 1 August 2003 – 1 August 2004. 3-day smoothing is performed on the data.

39. There are many other post-processing approaches • Neural nets Attempts to duplicate the complex interactions between neurons in the human brain.

41. They don’t do MOS!

42. Dynamic MOS Using Multiple Models • MOS equations are updated frequently, not static like the NWS. • Example: DiCast used by the Weather Channel

45. ForecastAdvisor.com

48. UW Bias Correction of WRF

49. Bayesian Model Averaging (BMA) • A good way to optimize the use of ensembles of forecasts to provide calibrated probabilistic guidance. • Can weight models and the variability by previous performance.