Winter road maintenance is one of the main tasks for the Washington State Department of Transportation. Anti-icing, that is, the preemptive application of chemicals, is often used to keep the roadways free of ice. Given the preventive nature of anti-icing, accurate predictions of road ice are needed. Currently, anti-icing decisions are usually based on deterministic weather forecasts. However, the costs of the two kinds of errors are highly asymmetric because the cost of a road closure due to ice is much greater than that of taking anti-icing measures. As a result, probabilistic forecasts are needed to optimize decision making.We propose two methods for forecasting the probability of ice formation. Starting with deterministic numerical weather predictions, we model temperature and precipitation using distributions centered around the bias-corrected forecasts. This produces a joint predictive probability distribution of temperature and precipitation, which then yields the probability of ice formation, defined here as the occurrence of precipitation when the temperature is below freezing. The first method assumes that temperatures, as well as precipitation, at different spatial locations are conditionally independent given the numerical weather predictions. The second method models the spatial dependence between forecast errors at different locations. The model parameters are estimated using a Bayesian approach via Markov chain Monte Carlo.We evaluate both methods by comparing their probabilistic forecasts with observations of ice formation for Interstate Highway 90 in Washington State for the 2003-2004 and 2004-2005 winter seasons. The use of the probabilistic forecasts reduces costs by about 50% when compared to deterministic forecasts. The spatial method improves the reliability of the forecasts, but does not result in further cost reduction when compared to the first method.