Abstract. Snow avalanche hazard is threatening people and infrastructure in all alpine regions with seasonal or permanent snow cover around the globe. Coping with this hazard is a big challenge and during the past centuries, different strategies were 10 developed. Today, in Switzerland, experienced avalanche engineers produce hazard maps with a very high reliability based on avalanche cadastre information, terrain analysis, climatological datasets and numerical modelling of the flow dynamics for selected avalanche tracks that might affect settlements. However, for regions outside the considered settlement areas such areawide hazard maps are not available mainly because of the too high cost, in Switzerland and in most mountain regions around the world. Therefore, hazard indication maps, even though they are less reliable and less detailed, are often the only spatial 15 planning tool available. To produce meaningful and cost-effective avalanche hazard indication maps over large regions (regional to national scale), automated release area delineation has to be combined with volume estimations and state-of-theart numerical avalanche simulations.In this paper we validate existing potential release area (PRA) delineation algorithms, published in peer-reviewed journals, 20 that are based on digital terrain models and their derivatives such as slope angle, aspect, roughness and curvature. For validation, we apply avalanche cadastre data from three different ski resorts in the vicinity of Davos, Switzerland, where experienced ski-patrol staff mapped most avalanches in detail since many years. After calculating the best fit input parameters for every tested algorithm, we compare their performance based on the reference datasets. Because all tested algorithms do not provide meaningful delineation between individual potential release areas (PRA), we propose a new algorithm based on object-25 based image analysis (OBIA). In combination with an automatic procedure to estimate the average release depth (d0), defining the avalanche release volume, this algorithm enables the numerical simulation of thousands of avalanches over large regions applying the well-established avalanche dynamics model RAMMS. We demonstrate this for the region of Davos for two hazard scenarios, frequent (10 -30 years return period) and extreme (100 -300 years return period). This approach opens the door for large scale avalanche hazard indication mapping in all regions where high quality and resolution digital terrain models and 30 snow data are available.Nat. Hazards Earth Syst. Sci. Discuss., https://doi