Research by Eleni and colleagues mentioned by the Independent

As part of a series of articles on International Bee Day 2021, The Independent mentioned work by Eleni on modelling bumblebee citizen science data. In particular, the article referred to results published in 2018, which suggested that even some common species of bumblebee are now in decline.

Recent analysis, using the new and freely available RShiny app, developed by research student Fabian Ketwaroo, suggest that the number of queens emerging from hibernation each year for species  Bombus pascuorum and Bombus hortorum continue to be stable in recent years compared to 2011, when BeeWalk data were first collected.

At the same time there is some indication of positive news for species Bombus lapidarius, also known as the red tailed bumblebee, as the number of queens emerging from hibernation is on the rise in recent years, with 2019 seeing a similar number of queens to 2011. (see plot above, showing the number of queens emerging each year in comparison to 2011).


Eleni gave invited talk during an BI-IBS event

Eleni was invited to present some of her work on modelling bumblebee citizen science data during a meeting organised by the British and Irish Region of the International Biometric Society on Advanced Topics in Spatial Sampling.

Eleni’s talk showcased a new Bayesian model, developed and implemented into a freely-available RShiny app  by research Fabian Ketwaroo and demonstrated it when modelling two species of UK bumblebees.


New paper: Removal modelling in ecology – A systematic review

SE@K members Oscar and Rachel have published a paper reviewing the current state-of-the-art in modelling removal data. This work is part of the the EPSRC project EP/S020470/1 “Modelling removal and re-introduction data for improved conservation”.

The paper is available in full here:

Removal modelling in ecology: A systematic review

Oscar Rodriguez de Rivera and Rachel McCrea

Removal models were proposed over 80 years ago as a tool to estimate unknown population size. More recently, they are used as an effective tool for management actions for the control of non desirable species, or for the evaluation of translocation management actions. Although the models have evolved over time, in essence, the protocol for data collection has remained similar: at each sampling occasion attempts are made to capture and remove individuals from the study area. Within this paper we review the literature of removal modelling and highlight the methodological developments for the analysis of removal data, in order to provide a unified resource for ecologists wishing to implement these approaches. Models for removal data have developed to better accommodate important features of the data and we discuss the shift in the required assumptions for the implementation of the models. The relative simplicity of this type of data and associated models mean that the method remains attractive and we discuss the potential future role of this technique.


New paper: A Test for the Underlying State-Structure of Hideen Markov Models – Partially Observed Capture-Recapture Data

Former SE@K PhD student Anita Jeyam, Rachel and Roger Pradel (Montpellier) have had a paper published presenting a new test for determining the underlying state structure of Hidden Markov models.  This is an exciting piece of work which provides the foundation for addressing a complex topic for general HMMs.

The full paper is available open access here:

A Test for the Underlying State-Structure of Hideen Markov Models – Partially Observed Capture-Recapture Data

Hidden Markov models (HMMs) are being widely used in the field of ecological modeling,
however determining the number of underlying states in an HMM remains a challenge.
Here we examine a special case of capture-recapture models for open populations,
where some animals are observed but it is not possible to ascertain their state (partial
observations), whilst the other animals’ states are assigned without error (complete
observations). We propose a mixture test of the underlying state structure generating
the partial observations, which assesses whether they are compatible with the set of
states observed in the complete observations. We demonstrate the good performance
of the test using simulation and through application to a data set of Canada Geese.


New paper: Modeling Recruitment of Birth Cohorts to the Breeding Population: A Hidden Markov Approach

Rachel McCrea and co-authors from St Andrews and Edinburgh have published a paper on modelling the recruitment process of gray seals.

The full paper is available open access here:

Modeling Recruitment of Birth Cohorts to the Breeding Population: A Hidden Markov Approach

Worthington, King, McCrea, Smout and Pomeroy

Long-term capture-recapture studies provide an opportunity to investigate the population dynamics of long-lived species through individual maturation and adulthood and/or
time. We consider capture-recapture data collected on cohorts of female gray seals
(Halichoerus grypus) born during the 1990s and later observed breeding on the Isle
of May, Firth of Forth, Scotland. Female gray seals can live for 30+ years but display
individual variability in their maturation rates and so recruit into the breeding population
across a range of ages. Understanding the partially hidden process by which individuals
transition from immature to breeding members, and in particular the identification of
any changes to this process through time, are important for understanding the factors
affecting the population dynamics of this species. Age-structured capture-recapture
models can explicitly relate recruitment, and other demographic parameters of interest,
to the age of individuals and/or time. To account for the monitoring of the seals from
several birth cohorts we consider an age-structured model that incorporates a specific
cohort-structure. Within this model we focus on the estimation of the distribution of
the age of recruitment to the breeding population at this colony. Understanding this
recruitment process, and identifying any changes or trends in this process, will offer
insight into individual year effects and give a more realistic recruitment profile for the
current UK gray seal population model. The use of the hidden Markov model provides an
intuitive framework following the evolution of the true underlying states of the individuals.
The model breaks down the different processes of the system: recruitment into the
breeding population; survival; and the associated observation process. This model
specification results in an explicit and compact expression for the model with associated
efficiency in model fitting. Further, this framework naturally leads to extensions to more
complex models, for example the separation of first-time from return breeders, through
relatively simple changes to the mathematical structure of the model.



New Paper: Parameter redundancy in Jolly‐Seber tag loss models

Diana along with Wei Cai, Stephanie Yurchak and Laura Cowen have published the paper:

Parameter redundancy in Jolly‐Seber tag loss models

in Ecology and Evolution


1. Capture–recapture experiments are conducted to estimate population parameters such as population size, survival rates, and capture rates. Typically, individuals are captured and given unique tags, then recaptured over several time periods with the assumption that these tags are not lost. However, for some populations, tag loss cannot be assumed negligible. The Jolly‐Seber tag loss model is used when the no‐tag‐loss assumption is invalid. Further, the model has been extended to incorporate group heterogeneity, which allows parameters to vary by group membership. Many mark–recapture models become overparameterized resulting in the inability to independently estimate parameters. This is known as parameter redundancy.
2. We investigate parameter redundancy using symbolic methods. Because of the complex structure of some tag loss models, the methods cannot always be applied directly. Instead, we develop a simple combination of parameters that can be used to investigate parameter redundancy in tag loss models.
3. The incorporation of tag loss and group heterogeneity into Jolly‐Seber models does not result in further parameter redundancies. Furthermore, using hybrid methods we studied the parameter redundancy caused by data through case studies and generated tag histories with different parameter values.
4. Smaller capture and survival rates are found to cause parameter redundancy in these models. These problems resolve when applied to large populations.


New Paper by Emily and Byron: Integrated modelling of insect population dynamics at two temporal scales

Emily and Byron, along with Marc Kery, Armin Coray, Michael Schaub and Bruno Baur have published the paper:

Integrated modelling of insect population dynamics at two temporal scales

in Ecological Modelling.


Population size of species with birth-pulse life-cycles varies both within and between seasons, but most population dynamics models assume that a population can be characterised adequately by a single number within a season. However, within-season dynamics can sometimes be too substantial to be ignored when modelling dynamics between seasons. Typical examples are insect populations or migratory animals. Numerous models for only between-season dynamics exist, but very few have combined dynamics at both temporal scales.

In a new approach, we extend appreciably the models of Dennis et al. (2016b): we show how to adapt them for a generation time  year and fit an integrated population model for multiple data types, by maximising a joint likelihood for population counts of unmarked individuals and capture–recapture data from a study with marked individuals. We illustrate the approach using annual monitoring data for the endangered flightless beetle Iberodorcadion fuliginator from 18 populations in the Upper Rhine Valley for 1998–2016, with a 2-year life cycle. Standard likelihood methods are used for model fitting and comparison, and a concentrated (profile) likelihood approach provides computational efficiency.

Additional information from the capture–recapture data makes the population model more robust and, importantly, enables true, rather than relative, abundance to be estimated. A dynamic stopover model provides estimates of both survival and phenology parameters within a season, and also of productivity between seasons. For  I. fuliginator, we demonstrate a population decline since 1998 and how this links with productivity, which is affected by temperature. A delayed mean emergence date in recent years is also shown.

A main point of interest is the focus on the two temporal scales at which perhaps most animal populations vary: in the short-term, a population is seldom truly closed within a single season, and in the long-term (between seasons) it never is. Hence our models may serve as a template for a general description of population dynamics in many species. This includes rare species with limited data sets, for which there is a general lack of population dynamic models, yet conservation actions may greatly benefit from this kind of models.


PhD Project: Modelling Citizen-science Data for Butterflies and Moths: Where and When are they Flying?

Eleni Matechou, Emily Dennis, Richard Fox, Byron Morgan and Diana Cole are offering the statistics PhD project: Modelling Citizen-science Data for Butterflies and Moths: Where and When are they Flying?
At a time of biodiversity loss, including widely reported insect declines, citizen science data play a vital role in measuring changes in species’ populations and distributions and in seeking to understand the pressures influencing such changes. Butterflies and moths respond quickly to habitat and climatic change, and hence are valuable biodiversity indicators. In the UK, millions of species occurrence records for Lepidoptera have been gathered by two large citizen science recording schemes, of which the full potential has not been fully realized.

Analysing recording data of this nature presents unique challenges relating to their vast quantity but also associated sampling biases. Using cutting edge modelling, this project will maximise these valuable datasets to enhance our understanding of species’ phenology (flight periods), distribution and range dynamics to help inform future conservation delivery and policy for UK butterflies and moths.

The candidate will undertake new statistical model developments applied to citizen science data. The research will involve:

  • Critically assessing sampling design to determine how much data are needed to reliably estimate species’ occurrence trends – can occupancy models be used for rare species with small ranges?
  • Modelling species’ phenology from citizen science data to provide new insights on variation over space and time.
  • Applying state-of-the-art variable selection techniques to better describe drivers of species’ range and distribution change through suitable spatial and environmental covariates.

This project has been shortlisted for Aries funding. More details can be found at: