Post Processing

1. 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. 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

6. 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.

7. 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.

10. Current “Operational” Gridded MOS

11. 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

14. 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.

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

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

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

21. 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.

22. 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?

23. 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).

24. National Weather Service locations used in the study.

25. 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

26. 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.

27. Temperature Comparisons

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

29. Large one-day temp changes MAE for each forecast type during periods of large temperature change (10F over 24-hr), 1 August 2003 – 1 August 2004. Includes data for all stations.

30. MAE for each forecast type during periods of large departure (20F) from daily climatological values, 1 August 2003 – 1 August 2004.

31. Number of days each forecast is the most accurate, all stations. In (a), tie situations are counted only when the most accurate temperatures are exactly equivalent. In (b), tie situations are cases when the most accurate temperatures are within 2F of each other. Looser Tie Definition

32. Number of days each forecast is the least accurate, all stations. In (a), tie situations are counted only when the least accurate temperatures are exactly equivalent. In (b), tie situations are cases when the least accurate temperatures are within 2F of each other. Looser Tie Definition

33. Highly correlated time series Time series of MAE of MAX-T for period one for all stations, 1 August 2003 – 1 August 2004. The mean temperature over all stations is shown with a dotted line. 3-day smoothing is performed on the data.

34. Coldspell Time series of bias in MAX-T for period one for all stations, 1 August 2003 – 1 August 2004. Mean temperature over all stations is shown with a dotted line. 3-day smoothing is performed on the data.

35. MAE for all stations, 1 August 2003 – 1 August 2004, sorted by geographic region. MOS Seems to have the most problems at high elevation stations.

36. Bias for all stations, 1 August 2003 – 1 August 2004, sorted by geographic region.

37. Precipitation Comparisons

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

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

40. Precipitation Brier Score for all stations, 1 August 2003 – 1 August 2004, sorted by geographic region.

41. Reliability diagrams for period 1 (a), period 2 (b), period 3 (c) and period 4 (d).

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

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

45. They don’t MOS!

46. UW Bias Correction of WRF

47. And many others…