Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species’ current and future distributions

Bioclimatic models are widely used tools for assessing potential responses of species to climate change. One commonly used model is BIOCLIM, which summarises up to 35 climatic parameters throughout a species’ known range, and assesses the climatic suitability of habitat under current and future climate scenarios. A criticism of BIOCLIM is that the use of all 35 parameters may lead to over-fitting of the model, which in turn may result in misrepresentations of species’ potential ranges and to the loss of biological reality. In this study, we investigated how different methods of combining climatic parameters in BIOCLIM influenced predictions of the current distributions of 25 Australian butterflies species. Distributions were modeled using three previously used methods of selecting climatic parameters: (i) the full set of 35 parameters, (ii) a customised selection of the most relevant parameters for individual species based on analysing histograms produced by BIOCLIM, which show the values for each parameter at all of the focal species known locations, and (iii) a subset of 8 parameters that may generally influence the distributions of butterflies. We also modeled distributions based on random selections of parameters. Further, we assessed the extent to which parameter choice influenced predictions of the magnitude and direction of range changes under two climate change scenarios for 2020. We found that the size of predicted distributions was negatively correlated with the number of parameters incorporated in the model, with progressive addition of parameters resulting in progressively narrower potential distributions. There was also redundancy amongst some parameters; distributions produced using all 35 parameters were on average half the size of distributions produced using only 6 parameters. The selection of parameters via histogram analysis was influenced, to an extent, by the number of location records for the focal species. Further, species inhabiting different biogeographical zones may have different sets of climatic parameters limiting their distributions; hence, the appropriateness of applying the same subset of parameters to all species may be reduced under these situations. Under future climates, most species were predicted to suffer range reductions regardless of the scenario used and the method of parameter selection. Although the size of predicted distributions varied considerably depending on the method of selecting parameters, there were no significant differences in the proportional change in range size between the three methods: under the worst-case scenario, species’ distributions decrease by an average of 12.6, 11.4, and 15.7%, using all parameters, the ‘customised set’, and the ‘general set’ of parameters, respectively.