Skip to main content

Model uncertainties in multispecies ecological models

We live in an increasingly uncertain world.  Therefore, when we model environmental processes of interest, it is vital to account for the inherent uncertainties in our analyses and ensure that this information is communicated to relevant parties.  Whilst the use of complex statistical models to estimate quantities of interest is becoming increasingly common in environmental sciences, one aspect of uncertainty that is frequently overlooked is that of model uncertainty.  Much of the research I conduct considers this additional aspect of uncertainty quantification; that is not just uncertainty in the quantities of interest, but also in the models that we use to estimate them.

An example of this is in a paper recently published in Ecology and Evolution (Swallow et al., 2016), which looks at how different species of birds that we commonly see in our gardens respond to the same environmental factors (or covariates).  Some of the species have declined rapidly over the past 40 years, whilst others have remained stable or even increased in number.  Possible drivers of these changes that have been suggested include increases in predators, changes in climate and availability of natural food sources.  Statistically speaking, we try to understand and quantify changes in observed numbers of birds by relating them to changes in measured environmental quantities that the birds will be subjected to, such as numbers of predators, weather variables, habitat quality etc.  Most previous analyses have modelled each of the species observed at many different geographical locations (or monitoring sites) independently of each other, and estimated the quantities of interest completely separately, despite the fact that all these species share the same environment and are subject to the same external influences.  So how do we go about accounting for the fact that similar species may share similar population drivers?

This essentially constitutes a model uncertainty problem – that is, which parameters should be shared across which species in our statistical model and which parameters should be distinct?

If we were to consider two different species and use two different environmental factors to explain changes in those species, say habitat type and average monthly temperature, there are four possible models to consider.  That is,

Model
Habitat type
Temperature
No parameters
1
Shared
Shared
2
2
Distinct
Shared
3
3
Shared
Distinct
3
4
Distinct
Distinct
4

This can easily be extended to a higher number of species and covariates.

There is also inevitably going to be some aspects of variability shown by some of the species that we cannot account for through the quantities we have measured.  We account for this using site-specific random effects, which explain variability that is linked to a specific monitoring site, but which is not accounted for by the environmental covariates in the model.  Again, we would usually assume this is a single quantity representing the discrepancy between what we have accounted for using our measured covariates and what is ‘left over’.  Following on from work of previous authors (Lahoz-Monfort et al., 2011), we again split this unexplained variation into two – unexplained variation that is common to all species and unexplained variation that is specific to a single species.  The ratio of these two quantities can give us a good idea of what measurements we may be missing.  Is it additional environmental factors that are wide-ranging in their effects or is it something relating to the specific ecology of an individual species?

In the paper, we apply our method to a large dataset spanning nearly 40 years, collected as part of the British Trust for Ornithology’s Garden Bird Feeding Survey.  We selected two groups of similar species commonly found in UK gardens during the winter.  For ecological reasons, we would expect the species within the two groups to show similar traits, so they act as ideal study species for detecting synchrony in responses to environmental factors.  Whilst most the results were consistent with those from single-species models (e.g. Swallow et al., 2015), studying the species at an ecosystem level also highlighted some additional relationships that it would be impossible to study under more simplistic models.  The results highlight that there is unsurprisingly a large degree of synchrony across many of these species, and that they share many of the traits and drivers of population change.  The synchronies observed in the results corresponded to both significant positive or negative relationships with covariates, as well as those species that collectively show no strong relationship with a given environmental factor.  There is, however, more to the story and some of the species showed strong differences in how they respond to external factors.  Highlighting these differences may offer important information on how best to halt or reverse population declines.

The results from our analyses showed the importance of considering model uncertainty in statistical analyses of this type, and that by incorporating relevant uncertainties, we can improve our understanding of the environmental processes of interest.  Incorporating more data into the analysis will help in further constraining common or shared parameters and reduce uncertainties in them.  It also allows us to guide and improve future data collection procedures if we can gain a better understanding of what is currently missing from our model.

Blog written by Dr Ben Swallow, a Postdoctoral Research Associate, studying Ecological and environmental statistics in the School of Chemistry.

References

Lahoz-Monfort, J. J., Morgan, B. J. T., Harris, M. P., Wanless, S., & Freeman, S. N. (2011). A capture-recapture model for exploring multi-species synchrony in survival. Methods in Ecology and Evolution, 2(1), 116–124.

Swallow, B., Buckland, S. T., King, R. and Toms, M. P. (2015). Bayesian hierarchical modelling of continuous non-negative longitudinal data with a spike at zero: An application to a study of birds visiting gardens in winter. Biometrical Journal, 58(2), 357–371

Swallow, B., King, R., Buckland, S. T. and Toms, M. P. (2016). Identifying multispecies synchrony in response to environmental covariates. Ecology and Evolution, 6(23), 8515–8525

Figure 1. Blue tits show a highly synchronous response with great tits, and to a lesser degree coal tits, to their surrounding environment.



Figure 2. Male house sparrow feeding on fat balls.  Whilst they show some synchrony in their response to environmental factors, they appear to be subject to a differing ecology to the other two species they were compared with.




Popular posts from this blog

Converting probabilities between time-intervals

This is the first in an irregular sequence of snippets about some of the slightly more technical aspects of uncertainty and risk assessment.  If you have a slightly more technical question, then please email me and I will try to answer it with a snippet. Suppose that an event has a probability of 0.015 (or 1.5%) of happening at least once in the next five years. Then the probability of the event happening at least once in the next year is 0.015 / 5 = 0.003 (or 0.3%), and the probability of it happening at least once in the next 20 years is 0.015 * 4 = 0.06 (or 6%). Here is the rule for scaling probabilities to different time intervals: if both probabilities (the original one and the new one) are no larger than 0.1 (or 10%), then simply multiply the original probability by the ratio of the new time-interval to the original time-interval, to find the new probability. This rule is an approximation which breaks down if either of the probabilities is greater than 0.1. For example

1-in-200 year events

You often read or hear references to the ‘1-in-200 year event’, or ‘200-year event’, or ‘event with a return period of 200 years’. Other popular horizons are 1-in-30 years and 1-in-10,000 years. This term applies to hazards which can occur over a range of magnitudes, like volcanic eruptions, earthquakes, tsunamis, space weather, and various hydro-meteorological hazards like floods, storms, hot or cold spells, and droughts. ‘1-in-200 years’ refers to a particular magnitude. In floods this might be represented as a contour on a map, showing an area that is inundated. If this contour is labelled as ‘1-in-200 years’ this means that the current rate of floods at least as large as this is 1/200 /yr, or 0.005 /yr. So if your house is inside the contour, there is currently a 0.005 (0.5%) chance of being flooded in the next year, and a 0.025 (2.5%) chance of being flooded in the next five years. The general definition is this: ‘1-in-200 year magnitude is x’ = ‘the current rate for eve

Coconuts and climate change

Before pursuing an MSc in Climate Change Science and Policy at the University of Bristol, I completed my undergraduate studies in Environmental Science at the University of Colombo, Sri Lanka. During my final year I carried out a research project that explored the impact of extreme weather events on coconut productivity across the three climatic zones of Sri Lanka. A few months ago, I managed to get a paper published and I thought it would be a good idea to share my findings on this platform. Climate change and crop productivity  There has been a growing concern about the impact of extreme weather events on crop production across the globe, Sri Lanka being no exception. Coconut is becoming a rare commodity in the country, due to several reasons including the changing climate. The price hike in coconuts over the last few years is a good indication of how climate change is affecting coconut productivity across the country. Most coconut trees are no longer bearing fruits and thos