In snow-fed catchments, it is crucial to monitor and model the snow water equivalent (SWE), particularly when simulating the melt water runoff. SWE distribution can, however, be highly heterogeneous, particularly in forested environments. Within these locations, scant studies have explored the spatiotemporal variability in SWE in relation with vegetation characteristics, with only few successful attempts. The aim of this paper is to fill this knowledge gap, through a detailed monitoring at nine locations within a 3.49 km 2 forested catchment in southern Québec, Canada (47 N, 71 W). The catchment receives an annual average of 633 mm of solid precipitation and is predominantly covered with balsam fir stands. Extracted from intensive field campaign and high-resolution LiDAR data, this study explores the effect of fine scale forest features (tree height, tree diameter, canopy density, leaf area index [LAI], tree density and gap fraction) on the spatiotemporal variability in the SWE distribution. A nested stratified random sampling design was adopted to quantify small-scale variability across the catchment and 1810 manual snow samples were collected throughout the consecutive winters of 2016-17 and 2017-18. This study explored the variability of SWE using coefficients of variation (CV) and relating to the LAI. We also present existing spatiotemporal differences in maximum snow depth across different stands and its relationship with average tree diameter. Furthermore, exploiting key vegetation characteristics, this paper explores different approaches to model SWE, such as multiple linear regression, binary regression tree and neural networks (NN).We were unable to establish any relationship between the CV of SWE and the LAI.However, we observed an increase in maximum snow depth with decreasing tree diameter, suggesting an association between these variables. NN modelling (Nash-Sutcliffe efficiency [NSE] = 0.71) revealed that, snow depth, snowpack age and forest characteristics (tree diameter and tree density) are key controlling variables on SWE.Using only variables that are deemed to be more readily available (snow depth, tree height, snowpack age and elevation), NN performance falls by only 7% (NSE = 0.66).

In the central part of the Hindu Kush Himalayan region, snowmelt is one of the main inputs that ensures the availability of surface water outside the monsoon period. A common approach for snowpack modeling is based on the degree day factor (DDF) method to represent the snowmelt rate. However, the important seasonal variability of the snow processes is usually not represented when using a DDF method, which can lead to large uncertainties for snowpack simulation. The SPOT-VGT and the MODIS-Terra sensors provide valuable information for snow detection over several years. The aim of this work was to use those data to parametrize the seasonal variability of the snow processes in the hydrological distributed snow model (HDSM), based on a DDF method. The satellite products were corrected and combined in order to implement a database of 8 day snow cover area (SCA) maps over the northern part of the Dudh Koshi watershed (Nepal) for the period 1998–2017. A revisited version of the snow module of the HDSM model was implemented so as to split it into two parameterizations depending on the seasonality. Corrected 8 day SCA maps retrieved from MODIS-Terra were used to calibrate the seasonal parameterization, through a stochastic method, over the period of study (2013–2016). The results demonstrate that the seasonal parameterization reduces the error in the simulated SCA and increases the correlation with the MODIS SCA. The two-set version of the model improved the yearly RMSE from 5.9% to 7.7% depending on the basin, compared to the one-set version. The correlation between the model and MODIS passes from 0.73 to 0.79 in winter for the larger basin, Phakding. This study shows that the use of a remote sensing product can improve the parameterization of the seasonal dynamics of snow processes in a model based on a DDF method.