High latitude regions, including the circumpolar boreal biome, are experiencing important changes in the availability of usable surface water because of climate change. In this context, an adequate representation of the land-atmosphere interaction is critical to ensure optimal management of current and future water resources, forest management, and climate prediction. However, the task is particularly intricate in high-latitude boreal forest, as land surface model faces several challenges due to the unique environmental conditions and ecological characteristics. The objective of this study is to quantify the impact of forest landscape heterogeneity, specifically stand leaf-area index (LAI), soil texture, and drainage regime, on surface water and energy balance in a small boreal high-latitude sub-catchment. To this end, hydrometeorological conditions at seventeen 20×20 m plots in a 1-km2 boreal forest sub-basin are simulated using the Canadian Land Surface Scheme (CLASS), a land surface model, at the point scale. The subplot-scale soil texture, drainage regime, and vegetation characteristics and type are based as closely as possible on field measurements and observations for the 17 plots. The model-driven experiment comprises two sets of simulations using CLASS, each employing the same model setup and run for the 17 experimental plots. The main set employs meteorological forcing from a local micrometeorological tower within the sub-basin to investigate the plot-to-plot variability of albedo, energy fluxes, and soil state variables. A second set of simulations is conducted using meteorological forcing from the ERA5-Land reanalysis, which spans from 1986 to 2022. This data provides a longer time series, enabling a more accurate representation of the interannual climatic variability in the sub-basin. The results of the main and secondary sets of CLASS simulations are used to assess the plot-to-plot and temporal variability of several key hydrometeorological variables by calculating a monthly spread. In brief, the following conclusions and broader implications can be drawn from the findings: i) The simulated total annual evapotranspiration remains relatively uniform between plots despite notable variation in its partitioning from plot to plot. ii) In the presence of a full snowpack, the albedo exhibits substantial heterogeneity at the subplot scale, linked to the canopy's LAI. iii) Local soil properties, drainage regime, and vegetation structure and type exhibit substantial influence on the plot-to-plot variability in soil water content. iv) When parameterized with localized observations and measurements, CLASS can represent and be responsive to the complex dynamics of energy and water fluxes at the plot scale within the heterogeneous surface of boreal forests.