“…A variety of methodological approaches have been developed over the last decades to generate SDMs, such as Neural Networks (Özesmi and Özesmi, 1999), Boosted Regression Trees (Elith et al, 2008;Wege et al, 2021), Maximum Entropy (Phillips et al, 2006), Generalized Linear Regression Model (Guisan et al, 2002), and Additive Regression Models (Swartzman et al, 1992;Austin, 2007). However, the statistical challenges using SDMs have increased as datasets have become more complex over time (Orúe et al, 2020;Lloret-Lloret et al, 2021). Recently, Bayesian hierarchical spatiotemporal models with the Integrated Nested Laplace Approximation (INLA) methodology has acquired an important role, since it is ideally suited for fitting complex spatiotemporal covariance structures (Rue et al, 2009;Silva et al, 2017).…”