Data error propagation in stacked bioclimatic envelope models

Stacking is the process of overlaying inferred species potential distributions for multiple species based on outputs of bioclimatic envelope models (BEMs). The approach can be used to investigate patterns and processes of species richness. If data limitations on individual species distributions are inevitable, but how do they affect inferences of patterns and processes of species richness? We investigate the influence of different data sources on estimated species richness gradients in China. We fitted BEMs using species distributions data for 334 bird species obtained from (1) global range maps, (2) regional checklists, (3) museum records and surveys, and (4) citizen science data using presence-only (Mahalanobis distance), presence-background (MAXENT), and presence-absence (GAM and BRT) BEMs. Individual species predictions were stacked to generate species richness gradients. Here, we show that different data sources and BEMs can generate spatially varying gradients of species richness. The environmental predictors that best explained species distributions also differed between data sources. Models using citizen-based data had the highest accuracy, whereas those using range data had the lowest accuracy. Potential richness patterns estimated by GAM and BRT models were robust to data uncertainty. When multiple data sets exist for the same region and taxa, we advise that explicit treatments of uncertainty, such as sensitivity analyses of the input data, should be conducted during the process of modeling.