Considering the overestimation of the maximum runout distance projected by the simulations of both case studies with the low-resolution DEMs, a calibration of the normal (EN) and tangential (ET) coefficients of restitution, as well as the rolling friction coefficient (AT) is required. In this section, the values of these coefficients are manually adjusted for each DEM scenario to obtain simulations that more closely resemble actual events. Specifically, to reduce the simulated runout distances, progressively lower restitution coefficients and higher friction coefficients are tested until the observation is matched. These adjustments compensate for different morphological roughness represented in each DEM.
4.2.1. Myloi Rockfall
After the coefficient adjustment procedure, it was determined that, for the Myloi case, only the Greek Cadastral DEM yielded reliable results, as the other DEMs (ALOS AW3D30 DEM, ASTER GDEM, SRTM30 DEM, SRTM90 DEM, and TanDEM_X) proved unsuitable for small-scale events like the Myloi rockfall, leading to unrealistic outputs. Regarding the UAV DEM, no new trials were conducted as the results with the initial coefficient values already align with the in situ observations. The forthcoming results showcase the transit frequency and maximum translational kinetic energy from the Greek Cadastral DEM, obtained with the revised coefficients shown in
Table 3.
The results from the Greek Cadastral DEM reveal that the transit frequency values were high, nearing 1000, in the rockfall source area and along the main rockfall path (
Figure 12a). In terms of the maximum translational kinetic energy, the HY-STONE calculations showed energy levels ranging from 0 to 500 kJ. The rockfall source area and the endpoint recorded minimal energies, whereas the main body exhibited significantly higher energies, up to 500 kJ (
Figure 12b). This simulation using the Greek Cadastral DEM accurately reflects the actual rockfall event.
Figure 6b and
Figure 12a from the HY-STONE simulations illustrate the effect of coefficient customization on rockfall modeling accuracy for the Myloi site by using the Greek Cadastral DEM at a 5 m resolution. In
Figure 12a, where the values of the coefficients were fine-adjusted to match the in situ observations, the depicted rockfall trajectories accurately terminate at the exact rock endpoints that were observed on-site, demonstrating a precise localization of rockfall impact. In contrast,
Figure 6b uses standard empirically derived coefficients, resulting in trajectories that extend beyond the actual observed endpoints, suggesting an overestimation of the impact area. This comparison highlights the importance of adjusting the simulation parameters to specific site conditions to ensure the accuracy and reliability of rockfall risk assessments. Furthermore,
Figure 7b and
Figure 12b illustrate the HY-STONE results for the maximum translation kinetic energy, also at the Myloi site, using the Greek Cadastral DEM at a 5 m resolution. As can be seen in
Figure 12b, utilizing the adjusted coefficients for a closer match to the in situ observations shows a precise distribution of kinetic energy with high values concentrated near the rockfall source, thus closely reflecting the actual rockfall paths. In contrast,
Figure 7b employs empirically derived coefficients, resulting in a more dispersed and extended kinetic energy distribution across the area.
4.2.2. Platiana Rockfall
For the Platiana rockfall, the coefficients EN, ET, and AT have been adjusted for each DEM scenario, except for the Greek Cadastral DEM, which was already calibrated in the previous section against the in situ observations. Presented next are the outcomes for the transit frequency and maximum translational kinetic energy from the simulations using the ALOS AW3D30 DEM, ASTER GDEM, SRTM30 DEM, SRTM90 DEM, and TanDEM_X using the optimal restitution and rolling friction coefficients for each scenario.
Table 4 and
Figure 13 showcase the optimal coefficients and simulation results in terms of the transit frequency and kinetic energies for the ALOS AW3D30 DEM.
According to
Figure 13a, the ALOS AW3D30 DEM exhibited a single corridor of trajectories that pass on the left side of the real rock endpoint and show only low transit frequencies close to the endpoint. The computation of maximum translational kinetic energy displayed a range of values, beginning with low energies (50–100 kJ), escalating to very high energies (500–1000 kJ), and finally ending again with low energies (50–100 kJ) before reaching the endpoint (
Figure 13b).
The comparison between
Figure 13a and
Figure 9b highlights the differences in the transit frequency before and after the parameter adjustment. While
Figure 13a depicts a very concentrated narrow path, suggestive of a well-defined rockfall trajectory,
Figure 9b shows a more dispersed and broad frequency distribution that extends beyond the actual observed endpoints, implying either different simulation settings or more complex terrain interactions. Similarly,
Figure 10b and
Figure 13b contrast in their portrayal of kinetic energy;
Figure 13b presents a linear decrease in energy from the source along a controlled path, whereas
Figure 10b depicts a diffuse energy pattern with values further from the rock endpoint.
Table 5 and
Figure 14 showcase the optimal coefficients and simulation results in terms of the transit frequency and kinetic energies for the ALOS AW3D30 DEM.
Regarding the transit frequency maps, the ASTER GDEM documented (
Figure 14a) high to low transit frequencies during the rock path crossing the rock stop point.
Concerning the kinetic energy calculations, the ASTER GDEM develops similarly to the ALOS AW3D30 DEM simulation, starting with low values, followed by an increase until the final part of the trajectories where the blocks stopped.
Figure 14a presents a single concentrated transit frequency along a narrow path, reflecting a precise and realistic depiction of rockfall trajectory. In contrast,
Figure 9c results in a broader and less defined distribution of transit frequencies, indicating multiple potential trajectories and a potential overestimation of the impact area. According to the kinetic energies,
Figure 10c shows a detailed gradient of kinetic energy with multiple zones ranging from low- to high-energy values, reflecting a nuanced understanding of energy dissipation along the trajectory. This gradient indicates a more refined simulation where kinetic energy decreases progressively from the source, offering a detailed perspective on energy distribution during rockfall. On the other hand,
Figure 14b exhibits a more simplified energy distribution, primarily concentrating high kinetic energy near the rockfall source with a sharp drop to lower values, suggesting a less complex interaction with the terrain.
Table 6 and
Figure 15 showcase the optimal coefficients and simulation results in terms of the transit frequency and kinetic energies for the SRTM30 DEM.
The SRTM30 DEM displayed medium transit frequencies over an extended section of the rock path, high frequencies in the source area, and no recorded values at the rock endpoint (
Figure 15c). The kinetic energy results delineated two rock paths: the left, with energy levels up to 500–1000 kJ, and the right, with lower values around 200 kJ (
Figure 15d).
Figure 9d,
Figure 10d, and
Figure 15a,b present the rockfall simulation results for the transit frequency and kinetic energies, correspondingly, using the SRTM30 DEM.
Figure 15a,b uses the adjusted coefficients, while
Figure 9d and
Figure 10d use coefficients on the highest resolution DEM.
Figure 15a shows a more diverse color spectrum, indicating a dynamic interaction with a varied terrain that ends at the left side of the rock endpoint, while
Figure 9d presents a more uniform and concentrated frequency, suggesting simpler terrain interactions or model settings that end after the rock endpoint.
Figure 10d displays a linear distribution of kinetic energy with clear high-energy zones concentrated near the rockfall source, suggesting a direct rockfall trajectory. In contrast,
Figure 15b shows a more uniform distribution, with fewer fluctuations in energy levels along the trajectory, indicating a potential oversimplification of the terrain’s impact on energy dissipation.
Table 7 and
Figure 16 showcase the optimal coefficients and simulation results in terms of the transit frequency and kinetic energies for the SRTM90 DEM.
The SRTM90 DEM shows a constant decrease in the transit frequencies along the slope (
Figure 16a), resembling the results of the ASTER GDEM and ALOS AW3D30 DEM. Regarding the kinetic energy (
Figure 16b), the SRTM90 DEM exhibited initial lower energies of 200–500 kJ at the source area, increasing to 500–1000 kJ along the slope, without any decrease in energy before the rock stop point (
Figure 16b).
The transit frequency results shown in
Figure 9e and
Figure 16a demonstrate a vertical layout with varying degrees of distribution;
Figure 9e shows a poor gradation of the frequency zones along the slope continuing after the rock endpoint, whereas
Figure 16a presents a more streamlined frequency layout ending at the rock endpoint. Similarly, the kinetic energy results in
Figure 10e and
Figure 16b display linear energy distributions along the rockfall path.
Figure 10e shows a complex pattern of low- and high-energy zones near the rockfall source and at the endpoint. On the other hand,
Figure 16b portrays a simpler and more uniform distribution of kinetic energy, highlighting the role of the model’s parameters in shaping the energy dissipation outcomes.
Table 8 and
Figure 17 showcase the optimal coefficients and simulation results in terms of the transit frequency and kinetic energies for the TanDEM_X.
According to the TanDEM_X rockfall simulation map, the transit frequency for the Platiana rockfall is high in the source area and low at the endpoint of the rockfall (
Figure 17a). The maximum translational kinetic energies show little variation along the slope, ranging in the 200–500 kJ class (
Figure 17b).
The comparison between
Figure 17a and
Figure 9f concluded that the simulation with the parameters optimized on the TanDEM_X presents a simplified and very linear representation of the transit frequency ending at the rock endpoint, which strictly follows a vertical pattern, indicating a highly targeted and specific rockfall path. Meanwhile,
Figure 9f displays a more complex distribution with a smaller spectrum of transit frequencies, suggesting more variability in rockfall behavior. Last but not least, the compared
Figure 17b and
Figure 10f sum up that both figures show a constant vertical distribution of kinetic energy values, but
Figure 17b depicts a less spread-out energy distribution with a more realistic ending point.
To explore how the spatial resolution of the DEMs influences the normal (E
N) and tangential (E
T) restitution coefficients and the rolling friction coefficient (A
T), a plot reporting the differences in the coefficients for the most representative land cover class (debris without vegetation) is shown for the Platiana case (
Figure 18). This plot clearly shows that the optimal normal (E
N) and tangential (E
T) coefficients decrease with increasing spatial resolution, whereas the rolling friction coefficient (A
T) increases with a higher spatial resolution.
The graph shows a roughly logarithmic relationship between the spatial resolution of DEMs and the optimal rockfall simulation coefficients. As the pixel size increases, indicative of coarser spatial resolution, the rolling friction coefficient (AT) exhibits an increasing trend. This suggests that larger pixel sizes, which represent less detailed terrain data, demand high values of rolling friction coefficient to halt rocks more rapidly on the slope.
In contrast, the restitution coefficients, EN and ET, show a decreasing trend with increasing pixel size. This pattern indicates that, in simulations, rocks lose more energy upon impact with coarser-resolution DEMs, resulting in shorter bounces and less travel distance, which is characteristic of more inelastic collisions.
The correlation depicted here highlights the significance of selecting DEMs with appropriate pixel sizes for rockfall modeling to ensure that the physical behavior of rockfalls is realistically replicated in the simulation outcomes. The variability of the coefficients for a certain resolution is due to slight differences in the morphological description of the different DEMs. Future research will expand the dataset to validate these findings further. However, the observed trend provides valuable insights into the behavior of each coefficient in relation to each spatial resolution.