Evaluation of the Experimental Climate Prediction Center’s fire danger forecasts with remote automated weather station observations
Hauss J. Reinbold A C , John O. Roads B and Timothy J. Brown AA Desert Research Institute, 2215 Raggio Parkway, Reno, NV 89512-1095, USA.
B Scripps Experimental Climate Prediction Center, University of California San Diego, 0224 La Jolla, CA 92093, USA.
C Corresponding author. Telephone: +1 775 673 7386; fax: +1 775 674 7007; email: hauss.reinbold@dri.edu
International Journal of Wildland Fire 14(1) 19-36 https://doi.org/10.1071/WF04042
Submitted: 24 August 2004 Accepted: 21 December 2004 Published: 7 March 2005
Abstract
The Scripps Experimental Climate Prediction Center has been routinely making regional forecasts of atmospheric elements and fire danger indices since 27 September 1997. This study evaluates these forecasts using selected remote automated weather station observations over the western USA. Bias and anomaly correlations are computed for daily 2-m maximum, minimum, average temperature, 2-m maximum, minimum and average relative humidity, precipitation and afternoon 10-m wind speed, and four National Fire Danger Rating System indices—ignition component, spread component, burning index and energy release component. Of the atmospheric elements, temperature generally correlates well, but relative humidity, precipitation and wind speed are less correlated. Fire danger indices have much lower correlations, but do show useful spatial structure in some areas such as Southern California, Arizona and Nevada.
Anderson BT , Roads J (2002) Regional simulation and of summertime precipitation over the Southwestern United States. Journal of Climate 15, 3321–3342.
| Crossref | GoogleScholarGoogle Scholar |
Chen SC, Roads J, Juang H , Kanamitsu M (1999) Global to regional simulation of California’s wintertime precipitation. Journal of Geophysical Research 104, 31 517–31 532.
Han J , Roads J (2004) US climate sensitivity simulated with the NCEP Regional Spectral Model. Climatic Change 62, 115–154.
| Crossref | GoogleScholarGoogle Scholar |
Juang H , Kanamitsu M (1994) The NMC nested regional spectral model. Monthly Weather Review 122, 3–26.
| Crossref | GoogleScholarGoogle Scholar |
Roads JO (2004) Experimental weekly to seasonal, global to regional US precipitation forecasts. Journal of Hydrology 288, 153–169.
| Crossref | GoogleScholarGoogle Scholar |
Appendix 1. Statistical measures
The purpose of forecast verification is to determine the quantitative accuracy of the forecast. As described previously by Roads et al. (2005) the statistical methods employed in this study as a means of forecast verification were bias, root mean square error, anomaly correlation and standard deviation. Due to space constraints, the root-mean square error and standard deviation are not discussed here, but these measures can be found in Reinbold (2003).
Bias is a simple calculation of forecast minus observation (Wilks 1995), or
where f is the forecast value and o is the value of the observation. When shown in graphical form, this calculation has the benefit of revealing under what situations the model is over- or under-forecasting and by how much. It is also useful in determining potential seasonal characteristics in the errors between the datasets.
Anomaly correlations are commonly used to evaluate extended forecasts. This measure is designed to reflect good forecasts in the pattern of an observed field, but does not effectively measure the magnitude of the values (Wilks 1995). There are two different equations, representing the two types of anomaly correlation used in this study. The first is for judging the spatial variation and correlation of the anomalies (Roads 2004; Roads et al. 2005) and is not shown here. The second is better described as temporal variations in spatial correlations (equation 3; Roads et al. 2005). Anomalies are first computed by taking the difference between the total forecast (either by region or for the entire Western USA) and the climatological monthly means. In other words,
where A is the anomaly, f is the forecast and Cf is the climatological mean for that forecast type (weekly, monthly or seasonal mean).
Given that A is a forecast anomaly of any type (weekly, monthly or seasonal mean) and that B is the validating anomaly from observation, the temporal variations in the spatial anomaly correlations (AC, sometimes known as pattern correlation) are calculated using
where M is the total number of RAWS in the current region (M = 262 for the western USA). Missing values in the observational anomalies for equation (3) reduce the anomaly summations for both datasets by the number missing (the missing RAWS values and matching interpolated validation or forecast values are removed).