Amphibian species specific landscape resistance maps.

For each amphibian species we used the weighted average Ensemble model suitability map and converted it to a resistance map using an exponential decay function following (Duflot et al., 2018) using the species-specific mpa threshold previously computed (Preau et al., 2020). The function assigns a resistance value of 1,000 when habitat suitability equals zero, and a value of one when habitat suitability is greater than or equal to the binarization threshold. The exponential shape strengthens the barrier effect of less favorable areas.