AbstractAccurate estimation of snow accumulation and melt is a critical part of decision‐making in snow‐dominated watersheds. In this study, we demonstrate a flexible methodology to couple a detailed snow model, Crocus, separately to two different land surface models (LSMs), Noah‐MP and Noah. The original LSMs and the coupled models (Noah‐MP‐Crocus and Noah‐Crocus) are used to simulate snow depth, snow water equivalent, and other water and energy states and fluxes. The results of simulations are compared against a wide range of independent gridded and point scale reference data sets. Our results show that coupling the detailed snow model, Crocus, with the LSMs improves the snow depth and snow water equivalent relative to independent observations. Overall, larger improvements are obtained with coupling Crocus to the Noah LSM, with the coupled Noah‐Crocus configuration reducing the RMSE and bias of snow depth from 2% to 12% and 57% to 75%, respectively, relative to Snow Data Assimilation System (SNODAS) and snow product from the University of Arizona. On the other hand, smaller improvements are obtained by coupling Crocus with Noah‐MP. The Coupled Noah‐MP‐Crocus reduces the snow depth bias but slightly degrades the RMSE of snow depth and snow water equivalent. The corresponding impacts in other water budget terms such as evapotranspiration, soil moisture, and streamflow, however, are mixed, pointing to the significant need to improve the coupling assumptions of these processes within land models. Overall, the interoperable coupling framework demonstrated here offers the opportunity to include more detailed snow physics and processes, and to advance data assimilation systems through improved exploitation of information from snow remote sensing instruments.