Model output from a circumpolar realisation of the Regional Ocean Modelling System (ROMS). Model was run at a horizontal resolution of 1/4 degree and 31 vertical levels. Spatial domain was circumpolar out to 30 degrees South. Forcing comes from prescribed salt and heat fluxes based on a derived climatology from Tamura et al (2008). For open water regions the Tamura data is blended with open-water heat, salt and surface stress fluxes from a monthly NCEP2 climatology.

A numerical ocean model based on the Regional Ocean Modelling System (ROMS) framework was run for the Sorsdal ice shelf region and included modifications for ice/ocean thermodynamics and mechanical pressure, following (Dinniman et al., 2007). The model domain was discretised on a polar stereographic grid with a uniform 2 km horizontal resolution. The vertical terrain-following coordinate had 31 vertical layers with a sigmoidal layer distribution to provide higher vertical resolution at the surface and bottom regions.

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 has 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. Here 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. ____ This dataset represents KRILLPODYM modelled estimates of Antarctic krill circumpolar biomass distribution for the final year of a 12-year spin up. Biomass distributions are given for each of the five key life stages outlined above. The accompanying background, model framework and initialisation description can be found in the following reference paper: Green, D. B., Titaud, O., Bestley, S., Corney, S. P., Hindell, M. A., Trebilco, R., Conchon, A. and Lehodey, P. in review. KRILLPODYM: a mechanistic, spatially resolved model of Antarctic krill distribution and abundance. - Frontiers in Marine Science