Physically based snow models provide valuable information on snow cover evolution and are therefore key to provide water availability projections. Yet, uncertainties related to snow modelling remain large as a result of differences in the representation of snow physics and meteorological forcing. While many studies focus on evaluating these uncertainties, no snow model comparison has been done in environments where sublimation is the main ablation process. This study evaluates a case study in the semi-arid Andes of Chile and aims to compare two snow models with different complexities, SNOWPACK and SnowModel, at a local point over one snow season and to evaluate their sensitivity relative to parameterisation and forcing. For that purpose, the two models are forced with (i) the most ideal set of input parameters, (ii) an ensemble of different physical parameterisations, and (iii) an ensemble of biased forcing. Results indicate large uncertainties depending on forcing, the snow roughness length z0, albedo parameterisation, and fresh snow density parameterisation. The uncertainty caused by the forcing is directly related to the bias chosen. Even though the models show significant differences in their physical complexity, the snow model choice is of least importance, as the sensitivity of both models to the forcing data was on the same order of magnitude and highly influenced by the precipitation uncertainties. The sublimation ratio ranges are in agreement for the two models: 36.4 % to 80.7 % for SnowModel and 36.3 % to 86.0 % for SNOWPACK, and are related to the albedo parameterisation and snow roughness length choice for the two models.

In the semiarid Andes of Chile, farmers and industry in the cordillera lowlands depend on water from snowmelt, as annual rainfall is insufficient to meet their needs. Despite the importance of snow cover for water resources in this region, understanding of snow depth distribution and snow mass balance is limited. Whilst the effect of wind on snow cover pattern distribution has been assessed, the relative importance of melt versus sublimation has only been studied at the point scale over one catchment. Analyzing relative ablation rates and evaluating uncertainties are critical for understanding snow depth sensitivity to variations in climate and simulating the evolution of the snowpack over a larger area and over time. Using a distributed snowpack model (SnowModel), this study aims to simulate melt and sublimation rates over the instrumented watershed of La Laguna (513 km², 3150-5630 a.s.l., 30° S, 70° W), during two hydrologically contrasting years (i.e., dry vs. wet). The model is calibrated and forced with meteorological data from nine Automatic Weather Stations (AWSs) located in the watershed and atmospheric simulation outputs from the Weather Research and Forecasting (WRF) model. Results of simulations indicate first a large uncertainty in sublimation-to-melt ratios depending on the forcing as the WRF data have a cold bias and overestimate precipitation in this region. These input differences cause a doubling of the sublimation-to-melt ratio using WRF forcing inputs compared to AWS. Therefore, the use of WRF model output in such environments must be carefully adjusted so as to reduce errors caused by inherent bias in the model data. For both input datasets, the simulations indicate a similar sublimation fraction for both study years, but ratios of sublimation to melt vary with elevation as melt rates decrease with elevation due to decreasing temperatures. Finally results indicate that snow persistence during the spring period decreases the ratio of sublimation due to higher melt rates.