KRILLPODYM: a mechanistic, spatially resolved model of Antarctic krill distribution and abundance

Abstract

Robust prediction of population responses to changing environments requires the integration of factors controlling population dynamics with processes affecting distribution. This is true everywhere but especially in polar pelagic environments. Biological cycles for many polar species are synchronised to extreme seasonality, while their distributions may be influenced by both the prevailing oceanic circulation and sea-ice distribution. Antarctic krill (krill, Euphausia superba) is one such species exhibiting a complex life history that is finely tuned to the extreme seasonality of the Southern Ocean. Dependencies on the timing of optimal seasonal conditions have led to concerns over the effects of future climate on krill’s population status, particularly given the species’ important role within Southern Ocean ecosystems. Under a changing climate, established correlations between environment and species may breakdown. Developing the capacity for predicting krill responses to climate change therefore requires methods that can explicitly consider the interplay between life history, biological conditions, and transport. The Spatial Ecosystem And Population Dynamics Model (SEAPODYM) is one such framework that integrates population and general circulation modelling to simulate the spatial dynamics of key organisms. Here, we describe a modification to SEAPODYM, creating a novel model – KRILLPODYM – that generates spatially resolved estimates of krill biomass and demographics. This new model consists of three major components: (1) an age-structured population consisting of five key life stages, each with multiple age classes, which undergo age-dependent growth and mortality, (2) six key habitats that mediate the production of larvae and life stage survival, and (3) spatial dynamics driven by both the underlying circulation of ocean currents and advection of sea-ice. We present the first results of KRILLPODYM, using published deterministic functions of population processes and habitat suitability rules. Initialising from a non-informative uniform density across the Southern Ocean our model independently develops a circumpolar population distribution of krill that approximates observations. The model framework lends itself to applied experiments aimed at resolving key population parameters, life-stage specific habitat requirements, and dominant transport regimes, ultimately informing sustainable fishery management.