P. 285-294The estimated potential of landscape metrics as a surrogate for biodiversity is strongly dependent on the spatial analytical unit used for evaluation. We assessed the relationship between terrestrial vertebrate species richness (total and taxonomic) and structural landscape heterogeneity, testing the impact of using different spatial analytical units in three mountain systems in Spain. Landscape heterogeneity was quantified through an additive partitioning of the Shannon diversity index of landscape classes. Both landscape heterogeneity and species richness were calculated using two spatial analytical unit approaches: eco-geographic vs. arbitrary (i.e., watersheds vs. square windows of different sizes 20 × 20 km, 50 × 50 km, 100 × 100 km). We predicted species richness on the basis of landscape heterogeneity by fitting separate linear models for each spatial analytical unit approach. The main results obtained showed that landscape heterogeneity influenced terrestrial vertebrate species richness. However, the emerging relationships were dependent on the spatial analytical unit approach. The eco-geographic approach showed significant relationships between landscape heterogeneity and total and taxonomic species richness in almost all cases (except mammals). Considering the arbitrary approach, landscape heterogeneity appeared as a predictor of species richness only for mammals and breeding birds and at the coarsest spatial scales. Our results claim for further consideration of eco-geographical spatial analytical unit approaches in biodiversity studies and show that the methods of this study offer a valuable cost-effective framework for biodiversity management and spatial modeling, with potential to be adapted to national and global applications.S