Snow avalanche hazard mapping has a long tradition in the European Alps. Hazard maps delineate areas of potential avalanche danger and are only available for selected areas where people and significant infrastructure are endangered. They have been created over generations, at specific sites, mainly based on avalanche activity in the past. For a large part of the area (90 % in the case of the canton of Grisons) only strongly generalized hazard indication maps are available (SilvaProtect), not showing impact information such as pressure. This is a problem when new territory with no or an incomplete historical record is to be developed. It is an even larger problem when trying to predict the effects of climate change at the state scale, where the historical record may no longer be valid. To close this gap, we develop an automated approach to generate spatially coherent hazard indication mapping based on a digital elevation model for the canton of Grisons (7105 km²) in the Swiss Alps. We calculate eight different scenarios with return periods ranging from frequent to very rare as well as with and without taking the protective effects of the forest into account, resulting in a total of approximately 2 million individual avalanche simulations. This approach combines the automated delineation of potential release areas, the calculation of release depths and the numerical simulation of the avalanche dynamics. We find that between 47 % (most frequent scenario) and 67 % (most extreme scenario) of the cantonal area can be affected by avalanches. Without forest, approximately 20 % more area would be endangered. This procedure can be applied worldwide, where high-spatial-resolution digital elevation models, detailed information on the forest and data on the snow climate are available, enabling reproducible hazard indication mapping also in regions where no avalanche hazard maps yet exist. This is invaluable for climate change studies. The simulation results are validated with official hazard maps, by assessments of avalanche experts, and by existing avalanche cadastres derived from manual mapping and mapping based on satellite datasets. The results for the canton of Grisons are now operationally applied in the daily hazard assessment work of the authorities. Based on these experiences, the proposed approach can be applied for further mountain regions.