Abstract The ensemble predictability of the January 2015 and 2016 East Coast winter storms is assessed, with model precipitation forecasts verified against observational datasets. Skill scores and reliability diagrams indicate that the large ensemble spread produced by operational forecasts was warranted given the actual forecast errors imposed by practical predictability limits. For the 2015 storm, uncertainties along the western edge’s sharp precipitation gradient are linked to position errors of the coastal low, which are traced to the positioning of the preceding 500-hPa wave pattern using the ensemble sensitivity technique. Predictability horizon diagrams indicate the forecast lead time in terms of initial detection, emergence of a signal, and convergence of solutions for an event. For the 2016 storm, the synoptic setup was detected at least 6 days in advance by global ensembles, whereas the predictability of mesoscale features is limited to hours. Convection-permitting WRF ensemble forecasts downscaled from the GEFS resolve mesoscale snowbands and demonstrate sensitivity to synoptic and mesoscale ensemble perturbations, as evidenced by changes in location and timing. Several perturbation techniques are compared, with stochastic techniques [the stochastic kinetic energy backscatter scheme (SKEBS) and stochastically perturbed parameterization tendency (SPPT)] and multiphysics configurations improving performance of both the ensemble mean and spread over the baseline initial conditions/boundary conditions (IC/BC) perturbation run. This study demonstrates the importance of ensembles and convective-allowing models for forecasting and decision support for east coast winter storms.