A Monte Carlo method to account for sampling error in multi-species indicators

Abstract The usefulness of biodiversity indicators strongly increases if accompanied by measures of uncertainty. In the case of indicators that combine population indices of species, however, the inclusion of the uncertainty of the species indices has shown to be hard to realize, usually due to imperfections in monitoring programmes. Missing values and time series of different lengths preclude the use of analytical approaches, whereas bootstrapping across sites requires the raw abundance data on the site level, which may not always be available. Sometimes bootstrapping across species rather than sites is opted for, but this approach ignores the uncertainty attached to species indices. We developed a method to account for sampling error of species indices in the calculation of multi-species indicators based on Monte Carlo simulation of annual species indices. The construction of confidence intervals enables various trend assessments, like testing for linear or smooth trends, testing for changes between two time points, testing the significance of a suspected change-point and testing for differences between two multi-species indicators. Here, we compare our method with conventional methods and illustrate the benefits of our approach using Dutch breeding bird indicators.