Probabilistic Modelling using Monte Carlo Simulation for Incorporating Uncertainty in Least Cost Path Results: a Roman Road Case Study. (DOI https://doi.org/10.31235/osf.io/mxas2)
This paper discusses the use of the Monte Carlo simulation for producing a probabilistic assessment of the accuracy of Least Cost path modelling. The author’s approach is based on generating 1000 realizations of random error fields representing the vertical error of the DEM (drawn from RMSE) and by spatially propagating this error throughout the DEM. This process is integrated into the R package leastcostpath for LCP modelling. Leastcostpath is the second designed R package of its type dealing with least Cost path modelling following Movecost (Alberti, 2019). The latter calculates accumulated slope-dependent anisotropic cost-surfaces and least-cost paths.
This paper offers an innovative approach to LCP analysis by modelling DEM uncertainty. It is well structured from a methodological and theoretical perspective. The abstract adequately reflects the content of the paper. The research question is clearly defined. The aims are stated and their significance is clearly explained. I therefore recommend that this paper be accepted. The following comments aim at providing recommendations for refining the preprint.
The Monte Carlo simulation requires substantial computational capacities for error propagation if the input DEM has a high spatial resolution (Gesch et al., 2020; Temme, Heuvelink, Schoorl, & Claessens, 2009) or when different cost functions are used comparatively for modelling the LCP. The number of modelled scenarios to be computed will moreover be increased if other factors, such as incorporating barriers with various impedance values, are included in modelling movement.
The following reference (Verhagen, Nuninger, & Groenhuijzen, 2019) could be added in the paragraph extending from lines 22-33 when discussing the role of the DEM, the slope and the key factors in LCP modelling.
Lines 60-63: It is worth highlighting that the Monte Carlo Simulation could also be applied to DEM with high accuracy and high spatial resolution. See for instance the study conducted by Gesch et al. (2020). As has been stated by Temme et al. (2009), “A high-resolution DEM may still have a greater uncertainty than a low resolution DEM if we are less certain of its attribute values”.
A study conducted by Herzog & Yépez (2015) on the impact of DEM on archaeological GIS studies could be added in lines 64-71, where the author highlights the rarity of archaeological studies discussing the vertical error and its impact on LCP modelling.
- Methodological Proposal: Temme et al. (2009) argue that simulated DEMs are not geomorphologically realistic as they have more local variation in altitude, therefore steeper slopes. In their study, the simulated DEMs were modified using sink removal algorithm and Monte Carlo analysis was performed on both unfilled and filled DEMs in order to assess the influence of sink removal on uncertainty propagation. Has the author taken into account this matter and its effect on LCP modelling?
Lines 96-97: Perhaps it would also be useful to explain in a few sentences the concept of spatial autocorrelation (see Temme et al., 2009) as it could help provide a better understanding of the usefulness of the neighborhood autocorrelation filter method.
Lines 203-205: It would be better to evaluate the RMSE of the SRTM DEM of the study area based on a number of reference points (datum points) rather than relying on consensus view of a minimum RMSE value equal to 9.73 m worldwide.
Lines 227-229: The author states that: “Effective slope, which takes into account the direction of descent, was computed in leastcostpath by calculating the difference in elevation between cells and their sixteen neighboring cells”. Is this based on the effective slope equation developed by Yu, Lee, & Munro-Stasiuk (2003)? In this case, the distance between cell centers should be also taken into account in the calculation of the slope value (See also Herzog 2014).
Lines 239-240: It is stated that the maximum distance from the known location of the High Street Roman road to the computed LCP from the south-to-north is 85% less than the maximum distance of the LCP calculated from north-to-south. In lines 292-294 the author also notes that the LCP from north-to-south is less accurate with 85% greater maximum distance than the LCP from south-to-north. Would it be possible to clarify how this percentage was calculated? Based on the numbers exposed in table 1 the maximum distance from the Least Cost Path calculated from North to South to the known location of the High Street Roman road of is 825.93 m while the maximum distance from the Least Cost Path calculated from South to North is 332.97 m. Thus, the percentage decrease from 825.93 to 332.97 is 59.68% while the percentage increase from 332.97 to 825.93 is 148.04%.
It would be preferable to add a scale and a north arrow to figures 6-8 since the author designates the computed LCP according to cardinal directions.
Alberti, G. (2019). movecost: An R package for calculating accumulated slope-dependent anisotropic cost-surfaces and least-cost paths. SoftwareX, 10, 100331. doi: https://doi.org/10.1016/j.softx.2019.100331.
Gesch, D., Palaseanu-Lovejoy, M., Danielson, J., Fletcher, C., Kottermair, M., Barbee, M., & Jalandoni, A. (2020). Inundation Exposure Assessment for Majuro Atoll, Republic of the Marshall Islands Using A High-Accuracy Digital Elevation Model. Remote Sensing, 12(1). doi: 10.3390/rs12010154.
Herzog, I., (2014). Least-cost Paths – Some Methodological Issues. IA. https://doi.org/10.11141/ia.36.5 461.
Herzog, I., & Yépez, A. (2015). The impact of the DEM on archaeological GIS studies: A case study in Ecudaor Paper presented at the Conference on Cultural Heritage and New Technologies, Vienna.
Temme, A. J. A. M., Heuvelink, G. B. M., Schoorl, J. M., & Claessens, L. (2009). Chapter 5 Geostatistical Simulation and Error Propagation in Geomorphometry. In T. Hengl & H. I. Reuter (Eds.), Developments in Soil Science (Vol. 33, pp. 121-140): Elsevier.
Verhagen, P., Nuninger, L., & Groenhuijzen, M. R. (2019). Modelling of Pathways and Movement Networks in Archaeology: An Overview of Current Approaches. In P. Verhagen, J. Joyce & M. R. Groenhuijzen (Eds.), Finding the Limits of the Limes: Modelling Demography, Economy and Transport on the Edge of the Roman Empire (pp. 217-249). Cham: Springer International Publishing.
Yu, C., Lee, J. A. Y., & Munro-Stasiuk, M. J. (2003). Research Article: Extensions to least-cost path algorithms for roadway planning. International Journal of Geographical Information Science, 17(4), 361-376. doi: 10.1080/1365881031000072645.