Abstract: This study leverages observations from NYSM to evaluate and improve representation of snow within the NWM and its associated land surface model. Distributed NWM simulations were ran and analyzed, forced by gridded meteorological analyses, and Noah-MP point simulations, forced by NYSM observations. Distributed NWM runs, with a baseline configuration, show substantial SWE biases caused by biases in meteorological forcing used, imperfect representation of snow processes, and mismatches between land cover in the model and NYSM station locations. Noah-MP point simulations, using baseline configuration, reveal a systematic positive bias in SWE accumulation. Noah-MP point simulations, with improved precipitation phase partitioning, reveal a systematic negative bias in SWE ablation rates. Sensitivity experiments highlight uncertain parameters within Noah-MP that strongly affect ablation rates and show particularly large sensitivity to snow albedo decay time-scale parameter, which modulates snow albedo decay rates. Distributed NWM experiments, with precipitation phase partitioning and TAU0 adjusted based on Noah-MP point simulation results, show qualitatively similar sensitivities. However, the distributed experiments do not show clear improvements when compared to SWE and streamflow observations. This is likely due to some combination of sources of bias in the baseline-distributed run and biases in other parameterized processes unrelated to snow in the NWM.