Mapping small mammal optimal habitats using satellite-derived proxy variables and species distribution models

Small mammal species play an important role influencing vegetation primary productivity and plant species composition, seed dispersal, soil structure, and as predator and/or prey species. Species which experience population dynamics cycles can, at high population phases, heavily impact agricultural sectors and promote rodent-borne disease transmission. To better understand the drivers behind small mammal distributions and abundances, and how these differ for individual species, it is necessary to characterise landscape variables important for the life cycles of the species in question. In this study, a suite of Earth observation derived metrics quantifying landscape characteristics and dynamics, and in-situ small mammal trapline and transect survey data, are used to generate random forest species distribution models for nine small mammal species for study sites in Narati, China and Sary Mogul, Kyrgyzstan. These species distribution models identify the important landscape proxy variables driving species abundance and distributions, in turn identifying the optimal conditions for each species. The observed relationships differed between species, with the number of landscape proxy variables identified as important for each species ranging from 3 for Microtus gregalis at Sary Mogul, to 26 for Ellobius tancrei at Narati. Results indicate that grasslands were predicted to hold higher abundances of Microtus obscurus, E. tancrei and Marmota baibacina, forest areas hold higher abundances of Myodes centralis and Sorex asper, with mixed forest-grassland boundary areas and areas close to watercourses predicted to hold higher abundances of Apodemus uralensis and Sicista tianshanica. Localised variability in vegetation and wetness conditions, as well as presence of certain habitat types, are also shown to influence these small mammal species abundances. Predictive application of the Random Forest (RF) models identified spatial hot-spots of high abundance, with model validation producing R2 values between 0.670 for M. gregalis transect data at Sary Mogul to 0.939 for E. tancrei transect data at Narati. This enhances previous work whereby optimal habitat was defined simply as presence of a given land cover type, and instead defines optimal habitat via a combination of important landscape dynamic variables, moving from a human-defined to species-defined perspective of optimal habitat. The species distribution models demonstrate differing distributions and abundances of host species across the study areas, utilising the strengths of Earth observation data to improve our understanding of landscape and ecological linkages to small mammal distributions and abundances.