Abstract. Glide-snow avalanches release at the ground–snow interface due to a loss in basal friction. They pose a threat to infrastructure because of the combination of unreliable mitigation measures, limited forecasting capabilities, and a lack of understanding of the release process. The aim of this study was to investigate the influence of spatial variability in basal friction and snowpack properties on the avalanche release area distribution and the release location. We developed a pseudo-3D, mechanical, threshold-based model that consists of many interacting snow columns on a uniform slope. Parameterizations in the model are based on our current understanding of glide-snow avalanche release. The model can reproduce the power law glide-snow avalanche release area distribution as observed on Dorfberg (Davos, Switzerland). A sensitivity analysis of the input parameters showed that the avalanche release area distribution was mostly influenced by the homogeneity (correlation length and variance) of the basal friction and whether the basal friction was reduced suddenly or in small increments. Larger release areas were modeled for a sudden decrease and a more homogeneous basal friction. The spatial variability of the snowpack parameters had little influence on the release area distribution. Extending the model to a real-world slope showed that the modeled location of avalanche releases qualitatively matched the observed locations. The model can help narrow down the length scales and timescales for field investigations. Simultaneously, it can grow in complexity with our increasing understanding of glide-snow avalanche release processes. Input parameters such as the basal friction or snowpack parameters could potentially all be connected to the liquid water content. This would allow for the use of meteorological measurements to drive the model. The model has the potential to help identify potentially dangerous conditions for large or numerous avalanches which would help improve glide-snow avalanche forecasting.