From a computational-archaeology perspective, this paper is especially valuable because it makes a large-scale archaeological argument through a chain of spatial models.
The workflow combines Neanderthal occurrence data with six environmental predictors, five modelling algorithms (GLM, GAM, random forest, GBM, and MaxNet) and an ensemble species-distribution model. The resulting suitability surfaces are then used in circuit-theory connectivity analysis to estimate potential movement pathways. This is a strong framework for testing whether climatic suitability became more fragmented through time.
However, the supplementary code also shows why model outputs must remain open to scrutiny. As supplied, the script creates a five-fold cross validation object with five repetitions, but the modelling call appears to use two repetitions and a 70/30 data split. The supplied loops also execute the LIG stage unless adapted. Suitable area is calculated by summing binary raster cells, and the connectivity model uses specific choices :: four-neighbour movement, raster aggregation, a fixed parameter, and an implicit suitability-to-conductance relationship.
These are not minor technical details. They can influence the resulting corridors and landscape metrics. The lesson is not that the model is invalid, it is that a corridor is a model-dependent hypothesis, not an observed Neanderthal route. Future work can add spatial-block validation, resistance-surface sensitivity tests, uncertainty propagation, independent archaeological validation, and explicit modelling of prey, water, seasonality, and demography.
Supplementary Data:
https://media.springernature.com/original/springer-static/esm/art%3A10.1038%2Fs43247-026-03960-8/MediaObjects/43247_2026_3960_MOESM2_ESM.pdf
#ComputationalArchaeology #GIS #SpeciesDistributionModels #ReproducibleResearch