Hierarchical models for describing space-for-time variations in insect population size and sex-ratio along a primary succession

Chronosequences of glacier retreat are useful for investigating primary successions over time periods that are longer than direct observation would permit. In this context, space-for-time substitution studies have been applied to assess the effects of climate change on invertebrate assemblages. However, population dynamics of insect species following retreating glaciers has been under-investigated until now due to difficulty in applying capture-recapture methods and correctly identifying species in the field. Removal sampling methods are commonly used, but imperfect detectability is rarely accounted for in the analytical framework. In this paper we study the effects of environmental drivers of spatial, and indirectly temporal, variation in population size and sex-ratio of cold-adapted insects through a hierarchical framework for abundance. We show the importance of a metapopulation design, where samples are replicated in space and time, to model data from small and scattered populations, typically present in habitats with climate mediated selective pressure like those along glacier forelands. This scattered distribution can influence the observation or sampling process and thus species detectability.