Accurate estimation of snowfall episodes several hours or even days in advance is essential to minimize risks to transport and other human activities. Every year, these episodes cause severe traffic problems on the northwestern Iberian Peninsula. In order to analyze the influence of different parameterization schemes, 15 snowfall days were analyzed with the Weather Research and Forecasting (WRF) model, defining three nested domains with resolutions of 27, 9, and 3 km. We implemented four microphysical parameterizations (WRF Single‐Moment 6‐class scheme, Goddard, Thompson, and Morrison) and two planetary boundary layer schemes (Yonsei University and Mellor‐Yamada‐Janjic), yielding eight distinct combinations. To validate model estimates, a network of 97 precipitation gauges was used, together with dichotomous data of snowfall presence/absence from snowplow requests to the emergency service of Spain and observatories of the Spanish Meteorological Agency. The results indicate that the most accurate setting of WRF for the study area was that using the Thompson microphysical parameterization and Mellor‐Yamada‐Janjic scheme, although the Thompson and Yonsei University combination had greater accuracy in determining the temporal distribution of precipitation over 1 day. Combining the eight deterministic members in an ensemble average improved results considerably. Further, the root mean square difference decreased markedly using a multiple linear regression as postprocessing. In addition, our method was able to provide mean ensemble precipitation and maximum expected precipitation,which can be very useful in the management of water resources. Finally, we developed an application that allows determination of the risk of snowfall above a certain threshold.