Abstract. Snow cover modeling remains a major challenge in climate and numerical weather prediction (NWP) models even in recent versions of high-resolution coupled surface–atmosphere (i.e., at kilometer scale) regional models. Evaluation of recent climate simulations, carried out as part of the WCRP-CORDEX Flagship Pilot Study on Convection (FPSCONV) with the CNRM-AROME convection-permitting regional climate model at 2.5 km horizontal resolution, has highlighted significant snow cover biases, severely limiting its potential in mountain regions. These biases, which are also found in AROME numerical weather prediction (NWP) model results, have multiple causes, involving atmospheric processes and their influence on input data to the land surface models in addition to deficiencies of the land surface model itself. Here we present improved configurations of the SURFEX-ISBA land surface model used in CNRM-AROME. We thoroughly evaluated these configurations on their ability to represent seasonal snow cover across the European Alps. Our evaluation was based on coupled simulations spanning the winters of 2018–2019 and 2019–2020, which were compared against remote sensing data and in situ observations. More specifically, the study tests the influence of various changes in the land surface configuration, such as the use of multi-layer soil and snow schemes, the division of the energy balance calculation by surface type within a grid cell (multiple patches), and new physiographic databases and parameter adjustments. Our findings indicate that using only more detailed individual components in the surface model did not improve the representation of snow cover due to limitations in the approach used to account for partial snow cover within a grid cell. These limitations are addressed in further configurations that highlight the importance, even at kilometer resolution, of taking into account the main subgrid surface heterogeneities and improving representations of interactions between fractional snow cover and vegetation. Ultimately, we introduce a land surface configuration, which substantially improves the representation of seasonal snow cover in the European Alps in coupled CNRM-AROME simulations. This holds promising potential for the use of such model configurations in climate simulations and numerical weather prediction both for AROME and other high-resolution climate models.