A global set of Fourier-transformed remotely sensed covariates for the description of abiotic niche in epidemiological studies of tick vector species

Background: Correlative modelling combines observations of species occurrence with environmental variables to capture the niche of organisms. It has been argued for the use of predictors that are ecologically relevant to the target species, instead of the automatic selection of variables. Without such biological background, the forced inclusion of numerous variables can produce models that are highly inflated and biologically irrelevant. The tendency in correlative modelling is to use environmental variables that are interpolated from climate stations, or monthly estimates of remotely sensed features. Methods: We produced a global dataset of abiotic variables based on the transformation by harmonic regression (time series Fourier transform) of monthly data derived from the MODIS series of satellites at a nominal resolution of 0.1 degrees. The dataset includes variables, such as day and night temperature or vegetation and water availability, which potentially could affect physiological processes and therefore are surrogates in tracking the abiotic niche. We tested the capacities of the dataset to describe the abiotic niche of parasitic organisms, applying it to discriminate five species of the globally distributed tick subgenus Boophilus and using more than 9,500 published records. Results: With an average reliability of 82%, the Fourier-transformed dataset outperformed the raw MODIS-derived monthly data for temperature and vegetation stress (62% of reliability) and other popular interpolated climate datasets, which had variable reliability (56%-65%). The transformed abiotic variables always had a collinearity of less than 3 (as measured by the variance inflation factor), in contrast with interpolated datasets, which had values as high as 300. Conclusions: The new dataset of transformed covariates could address the tracking of abiotic niches without inflation of the models arising from internal issues with the descriptive variables, which appear when variance inflation is higher than 10. The coefficients of the harmonic regressions can also be used to reconstruct the complete original time series, being an adequate complement for ecological, epidemiological, or phylogenetic studies. We provide the dataset as a free download under the GNU general public license as well as the scripts necessary to integrate other time series of data into the calculations of the harmonic coefficients.