The aim of this study was to delimit potential rockfall propagation zones based on simulated 2 m resolution rockfall trajectories using Rockyfor3D for block volume scenarios ranging from 0.05–30 m3, with explicit inclusion of the barrier effect of standing trees, for an area of approx. 7200 km2 in Switzerland and Liechtenstein. For the determination of the start cells, as well as the slope surface characteristics, we used the terrain morphology derived from a 1 m resolution digital terrain model, as well as the topographic landscape model geodataset of swisstopo and information from geological maps. The forest structure was defined by individual trees with their coordinates, diameters, and tree type (coniferous or broadleaved). These were generated from detected individual trees combined with generated trees on the basis of statistical relationships between the detected trees, remote sensing-based forest structure type definitions, and stem numbers from field inventory data. From the simulated rockfall propagation zones we delimited rockfall hazard indication zones (HIZ), as called by the practitioners (because they serve as a basis for the Swiss hazard index map), on the basis of the simulated reach probability rasters. As validation, 1554 mapped past rockfall events were used. The results of the more than 89 billion simulated trajectories showed that 94% of the mapped silent witnesses could be reproduced by the simulations and at least 82% are included in the delimited HIZ.