Weather forecasts and climate projections of precipitation phase and type in winter storms are challenging due to the complicated underlying microphysical and dynamical processes. In the Canadian numerical weather prediction model, explicit freezing rain (FR) at the surface is often overestimated during the winter season for situations in which snow is observed. For a case study simulated using this model with the Predicted Particle Properties (P3) microphysics scheme, the secondary ice production (SIP) process has a major impact on the surface precipitation type. Parameterized SIP substantially reduces FR due to increased collection of supercooled drops with ice particles formed by rime splintering. Hindcast simulations of 40 winter cases show that these results are systematic, and the decreased frequency of FR leads to improved forecast skill relative to observations. Thus, accounting for SIP in the model is critical for accurately simulating precipitation types.